要求
图像是彩色的;旋转后,原图像的四个角点应在新图像的四条边上(即新图像的大小是包含原图像旋转后的最小矩形);旋转的角度可以是 0~ 360度;函数的调用形式是 J=myrotate(I,angle);不允许使用 imrotate。
原理
原理部分仅仅是代码中所体现的,具体推导过程参看 百度文库 图像旋转
下图显示新图像与原图像的长宽关系,注意如果旋转角度是120度,则需要取绝对值,具体参见下面源代码。
下图显示原图像旋转后的点与原图像点的对应关系,以按此公式求得新图像每一点对应的原图像点。
旋转后,由于图像点可能是非整数,出现像素空洞问题,故需要进行插值法,下面选用最近邻插值法。
源代码
function img_rotate=my_img_rotate_chazhi(img1,angle)
% 旋转后图像=my_img_rotate(原图像,旋转角度) 0<旋转角度<360
% 查找新图对应的原图像素点
% 最邻近点插值方法
[h,w,d]=size(img1);
radian=angle/180*pi;
cos_val = cos(radian);
sin_val = sin(radian);
w2=round(abs(cos_val)*w+h*abs(sin_val));
h2=round(abs(cos_val)*h+w*abs(sin_val));
img_rotate = uint8(zeros(h2,w2,3)); %像素是整数
for x=1:w2
for y=1:h2
x0 = uint32(x*cos_val + y*sin_val -0.5*w2*cos_val-0.5*h2*sin_val+0.5*w);
y0= uint32(y*cos_val-x*sin_val+0.5*w2*sin_val-0.5*h2*cos_val+0.5*h);
x0=round(x0); %最邻近插值
y0=round(y0); %最邻近插值
if x0>0 && y0>0&& w >= x0&& h >= y0
img_rotate(y,x,:) = img1(y0,x0,:);
end
end
end
%I = imread('C:\Users\Desktop\love.jpg')
%I2=my_img_rotate(I,30);
%figure,imshow(I2);