圖像處理實驗3

1、如果您需要使用本文檔，請點擊下載按鈕下載!昆明理工大學(xué)信息工程與自動化學(xué)院學(xué)生實驗報告（ 20162017學(xué)年 第 一 學(xué)期 ）課程名稱：數(shù)字圖像基礎(chǔ) 開課實驗室： 2016年 月 日年級、專業(yè)、班學(xué)號姓名成績實驗項目名稱圖像增強（1）指導(dǎo)教師教師評語該同學(xué)是否了解實驗原理：a.了解b.基本了解c.不了解該同學(xué)的實驗?zāi)芰Γ篴.強 b.中等 c.差 該同學(xué)的實驗是否達(dá)到要求：a.達(dá)到b.基本達(dá)到c.未達(dá)到實驗報告是否規(guī)范：a.規(guī)范b.基本規(guī)范c.不規(guī)范實驗過程是否詳細(xì)記錄：a.詳細(xì)b.一般 c.沒有 教師簽名： 2016年 11月 16日一、實驗?zāi)康募皟?nèi)容目的：掌握和熟悉matlab編程環(huán)境及語

2、言；掌握直方圖統(tǒng)計的算法和用途。內(nèi)容：1. 調(diào)試教材p25頁例2.1輸出類似教材圖2.3的結(jié)果。2. 調(diào)試教材p33頁例2.4，編寫一個程序，分別使用imhist、bar、stem、plot四種方式顯示一幅灰度圖像的直方圖3. 調(diào)試教材p37頁例2.5。4. 直方圖均衡化的公式如下所示：根據(jù)上式及課堂所講直方圖均衡化原理及方法，自己寫一個matlab函數(shù)實現(xiàn)對灰度圖像的直方圖均衡化功能（類似于matlab提供的histeq函數(shù)）。（提示：實現(xiàn)中使用matlab函數(shù)cumsum（p38）可能會使程序簡單些）。1 / 15-1-如果您需要使用本文檔，請點擊下載按鈕下載!二、要求1. 描述直方圖的概

3、念并解釋直方圖均衡化原理。2. 程序結(jié)構(gòu)清晰，運行結(jié)果正確。3. 對于第1、2、3小題在實驗報告中給出所調(diào)試的程序，及其運行結(jié)果，對第4小題描述程序的設(shè)計、實現(xiàn)和結(jié)果，并對結(jié)果進行分析。2 / 15-2-如果您需要使用本文檔，請點擊下載按鈕下載!一、描述直方圖的概念并解釋直方圖均衡化原理。直方圖的概念：直方圖是用來量化曝光量，能夠使我們真實、直觀地看出照片曝光情況的一張二維坐標(biāo)系，其橫軸代表的是圖像中的亮度，從左向右，從全黑逐漸過渡到全白；縱軸代表的則是圖像中處于這個亮度范圍的像素的相對數(shù)量。當(dāng)直方圖中的黑色色塊偏向于左邊時，說明這張照片的整體色調(diào)偏暗，也可以理解為照片欠曝。而當(dāng)直方圖中的黑色

4、色塊集中在右邊時，說明這張照片整體色調(diào)偏亮，除非特殊構(gòu)圖需要，否則我們可以理解為照片過曝。直方圖均衡化原理：直方圖均衡化是圖像處理領(lǐng)域中利用圖像直方圖對對比度進行調(diào)整的方法。這種方法通常用來增加許多圖像的局部對比度，尤其是當(dāng)圖像的有用數(shù)據(jù)的對比度相當(dāng)接近的時候。通過這種方法，亮度可以更好地在直方圖上分布。這樣就可以用于增強局部的對比度而不影響整體的對比度，直方圖均衡化通過有效地擴展常用的亮度來實現(xiàn)這種功能。二、程序代碼以及運行結(jié)果第一題：代碼：f=imread('breast.tif');>> imshow(f)3 / 15-3-如果您需要使用本文檔，請點擊下載按鈕

5、下載!>> g1=imadjust(f,0 1,1 0);>> imshow(g1)>> g2=imadjust(f,0.5 0.75,0 1);>> imshow(g2)4 / 15-4-如果您需要使用本文檔，請點擊下載按鈕下載!>> g3=imadjust(f, , ,2);>> imshow(g3) >> g=imadjust(f,stretchlim(f), );>> imshow(g)>>g=imadjust(f,stretchlim(f),1 0);5 / 15-5-如果您需

6、要使用本文檔，請點擊下載按鈕下載!>> imshow(g)第二題：代碼： f=imread('pollen.tif');>> imhist(f);>> h=imhist(f,25);>> horz=linspace(0,255,25);>> bar(horz,h)6 / 15-6-如果您需要使用本文檔，請點擊下載按鈕下載!>> axis(0 255 0 60000)>> set(gca,'xtick',0:50:255)>> set(gca,'ytick

7、9;,0:20000:60000)>> stem(horz,h,'fill')>> axis(0 255 0 60000)>> set(gca,'xtick',0:50:255)>> set(gca,'ytick',0:20000:60000)>> hc=imhist(f);>> plot(hc)%use the default values.>> axis(0 255 0 15000)>> set(gca,'xtick',0:50:2

8、55)>> set(gca,'ytick',0:2000:15000)截圖:7 / 15-7-如果您需要使用本文檔，請點擊下載按鈕下載!第三題：代碼：f=imread('pollen.tif');>> imshow(f);>> figure,imhist(f)>> ylim('auto')>> g=histeq(f,256);>> figure,imshow(g);8 / 15-8-如果您需要使用本文檔，請點擊下載按鈕下載!>> figure,imhist(g);

9、>> ylim('auto')>>截圖：9 / 15-9-如果您需要使用本文檔，請點擊下載按鈕下載!第四題：代碼： function e = example( i )%untitled 此處顯示有關(guān)此函數(shù)的摘要% 此處顯示詳細(xì)說明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=

10、1:256 gpeq(i)=sum(gp(find(s2=i);endfigure(1);bar(0:255,gp,'b');figure(2);bar(0:255,gpeq,'b');10 / 15-10-如果您需要使用本文檔，請點擊下載按鈕下載!x=i;for i=0:255 x(find(i=i)=s2(i+1);endfigure(3);imshow(i);figure(4);imshow(x);end>> g=imread('pollen.tif');>> example(g) >>截圖：原圖像直方圖：11 / 15-11-如果您需要使用本文檔，請點