Cómo grabar un ráster como imagen renderizada con PyQGIS

Cuando se carga, por ejemplo, un raster DEM a la Map View de QGIS observamos un renderizado basado en diferentes tonos de grises (del blanco al negro) porque este tipo de imagenes son ‘Single Band Gray’. Si en la propiedades de capa asignamos otro tipo de renderizado, como por ejemplo ‘Single Band Pseudocolor’ (rampa de color Spectral con 5 clases), éste persistirá en las imágenes grabadas sólo si seleccionamos en “Save As…” el modo de salida “Rendered Image” en lugar de “Raw data”.

Sigue leyendo

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

Leer HDF4 data (productos modis) utilizando PyTHON con la librería GDAL

En un post anterior ya se había señalado el soporte GDAL para la lectura de archivos con formato hdf4/hdf5 de productos MODIS. Estos productos están constituidos, a su vez, por una serie de sub datasets que pueden ser accedidos por separado, por ejemplo, en un script de PyTHON.

Sigue leyendo

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

Cómo desplegar simultáneamente en QGIS polilíneas que se solapan mediante PyQGIS

Cuando se trata de representar rutas de medios de transporte es muy probable que se use la capa de vialidad como plantilla (vectoriales tipo polilínea). Esto traerá como consecuencia que los shapefiles se solapen en diferentes tramos de la ruta haciendo imposible la visualización simultánea en estos sino modificamos su simbología por defecto. Para ello, hay que recurrir a la opción ‘offset’ de la simbología que, en PyQGIS, se accede a través de incluirla explícitamente en las opciones de diccionario del método ‘createSimple’ de los objetos de la clase QgsLineSymbolV2.

Sigue leyendo

Publicado en PyQGIS, QGIS, SIG, Software Libre | 2 comentarios

Clases texturales de suelo utilizando el raster calculator de QGIS

Existen 12 clases texturales de suelo basadas en el contenido porcentual de arena (sand), limo (silt) y arcilla (clay). Un documento del Departamento de Agricultura de USA (USDA) resume, mediante expresiones matemáticas, la forma de determinar las clases texturales. La imagen siguiente, extracto de ese documento, contiene las referidas fórmulas:

Sigue leyendo

Publicado en QGIS, Software Libre | Deja un comentario

Renderizado de capas vectoriales basado en reglas utilizando PyQGIS

Con base en una pregunta realizada en http://gis.stackexchange.com/, a continuación expondré un ejemplo de código que puede ser ejecutado en la Python Console de QGIS que permite, en este caso, transformar el estilo ‘Single Symbol’ de un vectorial de línea con tres features, en un renderizado con base en reglas diferentes para cada uno de los features. El código es una adaptación de otro que aparece publicado en el Blog de snorfalorpagus dot net sobre la simbología de capas vectoriales en plugins para QGIS.

Sigue leyendo

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

Cómo alinear rásteres producidos a partir de la rasterización de capas vectoriales individuales en QGIS

Es un problema recurrente el hecho de que muchos usuarios reportan la imposibilidad de obtener resultados coherentes en el álgebra de mapas porque sus rásteres han sido obtenidos a partir de la rasterización de capas vectoriales individuales. La razón estriba en el hecho de que los rásteres así producidos no están perfectamente alineados porque la herramienta de rasterización, a pesar de que permite fijar la resolución de manera directa, no hace lo mismo con la extensión porque el softare, por default, la determina automáticamente según sus propias pautas.

Sigue leyendo

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

Cómo obtener rásteres de ‘prueba’ con lenguaje R

Normalmente, si estamos probando o desarrollando procedimientos para SIG mediante instrucciones sencillas o lenguajes de programación, puede ser conveniente tener muestras sintéticas de los objetos a manipular con el fin de verificar de manera rápida y fácil si las acciones sobre ellos producen los resultados esperados. En el caso de rásteres, aunque QGIS tiene a la disposición una serie de herramientas que nos podrían servir para muchas opciones, su versatilidad y potencia dista mucho de la que permite el empleo de la librería raster que se encuentra disponible en el lenguage R. Esta librería requiere a su vez la librería rgdal si se desean manipular los formatos ‘GTiff’. Por tanto, es conveniente tenerla también instalada.

Sigue leyendo

Publicado en Lenguaje R, SIG, Software Libre | 1 comentario