遥感图像变化检测算法

可以作为遥感图像变化检测学习的资源，有变化检测部分，有阈值分割部分

function [dispersionmeanw] = covw（xwvarargin）

% [covaw meanw] = covw（xwvarargin）
%
%   COVW Covariance matrix with weights.
%
%   COVW（WX） if X is a vector returns the variance.  For matrices
%   where each row is an observation and each column a variable
%   COVW（XW） is the variance-covariance matrix.  DIAG（COVW（XW）） is
%   a vector of variances for each column and SQRT（DIAG（COVW（XW）））
%   is a vector of standard deviations. 
%
%   Observations are weighted with vector of weights W which has all
%   nonnegative values.  COVW（XW） gives the weighted estimate of the
%   variance-covariance matrix.
%   
%   COVW（XW） normalizes by （N-1） where N is the number of
%   observations.
%
%   COVW（XW1） normalizes by N and produces the second
%   moment matrix of the observations about their mean.
%   COVW（XW0） is the same as COVW（XW）.
%
%   The weighted mean is removed from each column before calculating the
%   result.
%
%   See also COV CORRCOEF STD MEAN.

%   J. Little 5-5-86
%   Revised 6-9-88 LS 3-10-94 BJ
%   Copyright 1984-2002 The MathWorks Inc. 
%   $Revision: 5.16 $  $Date: 2002/06/05 17:06:38 $
%
%   Modified for weights
%   Allan Aasbjerg Nielsen
%   aa@imm.dtu.dk
%   10 Feb 2005

if nargin<2 error（‘Not enough input arguments.‘）; end
if nargin>3 error（‘Too many input arguments.‘）; end
if ndims（x）>2 error（‘x must be 1- or 2-D.‘）; end
w = w（:）; % w is now a column vector

[mn] = size（x）;

[mwnw] = size（w）;
if length（x）==（m*n）
    x = x（:）; % x is now a column vector
    [mn] = size（x）;
end
if m~=mw error（‘x and w must match.‘）; end 
if length（w（w<0））>0 error（‘Weights must be nonnegative.‘）; end

flag = 0;
% Check for covw（xw） or covw（xwflag）
if （nargin==3）
    flag = varargin{end};
end

if m==1  % Handle special case
    dispersion = 0;
else
    sumw = sum（w）;
    meanw = sum（repmat（w1n）.*x） / sumw;
    %xc = x - repmat（sum（repmat（w1n）.*x）m1） / sumw;  % Remove weighted mean
    xc = x - repmat（meanwm1）;  % Remove weighted mean
    % direct method uses way too much memory:
    %dispersion = xc‘ * diag（w） * xc / sumw;
    xc = repmat（sqrt（w）1n） .* xc;
    dispersion = xc‘ * xc / sumw;
    if flag==0
        dispersion = m/（m-1） * dispersion;
    end
end

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

     文件       2311  2016-10-31 15:10  Change Detection Code\covw.m

     文件       1145  2016-11-01 16:35  Change Detection Code\createfigure.m

     文件       1314  2016-11-14 16:18  Change Detection Code\CVADemo.m

     文件       4272  2016-11-15 11:16  Change Detection Code\DAcom.m

     文件      35124  2016-10-28 20:04  Change Detection Code\data\major-TestROI.roi

     文件     960000  2012-04-25 09:29  Change Detection Code\data\Taizhou\2000TM

     文件       1376  2012-04-25 09:29  Change Detection Code\data\Taizhou\2000TM.HDR

     文件     960000  2012-04-25 09:28  Change Detection Code\data\Taizhou\2003TM

     文件       1376  2012-04-25 09:28  Change Detection Code\data\Taizhou\2003TM.HDR

     文件      35124  2016-10-28 20:04  Change Detection Code\data\Taizhou\major-TestROI.roi

     文件     160000  2016-10-31 17:01  Change Detection Code\data\Taizhou\TaizhouChange

     文件        923  2016-10-31 17:01  Change Detection Code\data\Taizhou\TaizhouChange.hdr

     文件     161078  2016-10-31 17:04  Change Detection Code\data\Taizhou\TaizhouChange_blackWhite.bmp

     文件        222  2016-10-31 17:04  Change Detection Code\data\Taizhou\TaizhouChange_blackWhite.bpw

     文件     161078  2016-11-02 10:26  Change Detection Code\data\Taizhou\TaizhouChange_blackWhite_unchange.bmp

     文件        222  2016-11-02 10:26  Change Detection Code\data\Taizhou\TaizhouChange_blackWhite_unchange.bpw

     文件     160000  2016-11-02 16:22  Change Detection Code\data\Taizhou\TaizhouChange_overall

     文件        993  2016-11-02 16:22  Change Detection Code\data\Taizhou\TaizhouChange_overall.hdr

     文件        247  2005-01-27 12:14  Change Detection Code\eigen2.m

     文件       2147  2016-10-28 20:30  Change Detection Code\enviwrite.m

     文件       2139  2016-10-29 17:40  Change Detection Code\enviwrite2.m

     文件       3433  2016-10-18 20:07  Change Detection Code\freadenvi.m

     文件       7071  2016-11-14 15:16  Change Detection Code\IRMAD_Update.asv

     文件       7275  2016-11-15 11:19  Change Detection Code\IRMAD_Update.m

     文件        881  2016-11-04 15:50  Change Detection Code\KmeansMap.m

     文件       4120  2016-11-15 11:15  Change Detection Code\MADGet.m

     文件       1683  2016-09-28 11:13  Change Detection Code\PCA.m

     文件        910  2016-11-14 16:19  Change Detection Code\PCADemo.m

     文件     541222  2016-11-01 16:34  Change Detection Code\Test result\All Bands CCA EigValues Tendency 10 Iterations_Taizhou.bmp

     文件       9238  2016-11-01 16:34  Change Detection Code\Test result\All Bands CCA EigValues Tendency 10 Iterations_Taizhou.fig

............此处省略65个文件信息