數(shù)字圖像處理實驗報告63445298

1、數(shù)字圖像處理實驗報告一、實驗內(nèi)容及目的：已知一幅圖上有五個瓶子，裝有液體?；谶@幅圖片，設(shè)計一個檢測未完全裝滿的瓶子的解決方案。二、實驗思路分析：形態(tài)學(xué)操作，使圖像便于處理二值化原圖找出不符合要求的標(biāo)記塊顯示原圖標(biāo)記、計算三、實驗步驟1、打開matlab軟件，并將current folder改為所要處理的圖片fig1.jpg所在目錄。我的圖片所在目錄為e:12-13-2數(shù)字圖象處理di：2、讀取和顯示原始圖像，顯示原始圖像的直方圖：picop = imread(fig1.jpg); % 讀入圖像figure,imshow(picop),title(原始灰度圖像); % 顯示圖像figure,

2、imhist(picop), title(原始灰度圖像直方圖); % 顯示原始圖像的直方圖3、將圖像進(jìn)行二值化處理：由于原始圖像中白色比較明顯，瓶身的顏色和背景比較相近，所以直接進(jìn)行二值化處理，將未裝滿液體的信息提取出來。picb = picop; % 復(fù)制灰度圖像到picbpicbw = im2bw(picb,t/255); % 采用全局閾值進(jìn)行灰度圖像轉(zhuǎn)變?yōu)槎祱D像figure,imshow(picbw),title(全局閾值下二值圖像); % 顯示二值圖像4、進(jìn)行形態(tài)學(xué)操作，轉(zhuǎn)化為有利于處理的圖像：picmorphology = picbw ; % 復(fù)制準(zhǔn)備形態(tài)學(xué)處理se = strel

3、(square,10); % 結(jié)構(gòu)化元素fo = imopen(picmorphology,se); % 開操作figure,imshow(fo),title(開運算之后圖像); % 輸出開運算之后圖像5、標(biāo)記連通域，計算出各個白色連通區(qū)域（未裝灌區(qū)域）面積和質(zhì)心坐標(biāo)：圖像處理后五個白色區(qū)域面積為：2374;2739;8381;2739;1660五個質(zhì)心的坐標(biāo)數(shù)據(jù)：20.6251053074979;48.0433866891323;138.604600219058;46.6575392479007;256.266674621167;84.9748240066818;376.6046002190

4、58;46.6575392479007;484.618072289157;47.90843373493986、根據(jù)面積和質(zhì)心坐標(biāo)計算結(jié)果判斷那個裝灌不合格：通過沒有裝灌部分的面積大于2900或者質(zhì)心大于50則判斷為不合格，由步驟5的連通面積和質(zhì)心坐標(biāo)可得第三個瓶子不合格。四、實驗代碼clear all;close all;t = 200; % 全局閾值 200時效果要好一些，手動選出的值%step1 讀取和顯示原始圖像，顯示原始圖像的直方圖picop = imread(fig1.jpg); % 讀入圖像figure,imshow(picop),title(原始灰度圖像); % 顯示圖像figu

5、re, imhist(picop), title(原始灰度圖像直方圖); % 顯示原始圖像的直方圖% step2 轉(zhuǎn)化為二值圖像picb = picop; % 復(fù)制灰度圖像到picbpicbw = im2bw(picb,t/255); % 采用全局閾值進(jìn)行灰度圖像轉(zhuǎn)變?yōu)槎祱D像figure,imshow(picbw),title(全局閾值下二值圖像); % 顯示二值圖像% step3 進(jìn)行形態(tài)學(xué)操作，轉(zhuǎn)化為有利于處理的圖像picmorphology = picbw ; % 復(fù)制準(zhǔn)備形態(tài)學(xué)處理se = strel(square,10); % 結(jié)構(gòu)化元素fo = imopen(picmorphol

6、ogy,se); % 開操作figure,imshow(fo),title(開運算之后圖像); % 輸出開運算之后圖像% step4 計算出各個白色連通區(qū)域（未裝灌區(qū)域）面積和質(zhì)心坐標(biāo)l,num = bwlabel(fo,4); % 為白色連通區(qū)域做標(biāo)簽stats = regionprops(l,basic); % 利用 regionprops（）函數(shù)計算面積、質(zhì)心和邊界等s = stats.area; % 矩陣s 中保存所有白色連通區(qū)域的面積zx = stats.centroid; % zx保存白色連通區(qū)域質(zhì)心的坐標(biāo)% step5 根據(jù)面積和質(zhì)心坐標(biāo)計算結(jié)果判斷那個裝灌不合格for indexnum = 1:num % 有實際幾個連通區(qū)域就做幾個循環(huán),實際就是小于等于5個if ( (s(indexnum) 2900) | (zx(