[Gfoss] GRASS esportare layer 0

[Marco Guiducci]

----- Messaggio originale -----
> Da: Markus Neteler <neteler a osgeo.org>
> A: Marco Guiducci <marcoguiducci a ymail.com>
> Cc: "gfoss a lists.gfoss.it" <gfoss a lists.gfoss.it>
> Inviato: Martedì 2 Ottobre 2012 14:49
> Oggetto: Re: [Gfoss] GRASS esportare layer 0
> 
> 2012/10/2 Marco Guiducci <marcoguiducci a ymail.com>:
>> 
>> 
>>   Ciao,
>>   abbiamo bisogno di esportare nel formato shapefile tramite v.out.ogr da 
> linea di
>>   comando, il layer 0 di un dato vettoriale, ottenuto tramite v.in.ogr da un
>>   precedente shapefile.
>>   In pratica vogliamo travare i vuoti di uno shapefile tramite GRASS e
>>   riesportarli in shapefile.
> 
> Non conosco la parola "travare" :) Se significa mantenere, sarebbe -c
> di v.out.ogr.
> 
> Quale versione di GRASS usate?
> 

"travare" forse deriva, per contrazione, dal predicato francese travailler, e tutto sommato rende bene l'idea se lo traduciamo in travagliare. E' proprio un travaglio!
Venendo al serio: con v.in di uno shapefile ottengo i layer 0, 1 e 2.
Mentre non ho trovato problemi a riesportare i poligoni del layer 1 in overlay con il 2, mantenendo quindi i codici che mi distinguono i poligoni, non riesco a farlo con i poligoni del layer 0.
Riporto qui:

-------------------------
GRASS 6.4.2 (grosseto)> v.out.ogr -c input=sez331090_poligoni layer=0 dsn="c:\po
l_00.shp"
ATTENZIONE: 699 centroid(s) found, but not requested to be exported. Verify
            'type' parameter.
ATTENZIONE: 702 areas found, but not requested to be exported. Verify
            'type' parameter.
ATTENZIONE: No lines found, but requested to be exported. Will skip this
            geometry type.
Esportazione di 2475 geometrie in corso...
 100%
ATTENZIONE: 1776 elementi trovati senza categoria non considerati
v.out.ogr completo. 0 elementi scritti in '<sez331090_poligoni>
(ESRI_Shapefile).

GRASS 6.4.2 (grosseto)>
--------------

ora si da il caso che i 699 centroidi siano esattamente quelli del layer 1, mentre il layer 0 ha tre record (699+3=702)

Se ometto -c, mi esporta tutto il layer 1.
Il problema è nel parametro type? (come ho dato il comando ha esportato le linee)
Se metto type=area

-----
GRASS 6.4.2 (grosseto)> v.out.ogr -c type=area input=sez331090_poligoni layer=0
dsn="c:\pol_0000.shp"
Exporting 702 areas (may take some time)...
 100%
ATTENZIONE: 702 elementi trovati senza categoria non considerati
v.out.ogr completo. 0 elementi scritti in '<sez331090_poligoni>
(ESRI_Shapefile).
------

Da notare, infine, che il comando di QGis "Salva con nome" eseguito sullo 0_polygon lavora fine ;-)
Sono un po' disorientato. Forse è meglio se vado in Messico a trapanar....

Grazie e ciao
marco g.