求助大佬，焊接轨迹同心度测量

如上图，想要求左边二个焊接轨迹圆心的距离，请求大佬指点
下面是想要的效果图

1. *1.裁切roi
   
2. read_image (Image1, 'C:/Users/Annan/Downloads/焊接/557.png')
   
3. * read_image (Image1, 'C:/Users/Annan/Downloads/焊接/456.png')
   
4. * read_image (Image1, 'C:/Users/Annan/Downloads/焊接/253.png')
   
5. gen_circle (ROI_0, 273.938, 236.098, 165.243)
   
6. area_center (ROI_0, Area, RowC, ColumnC)
   
7. reduce_domain (Image1, ROI_0, ImageReduced)
   
8. gen_cross_contour_xld (Cross, RowC, ColumnC, 36, 0.785398)
   
9. dev_display (Image1)
   
10. dev_set_draw ('margin')
    
11. dev_display (Cross)
    
12. dev_display (ROI_0)
    
13. stop ()
    
14. *2.预处理
    
15. rgb1_to_gray (ImageReduced, GrayImage)
    
16. mean_image (GrayImage, ImageMean, 12, 12)
    
17. *3.动态阈值取边缘
    
18. dyn_threshold (GrayImage, ImageMean, RegionDynThresh, 20, 'dark')
    
19. connection (RegionDynThresh, ConnectedRegions2)
    
20. select_shape (ConnectedRegions2, SelectedRegions, 'outer_radius', 'and', 30, 200)
    
21. select_shape (SelectedRegions, SelectedRegions1, 'width', 'and', 10, 99999)
    
22. union1 (SelectedRegions1, RegionUnion)
    
23. *4.连接边缘
    
24. dilation_circle (RegionUnion, RegionDilation, 3.5)
    
25. *5.展开坐标系 去除毛刺干扰
    
26. area_center (RegionDilation, Area1, Row, Column)
    
27. get_image_size (ImageReduced, Width, Height)
    
28. OuterRadius:=200
    
29. InnerRadius:=0
    
30. polar_trans_region (RegionDilation, PolarTransRegion, Row, Column, 0, 6.28319, 0, OuterRadius, Width, Height, 'nearest_neighbor')
    
31. opening_circle (PolarTransRegion, RegionOpening, 7)
    
32. *6.逆变换
    
33. polar_trans_region_inv (RegionOpening, XYTransRegion, Row, Column, 0, 6.28319, 0, OuterRadius, Width, Height, Width, Height, 'nearest_neighbor')
    
34. *7.提取边缘区域骨架
    
35. connection (XYTransRegion, ConnectedRegions)
    
36. skeleton (ConnectedRegions, Skeleton)
    
37. *8.区域转轮廓
    
38. gen_contour_region_xld (Skeleton, Contours, 'center')
    
39. *9.轮廓拟合椭圆
    
40. fit_ellipse_contour_xld (Contours, 'fitzgibbon', -1, 0, 0, 200, 3, 2, Row1, Column1, Phi, Radius1, Radius2, StartPhi, EndPhi, PointOrder)
    
41. gen_ellipse (Ellipse, Row1, Column1, Phi, Radius1, Radius2)
    
42. *8.筛选两个圆
    
43. select_shape (Ellipse, SelectedRegions1, 'outer_radius', 'and', 95, 105)
    
44. select_shape (Ellipse, SelectedRegions2, 'outer_radius', 'and', 70, 80)
    
45. concat_obj (SelectedRegions1, SelectedRegions2, ObjectsConcat)
    
46. *9.求圆心并显示
    
47. area_center (ObjectsConcat, Area2, Row2, Column2)
    
48. gen_cross_contour_xld (Cross1, Row2, Column2, 36, 0.785398)
    
49. dev_display (Image1)
    
50. dev_set_draw ('margin')
    
51. dev_display (ObjectsConcat)
    
52. dev_display (Cross1)

复制代码

效果图

效果图

效果图

易安呀 发表于 2024-4-10 14:39

易安呀 发表于 2024-4-10 14:39

大佬好，mean_image算子中的12这个参数是根据哪里算出来的，不太明白

易安呀 发表于 2024-4-10 14:39

竟然还用到了极坐标变换，而且还是用区域转换的轮廓再拟合，搞得有点复杂了，这种图片直接用二维测量才是最佳方法

halcon2016 发表于 2024-4-11 13:20
竟然还用到了极坐标变换，而且还是用区域转换的轮廓再拟合，搞得有点复杂了，这种图片直接用二维测量才是 ...

通用算法 主要是思路 测量在halcon里好用 在opencv你不一定好实现。至于为什么要用到极坐标变换 是为了去除杂质 区域是圆形的 开运算会破坏原有的轮廓 所以需要转换后 在去除

aicanjia 发表于 2024-4-11 08:58
大佬好，mean_image算子中的12这个参数是根据哪里算出来的，不太明白

看你边缘的内宽大概是多少像素 一般设置为一半以上

学习一下

1. *1.裁切roi
   
2. 
   
3. read_image (Image1, './1.png')
   
4. 
   
5. * read_image (Image1, 'C:/Users/Annan/Downloads/焊接/557.png')
   
6. * read_image (Image1, 'C:/Users/Annan/Downloads/焊接/456.png')
   
7. * read_image (Image1, 'C:/Users/Annan/Downloads/焊接/253.png')
   
8. gen_circle (ROI_0, 273.938, 236.098, 165.243)
   
9. area_center (ROI_0, Area, RowC, ColumnC)
   
10. reduce_domain (Image1, ROI_0, ImageReduced)
    
11. gen_cross_contour_xld (Cross, RowC, ColumnC, 36, 0.785398)
    
12. dev_display (Image1)
    
13. dev_set_draw ('margin')
    
14. dev_display (Cross)
    
15. dev_display (ROI_0)
    
16. stop ()
    
17. *2.预处理
    
18. rgb1_to_gray (ImageReduced, GrayImage)
    
19. mean_image (GrayImage, ImageMean, 12, 12)
    
20. *3.动态阈值取边缘
    
21. dyn_threshold (GrayImage, ImageMean, RegionDynThresh, 20, 'dark')
    
22. connection (RegionDynThresh, ConnectedRegions2)
    
23. select_shape (ConnectedRegions2, SelectedRegions, 'outer_radius', 'and', 30, 200)
    
24. select_shape (SelectedRegions, SelectedRegions1, 'width', 'and', 10, 99999)
    
25. union1 (SelectedRegions1, RegionUnion)
    
26. *4.连接边缘
    
27. dilation_circle (RegionUnion, RegionDilation, 3.5)
    
28. 
    
29. *--------------方法1：blob------------------------
    
30. skeleton (RegionDilation, Skeleton1)
    
31. gen_contour_region_xld (Skeleton1, Contours1, 'center')
    
32. select_obj (Contours1, ObjectSelected, 1)
    
33. segment_contours_xld (ObjectSelected, ContoursSplit, 'lines_ellipses', 5, 4, 1)
    
34. select_shape_xld (ContoursSplit, SelectedXLD1, ['circularity','area'], 'and', [0.24704,3792.88], [1,50000])
    
35. 
    
36. fit_ellipse_contour_xld (SelectedXLD1, 'fhuber', -1, 0, 0, 200, 3, 2, \
    
37.                          Row1, Column1, Phi, Radius1, Radius2, StartPhi, EndPhi, PointOrder)
    
38. 
    
39. gen_ellipse_contour_xld (ContEllipse, Row1, Column1, Phi, Radius1, Radius2, \
    
40.                          StartPhi, EndPhi, PointOrder, 1.5)
    
41. 
    
42. select_shape_xld (ContEllipse, SelectedXLD, 'area', 'and', 25658.8, 50000)
    
43. 
    
44. 
    
45. select_obj (Contours1, ObjectSelected1, 2)
    
46. fit_ellipse_contour_xld (ObjectSelected1, 'fhuber', -1, 0, 0, 200, 3, 2, \
    
47.                          Row1, Column1, Phi, Radius1, Radius2, StartPhi, EndPhi, PointOrder)
    
48. 
    
49. gen_ellipse_contour_xld (ContEllipse1, Row1, Column1, Phi, Radius1, Radius2, \
    
50.                          StartPhi, EndPhi, PointOrder, 1.5)
    
51. *-----------------------------------------------------------------
    
52. 
    
53. 
    
54. *--------------------------方法2---------------------------
    
55. *5.展开坐标系 去除毛刺干扰
    
56. area_center (RegionDilation, Area1, Row, Column)
    
57. get_image_size (ImageReduced, Width, Height)
    
58. OuterRadius:=200
    
59. InnerRadius:=0
    
60. polar_trans_region (RegionDilation, PolarTransRegion, Row, Column, 0, 6.28319, 0, OuterRadius, Width, Height, 'nearest_neighbor')
    
61. opening_circle (PolarTransRegion, RegionOpening, 7)
    
62. *6.逆变换
    
63. polar_trans_region_inv (RegionOpening, XYTransRegion, Row, Column, 0, 6.28319, 0, OuterRadius, Width, Height, Width, Height, 'nearest_neighbor')
    
64. *7.提取边缘区域骨架
    
65. connection (XYTransRegion, ConnectedRegions)
    
66. skeleton (ConnectedRegions, Skeleton)
    
67. *8.区域转轮廓
    
68. gen_contour_region_xld (Skeleton, Contours, 'center')
    
69. *9.轮廓拟合椭圆
    
70. fit_ellipse_contour_xld (Contours, 'fitzgibbon', -1, 0, 0, 200, 3, 2, Row1, Column1, Phi, Radius1, Radius2, StartPhi, EndPhi, PointOrder)
    
71. gen_ellipse (Ellipse, Row1, Column1, Phi, Radius1, Radius2)
    
72. *8.筛选两个圆
    
73. select_shape (Ellipse, SelectedRegions1, 'outer_radius', 'and', 95, 105)
    
74. select_shape (Ellipse, SelectedRegions2, 'outer_radius', 'and', 70, 80)
    
75. concat_obj (SelectedRegions1, SelectedRegions2, ObjectsConcat)
    
76. 
    
77. *--------------------------------------------------------
    
78. 
    
79. *9.求圆心并显示
    
80. area_center (ObjectsConcat, Area2, Row2, Column2)
    
81. gen_cross_contour_xld (Cross1, Row2, Column2, 36, 0.785398)
    
82. dev_display (Image1)
    
83. dev_set_draw ('margin')
    
84. dev_display (ObjectsConcat)
    
85. dev_display (Cross1)

复制代码