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

hcube = hypercube('indian_pines.dat');

Вычислите значение спектральных индексов для каждого пикселя в кубе данных. По умолчанию в spectralIndices функция возвращает индекс простого отношения (SR), улучшенный индекс растительности (EVI) и модифицированный индекс коэффициента поглощения хлорофилла (MCARI).

indices = spectralIndices(hcube);

Смотрите имена индексов в выход struct indices. Считайте соответствующие изображения индекса, возвращенные на выходе.

indices.IndexName

ans = 
"Simple Ratio (SR)"

ans = 
"Enhanced Vegetation Index (EVI)"

ans = 
"Modified Chlorophyll Absorption Ratio Index (MCARI)"

srImg = indices(1).IndexImage;
eviImg = indices(2).IndexImage;
mcariImg = indices(3).IndexImage;

Оцените растянутое на контраст изображение RGB из исходного куба данных с помощью colorize функция.

rgbImg = colorize(hcube,'Method','RGB','ContrastStretching',true);

Отображение исходных и вычисленных индексных изображений. Значение индекса SR, больше 3, означает растительность. Индекс EVI идентифицирует плотную растительность, и типичное значение индекса EVI для здоровой растительности находится между 0,2 и 0,8. Индекс MCARI означает обилие хлорофилла в области.

fig = figure('Position',[0 0 800 700]);

axes1 = axes('Parent',fig,'Position',[0 0.54 0.42 0.42]);
imagesc(rgbImg,'Parent',axes1);
axis off
title('RGB Image of Data Cube')

axes2 = axes('Parent',fig,'Position',[0.5 0.54 0.45 0.42]);
imagesc(srImg,'Parent',axes2);
axis off
title('SR Image')
colorbar

axes3 = axes('Parent',fig,'Position',[0 0.035 0.45 0.42]);
imagesc(eviImg,'Parent',axes3);
axis off
title('EVI Image')
colorbar

axes4 = axes('Parent',fig,'Position',[0.5 0.035 0.45 0.42]);
imagesc(mcariImg,'Parent',axes4);
axis off
title('MCARI Image')
colorbar

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

hcube = hypercube('jasperRidge2_R198.img');

Вычислите значение MNDWI для каждого пикселя в кубе данных и считайте изображение индекса воды.

indices = spectralIndices(hcube,'MNDWI');
mndwiImg = indices.IndexImage;

Оцените растянутое на контраст изображение RGB из исходного куба данных с помощью colorize функция.

rgbImg = colorize(hcube,'Method','RGB','ContrastStretching',true);

Отображение оригинала и изображения MNDWI.

fig = figure('Position',[0 0 700 400]);
axes1 = axes('Parent',fig,'Position',[0 0.1 0.4 0.8]);
imshow(rgbImg,'Parent',axes1)
title('RGB Image of Data Cube')
axes2 = axes('Parent',fig,'Position',[0.45 0.15 0.47 0.7]);
imagesc(mndwiImg,'Parent',axes2)
colorbar
axis off
title('MNDWI Image')

Водные области обычно имеют значения MNDWI более 0,09. Выполните пороговое сохранение изображения MNDWI для сегментации водных областей. Задайте пороговое значение.

threshold = 0.09;

Сгенерируйте бинарное изображение со значением интенсивности 1 для пикселей с счетом, большей или равной указанному порогу. Все другие пиксели имеют значение 0. Области в бинарном изображении со значением 1 соответствуют областям воды в кубе данных со значениями MNDWI, больше порога.

bw = mndwiImg > threshold;

Наложите бинарное изображение на изображение RGB и отобразите наложенное изображение.

overlayImg = imoverlay(rgbImg,bw,[0 0 1]);
figure
imagesc(overlayImg)
axis off
title('Water Region Overlaid on RGB Image')