Matlab交互式图像配准实现

clear;
clc;
reference_img=imread('optical1.tif');%提取原始图片作为基准
target_img=imread('optical2.tif');%要求进行匹配的图像
subplot(2,3,1);%显示原始图像
imshow(reference_img);
title('原始图像');
subplot(2,3,2);%显示进行匹配的图像
imshow(target_img);
title('目标图像');
[Ix,Iy]=size(target_img);%取出图片大小
x=[160,103,102,209];%自己选定的特征点
y=[116,246,160,299];
u=[87,35,34,141];
v=[21,151,66,204];
Px=polyfit(x,u,1);%映射函数求解
Py=polyfit(y,v,1);
for i=1:Ix
    for j=1:Iy
        m=i*Px(1,1)+Px(1,2);%映射函数
        n=j*Py(1,1)+Py(1,2);    
        m_integer=floor(m);%对映射值进行取整
        n_integer=floor(n);
        m_decimal=m-m_integer;%对映射值取小数
        n_decimal=n-n_integer;
        if(ge(m_integer,Ix)||ge(n_integer,Iy)||lt(m_integer,0)||lt(n_integer,0)||m_integer==0||n_integer==0)%判断映射点是否满足条件
            result_img(i,j)=0;
            result_img3(i,j)=0;
        else
            result_img(i,j)=target_img(m_integer,n_integer);%没有进行插值