gvSIG blog

En anteriores entradas hemos visto cómo crear un script muy sencillito en gvSIG 2.0.

En esta entrada os traemos un vídeo en el que implementaremos un script que crea una capa shape de polígonos que se obtienen de aplicar un buffer a las geometrías de una capa de puntos. La distancia del buffer se obtiene de uno de los campos de la capa de puntos.

Los pasos que se siguen en el script son:

• Obtener la capa de entrada
• Crear la definición de datos (Schema) de la nueva capa
• Obtener los datos de creación de la nueva capa
• Crear la nueva capa
• Recorrer los fenómenos de la capa de entrada (puntos)
• Generar las nuevas geometrías
• Añadir las nuevas geometrías a la nueva capa
• Terminar la edición de la capa nueva
• Liberar el conjunto de fenómenos de la capa de entrada

El código con comentarios es el siguiente

from gvsig import *

def main():
  #Obtenemos la capa activa
  layer = currentLayer()

  #Creamos la definición de los datos
  schema = createSchema(layer.getSchema())
  schema.modify()

  #Obtenemos la proyeccion
  crs = currentProject().getProjectionCode()

  #Definimos la ruta de la nueva capa
  ruta = "/home/victor/carto/influencia.shp"

  #Creamos la capa
  newLayer = createShape(
    schema,
    ruta,
    CRS = crs,
    geometryType = SURFACE
  )

  #Obtenemos los fenomenos de la capa de entrada
  features = layer.features()

  #Recorremos los fenomenos 
  for feature in features:

    # Obtenemos el valor del area de influencia
    influencia = feature.get("Influencia")

    #Creamos la nueva geometria aplicando un buffer sobre la geometria del fenomeno de entrada
    geom = feature.geometry().buffer(influencia)

    # Obtenemos los valores de los atributos del fenomeno de entrada
    values = feature.getValues()

    # Insertamos la nueva geometria en la capa nueva
    values["GEOMETRY"] = geom
    newLayer.append(values)

  # Terminamos la edicion de la nueva capa guardando los cambios
  newLayer.commit()

  # Liberamos el conjunto de fenomenos de la capa de entrada
  features.dispose()

Debes tener en cuenta que para que funcione es necesario tener la última versión de gvSIG 2.0 y la extensión Scripting Framework instaladas. Además, la vista debe tener una capa activa y esta capa debe tener un campo llamado Influencia de tipo double. Y au mejor, podéis modificar el script para que se adapte a vuestros datos.

Espero que os haya gustado y recordad; nos vemos en las próximas Jornadas internacionales de gvSIG!

Una de las novedades que introduce la versión 2.0 de gvSIG es la inclusión de un pequeño entorno para el desarrollo y ejecución de scripts que interactúen con la aplicación. Estos scripts nos permiten automatizar pequeñas tareas o añadir alguna funcionalidad que echemos en falta y que seamos capaces de implementar.

El lenguaje elegido desde el proyecto gvSIG para dar soporte al scripting es Jython , que es una implemetación del lenguaje Python que se ‘ejecuta’ en la máquina virtual de Java (JVM), lo que permite que esté completamente integrado con gvSIG.

La instalación por defecto de gvSIG no lleva de base el plugin que nos aportara el soporte de scripting, así que lo primero que tendremos que hacer será instalarlo, puedes consultar cómo en la documentación del administrador de complementos de gvSIG .

El plugin de scripting aporta dos herramientas, por un lado un editor y organizador de nuestros scripts, Scripting Composer, y por otro un lanzador de estos, Scripting Launcher. Estas dos herramientas están disponibles en el menú Herramientas -> Scripting así como en dos botones en la barra de botones de gvSIG.

Como podéis ver en la imagen el Composer se divide en 3 zonas principales, un explorador de archivos, desde el que podemos ver los scripts que tenemos disponibles, un área de edición, y un área de notificaciones dónde podremos ver las incidencias que ocurran durante la ejecución de nuestros scripts.

Una de las ventajas que ofrece que la extensión de scripting esté integrada en gvSIG es que nos permite acceder a los documentos que tenemos cargados en la aplicación e interactuar con ellos, lo que junto con el refactoring de los APIS de gvSIG 2.0 y las librerías que se han creado, permite crear scripts de forma sencilla. Veámoslo con un ejemplo.

Supongamos que tenemos una capa cargada cuya estructura de datos tiene un campo llamado ELEVATION y nosotros queremos obtener los valores máximo y mínimo de este campo. Para ello debemos recorrer todas las features de la capa y comprobar los valores de este campo, almacenar el valor máximo y el mínimo y mostrarlos mediante una ventana.

from gvsig import *
from commonsdialog import *

def main(): 
  layer = currentLayer()
  if layer == None:
    msgbox("Deberá haber cargada y 
             seleccionada una capa.", "AVISO", 1)
    return

  emax = 0.0
  emin = 0.0

  for feature in layer.features():
    if feature.ELEVATION > emax :
      emax = feature.ELEVATION
    if feature.ELEVATION < emin or emin ==0.0:
      emin = feature.ELEVATION
  msgbox("Elevación máxima=%s y mínima=%s" % (emax, emin),
         "Elevation", 0)

El resultado de ejecutar el script es el que puede verse en la imágen siguiente.

Espero que os haya gustado y que os anime a empezar a programar vuestros propios scripts en Python y si es así, y os volvéis algo pythonicos tal vez os guste leer el zen de python o mejor aún leerlo ejecutando vuestro propio script en gvSIG. ¿Nos contáis cómo?