标签：血管 提取 lines gauss halcon 造影 dev true 线条

1.血管造影

血管造影，是一种介入检测方法，将显影剂注入血管里。因为X光无法穿透显影剂，血管造影正是利用这一特性，通过显影剂在X光下所显示的影像来诊断血管病变的

1.1 原始图

1.2 处理后的图

（1）提取血管轮廓线

 （2）直径测量

 

2.halcon实现代码

*把其他显示窗更新关闭，一般放在开头
dev_update_off ()
*读取图片
read_image (Angio, 'angio-part')
*获取图片长宽
get_image_size (Angio, Width, Height)
*关闭窗口
dev_close_window ()
*新增显示窗口句柄
dev_open_window (0, 0, Width, Height, 'black', WindowID)
*设置显示模板
set_display_font (WindowID, 14, 'mono', 'true', 'false')
*显示图片
dev_display (Angio)
*显示字体
disp_message (WindowID, 'Original image', 'window', 12, 12, 'black', 'true')

*在荧幕上显示暂停程序继续操作的信息。
disp_continue_message (WindowID, 'black', 'true')
stop ()

*根据要提取线的最大宽度以及对比度计算出lines_gauss算子输入的Sigma、Low、High值
calculate_lines_gauss_parameters (8, [12,0], Sigma, Low, High)

*检测线条以及其宽度
lines_gauss (Angio, Lines, Sigma, Low, High, 'dark', 'true', 'parabolic', 'true')

* 选择特定形状特征要求的xld轮廓或多边形
select_contours_xld (Lines, RelLines, 'length', 5.0, 999, 0, 0)

*设置多个输出颜色
dev_set_colored (12)

*显示图片
dev_display (Angio)

*显示检测的线条
dev_display (RelLines)

*显示消息
disp_message (WindowID, 'Extracted lines', 'window', 12, 12, 'black', 'true')

*在荧幕上显示暂停程序继续操作的信息。
disp_continue_message (WindowID, 'black', 'true')
stop ()

*count_obj算子是用来计算输入区域中连通域的个数。
count_obj (RelLines, Number)

*对xld进行排序， 以外接矩形的相应顶点（左上，右上，左下，右下）为排序点，true时升序，false降序
sort_contours_xld (RelLines, RelLines, 'lower_left', 'true', 'row')

*创建一个空的object
gen_empty_obj (PrintedLines)

*遍历所有检测的线
for I := 1 to Number by 1
    dev_display (Angio)
    dev_set_color ('white')
    dev_display (PrintedLines)
    
    *选择索引
    select_obj (RelLines, Line, I)
    dev_set_color ('green')
    dev_display (Line)
    
    *把两个区域连接在一起
    concat_obj (PrintedLines, Line, PrintedLines)

    *返回XLD轮廓(contour)的坐标
    get_contour_xld (Line, Row, Col)
    meanRow := sum(Row) / |Row|
    meanCol := sum(Col) / |Col|
    get_contour_attrib_xld (Line, 'width_left', WidthL)
    get_contour_attrib_xld (Line, 'width_right', WidthR)
    get_contour_attrib_xld (Line, 'contrast', Contrast)
    
    *
    meanContrast := sum(Contrast) / |Contrast|
    * To display the lines widths, the point at which the gray value drops
    * to 25% of the contrast between the line and the background will be
    * displayed.  This point is given by sqrt(3/4) for the parabolic line
    * model.
    Diameter := (WidthL + WidthR) * sqrt(0.75)
    Diam := sum(Diameter) / |Diameter|
    if (meanRow > Height - 50)
        disp_message (WindowID, 'diam: ' + Diam, 'image', Height - 70, meanCol, 'yellow', 'false')
    else
        disp_message (WindowID, 'diam: ' + Diam, 'image', meanRow, meanCol, 'yellow', 'false')
    endif
    disp_message (WindowID, 'Diameters of line segments', 'window', 135, 60, 'black', 'true')
    if (I < 5)
        disp_continue_message (WindowID, 'black', 'true')
        stop ()
    else
        wait_seconds (0.2)
    endif
endfor
dev_update_on ()

3.核心函数讲解

3.1 calculate_lines_gauss_parameters

calculate_lines_gauss_parameters( : : MaxLineWidth, Contrast : Sigma, Low, High)

（1）功能

根据要提取线的最大宽度以及对比度计算出lines_gauss算子输入的Sigma、Low、High值。

（2）参数列表

MaxLineWidth (input_control) ： lines_gauss要提取线条的最大宽度
Contrast (input_control) ：lines_gauss要提取线的对比度。Contrast 值不仅可以一个，也可以为两个。当值为两个时，数组中的第二个值是要提取线的最小对比度，并且其值不能大于第一个值。比如：[20,10]
Sigma (output_control) ：获取用于lines_gauss输入的Sigma值
Low (output_control)：获取用于lines_gauss输入的Low 值
High (output_control) ：获取用于lines_gauss输入的High 值

（3）详解

对于参数Contrast 需要说明的是，如果只选择一个值时，最小对比度将会默认为最大对比度的1/3，最小对比度越小，线条将会延伸到对比度较低的区域，即线条越长。反之，值越高，线条越短，但越突出。
lines_gauss算子用到的滞后阈值方法中的Low、High值是根据线的最大宽度以及对比度（两个高低对比度）计算得出的。

3.2 lines_gauss

（1）原型

lines_gauss(Image : Lines : Sigma, Low, High, LightDark, ExtractWidth, LineModel, CompleteJunctions : )

功能
检测线条以及其宽度。

（2）参数列表

Image (input_object) ：输入图像
Lines (output_object) ：检测线条（XLD）
Sigma (input_control) ：高斯滤波值
Low (input_control) ：滞后阈值分割的低阈值
High (input_control) ：滞后阈值分割的高阈值
LightDark (input_control) ：提取线条类型（ ‘dark’,‘light’）
ExtractWidth (input_control) ：是否提取线宽（‘false’,‘true’）
LineModel (input_control) ：用来调整线条位置和宽度的线模型（‘bar-shaped’, ‘gaussian’, ‘none’, ‘parabolic’）
CompleteJunctions (input_control) ：在断连的部分是否添加节点使线条连续（‘false’, ‘true’）

（3）详解

lines_gauss 主要功能是提取图像上的线条，提取的结果属于亚像素精度的XLD轮廓。

参数LightDark 其实就是让你选择你要提取线条的类型，暗色还是亮色，选’dark’, ‘light’。

参数ExtractWidth是否提取每条XLD轮廓线的线宽，true提取，false不提取。

参数LineModel 如果选择 除’none’以外的参数，lines_gauss算子会补偿非对称线条（即在线条的中心两侧有不同对比度的线条），来校正提取出的线条的位置和宽度。用于校准的线模型有’bar-shaped’, ‘gaussian’, ‘none’, ‘parabolic’四种。线模型’bar-shaped’适用于大多数应用；如果管状物体采用背光方式打光（比如血管的X光图像），提取线条可以使用’gaussian’, 'parabolic’模型；’parabolic’模型常用来提取边缘比较锐利的线条（比如采用背光照明图像中的线条）；’gaussian’则用来提取边缘相对不那么锐利的线条。只有将ExtractWidth设置为true时，参数LineModel才有意义。

由于微分几何算法的原因，线提取器（lines_gauss）不能提取出所有线段结合点，当CompleteJunctions 设置为true时，lines_gauss试图通过不同算法提取出那些被遗漏的节点。
 

标签：血管,提取,lines,gauss,halcon,造影,dev,true,线条
来源： https://blog.csdn.net/thehappysheep/article/details/120755120

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。