halcon-轴断面检测定位

前言

通常情况下轴检测时，通常会检测轴的各个阶段的长度。但是由于各种原因，在轴断面的区域现实不明显，无法正确提取，这时候需要根据轴断面的突出部分进行检测，但是由于部分轴的粗轴和细轴区域的宽度差距相当接近，所以就需要通过另外的处理，将轴的断面进行单独提取

1.halcon程序

* Image Acquisition 01: Code generated by Image Acquisition 01
list_files ('D:/2024Work/work/4.12断层检测', ['files','follow_links'], ImageFiles)
tuple_regexp_select (ImageFiles, ['\\.(tif|tiff|gif|bmp|jpg|jpeg|jp2|png|pcx|pgm|ppm|pbm|xwd|ima|hobj)$','ignore_case'], ImageFiles)
for Index := 0 to |ImageFiles| - 1 by 1
    read_image (Image, ImageFiles[Index])
    *gen_rectangle1 (ROI_0, 707.927, 631.228, 1072.01, 1398.08)
    *gen_rectangle1 (ROI_0, 735.583, 1488.93, 1089.37, 2298.93)
    gen_rectangle1 (ROI_0, 761.136, 824.068, 1021.81, 1090.71)
    reduce_domain (Image, ROI_0, ImageReduced)
    *矫正图像
    binary_threshold (ImageReduced, Region, 'max_separability', 'dark', UsedThreshold)
    erosion_circle (Region, RegionErosion,1)
    smallest_rectangle2 (RegionErosion, Row, Column, Phi2, Length1, Length2)
    vector_angle_to_rigid (Row, Column, Phi2, Row, Column, 1.5707963, HomMat2D)
    affine_trans_image (ImageReduced, ImageAffineTrans, HomMat2D, 'constant', 'false')
    reduce_domain (ImageAffineTrans, ROI_0, ImageReduced2)
    *获取断面下半部分
    binary_threshold (ImageReduced2, Region2, 'smooth_histo', 'dark', UsedThreshold1)
    *将上下断面进行联合
    shape_trans (Region2, RegionTrans, 'convex')
    *取区域的最大内接矩形
    inner_rectangle1 (RegionTrans, Row1, Column1, Row2, Column2)
    *创建矩形区域。根据实际情况，可以对内接矩形的左右列坐标进行修改，使得较粗的部分可以较多
    *较细的部分尽可能的少
    gen_rectangle1 (Rectangle1, Row1, Column1-4, Row2, Column2-3)
    *求取差值，将轴区域与内接矩形区域求差值，可得粗轴外轮廓
    difference (Region2, Rectangle1, RegionDifference)
    *剔除细轴多余的区域
    opening_circle (RegionDifference, RegionOpening1, 1)
    connection (RegionOpening1, ConnectedRegions)
    select_shape (ConnectedRegions, SelectedRegions, 'area', 'and', 200.96, 10000)
    union1 (SelectedRegions, RegionUnion)
    *进行横向闭运算，将粗轴进行相连接。
    closing_rectangle1 (RegionUnion, RegionClosing, 1000000000, 1)
    opening_circle (RegionClosing, RegionOpening, 5)
    connection (RegionOpening, ConnectedRegions1)    
    select_shape (ConnectedRegions1, SelectedRegions1, 'area', 'and', 10833.3, 363400.6)
    shape_trans (SelectedRegions1, RegionTrans1, 'convex')
    *获取断面上半部分
    *做差值提取到断面上半部分，剔除多余的噪点
    difference (Region2, RegionTrans1, RegionDifference2)
    connection (RegionDifference2, ConnectedRegions3)   
    select_shape (ConnectedRegions3, SelectedRegions3, 'area', 'and', 5196.22, 500000)
    opening_circle (SelectedRegions3, RegionOpening2, 3.5)
    *上下部分
    *对上下部分进行膨胀取交集
    dilation_circle (RegionOpening2, RegionDilation2, 1)
    dilation_circle (RegionTrans1, RegionDilation3, 3)
    intersection (RegionDilation2, RegionDilation3, RegionIntersection1)
    *结果显示
    skeleton (RegionIntersection1, Skeleton)
    junctions_skeleton (Skeleton, EndPoints, JuncPoints)
    get_region_points (EndPoints, Rows1, Columns1)
    gen_cross_contour_xld (Cross, Rows1[0], Columns1[0], 50, 0)
    gen_cross_contour_xld (Cross1, Rows1[1], Columns1[1], 50, 0)
    dev_display (ImageAffineTrans)
    dev_display (Cross)
    dev_display (Cross1)
    
stop ()
endfor

2.程序解析

2.1图像矫正

在进行检测前，需要将图像矫正在垂直状态，在几何图像中是最大内接矩形是不容易判断的，例如在圆中的最大内接矩形是无限个。所以halcon并没有提供带角度的最大内接矩形，只是提供了垂直角度下的最大内接矩形。所以，我们需要对图像进行矫正到垂直。

    *矫正图像
    binary_threshold (ImageReduced, Region, 'max_separability', 'dark', UsedThreshold)
    erosion_circle (Region, RegionErosion,1)
    smallest_rectangle2 (RegionErosion, Row, Column, Phi2, Length1, Length2)
    vector_angle_to_rigid (Row, Column, Phi2, Row, Column, 1.5707963, HomMat2D)
    affine_trans_image (ImageReduced, ImageAffineTrans, HomMat2D, 'constant', 'false')
    reduce_domain (ImageAffineTrans, ROI_0, ImageReduced2)

2.2获取断面粗轴部分

本文中以，断面的下半部分作为轴的粗轴部分。首先是将上下断面区域进行联合，取区域的内接矩形，做差值后，就可以提取出粗轴的部分。

*获取断面下半部分
    binary_threshold (ImageReduced2, Region2, 'smooth_histo', 'dark', UsedThreshold1)
    *将上下断面进行联合
    shape_trans (Region2, RegionTrans, 'convex')
    *取区域的最大内接矩形
    inner_rectangle1 (RegionTrans, Row1, Column1, Row2, Column2)
    *创建矩形区域。根据实际情况，可以对内接矩形的左右列坐标进行修改，使得较粗的部分可以较多
    *较细的部分尽可能的少
    gen_rectangle1 (Rectangle1, Row1, Column1-4, Row2, Column2-3)
    *求取差值，将轴区域与内接矩形区域求差值，可得粗轴外轮廓
    difference (Region2, Rectangle1, RegionDifference)
    *剔除细轴多余的区域
    opening_circle (RegionDifference, RegionOpening1, 1)
    connection (RegionOpening1, ConnectedRegions)
    select_shape (ConnectedRegions, SelectedRegions, 'area', 'and', 200.96, 10000)
    union1 (SelectedRegions, RegionUnion)
    *进行横向闭运算，将粗轴进行相连接。
    closing_rectangle1 (RegionUnion, RegionClosing, 1000000000, 1)
    opening_circle (RegionClosing, RegionOpening, 5)
    connection (RegionOpening, ConnectedRegions1)    
    select_shape (ConnectedRegions1, SelectedRegions1, 'area', 'and', 10833.3, 363400.6)
    shape_trans (SelectedRegions1, RegionTrans1, 'convex')

2.3获取断面细轴部分

将整个轴的区域减去粗轴的区域，获得的就是断面区域

    difference (Region2, RegionTrans1, RegionDifference2)
    connection (RegionDifference2, ConnectedRegions3)   
    select_shape (ConnectedRegions3, SelectedRegions3, 'area', 'and', 5196.22, 500000)
    opening_circle (SelectedRegions3, RegionOpening2, 3.5)

2.4获取断面并显示

对上下部分进行膨胀，膨胀后取交集，即可获取到断面区域的。可以根据实际情况条件膨胀的参数使得结果在所需要的位置

*上下部分
    *对上下部分进行膨胀取交集
    dilation_circle (RegionOpening2, RegionDilation2, 1)
    dilation_circle (RegionTrans1, RegionDilation3, 3)
    intersection (RegionDilation2, RegionDilation3, RegionIntersection1)
    *结果显示
    skeleton (RegionIntersection1, Skeleton)
    junctions_skeleton (Skeleton, EndPoints, JuncPoints)
    get_region_points (EndPoints, Rows1, Columns1)
    gen_cross_contour_xld (Cross, Rows1[0], Columns1[0], 50, 0)
    gen_cross_contour_xld (Cross1, Rows1[1], Columns1[1], 50, 0)
    dev_display (ImageAffineTrans)
    dev_display (Cross)
    dev_display (Cross1)

总结

断面检测的核心在于对区域的内接矩形和差值的运用。通过内接矩形剔除细轴的区域，在不断通过差值进行计算即可。