效果图

注意内部代码中有stop() 运行停止

可以将其全部注释，在不注释的情况下每运行一次，只能输出一次操作，需要一直点击运行，直到达到最终效果

dev_update_window ('off')
dev_close_window ()
*打开一个原点（左上角）在（0,0）的，宽：高：，黑色背景的图形窗口
dev_open_window (0, 0, 512, 512, 'black', WindowID)
read_image (Image, 'D:/Besktop/faimage/测试原图/23_50_12.bmp')
 
dev_display (Image)
set_display_font (WindowID, 22, 'mono', 'true', 'false')
disp_continue_message (WindowID, 'black', 'true')
stop ()
threshold (Image, Bright, 90, 255)
*画出最小外接矩形
shape_trans (Bright, rigontran, 'rectangle2')
dev_set_color ('green')
dev_set_line_width (2)
dev_set_draw ('margin')
dev_display (rigontran)
disp_continue_message (WindowID, 'black', 'true')
stop ()
*获得矩形内图像
reduce_domain (Image, rigontran, DieGrey)
threshold (DieGrey, Wires, 0, 50)
*填充区域内给定特征的孔
fill_up_shape (Wires, WiresFilled, 'area', 1, 100)
dev_display (Image)
*设置所画图形的形状，fill为填满
dev_set_draw ('margin')
dev_set_color ('red')
dev_display (WiresFilled)
disp_continue_message (WindowID, 'black', 'true')
stop ()
*开运算，除去孤立的小点、毛刺，总位置大小形状不变；输入控制半径
opening_circle (WiresFilled, Balls, 15.5)
dev_set_color ('green')
dev_display (Balls)
disp_continue_message (WindowID, 'black', 'true')
stop ()
connection (Balls, SingleBalls)
*根据特征选择区域；‘and’全部特征，‘or’特征之一；这里通过圆率进行筛选
select_shape (SingleBalls, IntermediateBalls, 'circularity', 'and', 0.85, 1.0)
*根据区域相关位置进行排序，这里选择‘column’，按照先列后行进行排序；递增
sort_region (IntermediateBalls, FinalBalls, 'first_point', 'true', 'column')
dev_display (Image)
dev_set_color('orange')
dev_display (FinalBalls)
disp_continue_message (WindowID, 'black', 'true')
stop ()
*计算给定区域的最小外接圆
smallest_circle (FinalBalls, Row, Column, Radius)
NumBalls := |Radius|
Diameter := 2 * Radius
meanDiameter := sum(Diameter) / NumBalls
mimDiameter := min(Diameter)
dev_display (Image)
*显示描述了圆心坐标及半径的圆
disp_circle (WindowID, Row, Column, Radius)
dev_set_color ('white')
 
*求数量用||
Row_length:=|Row|
Col_lenght:=|Column|
*遍历所有的圆心坐标
for i := 0 to Row_length-1 by 1 
    *圆心坐标
    R:=[Row[i],Column[i]]
    *用十字画出圆心位置
    disp_cross (WindowID, R[0], R[1], 10, 10)
    *disp_message(WindowID, 'Cc: ' + R, 'image',R[0],R[1], 'white', 'false')
endfor
stop()
 
for i := 0 to Row_length-1 by 1
    *圆心坐标
    R:=[Row[i],Column[i]]
    R1:=Row[i]
    W1:=Column[i]
    disp_message(WindowID, '(' + W1 +','+ R1 + ')', 'image',R[0],R[1], 'white', 'false')
    
endfor
stop()
 
for i := 0 to Row_length-1 by 1
    *圆心坐标
    R:=[Row[i],Column[i]]
    for k := 1 to Row_length-1 by 1
    W:=[Row[k],Column[k]]
    distance_pp(R[0],R[1],W[0],W[1], Distance)
    disp_line(WindowID,R[0],R[1],W[0],W[1])
    disp_message (WindowID, 'Dp: ' + Distance, 'image',(R[0]+W[0])/2,(R[1]+W[1])/2, 'white', 'false')
    endfor
endfor
 
 
for i := 1 to NumBalls by 1
    *判断余数是否为1
    if (fmod(i,2) != 1)
        disp_message (WindowID, 'D: ' + Diameter[i - 1], 'image', Row[i - 1] - 1.8 * Radius[i - 1], max([Column[i - 1] - 60,0]), 'white', 'false')
        *disp_message (WindowID, 'D: ' + Diameter[i - 1], 'image', Row[i - 1] - 2.0 * Radius[i - 1], max([Column[i - 1] - 60,0]), 'white', 'false')
    else
        disp_message (WindowID, 'D: ' + Diameter[i - 1], 'image', Row[i - 1] + 1.2 * Radius[i - 1], max([Column[i - 1] - 60,0]), 'white', 'false')
    endif
endfor
 
dev_set_color ('green')
dev_update_window ('on')
*disp_continue_message (WindowID, 'black', 'true')
stop ()
dev_close_window ()