前言
上一篇文章 ,我們講解了圖像處理中的閾值函數,這一篇文章我們來做膨脹和腐蝕函數。
?
膨脹與腐蝕
說概念可能很難解釋,我們來看圖,首先是原圖:
膨脹以后會變成這樣:
腐蝕以后則會變成這樣:
看起來可能有些莫名其妙,明明是膨脹,為什么字反而變細了,而明明是腐蝕,為什么字反而變粗了。
實際上,所謂膨脹應該指:
較亮色塊膨脹。
而所謂腐蝕應該指:
較亮色塊腐蝕。
上面圖里面,由于背景白色是較亮色塊,所以膨脹時就把黑色較暗色塊的字壓扁了……相反腐蝕時,字就吸水膨脹了……
用數學公式表示就是:
說白了就是在指定大小的內核里找最暗或者最亮的像素點,并用該點替換掉內核錨點上的像素。?
?
實現
首先我們來實現膨脹dilate函數。
var dilate = function (__src, __size, __borderType, __dst){ __src || error(arguments.callee, IS_UNDEFINED_OR_NULL /* {line} */ ); if (__src.type && __src.type == "CV_RGBA" ){ var width = __src.col, height = __src.row, size = __size || 3 , dst = __dst || new Mat(height, width, CV_RGBA), dstData = dst.data; var start = size >> 1 ; var withBorderMat = copyMakeBorder(__src, start, start, 0, 0 , __borderType), mData = withBorderMat.data, mWidth = withBorderMat.col; var newOffset, total, nowX, offsetY, offsetI, nowOffset, i, j; if (size & 1 === 0 ){ error(arguments.callee, UNSPPORT_SIZE /* {line} */ ); return __src; } for (i = height; i-- ;){ offsetI = i * width; for (j = width; j-- ;){ newOffset = 0 ; total = 0 ; for (y = size; y-- ;){ offsetY = (y + i) * mWidth * 4 ; for (x = size; x-- ;){ nowX = (x + j) * 4 ; nowOffset = offsetY + nowX; (mData[nowOffset] + mData[nowOffset + 1] + mData[nowOffset + 2] > total) && (total = mData[nowOffset] + mData[nowOffset + 1] + mData[nowOffset + 2]) && (newOffset = nowOffset); } } dstData[(j + offsetI) * 4] = mData[newOffset]; dstData[(j + offsetI) * 4 + 1] = mData[newOffset + 1 ]; dstData[(j + offsetI) * 4 + 2] = mData[newOffset + 2 ]; dstData[(j + offsetI) * 4 + 3] = mData[newOffset + 3 ]; } } } else { error(arguments.callee, UNSPPORT_DATA_TYPE /* {line} */ ); } return dst; };
(mData[nowOffset] + mData[nowOffset + 1] + mData[nowOffset + 2] > total) && (total = mData[nowOffset] + mData[nowOffset + 1] + mData[nowOffset + 2]) && (newOffset = nowOffset);
這行代碼里面,我們先用RGB的數值和與上一個最大值total對比,然后如果新的值比較大,就把新的值賦給total,并把新的偏移量newOffset賦值當前偏移量nowOffset。
然后整個內核大小的面積循環完了就可以找到一個最大的total以及對應的偏移量newOffset。就可以賦值了:
dstData[(j + offsetI) * 4] = mData[newOffset];
dstData[(j + offsetI) * 4 + 1] = mData[newOffset + 1];
dstData[(j + offsetI) * 4 + 2] = mData[newOffset + 2];
dstData[(j + offsetI) * 4 + 3] = mData[newOffset + 3];
那么腐蝕erode函數則相反,找一個最小的值就行了。
var erode = function (__src, __size, __borderType, __dst){ __src || error(arguments.callee, IS_UNDEFINED_OR_NULL /* {line} */ ); if (__src.type && __src.type == "CV_RGBA" ){ var width = __src.col, height = __src.row, size = __size || 3 , dst = __dst || new Mat(height, width, CV_RGBA), dstData = dst.data; var start = size >> 1 ; var withBorderMat = copyMakeBorder(__src, start, start, 0, 0 , __borderType), mData = withBorderMat.data, mWidth = withBorderMat.col; var newOffset, total, nowX, offsetY, offsetI, nowOffset, i, j; if (size & 1 === 0 ){ error(arguments.callee, UNSPPORT_SIZE /* {line} */ ); return __src; } for (i = height; i-- ;){ offsetI = i * width; for (j = width; j-- ;){ newOffset = 0 ; total = 765 ; for (y = size; y-- ;){ offsetY = (y + i) * mWidth * 4 ; for (x = size; x-- ;){ nowX = (x + j) * 4 ; nowOffset = offsetY + nowX; (mData[nowOffset] + mData[nowOffset + 1] + mData[nowOffset + 2] < total) && (total = mData[nowOffset] + mData[nowOffset + 1] + mData[nowOffset + 2]) && (newOffset = nowOffset); } } dstData[(j + offsetI) * 4] = mData[newOffset]; dstData[(j + offsetI) * 4 + 1] = mData[newOffset + 1 ]; dstData[(j + offsetI) * 4 + 2] = mData[newOffset + 2 ]; dstData[(j + offsetI) * 4 + 3] = mData[newOffset + 3 ]; } } } else { error(arguments.callee, UNSPPORT_DATA_TYPE /* {line} */ ); } return dst; };
?
效果
?
系列目錄
?
? 參考資料
更多文章、技術交流、商務合作、聯系博主
微信掃碼或搜索:z360901061

微信掃一掃加我為好友
QQ號聯系: 360901061
您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。
【本文對您有幫助就好】元
