9. Trabajando con capas: FLayer¶

Trabajando con capas FLayer

createShape(definition[, filename=None, geometryType=None, CRS=None, prefixname="tmpshp"])¶

Creates a new shape layer. If geometryType is None, will take geometry type from the definition. If parameter geometryType and the geometry type inside the definition are different, raises an error.

Parámetros:

definition – Layer schema :type definition: DefaultFeatureType
filename (str) – Full path
geometryType –

type geometryType:

GeometryType
CRS (str) – CRS Code :param str prefixname: first part of the temp name

La interfaz FLayer viene implementada en varias capas de uso común como son las capas vectoriales o raster.

9.1. Capas vectoriales: FLyrVect¶

Unas de las operaciones más importantes es el trabajo con capas vectoriales FLyrVect. Hemos querido simplificar al máximo posible el trabajo con ellas, sobretodo si eres un nuevo desarrollador.

Ya hemos visto anteriormente la generación de esquemas que pueden ser usados en capas vectoriales. La única condición de estos esquemas es que contengan un campo GEOMETRY con un tipo de geometría establecido.

9.1.1. Crear capa vectorial¶

Para la generación de capas hemos simplificado al máximo la función :py:func:createShape(definition). El único parámetro obligatorio será el esquema de la capa con su campo geometría.

Por ejemplo, lo mínimo necesario para crear una capa es:

from gvsig import *
from gvsig.geom import *

def main(*args):

        schema = createFeatureType() # DefaultFeatureType

        schema.append("GEOMETRY", "GEOMETRY")
        schema.get("GEOMETRY").setGeometryType(POINT, D2)
        shape = createShape(schema)

De esta forma, si la ruta o el nombre del fichero no es importante para nosotros, algo muy común a la hora de generar capas para resultados de geoprocesos o similares, la función se encargará de asignarle una ruta en una carpeta temporal y un nombre aleatorio (basado en códigos de tiempo).

Otra forma sería usar directamente la funcion createLayer():

layer = createLayer(schema=schema,
                    servertype="FilesystemExplorer",
                    layertype="Shape",
                    shpFile="/home/osc/temp/test1.shp",
                    CRS="EPSG:25830",
                    geometryType=geom.POINT
                    )

Así en el script anterior podemos añadir las líneas:

print "Nombre: ", shape.getName()
print "Ruta: ", shape.getDataStore().getFullName()

Que darán como resultado el nombre y ruta de la capa, similar a:

Nombre:  tmpshp-57afa1381035
Ruta:  C:\Users\Oscar\AppData\Local\Temp\tmp-andami\tmpshp-57afa1381035.shp

Si necesitas crear varias capas y necesitas poder diferenciarlas, por ejemplo, al enlazar una serie de operaciones, podemos establecer el parámetro prefixname el cual modificará la primera parte del nombre, pero seguirá creando todo el resto de ruta temporal y asegurándonos que será un nombre único:

from gvsig import *
from gvsig.geom import *

def main(*args):

        schema = createFeatureType() # DefaultFeatureType

        schema.append("GEOMETRY", "GEOMETRY")
        schema.get("GEOMETRY").setGeometryType(POINT, D2)
        shape1 = createShape(schema, prefixname="resultado")
        shape2 = createShape(schema, prefixname="errores")
        shape3 = createShape(schema, prefixname="puntos")

        print "Nombre: ", shape1.getName(), "\tRuta: ", shape1.getDataStore().getFullName()
        print "Nombre: ", shape2.getName(), "\tRuta: ", shape2.getDataStore().getFullName()
        print "Nombre: ", shape3.getName(), "\tRuta: ", shape3.getDataStore().getFullName()

Por consola muestra el nombre de las capas como indicábamos:

Nombre:  resultado-57afa4f612c9         Ruta:  C:\Users\Oscar\AppData\Local\Temp\tmp-andami\resultado-57afa4f612c9.shp
Nombre:  errores-57afa4f613a6   Ruta:  C:\Users\Oscar\AppData\Local\Temp\tmp-andami\errores-57afa4f613a6.shp
Nombre:  puntos-57afa4f61446    Ruta:  C:\Users\Oscar\AppData\Local\Temp\tmp-andami\puntos-57afa4f61446.shp

9.1.2. Modificar esquema de una capa¶

El siguiente script modificará el esquema de una capa. Para ello tenemos que crear un esquema nuevo basado en el anterior de la capa mediante createFeatureType(layer_schema), realizar las modificaciones y actualizar la capa:

from gvsig import *
from gvsig import geom

def main(*args):
        """Updating schema of existent layer"""

        layer = currentLayer()

        schema = layer.getSchema()
        newschema = createSchema(schema)
        newschema.append("ID2", "STRING")

        layer.edit()
        layer.update(newschema)
        layer.commit()

9.1.3. Operaciones con entidades¶

Una vez creada la capa nueva o accedido a una de ellas con currentLayer() o view.getLayer("Name"), podemos acceder a sus entidades mediante el método .features(), tal y como explicamos en la guía de Acceso a datos.

Lo siguiente que haremos es añadir datos a esta capa vectorial. Para ello ponemos la capa en modo de edición mediante layer.edit() y agregamos las entidades con layer.append(args):

from gvsig import *
from gvsig.geom import *

def main(*args):

        schema = createFeatureType() # DefaultFeatureType

        schema.append("ID", "INTEGER", 5)
        schema.append("NAME", "STRING", 10)
        schema.append("GEOMETRY", "GEOMETRY")
        schema.get("GEOMETRY").setGeometryType(POINT, D2)

        shape = createShape(schema, prefixname="resultado")


        print "Nombre: ", shape.getName(), "\tRuta: ", shape.getDataStore().getFullName()

        shape.edit()
        # Setting arguments
        shape.append(ID=1, NAME="Valencia", GEOMETRY=createPoint2D(10, 10))
        # Diccionary
        shape.append({"ID": 2, "NAME": "Paris", "GEOMETRY":createPoint2D(15, 15)})
        shape.commit()

        currentView().addLayer(shape)

Otro ejemplo añadiendo entidades, usando también la forma con la que se haría desde Java:

import gvsig
reload(gvsig)
from gvsig import *
from gvsig import geom

from org.gvsig.fmap.dal.feature import FeatureStore
def main(*args):

        # Creating new layer
        schema = createSchema()
        schema.append("ID", "INTEGER")
        schema.append("NAME", "STRING", 10)
        schema.append("GEOMETRY", "GEOMETRY")
        schema.get('GEOMETRY').setGeometryType(geom.POINT,geom.D2)

        layer = createShape(schema, prefixname="points_layer")

        # Insert with newfeature
        store = layer.getFeatureStore()

        newfeature = store.createNewFeature()
        newfeature.set("ID",1)
        newfeature.set("NAME","Feature1")
        newfeature.set("GEOMETRY", geom.createPoint(geom.D2, 1,2))

        layer.edit(FeatureStore.MODE_APPEND) #solo para capas recien creadas
        store.insert(newfeature)
        layer.commit()

        # Insert with append
        layer.edit()
        layer.append(ID=2,NAME='Feature2',GEOMETRY=geom.createPoint(geom.D2, 5, 3))

        layer.append({'ID':3,'NAME':'Feature2','GEOMETRY':geom.createPoint(geom.D2, 3, 3)})
        layer.append({'ID':4,'NAME':'Feature2','GEOMETRY':geom.createPoint(geom.D2, 2, 1)})
        layer.append({'ID':5,'NAME':'Feature3','GEOMETRY':geom.createPoint(geom.D2, 2, 6)})
        layer.append({'ID':6,'NAME':'Feature3','GEOMETRY':geom.createPoint(geom.D2, 6, 2)})
        layer.append({'ID':7,'NAME':'Feature3','GEOMETRY':geom.createPoint(geom.D2, 2, 7)})
        layer.commit()

        # Add layer to the view
        currentView().addLayer(layer)

        print "Features info"
        for l in layer.features():
                print l

Si al final del script anterior añadimos las siguientes líneas, veremos un ejemplo para eliminar entidades:

features = layer.features("ID < 6") #DefaultFeatureSet

layer.edit()
print type(layer)
print features, type(features)
for i in features:
        features.delete(i)

layer.commit()

Para modificar los valores de las entidades que contiene nuestra capa:

layer.edit()

for i in features:
        print i
        c = i.getEditable()
        c.set("NAME", "Modified_4")
        features.update(c)

layer.commit()

Se puede realizar copias de entidades (features) y poder modificarlas después en su capa original.

Ejemplo: Extraemos ciertas entidades de una capa que contiene un Campo1 de tipo Long. Estas entidades las copiamos a una lista. Después, modificamos estas entidades y las volvemos a modificar sobre la capa inicial:

from gvsig import *

def main(*args):
    layer = currentLayer()
    features = layer.features('Campo1>2',sortBy="Campo1",asc=True)
    lista = []
    for f in features:
        print f
        copia = f.getCopy()
        print type(copia)
        lista.append(copia)

    print len(lista)
    layer.edit()
    for i in lista:
        value = i.get('Campo1')+0.01
        i = i.getEditable()
        i.set('Campo1', value)
        print "new value", i.get('Campo1'), type(i)
        featureSet = layer.features()
        layer.features().update(i)
    layer.commit()

Tabla de Contenidos

Tema anterior

Próximo tema

Esta página

9. Trabajando con capas: FLayer¶

9.1. Capas vectoriales: FLyrVect¶

9.1.1. Crear capa vectorial¶

9.1.2. Modificar esquema de una capa¶

9.1.3. Operaciones con entidades¶