SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，.zip

SAR图像变化检测方法，包括了对数比和均值比两种较为经典的仿真。本方法从图像像素的角度出发，对图像像素的灰度值进行操作。分别才用了均值法和对数比法，然后。本代码适合新手使用，能让初学者从图像像素灰度值的角度去理解SAR图像形变监测的理论

% BFILTER2 Two dimensional bilateral filtering.
%    This function implements 2-D bilateral filtering using
%    the method outlined in:
%
%       C. Tomasi and R. Manduchi. Bilateral Filtering for 
%       Gray and Color Images. In Proceedings of the IEEE 
%       International Conference on Computer Vision 1998. 
%
%    B = bfilter2（AWSIGMA） performs 2-D bilateral filtering
%    for the grayscale or color image A. A should be a double
%    precision matrix of size NxMx1 or NxMx3 （i.e. grayscale
%    or color images respectively） with normalized values in
%    the closed interval [01]. The half-size of the Gaussian
%    bilateral filter window is defined by W. The standard
%    deviations of the bilateral filter are given by SIGMA
%    where the spatial-domain standard deviation is given by
%    SIGMA（1） and the intensity-domain standard deviation is
%    given by SIGMA（2）.
%
% Douglas R. Lanman Brown University September 2006.
% dlanman@brown.edu http://mesh.brown.edu/dlanman
% 我理解的意思是B=bfilter的第一个参数是待处理的图片，第二个参数是算子长度，SIGMA（1）是空间的公式里的sigma，SIGMA（2）是公式像素值的SIGMA

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Pre-process input and select appropriate filter.
% 下面几个if判断都是判断输入的参数是否正确
function B = bfilter2（Awsigma）

% Verify that the input image exists and is valid.
if ~exist（‘A‘‘var‘） || isempty（A）
   error（‘Input image A is undefined or invalid.‘）;
end
if ~isfloat（A） || ~sum（[13] == size（A3）） || ...
      min（A（:）） < 0 || max（A（:）） > 1
   error（[‘Input image A must be a double precision ‘...
          ‘matrix of size NxMx1 or NxMx3 on the closed ‘...
          ‘interval [01].‘]）;      
end

% Verify bilateral filter window size.
% nume1（w）判断数组w中元素个数
if ~exist（‘w‘‘var‘） || isempty（w） || ...
      numel（w） ~= 1 || w < 1
   w = 5;
end
% 将w朝正四舍五入
w = ceil（w）;

% Verify bilateral filter standard deviations.
if ~exist（‘sigma‘‘var‘） || isempty（sigma） || ...
      numel（sigma） ~= 2 || sigma（1） <= 0 || sigma（2） <= 0
   sigma = [3 0.1];
end

% Apply either grayscale or color bilateral filtering.
if size（A3） == 1
   B = bfltGray（Awsigma（1）sigma（2））;
else
   B = bfltColor（Awsigma（1）sigma（2））;
end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Implements bilateral filtering for grayscale images.
function B = bfltGray（Awsigma_dsigma_r）

% Pre-compute Gaussian distance weights.
[XY] = meshgrid（-w:w-w:w）;
G = exp（-（X.^2+Y.^2）/（2*sigma_d^2））;
% Create waitbar.
h = waitbar（0‘Applying bilateral filter...‘）;
set（h‘Name‘‘Bilateral Filter Progress‘）;

% Apply bilateral filter.
dim = size（A）;
B = zeros（dim）;
for i = 1:dim（1）
   for j = 1:dim（2）
      
         % Extract local region.
         iMin = max（i-w1）;
         iMax = min（i+wdim（1））;
         jMin = max（j-w1）;
         jMax = min（j+wdim（2））;
         I = A（iMin:iMaxjMin:jMax）;
      
         % Compute Gaussian intensity weights.
         H = exp（-（I-A（ij））.^2/（2*sigma_r^2））;
      
         % Calculate bilateral filter response.
      

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-08-29 20:42  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\
     文件        4905  2018-01-16 12:11  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\bfilter2.m
     文件           0  2019-08-27 17:23  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\cols2.m
     文件        2863  2019-08-26 20:06  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\LOG_D_image.m
     文件        3485  2019-08-26 20:03  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\Maxture.m
     文件        3504  2019-08-26 20:00  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\MEAN_D_image.m
     文件      305254  2016-04-15 16:26  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\sar1gray.bmp
     文件      305254  2016-04-15 16:26  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\sar2gray.bmp
     文件      305254  2016-04-15 16:26  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\sarreference.bmp
     文件           0  2019-08-26 17:51  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\Xm.m
     文件      103278  2019-08-27 17:41  SAR图像变化检测方法，包括了对数比和均值比两种较为经典的，\Xm2.bmp