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