removeBands

データキューブからのスペクトルバンドの削除

R2020a 以降

構文

newhcube = removeBands(hcube,'Wavelength',wlrange)

newhcube = removeBands(hcube,'BandNumber',band)

説明

newhcube = removeBands(hcube,'Wavelength',wlrange) は指定した波長範囲内のスペクトルバンドをデータキューブから削除します。この関数は、元のデータキューブから、残りの波長、そのメタデータ情報、対応するスペクトルバンドを含む新しい hypercube オブジェクトを返します。

例

newhcube = removeBands(hcube,'BandNumber',band) は指定されたスペクトルバンド番号を持つスペクトルバンドをハイパースペクトルデータキューブから削除します。

メモ

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

例

すべて折りたたむ

指定した波長範囲にあるバンドの削除

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

Image Processing Toolbox Hyperspectral Imaging Library Image Processing Toolbox Hyperspectral Imaging Library

ライブスクリプトを開く

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

hcube = hypercube('paviaU.dat');

hypercube オブジェクトのプロパティを調べます。

hcube

hcube = 
  hypercube with properties:

      DataCube: [610×340×103 double]
    Wavelength: [103×1 double]
      Metadata: [1×1 struct]

ハイパースペクトルデータキューブのスペクトル波長範囲を求めます。

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:

      DataCube: [610×340×42 double]
    Wavelength: [42×1 double]
      Metadata: [1×1 struct]

元の波長値と新しい波長値をプロットします。

figure
plot(hcube.Wavelength,'o')
hold on
plot(newhcube.Wavelength,'or')
xlabel('Band Number')
ylabel('Wavelength')
legend('Original Values','New Values','Location','SouthEast')

データキューブからの情報のないバンドの削除

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

Image Processing Toolbox Hyperspectral Imaging Library Image Processing Toolbox Hyperspectral Imaging Library

ライブスクリプトを開く

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

hcube = hypercube('paviaU.dat')

hcube = 
  hypercube with properties:

      DataCube: [610×340×103 double]
    Wavelength: [103×1 double]
      Metadata: [1×1 struct]

関数 ppi を使用して、ハイパースペクトルデータキューブのスペクトル的に異なる 5 つのエンドメンバーを計算します。

endmembers = fippi(hcube,5);

エンドメンバースペクトルに基づいて、入力データキューブの最も情報量の多いバンドを 10 個求めます。

[~,informativeband] = selectBands(hcube,endmembers,'NumberOfBands',10);

情報量の多いバンドのバンド番号を使用して、データキューブの情報量のないバンドのバンド番号を求めます。

band = setdiff(1:size(hcube.DataCube,3),informativeband);

ハイパースペクトルデータキューブから情報量のないバンドを削除します。この関数は、最も情報量の多いバンドのみを含む、新しい hypercube オブジェクトを返します。

newhcube = removeBands(hcube,'BandNumber',band)

newhcube = 
  hypercube with properties:

      DataCube: [610×340×10 double]
    Wavelength: [10×1 double]
      Metadata: [1×1 struct]

元の波長値と新しい波長値をプロットします。

figure
plot(hcube.Wavelength,'o')
hold on
plot(newhcube.Wavelength,'or')
xlabel('Band Number')
ylabel('Wavelength')
legend('Original Values','New Values','Location','SouthEast')

入力引数

すべて折りたたむ

`hcube` — 入力ハイパースペクトルデータ
`hypercube` オブジェクト

入力ハイパースペクトルデータ。hypercube オブジェクトとして指定します。hypercube オブジェクトの DataCube プロパティにはハイパースペクトルデータキューブが格納されます。

`wlrange` — 削除する波長範囲
K 行 2 列の行列

削除する波長範囲。K 行 2 列の行列として指定します。K は入力データから削除する波長範囲の数です。各行の形式は [W_min W_max] です。W_min および W_max は削除する範囲の最小波長と最大波長です。指定された波長範囲の少なくとも 1 つは、入力 hypercube オブジェクト内の少なくとも 1 つのスペクトルバンドの波長値と重なり合っていなければなりません。