[Gvsig_usuarios] coordenadas iniciales y finales de polígonos

[Mario Carrera]

Hola Edu,

hemos estado haciendo pruebas con el shp que nos has enviado y 
efectivamente da error dicho geoproceso, incluso hemos probado con la 
versión 1.12 por si hubo algún error en la migración y salta el mismo 
error.

Hemos estado probando también otras posibles soluciones, con otros 
geoprocesos, pero no sale una solución buena.

Para extraer las coordenadas de un polígono sí que hay un script, pero 
no creo que te de mucho juego para obtener lo que quieres. Te adjunto el 
script por si quieres hacer pruebas.

Aquí puedes ver cómo crearlo: 
https://www.youtube.com/watch?v=7c_6KetDOAM (pero en tu caso debes 
copiar el texto que te adjunto)

Saludos,
   Mario



El 26/07/17 a las 11:30, Mario Carrera escribió:
>
> Hola Edu,
>
> la forma óptima sería con el de "eje de polígono". ¿Puedes adjuntarnos 
> un shp con algunos polígonos con los que hayas probado? Así haremos 
> zquí alguna prueba.
>
> Saludos,
>  Mario
>
>
>
> El 26/07/17 a las 09:37, Edu Requena Cerdan escribió:
>> Buenos días,
>>
>> Tengo una duda sobre como se obtienen las coordenadas finales e 
>> iniciales de unos polígonos. Para que se hagan una idea, los 
>> polígonos hacen referencia a caminos. He intentado varios procesos:
>>
>>
>> - en primer lugar intenté obtener el eje de los polígonos por medio 
>> de la herramienta "eje de polígono (medial axis)", pero no obtuve 
>> ningún resultado.
>>
>>
>> -más tarde, converti los polígonos a lineas por medio de la 
>> herramienta "convertir polígonos en polilíneas".
>>
>>
>> -seguidamente, probé la herramienta "muestrear extremos de líneas" 
>> pero no me permite seleccionar ese proceso.
>>
>>
>> -por otro lado,y lo único que me ha servido, ha sido la herramienta 
>> "convertir líneas en puntos equiespaciados", pero como es de 
>> esperar obtengo muchos registros según la tolerancia aplicada.
>>
>>
>> No se si habrá alguna manera mejor o mas sencilla para crear una capa 
>> de puntos sobre el inicio y final de capa polígono, ya que son muchos 
>> registros para hacerlo de manera individualizada.
>>
>>
>> Agradezco vuestra atención.
>>
>>
>>  Un saludo, Eduardo.
>>
>>
>> <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> 
>> 	Libre de virus. www.avg.com 
>> <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> 
>>
>>
>>
>>
>> _______________________________________________
>> gvSIG_usuarios mailing list
>> gvSIG_usuarios en 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:
>>
>> https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios
>
>
>
> _______________________________________________
> gvSIG_usuarios mailing list
> gvSIG_usuarios en 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:
>
> https://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://listserv.gva.es/pipermail/gvsig_usuarios/attachments/20170726/23333196/attachment.html>
------------ próxima parte ------------
from gvsig import *
from gvsig.commonsdialog import *
 
def main(*args):
    """Read wkt"""
    sel = currentLayer().getSelection()
    pfile = str(saveFileDialog("Seleccionar fichero texto de salida")[0])
    f = open(pfile,'w')
    for s in sel:
        f.write("\n\n================\n\n")
        f.write("\n"+str(s.getValues())+"\n\n")
        g = str(s.geometry().convertToWKT())
        f.write(g)
    f.close()
    print pfile