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))
