透视变换源代码_matlab；图像处理_透视变换_图像变换

%程序参考https://www.cnblogs.com/tiandsp/archive/2012/12/16/2820916.html clear all; close all; clc; %读取并显示原图像。 img= imread('class4-透视变换.jpg'); imshow(img); [N M layers] = size(img); %M:图宽，N:图高。Matlab 的image里面 x在后，y在前！ %取四个点，依次是左上，右上，右下，左下,这里我们取的是牌照的四个角 dot=ginput(); x=[dot(1,1) dot(2,1) dot(3,1) dot(4,1)]; %四个原顶点x坐标 y=[dot(1,2) dot(2,2) dot(3,2) dot(4,2)]; %四个原顶点y坐标 line([x x(1)],[y y(1)],'color','red','linewidth',3); %画出矩形 %设置变换后的新顶点，存在向量X和Y中。 %这里我们做成矩形,也可以做成其它形状。 w=round(sqrt((x(1)-x(2))^2+(y(1)-y(2))^2)); %从原四边形获得新矩形宽 h=round(sqrt((x(1)-x(4))^2+(y(1)-y(4))^2)); %从原四边形获得新矩形高 X=[x(1) x(1)+w x(1)+w x(1)]; %变换后的新顶点x坐标 Y=[y(1) y(1) y(1)+h y(1)+h]; %变换后的新顶点y坐标 %联立解方程组：A Tv = B，其中Tv是以向量形式表示的透视变换。 B=[X(1) Y(1) X(2) Y(2) X(3) Y(3) X(4) Y(4)]'; %变换后的四个顶点，方程组右边的值 A=[x(1) y(1) 1 0 0 0 -X(1)*x(1) -X(1)*y(1); 0 0 0 x(1) y(1) 1 -Y(1)*x(1) -Y(1)*y(1); x(2) y(2) 1 0 0 0 -X(2)*x(2) -X(2)*y(2); 0 0 0 x(2) y(2) 1 -Y(2)*x(2) -Y(2)*y(2); x(3) y(3) 1 0 0 0 -X(3)*x(3) -X(3)*y(3); 0 0 0 x(3) y(3) 1 -Y(3)*x(3) -Y(3)*y(3); x(4) y(4) 1 0 0 0 -X(4)*x(4) -X(4)*y(4); 0 0 0 x(4) y(4) 1 -Y(4)*x(4) -Y(4)*y(4)]; Tv=inv(A)*B; %求方程组的解 Tm=[Tv(1) Tv(2) Tv(3); %将Tv写成矩阵形式Tm。 Tv(4) Tv(5) Tv(6); Tv(7) Tv(8) 1 ]; inv_Tm=inv(Tm); %逆变换，用于计算变换后图像。 %计算变换后的图像大小。 pix1=Tm*[1 1 1]'; pix1 = pix1/pix1(3); %变换后图像左上点 pix2=Tm*[M 1 1]'; pix2 = pix2/pix2(3); %变换后图像右上点 pix3=Tm*[1 N 1]'; pix3 = pix3/pix3(3); %变换后图像左下点 pix4=Tm*[M N 1]'; pix4 = pix4/pix4(3); %变换后图像右下点 min_x=round(min([pix1(1) pix2(1) pix3(1) pix4(1)])); %变换后x坐标最小值 min_y=round(min([pix1(2) pix2(2) pix3(2) pix4(2)])); %变换后y坐标最小值 width=round(max([pix1(1) pix2(1) pix3(1) pix4(1)])-min_x); %变换后图像的宽度 height=round(max([pix1(2) pix2(2) pix3(2) pix4(2)])-min_y); %变换后图像的高度 width = min(5000,width); %转换后的图像可能会非常大，我们必须设置宽度和高度上限。 height= min(5000,height); %计算变换后图像。 imgn=zeros(height,width,layers); %Matlab 的image里面 x在后，y在前！ for i = min_x:width+min_x-1 %从变换图像中反向寻找原图像的点，以免出现空洞，和旋转放大原理一样 for j = min_y:height+min_y-1 pix=inv_Tm*[i j 1]'; %求原图像中坐标. pix=round(pix/pix(3)); if pix(1)>=1 && pix(2)>=1 && pix(1)<=M && pix(2)<=N for k=1:layers %Matlab 的image里面 x在后，y在前！ imgn(j+1-min_y, i+1-min_x,k)=img(pix(2), pix(1),k); %最邻近插值,也可以用双线性或双立方插值 end end end end %显示变换后图像。 figure; imshow(uint8(imgn)); line([X X(1)]+1-min_x,[Y Y(1)]+1-min_y,'color','red','linewidth',3); %画出变换后矩形