1、昆明理工大學信息工程與自動化學院學生實驗報告（ 20162017學年 第 一 學期 ）課程名稱：數字圖像基礎 開課實驗室： 2016年 月 日年級、專業、班學號姓名成績實驗項目名稱圖像增強（1）指導教師教師評語該同學是否了解實驗原理：A.了解B.基本了解C.不了解該同學的實驗能力：A.強 B.中等 C.差 該同學的實驗是否達到要求：A.達到B.基本達到C.未達到實驗報告是否規范：A.規范B.基本規范C.不規范實驗過程是否詳細記錄：A.詳細B.一般 C.沒有 教師簽名： 2016年 11月 16日一、實驗目的及內容目的：掌握和熟悉Matlab編程環境及語言；掌握直方圖統計的算法和用途。內容：1.

2、 調試教材P25頁例2.1輸出類似教材圖2.3的結果。2. 調試教材P33頁例2.4，編寫一個程序，分別使用imhist、bar、stem、plot四種方式顯示一幅灰度圖像的直方圖3. 調試教材P37頁例2.5。4. 直方圖均衡化的公式如下所示：根據上式及課堂所講直方圖均衡化原理及方法，自己寫一個Matlab函數實現對灰度圖像的直方圖均衡化功能（類似于Matlab提供的histeq函數）。（提示：實現中使用Matlab函數cumsum（P38）可能會使程序簡單些）。二、要求1. 描述直方圖的概念并解釋直方圖均衡化原理。2. 程序結構清晰，運行結果正確。3. 對于第1、2、3小題在實驗報告中給出

3、所調試的程序，及其運行結果，對第4小題描述程序的設計、實現和結果，并對結果進行分析。一、描述直方圖的概念并解釋直方圖均衡化原理。直方圖的概念：直方圖是用來量化曝光量，能夠使我們真實、直觀地看出照片曝光情況的一張二維坐標系，其橫軸代表的是圖像中的亮度，從左向右，從全黑逐漸過渡到全白；縱軸代表的則是圖像中處于這個亮度范圍的像素的相對數量。當直方圖中的黑色色塊偏向于左邊時，說明這張照片的整體色調偏暗，也可以理解為照片欠曝。而當直方圖中的黑色色塊集中在右邊時，說明這張照片整體色調偏亮，除非特殊構圖需要，否則我們可以理解為照片過曝。直方圖均衡化原理：直方圖均衡化是圖像處理領域中利用圖像直方圖對對比度進行

4、調整的方法。這種方法通常用來增加許多圖像的局部對比度，尤其是當圖像的有用數據的對比度相當接近的時候。通過這種方法，亮度可以更好地在直方圖上分布。這樣就可以用于增強局部的對比度而不影響整體的對比度，直方圖均衡化通過有效地擴展常用的亮度來實現這種功能。二、程序代碼以及運行結果第一題：代碼：f=imread(breast.tif); imshow(f) g1=imadjust(f,0 1,1 0); imshow(g1) g2=imadjust(f,0.5 0.75,0 1); imshow(g2) g3=imadjust(f, , ,2); imshow(g3) g=imadjust(f,stre

5、tchlim(f), ); imshow(g)g=imadjust(f,stretchlim(f),1 0); imshow(g)第二題：代碼： f=imread(pollen.tif); imhist(f); h=imhist(f,25); horz=linspace(0,255,25); bar(horz,h) axis(0 255 0 60000) set(gca,xtick,0:50:255) set(gca,ytick,0:20000:60000) stem(horz,h,fill) axis(0 255 0 60000) set(gca,xtick,0:50:255) set(gc

6、a,ytick,0:20000:60000) hc=imhist(f); plot(hc)%Use the default values. axis(0 255 0 15000) set(gca,xtick,0:50:255) set(gca,ytick,0:2000:15000)截圖:第三題：代碼：f=imread(pollen.tif); imshow(f); figure,imhist(f) ylim(auto) g=histeq(f,256); figure,imshow(g); figure,imhist(g); ylim(auto)截圖：第四題：代碼： function E = E

7、xample( I )%UNTITLED 此處顯示有關此函數的摘要% 此處顯示詳細說明m,n=size(I);GP=zeros(1,256);for k=0:255 GP(k+1)=length(find(I=k)/(m*n);endS1=zeros(1,256);for i=1:256 for j=1:i S1(i)=GP(j)+S1(i); endendS2=round(S1*256)+0.5);for i=1:256 GPeq(i)=sum(GP(find(S2=i);endfigure(1);bar(0:255,GP,b);figure(2);bar(0:255,GPeq,b);X=I;for i=0:255 X(find(I=i)=S2(i+1);e