Hola amigos,
estoy haciendo mis primeros pinitos con el python y gvsig. Me ha surgido un
problema con un script que estoy haciendo.
Es muy simple, consiste en intersectar los polígonos de una capa y generar
un shapefile nuevo con las áreas de solape, para identificar los errores.
Necesito generar esta capa, por lo que no me vale que me señale los solapes
como hace la herramienta de topología. Además es una excusa perfecta para
practicar con gvsig el tema de la programación.
El script: es el siguiente;
El problema me surge una vez que tengo las intersecciones creadas y quiero
generar el nuevo shape. He probado pero no lo consigo.
¿Alguién me puede orientar?
Muchas gracias.
from gvsig import *
from geom import *
def main():
path = "D:\\JJV_2E\\CURSO\\salida.shp"
layer = currentLayer() # active layer
### proceso
schema = layer.getSchema()
newSchema = createSchema(schema)
CRS = currentView().getProjection()
output = createShape(newSchema, path, CRS=CRS,geometryType=POLYGON)
features = layer.features()
for feature in features:
g1 = feature.geometry()
for feature2 in features:
values = feature.getValues()
g2 = feature.geometry()
intersect = g1.intersection (g2)
_______________________________________________
gvSIG_desarrolladores mailing list
gvSIG_desarrolladores@listserv.gva.es
Para ver histórico de mensajes, editar sus preferencias de usuario o darse de
baja en esta lista, acuda a la siguiente dirección:
http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_desarrolladores
