ChangeDetectionPCAkMeans 用于变化检测的方法

用于变化检测的方法，是基于PCA+k-means算法实现的

function  [FAMAOECA]=DAcom（refImagetestImage）
% DAcom.m
% Usage: [FAMAOE]=DAcom（refImagetestImage）
% Inputs:
%   refImage=reference map
%   testImage=detection map
% Outputs:
%   FA=False alarms
%   MA=Missed alarms
%   OE=Overall Error
% July152009
% Copyright（C） 2008-2009 by Fan
%-------------------------------------------------------------------------
if isempty（refImage）
    error（‘!!!Not exist reference map‘）;
end
refImage（find（refImage>=128））=255; %% reference image binaryzation 
refImage（find（refImage<128））=0;
RI=refImage（:）;TI=testImage（:）; 
aa=find（RI==0&TI~=0）;bb=find（RI~=0&TI==0）; %% Calculation of indicators
FA=numel（aa）;MA=numel（bb）;OE=FA+MA; 
[mn]=size（testImage）;
CA=1-OE/（（m+4）*（n+4））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-12-06 16:45  ChangeDetectionPCAkMeans\
     目录           0  2012-12-11 10:55  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\
     文件      272158  2008-02-01 09:47  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\04.bmp
     文件      272158  2008-02-01 14:35  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\045.bmp
     文件      272158  2008-02-01 09:48  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\05.bmp
     文件      125662  2012-09-17 15:01  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\1.bmp
     文件      125934  2012-09-17 15:01  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\2.bmp
     文件      372158  2012-09-17 15:01  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\3.bmp
     文件       91078  2009-08-24 16:12  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\BangladeshT1.bmp
     文件       91078  2009-08-24 16:12  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\BangladeshT2.bmp
     目录           0  2012-12-06 16:45  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\
     目录           0  2012-12-09 00:19  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\
     文件      272158  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\04.bmp
     文件      272158  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\045.bmp
     文件      272158  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\05.bmp
     文件      125662  2012-12-06 16:56  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\1.bmp
     文件      125934  2012-12-06 16:56  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\2.bmp
     文件      372158  2012-12-06 16:56  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\3.bmp
     文件      125662  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\A1.bmp
     文件      125934  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\A2.bmp
     文件      372158  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\A3.bmp
     文件      786486  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\B1.bmp
     文件           0  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\B2.bmp
     文件      263222  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\B3.bmp
     文件       91078  2009-08-24 16:12  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\BangladeshT1.bmp
     文件       91078  2009-08-24 16:12  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\BangladeshT2.bmp
     文件      370854  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\C1.bmp
     文件      370854  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\C2.bmp
     文件      370854  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\C3.bmp
     文件      405598  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\D1.bmp
     文件      405598  2012-12-09 00:09  ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\ChangeDetectionPCAkMeans\D2.bmp
............此处省略40个文件信息