SAR图像变化检测代码，基于差异图和比值图融合的.zip

SAR图像变化检测代码，基于差异图和比值图融合的代码。本方法从图像像素的角度出发，对图像像素的灰度值进行操作。分别才用了差值法和比值法，然后在各自的基础上对代码进行融合，这样就避免了差值法和比值法各种带来的缺点。本代码适合新手使用，能让初学者从图像像素灰度值的角度去理解SAR图像形变监测的理论。

clear all;
close all;

f=imread（‘1999.04.bmp‘）;
g=imread（‘1999.05.bmp‘）;
figureimshow（f）;
figureimshow（g）;
f=double（rgb2gray（f））; %意思是将真彩色图像RGB转换为灰度强度图像I
g=double（rgb2gray（g））;

[ab]=size（f）;        %读取格式
A1=abs（f-g）;
figureimshow（uint8（A1））
A=255-abs（f-g）; %A为差值法做出来图像变换信息黑白互换
figureimshow（uint8（A））
B1=g./f;
[cd]=size（B1）;
for i=1:c
    for j=1:d
        if B1（ij）>1.3 || B1（ij）<0.7
            B2（ij）=1;
        else B2（ij）=0;
        end
    end
end
figureimshow（B2）  

B=255*min（f./gg./f）; %B为比值法求出来的图像变换信息
C=A.*B/max（max（B））;   %将差值法和比值法融合得到的C
figureimshow（uint8（C））
[ab]=size（C）;

for i=1:a
    for j=1:b
        if C（ij）>0 && C（ij）<100
            C1（ij）=1;
        else C1（ij）=0;
        end
    end
end

figureimshow（uint8（C））;%unit格式，就是二值图像无符号8位二进制整型
figureimshow（C1）;  %融合出来的图像差异
fim=uint8（C）;

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

     文件     272158  2008-02-01 09:47  1999.04.bmp

     文件     272158  2008-02-01 09:48  1999.05.bmp

     文件        416  2010-04-08 11:19  fusion.m

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

               544732                    3