圖像處理實(shí)驗(yàn)3

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

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

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

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

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

6、要使用本文檔，請(qǐng)點(diǎn)擊下載按鈕下載!>> imshow(g)第二題：代碼： f=imread('pollen.tif');>> imhist(f);>> h=imhist(f,25);>> horz=linspace(0,255,25);>> bar(horz,h)6 / 15-6-如果您需要使用本文檔，請(qǐng)點(diǎn)擊下載按鈕下載!>> 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-如果您需要使用本文檔，請(qǐng)點(diǎn)擊下載按鈕下載!第三題：代碼：f=imread('pollen.tif');>> imshow(f);>> figure,imhist(f)>> ylim('auto')>> g=histeq(f,256);>> figure,imshow(g);8 / 15-8-如果您需要使用本文檔，請(qǐng)點(diǎn)擊下載按鈕下載!>> figure,imhist(g);

9、>> ylim('auto')>>截圖：9 / 15-9-如果您需要使用本文檔，請(qǐng)點(diǎn)擊下載按鈕下載!第四題：代碼： 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-如果您需要使用本文檔，請(qǐng)點(diǎn)擊下載按鈕下載!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-如果您需要使用本文檔，請(qǐng)點(diǎn)