直方圖均衡化及空域濾波

1、2012-2013學年第一學期醫(yī)學圖像處理實驗報告班級學號姓名實驗時間2013-04實驗地點4-601實驗成績實驗題目直方圖均衡化及空域濾波一、實驗?zāi)康?掌握灰度直方圖的概念及直力圖均衡化過程；2掌握空域濾波的基本原理及方法；3掌握用MATLAB語言進行圖像的空域濾波的方法。二、實驗內(nèi)容1、讀入boats720x576.tif圖像，用histeq函數(shù)將原始圖像的灰度直方圖均衡化，同時觀察均衡化后的圖像與前面圖像的差別，均衡化后的灰度直方圖與前面的灰度直方圖的區(qū)別。2、讀取一幅灰度圖像。給讀出的圖像加入椒鹽噪聲，采用大小相同的模板對加有噪聲的圖像分別進行均值濾波和中值濾波，用一個圖像處理對話框顯

2、示原圖像、加有噪聲的圖像、均值濾波的圖像及中值濾波的圖像。比較結(jié)果。3、給灰度圖像加入高斯噪聲，重復(fù)均值濾波和中值濾波，比較結(jié)果。三、實驗環(huán)境1計算機； 2MATLAB程序；四、實驗過程及結(jié)果1、實驗步驟（列出具體每一步）（1）、讀入boats720x576.tif圖像，用histeq函數(shù)將原始圖像的灰度直方圖均衡化，同時觀察均衡化后的圖像與前面圖像的差別，均衡化后的灰度直方圖與前面的灰度直方圖的區(qū)別。>> 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）、讀取一幅灰度圖像。給讀出的圖像加入椒鹽噪聲，采用大小相同的模板對加有噪聲的圖像分別進行均值濾波和中值濾波，用一個圖像處理對話框顯示原圖像、加有噪聲的圖像、均值濾波的圖像及中值濾波的圖像。比較結(jié)果。>> 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）、給灰度圖像加入高斯噪聲，重復(fù)均值濾波和中值濾波，比較結(jié)果。>> 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、實驗結(jié)果(1)、