Scripting en gvSIG: Coordenadas de los extremos de líneas en la tabla

El otro día un compañero en las listas de usuarios de gvsig de nos planteó una necesidad que tenía: añadir a la tabla de datos 4 campos nuevos que contuvieran las coordenadas X, Y de inicio y fin de las lineas que tenía su shape.

La solución sin usar programación la podéis seguir en dos emails[1] [2], explicada gracias a Gustavo Agüero Córdoba . 

Me parecía un buen ejemplo sencillo para exponer y aquí lo tenéis, simplemente tenéis que tener la capa activa donde queréis realizar los cambios,  y ejecutarlo. Sencillo pero que te puede ser muy útil, o servirte de base para algo que utilices.

import gvsig
from geom import *

def main():
    """Obtiene los vertices finales de 
        lineas y los agrega a tabla"""
    #Nuevo nombre para los campos X,Y
    campoX = "X_final"
    campoY = "Y_final"
    campoXi = "X_inicial"
    campoYi = "Y_inicial"
    layer = gvsig.currentLayer()
    
    schema = layer.getSchema()
    schema.modify()
    schema.append(campoXi,"DOUBLE", 30)
    schema.append(campoYi,"DOUBLE", 30)
    schema.append(campoX, "DOUBLE", 30)
    schema.append(campoY, "DOUBLE", 30)
    layer.edit()
    layer.updateSchema(schema)
    layer.commit()
    features = layer.features()
    for line in features:
        geom = line.geometry()
        vertex = geom.getVertex((geom.numVertices-1))
        vertexi = geom.getVertex(0)
        x,y = vertex.getX(), vertex.getY()
        xi,yi = vertexi.getX(), vertexi.getY()
        line.edit()
        line.set(campoX, x)
        line.set(campoY, y)
        line.set(campoXi, xi)
        line.set(campoYi, yi)
        layer.update(line)
        
    layer.commit()