例程功能

 提取连续的边缘

代码如下

dev_update_off ()
dev_close_window ()
* ****
* step: acquire image
* ****
read_image (Image, 'mreut')
get_image_size (Image, Width, Height)
dev_open_window_fit_image (Image, 0, 0, Width, Height, WindowID)
set_display_font (WindowID, 12, 'mono', 'true', 'false')
dev_set_draw ('margin')
dev_set_line_width (3)
dev_display (Image)
disp_continue_message (WindowID, 'black', 'true')
stop ()
* ****
* step: filter image
* ****
edges_image (Image, ImaAmp, ImaDir, 'lanser2', 0.5, 'nms', 20, 40)
dev_display (ImaAmp)
disp_continue_message (WindowID, 'black', 'true')
stop ()
* ****
* step: extract edges
* ****
threshold (ImaAmp, Region, 1, 255)
connection (Region, ConnectedRegions)
dev_clear_window ()
dev_set_colored (12)
dev_display (ConnectedRegions)
disp_continue_message (WindowID, 'black', 'true')
stop ()
* ****
* step: process edges
* ****
dev_clear_window ()
count_obj (ConnectedRegions, Number)
gen_empty_obj (XLDContours)
for i := 1 to Number by 1
    select_obj (ConnectedRegions, SingleEdgeObject, i)
    split_skeleton_lines (SingleEdgeObject, 2, BeginRow, BeginCol, EndRow, EndCol)
    for k := 0 to |BeginRow| - 1 by 1
        gen_contour_polygon_xld (Contour, [BeginRow[k],EndRow[k]], [BeginCol[k],EndCol[k]])
        concat_obj (XLDContours, Contour, XLDContours)
    endfor
endfor
dev_display (XLDContours)

要点

  1. 明确边缘检测的一般步骤:获取图像 ~> 设置ROI ~> 图像滤波 ~> 提取边缘 ~> 处理边缘 ~> 结果可视化;
  2. 基于边缘振幅和方向对图像进行分割;
	read_image (Image, 'mreut') # 读图
	edges_image (Image, ImaAmp, ImaDir, 'lanser2', 0.5, 'nms', 20, 40) # 阈值分割
	threshold (ImaAmp, Region, 1, 255) # 灰度值分割
	connection (Region, ConnectedRegions) # 打散

edges_image——使用 Deriche、Lanser、Shen 或 Canny 过滤器提取边缘;
Image ——入参,待处理图像;
ImaAmp ——出参,边缘幅度(梯度)分割结果;
ImaDir ——出参,边缘方向分割结果;
Filter ——入参,过滤器;
Alpha ——入参,滤镜参数(值越小,平滑效果越强,细节越少,canny算子相反);
NMS ——入参,极大值抑制;
Low ——入参,滞后阈值分割的低阈值(不需要则置negative);
High ——入参,滞后阈值分割的高阈值(不需要则置negative);
3. 将区域数据转换为数值(用线段近似表示区域)

	count_obj (ConnectedRegions, Number)
	gen_empty_obj (XLDContours)
	for i := 1 to Number by 1
		select_obj (ConnectedRegions, SingleEdgeObject, i)
		split_skeleton_lines (SingleEdgeObject, 2, BeginRow, BeginCol, EndRow, \
		EndCol)
		for k := 0 to |BeginRow| - 1 by 1
			gen_contour_polygon_xld (Contour, [BeginRow[k],EndRow[k]], \
			[BeginCol[k],EndCol[k]])
			concat_obj (XLDContours, Contour, XLDContours)
		endfor
	endfor
	dev_display (XLDContours)

 注意在for循环外创建了空对象XLDContours,在循环内部对区域进行逐个处理,将处理的结果通过concat_obj放入XLDContours中,处理结束统一显示,这种处理形式很有意思,挺常用的;
split_skeleton_lines——用一条一像素宽的无分支的线分割;
SkeletonRegion ——入参,输入的线;
MaxDistance ——入参,线指向连接两个端点的线段的最大距离(没太明白);
BeginRow ——出参,输出直线的起始行值(纵坐标);
BeginCol ——出参,输出直线的起始列值(横坐标);
EndRow ——出参,输出直线的终止行值(纵坐标);
EndCol ——出参,输出直线的终止列值(横坐标);

01-25 11:52