光谱角matlab代码（Spectral Angle Mapper,SAM）

以下是使用 Matlab 计算光谱角（Spectral Angle Mapper，SAM）的示例程序，你可以根据自己的需求进行修改和调整。

% Load reference and test spectra
ref_spec = double(imread('ref_spectrum.tif')); % reference spectrum
test_spec = double(imread('test_spectrum.tif')); % test spectrum

% Normalize spectra
ref_spec = ref_spec./sum(ref_spec(:));
test_spec = test_spec./sum(test_spec(:));

% Calculate spectral angle
numerator = dot(ref_spec,test_spec);
denominator = norm(ref_spec)*norm(test_spec);
theta = acosd(numerator/denominator);

% Output
fprintf('Spectral angle: %f degrees\n', theta);

注意，该示例程序只适用于单波段光谱。如果你需要计算多波段光谱角，请参考相关文献或使用现有的光谱角计算工具。此外，你需要根据实际情况调整输入和输出的数据格式和文件路径。