OpenCvSharp学习笔记11--一文搞懂击中击不中(HitMiss)

目的

学会使用形态学操作(MorphologyEx)中的击中击不中(MorphTypes.HitMiss)算法。

原理

OpenCV官网详细讲解击中击不中的算法，这里简单说明如下。

假设源图像为A，结构元素(核)B1和结构元素B2

击中击不中算法的运算如下：

1、使用结构元素B1对源图像A进行腐蚀操作；

2、使用结构元素B2对源图像A的互补图(取反)进行腐蚀操作(Erode)；

3、将步骤1与步骤2的结果进行和操作(AND)，即为输出结果。

可以将B1和B2组合为结构元素B，如下

如果上面的话理解起来有些晦涩的话，我个人的理解是：定义一个待匹配的核，其中核内元素值为1，表示该位置需要匹配前景（白色）；若为-1，表示该位置需要匹配背景（黑色）；若为0，表示任意(前景、背景皆可)。使用该核对源图像扫描后，若匹配上述规则，在锚点位置记为255(非0)，若不匹配，则锚点位置记为0，最后得到的结果就是输出图像。

注意，在外推边框没有指定值时，左上角和右下角也出现匹配情况。

将上图放大50倍后显示效果如下

另一个核

 再换一个核

 

通上面三个示例，应该不难明白击中击不中算法是如何匹配的了。 

本人将该算法应用到实际生产中的十字架定位，效果很好。

最后，在使用时HitMiss算法，建议还是指定边框外推像素的默认值，避免多余匹配点。

源码示例

public void Run(ParamBase paramBase)
{
    //源图像
    using Mat inputImage = new Mat(8, 8, MatType.CV_8UC1, new byte[] {0,   0,   0,   0,  0,   0, 255,   0,
                                                                      0, 255, 255, 255,  0,   0,   0, 255,
                                                                      0, 255, 255, 255,  0,   0,   0,   0,
                                                                      0, 255, 255, 255,  0, 255,   0,   0,
                                                                      0,   0, 255,   0,  0,   0,   0,   0,
                                                                      0,   0, 255,   0,  0, 255, 255,   0,
                                                                      0, 255,   0, 255,  0,   0, 255,   0,
                                                                      0, 255, 255,   0,  0, 255,   0, 255 });

    //待匹配的核
    //中间为黑(-1),上、下、左、右为白(1)，四个角任意(0)
    Mat kernel = new Mat(3, 3, MatType.CV_32SC1, new int[] {0, 1, 0,
                                                            1,-1, 1,
                                                            0, 1, 0});

    TestHitMiss(inputImage.Clone(), kernel);

    //核2
    kernel = new Mat(3, 3, MatType.CV_32SC1, new int[] {0,-1,-1,
                                                        1, 1,-1,
                                                        0, 1, 0});
    TestHitMiss(inputImage.Clone(), kernel);

    //核3
    kernel = new Mat(3, 3, MatType.CV_32SC1, new int[] {-1,-1,0,
                                                        -1, 1,0,
                                                        -1,-1,0});
    TestHitMiss(inputImage.Clone(), kernel);
}

/// <summary>
/// 测试HitMiss算法
/// </summary>
/// <param name="inputImage">源图像</param>
/// <param name="kernel">待匹配的核</param>
private void TestHitMiss(Mat inputImage,Mat kernel)
{
    using Mat outputImage1 = new Mat();
    //使用形态学的击中击不中算法，查找上面的核。如果匹配，在锚点位置记为255，否则记为0
    Cv2.MorphologyEx(inputImage, outputImage1, MorphTypes.HitMiss, kernel);

    using Mat outputImage2 = new Mat();
    //注意，默认的borderValue为Scalar.All(double.MaxValue)，这里改为0
    Cv2.MorphologyEx(inputImage, outputImage2, MorphTypes.HitMiss, kernel, borderType: BorderTypes.Constant, borderValue: Scalar.All(0));

    //放大倍数
    const int rate = 50;
    kernel = (kernel + 1) * 127;//将-1显示为黑，0显示为灰，1显示为白
    kernel.ConvertTo(kernel, MatType.CV_8U);

    ResizeAndShow(kernel, rate, "kernel");
    ResizeAndShow(inputImage, rate, "Original");
    ResizeAndShow(outputImage1, rate, "Hit or Miss without Set BorderValue");
    //建议在做HitMiss时，给外框指定值
    ResizeAndShow(outputImage2, rate, "Hit or Miss Set BorderValue=0");
    Utils.WaitDestoryAll();
}

/// <summary>
/// 放大显示
/// </summary>
/// <param name="mat"></param>
/// <param name="rate"></param>
/// <param name="winName"></param>
private void ResizeAndShow(Mat mat, int rate, string winName)
{
    Cv2.Resize(mat, mat, new Size(), rate, rate, InterpolationFlags.Nearest);
    Cv2.ImShow(winName, mat);
}

OpenCvSharp函数示例(目录)

参考

OpenCV: Hit-or-Miss