calcolo aree poligoni "minori" incuse in poligoni "maggiori"

Ciao a tutti/e

Ho un problema che in linea teorica non mi sembra
particolarmente complesso ma per il quale sto trovando dei problemi...


ho due file vettoriali:

(1) fogli_catastali
(2)
corpi_idrici_sotterranei

Fra i vari attributi, il file
"fogli_catastali" contiene una colonna con un codice univoco per
ciascun foglio catastale (Com_fogl) ed una colonna con l'area di
ciascun foglio catastale (Area).
Mentre il file
"corpi_idrici_sotterranei" contiene una colonna in cui sono indicati i
nomi dei corpi idrici (COD_C_ID)

Quello che vorrei fare è sapere quali
fogli catastali entrano in ciascun corpo idrico ed in che quantità (per
esempio nel corpo idrico X entra tutto il foglio catastale "a" e tot m2
del foglio catastale "b").

Quello che pensavo di fare è:

(1) un
overlay fra le due mappe: v.overlay (a=fogli_catastali,
b=corpi_idrici_sotterranei, operatore=and) [1]
(2) aggiungere una
colonna al file overlay: v.db.addcol [2]
(3) calcolare l'area nella
nuova colonna: v.to.db [3]

Qualcosa però non ha funzionato in quanto
c'è una grossa discrepanza fra le aree calcolate e le aree dei fogli
catastali...

dove posso aver sbagliato??
grazie mille!!
marco

[1] v.
overlay ainput=fogli_catastali@Umbria atype=area alayer=1
binput=corpi_idrici_sotterranei@Umbria btype=area blayer=1
output=overlay_cat_loc operator=and olayer=1,0,0
[2] v.db.addcol
map=overlay_cat_loc@Umbria layer=1 'columns=area_new double precision'

[3] v.to.db map=overlay_cat_loc@Umbria layer=1 qlayer=1 option=area
units=meters columns=area_new



Con Tiscali Photo puoi creare splendidi calendari con le tue foto direttamente via WEB.Stampali o condividili con chi vuoi tu!http://photo.tiscali.it
_______________________________________________
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.
630 iscritti al 1.12.2012

[giohappy]

Domanda banale, in che sistema di riferimento stai lavorando? Le aree precalcolate, presenti nella relativa colonna, sono state calcolate nello stesso SR?

Che entità hanno le discrepanze?

Ultima domanda, hai provato a calcolare l'area direttamente sul layer dei catastali, per confrontarle con quelle precalcolate?

giovanni

Il giorno 17 dicembre 2012 10:59, [hidden email] <[hidden email]> ha scritto:

Ciao a tutti/e

Ho un problema che in linea teorica non mi sembra
particolarmente complesso ma per il quale sto trovando dei problemi...


ho due file vettoriali:

(1) fogli_catastali
(2)
corpi_idrici_sotterranei

Fra i vari attributi, il file
"fogli_catastali" contiene una colonna con un codice univoco per
ciascun foglio catastale (Com_fogl) ed una colonna con l'area di
ciascun foglio catastale (Area).
Mentre il file
"corpi_idrici_sotterranei" contiene una colonna in cui sono indicati i
nomi dei corpi idrici (COD_C_ID)

Quello che vorrei fare è sapere quali
fogli catastali entrano in ciascun corpo idrico ed in che quantità (per
esempio nel corpo idrico X entra tutto il foglio catastale "a" e tot m2
del foglio catastale "b").

Quello che pensavo di fare è:

(1) un
overlay fra le due mappe: v.overlay (a=fogli_catastali,
b=corpi_idrici_sotterranei, operatore=and) [1]
(2) aggiungere una
colonna al file overlay: v.db.addcol [2]
(3) calcolare l'area nella
nuova colonna: v.to.db [3]

Qualcosa però non ha funzionato in quanto
c'è una grossa discrepanza fra le aree calcolate e le aree dei fogli
catastali...

dove posso aver sbagliato??
grazie mille!!
marco

[1] v.
overlay ainput=fogli_catastali@Umbria atype=area alayer=1
binput=corpi_idrici_sotterranei@Umbria btype=area blayer=1
output=overlay_cat_loc operator=and olayer=1,0,0
[2] v.db.addcol
map=overlay_cat_loc@Umbria layer=1 'columns=area_new double precision'

[3] v.to.db map=overlay_cat_loc@Umbria layer=1 qlayer=1 option=area
units=meters columns=area_new



Con Tiscali Photo puoi creare splendidi calendari con le tue foto direttamente via WEB.Stampali o condividili con chi vuoi tu!http://photo.tiscali.it
_______________________________________________
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.
630 iscritti al 1.12.2012

--
--
Giovanni Allegri
website: http://giovanniallegri.it

GEO+ geomatica in Italia http://bit.ly/GEOplus

_______________________________________________
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.
630 iscritti al 1.12.2012