Python para Geoprocesamiento en SIG

• Introducción
• Relación en Python y SIG
• APIs Python para Geoprocesamiento
• Ejemplos de scripting con python para gvSIG y QGIS

El uso de Python en sistemas de información geográfica requiere por lo menos las siguientes bases de conocimiento:

• La sintaxis de Python
• Las bases de SIG:
  • Las estructuras de los datos espaciales; y,
  • La gestión de proyectos.
• Los 'Cálculos' relacionados
• La API que se ha de utilizar

Aqui dejare que los expertos (La gente de Python Ecuador) les expliquen los detalles de la sintaxis de este lenguaje.

El uso principal de python en Sistemas de información geográfica está relacionado con “Scripting”.

Recientemente se han introducido intérpretes python para la ampliación de algoritmos y herramientas desde la interfaz de los SIG. Algunos permiten combinar Python con otros lenguajes (ej. bash, C++, Fortan, java, …)

Necesitamos el módulo de scripting de gvSIG que se importará con las siguientes líneas:

from gvsig import *
from commonsdialog import *

Luego, empieza el juego de las líneas de código. Una buena práctica es el reciclaje del código, para los cual definimos las Clases y métodos necesarios.

def cal_suma(mapContext, layer, fieldname):
  canv = mapContext.getViewPort()
  rect = canv.getEnvelope().getGeometry()
  tabla = layer.getFeatureStore()
  ents = tabla.getFeatureSet().iterator()
  suma = 0
  for ent in ents:
    g = ent.getDefaultGeometry()
    if rect.intersects(g) :
      suma += int(ent.get(fieldname))
  print "Sumatorio de %s: %s" % (fieldname, suma)

De esta manera obtenemos la vista y su contexto

mapContext = currentView().getMapContext()

Por otro lado así obtenemos la capa activa y solo es cuestión de decir cual campo..

layer = currentLayer()
fieldname = 'Total'

Solo queda llamar a la función para hacer

mapContext = currentView().getMapContext()

…a partir de esto queda poner algunos “adornos” para hacerlo mas interactivo. Para eso usamos el modulo de commonsdialog

El resultado les mostraré en seguida…

De la API de PyQgis importaremos solo las clases RasterCalculator y RasterCalculatorEntry

from qgis.analysis import QgsRasterCalculator, QgsRasterCalculatorEntry

De la API de PyQgis importaremos solo las clases RasterCalculator y RasterCalculatorEntry

# Función para registrar las capas de entrada

def rasterCalcEntry(EntLst,band,refNam,bNum):
    raster = QgsRasterCalculatorEntry()
    raster.ref = refNam
    raster.raster = band
    raster.bandNumber = bNum
    EntLst.append(raster)

De la API de PyQgis importaremos solo las clases RasterCalculator y RasterCalculatorEntry

# Definimos el método que calcule el NDVI

def ndvi(R, IRC, output):
    Ent = []
    rasterCalcEntry(Ent, R, 'R@1', 1)
    rasterCalcEntry(Ent, IRC, 'IRC@1', 1)
    calc = QgsRasterCalculator(
        '(IRC@1-R@1)/(IRC@1+R@1)', 
        output, 'GTiff', band3.extent(), 
        band3.width(), band3.height(),Ent)
    calc.processCalculation()

Finalmente definimos las entradas y llamamos a las funciones definidas:

# Haremos un listado de capas disponibles 
# en el proyecto actual

lyrs = iface.mapCanvas().layers()

# Definirmos la ruta y nombre de salida

output = '../datos/geodata/Landsat/ndvi.tif'

# Llamamos a la función

ndvi(lyrs[1], lyrs[0], output)

Ahora lo veremos en acción…

Consolidaremos los valores de varios campos en uno solo.

Solo escribiremos nuestra función y la ejecutaremos desde la calculadora de campos.

def consolide(*args):
    value0 = "-"
    for arg in args:
        val = str(arg).strip()
        if val == "-":
            continue
        elif val == "":
            value0 += " "
        else:
            value0 += val
    return value0.strip("-")

• El scripting con Python permite implementar tareas automáticas que pueden dar mayor valor agregado a los análisis de información goegráfica.
• Además se pueden ampliar las funcionalidades de los diferentes software.
• Con la tendencia de los nuevos análisis en aplicaciones web, se han desarrollado librerías para implementar geoprocesos y herramientas de análisis geoespacial en línea.