xuanzhuan.rar_MATLAB 图像旋转_图像旋转_图像转换_旋转_灰度图像转换

a=imread('d:/1.jpg'); a=rgb2gray(a); a=double(a); sizea=size(a); m=sizea(1,1);%得到行数 n=sizea(1,2);%得到列数 b=a; angle=input('degree: '); angle=angle*pi/180;%弧度化 i=zeros(m,n); j=i; for p=1:m i(p,:)=p; end for p=1:n j(:,p)=p; end%i是行坐标，j是列坐标 it=i*cos(angle)-j*sin(angle); jt=i*sin(angle)+j*cos(angle);%旋转公式 it=round(it); jt=round(jt);%取整 imin=min(min(it)); jmin=min(min(jt)); if (imin<=0) it=it+(-imin+1); end if (jmin<=0) jt=jt+(-jmin+1); end%画布扩大，因为坐标只能是正整数 imax=max(max(it)); jmax=max(max(jt)); %b=ones(imax,jmax); for p=1:imax for q=1:jmax b(p,q)=255;%将新矩阵全填充为255 end end for p=1:m for q=1:n b(it(p,q),jt(p,q))=a(p,q);%变换到新图像 end end for p=2:(imax-1) for q=2:(jmax-1) if( b(p,q)==255 ) b(p,q)=round((double(b(p-1,q))+b(p+1,q)+b(p,q+1)+b(p,q-1))/4); end end end a=uint8(a); b=uint8(b); subplot(1,2,1);imshow(a);title('原图像'); subplot(1,2,2);imshow(b);title('处理后图像');