Landsat8 DNからTOAへの反射率(マイナス値)


8

大気補正を手動で行う方法を学んでいるため、ENVIで利用可能な自動化された方法を使用していません

Landsat 8シーンのメタデータ(バンド1)は次のとおりです。

RADIANCE_MULT_BAND_1 = 1.2344E-02

RADIANCE_ADD_BAND_1 = -61.72016

REFLECTANCE_MULT_BAND_1 = 2.0000E-05

REFLECTANCE_ADD_BAND_1 = -0.100000

SUN_ELEVATION = 48.12237303

したがって、バンド数学を使用してTOA放射輝度を計算する場合、次のように入力します。

(0.012344 * B1)-61.72016

そして、それは一種の正しいように見えるスペクトル統計をもたらします SpectralStats

次のバンドの数式を使用すると、

TOA反射率(太陽高度調整を除く)

(B1 * 0.00002)-0.1

TOA反射率(太陽高度調整を含む)

((B1 * 0.00002)-0.1)/0.74457226676389733207607359928648

または、生のDNからTOA反射率(太陽高度調整を含む)まで直接計算します。

((((0.012344 * B1)-61.72016)* 0.00002)-0.1)/0.74457226676389733207607359928648

私はこのようなものになってしまいます:

TOA反射率

正しくないように見えます。

私はググっていて、私の検索語がごみであるか、手動バンドの計算を使ってこれをしている人はいません。

[編集]

以下のMikkelの解答と修正されたバンド数学計算(中間ステップを省略)を使用すると、次の統計が得られました。

ここに画像の説明を入力してください

ややリアルに見えます。

各バンドに13.4305を追加する必要があるか、または負の最小値を無視する必要があるかどうか疑問に思いますか?

[編集を終了]

[編集#2]

Mikkelは再び注目され、-13.4305の値は画像の「データがない」領域に分離されたことを示唆しています。

ここに画像の説明を入力してください

ご覧のとおり、シーンの北端のみに孤立しています。

[編集#2を終了]

回答:


7

あなたのアプローチの一つの欠陥。輝きを得るためにDNを通過する必要はありません。
DNから反射率に直接進むことができます。
だけに固執する((B1*0.00002)-0.1)/0.74457226676389733207607359928648


これはかなりうまくいきました!TOA放射輝度/ TOA反射率は除外されると説明されたことがないと思います。太陽の上昇は、補正シーケンスでは必要な手順ではありませんでした。
kyuss

上記の編集された投稿によれば、各バンドに13.4305を追加してmin = 0にする必要がありますか、それとも不要ですか?
kyuss

@kyuss、-0.134305がどこにあるか確認してみてください。画像の端の周りのデータがない領域だと思います。
Mikkel Lydholm Rasmussen、2015

素晴らしい、私は何もチェックしなかったが、負の値がなかった南の境界線。あなたのすべての助けに感謝します。
kyuss

@kyussは、プロセスの次のステップは「単純に」大気補正を行うことです。
Mikkel Lydholm Rasmussen、2015

2

これが大気補正の物理学に関する良い読みです、そして問題を助けるかもしれません...私はボタンを押す男です、私はあなたの痛みを感じます、あなたが外に出たら、あなたはその美しいQUACまたはFLAASHボタンをクリックするだけです。学界の。

http://modis.gsfc.nasa.gov/data/atbd/atbd_mod08.pdf


スティーブンに感謝します。私は-あなたが推論したように-私の研究中に多くを学びましたが、そのmodis文書は、放射伝達に関する私の理解のギャップを埋めるのに役立つかもしれません。
kyuss
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.