fastInScene

シーン内高速大気補正の実施

ページ内をすべて折りたたむ

    構文

    newspcube = fastInScene(spcube)

    説明

    newspcube = fastInScene(spcube) は、入力のハイパースペクトル データ内またはマルチスペクトル データ内の大気上端 (TOA) の反射率値から表面反射率値を計算することにより、大気補正されたデータを返します。補正パラメーターの推定にシーン内の特徴を使用する高速大気補正を実施する場合は、この関数を使用します。

    スペクトル データが次に該当する場合は、シーン内高速手法により最良の補正結果が得られます。

    • ラジオメトリック キャリブレーションを実施済みの場合

    • 明るさ条件が均一な場合

    • 大きな水塊、雲、または雲の影が含まれていない場合

    • ベースライン スペクトルの近似値を計算するのに十分な量の暗いピクセルを含む場合

    • 土、水、植物、人工建造物などの異種混合領域を含む場合。この手法では、複数の異なるエンドメンバー スペクトルの平均反射スペクトルはシーンに依存しないと仮定します。

    メモ

    この関数には、Hyperspectral Imaging Library for Image Processing Toolbox™ が必要です。Hyperspectral Imaging Library for Image Processing Toolbox はアドオン エクスプローラーからインストールできます。アドオンのインストールの詳細については、アドオンの入手と管理を参照してください。

    Hyperspectral Imaging Library for Image Processing Toolbox は、MATLAB® Online™ および MATLAB Mobile™ によってサポートされないため、デスクトップの MATLAB が必要です。

    例

    例

    すべて折りたたむ

    この例では次を使用します。

    ライブ スクリプトを開く

    ハイパースペクトル データをワークスペースに読み取ります。

    inputCube = imhypercube("EO1H0440342002212110PY_cropped.dat");

    ハイパースペクトル データ キューブから S/N 比 (SNR) の低いバンドを削除します。

    inputCube = removeBands(inputCube,BandNumber=find(~inputCube.Metadata.BadBands));

    デジタル値 (DN) を大気上端 (TOA) の反射率値に変換します。

    inputCube = dn2reflectance(inputCube);

    シーン内の特徴に基づき、入力ハイパースペクトル データから大気効果を削除します。

    correctedCube = fastInScene(inputCube);

    入力されたハイパースペクトル データおよびシーン内高速大気補正後のハイパースペクトル データのフォールスカラー イメージを表示します。

    colorImg = colorize(inputCube);
colorImg_corrected = colorize(correctedCube);
figure
imagesc([colorImg colorImg_corrected])

    Figure contains an axes object. The axes object contains an object of type image.

    入力引数

    すべて折りたたむ

    入力スペクトル データ。hypercube オブジェクトまたは multicube オブジェクトとして指定します。より良い結果を得るには、入力値が TOA の反射率値でなければなりません。入力値がデジタル値の場合は、関数 dn2reflectance を使用して TOA の反射率値を推定します。

    spcube が multicube オブジェクトである場合、そのすべてのスペクトル バンドのデータ解像度は同じでなければなりません。multicube オブジェクトのスペクトル バンドの解像度がすべて同じではない場合、resampleBands 関数を使用して帯域を再サンプリングするか、selectBands 関数を使用して均一な解像度の帯域を選択します。

    出力引数

    すべて折りたたむ

    出力スペクトル データ。hypercube オブジェクトまたは multicube オブジェクトとして返されます。出力で返されるデータ キューブのピクセル値は、表面反射率値を指定します。

    制限

    この関数は、パフォーマンスが既に最適化されているため、parfor ループをサポートしません。 (R2023a 以降)

    参照

    [1] Bernstein, L.S., S.M. Adler-Golden, R.L. Sundberg, R.Y. Levine, T.C. Perkins, A. Berk, A.J. Ratkowski, G. Felde, and M.L. Hoke. “A New Method for Atmospheric Correction and Aerosol Optical Property Retrieval for VIS-SWIR Multi- and Hyperspectral Imaging Sensors: QUAC (QUick Atmospheric Correction).” In Proceedings. 2005 IEEE International Geoscience and Remote Sensing Symposium, 2005. IGARSS ’05., 5:3549–52. Seoul, Korea: IEEE, 2005. https://doi.org/10.1109/IGARSS.2005.1526613.

    バージョン履歴

    R2020b で導入

    すべて展開する

    参考

    | | |