1 cl;
 2 H=1;                        %索引pix中第一个元素，即高度
 3 W=2;                        %索引pix中第二个元素，即宽度
 4 jiaodu=45;                       %要旋转的角度，旋转方向为顺时针
 5 img=imread(\'Corner.png\');       %这里v为原图像的高度，u为原图像的宽度
 6 imshow(img);                    %这里y为变换后图像的高度，x为变换后图像的宽度
 7 [v u]=size(img);
 8 
 9 theta=jiaodu/180*pi;
10 rot=[cos(theta) -sin(theta) 0;sin(theta) cos(theta) 0;0 0 1]; 
11 inv_rot=inv(rot);
12 pix1=[1 1 1]*rot;               %变换后图像左上点的坐标
13 pix2=[1 u 1]*rot;               %变换后图像右上点的坐标
14 pix3=[v 1 1]*rot;               %变换后图像左下点的坐标
15 pix4=[v u 1]*rot;               %变换后图像右下点的坐标
16 
17 height=round(max([abs(pix1(H)-pix4(H))+0.5 abs(pix2(H)-pix3(H))+0.5]));     %变换后图像的高度
18 width=round(max([abs(pix1(W)-pix4(W))+0.5 abs(pix2(W)-pix3(W))+0.5]));      %变换后图像的宽度
19 imgn=zeros(height,width);
20 
21 delta_y=abs(min([pix1(H) pix2(H) pix3(H) pix4(H)]));            %取得y方向的负轴超出的偏移量
22 delta_x=abs(min([pix1(W) pix2(W) pix3(W) pix4(W)]));            %取得x方向的负轴超出的偏移量
23 
24 for y=1-delta_y:height-delta_y
25     for x=1-delta_x:width-delta_x
26         pix=[y x 1]*inv_rot;                                %用变换后图像的点的坐标去寻找原图像点的坐标，                                         
27                                                             %否则有些变换后的图像的像素点无法完全填充
28         if pix(H)>=0.5 && pix(W)>=0.5 && pix(H)<=v && pix(W)<=u
29             imgn(y+delta_y,x+delta_x)=img(round(pix(H)),round(pix(W)));
30         end   
31         
32     end
33 end
34 
35 figure,imshow(uint8(imgn))