insar图像配准及滤波

function [moverowmovecol]=rcorr（im11im22）
%功能：对输入的im2矩阵进行平移，并返回im1im2相关系数最大时im2的平移量.
rmove=zeros（100100）;%平移所对应的相关系数
for i1=-2:2
    for j=-20:2 %负号代表左移和上移
        %%  获得一次平移后pq的im2矩阵im2move
        if  （i1>=0&&j>=0） %下移i行右移j列
            im1move=zeros（50-i150-j）;%主窗口
            im2move=zeros（50-i150-j）;%移动窗口
            im1move=im11（（i1+1）:50（j+1）:50）;
            im2move=im22（1:50-i11:50-j）;
        else
            if  （i1>=0&&j<0） %下移i行左移j列
                im1move=zeros（50-i150+j）;%主窗口
                im2move=zeros（50-i150+j）;%移动窗口
                im1move=im11（（i1+1）:501:50+j）;
                im2move=im22（1:50-i1（-j+1）:50）;
            else
                if  （i1<0&&j>=0）%上移i行右移j列
                    im1move=zeros（50+i150-j）;%主窗口
                    im2move=zeros（50+i150-j）;%移动窗口
                    im1move=im11（1:50+i1（j+1）:50）;
                    im2move=im22（（-i1+1）:501:50-j）;
                else
                    
                    im1move=zeros（50+i150+j）;%主窗口
                    im2move=zeros（50+i150+j）;%移动窗口
                    im1move=im11（1:50+i11:50+j）;
                    im2move=im22（（-i1+1）:50（-j+1）:50）;
                end
            end
        end
        
        %%将该次平移得到的矩阵进行求相关系数并存储到rmove中
        im=im1move.*conj（im2move）;
        [mn]=size（im1move）;
        c=1;d=1; Eim1=0;Eim2=0;w=0;
        for c=1:m
            for d=1:n
                w=im（cd）+w;
                Eim1= Eim1+im1move（cd）*conj（im1move（cd））;
                Eim2=Eim2+im2move（cd）*conj（im2move（cd））;
            end
        end
        rmove（i1+10j+10）=abs（w）/sqrt（Eim1*Eim2）;%由平移得到的相关系数矩阵
    end
end
[xy]=find（rmove==max（max（rmove）））;%返回最大move对应的坐标
moverow=x-10;
movecol=y-10;

end
   
  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-12-17 10:56  insar图像配准及滤波
     文件     9286147  2018-11-29 16:52  insar图像配准及滤波im1.mat
     文件    15357756  2018-11-29 16:52  insar图像配准及滤波im2.mat
     文件       27843  2018-12-14 18:34  insar图像配准及滤波mmmove.mat
     文件       28022  2018-12-14 18:35  insar图像配准及滤波
nmove.mat
     文件        1950  2018-12-14 23:34  insar图像配准及滤波
corr.m
     文件        5623  2018-12-15 09:59  insar图像配准及滤波sar_phase_main.m
     文件    34132971  2018-12-17 10:50  insar图像配准及滤波第四次作业.docx
     文件         170  2019-05-08 22:30  说明.txt