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

1、昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院學(xué)生實(shí)驗(yàn)報(bào)告（ 20162017學(xué)年 第 一 學(xué)期 ）課程名稱：數(shù)字圖像基礎(chǔ) 開(kāi)課實(shí)驗(yàn)室： 2016年 月 日年級(jí)、專業(yè)、班學(xué)號(hào)姓名成績(jī)實(shí)驗(yàn)項(xiàng)目名稱圖像增強(qiáng)（1）指導(dǎo)教師教師評(píng)語(yǔ)該同學(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)過(guò)程是否詳細(xì)記錄：A.詳細(xì)B.一般 C.沒(méi)有 教師簽名： 2016年 11月 16日一、實(shí)驗(yàn)?zāi)康募皟?nèi)容目的：掌握和熟悉Matlab編程環(huán)境及語(yǔ)言；掌握直方圖統(tǒng)計(jì)的算法和用途。內(nèi)容：1.

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

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

4、調(diào)整的方法。這種方法通常用來(lái)增加許多圖像的局部對(duì)比度，尤其是當(dāng)圖像的有用數(shù)據(jù)的對(duì)比度相當(dāng)接近的時(shí)候。通過(guò)這種方法，亮度可以更好地在直方圖上分布。這樣就可以用于增強(qiáng)局部的對(duì)比度而不影響整體的對(duì)比度，直方圖均衡化通過(guò)有效地?cái)U(kuò)展常用的亮度來(lái)實(shí)現(xiàn)這種功能。二、程序代碼以及運(yùn)行結(jié)果第一題：代碼：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 此處顯示有關(guān)此函數(shù)的摘要% 此處顯示詳細(xì)說(shuō)明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