图像配准

作业三图像配准

一、实验题目：

图(a)是一幅参考图像，图(b)是因垂直和水平切变产生几何畸变的输入图像，选用合适的约束点配准这两幅图像。

注：参考图像已提供，由参考图像，可通过水平偏移20像素和垂直偏移200像素得到输入图像（图(b)）

(a) 参考图像(b) 输入图像

二、实验代码：

function Main()

baseIm= imread('Fig0237.tif');%读入基准图像

[inputIm,map] = imread('t.tif');%读入待配准图像

A=[621 257];B=[403 443];C=[644 823];D=[472 606];%select_point函数选取坐标点

A1=[620 75];B1=[393 329];C1=[622 640];D1=[456 471];

input_points=[A;B;C;D];

base_points=[A1;B1;C1;D1];

%input_points_corr = cpcorr(input_points,base_points,input,base);

%%选择变换类型，不同的变换类型配准精度不一，选择最优配准即可

type ='affine';% 'similarity'%'affine' %'piecewise

linear'%'nonreflectivesimilarity'%'piecewise linear'%;

txpz(baseIm,inputIm,input_points,base_points,type);

End

function txpz(baseIm,inputIm,input_points,base_points,type)

%%求出从待配准图到基准图的坐标变换关系

T = cp2tform(input_points,base_points,type);

%%将待配准图放到基准图的坐标系下，xRange是横坐标的取值范围，yRange是纵坐标的取值范围,

%%Reg是配准后的图像

[Reg xRange yRange] =imtransform(inputIm,T);

[m1,n1,p1] = size(baseIm);

[m2,n2,p2] = size(Reg);

%从配准后的大图中截取与基准图等大的局部图像

rowStart = ceil(-yRange(1));

rowEnd = floor(-yRange(1)+m1);

colStart = ceil(-xRange(1));

colEnd = floor(-xRange(1)+n1);

regIm = Reg(rowStart:rowEnd,colStart:colEnd,:);

%写图像

imwrite(regIm,['name' type '.bmp']);

%显示配准结果图

figure;

subplot(131);imshow(baseIm);title('BaseImage');

subplot(132);imshow(inputIm);title('InputImage');

subplot(133);imshow(regIm);title(['Registration Image using ' type]); end

三、实验结果：

（图一）

（图二）