利用C#版OpenCV实现圆心求取实例代码

前言

opencvsharp是opencv的.net wrapper，是一名日本工程师开发的，项目地址为：。

该源码是 bsd开放协议，bsd开源协议是一个给于使用者很大*的协议。基本上使用者可以”为所欲为”,可以*的使用，修改源代码，也可以将修改后的代码作为开源或者专有软件再发布或商业化销售。

1.opencvsharp的下载

可以直接从上面的github上下载源码，自行编译引用；

也可用vs中的nuget包管理器下载；

打开【工具】->【库程序包管理器】->【管理解决方案的nuget程序包】，在其中搜索opencvsharp，选择合适的点击【安装】（最好安装最新的）。一直等待完成。

我的环境是vs2017，下载之后最好直接拷贝opencvsharp系列dll到项目中引用即可。

2.拟合圆并求取圆心

本次举例比较贴近实际，我们求如下原始图片的中间部分圆的圆心，选取的图片故意只留了一半：

源码如下：

using opencvsharp;

namespace sampleswinform
{
  public class circlefit 
  {
    public void run()
    {
      //读取图片
      var img = cv2.imread("data/image/c1.bmp");
      //显示图片
      //cv2.imshow("input image", img);
      
      //转换成灰度图
      mat gray = img.cvtcolor(colorconversioncodes.bgr2gray);
      
      //阈值二值化操作 阈值参数可以用一些可视化工具来调试得到
      mat thresholdimg = gray.threshold(11, 255, thresholdtypes.binary);
      cv2.imshow("threshold", thresholdimg);

      //降噪
      //方法一：高斯变化
      //mat gaussimg= thresholdimg.gaussianblur(new size(5, 5), 0.8);
      //cv2.imshow("gaussianblur", gaussimg);
      //方法二：中值滤波降噪
      mat medianimg = thresholdimg.medianblur(5);
      cv2.imshow("medianblur", medianimg);

      //方法三：膨胀+腐蚀
      ////膨胀处理
      //mat kernel = new mat(15, 15, mattype.cv_8uc1);
      //mat dilateimg = thresholdimg.dilate(kernel);
      ////腐蚀处理
      //mat binary = dilateimg.erode(kernel);
      ////显示中间结果
      //cv2.imshow("dilate & erode", binary);

      //设置感兴趣的区域
      int x = 150, y = 100, w = 294, h = 337;
      rect roi = new rect(x, y, w, h);
      mat roiimg = new mat(medianimg, roi);
      //cv2.imshow("roi image", roiimg);

      //寻找图像轮廓
      point[][] contours;
      hierarchyindex[] hierachy;
      cv2.findcontours(roiimg, out contours, out hierachy, retrievalmodes.list, contourapproximationmodes.approxtc89kcos);
      //根据找到的轮廓点，拟合椭圆
      for (int i = 0; i < contours.length; i++)
      {
        //拟合函数必须至少5个点，少于则不拟合
        if (contours[i].length < 5) continue;
        //椭圆拟合
        var rrt = cv2.fitellipse(contours[i]);

        //roi复原
        rrt.center.x += x;
        rrt.center.y += y;

        //画椭圆
        cv2.ellipse(img, rrt, new scalar(0, 0, 255), 2, linetypes.antialias);
        //画圆心
        cv2.circle(img, (int)(rrt.center.x), (int)(rrt.center.y), 4, new scalar(255, 0, 0), -1, linetypes.link8, 0);
      }
      
      cv2.imshow("fit circle", img);
    }


   }
}

中间处理过程效果图如下：

分别是阈值分割后，中值滤波后，拟合圆之后

总结

以上就是这篇文章的全部内容了，希望本文的内容对大家的学习或者工作具有一定的参考学习价值，谢谢大家对的支持。