數(shù)字圖像處理實(shí)驗(yàn)報(bào)告63445298

1、數(shù)字圖像處理實(shí)驗(yàn)報(bào)告一、實(shí)驗(yàn)內(nèi)容及目的：已知一幅圖上有五個(gè)瓶子，裝有液體?；谶@幅圖片，設(shè)計(jì)一個(gè)檢測未完全裝滿的瓶子的解決方案。二、實(shí)驗(yàn)思路分析：形態(tài)學(xué)操作，使圖像便于處理二值化原圖找出不符合要求的標(biāo)記塊顯示原圖標(biāo)記、計(jì)算三、實(shí)驗(yàn)步驟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(開運(yùn)算之后圖像); % 輸出開運(yùn)算之后圖像5、標(biāo)記連通域，計(jì)算出各個(gè)白色連通區(qū)域（未裝灌區(qū)域）面積和質(zhì)心坐標(biāo)：圖像處理后五個(gè)白色區(qū)域面積為：2374;2739;8381;2739;1660五個(gè)質(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)計(jì)算結(jié)果判斷那個(gè)裝灌不合格：通過沒有裝灌部分的面積大于2900或者質(zhì)心大于50則判斷為不合格，由步驟5的連通面積和質(zhì)心坐標(biāo)可得第三個(gè)瓶子不合格。四、實(shí)驗(yàn)代碼clear all;close all;t = 200; % 全局閾值 200時(shí)效果要好一些，手動選出的值%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(開運(yùn)算之后圖像); % 輸出開運(yùn)算之后圖像% step4 計(jì)算出各個(gè)白色連通區(qū)域（未裝灌區(qū)域）面積和質(zhì)心坐標(biāo)l,num = bwlabel(fo,4); % 為白色連通區(qū)域做標(biāo)簽stats = regionprops(l,basic); % 利用 regionprops（）函數(shù)計(jì)算面積、質(zhì)心和邊界等s = stats.area; % 矩陣s 中保存所有白色連通區(qū)域的面積zx = stats.centroid; % zx保存白色連通區(qū)域質(zhì)心的坐標(biāo)% step5 根據(jù)面積和質(zhì)心坐標(biāo)計(jì)算結(jié)果判斷那個(gè)裝灌不合格for indexnum = 1:num % 有實(shí)際幾個(gè)連通區(qū)域就做幾個(gè)循環(huán),實(shí)際就是小于等于5個(gè)if ( (s(indexnum) 2900) | (zx(