1.示例图
(1)原始图
(2)结果图
2
2.halcon 代码
* vessel.hdev: Segmentation & Measurement of a blood vessel
*
dev_update_window ('off')
dev_close_window ()
dev_open_window (0, 0, 512, 512, 'black', WindowID)
set_display_font (WindowID, 14, 'mono', 'true', 'false')
read_image (Image, 'vessel')
dev_display (Image)
disp_continue_message (WindowID, 'black', 'true')
stop ()
dev_set_draw ('fill')
dev_set_color ('green')
dev_set_line_width (3)
texture_laws (Image, Texture, 'el', 2, 5)
dev_display (Texture)
disp_continue_message (WindowID, 'black', 'true')
stop ()
mean_image (Texture, Energy, 211, 61)
dev_display (Energy)
disp_continue_message (WindowID, 'black', 'true')
stop ()
dev_display (Image)
binary_threshold (Energy, Vessel, 'smooth_histo', 'dark', UsedThreshold)
dev_display (Vessel)
disp_continue_message (WindowID, 'black', 'true')
stop ()
dev_set_draw ('margin')
get_region_contour (Vessel, VesselRows, VesselColumns)
gen_region_points (VesselBoundary, VesselRows, VesselColumns)
dev_display (Image)
clip_region (VesselBoundary, Boundaries, 0, 2, 511, 509)
dev_display (Boundaries)
smallest_rectangle1 (Boundaries, UpperRow, Column1, LowerRow, Column2)
disp_continue_message (WindowID, 'black', 'true')
stop ()
*
* Measurement of Diameter
*
Diameter := []
dev_set_color ('red')
dev_update_var ('off')
dev_update_time ('off')
dev_update_pc ('off')
for i := 2 to 509 by 1
gen_region_line (ScanLine, UpperRow, i, LowerRow, i)
dev_display (ScanLine)
intersection (Boundaries, ScanLine, BoundaryPoints)
get_region_points (BoundaryPoints, Rows, Columns)
Diameter := [Diameter,abs(Rows[|Rows| - 1] - Rows[0])]
endfor
dev_update_pc ('on')
dev_update_time ('on')
dev_update_var ('on')
Diameter := Diameter
*
* statistics
*
meanDiameter := sum(Diameter) / |Diameter|
minDiameter := min(Diameter)
dev_set_color ('green')
dev_display (Image)
dev_display (Boundaries)
dev_set_color ('yellow')
get_image_size (Image, Width, Height)
get_string_extents (WindowID, 'meanDiam: ' + meanDiameter + ', minDiameter: ' + minDiameter, Ascent, Descent, TxtWidth, TxtHeight)
disp_message (WindowID, 'meanDiam: ' + meanDiameter + ', minDiameter: ' + minDiameter, 'image', (UpperRow + LowerRow) / 2, (Width - TxtWidth) / 2, 'yellow', 'false')
* get_mbutton (WindowID, Row, Column, Button)
dev_set_draw ('fill')
dev_set_line_width (1)
dev_update_window ('on')