1 close all;
 2 clc;
 3 H=1;                        %索引pix中第一個元素，即高度
 4 W=2;                        %索引pix中第二個元素，即寬度
 5 left_right=0.3;               %抬起左邊或右邊時值為0-1之間，不抬起時為0
 6 up_down=0;                %抬起上邊或下邊時值為0-1之間，不抬起時為0
 7 
 8 img=imread('lena.jpg');       %這裡v為原影象的高度，u為原影象的寬度
 9 imshow(img);                    %這裡y為變換後圖像的高度，x為變換後圖像的寬度
 
10 %img=flipud(img);           %註釋掉為抬起下邊，沒註釋掉為抬起上邊
11 %img=fliplr(img);           %註釋掉為抬起右邊，沒註釋掉為抬起左邊
12 [v u]=size(img);
13 
14 
15 a=1;b=up_down;c=0;
16 d=left_right;e=1;f=0;
17 g=up_down/v;h=left_right/u;i=1;
18 rot=[a b c;d e f;g h i];
19 
20 pix1=[1 1 1]*rot./(g+h+i);                 %變換後圖像左上點的座標
21 pix2=[1
 
 u 1]*rot./(g*v+h+i);               %變換後圖像右上點的座標
22 pix3=[v 1 1]*rot./(g+h*u+i);               %變換後圖像左下點的座標
23 pix4=[v u 1]*rot./(g*v+h*u+i);             %變換後圖像右下點的座標
24 
25 height=round(max([abs(pix1(H)-pix3(H))+0.5 abs(pix2(H)-pix3(H))+0.5 ...
26                   abs(pix1(H)-pix4(H))+0.5 abs(pix2(H)-pix4(H))+0.5
 
]));     %變換後圖像的高度
27 
28 width=round(max([abs(pix1(W)-pix2(W))+0.5 abs(pix3(W)-pix2(W))+0.5 ...
29                  abs(pix1(W)-pix4(W))+0.5 abs(pix3(W)-pix4(W))+0.5]));      %變換後圖像的寬度
30 imgn=zeros(height,width);
31 
32 delta_y=abs(min([pix1(H)-0.5 pix2(H)-0.5 pix3(H)-0.5 pix4(H)-0.5]));            %取得y方向的負軸超出的偏移量
33 delta_x=abs(min([pix1(W)-0.5 pix2(W)-0.5 pix3(W)-0.5 pix4(W)-0.5]));            %取得x方向的負軸超出的偏移量
34 
35 for y=1-floor(delta_y):height-floor(delta_y)
36     for x=1-floor(delta_x):width-floor(delta_x)
37         pix=[y x 1]/rot*(g*y+h*x+i);                                %用變換後圖像的點的座標去尋找原影象點的座標，                                         
38                                                             %否則有些變換後的影象的畫素點無法完全填充
39         if pix(H)>=0.5 && pix(W)>=0.5 && pix(H)<=v && pix(W)<=u
40             imgn(y+floor(delta_y),x+floor(delta_x))=img(round(pix(H)),round(pix(W)));
41         end   
42         
43     end
44 end
45 figure,imshow(uint8(imgn));
46 
47 %%二次變換,暫時只想起來這種笨方法，一次變換好像沒法進行透視變換
48 img=imgn;
49 [v u]=size(img);
50 a=1;b=-b/2;c=0;
51 d=-d/2;e=1;f=0;
52 g=0;h=0;i=1;
53 rot=[a b c;d e f;g h i];
54 
55 pix1=[1 1 1]*rot./(g+h+i);                 %變換後圖像左上點的座標
56 pix2=[1 u 1]*rot./(g*v+h+i);               %變換後圖像右上點的座標
57 pix3=[v 1 1]*rot./(g+h*u+i);               %變換後圖像左下點的座標
58 pix4=[v u 1]*rot./(g*v+h*u+i);             %變換後圖像右下點的座標
59 
60 height=round(max([abs(pix1(H)-pix3(H))+0.5 abs(pix2(H)-pix3(H))+0.5 ...
61                   abs(pix1(H)-pix4(H))+0.5 abs(pix2(H)-pix4(H))+0.5]));     %變換後圖像的高度
62 
63 width=round(max([abs(pix1(W)-pix2(W))+0.5 abs(pix3(W)-pix2(W))+0.5 ...
64                  abs(pix1(W)-pix4(W))+0.5 abs(pix3(W)-pix4(W))+0.5]));      %變換後圖像的寬度
65 imgn=zeros(height,width);
66 
67 delta_y=abs(min([pix1(H)-0.5 pix2(H)-0.5 pix3(H)-0.5 pix4(H)-0.5]));            %取得y方向的負軸超出的偏移量
68 delta_x=abs(min([pix1(W)-0.5 pix2(W)-0.5 pix3(W)-0.5 pix4(W)-0.5]));            %取得x方向的負軸超出的偏移量
69 
70 for y=1-floor(delta_y):height-floor(delta_y)
71     for x=1-floor(delta_x):width-floor(delta_x)
72         pix=[y x 1]/rot*(g*y+h*x+i);                                %用變換後圖像的點的座標去尋找原影象點的座標，                                         
73                                                             %否則有些變換後的影象的畫素點無法完全填充
74         if pix(H)>=0.5 && pix(W)>=0.5 && pix(H)<=v && pix(W)<=u
75             imgn(y+floor(delta_y),x+floor(delta_x))=img(round(pix(H)),round(pix(W)));
76         end   
77         
78     end
79 end
80 %imgn=flipud(imgn);             %註釋掉為抬起下邊，沒註釋掉為抬起上邊
81 %imgn=fliplr(imgn);             %註釋掉為抬起右邊，沒註釋掉為抬起左邊
82 figure,imshow(uint8(imgn));

圖片福利：