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

hcube = hypercube('indian_pines.dat');

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

indices = spectralIndices(hcube);

Проверьте имена индексов в структуре вывода 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')