Buongiorno a tutti,
in alternativa dell'uso di "Raster volume" stavo provando a fare il calcolo
di volumi utilizzando il plugin "statistiche zonali".
L'idea è quella di calcolare, per una serie di aree racchiuse neri
rispettivi poligoni, le statistiche zonali sul DEM ed utilizzare
SUM*(superficie del pixel) per ottenere il volume.
Vi condivido alcuni problemi:
- la formula funziona solo se tutti i dati di elevazione sono positivi. Per
capirci, immaginate di lavorare in una zona di poco al di sotto della quota
"zero", che ci siano alcuni punti del DEM postivi ed altri negativi e che il
volume vada calcolato a partire alla quota -2.00. Come procedereste senza
usare "Raster Volume" e quindi senza "tagliare con maschere" l'intero DEM?
- ho notato che "statistiche zonali" escluda dal calcolo non solo i pixel a
contatto con poligono (quelli perimetrali all'area) ma che l'esclusione vada
addirittura più all'interno come se ci fosse una specie di offset. Vi
risulta?
Vi allego un DEM esempio di soli 4 pixel 1m x 1m, ed il relativo vettore
agganciato al perimetro dell'area
https://www.dropbox.com/s/pmfw11yg5z2ogqo/Area.sqlite?dl=0
<https://www.dropbox.com/s/pmfw11yg5z2ogqo/Area.sqlite?dl=0>
https://www.dropbox.com/s/3w429j7s6guzjeb/DEM_points.tif?dl=0
<https://www.dropbox.com/s/3w429j7s6guzjeb/DEM_points.tif?dl=0>
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
Luca Tafuro wrote
Buongiorno a tutti,
in alternativa dell'uso di "Raster volume" stavo provando a fare il
calcolo
di volumi utilizzando il plugin "statistiche zonali".
L'idea è quella di calcolare, per una serie di aree racchiuse neri
rispettivi poligoni, le statistiche zonali sul DEM ed utilizzare
SUM*(superficie del pixel) per ottenere il volume.
Vi condivido alcuni problemi:
- la formula funziona solo se tutti i dati di elevazione sono positivi.
Per
capirci, immaginate di lavorare in una zona di poco al di sotto della
quota
"zero", che ci siano alcuni punti del DEM postivi ed altri negativi e che
il
volume vada calcolato a partire alla quota -2.00. Come procedereste senza
usare "Raster Volume" e quindi senza "tagliare con maschere" l'intero DEM?
- ho notato che "statistiche zonali" escluda dal calcolo non solo i pixel
a
contatto con poligono (quelli perimetrali all'area) ma che l'esclusione
vada
addirittura più all'interno come se ci fosse una specie di offset. Vi
risulta?
Vi allego un DEM esempio di soli 4 pixel 1m x 1m, ed il relativo vettore
agganciato al perimetro dell'area
Ciao Luca;
ho fatto una rapida prova:
usando il raster calculator, applico una sorta di ABS (valore assoluto) a
tutti i valori dei pixel con la seguente formula:
( "DEM_points@1"<0 )*-"DEM_points@1" +( "DEM_points@1">0 ) *"DEM_points@1"
successivamente utilizzo 'zonal statistics' ottenendo, come sum, 6, credo
sia iòl risultato corretto;
per quanto riguarda il secondo problema che segnali NON si è verificato in
questo caso.
saluti
-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
Ciao Salvatore.
Otti, appena posso provo.
In merito al problema, hai usato i due file che ho condiviso o con altro
dataset?
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
Ciao Salvatore,
grazie per la dritta con Raster Calculator, ottima idea.
Invece, in merito ai problemi con "statistiche zonali" hai provato con il
mio dataset o con file tuoi?
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
Luca Tafuro wrote
Ciao Salvatore,
grazie per la dritta con Raster Calculator, ottima idea.
Invece, in merito ai problemi con "statistiche zonali" hai provato con il
mio dataset o con file tuoi?
Con i tuoi dati condivisi.
-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html