一种简单的RX算法程序

一种典型的RX算法程序，可与PCA,ICA等程序联合使用，用于目标检测。

function outputimage=RX（A）
[mnq]=size（A）;
newimage= reshape（ Am*nq）;
lastimage=newimage‘;
meanRow=mean（lastimage2）;%求每一行的均值
% meanRow=round（meanRow）;%取整
% C=zeros（qq）;
% for i=1:m*n 
%      lastimage（:i）= lastimage（:i）-meanRow;
%     
% C= lastimage*lastimage‘+C; 
% Csum=C;
% end
% Wcov=Csum/（m*n）;%协方差矩阵

for i=1:m*n 
     lastimage（:i）= lastimage（:i）-meanRow;
end
Wcov=cov（lastimage‘）;
invWcov=inv（Wcov）/（m*n-1）;


rx=zeros（1m*n）;
for i=1:m*n
    rx（i）=lastimage（:i）‘* invWcov *lastimage（:i）;  
end
outputimage=reshape（rxmn）;
outputimage=mat2gray（outputimage）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        627  2014-01-02 09:26  RX.m

----------- ---------  ---------- -----  ----

                  627                    1