Wie landeSAT 8 in GRASS pansharpen?

9

Nachdem ich QGIS einige Zeit für das Mapping verwendet hatte, war ich gezwungen, mich mit GRASS zu befassen, da das Pansharpening in QGIS nicht direkt ausführbar ist. Nachdem ich die Struktur mehr oder weniger verstanden hatte, führte ich einige Versuche mit i.fusion.brovey durch. Ich konnte die Raster mit einer Auflösung von 15 m erzeugen, aber beim Laden in QGIS kann ich den schönen Farbkontrast, wie er in den ursprünglichen 30 m-Rastern vorhanden ist, nicht erreichen / wiederherstellen.

In einem früheren Beitrag hat jemand dies bestätigt und vorgeschlagen, eine andere Funktion für Landsat zu verwenden: 'i.his.rgb' (nach der Verwendung von i.landsat.toar, die ich nicht verwende). Er beginnt, min / max in 1-255 umzuwandeln, bevor er i.his.rgb verwendet ... Meine ursprünglichen Raster haben einen Bereich von 0-65535. Für mich hat es nicht geklappt ...

Ich höre auf, nachdem ich einige Tage lang darum gekämpft habe, nur einen Pansharp auszuführen. Ich fühle mich festgefahren ... Warum sollte die Anwendung von i.pansharp für Landsat nicht funktionieren, wie von anderen vorgeschlagen? Kann jemand einen einfacheren Workflow vorschlagen, um das Pansharpening von Landsat 8-Bildern durchzuführen? Ich freue mich auf einige Ratschläge, danke im Voraus ...

Arnoud Cuppens
quelle

Antworten:

7

Schauen Sie sich i.landsat.rgb an - Führt einen automatischen Farbausgleich für LANDSAT-Bilder durch, wahrscheinlich bevor Sie das Pansharpening ausführen. Sie können auch in Betracht ziehen, die digitalen Zahlen der einzelnen Kanäle mit i.landsat.toar in Strahlung oder Reflexion der Atmosphäre umzuwandeln . Siehe auch http://grasswiki.osgeo.org/wiki/LANDSAT

Übrigens: Ein Bereich von 0-65535 für ursprüngliche LANDSAT-Daten ist ungerade, da sie als 8-Bit-Kanäle (0-255) geliefert werden. Überprüfen Sie also zuerst Ihre Eingabedaten.

Update: 10/2013

Eine vereinfachte Vorgehensweise mit dem aktuellen GRASS GIS 7 finden Sie auf der Seite Geben Sie hier die Bildbeschreibung einVerarbeiten von Landsat 8-Daten in GRASS GIS 7: RGB-Verbundwerkstoffe und Pan-Schärfen

markusN
quelle
Ich bin ursprünglich auf den Artikel gestoßen, den Sie von neteler.org verlinken, aber die Anweisungen haben überhaupt nicht funktioniert (aufgrund des Gleitkommaformats der mit generierten Raster i.landsat.toar, sodass ich den Ratschlägen von hier gefolgt bin
Juan
3
Beachten Sie, dass dieser Rat vor vielen Jahren gegeben wurde und nicht den neuesten GRASS GIS-Code widerspiegelt. Übrigens: Ein wirklich großartiges Pansharpening-Modul ist das neue Addon: i.fusion.hpf - Zusammenführen von hochauflösenden panchromatischen und niedrigauflösenden Multispektraldaten basierend auf der Hochpassfilter-
markusN
4

Danke Michal und Markus. Endlich konnte ich das Pfannenschärfen mit der angegebenen Funktion brov durchführen. Einmal in GRASS ausgeführt, habe ich die RGB-Raster mit dem Datentyp Uint16 nach gdal GTiff exportiert. Markus, ich habe mehrere Landsat 8-Bilder heruntergeladen und verwendet (tiff von jedem der 11 Bänder), um schöne visuelle Ergebnisse zu erzielen, sowohl in natürlicher als auch in falscher Farbe. Meines Wissens kommen die ursprünglichen Landsat-Tiffs als Uint16 mit Werten zwischen 0 und 65535 (oder fehlt mir etwas?) ... Danke für den GRASS LANDSAT-Link ...

Arnoud Cuppens
quelle
3
Ja, Landsat 8-Bilder werden als Uint16-Typen gespeichert.
Webrian
3

Das Schärfen von FIHS oder Brovery Transform ist recht einfach, wenn Sie den Formeln eines Artikels wie diesem folgen: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Angenommen, Sie verwenden GRASS GIS und Pansharpening Landsat 8.

Stellen Sie Ihre Region mit Pan 8 ein - wichtig, um eine Auflösung von 15 m zu haben, die in diesem Fall eingestellt wird.

Für Brovery ohne Gewichte geben Sie dies in den Kartenrechner für Blau, Grün und Rot ein:

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

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

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

Führen Sie r.colors für jede der oben genannten Ausgaben aus und setzen Sie den Histogrammausgleich und die Farbtabelle auf grau.

Führen Sie r.composite aus und wählen Sie die entsprechenden Bänder aus (rot = 4, grün = 3, blau = 2).

Schauen Sie sich das Bild an.

Offensichtlich könnte dieser Prozess unter stärkerer Berücksichtigung der atmosphärischen Korrektur, der Farbdehnung und der Bandgewichte durchgeführt werden, wie im obigen Artikel erörtert.

Ben
quelle
2

Sie können auch Monteverdi ausprobieren , um das Pansharpening durchzuführen. Dieses Tool kann reguläre GeoTIFFs als Eingabe verwenden (separate Multispektralbänder müssen jedoch zu einem Multiband-Raster zusammengeführt werden) und bietet eine recht gute Leistung.

Michal Mackiewicz
quelle
2

Nachdem ich mit diesem Problem zu kämpfen hatte, fand ich eine einfache Lösung, mit der Sie QGIS (2.18.21 ist meine Version) mithilfe von SCP pan-schärfen können. Das Update für SCP mit dieser Methode wird hier beschrieben: https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Das ist einfach und lief wie ein Zauber! Ich habe die endgültige Ausgabe überprüft und die Pixel haben eine Auflösung von 15 m - passend zum panchromatischen 8-Band. Wenn jemand neugierig ist, verwendet er die im Handbuch beschriebene Brovey-Transformationsmethode (siehe https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

Schritt 1 für das SCP-Pansharpen-Tool Schritt 3-4 für SCP Pansharpen

Mark Thompson
quelle
1

Willkommen auf der Seite Janos. Im Idealfall sollten Antworten wie Ihre eine Beschreibung enthalten, warum sie gute Methoden sind. ..

Für das Pansharpening von Landsat 8-Bildern ist die Methode http://www.geosage.com/highview/download.html einfacher . Mein Rat, probieren Sie es aus. Nur ein Klick ....Spektraltransformator

János Breznai
quelle
-2

Für pansharpenig habe ich einen GUI-Spektraltransformator für Landsat 8-Bilder verwendet. Geben Sie hier die Bildbeschreibung ein Http://www.geosage.com/highview/download.html Sehr gut. BR Janos

János Breznai
quelle
2
Willkommen auf der Seite Janos. Im Idealfall sollten Antworten wie Ihre eine Beschreibung enthalten, warum sie gute Methoden sind. Nur-Link-Antworten werden in der Regel ungültig, wenn der Link nicht mehr angezeigt wird.
Aaron