資源描述：
《圖像空域增強.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。
1、本科生實驗報告（二）姓名:學(xué)院:專業(yè):信息安全班級:信息安全實驗課程名稱:數(shù)字圖像處理與分析實驗日期:2013年4月25日指導(dǎo)教師及職稱:王云峰實驗成績:開課時間：2012-2013學(xué)年第2學(xué)期甘肅政法學(xué)院實驗管理中心印制實驗題目圖像空域增強小組合作否姓名班級信息安全學(xué)號一、實驗?zāi)康?、掌握灰度變換的基本原理。2、掌握直方圖處理的基本原理。3、掌握Matlab中灰度變換和直方圖處理的實現(xiàn)方法。4、理解圖像濾波的基本定義及目的；5、掌握空域濾波的基本原理及方法；6、掌握用MATLAB語言進行圖像的空域濾波的方法。二．實驗環(huán)
2、境裝有的MATLAB軟件的計算機和一些圖片。三、實驗內(nèi)容與步驟圖象增強是通過某種技術(shù)有選擇地突出對某一具體應(yīng)用有用的信息，削弱或抑制一些無用的信息。圖象增強按增強處理所在空間不同分為空域增強方法和頻域增強方法?？沼蛟鰪?直接在圖像所在的二維空間進行處理，即直接對每一像素的灰度值進行處理?？臻g域增強方法主要有灰度變換增強、直方圖增強、圖像平滑和圖像銳化等。1.灰度變換灰度變換:基于點操作，將每一個像素的灰度值按照一定的數(shù)學(xué)變換公式轉(zhuǎn)換為一個新的灰度值。常用的有:對比度增強、直方圖均衡化等方法。分段線性灰度變換(1)程序源代
3、碼：closeallclearclcdisp('====E4_4_2.m====');I=imread('rice.png');subplot(3,3,1),imshow(I),ylabel('原圖像');subplot(3,3,2),imhist(I);%方法1-系統(tǒng)函數(shù)%把灰度值范圍[80,160]映射到[20,220]，灰度值小于80的映射為20，灰度值大于160的映射為220J=imadjust(I,[80/255160/255],[20/255220/255]);subplot(3,3,4),imshow(J)
4、,ylabel('變換圖像(方法1)');subplot(3,3,5),imhist(J);%方法2-編程實現(xiàn)%把灰度值范圍[80,160]映射到[20,220]，灰度值小于80的映射為20，灰度值大于160的映射為220f0=0;g0=20;%分段曲線的第1個點f1=80;g1=20;%分段曲線的第2個點f2=160;g2=220;%分段曲線的第3個點f3=255;g3=220;%分段曲線的第4個點subplot(3,3,9),plot([f0,f1,f2,f3],[g0,g1,g2,g3]),xlabel('f'),
5、ylabel('g'),axis([02550255]);%繪制變換曲線r1=(g1-g0)/(f1-f0);%曲線1的斜率b1=g0-r1*f0;%曲線1的截距r2=(g2-g1)/(f2-f1);%曲線2的斜率b2=g1-r2*f1;%曲線2的截距r3=(g3-g2)/(f3-f2);%曲線3的斜率b3=g2-r3*f2;%曲線3的截距[m,n]=size(I);K=double(I);fori=1:mforj=1:nf=K(i,j);g(i,j)=0;17if(f>=f0)&(f<=f1)g(i,j)=r1*f+b
6、1;%曲線1的方程y=r1*x+b1elseif(f>=f1)&(f<=f2)g(i,j)=r2*f+b2;%曲線2的方程y=r2*x+b2elseif(f>=f2)&(f<=f3)g(i,j)=r3*f+b3;%曲線3的方程y=r3*x+b3endendendendendsubplot(3,3,7),imshow(uint8(g)),ylabel('變換圖像(方法2)');subplot(3,3,8),imhist(uint8(g));(2)實驗結(jié)果灰度變換（對數(shù)變換）(1)程序源代碼：closeallclearclc
7、disp('====E4_4_6.m====');I=imread('rice.png');subplot(3,3,1),imshow(I),ylabel('原圖像');subplot(3,3,2),imhist(I);%對數(shù)變換c=255/log(1+255);%語句1x=0:1:255;y=c*log(1+x);subplot(3,3,6),plot(x,y),axistight,xlabel('f'),ylabel('g');%繪制變換曲線[m,n]=size(I);K=double(I);g=c*log(K+1)
8、;%{%等價的循環(huán)語句實現(xiàn)fori=1:mforj=1:ng(i,j)=0;g(i,j)=c*log(K(i,j)+1);endend%}subplot(3,3,4),imshow(uint8(g)),ylabel('對數(shù)變換');subplot(3,3,5),imhist(uint8(g));%###系統(tǒng)函數(shù)調(diào)用#
