Считывание гиперспектральных данных в рабочую область.

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')

Считывание гиперспектральных данных в рабочую область.

hcube = hypercube('paviaU.dat')

hcube = 
  hypercube with properties:

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

Вычислите пять спектрально различающихся конечных элементов гиперспектрального куба данных с помощью ppi функция.

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')