数字图像 代码.docx

function?figure_rotate%?构造原始图像I=zeros(256,256);I(88:168,124:132)=1;%?subplot(1,2,1)imshow(I)%?求原始图像的频谱J=fft2(I);F=abs(J);J1=fftshift(F);?figure%subplot(1,2,2)imshow(J1,[5?50])?%?灰度级低于5的显示黑色，高于50的显示白色，在[5?50]范围内的保持其灰度级不变%imshow(J1)%?对一副图像进行旋转J=imrotate(I,90,bilinear,crop);figureimshow(J)%?求旋转后图像的频谱J1=fft2(J);F=abs(J1);J2=fftshift(F);figureimshow(J2,[5?50])2015/5/7?11:30:04曾曙?2015/5/7?11:30:04function?zhifangtuclear?all?close?all%?I=imread(cameraman.tif);????%??lena.png%?I=imread(lena.jpg);????%??lena.pngI=imread(car0.jpg);????%??lena.pngI=rgb2gray(I);??%?将RGB图像转换为灰度图像imshow(I)figure(2)imhist(I,64)?%?显示一副图像的直方图，只能用于灰度图像J=histeq(I);?%?直方图均衡化figure(3)imshow(J)figure(4)imhist(J,64)2015/5/7?11:33:07曾曙?2015/5/7?11:33:07I=rgb2gray(I);??%?将RGB图像转换为灰度图像2015/5/7?11:48:051班-黄?2015/5/7?11:48:052015/5/8?15:42:50曾曙?2015/5/8?15:42:50由于L2621机位有限，周四上午?3~4节（10~14周）在L2621的上机改为?周二上午?3~4节（11~14周）在L2616上。同时周四上午?3~4节（11~14周）的理论课照常在J3216上。请相互转告下。曾曙?2015/5/8?15:43:28由于L2621机位有限，周四上午?3~4节（10~14周）在L2621的上机改为?周二上午?3~4节（11~14周）在L2616上。同时周四上午?3~4节（11~14周）的理论课照常在J3216上。请相互转告下。2015/5/8?18:25:16落花雨~?2015/5/8?18:25:162015/5/12?11:00:54曾曙?2015/5/12?11:00:54function?lvbo%?采用领域平均、中值滤波等方法对一副添加了噪声的图像进行滤波处理clear?allclose?allI=imread(pout.tif);%?subplot(2,2,1);imshow(I);%?xlabel((a)原始图像);J0=imnoise(I,salt??pepper,0.05);?%添加椒盐噪声J=imnoise(I,gaussian,0,0.025);?%添加高斯噪声subplot(2,2,1);imshow(J);xlabel(（a）添加高斯噪声);J1=filter2(fspecial(average,3),J)/255;?%3*3模板均值滤波subplot(2,2,2);imshow(J1);xlabel((b)3*3模板均值滤波去噪);J2=filter2(fspecial(average,5),J)/255;?%5*5模板均值滤波subplot(2,2,3);imshow(J2);xlabel((c)5*5模板均值滤波去噪);J3=filter2(fspecial(average,7),J)/255;?%7*7模板均值滤波subplot(2,2,4);imshow(J3);xlabel((d)7*7模板均值滤波去噪);figure(2)subplot(2,2,1);imshow(J0);xlabel((a)添加椒盐噪声);subplot(2,2,2);imshow(J);xlabel((b)添加高斯噪声);J4=medfilt2(J0,[5,5]);??%?添加椒盐噪声图像进行5×5方形窗口中值滤波J5=medfilt2(J,[5,5]);??%?添加高斯噪声图像进行5×5方形窗口中值滤波subplot(2,2,3);imshow(J4);xlabel((c)椒盐噪声图像中值滤波);subplot(2,2,4);imshow(J5);xlabel((d)高斯噪声图像中值滤波);2015/5/12?11:17:21曾曙?2015/5/12?11