新·MATLAB图像几何处理.docVIP

第16周上机实验报告 学号：912104520120 姓名：鲍春霖 1、 内容 （1）写出以图像中心点作为支点旋转任意角度，图像仍保持完整的旋转程序。 （2）写出N*N模板的对称近邻平滑滤波算法的具体实现程序。 （3）PS上机运行旋转操作。 2、 原理clear all; A=imread(D:\cameraman.tif) [row,col]=size(A); B=A; C=A; for x=1:row for y=1:col B(x,y)=A(x,col+1-y); end end a=(-pi/6); xd=row*cos(a)-col*sin(a); yd=col*cos(a)-row*sin(a); for i=1:xd for j=1:yd m0=(i*cos(a)+(j+col*sin(a))*sin(a)); n0=(-i*sin(a)+(j+col*sin(a))*cos(a)); m=floor(m0); n=floor(n0); p=m0-m; q=n0-n; C(i,j)=0; if m=1n=1m256n256 C(i,j)=(1-q)*((1-p)*B(m,n)+p*B(m+1,n))+q*((1-p)*B(m,n+1)+p*B(m+1,n+1)); end end end subplot(1,2,1),imshow(A),title(?-í?); subplot(1,2,2),imshow(C),title(Dy×a?μ??oóí???); 我们通过老师给予的其中一个赋值程序，得到了在顺时针时候成立的一个程序。顺时针的时候，画布的扩展需要纵坐标的改变。而这个程序并不适用于旋转角度大于零，即逆时针的旋转，为此做出了一个优化，使得程序变得更加完整。 C(i,j)=(1-q)*((1-p)*B(m,n)+p*B(m+1,n))+q*((1-p)*B(m,n+1)+p*B(m+1,n+1)); 该语句是一个插值语句，语句作用是为了提升旋转后的画面质量。 （2）对称近邻平滑滤波 本程序作用在于做出一个N*N模板的对称近邻平滑滤波算法的具体实现程序。对称算法的要求就是当N为级数时候，有（N*N-1）/2对像素点，然后通过这些像素点与变换点的比较，选择这些像素点对中大小接近需要变换点的数值，然后对这些选择出来的（N*N-1）/2个像素点的灰度值求平均值，从而得出来需变换点的像素点灰度值。 3、 方法及程序 （1）旋转程序 clear all; A=imread(D:\cameraman.tif) [row,col]=size(A); B=A; C=A; for x=1:row for y=1:col B(x,y)=A(x,col+1-y); end end a=(-pi/6); if (a=0) xd=row*cos(a)-col*sin(a); yd=col*cos(a)-row*sin(a); else xd=row*cos(a)+col*sin(a); yd=col*cos(a)+row*sin(a); end for i=1:xd for j=1:yd if (a=0) m0=(i*cos(a)+(j+col*sin(a))*sin(a)); n0=(-i*sin(a)+(j+col*sin(a))*cos(a)); else m0=((i-row*sin(a))*cos(a)+(j)*sin(a)); n0=(-(i-row*sin(a))*sin(a)+(j)*cos(a)); end m=floor(m0); n=floor(n0); p=m0-m; q=n0-n; C(i,j)=0; if m=1n=1m256n256 C(i,j)=(1-q)*((1-p)*B(m,n)+p*B(m+1,n))+q*((1-p)*B(m,n+1)+p*B(m+1,n+1));