直方图均衡化及空域滤波

1、2012-2013学年第一学期医学图像处理实验报告班级学号姓名实验时间2013-04实验地点4-601实验成绩实验题目直方图均衡化及空域滤波一、实验目的1掌握灰度直方图的概念及直力图均衡化过程；2掌握空域滤波的基本原理及方法；3掌握用MATLAB语言进行图像的空域滤波的方法。二、实验内容1、读入boats720x576.tif图像，用histeq函数将原始图像的灰度直方图均衡化，同时观察均衡化后的图像与前面图像的差别，均衡化后的灰度直方图与前面的灰度直方图的区别。2、读取一幅灰度图像。给读出的图像加入椒盐噪声，采用大小相同的模板对加有噪声的图像分别进行均值滤波和中值滤波，用一个图像处理对话框显

2、示原图像、加有噪声的图像、均值滤波的图像及中值滤波的图像。比较结果。3、给灰度图像加入高斯噪声，重复均值滤波和中值滤波，比较结果。三、实验环境1计算机； 2MATLAB程序；四、实验过程及结果1、实验步骤（列出具体每一步）（1）、读入boats720x576.tif图像，用histeq函数将原始图像的灰度直方图均衡化，同时观察均衡化后的图像与前面图像的差别，均衡化后的灰度直方图与前面的灰度直方图的区别。>> I=imread('pout.tif'); % 读入原图像>> J=histeq(I); %对原图像进行直方图均衡化处理>> %对直方图

3、均衡化后图像进行屏幕控制；作一幅子图，并排两幅图的第1幅>> subplot(2,2,1);imshow(I); %显示原图像>> title('原图像'); %给原图像加标题名>> subplot(2,2,2);imshow(J); >> title('直方图均衡化后的图像') ; >> subplot(2,2,3); >> imhist(I,64); %将原图像直方图显示为64级灰度>> title('原图像直方图') ; %给原图像直方图加标题名>&g

4、t; subplot(2,2,4); >> imhist(J,64) ; %将均衡化后图像的直方图显示为64级灰度>> title('均衡变换后的直方图') ; %给均衡化后图像直方图加标题名（2）、读取一幅灰度图像。给读出的图像加入椒盐噪声，采用大小相同的模板对加有噪声的图像分别进行均值滤波和中值滤波，用一个图像处理对话框显示原图像、加有噪声的图像、均值滤波的图像及中值滤波的图像。比较结果。>> I=imread('pout.tif');>> J=imnoise(I, 'salt & pepper

5、', 0.02);>> h=fspecial('average', 3);>> I2=filter2(h, J);>> I3=medfilt2(J, 3,3);>> subplot(3,3,1), imshow(I), title('原图像');>> subplot(3,3,5), imshow(J), title('加噪声图像');>> subplot(3,3,7), imshow(I2, ), title('加均值滤波后图像');>>

6、subplot(3,3,6), imshow(I3), title('加中值滤波后图像');（3）、给灰度图像加入高斯噪声，重复均值滤波和中值滤波，比较结果。>> I=imread('pout.tif');>> J=imnoise(I, 'gaussian', 0.02);>> I2=medfilt2(J, 10,10);>> h=fspecial('average', 10);>> I3=filter2(h, J);>> subplot(2,2,1), imshow(I), title('原图像');>> subplot(2,2,2), imshow(J), title('加噪声图像');>> subplot(2,2,4), imshow(I3, ), title('均值滤波后图像');>> subplot(2,2,3), imshow(I2), title('中值滤波后图像');2、遇到的主要问题及其解决方案 编程时容易把公式记混淆了3、实验结果(1)、