Creando buffers cuadrados alrededor de puntos con PyQGIS

Crear buffers cuadrados alrededor de puntos usando PyQGIS no es muy difícil. En la clase QgsGeometry existe el método (‘buffer’) y para ello sólo hay que seleccionar previamente todos los rasgos. Sin embargo, el resultado inicial es el ya conocido buffer circular.

Seguir leyendo

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

Determinando el punto de un “Gran Círculo” con la menor distancia a un punto exterior

Cuando se trabaja con coordenadas proyectadas referidas a pequeñas distancias, la separación más corta entre dos puntos (expresada mediante una línea adecuadamente densificada) probablemente no difiera visualmente de una recta. Sin embargo, cuando éstas se incrementan, como la representación de la superficie de la tierra corresponde a la de un elipsoide, la menor distancia entre dos puntos, en la mayoría de las proyecciones que comunmente usamos, debería visualizarse como una curva; no como una recta. Estas curvas corresponden a lo que denominamos “Grandes Círculos”.

Seguir leyendo

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

Estadística zonal con puntos espaciados regularmente dentro de los rasgos (features) de un vectorial tipo polígono con PyQGIS

En el post anterior se refirió la manera de cómo crear puntos espaciados regularmente dentro de los rasgos (features) de un vectorial tipo polígono con PyQGIS. Vamos a aprovechar este hecho para, con un ráster base como capa subyacente, tome los valores de la réticula de éste y, mediante el método ‘identify’ de QgsRasterDataProvider muestree todos los values por feature con el fin de determinar parte de sus estadísticas locales (en este caso número de valores, media y desviación estándar).

Seguir leyendo

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

Creando puntos espaciados regularmente dentro de los rasgos (features) de un vectorial tipo polígono con PyQGIS

La creación o generación de puntos regularmente espaciados es un procedimiento que puede realizarse utilizando las herramientas que pone a su disposición PyQGIS. Si a eso le añadimos la posibilidad de que los puntos esten sólo circunscritos dentro de los rasgos de un archivo vetorial tipo polígono y que el espaciamiento sea el de un ráster base, entonces podrían muestrearse los valores del ráster para hacer estadística zonal.

Seguir leyendo

Publicado en PyQGIS, SIG, Software Libre | 2 comentarios

Cortando (clipping) un ráster mediante un objeto de la clase QgsRasterBlock en PyQGIS

En el post anterior se desarrolló un script para seleccionar todos los valores de un ráster mediante un objeto de la clase QgsRasterBlock en PyQGIS. Cualquier sub región de dicho ráster puede ser establecida simplemente con indicar su vértice superior izquierdo y su vértice inferior derecho. Por tanto, es lógico suponer que, una vez seleccionada dicha región, pueda ser escrita como un nuevo ráster. Esto sería equivalente a lo que produce la herramienta ‘Clipper’.

Seguir leyendo

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

Obtener todos los valores de un ráster mediante un objeto de la clase QgsRasterBlock en PyQGIS

En otros posts se ha usado el método ‘identify’ de QgsRasterDataProvider para acceder a valores individuales o a todos ellos de un ráster. Tiene la particularidad de que, al final, permite generar una capa de puntos con patrón de retícula si el muestreo tiene como objetivo todos los valores. Sin embargo, existe un método de QgsRasterDataProvider más eficiente: ‘block’. Este genera un objeto de la clase QgsRasterBlock que, a través de su método ‘value’, hace posible la obtención de todos ellos en un sólo paso o modificarlos mediante el método ‘setValue’.

Seguir leyendo

Publicado en PyQGIS, SIG, Software Libre | 1 Comentario

Rasterizando por valores mínimos, máximos o promedio de un campo en un vectorial de puntos

La situación que se tiene en la imagen siguiente es la de un pequeño ráster de 3 x 3 (resolución 73.9887 m x 73.9887 m) en cuya extensión existen 17 puntos generados aleatoriamente y donde en algunas celdas se observa que hay más de un punto. Además, en la tabla atributiva del vectorial de puntos se observa que al seleccionar los dos primeros, estos aparecen (en amarillo) en celdas no contiguas en el Map Canvas reflejando que no existe ninguna sistematización en la asignación de los id.

Seguir leyendo

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