EmguCV学习笔记 VB.Net 6.S 特别示例

版权声明：本文为博主原创文章，转载请在显著位置标明本文出处以及作者网名，未经作者允许不得用于商业目的。

EmguCV是一个基于OpenCV的开源免费的跨平台计算机视觉库,它向C#和VB.NET开发者提供了OpenCV库的大部分功能。

教程VB.net版本请访问：EmguCV学习笔记 VB.Net 目录-CSDN博客

教程C#版本请访问：EmguCV学习笔记 C# 目录-CSDN博客

笔者的博客网址：https://blog.csdn.net/uruseibest

教程配套文件及相关说明以及如何获得pdf教程和代码，请移步：EmguCV学习笔记

学习VB.Net知识，请移步： vb.net 教程 目录_vb中如何用datagridview-CSDN博客

 学习C#知识，请移步：C# 教程 目录_c#教程目录-CSDN博客

 

6.S 特别示例

6.S.1图像中的圆半径

本示例来自一名网友咨询如何获得图像中圆形的半径，其中有两个十字作为标记，十字之间距离为100mm。如下图：

图6-28 根据给出的定位标记求圆半径

【代码位置：frmChapter6_S1】Button1_Click、PointFToPoint

    Private Sub Button1_Click(sender As Object, e As EventArgs) Handles Button1.Click

        Dim msrc As New Mat("C:\learnEmgucv\celiang.jpg", ImreadModes.Color)

        Dim mgray As New Mat()

        CvInvoke.CvtColor(msrc, mgray, ColorConversion.Bgr2Gray)

        Dim kernel As New Mat

        kernel = CvInvoke.GetStructuringElement(ElementShape.Cross, New Drawing.Size(3, 3), New Point(-1, -1))

        Dim merode As New Mat        '

        '这里使用了2次迭代

        CvInvoke.Dilate(mgray, merode, kernel, New Point(-1, -1), 1, BorderType.Constant, Nothing)

        CvInvoke.Threshold(merode, merode, 200, 255, ThresholdType.BinaryInv)

        ImageBox1.Image = merode

        '获得所有轮廓

        Dim contours As New VectorOfVectorOfPoint

        Dim hierarchy As New VectorOfRect

        CvInvoke.FindContours(merode, contours, hierarchy, RetrType.List, ChainApproxMethod.ChainApproxSimple)

        Dim m2 As New Mat(merode.Size, DepthType.Cv8U, 1)

        m2.SetTo(New MCvScalar(0))

        '圆轮廓

        Dim contourCircle As VectorOfPoint

        '圆轮廓的周长

        Dim perimeter As Double

        '绘制轮廓

        For i As Integer = 0 To contours.Size - 1

            Dim carea As VectorOfPoint = contours(i)

            '获得轮廓面积

            Dim area As Double = CvInvoke.ContourArea(carea, False)

            '符合条件时，绘制轮廓，排除圆形，只保留十字线

            '本图中圆形面积为2449，直线面积为8，需要根据实际情况调整

            If area < 200 Then

                CvInvoke.DrawContours(m2, contours, i, New MCvScalar(255), 0.4)

            Else

                '得到圆形，图像中只有三个轮廓：2个交叉十字线段、1个圆形

                '这里简化操作，否则在多个轮廓情况下，应获取最大面积的轮廓判断为圆形

                contourCircle = contours(i)

                '获取轮廓周长

                perimeter = CvInvoke.ArcLength(contourCircle, True)

            End If

        Next

        ImageBox1.Image = m2

        '使用HoughLinesP方法检测图像中的直线，并将其绘制到图像

        '因为本图中十字线上的线段较短，所以这里阈值设置很小

        Dim lines As LineSegment2D() = CvInvoke.HoughLinesP(m2, 1, Math.PI / 180, 5, 5, 80)

        Dim m3 As New Mat(merode.Size, DepthType.Cv8U, 3)

        m3.SetTo(New MCvScalar(0, 0, 0))

        For Each line As LineSegment2D In lines

            CvInvoke.Line(m3, line.P1, line.P2, New MCvScalar(0, 255, 0), 2)

        Next

        ImageBox2.Image = m3

        '对直线进行分类，将其分为垂直和水平两类：

        Dim verticalLines As New List(Of LineSegment2D)

        Dim horizontalLines As New List(Of LineSegment2D)

        '计算每条直线的倾斜角度来进行分类，

        '将倾斜角度在60 - 120度之间的直线划分为垂直类，

        '将倾斜角度在30 - 150度之间的直线划分为水平类。

        For Each line As LineSegment2D In lines

            Dim angle As Double = Math.Atan2(line.P2.Y - line.P1.Y, line.P2.X - line.P1.X) * 180 / Math.PI

            If angle < 0 Then angle += 180

            If angle > 60 AndAlso angle < 120 Then

                verticalLines.Add(line)

            ElseIf angle > 150 OrElse angle < 30 Then

                horizontalLines.Add(line)

            End If

        Next

        '对垂直和水平直线进行匹配，并计算十字中心点的位置：

        Dim intersections As New List(Of PointF)

        '得到两个相交点

        For Each verticalLine As LineSegment2D In verticalLines

            For Each horizontalLine As LineSegment2D In horizontalLines

                '基于图像中两条直线真实相交，

                '如果垂直线的中点X坐标在水平线两端点X坐标之间

                '那么，这条垂直线段和这条水平线段相交

                Dim centerX As Single = (verticalLine.P1.X + verticalLine.P2.X) / 2

                If horizontalLine.P1.X < horizontalLine.P2.X Then

                    If centerX > horizontalLine.P1.X And centerX < horizontalLine.P2.X Then

                        Dim intersectionPoint As New PointF(

                            (horizontalLine.P1.X + horizontalLine.P2.X + verticalLine.P1.X + verticalLine.P2.X) / 4,

                            (horizontalLine.P1.Y + horizontalLine.P2.Y + verticalLine.P1.Y + verticalLine.P2.Y) / 4

                        )

                        intersections.Add(intersectionPoint)

                    End If

                Else

                    If centerX > horizontalLine.P2.X And centerX < horizontalLine.P1.X Then

                        Dim intersectionPoint As New PointF(

                            (horizontalLine.P1.X + horizontalLine.P2.X + verticalLine.P1.X + verticalLine.P2.X) / 4,

                            (horizontalLine.P1.Y + horizontalLine.P2.Y + verticalLine.P1.Y + verticalLine.P2.Y) / 4

                        )

                        intersections.Add(intersectionPoint)

                    End If

                End If

            Next

        Next

        If intersections.Count <> 2 Then

            MessageBox.Show("未能获得两个十字线的交叉点")

            Exit Sub

        End If

        CvInvoke.Line(msrc, PointFToPoint(intersections(0)), PointFToPoint(intersections(1)), New MCvScalar(0, 255, 0), 2)

        CvInvoke.Imshow("m3", msrc)

        '计算两个交点的距离

        Dim distance As Double = Math.Sqrt(

            (intersections(0).X - intersections(1).X) ^ 2 +

            (intersections(0).Y - intersections(1).Y) ^ 2

        )

        '实际中两交点距离为100毫米，计算相应比例

        Dim proportion As Double = 100 / distance

        '以下是基于最小外接圆来计算实际圆半径

        Dim cf As CircleF

        cf = CvInvoke.MinEnclosingCircle(contourCircle)

        '获得外接圆形

        CvInvoke.Circle(msrc, New Point(CInt(cf.Center.X), CInt(cf.Center.Y)), cf.Radius, New MCvScalar(0, 0, 255), 2)

        CvInvoke.Imshow("m4", msrc)

        '实际圆半径

        Dim realradius1 As Double

        realradius1 = proportion * cf.Radius

        '以下是基于轮廓周长来计算实际圆半径

        '实际圆周长

        Dim realperimeter As Double = perimeter * proportion

        '图像中的圆半径

        Dim radius As Double

        radius = (perimeter / Math.PI) / 2

        '实际圆半径

        Dim realradius2 As Double

        realradius2 = proportion * radius

        MessageBox.Show("最小外接圆来计算实际圆半径：" & realradius1 & ControlChars.CrLf &

                        "基于轮廓周长来计算实际圆半径：" & realradius2)

    End Sub

输出结果如下图所示：

图6-29 求出的圆半径