opencv-图像形态学之开运算、闭运算、形态学梯度、顶帽、黑帽合辑

转自：https://blog.csdn.net/poem_qianmo/article/details/24599073

1.1 开运算（Opening Operation）

开运算（Opening Operation），其实就是先腐蚀后膨胀的过程。其数学表达式如下：

开运算可以用来消除小物体、在纤细点处分离物体、平滑较大物体的边界的同时并不明显改变其面积。效果图是这样的：

1.2 闭运算(Closing Operation)

先膨胀后腐蚀的过程称为闭运算(Closing Operation)，其数学表达式如下：

 

闭运算能够排除小型黑洞(黑色区域)。效果图如下所示：

 

1.3 形态学梯度（MorphologicalGradient）

形态学梯度（Morphological Gradient）为膨胀图与腐蚀图之差，数学表达式如下：

对二值图像进行这一操作可以将团块（blob）的边缘突出出来。我们可以用形态学梯度来保留物体的边缘轮廓，如下所示:

 

1.4 顶帽（Top Hat）

顶帽运算（Top Hat）又常常被译为”礼帽“运算。为原图像与上文刚刚介绍的“开运算“的结果图之差，数学表达式如下：

因为开运算带来的结果是放大了裂缝或者局部低亮度的区域，因此，从原图中减去开运算后的图，得到的效果图突出了比原图轮廓周围的区域更明亮的区域，且这一操作和选择的核的大小相关。

顶帽运算往往用来分离比邻近点亮一些的斑块。当一幅图像具有大幅的背景的时候，而微小物品比较有规律的情况下，可以使用顶帽运算进行背景提取。

如下所示:

1.5 黑帽（Black Hat）

黑帽（Black Hat）运算为”闭运算“的结果图与原图像之差。数学表达式为：

黑帽运算后的效果图突出了比原图轮廓周围的区域更暗的区域，且这一操作和选择的核的大小相关。

所以，黑帽运算用来分离比邻近点暗一些的斑块。非常完美的轮廓效果图：

 

 

二、深入——OpenCV源码分析溯源

本文的主角是OpenCV中的morphologyEx函数，它利用基本的膨胀和腐蚀技术，来执行更加高级的形态学变换，如开闭运算，形态学梯度，“顶帽”、“黑帽”等等。这一节我们来一起看一下morphologyEx函数的源代码。

//-----------------------------------【erode（）函数中文注释版源代码】----------------------------
 
// 说明：以下代码为来自于计算机开源视觉库OpenCV的官方源代码
 
// OpenCV源代码版本：2.4.8
 
// 源码路径：…opencvsourcesmodulesimgprocsrcmorph.cpp
 
// 源文件中如下代码的起始行数：1369行
 
// 中文注释by浅墨
 
//--------------------------------------------------------------------------------------------------------
 
void cv::morphologyEx( InputArray _src,OutputArray _dst, int op,
 
InputArray kernel, Pointanchor, int iterations,
 
int borderType, constScalar& borderValue )
 
{
 
//拷贝Mat数据到临时变量
 
Mat src = _src.getMat(), temp;
 
_dst.create(src.size(), src.type());
 
Mat dst = _dst.getMat();
 
 
 
//一个大switch，根据不同的标识符取不同的操作
 
switch( op )
 
{
 
case MORPH_ERODE:
 
erode( src, dst, kernel, anchor, iterations, borderType, borderValue );
 
break;
 
case MORPH_DILATE:
 
dilate( src, dst, kernel, anchor, iterations, borderType, borderValue );
 
break;
 
case MORPH_OPEN:
 
erode( src, dst, kernel, anchor, iterations, borderType, borderValue );
 
dilate( dst, dst, kernel, anchor, iterations, borderType, borderValue );
 
break;
 
case CV_MOP_CLOSE:
 
dilate( src, dst, kernel, anchor, iterations, borderType, borderValue );
 
erode( dst, dst, kernel, anchor, iterations, borderType, borderValue );
 
break;
 
case CV_MOP_GRADIENT:
 
erode( src, temp, kernel, anchor, iterations, borderType, borderValue );
 
dilate( src, dst, kernel, anchor, iterations, borderType, borderValue );
 
dst -= temp;
 
break;
 
case CV_MOP_TOPHAT:
 
if( src.data != dst.data )
 
temp = dst;
 
erode( src, temp, kernel, anchor, iterations, borderType, borderValue );
 
dilate( temp, temp, kernel, anchor,iterations, borderType, borderValue );
 
dst = src - temp;
 
break;
 
case CV_MOP_BLACKHAT:
 
if( src.data != dst.data )
 
temp = dst;
 
dilate( src, temp, kernel, anchor, iterations, borderType, borderValue);
 
erode( temp, temp, kernel, anchor, iterations, borderType, borderValue);
 
dst = temp - src;
 
break;
 
default:
 
CV_Error( CV_StsBadArg, "unknown morphological operation" );
 
}
 
}

看上面的源码可以发现，其实morphologyEx函数其实就是内部一个大switch而已。根据不同的标识符取不同的操作。比如开运算MORPH_OPEN，按我们上文中讲解的数学表达式，就是先腐蚀后膨胀，即依次调用erode和dilate函数，为非常简明干净的代码。

三、浅出——API函数快速上手

3.1 morphologyEx函数详解

上面我们已经讲到，morphologyEx函数利用基本的膨胀和腐蚀技术，来执行更加高级形态学变换，如开闭运算，形态学梯度，“顶帽”、“黑帽”等等。这一节我们来了解它的参数意义和使用方法。

C++: void morphologyEx(
 
InputArray src,
 
OutputArray dst,
 
int op,
 
InputArraykernel,
 
Pointanchor=Point(-1,-1),
 
intiterations=1,
 
intborderType=BORDER_CONSTANT,
 
constScalar& borderValue=morphologyDefaultBorderValue() );

• 第一个参数，InputArray类型的src，输入图像，即源图像，填Mat类的对象即可。图像位深应该为以下五种之一：CV_8U, CV_16U,CV_16S, CV_32F 或CV_64F。
• 第二个参数，OutputArray类型的dst，即目标图像，函数的输出参数，需要和源图片有一样的尺寸和类型。
• 第三个参数，int类型的op，表示形态学运算的类型，可以是如下之一的标识符：

• • MORPH_OPEN – 开运算（Opening operation）
  • MORPH_CLOSE – 闭运算（Closing operation）
  • MORPH_GRADIENT -形态学梯度（Morphological gradient）
  • MORPH_TOPHAT - “顶帽”（“Top hat”）
  • MORPH_BLACKHAT - “黑帽”（“Black hat“）

 

另有CV版本的标识符也可选择，如CV_MOP_CLOSE，CV_MOP_GRADIENT，CV_MOP_TOPHAT，CV_MOP_BLACKHAT，这应该是OpenCV1.0系列版本遗留下来的标识符，和上面的“MORPH_OPEN”一样的效果。

• 第四个参数，InputArray类型的kernel，形态学运算的内核。若为NULL时，表示的是使用参考点位于中心3x3的核。我们一般使用函数 getStructuringElement配合这个参数的使用。getStructuringElement函数会返回指定形状和尺寸的结构元素（内核矩阵）。关于getStructuringElement我们上篇文章中讲过了，这里为了大家参阅方便，再写一遍：

其中，getStructuringElement函数的第一个参数表示内核的形状，我们可以选择如下三种形状之一:

• • 矩形: MORPH_RECT
  • 交叉形: MORPH_CROSS
  • 椭圆形: MORPH_ELLIPSE

而getStructuringElement函数的第二和第三个参数分别是内核的尺寸以及锚点的位置。

我们一般在调用erode以及dilate函数之前，先定义一个Mat类型的变量来获得getStructuringElement函数的返回值。对于锚点的位置，有默认值Point(-1,-1)，表示锚点位于中心。且需要注意，十字形的element形状唯一依赖于锚点的位置。而在其他情况下，锚点只是影响了形态学运算结果的偏移。

getStructuringElement函数相关的调用示例代码如下：

int g_nStructElementSize = 3; //结构元素(内核矩阵)的尺寸
 
//获取自定义核
 
Mat element =getStructuringElement(MORPH_RECT,
 
Size(2*g_nStructElementSize+1,2*g_nStructElementSize+1),
 
Point(g_nStructElementSize, g_nStructElementSize ));

调用这样之后，我们便可以在接下来调用erode、dilate或morphologyEx函数时，kernel参数填保存getStructuringElement返回值的Mat类型变量。对应于我们上面的示例，就是填element变量。

• 第五个参数，Point类型的anchor，锚的位置，其有默认值（-1，-1），表示锚位于中心。
• 第六个参数，int类型的iterations，迭代使用函数的次数，默认值为1。
• 第七个参数，int类型的borderType，用于推断图像外部像素的某种边界模式。注意它有默认值BORDER_ CONSTANT。
• 第八个参数，const Scalar&类型的borderValue，当边界为常数时的边界值，有默认值morphologyDefaultBorderValue()，一般我们不用去管他。需要用到它时，可以看官方文档中的createMorphologyFilter()函数得到更详细的解释。

 

其中的这些操作都可以进行就地（in-place）操作。且对于多通道图像，每一个通道都是单独进行操作。