matlab練習程序（圖像錯切）

cl;
H=1;                            %索引pix中第一個元素，即高度
W=2;                            %索引pix中第二個元素，即寬度
jiaodu=0;                       %要旋轉的角度，旋轉方向為順時針
img=imread('Corner.png');         %這里v為原圖像的高度，u為原圖像的寬度
imshow(img);                    %這里y為變換后圖像的高度，x為變換后圖像的寬度
[v u]=size(img);
h=v;                            %原圖像的高度
w=u;                            %原圖像的寬度
pro=w/h;                        %寬高比

%雖然這段程序是本人寫的，不過這個仿射變換系數還是不太好表述，暫時是本人現在的理解。
%仿射變換系數,變換后圖像的傾斜強度。取值區間為（0,1)為圖像向右上錯切；取0時不變換圖像；取1時為垂直錯切，不顯圖像；大於1時向左下錯切圖像
%區間為（-1,0）為向左上錯切圖像；為-1時垂直錯切，不顯圖像；小於-1時為右下錯切圖像。 系數過大或過小處理時間會很長，再大/再小就會報錯
aff_cof=0.6;          

theta=jiaodu/180*pi;
rot1=[cos(theta) -sin(theta) 0;sin(theta) cos(theta) 0;0 0 1];      %旋轉變換矩陣
rot2=[1 aff_cof*pro 0;aff_cof*(1/pro) 1 0;0 0 1];                   %錯切變換矩陣
rot=rot2*rot1;                                                      %聯合變換矩陣
inv_rot=inv(rot);
pix1=[1 1 1]*rot;                                                   %變換后圖像左上點的坐標
pix2=[1 u 1]*rot;                                                   %變換后圖像右上點的坐標
pix3=[v 1 1]*rot;                                                   %變換后圖像左下點的坐標
pix4=[v u 1]*rot;                                                   %變換后圖像右下點的坐標

height=round(max([abs(pix1(H)-pix4(H))+0.5 abs(pix2(H)-pix3(H))+0.5]));     %變換后圖像的高度
width=round(max([abs(pix1(W)-pix4(W))+0.5 abs(pix2(W)-pix3(W))+0.5]));      %變換后圖像的寬度
imgn=zeros(height,width);

delta_y=abs(min([pix1(H) pix2(H) pix3(H) pix4(H)]));                        %取得y方向的負軸超出的偏移量
delta_x=abs(min([pix1(W) pix2(W) pix3(W) pix4(W)]));                        %取得x方向的負軸超出的偏移量

for y=1-delta_y:height-delta_y
    for x=1-delta_x:width-delta_x
        pix=[y x 1]*inv_rot;                                                %用變換后圖像的點的坐標去尋找原圖像點的坐標                                        
                                                                            %否則有些變換后的圖像的像素點無法完全填充
        if pix(H)>=0.5 && pix(W)>=0.5 && pix(H)<=v && pix(W)<=u
            imgn(y+delta_y,x+delta_x)=img(round(pix(H)),round(pix(W)));
        end   
        
    end
end

figure,imshow(uint8(imgn))