高光谱异常探测算法RXD的matlab实现

function Result = RXD（HIM）
% By Yulei Wang （last modified on 03/18/2014）
% RXD: anomaly detector
%       RXD performs the commonly used RX anomaly detector based on
%       covariance matrix
%       This funcion is a global RX detector in which covariance matrix is
%       calculated only once using all the data points
% Usage:
%   Result = RXD（HIM）
% Inputs:
%   HIM - 3D data matrix （XX*YY*bnd） which is the original hyperspectral
%         imagery data
% Outputs:
%   Result - Detector output （XX*YY）

% ------------------------ RX Algorithm --------------------------------------
display（‘start RXD algorithm‘）;%显示
t1 = cputime;
[XXYYbnd] = size（HIM）;% XX为图像的长度 YY为图像的高度 bnd为图像的波段数
r = reshape（HIMXX*YYbnd）‘;%将HTM中各个像素点的像素值放入到r（行数bnd 列数xx*yy）矩阵中