Librerías liblas y laspy para manejar datos LIDAR

Los archivos tipo LAS corresponden a un formato de archivo binario (con licencia pública) usado para el intercambio de nubes de puntos tridimensionales y que ha sido desarrollado, en un principio, para el intercambio de datos LIDAR (Light Detection and Ranging); los cuales son obtenidos por el barrido de un LASER aerotransportado. Además, este formato soporta el intercambio de cualquier tupla de datos tridimensionales x,y,z que a bien tenga considerarse.

Seguir leyendo

Publicado en LIDAR, SIG, Software Libre | Deja un comentario

Polígonos de Thiessen con el algoritmo de Fortune en PyQGIS

En el post pasado se consideró un posible bug de QGIS que no hacía posible obtener los polígonos de Thiessen (o Voronoi como también se les conoce) cuando el número de puntos de control era bajo. En la búsqueda de soluciones se corroboró que ArcMap no presentaba el problema y también se encontró una opción programada para Python que empleaba el algoritmo de Fortune.

Seguir leyendo

Publicado en PyQGIS, QGIS, SIG, Software Libre | Deja un comentario

Polígonos de Thiessen en ArcMap y QGIS

Los polígonos de Thiessen (o Voronoi) se obtienen por un método de interpolación que se basa en la distancia euclidiana y resultan de unir puntos entre sí mediante las mediatrices de los segmento de unión. Las intersecciones de estas mediatrices determinan una serie de polígonos en un espacio bidimensional alrededor de un conjunto de puntos de control, de manera que el perímetro de los polígonos generados sea equidistante a los puntos vecinos y designan su área de influencia.

Seguir leyendo

Publicado en ArcGis10.1, QGIS | 2 comentarios

GeoJSON como objetos python y las librerías json y OGR

GeoJSON es un formato estándar abierto diseñado para representar elementos geográficos sencillos, junto con sus atributos no espaciales, basado en JavaScript Object Notation. El formato es ampliamente utilizado en aplicaciones de cartografía en entornos web al permitir el intercambio de datos de manera rápida, ligera y sencilla.

Seguir leyendo

Publicado en Código Python, PyQGIS, QGIS, SIG, Software Libre | Deja un comentario

Asignación de valores ráster más cercanos a puntos sobre ‘None’ values con PyQGIS

En el post anterior se consideró la conformación de una memory layer de puntos en las celdas del perímetro de una máscara ráster. La máscara es equivalente a “recortar” un ráster por una geometría tipo polígono y asigna como valores ‘nodata’ a los píxeles que se encuentran más allá del perímetro de la geometría.

Seguir leyendo

Publicado en PyQGIS, QGIS, SIG, Software Libre | Deja un comentario

Shapefile de puntos en las celdas del perímetro de una máscara ráster

Supongamos que tenemos rasgos (features) de un archivo vectorial sobre celdas que albergan ‘nodata values’ o estén más allá de la extensión de un ráster y, aún así, queremos asociarlas a algún valor del ráster. Lo más lógico es que sea el que corresponda a la celda de éste más cercana al punto.

Seguir leyendo

Publicado en Código Python, GDAL, SIG, Software Libre | 1 Comentario

Distancia entre el centroide y el punto más lejano de un polígono con PyQGIS

Encontrar la distancia entre el centroide y el punto más lejano de un polígono no es una tarea difícil con PyQGIS; aún si éste está constituido por numerosos features y el cálculo tiene que realizarse para todos. Para ello, se guardan sendas listas con los centroides y los polígonos y, usando un valor de prueba para la distancia máxima que inicialmente es cero, se determinan todas ellas para cada centroide y los puntos respectivos del polígono; seleccionando la mayor, el punto respectivo y su índice. Estos puntos, conjuntamente con el centroide, se usan para configurar una lista de QgsPoint con la sintaxis propia de polilíneas para ser conformadas como una memory layer.

Seguir leyendo

Publicado en PyQGIS, QGIS, SIG, Software Libre | Deja un comentario