武汉科技大学图像处理实验报告.docVIP

实验报告要求 每个程序给出图像效果，每个子图像上方标注该图像简要注释。程序代码进行必要的修改完善。 第一次实验报告内容如下： 实验一：图像直方图均衡化与图像滤波 【实验目的】 了解图像直方图的基本数学原理,会编写相应的m函数文件； 掌握matlab图像直方图均衡化、规定化函数.理解图像直方图处理的意义； 掌握噪声模拟和图像滤波函数的使用方法； 理解图像均值滤波、中值滤波的基本原理,掌握Matlab图像滤波基本函数的使用方法。 【实验内容】 自己编写m文件代码，实现求灰度图像的灰度直方图，并与imhist函数的结果比较。 I=imread(rice.bmp); imshow(I) figure,imhist(I,64) figure,imhist(I,128) 自己编写m文件代码，实现某幅图像的3*3的均值滤波。 I=imread(eight.tif); J=imnoise(I,salt pepper,0.02); %对指定的图像加入椒盐噪声 subplot(1,2,1),imshow(J);title(含噪声的图像); K1=filter2(fspecial(average,3),J)/255; %进行3*3模板的均值滤波 subplot(1,2,2),imshow(K1);title(3*3模板的均值滤波后的图像); 读入一张常见的图象处理图像(譬如: Lena, house等图像)，应用IPT的函数对图像作模糊化和添加噪声处理；然后, 参照课本和课堂讲解的例题中的方法对退化的图像进行复原。请写下实验程序代码和必要的代码注释, 整理实验报告. K1=imread(rice.bmp); K2=filter2(fspecial(average,3),K1)/255; %进行3*3模板的均值滤波 K3=imnoise(K2,salt pepper,0.02); %对指定的图像加入椒盐噪声 K4=medfilt2(K3,[3 3]); %进行3*3模板的中值滤波 subplot(2,2,1),imshow(K1);title(原图像); subplot(2,2,2),imshow(K2);title(3*3均值滤波后的图像，模糊化); subplot(2,2,3),imshow(K3);title(加入椒盐噪声后的图像); subplot(2,2,4),imshow(K4);title(进行3*3模板的中值滤波后的图像); 第二次实验报告内容如下： 实验二：图像频域变换及应用 【实验目的】 了解傅立叶变换、离散余弦变换及Radon变换在图像处理中的应用； 掌握Matlab图像频域变换的基本函数的使用； 理解图像频域变换的意义； 【实验内容】 1 绘制一个二值图像矩阵，并将其傅里叶函数可视化。 说明：这个例子将创建一个矩阵，然后用一个二进制图像显示。再对这个图像矩阵做傅里叶变换，显示其幅值。 %创建矩阵 f=zeros(40,40); f(5:24,13:17)=1; subplot(1,2,1),imshow(f);title(矩阵的显示); %进行图像f的傅里叶变换 F=fft2(f); %对幅值取对数 F2=log(abs(F)); subplot(1,2,2),imshow(F2,[-1 5],notruesize);title(傅里叶变换显示); 2 调试如下程序。分析该程序的功能，写上必要注释。 figure(1); f=imread(room.bmp); f=rgb2gray(f); subplot(2,2,1),imshow(f); title(原图像); F=fft2(f); S=fftshift(log(1+abs(F))); S=gscale(S); subplot(2,2,2),imshow(S); title(图像1); h=fspecial(sobel); PQ=paddedsize(size(f)); H=freqz2(h,PQ(1),PQ(2)); H1=double(ifftshift(H)); subplot(2,2,3),imshow(abs(H),[]); title(图像2); subplot(2,2,4),imshow(abs(H1),[]); title(图像3); gs=imfilter(double(f),h); %generate the filtered images in the spatial domain. gf=dftfilt(f,H1); %generate the filtered images in the spatial domain. figure(2); subplot(3,2,1),imshow(gs,[]); t