[Gfoss] Riproiezione di un raster

[ludovico frate]

Come hai detto giustamente, "Mettendo da parte i perchè, si è scelto di
lavorare in questo modo" perché lo trovo davvero poco efficace.
Comunque, hai controllato che il resampling method sia lo stesso?


Il giorno mar 17 mar 2020 alle ore 09:52 Massimiliano Moraca <
info a massimilianomoraca.it> ha scritto:

> Buongiorno a tutti, scrivo per esporvi alcuni dubbi che mi sono venuti
> relativamente alla riproiezione di un raster. In ufficio, per alcuni indici
> di vegetazione, abbiamo la necessità di visualizzare in un diagramma a
> barre
> quanta superficie del raster in esame rientra in un certo range di valori.
> Mi spiego meglio: posti tre range x, y e z, vogliamo sapere quanti pixel
> rientrano in ognuno dei range e quale è la superficie totale di ognuno dei
> range. Ad esempio: nel range x rientrano 20 pixel, la dimensione di ogni
> pixel è 10x10 metri, quindi la superficie totale di x è 10x10x20=2000mq
>
> Mettendo da parte i perchè, si è scelto di lavorare in questo modo:
> 1 - vettorializzazione del raster dell'indice di vegetazione
> 2 - calcolo della superficie di ogni poligono del vettore(ogni poligono
> corrisponderà al singolo pixel del raster al punto 1)
> 3 - raggruppamento delle ricorrenze dei pixel secondo i range e sommatoria
> delle aree
>
> Il dato di partenza è in 4326 con dimensione del singolo pixel pari a
> 0.000112148462163201668 x 0.000112148462163201668
>
> Uno sviluppatore si sta occupando di creare il flusso di lavoro in python,
> io con QGIS ho fatto una verifica sui suoi risultati partendo da uno degli
> indici calcolati.
>
> Usando /Polygonize(raster to vector)/ ho ottenuto un vettore con una
> colonna
> DN che contiene il valore di indice per ogni poligono. Ho a questo punto
> aggiunto una colonna che ho chiamato area_mq con cui tramite field
> calculator ho calcolato la superficie di ogni poligono: *area( transform(
> $geometry, 'EPSG:4326', 'EPSG:3857'))*. Il risultato è pari a circa 206 mq
> per ogni pixel.
>
> Lo sviluppatore usando rasterio ha riproiettato il raster in 3857 e si è
> poi
> calcolato l'estensione di ogni pixel che gli risulta pari a circa 175 mq
> con
> pixel di dimensioni 13.26663473960738138 x 13.26663473960738138 m.
>
> Per capire il perchè di queste differenze ho esportato il mio raster in
> 3857(Export -> Save as..) trovandomi con pixel di dimensioni
> 12.48430981132078799 x 16.5115487603316744 e quindi la superficie di ogni
> pixel è circa 206 mq. Mi trovo con il risultato che mi sono calcolato in
> precedenza ma non con il risultato dello sviluppatore.
>
> Secondo test: ho usato /Warp/ per riproiettare il raster e mi trovo con le
> dimensioni in pixel dello sviluppatore.
>
> Ora, posto che ho sempre riproiettato un raster nello stesso modo con cui
> si
> riproietta un vettore non avendo mai riscontrato problemi di
> riposizionamento vorrei capire il perchè di quello che sembra essere un mio
> errore metodologico.
>
> -----
> Consulente GIS,  Formatore, Blogger e Ciclista Urbano
> email: info a massimilianomoraca.it
> cell: 333 5949583 (lun-ven, 9.00-18.00)
> website: massimilianomoraca.it
> --
> Sent from:
> http://gfoss-geographic-free-and-open-source-software-italian-mailing.3056002.n2.nabble.com/
> _______________________________________________
> Gfoss a lists.gfoss.it
> 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.
> 764 iscritti al 23/08/2019



-- 
*Dott. For. Ludovico Frate, Ph.D.*

*Via Kennedy 80, 86170 - Isernia - ITALIA.*
*Via Montalto 17, 86087 - Rionero Sannitico (IS) - ITALIA.*
Studio Tecnico Forestale e GIS *FORGIS* <http://www.forgis.it/>
*Collaboratore presso Studio Associato Ecoview <http://www.ecoview.it/>*
*Collaboratore presso Unimol EnviXlab <http://envixlab.unimol.it/>*

*Docente GIS presso lezionigis.it <https://www.lezionigis.it>*
*Cel: ++39 3333767557*

*P.IVA 00960030948E-mail: *frateludovico a gmail.com*|*
ludovicofrate a hotmail.it
*Research Gate
<https://www.researchgate.net/profile/Ludovico_Frate?ev=hdr_xprf&_sg=lGaI2daIBO6kPmtye069ckFfDExcPoNVKJHrqvQEPQmsmnHolvnRZzPQdcyxs1g7>*
|*Google Scholar
<https://scholar.google.it/citations?user=wGFhBrkAAAAJ&hl=it>*|*LinkedIn
<https://www.linkedin.com/in/ludovico-frate-a387aa57?trk=hp-identity-name>*