帮我写一段完整详细的matlab提取高光谱反射率敏感波段的代码

在MATLAB中，提取高光谱反射率敏感波段通常涉及到图像处理和数据分析。以下是一个简单的示例代码，用于从高光谱影像中选择几个特定波段作为敏感区域。假设我们已经有了一个三维数组`HSI`，其中第一维代表波长，第二、三维分别代表空间维度。

% 假设HSI是高光谱影像数据，例如Bands是波长信息，Data是RGB或全波段数据
Bands = HSI(:, 1); % 获取波段信息
Data = HSI(:, 2:end, :); % 提取波段数据

% 定义感兴趣的敏感波段范围，这里以红、绿、蓝波段为例
Wavelengths_of_interest = [400, 570]; % 红色波段
Green_Wavelengths = [500, 680];   % 绿色波段
Blue_Wavelengths = [450, 530];   % 蓝色波段

% 找到对应波段的索引
Red Indices = find(Bands >= Wavelengths_of_interest(1) & Bands <= Wavelengths_of_interest(2));
Green Indices = find(Bands >= Green_Wavelengths(1) & Bands <= Green_Wavelengths(2));
Blue Indices = find(Bands >= Blue_Wavelengths(1) & Bands <= Blue_Wavelengths(2));

% 提取敏感波段的数据
Sensitive_Bands_R = Data(:, Red_Indices, :);
Sensitive_Bands_G = Data(:, Green_Indices, :);
Sensitive_Bands_B = Data(:, Blue_Indices, :);

% 结果存储
Sensitive Bands = [Sensitive_Bands_R; Sensitive_Bands_G; Sensitive_Bands_B];

% 可视化或进一步分析
figure;
subplot(1, 3, 1), imshow(Sensitive_Bands_R, []);
title('Sensitive Red Band');
subplot(1, 3, 2), imshow(Sensitive_Bands_G, []);
title('Sensitive Green Band');
subplot(1, 3, 3), imshow(Sensitive_Bands_B, []);
title('Sensitive Blue Band');

%