% 直角坐标系实现图像旋转 im=rgb2gray(imread(\'robot.jpg\')); figure,imshow(im); [row,col]=size(im); a=pi/6; % 逆时针旋转30度 %得到行的初始矩阵 r1=repmat(1:row,row,1); i=r1\'; %得到列的初始矩阵 j=repmat(1:col,col,1); i1=zeros(row,col); for m=1:row for n=1:col i1(m,n)=round(m*cos(a)-n*sin(a)); %获得旋转后的i像素值 并且取整 end end j1=zeros(row,col); for m=1:row for n=1:col j1(m,n)=round(m*sin(a)+n*cos(a)); %获得旋转后的j像素值 并且取整 end end %判断旋转后的i,j是否小于0 坐标调整 if min(min(i1))<min(min(j1)) if min(min(i1))<0 k=abs(min(min(i1))); new_min=min(min(i1))+(k+1); A=i1+ones(row,col)*(k+1); end elseif min(min(j1))<0 k=abs(min(min(j1))); new_min=min(min(j1))+(k+1); A=j1+ones(row,col)*(k+1); end new_max_i=max(max(A)); %最大行 new_max_j=max(max(j1)); %最大列 G=zeros(new_max_i,new_max_j); %构造空白的新画布 %接下来 进项数据填充 for m=1:row for n=1:col G(A(m,n),j1(m,n))=im(m,n); end end figure,imshow(uint8(G));