Jak pansharpen LANDSAT 8 w GRASS?

Korzystając z QGIS od dłuższego czasu do mapowania, zostałem „zmuszony” do wkopania się w GRASS, ponieważ funkcja ostrzenia pansów nie jest bezpośrednio wykonywalna w QGIS. Po zrozumieniu mniej więcej tej struktury przeprowadziłem kilka prób z i.fusion.brovey. Byłem w stanie wygenerować rastry o rozdzielczości 15 m, ale podczas ładowania do QGIS nie mogę uzyskać / odzyskać ładnego kontrastu kolorów, jaki występuje w oryginalnych rastrach 30 m.

We wcześniejszym poście ktoś to potwierdził i zaproponował użycie innej funkcji dla landsat: „i.his.rgb” (po użyciu i.landsat.toar, której nie używam). Zaczyna konwertować min / max na 1-255 przed użyciem i.his.rgb ... Moje oryginalne rastry mają zakres 0-65535. Nie wyszło mi to ...

Rezygnuję po kilku dniach zmagań, aby po prostu wykonać pansharp, czuję się utknięty ... Dlaczego po prostu zastosowanie i.pansharp nie działa dla Landsata, jak sugerują inni? Czy ktoś może zasugerować prostszy obieg pracy w celu wykonania ostrzenia obrazu Landsat 8? Czekamy na porady, z góry dzięki ...

Spójrz na i.landsat.rgb - Wykonuje automatyczne równoważenie kolorów dla obrazów LANDSAT, prawdopodobnie przed uruchomieniem funkcji ostrzenia. Można również rozważyć konwersję cyfrowej liczby poszczególnych kanałów na najwyższy poziom jasności lub odbicia za pomocą i.landsat.toar . Zobacz także http://grasswiki.osgeo.org/wiki/LANDSAT

BTW: posiadanie zakresu 0-65535 dla oryginalnych danych LANDSAT jest nieparzyste, ponieważ są one dostarczane jako kanały 8-bitowe (0-255). Najpierw sprawdź swoje dane wejściowe.

Aktualizacja: 10/2013

Uproszczona procedura z obecnym GRASS GIS 7 znajduje się na stronie Przetwarzanie danych Landsat 8 w GRASS GIS 7: Kompozyty RGB i wyostrzanie pan

Dzięki Michał i Markus. W końcu udało mi się wykonać panoramowanie ze wskazaną funkcją brov. Po wykonaniu w GRASS wyeksportowałem rastry RGB do gdal GTiff z typem danych Uint16. Markus, pobierałem i korzystam z kilku zdjęć Landsat 8 (tiff każdego z 11 pasm), uzyskując ładne efekty wizualne, zarówno w naturalnym / fałszywym kolorze. O ile mi wiadomo, oryginalne tiffy typu landsat mają postać Uint16 z wartościami od 0 do 65535 (czy czegoś mi brakuje?) ... Dzięki za link GRASS LANDSAT ...

Pan-ostrzenie FIHS lub Brovery Transform jest dość proste, jeśli postępujesz zgodnie z formułami z artykułu takiego jak ten: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Zakładając, że używasz GRASS GIS i Pansharpening Landsat 8.

Ustaw swój region za pomocą panoramy 8 tif - ważne, aby mieć rozdzielczość 15 m, która zostanie ustawiona, jeśli to zrobisz.

W przypadku Brovery bez ciężarów wprowadź to do Kalkulatora mapy odpowiednio dla niebieskiego, zielonego i czerwonego:

pasmo 8 / (pasmo 2 + pasmo 3 + pasmo 4) / 3 * pasmo 2

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

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

Uruchom r. Kolory dla każdego z powyższych wyjść i ustaw korektę histogramu oraz tabelę kolorów na szarą.

Uruchom r.composite, wybierz odpowiednie pasma (czerwony = 4, zielony = 3, niebieski = 2)

Sprawdź obraz.

Oczywiście proces ten można wykonać z większym uwzględnieniem korekcji atmosferycznej, rozciągnięcia koloru i ciężaru pasma, jak omówiono w powyższym artykule.

Możesz także wypróbować Monteverdiego, aby wykonać ostrzenie. To narzędzie może pobierać zwykłe GeoTIFF jako dane wejściowe (jednak oddzielne pasma wielospektralne muszą zostać połączone w jeden wielopasmowy raster) i działa dość dobrze.

Po zmaganiu się z tym problemem znalazłem proste rozwiązanie, które pozwala na panoramowanie w QGIS (2.18.21 to moja wersja) za pomocą SCP. Aktualizacja SCP obejmująca tę metodę została opisana tutaj https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Jest to łatwe do zrobienia i działało jak urok! Sprawdziłem końcową moc wyjściową, a rozdzielczość pikseli wynosiła 15 m - pasując do panchromatycznego pasma 8. Jeśli ktoś jest ciekawy - korzysta z metody Brovey Transform zgodnie z opisem w instrukcji (patrz https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

Witamy na stronie Janos. Idealne odpowiedzi, takie jak Twoja, powinny zawierać opis, dlaczego są to dobre metody. ..

W przypadku podglądu zdjęć Landsat 8 łatwiejszą metodą jest http://www.geosage.com/highview/download.html . Moja rada, spróbuj. Tylko jedno kliknięcie ....

Do pansharpenig użyłem transformatora GUI do zdjęć Landsat 8 http://www.geosage.com/highview/download.html Bardzo dobrze. BR Janos