一种冠脉图像造影方法和装置的制造方法
【专利摘要】本发明提供了一种冠脉图像造影方法和装置，所述方法包括：获取病人的冠脉计算机断层摄影血管造影术CTA图像；对得到的冠脉CTA图像进行三维重建，得到三维冠脉图像；获取医生在三维冠脉图像中选择的感兴趣血管段；通过计算确定感兴趣血管段的最佳造影角度组合；根据最佳造影角度组合对三维冠脉图像进行造影，得到最佳造影角度组合下的冠脉造影图像，并将得到的冠脉造影图像进行展示。通过本发明可以提高冠脉造影的准确性。
【专利说明】
-种冠脉图像造影方法和装置
技术领域
[0001] 本发明设及图像处理领域，具体而言，设及一种冠脉图像造影方法和装置。
【背景技术】
[0002] 目前，医学影像是指为了医疗或医学研究，对人体或人体某部分，W非侵入方式取 得的人体内部组织影像。现有的医学影像主要包括计算机断层摄影血管造影术(Computed Tomographic Angiography，CTA)和数字减影血管造影技术（Digital Subtraction Angiography, DSA)等。当医生在给病人的冠状动脉（简称冠脉）病症进行诊断治疗的过程 中，会先获取病人的冠脉造影图像，并通过得到的冠脉造影图像来对病人的冠脉病症进行 诊断治疗。
[0003] 现有的冠脉造影图像的造影过程是利用DSA，先将作为造影剂的舰对比剂注入冠 脉开口，依靠舰对比剂流经的方向和充盈的冠脉，然后借助X线辅助采集病人的冠脉中被舰 对比剂充盈的冠脉的造影图像。
[0004] 由于病人的屯、脏冠脉形状不同，且冠脉的走线复杂，在利用默认体位或者医生经 验选择的造影体位对病人冠脉进行造影的过程中，仅通过一次造影往往无法充分暴露病人 冠脉的病变位置，只能对病人的冠脉进行多次造影后才能明确病人冠脉的病变位置，造成 冠脉造影的准确性较差，而且对病人进行多次冠脉造影，会对病人增加过多的射线福射量 和对比剂用量，给病人的身体带来不必要的损害。
【发明内容】
[0005] 有鉴于此，本发明实施例的目的在于提供一种冠脉图像造影方法和装置，W提高 冠脉造影的准确性。
[0006] 第一方面，本发明实施例提供了一种冠脉图像造影方法，包括：
[0007] 获取病人的冠脉计算机断层摄影血管造影术CTA图像；
[000引对得到的所述冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图像；
[0009] 获取医生在所述Ξ维冠脉图像中选择的感兴趣血管段；
[0010] 通过计算确定所述感兴趣血管段的最佳造影角度组合，所述最佳造影角度组合包 括最佳造影主角度和最佳造影副角度；
[0011] 根据所述最佳造影主角度和所述最佳造影副角度对所述Ξ维冠脉图像进行造影， 得到所述最佳造影角度组合下的冠脉造影图像，并将得到的所述冠脉造影图像进行展示。
[0012] 结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，所 述对得到的所述冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图像，包括：
[0013] 对得到的所述冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图像；
[0014] 从得到的所述Ξ维冠脉图像中，确定出冠脉CTA图像中屯、线；
[0015] 将确定出所述冠脉CTA图像中屯、线的所述Ξ维冠脉图像展示给所述医生。
[0016] 结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所 述从得到的所述Ξ维冠脉图像中，确定出冠脉CTA图像中屯、线，包括：
[0017] 利用霍夫圆检测算法，从所述冠脉CTA图像中确定出升主动脉图像的中屯、；
[0018] 根据确定的所述升主动脉图像的中屯、，从所述冠脉CTA图像中提取出冠脉图像；
[0019] 从提取出的所述冠脉图像中，确定出冠脉CT姻像中屯、线。
[0020] 结合第一方面，本发明实施例提供了第一方面的第Ξ种可能的实施方式，其中，确 定所述感兴趣血管段的最佳造影角度组合，所述最佳造影角度组合包括最佳造影主角度和 最佳造影副角度，包括：
[0021] 随机生成预设数量的造影角度组合，其中，每个所述造影角度组合均包括一造影 主角度和一造影副角度；
[0022] 根据预设的迭代次数，通过遗传算法对随机生成的所述造影角度组合进行迭代计 算，得到多种造影角度组合下所述感兴趣血管段的模拟冠脉造影图像面积值；
[0023] 确定所述感兴趣血管段的最大模拟冠脉造影图像面积值对应的造影角度组合作 为所述感兴趣血管段的最佳造影角度组合。
[0024] 结合第一方面，本发明实施例提供了第一方面的第四种可能的实施方式，其中，根 据预设的迭代次数，通过遗传算法对随机生成的所述造影角度组合进行迭代计算，得到多 种造影角度组合下所述感兴趣血管段的冠脉造影图像面积值，包括：
[0025] 根据所述造影角度组合中记录的造影主角度和造影副角度，对所述Ξ维冠脉图像 进行模拟造影，分别得到所述造影角度组合下所述Ξ维冠脉图像中所述感兴趣血管段和其 他血管段的模拟造影图像区域；
[00%]根据公式
确定所述感兴趣血管 段的模拟造影图像区域中各像素点的参考值f(i，j)，其中，select(i，j) = l表示像素点在 所述感兴趣血管段上，unselect (i，j ) = 1表示像素点在其他血管段上，i表示所述感兴趣血 管段的模拟造影图像区域中各像素点的横坐标，j表示所述感兴趣血管段的模拟造影图像 区域中各像素点的纵坐标；
[0027] 根据确定出的所述感兴趣血管段的模拟造影图像区域中各像素点的参考值，得到 所述造影角度组合下所述感兴趣血管段的冠脉造影图像面积值。
[0028] 第二方面，本发明实施例提供了一种冠脉图像造影装置，包括：
[0029] 图像获取模块，用于获取病人的冠脉计算机断层摄影血管造影术CTA图像；
[0030] Ξ维重建模块，用于对得到的所述冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图 像；
[0031] 血管获取模块，用于获取医生在所述Ξ维冠脉图像中选择的感兴趣血管段；
[0032] 计算模块，用于通过计算确定所述感兴趣血管段的最佳造影角度组合，所述最佳 造影角度组合包括最佳造影主角度和最佳造影副角度；
[0033] 造影模块，用于根据所述最佳造影主角度和所述最佳造影副角度对所述Ξ维冠脉 图像进行造影，得到所述最佳造影角度组合下的冠脉造影图像，并将得到的所述冠脉造影 图像进行展示。
[0034] 结合第二方面，本发明实施例提供了第二方面的第一种可能的实施方式，其中，所 述Ξ维重建模块，包括：
[0035] Ξ维重建单元，用于对得到的所述冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图 像；
[0036] 图像中屯、线确定单元，用于从得到的所述Ξ维冠脉图像中，确定出冠脉CTA图像中 屯、线；
[0037] 展示单元，用于将确定出所述冠脉CTA图像中屯、线的所述Ξ维冠脉图像展示给所 述医生。
[0038] 结合第二方面，本发明实施例提供了第二方面的第二种可能的实施方式，其中，所 述图像中屯、线确定单元，包括：
[0039] 计算子单元，用于利用霍夫圆检测算法，从所述冠脉CTA图像中确定出升主动脉图 像的中屯、；
[0040] 提取子单元，用于根据确定的所述升主动脉图像的中屯、，从所述冠脉CTA图像中提 取出冠脉图像；
[0041] 图像中屯、线确定子单元，用于从提取出的所述冠脉图像中，确定出冠脉CTA图像中 屯、线。
[0042] 结合第二方面，本发明实施例提供了第二方面的第Ξ种可能的实施方式，其中，所 述计算模块，包括：
[0043] 生成单元，用于随机生成预设数量的造影角度组合，其中，每个所述造影角度组合 均包括一造影主角度和一造影副角度；
[0044] 计算单元，用于根据预设的迭代次数，通过遗传算法对随机生成的所述造影角度 组合进行迭代计算，得到多种造影角度组合下所述感兴趣血管段的模拟冠脉造影图像面积 值；
[0045] 确定单元，用于确定所述感兴趣血管段的最大模拟冠脉造影图像面积值对应的造 影角度组合作为所述感兴趣血管段的最佳造影角度组合。
[0046] 结合第二方面，本发明实施例提供了第二方面的第四种可能的实施方式，其中，所 述计算单元，包括：
[0047] 模拟造影单元，用于根据所述造影角度组合中记录的造影主角度和造影副角度， 对所述Ξ维冠脉图像进行模拟造影，分别得到所述造影角度组合下所述Ξ维冠脉图像中所 述感兴趣血管段和其他血管段的模拟造影图像区域；
[0048] 参考值计算单元，用于根据公式
[0049]
确定所述感兴趣血管段的模拟 造影图像区域中各像素点的参考值f(i，j)，其中，select(i，j) = l表示像素点在所述感兴 趣血管段上，unselect (i，j ) = 1表示像素点在其他血管段上，i表示所述感兴趣血管段的模 拟造影图像区域中各像素点的横坐标，j表示所述感兴趣血管段的模拟造影图像区域中各 像素点的纵坐标；
[0050] 冠脉造影图像面积值计算单元，用于根据确定出的所述感兴趣血管段的模拟造影 图像区域中各像素点的参考值，得到所述造影角度组合下所述感兴趣血管段的冠脉造影图 像面积值。
[0051] 本发明实施例提供的一种冠脉图像造影方法和装置，通过计算确定医生选择的感 兴趣血管段的最佳造影主角度和最佳造影副角度，然后根据确定的最佳造影主角度和最佳 造影副角度对病人的冠脉进行造影，得到病人的冠脉造影图像，与现有技术中需对病人的 冠脉进行多次造影后才能明确病人冠脉病变位置相比，通过最佳造影角度组合对病人的冠 脉只需进行一次造影，就可W确定病人的冠脉的病变位置，提高了冠脉造影的准确性，而 且，由于减少了冠脉造影次数，所W减少了病人得到的射线福射量和对比剂用量，尽可能降 低了冠脉造影给病人身体所带来的损害。
[0052] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合 所附附图，作详细说明如下。
【附图说明】
[0053] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附 图作简单地介绍，应当理解，W下附图仅示出了本发明的某些实施例，因此不应被看作是对 范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W根据运 些附图获得其他相关的附图。
[0054] 图1示出了本发明实施例1所提供的一种冠脉图像造影方法的流程图；
[0055] 图2示出了本发明实施例2所提供的一种冠脉图像造影装置的结构示意图。
【具体实施方式】
[0056] 下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整 地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在 此处附图中描述和示出的本发明实施例的组件可各种不同的配置来布置和设计。因 此，W下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的 范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做 出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0057] 相关技术中，在对病人的冠脉进行冠冕造影的过程中，由于病人的屯、脏冠脉形状 不同，且冠脉的走线复杂，在利用默认体位或者医生经验选择的造影体位对病人冠脉进行 造影的过程中，仅通过一次造影往往无法充分暴露病人冠脉的病变位置，只能对病人的冠 脉进行多次造影后才能明确病人冠脉的病变位置，造成冠脉造影的准确性较差，而且对病 人进行多次冠脉造影，会对病人增加过多的射线福射量和对比剂用量，给病人的身体带来 不必要的损害。
[0化引实施例1
[0059]本实施例提供了一种冠脉图像造影方法，本发明实施例的执行主体为包括服务器 和CTA机的医学影像系统，其中，在医学影像系统中，服务器与CTA机进行数据交互，医生先 通过服务器向CTA机发出指令，得到病人的冠脉CTA图像，然后服务器会对获取的病人的冠 脉CTA图像进行处理，确定医生选择的感兴趣血管段的最佳造影主角度和最佳造影副角度， 然后根据确定的最佳造影主角度和最佳造影副角度对病人的冠脉进行造影，得到病人的冠 脉造影图像。
[0060] 其中，服务器包括显示设备和鼠标、键盘等输入设备，w和医生进行数据交互；而 且，服务器可W使用现有的任何型号的计算终端，运里不再一一寶述。
[0061] 参见图1，本实施例提出的一种冠脉图像造影方法，包括W下步骤：
[0062] 步骤100、获取病人的冠脉CTA图像。
[0063] 在步骤100中，服务器获取病人的冠脉CTA图像包括W下步骤(1)至步骤(3):
[0064] (1)监测医生触发的病人冠脉CTA图像获取指令；
[0065] (2)向CTA机发出CTA图像获取指令，W使得CTA机获取病人的冠脉CTA图像；
[0066] (3)获取CTA机得到的病人冠脉CTA图像。
[0067] 步骤102、对得到的冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图像。
[0068] 医生可W通过Ξ维重建图像清楚的看到每条冠脉的走势与状态，从而帮助医生对 病人整体的病情有更清楚的了解，医生也能更加清楚的对病变冠脉的狭窄程度与狭窄位置 进行判断，从而选择感兴趣血管段。
[0069] 具体地，步骤102包括W下步骤(1)至步骤(3):
[0070] (1)对得到的冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图像；
[0071] (2)从得到的Ξ维冠脉图像中，确定出冠脉CTA图像中屯、线；
[0072] (3)将确定出冠脉CTA图像中屯、线的Ξ维冠脉图像展示给医生。
[0073] 其中，服务器对得到的冠脉CTA图像通过现有的任何Ξ维重建算法，对CT冠脉造影 图像进行Ξ维重建，运里不再一一寶述。
[0074] 上述步骤2具体包括W下步骤(21)至步骤(23):
[0075] (21)利用霍夫圆检测算法，从冠脉CTA图像中确定出升主动脉图像的中屯、；
[0076] (22)根据确定的升主动脉图像的中屯、，从冠脉CTA图像中提取出冠脉图像；
[0077] (23)从提取出的冠脉图像中，确定出冠脉CTA图像中屯、线。
[0078] 在步骤21中，利用升主动脉圆柱状特性，利用霍夫圆检测算法从病人序列图像初 始层逐层对整个图像进行圆形检测，当其中某一个圆屯、在连续几层上位置变化范围很小的 时候，服务器可W确定该圆屯、即为升主动脉的中屯、。
[0079] 步骤22具体包括W下步骤(221)至步骤(223):
[0080] (221)根据升主动脉图像的中屯、，利用区域增长算法提取升主动脉图像；
[0081] (222)利用膨胀腐蚀算法对Ξ维冠脉图像中的冠脉图像进行增强；
[0082] (223)利用升主动脉图像，对增强后的冠脉图像进行区域增长，提取出冠脉图像。
[0083] 在步骤(223)中，在增强后的图像中利用已经提取出的升主动脉最外表面上的所 有点作为种子点，在增强图像中做区域增长，其中增长长度最长且不同的两条血管即为冠 脉左右两只。升主动脉最外表面的提取方法为，利用形态学算法对升主动脉进行腐蚀操作， 然后利用原升主动脉图像减去腐蚀后的升主动脉图像即可得出。
[0084] 针对血管闭塞严重的血管，因为闭塞部分血管显影不明显，可能无法自动提取出， 运里需要医生针对闭塞严重血管进行手动提取。手动提取过程可W利用血管CT值局部相似 性，W多个种子点在小范围内对血管进行提取从而保证血管的连续性。
[0085] 步骤23具体包括W下步骤(231)至步骤(234):
[0086] (231)对冠脉图像进行二值化操作，得到冠脉的二值图像；
[0087] (232)对冠脉的二值图像进行Ξ维细化操作，提取血管的粗略中屯、线；
[0088] (233)对冠脉的粗略中屯、线进行连通域检测，确定出冠脉的端点；
[0089] (234)计算升主动脉图像的中屯、到冠脉的端点的最短路径，W得到的最短路径作 为冠脉CTA图像中屯、线。
[0090] 在步骤233中，连通域检测结果表示每一点在1个像素的Ξ维邻域内相邻点的个 数。当结果为1时，该像素点即为冠脉的端点。统计连通域检测所有结果为1且属于冠脉上的 点即为每条冠脉分支的端点。
[0091] 在步骤234中，根据欧式距离变换对提取出的冠脉的二值图像进行距离变换，在距 离场中分别计算升主动脉几何中点到冠脉末梢点的最短路径。每个末梢点与升主动脉的几 何中屯、所构成的最短路径记为该条血管所对应的中屯、线。
[0092] 在步骤3中，服务器通过预先安装的显示设备向医生展示确定出冠脉CTA图像中屯、 线的Ξ维冠脉图像。
[0093] 步骤104、获取医生在Ξ维冠脉图像中选择的感兴趣血管段。
[0094] 在步骤104中，医生会通过服务器的输入设备在服务器的显示设备上显示的Ξ维 冠脉图像中，在确定出的血管中屯、线上设置感兴趣血管段的起始点和终止点，来选择Ξ维 冠脉图像中的感兴趣血管段。
[00M]步骤106、通过计算确定感兴趣血管段的最佳造影角度组合，最佳造影角度组合包 括最佳造影主角度和最佳造影副角度。
[0096] 其中，投影主角度是投影光源在冠脉CTA图像中屯、线方向上旋转的角度，投影副角 度是投影光源在垂直于冠脉CTA图像中屯、线方向上旋转的角度。
[0097] 投影光源就是设置在CTA机上，用于对冠脉进行造影得到冠脉造影图像的射线源。
[0098] 上述感兴趣血管段的最佳造影角度组合应符合W下条件：
[0099] (1)最佳造影角度组合应最大可能地使CTA机在对冠脉造影时垂直于感兴趣血管 段，即应满足最小投影距离；
[0100] (2)在最佳造影角度组合下，上述感兴趣血管段应尽可能没有被其他血管遮盖，即 满足感兴趣血管段投影面积的最小遮盖。
[0101] 步骤108、根据最佳造影主角度和最佳造影副角度对Ξ维冠脉图像进行造影，得到 最佳造影角度组合下的冠脉造影图像，并将得到的冠脉造影图像进行展示。
[0102] 具体地，步骤108包括W下步骤(1)至步骤(3):
[0103] (1)根据得到的最佳造影主角度和最佳造影副角度，调整CTA机的射线源的造影角 度，使得CTA机在最佳造影角度组合下对病人的冠脉进行造影，得到最佳造影角度组合下的 冠脉造影图像；
[0104] (2)获取CTA机在最佳造影角度组合下投影所得到的冠脉造影图像；
[0105] (3)将得到的冠脉造影图像进行展示。
[0106] 综上所述，本实施例提供的一种冠脉图像造影方法，通过计算确定医生选择的感 兴趣血管段的最佳造影主角度和最佳造影副角度，然后根据确定的最佳造影主角度和最佳 造影副角度对病人的冠脉进行造影，得到病人的冠脉造影图像，与现有技术中需对病人的 冠脉进行多次造影后才能明确病人冠脉病变位置相比，通过最佳造影角度组合对病人的冠 脉只需进行一次造影，就可W确定病人的冠脉的病变位置，提高了冠脉造影的准确性，而 且，由于减少了冠脉造影次数，所W减少了病人得到的射线福射量和对比剂用量，尽可能降 低了冠脉造影给病人身体所带来的损害。
[0107] 为了得到感兴趣血管段的最佳造影角度组合，确定感兴趣血管段的最佳造影角度 组合，包括W下步骤(1)至步骤(3):
[0108] (1)随机生成预设数量的造影角度组合，其中，每个造影角度组合均包括一造影主 角度和一造影副角度；
[0109] (2)根据预设的迭代次数，通过遗传算法对随机生成的造影角度组合进行迭代计 算，得到多种造影角度组合下感兴趣血管段的模拟冠脉造影图像面积值；
[0110] (3)确定感兴趣血管段的最大模拟冠脉造影图像面积值对应的造影角度组合作为 感兴趣血管段的最佳造影角度组合。
[0111] 在上述步骤1中，在造影主角度和造影副角度的预设角度范围内，随机生成预设数 量的造影角度组合。
[0112] 在上述步骤2中，计算任一造影角度组合下冠脉造影图像面积值的过程包括如下 步骤(21)至步骤(23):
[0113] (21)根据造影角度组合中记录的造影主角度和造影副角度，对Ξ维冠脉图像进行 模拟造影，分别得到造影角度组合下Ξ维冠脉图像中感兴趣血管段和其他血管段的模拟造 影图像区域；
[0114] (22)根据公式
确定感兴趣血管 段的模拟造影图像区域中各像素点的参考值f(i，j)，其中，select(i，j) = l表示像素点在 感兴趣血管段上，unselect(i，j ) = 1表示像素点在其他血管段上，i表示感兴趣血管段的模 拟造影图像区域中各像素点的横坐标，j表示感兴趣血管段的模拟造影图像区域中各像素 点的纵坐标；
[0115] (23)根据确定出的感兴趣血管段的模拟造影图像区域中各像素点的参考值，得到 造影角度组合下感兴趣血管段的冠脉造影图像面积值。
[0116] 在步骤21中，对感兴趣血管段中的模拟造影图像区域中的像素点设置预定的感兴 趣血管段标识，对其他血管段的模拟造影图像区域中的像素点设置预定的其他血管段标 识。如果感兴趣血管段中的模拟造影图像区域中的某一像素点被其他血管段的模拟造影图 像区域覆盖时，那么该像素点的就即拥有感兴趣血管段标识，又拥有其他血管段标识。
[0117] 在步骤22中，根据感兴趣血管段中的模拟造影图像区域中各像素点中所具有的标 识，确定感兴趣血管段的模拟造影图像区域中各像素点的参考值f(i，j)。
[0118] 如果感兴趣血管段中的模拟造影图像区域中的某一像素点即拥有感兴趣血管段 标识，又拥有其他血管段标识，那么该像素点的f(i，j) = 〇;若只拥有感兴趣血管段标识，那 么该像素点的f (i，j) = 1;若只拥有其他血管段标识，那么该像素点的f (i，j) = 0。
[0119] 在步骤23中，任一影角度组合下感兴趣血管段的冠脉造影图像面积值F=Ef(i， j )，如果F值越大，说明该F值所对应的造影角度组合就越符合上述感兴趣血管段的最佳造 影角度组合所满足的条件。
[0120] 在一个具体的实施方式中，在上述步骤2中，遗传算法的基本运算过程包括如下步 骤(a)至步骤(f):
[0121] (a)初始化:设置进化代数计数器t = 0，设置最大进化代数Τ，随机生成Μ个体作为 初始群体P(〇);
[0122] (b)个体评价:计算群体P(t)中各个个体的适应度。
[0123] (C)选择运算:将选择算子作用于群体。选择的目的是把优化的个体直接遗传到下 一代或通过配对交叉产生新的个体再遗传到下一代。选择操作是建立在群体中个体的适应 度评估基础上的。
[0124] (d)交叉运算:将交叉算子作用于群体。遗传算法中起核屯、作用的就是交叉算子。
[0125] (e)变异运算:将变异算子作用于群体。即是对群体中的个体串的某些基因座上的 基因值作变动。
[01%]群体P(t)经过选择、交叉、变异运算之后得到下一代群体P(t+1)。
[0127] (f)终止条件判断:若t = T，则W进化过程中所得到的具有最大适应度个体作为最 优解输出，终止计算。
[0128] 优选地，在上述步骤中，遗传算法使用的具体参数如下：
[0129] 设置每个投影的主角度与副角度的组合为一个个体，并在血管机机械约束范围内 随机设置30组不同的个体;个体的适应度为视角的代价函数优化迭代次数设置为200;交叉 概率设置为0.7;变异概率是0.05。
[0130] 综上所述，通过迭代算法，结合血管间遮盖的影响及造影系统机构的实际限制，确 定感兴趣血管段的最佳造影角度组合，从而在根据最佳造影角度组合对感兴趣血管段进行 造影时，提高了后续冠脉介入诊疗的便利性;且采用该最佳投射角度进行冠脉造影时，使得 造射角度选择的时间也相应缩短，节省了成本，提高了效率。
[0131] 实施例2
[0132] 参见图2,本实施例提供一种冠脉图像造影装置，用于执行上述的冠脉图像造影方 法，包括：
[0133] 图像获取模块200,用于获取病人的冠脉CTA图像;Ξ维重建模块202,用于对得到 的冠脉CTA图像进行Ξ维重建，得到Ξ维冠脉图像;血管获取模块204,用于获取医生在Ξ维 冠脉图像中选择的感兴趣血管段;计算模块206，用于通过计算确定感兴趣血管段的最佳造 影角度组合，最佳造影角度组合包括最佳造影主角度和最佳造影副角度;造影模块208，用 于根据最佳造影主角度和最佳造影副角度对Ξ维冠脉图像进行造影，得到最佳造影角度组 合下的冠脉造影图像，并将得到的冠脉造影图像进行展示。
[0134] 具体地，Ξ维重建模块202,包括:Ξ维重建单元，用于对得到的冠脉CTA图像进行 Ξ维重建，得到Ξ维冠脉图像；图像中屯、线确定单元，用于从得到的Ξ维冠脉图像中，确定 出冠脉CTA图像中屯、线；展示单元，用于将确定出冠脉CTA图像中屯、线的Ξ维冠脉图像展示 给医生。
[0135] 进一步地，图像中屯、线确定单元，包括:计算子单元，用于利用霍夫圆检测算法，从 冠脉CTA图像中确定出升主动脉图像的中屯、;提取子单元，用于根据确定的升主动脉图像的 中屯、，从冠脉CTA图像中提取出冠脉图像；图像中屯、线确定子单元，用于从提取出的冠脉图 像中，确定出冠脉CTA图像中屯、线。
[0136] 为了得到感兴趣血管段的最佳造影角度组合，计算模块206,包括:生成单元，用于 随机生成预设数量的造影角度组合，其中，每个造影角度组合均包括一造影主角度和一造 影副角度;计算单元，用于根据预设的迭代次数，通过遗传算法对随机生成的造影角度组合 进行迭代计算，得到多种造影角度组合下感兴趣血管段的模拟冠脉造影图像面积值;确定 单元，用于确定感兴趣血管段的最大模拟冠脉造影图像面积值对应的造影角度组合作为感 兴趣血管段的最佳造影角度组合。
[0137] 其中，计算单元，包括：
[0138] 模拟造影单元，用于根据造影角度组合中记录的造影主角度和造影副角度，对Ξ 维冠脉图像进行模拟造影，分别得到造影角度组合下Ξ维冠脉图像中感兴趣血管段和其他 血管段的模拟造影图像区域；
[0139] 参考值计算单元，用于根据公式
[0140]
确定感兴趣血管段的模拟造影 图像区域中各像素点的参考值f(i，j)，其中，select(i，j) = l表示像素点在感兴趣血管段 上，unselect(i，j ) = 1表示像素点在其他血管段上，i表示感兴趣血管段的模拟造影图像区 域中各像素点的横坐标，j表示感兴趣血管段的模拟造影图像区域中各像素点的纵坐标；
[0141] 冠脉造影图像面积值计算单元，用于根据确定出的感兴趣血管段的模拟造影图像 区域中各像素点的参考值，得到造影角度组合下所述感兴趣血管段的冠脉造影图像面积 值。
[0142] 综上所述，通过迭代算法，结合血管间遮盖的影响及造影系统机构的实际限制，确 定感兴趣血管段的最佳造影角度组合，从而在根据最佳造影角度组合对感兴趣血管段进行 造影时，提高了后续冠脉介入诊疗的便利性;且采用该最佳投射角度进行冠脉造影时，使得 造射角度选择的时间也相应缩短，节省了成本，提高了效率。
[0143] 综上所述，本实施例提供的一种冠脉图像造影装置，通过计算确定医生选择的感 兴趣血管段的最佳造影主角度和最佳造影副角度，然后根据确定的最佳造影主角度和最佳 造影副角度对病人的冠脉进行造影，得到病人的冠脉造影图像，与现有技术中需对病人的 冠脉进行多次造影后才能明确病人冠脉病变位置相比，通过最佳造影角度组合对病人的冠 脉只需进行一次造影，就可W确定病人的冠脉的病变位置，提高了冠脉造影的准确性，而 且，由于减少了冠脉造影次数，所W减少了病人得到的射线福射量和对比剂用量，尽可能降 低了冠脉造影给病人身体所带来的损害。
[0144] 本发明实施例所提供的进行冠脉图像造影方法的计算机程序产品，包括存储了程 序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所 述的方法，具体实现可参见方法实施例，在此不再寶述。
[0145] 所属领域的技术人员可W清楚地了解到，为描述的方便和简洁，上述描述的系统、 装置和单元的具体工作过程，可W参考前述方法实施例中的对应过程，在此不再寶述。
[0146] 在本申请所提供的几个实施例中，应该理解到，所掲露的系统、装置和方法，可W 通过其它的方式实现。W上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分， 仅仅为一种逻辑功能划分，实际实现时可W有另外的划分方式，又例如，多个单元或组件可 W结合或者可W集成到另一个系统，或一些特征可W忽略，或不执行。另一点，所显示或讨 论的相互之间的禪合或直接禪合或通信连接可w是通过一些通信接口，装置或单元的间接 禪合或通信连接，可W是电性，机械或其它的形式。
[0147] 所述作为分离部件说明的单元可W是或者也可W不是物理上分开的，作为单元显 示的部件可W是或者也可W不是物理单元，即可W位于一个地方，或者也可W分布到多个 网络单元上。可W根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目 的。
[0148] 另外，在本发明各个实施例中的各功能单元可W集成在一个处理单元中，也可W 是各个单元单独物理存在，也可W两个或两个W上单元集成在一个单元中。
[0149] 所述功能如果W软件功能单元的形式实现并作为独立的产品销售或使用时，可W 存储在一个计算机可读取存储介质中。基于运样的理解，本发明的技术方案本质上或者说 对现有技术做出贡献的部分或者该技术方案的部分可软件产品的形式体现出来，该计 算机软件产品存储在一个存储介质中，包括若干指令用W使得一台计算机设备(可W是个 人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。 而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memo巧）、随机存取存 储器(RAM,Random Access Memory)、磁碟或者光盘等各种可W存储程序代码的介质。
[0150] W上所述，仅为本发明的【具体实施方式】，但本发明的保护范围并不局限于此，任何 熟悉本技术领域的技术人员在本发明掲露的技术范围内，可轻易想到变化或替换，都应涵 盖在本发明的保护范围之内。因此，本发明的保护范围应所述W权利要求的保护范围为准。
【主权项】
1. 一种冠脉图像造影方法，其特征在于，包括： 获取病人的冠脉计算机断层摄影血管造影术CTA图像； 对得到的所述冠脉CTA图像进行三维重建，得到三维冠脉图像； 获取医生在所述三维冠脉图像中选择的感兴趣血管段； 通过计算确定所述感兴趣血管段的最佳造影角度组合，所述最佳造影角度组合包括最 佳造影主角度和最佳造影副角度； 根据所述最佳造影主角度和所述最佳造影副角度对所述三维冠脉图像进行造影，得到 所述最佳造影角度组合下的冠脉造影图像，并将得到的所述冠脉造影图像进行展示。2. 根据权利要求1所述的方法，其特征在于，所述对得到的所述冠脉CTA图像进行三维 重建，得到三维冠脉图像，包括： 对得到的所述冠脉CTA图像进行三维重建，得到三维冠脉图像； 从得到的所述三维冠脉图像中，确定出冠脉CTA图像中心线； 将确定出所述冠脉CTA图像中心线的所述三维冠脉图像展示给所述医生。3. 根据权利要求2所述的方法，其特征在于，所述从得到的所述三维冠脉图像中，确定 出冠脉CTA图像中心线，包括： 利用霍夫圆检测算法，从所述冠脉CTA图像中确定出升主动脉图像的中心； 根据确定的所述升主动脉图像的中心，从所述冠脉CTA图像中提取出冠脉图像； 从提取出的所述冠脉图像中，确定出冠脉CTA图像中心线。4. 根据权利要求1所述的方法，其特征在于，确定所述感兴趣血管段的最佳造影角度组 合，所述最佳造影角度组合包括最佳造影主角度和最佳造影副角度，包括： 随机生成预设数量的造影角度组合，其中，每个所述造影角度组合均包括一造影主角 度和一造影副角度； 根据预设的迭代次数，通过遗传算法对随机生成的所述造影角度组合进行迭代计算， 得到多种造影角度组合下所述感兴趣血管段的模拟冠脉造影图像面积值； 确定所述感兴趣血管段的最大模拟冠脉造影图像面积值对应的造影角度组合作为所 述感兴趣血管段的最佳造影角度组合。5. 根据权利要求4所述的方法，其特征在于，根据预设的迭代次数，通过遗传算法对随 机生成的所述造影角度组合进行迭代计算，得到多种造影角度组合下所述感兴趣血管段的 冠脉造影图像面积值，包括： 根据所述造影角度组合中记录的造影主角度和造影副角度，对所述三维冠脉图像进行 模拟造影，分别得到所述造影角度组合下所述三维冠脉图像中所述感兴趣血管段和其他血 管段的模拟诰影图像冈域： 根据公确定所述感兴趣血管段的模 拟造影图像区域中各像素点的参考值f(i，j)，其中，select(i，j) = l表示像素点在所述感 兴趣血管段上，unselect( i，j) = 1表示像素点在其他血管段上，i表示所述感兴趣血管段的 模拟造影图像区域中各像素点的横坐标，j表示所述感兴趣血管段的模拟造影图像区域中 各像素点的纵坐标； 根据确定出的所述感兴趣血管段的模拟造影图像区域中各像素点的参考值，得到所述 造影角度组合下所述感兴趣血管段的冠脉造影图像面积值。6. -种冠脉图像造影装置，其特征在于，包括： 图像获取模块，用于获取病人的冠脉计算机断层摄影血管造影术CTA图像； 三维重建模块，用于对得到的所述冠脉CTA图像进行三维重建，得到三维冠脉图像； 血管获取模块，用于获取医生在所述三维冠脉图像中选择的感兴趣血管段； 计算模块，用于通过计算确定所述感兴趣血管段的最佳造影角度组合，所述最佳造影 角度组合包括最佳造影主角度和最佳造影副角度； 造影模块，用于根据所述最佳造影主角度和所述最佳造影副角度对所述三维冠脉图像 进行造影，得到所述最佳造影角度组合下的冠脉造影图像，并将得到的所述冠脉造影图像 进行展示。7. 根据权利要求6所述的装置，其特征在于，所述三维重建模块，包括： 三维重建单元，用于对得到的所述冠脉CTA图像进行三维重建，得到三维冠脉图像； 图像中心线确定单元，用于从得到的所述三维冠脉图像中，确定出冠脉CTA图像中心 线； 展示单元，用于将确定出所述冠脉CTA图像中心线的所述三维冠脉图像展示给所述医 生。8. 根据权利要求7所述的装置，其特征在于，所述图像中心线确定单元，包括： 计算子单元，用于利用霍夫圆检测算法，从所述冠脉CTA图像中确定出升主动脉图像的 中心； 提取子单元，用于根据确定的所述升主动脉图像的中心，从所述冠脉CTA图像中提取出 冠脉图像； 图像中心线确定子单元，用于从提取出的所述冠脉图像中，确定出冠脉CTA图像中心 线。9. 根据权利要求6所述的装置，其特征在于，所述计算模块，包括： 生成单元，用于随机生成预设数量的造影角度组合，其中，每个所述造影角度组合均包 括一造影主角度和一造影副角度； 计算单元，用于根据预设的迭代次数，通过遗传算法对随机生成的所述造影角度组合 进行迭代计算，得到多种造影角度组合下所述感兴趣血管段的模拟冠脉造影图像面积值； 确定单元，用于确定所述感兴趣血管段的最大模拟冠脉造影图像面积值对应的造影角 度组合作为所述感兴趣血管段的最佳造影角度组合。10. 根据权利要求9所述的装置，其特征在于，所述计算单元，包括： 模拟造影单元，用于根据所述造影角度组合中记录的造影主角度和造影副角度，对所 述三维冠脉图像进行模拟造影，分别得到所述造影角度组合下所述三维冠脉图像中所述感 兴趣血管段和其他血管段的模拟造影图像区域； 参考值计算单元，用于根据公式确定所述感兴趣血管段的模拟造影图 像区域中各像素点的参考值f(i，j)，其中，seleCt(i，j) = l表示像素点在所述感兴趣血管 段上，unselect( i，j) = 1表示像素点在其他血管段上，i表示所述感兴趣血管段的模拟造影 图像区域中各像素点的横坐标，j表示所述感兴趣血管段的模拟造影图像区域中各像素点 的纵坐标； 冠脉造影图像面积值计算单元，用于根据确定出的所述感兴趣血管段的模拟造影图像 区域中各像素点的参考值，得到所述造影角度组合下所述感兴趣血管段的冠脉造影图像面 积值。
【文档编号】A61B6/00GK105877767SQ201610201535
【公开日】2016年8月24日
【申请日】2016年3月31日
【发明人】陈韵岱, 杨俊杰, 郭军, 田峰, 王民汉, 周宇, 李晓光
【申请人】北京思创贯宇科技开发有限公司