ArcMapとQGISでSentinel 2を使用するNDVIに違いがあるのはなぜですか？

Sentinel 2バンド4（赤）と8（NIR）を使用してシーンのNDVIを計算しようとしています。最初はArcMapで行い、次にQGISで行いました。どちらのソフトウェアでも、ラスター計算機を使用し、基本的なNDVI式=（NIR-Red）/（NIR + Red）を含めました。

ArcMapの範囲は0〜9221（これは適用されたストレッチによって異なります）でしたが、QGISの場合は0.0258〜0.5255のより妥当な値が得られました。

Sentinel 2はTOAの反射率であり、スケーリング係数が10000であることを簡単に読みました。これは、画像と一緒にダウンロードされるxmlファイルにあります。ダウンロードしたxmlファイルに、このスケーリングファクターが記載されていません。

なぜこの違いがあるのですか？QGISで計算はどのように行われますか？

— tg110
ソース

ラスター計算機またはツール/プラグインを使用して、NDVIをどのように計算しますか？Arcmapの計算で問題が発生したことは明らかです。投稿を編集して、計算方法を正確に説明できますか？

— HDunn

@ HDunn-私はあなたが提案したように詳細を含めました。プラグインは使用せず、ArcMapとQGISの両方でラスター計算機を使用しました。

— tg110 2016

それらを8/16ビットではなく浮動小数点として処理するようにArcmapを強制してみてください。そのためには、引数の1つに1.0を掛けます。.0に注意してください。（（NIR-Red）* 1.0）/（NIR + Red）のようなもの。

— HDunn

@ HDunn-うまくいきました！しかし、ArcMapの値の範囲-0.75468〜0.877372は、QGISを使用して取得した範囲（0.0258〜0.5255）とは大きく異なります。どちらが正しいですか。そして、なぜこの違いがあるのですか？

— tg110 2016

ストレッチは似ていますか？統計の計算を試すか、両方のプログラムで個々のピクセルを識別してみてください。同じであるか、または桁を丸めて正確に丸める必要があります

— HDunn

回答:

Arcmapは、ラスタータイプを保持し、16ビット整数ラスターから計算したときに値を小数/浮動小数点として扱わないことで悪名高くなっています。

Arcmapがラスター計算結果を浮動小数点として扱うように強制するには、引数の1つに1.0を掛けるか、式の定数または実数に単に「.0」を追加する必要があります。

あなたの場合に(NIR-Red)*1.0/(NIR+Red)は十分でしょう。

— HDunn
ソース

@ HDunn-同じストレッチを適用し、統計も計算しましたQGIS-最小= -0.684、最大= 0.827、平均= 0.252、標準偏差= 0.113 ArcMap-最小= -0.7546、最大= 0.87737、平均= 0.252、標準偏差= 0.113つまりほぼ同じ。なぜまだ微妙な違いがあるのかしら。

— tg110

使用しているそれぞれのバンドの初期ヒストグラムを確認しましたか？各ソフトウェアが異なるno-data値を解釈して、結果を歪める可能性があります。

— アレクサンダー
ソース

@アレクサンダー-それは私が持っていた別の質問です。QGISのNIRバンドの値は527-3234ですが、ArcMapでは0-9212です。方法と理由がわかりません。

— tg110 2016

これが問題の原因である可能性があります。ヒストグラムを表示すると、データのない値がわかります。これは、ヒストグラムの最初または最後に表示されるため、多くの場合、ピクセル数が多いためです。あなたがする必要があるのは、ヒストグラムクリップであり、データなしの値を切り取り、インデックスを再実行して、結果をお知らせください:)

— Alexander

@ Alexander-ヒストグラムクリップの方法がわかりません。どうすればいいですか？

— tg110

hmmm-私はerdasでそれを行う方法を知っています。これは、Rescale関数を使用することです。ただし、Arcmapでバンド計算機を使用せずにNDVIを試してみたことはありますか？：このリンクをチェックdesktop.arcgis.com/en/arcmap/10.3/manage-data/raster-and-images/...

— アレクサンダー

基本的に、画像分析ウィンドウでヒストグラムを再スケーリングできます（多くの場合、デフォルトでデータなしが除外されます）-画像分析オプションメニューに移動して、正しいバンドが割り当てられていることを確認してください。後でndviを試して、機能するかどうかを確認してください。

— Alexander