spectraMatch.rar_spectral angle_光谱匹配 matlab_光谱匹配程序_匹配 matlab_角度匹

function sid_sa=SID_SA(t,r) % t、r 都是 Nx1 矩阵,p为参数 'tan' or 'sin' sa=acos((t'*r)/(norm(t)*norm(r))); % SA Spectral Angle tm=t./sum(t); rm=r./sum(r); sid=(tm-rm)'*(log10(tm)-log10(rm)); sid_sa=sid*sin(sa); % elseif(strcmp(p,'sin')) % sid_sa=sid*sin(sa); % else % error('You can only use ''tan'' or ''sin''');