Home GeoForum Forums Geologia & Informatica GIS, CAD e grafica Riclassificare raster sulla base di poligoni

Salve a tutti

Sto lavorando su un'area piuttosto ampia di cui possiedo il particellare ed una mappa NDVI da me calcolata;

quello che ora dovrei fare e riclassificare in 3 classi ciascun poligono separatamente usando un equal interval (prndendo MIN e MAX dio ciascun poligono e dividendo per 3);

l'unico modo cervellotico trovato fino ad ora è quello di creare ogni volta una maschera differente (una per ciascun campo) con cui clippare l'NDVI prima di riclassificarlo;
In questa volta ottengo uno strato per ogni poligono, con la giusta riclassificazione..

Conoscete per caso un modo automatico per poter riclassificare ciascuna particella automaticamente in 3 classi senza ripetere la medesima procedura particella per particella??

Grazie mille

Roberto

PS conosco e dispongo di QGis ma anche di ArcMap nel caso fosse necessario..

Ehm, sarebbe meglio che ti spiegassi con un esempio grafico, schermate, ecc.

Provo a chiarire meglio con l'ausilio di questa immagine..



Quello che vedete in toni di grigio è un NDVI, i poligoni in rosso sono il particellare;

Quello che dovrei fare è prendere ogni singolo poligono e riclassificarlo mediante "equal interval" in tre classi (leggendo il valore max e min di ndvi di quel poligono, dividendolo per tre e riclassificando i punti di quel poligono in tre classi).

In generale con i vari softwares riesco solo a calcolare tre classi uniche ed uguali per tutti i poligoni del particellare; i vari metodi che ho trovato leggono max e min dell'intera immagine, dividono per tre ed applicano quegli intervalli a tutti i pixels dei vari poligoni, non ragionano in relativo come servirebbe a me..

come dicevo l'unico modo che ho trovato per farlo e mascherare di volta in volta ogni singolo poligono e poi riclassificarlo; a quel punto i programmi leggono comunque max e min dell'intera immagine, ma corrispondendo questa con un solo campo mi calcolano gli intervalli giusti..

spero di essermi spiegato meglio

Grazie anticipatamente a tutti voi

Spiegaci meglio il suddetto concetto con un esempio numerico, grazie.

Provo a chiarire..

all'interno di un poligono (vettoriale) ricadono per es 90 pixels della mappa NDVI ciascuno con un proprio valore di DN;

questi 90 pixel devono essere riclassificati in 3 classi usando il criterio equal interval;

per farlo si legge il valore del pixel più alto di NDVI ricadente nel poligono (es 0.8), si legge quello del pixel più basso (es 0,2), si calcola un range (0.8-0.2=0,6) e si divide il range per 3 (0,6/3=0,2) trovando l'ampiezza delle classi.

a questo punto si dividono i pixel di quel determinato poligono in tre classi in base ai valori calcolati:
classe 1 da 0.2 a 0.4 (0,2+ampiezza classe calolata di 0.2)
classe 2 da 0.41 a 0.6 (0.4 +ampiezza classe calolata di 0.2)
classe 3 da 0.61 a 0.8 (0.6 +ampiezza classe calolata di 0.2)

più lungo a dirsi che a farsi

spero di aver chiarito..

a me serve farlo poligono per poligono e non applicando un unico max e min calcolato sull'intera immagine ed applicato a tutti i poligoni

Grazie per la pazienza

Una cosa mi manca: con quale criterio vorresti dividere in 3 parti il mini-raster sotteso ad ogni poligono, verticalmente, orizzontalmente, come..?

devo riclassificare in base al valore di ogni pixel che ricade in ogni poligono, assegnando ad ogni singolo pixel ad una delle tre classi calcolate come sopra.. non si tratta di una divisione geometrica ma di una riclassificazione di ciascun pixel che ricade in ogni poligono..

alla fine dovrei ottenere un risultato tipo questo:



(qui avevo riclassificato in 6 classi invece che 3 ma usandop lo stesso criterio..)

Spiacente, almeno in Qgis non riesco a trovare un algoritmo e/o plugin che "legga" automaticamente i valori min/max di un raster, per poi compiere l'operazione:

min ÷ (min+(max-min)/3) => classe A
(min+(max-min)/3) ÷ 2*(min+(max-min)/3) => classe B
2*(min+(max-min)/3) ÷ 3*(min+(max-min)/3) => classe C

Se trovi una soluzione, sempre in Qgis, hai una cena pagata...