圖像配準(zhǔn)matlab代碼

registration.mclearclcimg=imread('.\image\test1.tif');figure(1)imshow(img)disp(impixelinfo%四個(gè)頂點(diǎn)的坐標(biāo)x1=0;y1=0;x2=688;y2=200;x3=20;y3=688;x4=708;y4=888;rs=size(img);rows=rs(1); %行數(shù)columns=rs(2); %列%[rows,columns]=size(img);%獲取原圖的大小R=round(((y3-y1)+(y4-y2))./2);%行C=round(((x2-x1)+(x4-x3))./2);%列D1=rows-R;%垂直偏移像素?cái)?shù)D2=columns-C;%水平偏移像素?cái)?shù)img1=zeros(R,C);fori=1:Rd2=round(D2*i./R);forj=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=imread('.\image\test.tif');figure(1)imshow(img)%計(jì)算出圖像的大小rs=size(img);rows=rs(1); %行數(shù)columns=rs(2); %列20200img1=zeros(rows+200,columns+20);fori=1:rowsx=round(20*i./rows);forj=1:columnsy=round(200*j./columns);img1(i+y,j+x)=floor(img(i,j));endendimwrite(uint8(img1),strcat('.\image\','test1.tif'),'tif');figure(2)imshow(img1)WORK.mclearclcimg= org=imread('.\image\test.tif');figure('name',imshow(org)impixelinfofigure('name',imshow(img)disp(impixelinfo[r,c]=size(org);u1=619;v1=150;u2=574;v2=605;u3=130;v3=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=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);fori=1:rforj=1:cx=round(c1*i+c2*j+c3*i*j+c4);y=round(c5*i+c6*j+c7*i*j+c8);ifx>rre