图像配准matlab代码

1、registration.mclearclcimg = imread(.imagetest1.tif);figure(1)imshow(img)disp(获取原图像四个顶点在该图中的坐标);impixelinfo%四个顶点的坐标x1 = 0;y1 = 0;x2 = 688;y2 = 200;x3 = 20;y3 = 688;x4 = 708;y4 = 888;rs = size(img);rows = rs(1); %行数columns = rs(2); % 列数%rows , columns = size(img);%获取原图的大小R = round(y3-y1)+(y4-y2)./2);

2、%行C = round(x2-x1)+(x4-x3)./2); %列D1 = rows -R; %垂直偏移像素数D2 = columns -C; %水平偏移像素数img1 = zeros(R,C);for i = 1:Rd2 = round(D2*i./R);for j = 1:Cd1 = round(D1*j./C);img1(i,j) = floor(img(i+d1,j+d2);endendfigure(2)imshow(img1)imwrite(uint8(img1),strcat(.image,result.tif),tif);transform.mclearclc%读取图像img

3、= imread(.imagetest.tif);figure(1)imshow(img)%计算出图像的大小rs = size(img);rows = rs(1); %行数columns = rs(2); % 列数%将图像水平偏移20个像素，垂直偏移200个像素img1 = zeros(rows+200,columns+20);for i = 1:rowsx = round(20 * i ./ rows);for j = 1:columnsy = round(200 * j ./columns);img1(i+y,j+x) = floor(img(i,j);endendimwrite(uint

4、8(img1),strcat(.image,test1.tif),tif);figure(2)imshow(img1)WORK.mclearclcimg = imread(.imagetest1.tif);org = imread(.imagetest.tif);figure(name,原图)imshow(org)impixelinfofigure(name,偏移之后的图像)imshow(img)disp(获取原图像四个顶点在该图中的坐标);impixelinfor,c = size(org);u1 = 619;v1 = 150;u2 = 574;v2 = 605;u3 = 130;v3 =

5、234;u4 = 58;v4 = 538;x1 = 622;y1 = 329;x2 = 590;y2 = 773;x3 = 139;y3 = 275;x4 = 75;y4 = 557;A = u1,v1,u1*v1,1,0,0,0,0;.0,0,0,0,u1,v1,u1*v1,1;.u2,v2,u2*v2,1,0,0,0,0;.0,0,0,0,u2,v2,u2*v2,1;.u3,v3,u3*v3,1,0,0,0,0;.0,0,0,0,u3,v3,u3*v3,1;.u4,v4,u4*v4,1,0,0,0,0;.0,0,0,0,u4,v4,u4*v4,1;B = x1;y1;x2;y2;x3;y3;x4;y4;%C = c1;c2;c3;c4;c5;c6;c7;c8;C = inv(A)*B;%c1,c2,c3,c4,c5,c6,c7,c8 = C;c1 = C(1,1);c2 = C(2,1);c3 = C(3,1);c4 = C(4,1);c5 = C(5,1);c6 = C(6,1);c7 = C(7,1);c8 = C(8,1);I = zeros(r,c);rr,cc = size(img);for i =1:rfor j = 1:cx = round(c1*i+c2*j+c3*i*j+c4);y = round(c5*i+c6*j+c7*i*j+