removeBands
スペクトル イメージからのスペクトル バンドの削除
説明
newhcube = removeBands(hcube,BandNumber=band)
メモ
この関数には、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 が必要です。
例
ハイパースペクトル データをワークスペースに読み取ります。
hcube = imhypercube("paviaU.dat");hypercube オブジェクトのプロパティを調べます。
hcube
hcube =
hypercube with properties:
ImageDims: "[610x340x103 double]"
Wavelength: [103×1 double]
Metadata: [1×1 struct]
BlockSize: [610 340]
ハイパースペクトル データ キューブのスペクトル波長範囲を求めます。
range = [min(hcube.Wavelength) max(hcube.Wavelength)]
range = 1×2
430 838
ハイパースペクトル データ キューブから削除する波長の範囲を指定します。
wlrange = [410 450; 620 850];
指定した波長範囲にあるスペクトル バンドを削除します。この関数は、削除されたバンドを含まない、新しい hypercube オブジェクトを返します。
newhcube = removeBands(hcube,Wavelength=wlrange)
newhcube =
hypercube with properties:
ImageDims: "[610x340x42 double]"
Wavelength: [42×1 double]
Metadata: [1×1 struct]
BlockSize: [610 340]
元の波長値と新しい波長値をプロットします。
figure plot(hcube.Wavelength,"o") hold on plot(newhcube.Wavelength,"or") xlabel("Band Number") ylabel("Wavelength") legend("Original Values","New Values",Location="SouthEast")
ハイパースペクトル データをワークスペースに読み取ります。
hcube = imhypercube("paviaU.dat")hcube =
hypercube with properties:
ImageDims: "[610x340x103 double]"
Wavelength: [103×1 double]
Metadata: [1×1 struct]
BlockSize: [610 340]
関数 ppi を使用して、ハイパースペクトル データ キューブのスペクトル的に異なる 5 つのエンドメンバーを計算します。
endmembers = fippi(hcube,5);
エンドメンバー スペクトルに基づいて、入力データ キューブの最も情報量の多いバンドを 10 個求めます。
[~,informativeband] = selectBands(hcube,endmembers,NumberOfBands=10);
情報量の多いバンドのバンド番号を使用して、データ キューブの情報量のないバンドのバンド番号を求めます。
datacube = gather(hcube); band = setdiff(1:size(datacube,3),informativeband);
ハイパースペクトル データ キューブから情報量のないバンドを削除します。この関数は、最も情報量の多いバンドのみを含む、新しい hypercube オブジェクトを返します。
newhcube = removeBands(hcube,BandNumber=band)
newhcube =
hypercube with properties:
ImageDims: "[610x340x10 double]"
Wavelength: [10×1 double]
Metadata: [1×1 struct]
BlockSize: [610 340]
元の波長値と新しい波長値をプロットします。
figure plot(hcube.Wavelength,"o") hold on plot(newhcube.Wavelength,"or") xlabel("Band Number") ylabel("Wavelength") legend("Original Values","New Values",Location="SouthEast")
入力引数
出力引数
バージョン履歴
R2020a で導入
参考
hypercube | selectBands | ppi | fippi | nfindr
