circles.hdev

read_image (Image, 'double_circle') * * Init window dev_close_window () get_image_size (Image, Width, Height) dev_open_window (0, 0, Width, Height, 'black', WindowHandle) * * Segment a region containing the edges *快速阈值分割 fast_threshold (Image, Region, 0, 120, 7) *使用形态学运算来计算区域的边界，这里就是将图片边界提取出来 boundary (Region, RegionBorder, 'inner') *在原区域的基础上做一个最小外接矩形得到区域1，然后让刚才那个矩形上下左右都以5的像素内缩，得到矩形2， *现在形成了一个矩形环，用区域1减去这个矩形环，得到剩下的区域就是目标区域 clip_region_rel (RegionBorder, RegionClipped, 5, 5, 5, 5) *使用半径为2.5的圆形元素对区域进行膨胀 dilation_circle (RegionClipped, RegionDilation, 2.5) *分离出目标区域（从原图中分离），不是分割，减小了图像的定义域，减小了计算 reduce_domain (Image, RegionDilation, ImageReduced) * * In the subdomain of the image containing the edges, * extract subpixel precise edges. *过滤器提取亚像素精确边缘 edges_sub_pix (ImageReduced, Edges, 'canny', 2, 20, 60) *分割XLD轮廓为线段和圆弧或椭圆弧 *segment_contours_xld(Contours : ContoursSplit : *Mode, SmoothCont, MaxLineDist1, MaxLineDist2 : ) *Contours :输入轮廓；ContoursSplit：分割后的轮廓；Mode：轮廓分割模式； *SmoothCont：用于平滑轮廓的点数；MaxLineDist1：轮廓与近似线之间的最大距离（第一次迭代）； *MaxLineDist2 ：轮廓与近似线之间的最大距离（第二次迭代） segment_contours_xld (Edges, ContoursSplit, 'lines_circles', 5, 4, 3) *将目标区域分开成单个单个的区域并标记序号 count_obj (ContoursSplit, Number) dev_display (Image) dev_set_draw ('margin') dev_set_color ('white') dev_update_window ('off') for I := 1 to Number by 1 *拿出其中一个区域 select_obj (ContoursSplit, ObjectSelected, I) *返回XLD轮廓的全局属性值’cont_approx’表示:对于’cont_approx’=-1，轮廓近似线段; *‘cont_approx’=0，曲线近似椭圆;‘cont_approx’=1，曲线近似圆弧。 get_contour_global_attrib_xld (ObjectSelected, 'cont_approx', Attrib) * Fit a circle to the line segment that are arcs of a circle *大于0表示为圆弧段，用圆弧拟合 if (Attrib > 0) *拟合圆 fit_circle_contour_xld (ObjectSelected, 'ahuber', -1, 2, 0, 3, 2, Row, Column, Radius, StartPhi, EndPhi, PointOrder) *将拟合的圆画出来 gen_circle_contour_xld (ContCircle, Row, Column, Radius, 0, rad(360), 'positive', 1.0) dev_display (ContCircle) endif endfor dev_set_colored (12) dev_set_line_width (3) dev_display (ContoursSplit)