資源描述:
《空域圖像增強(qiáng)報(bào)告》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。
1、UniverseofEledronicSoenoeandTedinclcgyofChina實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程:光電圖像外理名:學(xué)號:實(shí)驗(yàn)地點(diǎn):指導(dǎo)老師:實(shí)驗(yàn)時(shí)間:一.實(shí)驗(yàn)名稱:空間圖象增強(qiáng)(一)一.實(shí)驗(yàn)?zāi)康?.熟悉和掌握利用matlabI具進(jìn)行數(shù)字圖像的讀、寫、顯示、像素處理等數(shù)字圖像處理的基本步驟和流程。2.熟練掌握各種空間域圖像增強(qiáng)的基木原理及方法。3.熟悉通過灰度變換方式進(jìn)行圖像增強(qiáng)的基本原理、方法和實(shí)現(xiàn)。4.熟悉直方圖均衡化的基本原理、方法和實(shí)現(xiàn)。二.實(shí)驗(yàn)原理(一)數(shù)字圖像的灰度變換1、線性灰度變換令原圖像灰度范圍由[a,b]線性變換后圖像灰度范
2、圍[c,d],以便于提升原圖像的對比度。線性灰度變換公式如下:df(x,y)>bg(x,y)=+ca(x,y)5b(1)b-aC/(X,刃Va其中,./(x,y)為原始圖像,灰度范圍為[a,b],g(x,y)為增強(qiáng)后的數(shù)字圖像,灰度范圍為[C,d]O注:實(shí)驗(yàn)中[a,b]可由實(shí)際圖像來確定,[c,d]可視具體情況人為給定。2、非線性灰度變換(對數(shù)變換或幕律變換選作其一)非線性拉仲不是對圖像的灰度值進(jìn)行擴(kuò)展,而是有選擇地對某一灰度范圍進(jìn)行擴(kuò)展,其他范圍的灰度值則可能被壓縮。常用的非線性變換:對數(shù)變換和指數(shù)變換。即對數(shù)變換公式為:g(x,y)=clog(
3、l+/(兀,y))幕律(伽馬)變換為:g(x,y)=c/(x,y)r(二)直方圖處理對圖像2種定義下的直方圖統(tǒng)計(jì)統(tǒng)計(jì),并分別畫出兩種定義下的直方圖。直方圖均衡化主要用于増強(qiáng)動態(tài)范圍偏小的圖像的反差。該方法的基本思想是把原始的直方圖變換為均勻分布的形狀,這樣就增加了像素灰度值的動態(tài)范圍,從而達(dá)到增強(qiáng)圖像整體對比度的效果。!1!實(shí)驗(yàn)步驟(一)數(shù)字圖像的灰度變換1、線性灰度變換1)讀取一幅對比度低的灰度圖像并顯示。2)以m文件形式編寫matlab代碼,實(shí)現(xiàn)數(shù)字圖像的灰度范圍由[a,b倒[c,d]的線性拉伸,以便于提升原圖像的對比度。線性灰度變換公式如下:d/
4、(x,y)>bg(x,y)=<-Q]+ca(x,y)
5、(2)幕律(伽馬)變換為:3)顯示變換后的圖像,并利用matlabI具畫出對數(shù)或幕律變換曲線。4)以另一個文件名形式保存變換后的圖像。(二)直方圖處理1)讀取1幅灰度圖像并顯示原始圖像。2)編寫m文件實(shí)現(xiàn)對輸入圖像2種定義下的直方圖統(tǒng)計(jì)統(tǒng)計(jì)(即各個灰度級出現(xiàn)的次數(shù)及概率分布),并分別畫出兩種定義下的直方圖。3)對輸入灰度圖像進(jìn)行直方圖均衡化處理,分別顯示均衡處理前后的圖像和對應(yīng)的直方圖。五.實(shí)驗(yàn)結(jié)果及分析(一)數(shù)字圖像的灰度變換1、線性灰度變換實(shí)驗(yàn)結(jié)果:3000250020001500100050000501001502002502、非線性灰度變換(對
6、數(shù)變換或幕律變換選作其一)12000100008000600040002000050010000100200(二)直方圖處理原始圖像均衡化后的圖俟0100200結(jié)果分析:實(shí)驗(yàn)結(jié)果與預(yù)期一致六.程序源代碼數(shù)字圖像的灰度變換1、線性灰度變換程序源代碼:Il=imread(,Tl.jpg');ll=rgb2gray(II);subplot(221);imshow(Il);a二min(min仃1));b二max(max(II));B=(250-10)/double((b-a))*(11-a)+10;subplot(222);imshow(B);subplot(
7、223);imhist(II);subplot(224);imhist(B);2、非線性灰度變換(對數(shù)變換或幕律變換選作其一)subplot(231)I=imread(,T2.jpg');I=rgb2gray(1);imshow(I);a二min(min(I));b=max(max(T));c二2,r=0.8G=uint8(12*log(1+double(1)));subplot(232)imshow(G);subplot(233)z=0:1:1000;plot(z,12*log(double(z)));subplot(234)imhist(I);su
8、bplot(235)imhist(G);(二)直方圖處理I=imrcad(,pout?tif'