GRASSでLANDSAT 8をパンシャープンする方法は?


9

かなり長い間QGISをマッピングに使用していたため、パンシャープンはQGISで直接実行できないため、GRASSを掘り下げることを強制されました。構造をある程度理解した後、私はi.fusion.broveyでいくつかの試行を行いました。15mの解像度でラスターを生成することはできましたが、QGISにロードするときに、元の30mラスターに存在するような良好な色のコントラストを実現/回復できません。

以前の記事誰かことを確認し、LANDSATのために別の関数を使用することが提案されている:「i.his.rgb」(i.landsat.toarを使用した後、私は使用しません)。彼はi.his.rgbを使用する前に最小/最大を1〜255に変換し始めます...私の元のラスターの範囲は0〜65535です。うまくいかなかった...

Pansharpを実行するのに苦労して数日後に終了すると、行き詰まってしまいます...他の人が提案しているように、Landsatでi.pansharpを適用するだけではうまくいかないのはなぜですか?誰かがLandsat 8画像のパンシャープンを実行するより簡単なワークフローを提案できますか?いくつかのアドバイスを楽しみにして、事前に感謝します...

回答:


7

i.landsat.rgbを見て ください -おそらくパンシャープンを実行する前に、LANDSAT画像の色の自動バランスを実行します。また、i.landsat.toarを使用して、個々のチャネルのデジタル番号を大気圏外の輝度または反射率に変換することを検討することもできます。http://grasswiki.osgeo.org/wiki/LANDSAT参照してください

ところで:元のLANDSATデータの範囲が0〜65535であるのは、8ビットチャネル(0〜255)として配信されるため、奇妙です。したがって、最初に入力データを確認してください。

更新:2013年10月

現在のGRASS GIS 7での簡略化された手順については、GRASS GIS 7でのLandsat 8データの処理:RGBコンポジットとパンシャープニングのページを参照してくださいここに画像の説明を入力してください


最初にneteler.orgからリンクしている記事に出くわしましたが、手順はまったく機能しませんでした(で生成されたラスターの浮動小数点形式のi.landsat.toarため、私はここからのアドバイスに従い
Juan

3
このアドバイスは何年も前に行われたものであり、最新のGRASS GISコードを反映していないことに注意してください。ところで:本当に素晴らしいpansharpeningモジュールは、新しいアドオンです:i.fusion.hpf -ハイパスフィルタの追加技術に基づいて、高解像度のパンクロマチックと低解像度のマルチスペクトルデータの融合は、参照grass.osgeo.org/grass70/manuals/をaddons / i.fusion.hpf.html
markusN

4

MichalとMarkusに感謝します。最後に、示された関数brovでパンシャープニングを行うことができました。GRASSで実行したら、RGBラスターをデータ型Uint16としてgdal GTiffにエクスポートしました。マーカス、私はいくつかのランドサット8画像(11バンドのそれぞれのティフ)をダウンロードして使用しており、自然/偽色の両方で素晴らしい視覚的結果を得ています。私の知る限りでは、元のlandsat tiffは0から65535の範囲の値を持つUint16として提供されます(または何か不足していますか?)... GRASS LANDSATリンクに感謝します...


3
はい、Landsat 8画像はUint16タイプとして保存されます。
webrian 2013

3

次のような記事の公式に従うと、パンシャープンFIHSまたはBrovery変換は非常に簡単です。http//www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

GRASS GISを使用してLandsat 8をパンシャープンしていると仮定します。

パン8 tifを使用して地域を設定します-これを行うと設定される解像度15mが重要です。

重みのないBroveryの場合、これをそれぞれマップ計算機に青、緑、赤について入力します。

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

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

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

上記の各出力に対してr.colorsを実行し、ヒストグラム均等化とカラーテーブルを灰色に設定します。

r.compositeを実行し、適切なバンドを選択します(赤= 4、緑= 3、青= 2)

画像を確認してください。

明らかに、このプロセスは、上記の記事で説明したように、大気補正、カラーストレッチ、およびバンドの重みを考慮して行うことができます。


2

また、モンテヴェルディを使用してパンシャープンを実行することもできます。このツールは、通常のGeoTIFFを入力として使用できます(ただし、個別のマルチスペクトルバンドを1つのマルチバンドラスタにマージする必要があります)。


2

この問題に取り組んだ後、SCPを使用してQGIS(2.18.21が私のバージョンです)でパンシャープンできる単純なソリューションを見つけました。この方法を組み込んだSCPの更新は、https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.htmlで説明されています

これは簡単に実行でき、魅力的に実行されました!私は最終出力と15 mの解像度でのピクセル測定を確認しました-パンクロマティック8バンドと一致しています。気になる人がいる場合-マニュアルに記載されているBrovey Transformメソッドを使用します(https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdfを参照)。

SCPパンシャープンツールのステップ1 SCPパンシャープンのステップ3-4


1

サイトJanosへようこそ。理想的には、あなたのような回答には、なぜそれらが優れた方法であるかについての説明が必要です。..

Landsat 8画像のパンシャープンの場合、より簡単な方法はhttp://www.geosage.com/highview/download.htmlです。私のアドバイス、それを試してみてください。ワンクリック...。スペクトル変換器


-2

pansharpenigの場合、Landsat 8の写真にGUIスペクトルトランスフォーマーを使用しましたここに画像の説明を入力してください http://www.geosage.com/highview/download.html 非常に良い BRジャノス


2
サイトJanosへようこそ。理想的には、あなたのような回答には、なぜそれらが優れた方法であるかについての説明が必要です。リンクがなくなると、リンクのみの回答は機能しなくなる傾向があります。
アーロン
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.