Evaluando áreas de intersección con QgsDistanceArea y coordenadas geográficas

Con el uso combinado de las clases QgsDistanceArea y QgsGeometry de PyQGIS es posible determinar áreas de intersección cuyas coordenadas vienen expresadas como geográficas. En la imagen siguiente, se pueden observar tres capas tipo polígono, proyectadas en WGS84, donde la intersección corresponde a la capa de color amarillo ocre.

Seguir leyendo

Publicado en PyQGIS, SIG, Software Libre | 1 Comentario

Polígonos con “agujeros” en PyQGIS

Para considerar los polígonos con agujeros en PyQGIS sólo basta con tener una lista de listas. Sin embargo, a pesar de que puedan ser visualizados en el Map Canvas de QGIS, eso no significa que la geometría sea válida. Por ello, es conveniente incluir el método ‘isGeosValid’ de QgsGeometry para tal fin.

Seguir leyendo

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

Intersecciones con shapely en shapefiles con número de features elevado

Una de las desventajas que tienen los lenguajes interpretados es la velocidad de proceso y eso de pone de manifiesto cuando nos toca hacerlo con shapefiles con gran número de features; tal como en el del ejemplo siguiente. El shapefile tiene 20.730 features y el número de combinaciones, tomadas de 2 en 2 para evaluar las primeras intersecciones, es de 214.856.085.

Seguir leyendo

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

Usando Fiona y Shapely en la Python Console de QGIS (GNU/Linux Debian)

Hace unos días publiqué un primer post acerca de ciertos módulos complementarios (GeoPandas) a la labor que puede realizarse con PyQGIS. Hoy me tocará escribir sobre Fiona y Shapely.

Seguir leyendo

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

Cómo producir rejillas (grid) triangulares mediante PyQGIS

En posts pasados (1,2) se consideró la manera de producir rejillas hexagonales mediante PyQGIS y en este le toca el turno a las triangulares. Aunque pueda parecer extraño, a menor número de lados mayor complejidad.

Seguir leyendo

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

Cortando una rejilla (grid) hexagonal creada con PyQGIS

En el post anterior se consideró la creación de una rejilla hexagonal. Esta contenía unos “salientes” laterales y en las partes superior e inferior de la rejilla. Además, el punto de referencia usado para su confección luce desplazado por la necesidad de encajar las filas de la manera adecuada.

Seguir leyendo

Publicado en PyQGIS, SIG, Software Libre | 1 Comentario

Cómo producir rejillas (grid) hexagonales mediante PyQGIS

En el post anterior se consideró la manera de obtener “buffers” con diferentes formas poligonales alrededor de puntos (centroides). En este post se va a aprovechar este hecho para generar una rejilla hexagonal a partir de un punto de referencia dado.

Seguir leyendo

Publicado en PyQGIS, SIG, Software Libre | 2 comentarios