[QGIS-it-user] Statistiche dati raster su poligoni di un vettoriale

Buongiorno, avrei bisogno di elaborare dei calcoli su un modello digitale del
terreno in modo da avere, per ogni poligono di un layer vettoriale, alcuni
parametri delle quote quali, almeno: max; min; media.
Mi sapete indicare quale sia lo strumento adatto per favore?
Grazie

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.

Ciao Lorenzo,
il plugin già preinstallato Statistiche Zonali (lo trovi anche in
processing) dovrebbe fare al caso tuo.

Saluti

Matteo

ghtmtt wrote

il plugin già preinstallato Statistiche Zonali (lo trovi anche in
processing) dovrebbe fare al caso tuo.
Saluti
Matteo

Grazie mille Matteo, ho usato quello del processing (Qgis 2.8.2 su win) che
da più valori, anche massimi, minimi range ecc.
Solo che ho un problema (che ho risolto a mano) con le virgole, facendo il
copia/incolla dalla tabella excel riconosce il punto di Qgis come migliaia
ed i valori, che evidentemente sono preimpostati con un numero alto di
decimali, risultano sballati (dell'ordine di grandezza delle centinaia di
milioni), sapresti per caso come risolvere? grazie ancora

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257668.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.

Ciao Lorenzo,

mmm allora il fatto dei punti e delle virgole è sempre un grosso
problema (almeno per noi italiani e che usiamo le cose al contrario..)

secondo me, da come descrivi il problema, dovresti cambiare le
impostazioni in excel (non so da dove).. probabilmente c'è qualche
conflitto fra QGIS (che vede i punti come separatore di decimali) e excel..
cambia le impostazioni della lingua in excel e dovresti trovare qualche
opzione..

Saluti

Matteo

Ciao Lorenzo,

ma perché non usi OpenOffice o LibreOffice? In Excel 2007 e 2010: File - Opzioni - Impostazioni avanzate - Utilizza separatori di sistema - Separatore decimale / Separatore delle migliaia. Spero esserti stata d’aiuto.

Saluti,
Gabriela


Da: matteo matteo.ghetta@gmail.com
A: qgis-it-user@lists.osgeo.org
Inviato: Lunedì 21 Marzo 2016 17:58
Oggetto: Re: [QGIS-it-user] Statistiche dati raster su poligoni di un vettoriale

Ciao Lorenzo,

mmm allora il fatto dei punti e delle virgole è sempre un grosso
problema (almeno per noi italiani e che usiamo le cose al contrario…)

secondo me, da come descrivi il problema, dovresti cambiare le
impostazioni in excel (non so da dove)… probabilmente c’è qualche
conflitto fra QGIS (che vede i punti come separatore di decimali) e excel…
cambia le impostazioni della lingua in excel e dovresti trovare qualche
opzione…

Saluti

Matteo


QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-it-user

Grazie mille per gli aiuti tempestivi vediamo un po' di ricapitolare:
innanzitutto dopo il primo colpo andato a buon fine ho ripetuto l'analisi
con lo stesso strumento di geoprocessing (grid statistic for polygon), ma mi
restituiva un unico record (creando un file temporaneo) con sette colone con
numeri che non saprei dove sia andato a prendere, quindi ho installato essen
e con il plugin nativo "statstiche zonali" me la sono cavata, tra l'altro è
comodissimo l'aggiungere le statistiche alla tabella del vettoriale
aggiungendo un prefisso a seconda del raster utilizzato.

Veniamo al problema delle virgole:

Gabriela Osaci Costache wrote

Ciao Lorenzo,
ma perché non usi OpenOffice o LibreOffice? In Excel 2007 e 2010: File -
Opzioni - Impostazioni avanzate - Utilizza separatori di sistema -
Separatore decimale / Separatore delle migliaia. Spero esserti stata
d'aiuto.
Saluti,Gabriela

Devo prendere il treno e non ho tempo di scaricare openoffice, ma proverò ed
aggiorno il topic

Gabriela Osaci Costache wrote

      Da: matteo <

matteo.ghetta@

>
A:

qgis-it-user@.osgeo

Inviato: Lunedì 21 Marzo 2016 17:58
Oggetto: Re: [QGIS-it-user] Statistiche dati raster su poligoni di un
vettoriale
   
Ciao Lorenzo,

mmm allora il fatto dei punti e delle virgole è sempre un grosso
problema (almeno per noi italiani e che usiamo le cose al contrario..)

secondo me, da come descrivi il problema, dovresti cambiare le
impostazioni in excel (non so da dove).. probabilmente c'è qualche
conflitto fra QGIS (che vede i punti come separatore di decimali) e
excel..
cambia le impostazioni della lingua in excel e dovresti trovare qualche
opzione..

Saluti

Matteo

Le impostazioni di office per i separatori dei decimali e delle migliaia
(opzioni di office-->impostazioni avanzate) le ho provate tutte (lasciando
la spunta, togliendola, mettendo il punto come separatore dei decimali
ecc.), ma non me la sono cavata.
Devo dire inoltre che i valori generati per le statistiche sulle pendenze mi
lasciano alquantro allibito, per ora me le tengo così, ma spero di riuscire
a fare qualche controprova.

Grazie ancora a tutti e buona serata.

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257697.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.

Buongiorno,

riguardo la virgola non credo bastino le impostazioni di excel ma devi modificare le impostazioni di sistema di windows. su W7

pannello di controllo->Paese e lingua->impostazioni aggiuntive.

Personalmente di sconsiglio di farlo perchè questo influenzerebbe tutto (chi lavora con geomedia è costretto a farlo ad esempio).

Per ovviare al problema puoi fare questo: aprire il file dbf con calc (ti sconsiglio excel anche io) e fai un trova sostituisci dei punti e delle virgole.

Un altra cosa che puoi fare è salvare il dato come csv, editarlo con un editor di testo (notepad++) e rifare il join con il vettoriale originale utilizzando un campo univoco.
Io ho provato a fare la tua operazione su un vettoriale che sta su db postgresql+postgis e il problema delle virgole non si pone. Forse è lo stesso se usi un db spatialite ma non l’ho testato.

Riguardo ai valori della pendenza che ottieni i problemi potrebbero essere diversi:

-la pendenza deve essere espressa in percentuale non in gradi (fare la media fra 360 e 1 grado, che è quasi lo stesso, ti comporta che la pendenza ti vienei di 180,5)

-il dem da cui parti ha degli errori (tutti i dem hanno errori) Per eliminare questo problema l’unica cosa che puoi fare è un analisi degli outleir ed escluderli ma qui si va a toccare una geostatistica un po’ più complessa.

-la risoluzione del raster non è sufficiente per il lavoro che vuoi fare. Il raster è fatto di quadrati e solo se il centroide di questi quadrati cade dentro il tuo poligono viene preso in considerazione. Ti accorgi di questo guardando l’area che calcola sul raster partendo dal vettoriale che non è mai identica all’area calcolata solo sul vettoriale. Se il tuo problema è questo l’unica soluzione è trovare un dem a risoluzione più alta (cosa che mi rendo conto non è proprio banale). Se il poligono più piccolo del tuo vettoriale ha una superficie inferiore a quella del pixel del raster allora non puoi fare il lavoro che stai facendo perchè quel poligono potrebbe benissimo avere un valore di pendenza nullo. Stessa cosa se i tuoi poligoni sono ad esempio lunghi e stretti, se cadono a cavallo fra due pixel (senza perciò passare sul centroide) il valore risulta nullo.

Spero di essere stato d’aiuto.

Saluti

Sergio

···

Il giorno 21 marzo 2016 17:59, Lorenzo Fusco <lorenzofusco@gmail.com> ha scritto:

Grazie mille per gli aiuti tempestivi vediamo un po’ di ricapitolare:
innanzitutto dopo il primo colpo andato a buon fine ho ripetuto l’analisi
con lo stesso strumento di geoprocessing (grid statistic for polygon), ma mi
restituiva un unico record (creando un file temporaneo) con sette colone con
numeri che non saprei dove sia andato a prendere, quindi ho installato essen
e con il plugin nativo “statstiche zonali” me la sono cavata, tra l’altro è
comodissimo l’aggiungere le statistiche alla tabella del vettoriale
aggiungendo un prefisso a seconda del raster utilizzato.

Veniamo al problema delle virgole:

Gabriela Osaci Costache wrote

Ciao Lorenzo,
ma perché non usi OpenOffice o LibreOffice? In Excel 2007 e 2010: File -
Opzioni - Impostazioni avanzate - Utilizza separatori di sistema -
Separatore decimale / Separatore delle migliaia. Spero esserti stata
d’aiuto.
Saluti,Gabriela

Devo prendere il treno e non ho tempo di scaricare openoffice, ma proverò ed
aggiorno il topic

Gabriela Osaci Costache wrote

Da: matteo <

matteo.ghetta@

>
A:

qgis-it-user@.osgeo

Inviato: Lunedì 21 Marzo 2016 17:58
Oggetto: Re: [QGIS-it-user] Statistiche dati raster su poligoni di un
vettoriale

Ciao Lorenzo,

mmm allora il fatto dei punti e delle virgole è sempre un grosso
problema (almeno per noi italiani e che usiamo le cose al contrario…)

secondo me, da come descrivi il problema, dovresti cambiare le
impostazioni in excel (non so da dove)… probabilmente c’è qualche
conflitto fra QGIS (che vede i punti come separatore di decimali) e
excel…
cambia le impostazioni della lingua in excel e dovresti trovare qualche
opzione…

Saluti

Matteo

Le impostazioni di office per i separatori dei decimali e delle migliaia
(opzioni di office–>impostazioni avanzate) le ho provate tutte (lasciando
la spunta, togliendola, mettendo il punto come separatore dei decimali
ecc.), ma non me la sono cavata.
Devo dire inoltre che i valori generati per le statistiche sulle pendenze mi
lasciano alquantro allibito, per ora me le tengo così, ma spero di riuscire
a fare qualche controprova.

Grazie ancora a tutti e buona serata.


View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257697.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.


QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-it-user

Sergio Gollino

Questo messaggio, comprensivo di eventuali allegati, è ad uso esclusivo del destinatario e potrebbe contenere informazioni riservate; se è stato recapitato per errore ci scusiamo per l’accaduto e Vi invitiamo cortesemente a darcene notizia provvedendo alla sua distruzione.

Vi ricordiamo che la diffusione, l’utilizzo e/o la conservazione dei dati ricevuti per errore costituiscono violazione alle disposizioni del D.Lgs. n. 196/2003 “Codice in materia di protezione dei dati personali”.

Sergio Gollino wrote

[...]
Il raster è fatto di quadrati e solo se il centroide di questi quadrati
cade dentro il tuo poligono viene preso in considerazione.

ma se usi la funzione di saga shapes-grid-->add grid values to point oppure
add grid values to shapes, puoi utilizzare una funzione di interpolazione
(bilineare in genere) ed il valore campionato non è quello del centroide ma,
appunto, un'interpolazione della superficie.

s.

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257752.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.

Il valore è sempre quello del centroide. L’interpolazione non fa altro che aumentare il numero di pixel creando dei pixel più piccoli (e di conseguenza la possibilità di cadere fuori dal centroide è minore) il cui valore dipende dai pixel vicini (nel caso del bilineare le 4 celle che toccano il pixel).
Il problema è sempre l’analisi che vogliamo fare.
L’interpolazione è un operazione che falsa il dato e per determinate analisi questo non è accettabile. è un po’ come prendere una foto a 1 megapixel, ingrandirla a 10megapixel e pretendere di vedere i dettagli della foto come se fosse un originale da 10megapixel. L’informazione che non c’è non può sempre essere dedotta e solo perchè i GIS ci permettono di fare un operazione non significa che si debba farla.

sergio

···

Il giorno 22 marzo 2016 08:57, skampus <skampus@gmail.com> ha scritto:

Sergio Gollino wrote

[…]
Il raster è fatto di quadrati e solo se il centroide di questi quadrati
cade dentro il tuo poligono viene preso in considerazione.

ma se usi la funzione di saga shapes-grid–>add grid values to point oppure
add grid values to shapes, puoi utilizzare una funzione di interpolazione
(bilineare in genere) ed il valore campionato non è quello del centroide ma,
appunto, un’interpolazione della superficie.

s.


View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257752.html

Sent from the QGIS Italian User mailing list archive at Nabble.com.


QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-it-user

Sergio Gollino

Questo messaggio, comprensivo di eventuali allegati, è ad uso esclusivo del destinatario e potrebbe contenere informazioni riservate; se è stato recapitato per errore ci scusiamo per l’accaduto e Vi invitiamo cortesemente a darcene notizia provvedendo alla sua distruzione.

Vi ricordiamo che la diffusione, l’utilizzo e/o la conservazione dei dati ricevuti per errore costituiscono violazione alle disposizioni del D.Lgs. n. 196/2003 “Codice in materia di protezione dei dati personali”.

assolutamente corretto e d'accordo, ma anche il dem (sia dtm sia dsm) è
frutto di un'interpolazione che regolarizza una nuvola di punti in una
griglia avente celle di tot metri quadri.
dunque è così improprio stimare "che cosa succede" tra il centroide di una
cella e quello della cella adiacente?
ovviamente non parlo di celle 1 x 1 km2...

dobbiamo considerare una superficie grid come qbert?
:slight_smile:

http://www.pidownload.it/punto/trovafiles/pcbert.png

s.

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257755.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.

Si ma è un nuvola di punti raccolta con una densità tale da permettere di ottenere una griglia con un determinato passo. Se io ho un dem di pixel 20, 40 o 75 (vedi portale nazionale) e ho bisogno di conoscere l’andamento della pendenza sulla zona del carso triestino/isontino ho poco da interpolare perchè li è pieno di doline con un interpolazione otterrei che quella zona è piatta o quasi (cosa molto lontana dalla verità).

In questo caso cosa succedere tra un centroide e l’altro mi interessa eccome e dovrò andare a cercare un dato diverso (lidar ad esempio).

Io non conosco l’analisi che vuole fare Lorenzo e magari lui di cosa succedere fra i due pixel se ne frega. Il mio intento è solo quello di far notare una possibile criticità dell’analisi, non posso dire se questo è il suo caso o meno.

Sergio

···

Il giorno 22 marzo 2016 09:25, skampus <skampus@gmail.com> ha scritto:

assolutamente corretto e d’accordo, ma anche il dem (sia dtm sia dsm) è
frutto di un’interpolazione che regolarizza una nuvola di punti in una
griglia avente celle di tot metri quadri.
dunque è così improprio stimare “che cosa succede” tra il centroide di una
cella e quello della cella adiacente?
ovviamente non parlo di celle 1 x 1 km2…

dobbiamo considerare una superficie grid come qbert?
:slight_smile:

http://www.pidownload.it/punto/trovafiles/pcbert.png

s.


View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257755.html

Sent from the QGIS Italian User mailing list archive at Nabble.com.


QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-it-user

Sergio Gollino

Questo messaggio, comprensivo di eventuali allegati, è ad uso esclusivo del destinatario e potrebbe contenere informazioni riservate; se è stato recapitato per errore ci scusiamo per l’accaduto e Vi invitiamo cortesemente a darcene notizia provvedendo alla sua distruzione.

Vi ricordiamo che la diffusione, l’utilizzo e/o la conservazione dei dati ricevuti per errore costituiscono violazione alle disposizioni del D.Lgs. n. 196/2003 “Codice in materia di protezione dei dati personali”.

scusami, ho colto solo adesso il senso del tuo discorso.
concordo al 100%

s.

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257757.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.

Non c’è nulla di cui scursarsi, anzi grazie per avermi aiutato a chiarire il concetto!

Ciao!

Sergio

···

Il giorno 22 marzo 2016 09:41, skampus <skampus@gmail.com> ha scritto:

scusami, ho colto solo adesso il senso del tuo discorso.
concordo al 100%

s.


View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257757.html

Sent from the QGIS Italian User mailing list archive at Nabble.com.


QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-it-user

Sergio Gollino

Questo messaggio, comprensivo di eventuali allegati, è ad uso esclusivo del destinatario e potrebbe contenere informazioni riservate; se è stato recapitato per errore ci scusiamo per l’accaduto e Vi invitiamo cortesemente a darcene notizia provvedendo alla sua distruzione.

Vi ricordiamo che la diffusione, l’utilizzo e/o la conservazione dei dati ricevuti per errore costituiscono violazione alle disposizioni del D.Lgs. n. 196/2003 “Codice in materia di protezione dei dati personali”.

Buongiorno ed ancora grazie per il supporto, ho risolto (evito di quotare ma
riassumo io i concetti, magari chi si trova a leggere fa prima)

Il problema delle virgole/punti si risolve impostando in calc di open
office, la lingua inglese (USA) che poi voglio sincerarmi con una ricerca
che non siano loro a fare le cose al contrario (vedi guida a sn, once,
galloni, pinte, pollici ecc, il Sistema Internazionale non ne adotta nemmeno
una).

Per quanto riguarda la discussione sul DEM e le interpolazioni in generale
posso dire che da un lato la zona di studio non ha queste asperità eccessive
per cui il dem a 20m (scaricato dal pcn quando non c'erano limitazioni
spaziali) va più che bene per il sito, oltre che per lo scopo dell'analisi
che si prefigge di sintetizzare i valori medi, massimi e minimi di quota,
pendenza ed esposizione prevalente di alcune particelle forestali nelle
schede descrittive, non devo costruire alcun ponte. Inoltre il dem ha pixel
lato 20 per cui 400 mq di superficie mentre le superfici dei poligoni sono
nell'ordine della decina di ettari, per cui va bene così, resta da sistemare
la questione sulle pendenze su cui non sono proprio convinto che sia meglio
quella % rispetto ai gradi visto che (se non ricordo male dall'ultima volta
che ho indagato sul concetto) l'andamento della curva della pendenza
espressa in % è di tipo logaritmico per cui all'aumentare del valore
percentuale aumenta moltissimo il valore pendenza che tende ad infinito su
100%. Farò delle prove se trovo il tempo

Grazie ancora e buona giornata a tutti.

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Statistiche-dati-raster-su-poligoni-di-un-vettoriale-tp5257664p5257775.html
Sent from the QGIS Italian User mailing list archive at Nabble.com.