[Talk-it] Confini comunali per ritaglio shape

[marco bra]

Per ottenere gli shp dagli osm dei vari boundary derivanti dagli admin_levels:

- Per un singolo comune, Arenzano in questo caso:
wget -O b_arenzano.osm
'http://overpass.osm.rambler.ru/cgi/interpreter?data=area[name="Arenzano"][admin_level=8];relation(area)["admin_level"="8"];(._;>;);out
meta;'
ogr2ogr -overwrite -sql "select name as nomecomu from multipolygons"
b_arenzano.shp b_arenzano.osm multipolygons
ogr2ogr -overwrite boundary_arenzano.shp arenzano.osm multipolygons

- Per il boundary di una intera regione, Liguria in questo caso:
wget -O b_liguria.osm
'http://overpass.osm.rambler.ru/cgi/interpreter?data=area[name="Liguria"][admin_level=4];relation(area)["admin_level"="4"];(._;>;);out
meta;'; ogr2ogr -overwrite
ogr2ogr -overwrite  b_liguria.shp b_liguria.osm multipolygons


- Per boundary di tutti i comuni liguri
wget -O b_com_liguria.osm
'http://overpass.osm.rambler.ru/cgi/interpreter?data=area[name="Liguria"][admin_level=4];relation(area)["admin_level"="8"];(._;>;);out
meta;'
ogr2ogr -overwrite -sql "select name as nomecomu from multipolygons"
b_com_liguria.shp b_com_liguria.osm

Con qgis:

-  verificare eventuali sovrapposizioni fra i layer edifici ed i
confini degli shape "di taglio".

Poi supponendo di avere tutti gli edifici della Liguria e li volessimo
dividere per comune
per prima cosa assegniamo una nuova colonna "nomecomu" ai buildings
che contenga il nome del comune di appartenenza, tramite:

- Vettore → Strumenti di gestione dati → Unisci per posizione unendo i
buildings.shp  con b_com_liguria.shp

creeremo in tal modo un nuovo shape ad esempio buildings_liguria_per_comune.shp
dove all'interno troveremo tutti gli attributi dei buildings.shp ed in
più  la colonna "nomecomu" del comune di appartenenza.

- poi tramite  Vettore → Strumenti di gestione dati → Dividi vettore
potremo spezzare il file shape buildings_liguria_per_comune.shp usando
come campo ID unico il campo "nomecomu"
verranno creati i singoli files shape degli edifici per ogni comune...
nella cartella desiderata.

Ciao
Marco /mcheck/mcheckimport

Il 02 marzo 2014 18:58, Maurizio Napolitano <napoogle a gmail.com> ha scritto:
> secondo me fai prima a prendere gli shape file dei confini comunali dell'ISTAT.
> Cmq, ,interrogando per area di relation, puoi fare le query con overpass-api per
> vedere quali sono gli edifici contenuti in un comune
>
> 2014-03-02 17:09 GMT+01:00 Leonardo Frassetto <kinetocore86 a gmail.com>:
>> Consigli su come impostare la query su josm per scaricare i confini?
>>
>> Il 02/mar/2014 17:02 "sabas88" <sabas88 a gmail.com> ha scritto:
>>
>>>
>>>
>>>
>>> Il giorno 02 marzo 2014 13:58, Leonardo Frassetto <kinetocore86 a gmail.com>
>>> ha scritto:
>>>>
>>>> Domanda per gli smanettoni di GIS e dati OSM: è possibile ricavare i
>>>> confini comunali da osm, esportarli in qualche  formato compatibile con qgis
>>>> e usarli come area di selezione per selezionare gli edifici di un file shp?
>>>> Vorrei tentare di separare in piccole parti un file molto molto grosso .
>>>> Grazie!
>>>>
>>>>
>>>
>>> Scarica i confini comunali e salvali come geojson, qgis dovrebbe leggerlo.
>>> Al limite i confini copiali su un nuovo livello e falli diventare un
>>> poligono unico...
>>>
>>> Ciao,
>>> Stefano
>>>>
>>>> _______________________________________________
>>>> Talk-it mailing list
>>>> Talk-it a openstreetmap.org
>>>> https://lists.openstreetmap.org/listinfo/talk-it
>>>>
>>>
>>>
>>> _______________________________________________
>>> Talk-it mailing list
>>> Talk-it a openstreetmap.org
>>> https://lists.openstreetmap.org/listinfo/talk-it
>>>
>>
>> _______________________________________________
>> Talk-it mailing list
>> Talk-it a openstreetmap.org
>> https://lists.openstreetmap.org/listinfo/talk-it
>>
>
>
>
> --
> Maurizio "Napo" Napolitano
> http://de.straba.us
>
> _______________________________________________
> Talk-it mailing list
> Talk-it a openstreetmap.org
> https://lists.openstreetmap.org/listinfo/talk-it



-- 
Linux Infinite Freedom

I'm writing from this place:
http://www.openstreetmap.org/?lat=44.39945&lon=8.6798&zoom=15&layers=M