Come pansharpen LANDSAT 8 in GRASS?

Avendo usato QGIS per un bel po 'di tempo per la mappatura, sono stato "costretto" a scavare in GRASS poiché il pansharpening non è direttamente eseguibile in QGIS. Dopo aver compreso più o meno la struttura, ho fatto alcune prove con i.fusion.brovey. Sono stato in grado di generare i raster con una risoluzione di 15 m, ma durante il caricamento in QGIS non riesco a ottenere / recuperare il piacevole contrasto cromatico come è presente nei raster originali da 30 m.

In un post precedente qualcuno lo ha confermato e ha proposto di utilizzare un'altra funzione per landsat: 'i.his.rgb' (dopo aver usato i.landsat.toar, che non uso). Inizia a convertire min / max in 1-255 prima di usare i.his.rgb ... I miei raster originali hanno un intervallo compreso tra 0 e 65535. Non ha funzionato per me ...

Ho smesso dopo alcuni giorni di difficoltà a eseguire solo una pansharp, mi sento bloccato ... Perché l'applicazione di i.pansharp non dovrebbe funzionare per Landsat come suggerito da altri? Qualcuno può suggerire un flusso di lavoro più semplice per eseguire il pansharpening delle immagini di Landsat 8? In attesa di qualche consiglio, grazie in anticipo ...

Dai un'occhiata a i.landsat.rgb - Esegue il bilanciamento automatico dei colori per le immagini LANDSAT, probabilmente prima di eseguire la pansharpening. Con i.landsat.toar puoi anche considerare di convertire i numeri digitali dei singoli canali in radiosità o riflettenza nella parte alta dell'atmosfera . Vedi anche http://grasswiki.osgeo.org/wiki/LANDSAT

A proposito: avere un intervallo di 0-65535 per i dati LANDSAT originali è dispari poiché vengono forniti come canali a 8 bit (0-255). Quindi controlla prima i tuoi dati di input.

Aggiornamento: 10/2013

Per una procedura semplificata con l'attuale GRASS GIS 7, vedere la pagina Elaborazione dei dati Landsat 8 in GRASS GIS 7: compositi RGB e nitidezza della panoramica

Grazie Michal e Markus. Finalmente sono stato in grado di eseguire il pansharpening con la funzione indicata brov. Una volta eseguito in GRASS, ho esportato i raster RGB su gdal GTiff con il tipo di dati Uint16. Markus, ho scaricato e utilizzato diverse immagini Landsat 8 (tiff di ciascuna delle 11 bande), ottenendo risultati visivi piacevoli, sia in colore naturale / falso. Per quanto ne so, i tiffs di landat originali arrivano come Uint16 con valori che vanno da 0 a 65535 (o mi sto perdendo qualcosa?) ... Grazie per il link GRASS LANDSAT ...

Il FIHS o la trasformazione di Brovery sono abbastanza semplici se segui le formule di un articolo come questo: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Supponendo che si stia utilizzando GRASS GIS e Pansharpening Landsat 8.

Imposta la tua regione usando il tif di pan 8 - importante per avere una risoluzione di 15m che verrà impostata se lo fai.

Per Brovery senza pesi inseriscilo nel Calcolatore mappe rispettivamente per blu, verde e rosso:

band8 / (band2 + band3 + band4) / 3 * band2

band8 / (band2 + band3 + band4) / 3 * band3

band8 / (band2 + band3 + band4) / 3 * band4

Esegui r.colors per ognuna delle uscite precedenti e imposta l'equalizzazione dell'istogramma e la tabella dei colori su grigio.

Esegui r.composite, seleziona le bande appropriate (rosso = 4, verde = 3, blu = 2)

Guarda l'immagine.

Ovviamente questo processo potrebbe essere fatto con maggiore considerazione per la correzione atmosferica, l'allungamento del colore e i pesi della banda, come discusso nell'articolo precedente.

Puoi anche provare Monteverdi per eseguire il pansharpening. Questo strumento può utilizzare GeoTIFF regolari come input (tuttavia, le bande multispettrali separate devono essere unite in un raster multibanda) e si comporta abbastanza bene.

Dopo aver affrontato questo problema, ho trovato una soluzione semplice che ti consente di mettere a fuoco in QGIS (2.18.21 è la mia versione) usando SCP. L'aggiornamento a SCP che incorpora questo metodo è descritto qui https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Questo è facile da fare e ha funzionato come un fascino! Ho controllato l'output finale e i pixel misurano con una risoluzione di 15 m, abbinando la banda panchromatic 8. Se qualcuno è curioso, utilizza il metodo di trasformazione di Brovey come descritto nel manuale (vedere https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

Benvenuti nel sito Janos. Idealmente le risposte come la tua dovrebbero avere una descrizione del perché sono buoni metodi. ..

Per la pansharpening delle immagini di Landsat 8 il metodo più semplice è il http://www.geosage.com/highview/download.html . Il mio consiglio, provalo. Solo un clic ....

Per pansharpenig ho usato il trasformatore spettrale GUI per le immagini di Landsat 8 http://www.geosage.com/highview/download.html Molto buono. BR Janos