【OpenCV入門教程之十】 形態(tài)學圖像處理(一):膨脹與腐蝕
一向族、理論與概念講解——從現(xiàn)象到本質(zhì)
形態(tài)學(morphology)一詞通常表示生物學的一個分支,該分支主要研究動植物的形態(tài)和結(jié)構(gòu)狸剃。而我們圖像處理中指的形態(tài)學怯屉,往往表示的是數(shù)學形態(tài)學抛猫。
數(shù)學形態(tài)學(Mathematical morphology) 是一門建立在格論和拓撲學基礎之上的圖像分析學科贼陶,是數(shù)學形態(tài)學圖像處理的基本理論刃泡。
其基本的運算包括:二值腐蝕和膨脹、二值開閉運算碉怔、骨架抽取烘贴、極限腐蝕、擊中擊不中變換撮胧、形態(tài)學梯度桨踪、Top-hat變換、顆粒分析芹啥、流域變換锻离、灰值腐蝕和膨脹、灰值開閉運算墓怀、灰值形態(tài)學梯度等汽纠。
最基本的形態(tài)學操作有二種:膨脹與腐蝕(Dilation與Erosion)。膨脹與腐蝕能實現(xiàn)多種多樣的功能捺疼,主要如下:
- 消除噪聲
- 分割(isolate)出獨立的圖像元素疏虫,在圖像中連接(join)相鄰的元素。
- 尋找圖像中的明顯的極大值區(qū)域或極小值區(qū)域
- 求出圖像的梯度
在進行腐蝕和膨脹的講解之前啤呼,首先需要注意,腐蝕和膨脹是對白色部分(高亮部分)而言的呢袱,不是黑色部分官扣。膨脹就是圖像中的高亮部分進行膨脹,“領域擴張”羞福,效果圖擁有比原圖更大的高亮區(qū)域惕蹄。腐蝕就是原圖中的高亮部分被腐蝕,“領域被蠶食”,效果圖擁有比原圖更小的高亮區(qū)域卖陵。
1.1 結(jié)構(gòu)元和腐蝕
1.2 膨脹
按數(shù)學方面來說遭顶,膨脹或者腐蝕操作就是將圖像(或圖像的一部分區(qū)域,我們稱之為A)與核(我們稱之為B)進行卷積泪蔫。核可以是任何的形狀和大小棒旗,它擁有一個單獨定義出來的參考點,我們稱其為錨點(anchorpoint)撩荣。多數(shù)情況下铣揉,核是一個小的中間帶有參考點和實心正方形或者圓盤,其實餐曹,我們可以把核視為模板或者掩碼逛拱。
而膨脹就是求局部最大值的操作,核B與圖形卷積台猴,即計算核B覆蓋的區(qū)域(體現(xiàn)局部)的像素點的最大值朽合,并把這個最大值賦值給參考點指定的像素。這樣就會使圖像中的高亮區(qū)域逐漸增長饱狂。
注意:其實右圖要比左圖大了一圈
膨脹可以理解為B的中心(錨點)沿著A的外邊界走了一圈旁舰。膨脹是對高亮部分而言,A區(qū)域之外的部分 < A的高亮像素嗡官,所里外面被里面取代箭窜。
效果圖,高亮部分膨脹
膨脹的數(shù)學表達式:
用(x, y)
周邊區(qū)域(x+x', y+y')
內(nèi)的最大值代替(x, y)
的值衍腥。
1.3 腐蝕
腐蝕與膨脹是相反的操作磺樱,腐蝕是求局部最小值。
可與膨脹對比理解婆咸。
注意:其實右圖要比左圖小了一圈
腐蝕可以理解為B的中心(錨點)沿著A的內(nèi)邊界走了一圈竹捉。腐蝕也是對高亮部分而言,A區(qū)域之外的部分 < A的高亮像素尚骄,所里里面被外面取代块差。
A中能完全包含B的像素被留下來了。
效果圖倔丈,高亮部分被腐蝕
二憨闰、函數(shù)和實例
2.1 函數(shù)源碼
opencv\sources\modules\imgproc\src\morph.cpp
// 腐蝕
void cv::erode( InputArray src, OutputArray dst, InputArray kernel,
Point anchor, int iterations,
int borderType, const Scalar& borderValue )
{
CV_INSTRUMENT_REGION()
morphOp( MORPH_ERODE, src, dst, kernel, anchor, iterations, borderType, borderValue );
}
// 膨脹
void cv::dilate( InputArray src, OutputArray dst, InputArray kernel,
Point anchor, int iterations,
int borderType, const Scalar& borderValue )
{
CV_INSTRUMENT_REGION()
morphOp( MORPH_DILATE, src, dst, kernel, anchor, iterations, borderType, borderValue );
}
erode和dilate這兩個函數(shù)內(nèi)部就是調(diào)用了一下morphOp,只是調(diào)用morphOp時需五,第一個參數(shù)標識符不同鹉动,一個為MORPH_ERODE(腐蝕),一個為MORPH_DILATE(膨脹)宏邮。
這些函數(shù)在imgproc.hpp
中后面的參數(shù)是設置了默認值泽示。
void dilate( InputArray src, OutputArray dst, InputArray kernel,
Point anchor = Point(-1,-1), int iterations = 1,
int borderType = BORDER_CONSTANT,
const Scalar& borderValue = morphologyDefaultBorderValue() );
- 第三個參數(shù)缸血,InputArray類型的kernel,膨脹操作的核械筛。若為NULL時捎泻,表示的是使用參考點位于中心3x3的核。我們一般使用函數(shù)
getStructuringElement
返回指定形狀和尺寸的 結(jié)構(gòu)元(SE)埋哟。
Mat getStructuringElement(int shape, Size ksize, Point anchor = Point(-1,-1));
# 第一個參數(shù)表示內(nèi)核的形狀笆豁,我們可以選擇如下三種形狀之一
# 矩形: MORPH_RECT
# 交叉形: MORPH_CROSS
# 橢圓形: MORPH_ELLIPSE
#
# 第二和第三個參數(shù)分別是內(nèi)核的尺寸以及錨點的位置。
# 我們一般在調(diào)用erode以及dilate函數(shù)之前定欧,先定義一個Mat類型的變量來獲得getStructuringElement函數(shù)的返回值渔呵。
# 對于錨點的位置,有默認值Point(-1,-1)砍鸠,表示錨點位于中心扩氢。
# 且需要注意,交叉形的element形狀唯一依賴于錨點的位置爷辱。
# 而在其他情況下录豺,錨點只是影響了形態(tài)學運算結(jié)果的偏移。
- 第四個參數(shù)饭弓,Point類型的anchor双饥,錨的位置,默認值(-1弟断,-1)咏花,表示錨位于中心。
- 第五個參數(shù)阀趴,int類型的iterations昏翰,迭代使用erode()函數(shù)的次數(shù),默認值為1刘急。
- 第六個參數(shù)棚菊,int類型的borderType,用于推斷圖像外部像素的某種邊界模式叔汁。注意它有默認值BORDER_DEFAULT统求。
- 第七個參數(shù),const Scalar&類型的borderValue据块,當邊界為常數(shù)時的邊界值码邻,有默認值morphologyDefaultBorderValue(),一般我們不用去管他瑰钮。需要用到它時冒滩,可以看官方文檔中的createMorphologyFilter()函數(shù)得到更詳細的解釋。
2.2 實例
#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/highgui/highgui.hpp"
#include <iostream>
using namespace cv;
using namespace std;
int main() {
Mat src = imread("../pics/pig.jpg");
namedWindow("src");
imshow("src", src);
// 獲得結(jié)構(gòu)元
Mat se = getStructuringElement(MORPH_RECT, Size(10, 10));
Mat dst;
dilate(src, dst, se); // 后4個參數(shù)使用默認值
namedWindow("膨脹");
imshow("膨脹", dst);
waitKey(0);
}
Size(10, 10)
將膨脹代碼給為腐蝕
erode(src, dst, se);
三浪谴、總結(jié)
雖然膨脹和腐蝕是相反的操作开睡,但是如果用同樣的SE連續(xù)執(zhí)行2個操作,也不一定能恢復原圖苟耻。其實就是開閉操作篇恒。
先腐蝕再膨脹 其實就是開操作??
開操作一般會平滑物體的輪廓,斷開較窄的狹頸并消除細的突出物凶杖。
erode(src, dst, se);
namedWindow("先腐蝕");
imshow("先腐蝕", dst);
dilate(dst, dst, se);
namedWindow("再膨脹");
imshow("再膨脹", dst);
先膨脹再腐蝕 其實就是閉操作??
閉操作也會平滑輪廓的一部分胁艰,但與開操作相反,通常會彌合較窄的間斷和細長的溝壑智蝠,消除小的孔洞腾么,填補輪廓線中的斷裂。
dilate(src, dst, se);
namedWindow("先膨脹");
imshow("先膨脹", dst);
erode(dst, dst, se);
namedWindow("再腐蝕");
imshow("再腐蝕", dst);