Application: gvSIG desktop » gvSIG scripting

from gvsig import *

from simpleui import *

from commonsdialog import *

"""

Este script se recorre todas la features de la capa activa, y asumiendo

que tiene un campo ELEVATION, calcula la elevacion maxima y minima de 

estas, y con esto:

- Presentar un mensaje al usuario informando de cuales son

- Seleccionar las geometrias asociadas a la elevacion maxima y minima vista.

En la carpeta data puede encontrar la capa cv05_3d_polygons usada

para realizar las pruebas con este scripts.

"""

def main(): 

  #

  # Lo primero invocaremos a la funcion currentlayer() para obtener 

  # la capa activa del documento vista activo.

  # Si no hubiese un documento vista o este no tuviese una capa activa

  # currentlater() devolver? None.

  layer = currentLayer()

  if layer == None:

    # Si no somos capaces de recuperar la capa corriente, presentamos

    # un mensaje y terminamos.

    msgbox("Debera haber cargada y seleccionada una capa.")

    return

  emax = 0.0

  emin = 0.0

  fmax = None

  fmin = None

  #

  # Ahora recorreremos todas las features asociadas a la capa activa.

  for feature in layer.features():

    # 

    # Comprobaremos si la elevacion de la feature corriente es mayor

    # que la maxima calculada hasta el momento, y de ser asi?, nos la

    # guardamos como maxima.

    if feature.ELEVATION > emax :

      try:

        emax = feature.ELEVATION

        fmax = feature.getCopy()

      except Exception, ex:

        print repr(ex)

    #

    # Lo mismo pero con la elevacion minima

    if feature.ELEVATION < emin or emin == 0.0:

      emin = feature.ELEVATION

      fmin = feature.getCopy()

  #

  # Ahora seleccionamos las features en la capa.

  if fmax!=None:

    layer.getSelection().select(fmax)

  if fmin!=None:

    layer.getSelection().select(fmin)

  #

  # Por ultimo presentamos las elevaciones calculadas en un cuadro de dialogo.

  msgbox("Elevacion maxima=%s, elevacion m?inima=%s" % (emax, emin))