<div dir="ltr">Te faltan al menos dos cosas. Escribo de cabeza así que puede haber algún error:<div><br></div><div>* Comitear los cambios a la capa destino. Fuera del for tienes que hacer un: output.commit()</div><div>* Añadir la nueva geometría variable "intersect" en tu casa a la capa destino: output.append({"GEOMETRY":intersect})<br></div><div><br></div><div>* Además sólo deberías crear una nueva geometría cuando realmente se produzca una intersección</div><div><br></div><div>if g1.intersects(g2):</div><div> intersect = g1.intersection (g2)</div><div> output.append({"GEOMETRY":intersect})</div><div><br></div><div>Y por último no estoy seguro de que el append permita crear un nuevo feature sólo pasándole el valor de GEOMETRY, si no es así, el código sería por ejemplo:</div><div><br></div><div><div>if g1.intersects(g2):</div><div> values = feature.getValues()</div><div> intersect = g1.intersection (g2)</div><div> values["GEOMETRY"] = intersect<br></div></div><div> output.append({"GEOMETRY":intersect})<br></div><div><br></div><div>Esto te crearía el nuevo registro con los valores de una de los registros que intersecan con geometría sólo la parte que interseca.</div><div><br></div><div>Échale un ojo a este script que te puede dar alguna pista:</div><div><br></div><div><a href="http://outreach.gvsig.org/scripts/buffer-intersection-condition-and-new-field-spatial-info">http://outreach.gvsig.org/scripts/buffer-intersection-condition-and-new-field-spatial-info</a><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">El 20 de abril de 2015, 12:24, Juan Pedro Ruiz Castellano <span dir="ltr"><<a href="mailto:juanpedro.ruizcastellano@gmail.com" target="_blank">juanpedro.ruizcastellano@gmail.com</a>></span> escribió:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hola amigos,<br>
estoy haciendo mis primeros pinitos con el python y gvsig. Me ha surgido un problema con un script que estoy haciendo.<br>
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.<br>
<br>
El script: es el siguiente;<br>
El problema me surge una vez que tengo las intersecciones creadas y
quiero generar el nuevo shape. He probado pero no lo consigo.<br>
<br>
¿Alguién me puede orientar?<br>
<br>
Muchas gracias.<br>
<br>
from gvsig import *<br>
from geom import *<br>
<br>
def main():<br>
path = "D:\\JJV_2E\\CURSO\\salida.shp"<br>
layer = currentLayer() # active layer<br>
### proceso<br>
schema = layer.getSchema()<br>
newSchema = createSchema(schema)<br>
CRS = currentView().getProjection()<br>
output = createShape(newSchema, path, CRS=CRS,geometryType=POLYGON)<br>
features = layer.features()<br>
for feature in features:<br>
g1 = feature.geometry()<br>
for feature2 in features:<br>
values = feature.getValues()<br>
g2 = feature.geometry()<br>
intersect = g1.intersection (g2)<br>
</div>
<br>_______________________________________________<br>
gvSIG_desarrolladores mailing list<br>
<a href="mailto:gvSIG_desarrolladores@listserv.gva.es">gvSIG_desarrolladores@listserv.gva.es</a><br>
Para ver histórico de mensajes, editar sus preferencias de usuario o darse de baja en esta lista, acuda a la siguiente dirección: <a href="http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_desarrolladores" target="_blank">http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_desarrolladores</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>Francisco Puga</div><div>iCarto | Innovación, Cooperación, Cartografía y Territorio S.L.</div><div><a href="http://www.icarto.es/" target="_blank">http://www.icarto.es/</a></div><div><br></div><div>c/ Rafael Alberti nº 13 – 1º D</div><div>15008 A Coruña</div><div>Galicia (Spain)</div><div>+34 881927808</div><div><br></div><div>Este correo electrónico contiene información estrictamente confidencial y es de uso exclusivo del destinatario, quedando prohibida a cualquier otra persona su revelación, copia, distribución, o el ejercicio de cualquier acción relativa a su contenido. Si ha recibido este mensaje por error, por favor conteste a su remitente mediante correo electrónico y proceda a borrarlo de su sistema.</div><div><br></div><div>Sus datos personales serán tratados de forma confidencial y no serán cedidos a terceros ajenos a ICARTO. En cualquier caso, podrá ejercer los derecho de oposición, acceso, rectificación y cancelación de acuerdo con lo establecido en la Ley Orgánica 15/99, de 13 de diciembre, de Protección de Datos de Carácter Personal dirigiéndose a Innovación, Cooperación, Cartografía e Territorio, SL. (ICARTO) en la dirección postal a C/ Rafael Alberti, nº 13, 1ºD, 15.008 – (A Coruña).</div></div></div>
</div>