Temperatura de brillo a partir de imágenes Landsat 8

En otros posts he escrito acerca de la determinación de la temperatura de brillo usando la banda térmica (b6) de imágenes Landsat 5. Hoy voy a describir cómo hacer lo mismo pero usando imágenes Landsat 8. Los recursos que pone a la disposición el USGS (U.S. Geological Survey) para el proceso de downloading permite a los usuarios registrados ubicar fácilmente el dataset completo (11 bandas más el archivo *.txt de la metadata) si se cuenta con un shapefile de la zona de interés. Por otra parte, este archivo de metadata es muy explícito en cuanto a los parámetros que hay que seleccionar e, incluso, existen tutoriales bastante ilustrativos acerca de como hay que proceder publicados por el propio USGS y el Blog de ESRI.

Sigue leyendo

Publicado en ArcGis10.1, Teledetección | 1 comentario

Subclassing una Class en PyQGIS

A menudo es necesario crear una clase que herede el comportamiento de una clase existente. Esto se logra mediante el “subclassing”; lo que resulta en una nueva clase que hereda todos los rasgos de la clase parental. Uno de los requerimientos más deseados en el comportamiento de QGIS es que permita manejar las geometrías 3D. En un artículo en GeoTux se afirma que QGIS o PyQGIS no saben tratar las geometrías 3D porque las consideran como shapefiles 2D. En realidad si “sabe” pero hay que construir la clase desde “cero” mediante subclassing de la clase QgsPoint.

Sigue leyendo

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

Lista de colores ampliada para una Color Ramp Shader en PyQGIS

En el post pasado se construyó una rampa de pseudocolor empleando los métodos de clase de PyQGIS. Las principales clases usadas fueron la QgsColorRampShader, la QgsRasterShader y la QgsSingleBandPseudoColorRenderer. En este post se utiliza ahora una lista con tres colores (en lugar de 2), que se seleccionan a partir de un diccionario, para colorizar los pixeles de un raster ‘singlebandgray’ con base en tres intervalos iguales que se han generado a partir de los valores mínimo y máximo del raster.

Sigue leyendo

Publicado en Código Python, PyQGIS, SIG, Software Libre | 1 comentario

Rampa de pseudocolor mediante PyQGIS

En este post voy a considerar como establecer una rampa de pseudocolor mediante PyQGIS. A pesar de lo sencillo del script, me ha costado un poco encontrar la solución debido a que el Manual Oficial del 2014 (PyQGIS Cook Book) mantiene aún sin actualizar un considerable número de métodos de clase y sus ejemplos no funcionan. El problema es especialmente notable en lo que respecta a las Raster Layers.

Sigue leyendo

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

Realce por contraste (Contrast Enhancement) con Raster Layers en PyQGIS

El realce por contraste en PyQGIS se alcanza a través de la data, atributos y métodos de la clase QgsContrastEnhancement. Para ello, es necesario usar el constructor de la clase pasando como argumento el tipo de banda usada el cual se obtiene, si se está trabajando con Raster Layers del tipo Single Band Gray, con el método ‘dataType’ de la clase QgsSingleBandGrayRenderer. Vamos a cargar en la Map View de QGIS un raster de este tipo, nos aseguramos que es capa activa y creamos una referencia al mismo a través de iface. También extraemos el código entero de la banda a través del ‘renderer’ así como su ‘dataType’.

Sigue leyendo

Publicado en Código Python, PyQGIS, SIG, Software Libre | 1 comentario

Clase QgsDistanceArea en PyQGIS

El PyQGIS Cook Book sugiere emplear la clase QgsDistanceArea cuando se requiera más precisión a la hora de determinar áreas y distancias, en lugar de emplear QgsGeometry o QgsPoint, porque estas últimas están basadas en el plano; no en el esferoide. Sin embargo, los ejemplos que refiere el PyQGIS Cook Book no han sido actualizados en cuanto a los métodos de esta clase (QgsDistanceArea). Para ejemplificar su uso se va a emplear un raster (ver imagen) donde, sólo a modo ilustrativo, se le coloca un shapefile de puntos para visualizar aquellos que van a ser manipulados desde la Python Console.

Sigue leyendo

Publicado en Código Python, PyQGIS, SIG, Software Libre | 1 comentario

Estadísticas de un raster “Single Band” usando PyQGIS

PyQGIS tiene más de 500 clases y, en consecuencia, miles de métodos de clase listados en la documentación. Por tanto, no es difícil “perderse” explorando entre ellos; sobre todo cuando los objetos de diferentes clases tienen métodos coincidentes. En este post se van determinar las estadísticas (y otra información relevante) de un raster mediante PyQGIS. Para ello, vamos a acceder a los objetos de diferentes clases utilizando los métodos adecuados.

Sigue leyendo

Publicado en PyQGIS, SIG, Software Libre | 1 comentario