CN114418947B - 心脏冠脉血管的命名方法及装置 - Google Patents
心脏冠脉血管的命名方法及装置 Download PDFInfo
- Publication number
- CN114418947B CN114418947B CN202111545244.4A CN202111545244A CN114418947B CN 114418947 B CN114418947 B CN 114418947B CN 202111545244 A CN202111545244 A CN 202111545244A CN 114418947 B CN114418947 B CN 114418947B
- Authority
- CN
- China
- Prior art keywords
- point
- line
- intersection
- ventricular
- atrial
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 210000004204 blood vessel Anatomy 0.000 title abstract description 26
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 107
- 210000003484 anatomy Anatomy 0.000 claims abstract description 89
- 230000002861 ventricular Effects 0.000 claims abstract description 72
- 230000001746 atrial effect Effects 0.000 claims abstract description 70
- 230000000747 cardiac effect Effects 0.000 claims abstract description 66
- 230000011218 segmentation Effects 0.000 claims abstract description 54
- 210000002837 heart atrium Anatomy 0.000 claims abstract description 43
- 210000001308 heart ventricle Anatomy 0.000 claims abstract description 17
- 239000013598 vector Substances 0.000 claims description 102
- 210000005240 left ventricle Anatomy 0.000 claims description 22
- 210000005241 right ventricle Anatomy 0.000 claims description 18
- 210000005246 left atrium Anatomy 0.000 claims description 12
- 210000005245 right atrium Anatomy 0.000 claims description 10
- PPTYJKAXVCCBDU-UHFFFAOYSA-N Rohypnol Chemical compound N=1CC(=O)N(C)C2=CC=C([N+]([O-])=O)C=C2C=1C1=CC=CC=C1F PPTYJKAXVCCBDU-UHFFFAOYSA-N 0.000 claims description 6
- 210000000596 ventricular septum Anatomy 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 16
- 238000004590 computer program Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010968 computed tomography angiography Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 238000002583 angiography Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请提供了一种心脏冠脉血管的命名方法、装置、计算机存储介质及处理器。其中,该方法包括:获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;获取心脏心房心室分割模型;根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本申请解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
Description
技术领域
本申请涉及图像处理领域,具体而言,涉及一种心脏冠脉血管的命名方法、装置、计算机存储介质及处理器。
背景技术
在对心脏冠脉血管进行命名时,研究学者一般使用AI命名的方法。这种方法的局限性在于,没有参考实际的解剖结构,同时需要大量数据集进行大规模耗时的训练,而且AI无法解释结果的准确性。又有研究学者使用三维模板匹配的命名方法,这种方法的局限性在于,需要提供合适的模板,模板的选择对最终结果的影响非常大,同时,匹配的精度低,无法覆盖大规模的数据,无法应对血管的变异情况。又有研究学者使用二维树形结构血管命名的方法,该方法依赖三维血管树抽象化为二维建模,同样无法应对血管变异情况。因此,有待提出一种技术方案能准确快速地对心脏冠脉血管进行命名的方法。
发明内容
本申请的主要目的在于提供一种心脏冠脉血管的命名方法、装置、计算机存储介质及处理器,以解决现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
为了实现上述目的,根据本申请的一个方面,提供了一种心脏冠脉血管的命名方法,该方法包括:获取心脏冠脉血管的3D影像,且从所述3D影像中提取出所述心脏冠脉血管的中心线树;获取心脏心房心室分割模型;根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,所述解剖学结构至少包括房室沟交线和房室沟平面;根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名。
进一步地,所述心脏心房心室分割模型包括右心房右心室模型和左心房左心室模型,根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,所述解剖学结构至少包括房室沟交线和房室沟平面,该方法包括:将所述右心房右心室模型扩大预定倍数,使得右心房和右心室产生第一交集;根据所述第一交集确定右房室沟交线;根据所述右房室沟交线确定右房室沟平面;将所述左心房左心室模型扩大所述预定倍数,使得左心房和左心室产生第二交集;根据所述第二交集确定左房室沟交线;根据所述左房室沟交线确定左房室沟平面。
进一步地,根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,该方法还包括:获取所述右房室沟交线的第一中心点和所述左房室沟交线的第二中心点;获取所述第一中心点和所述第二中心点的第一连线的第一法向量;获取所述右房室沟平面的第三中心点和所述左房室沟平面的第四中心点;获取所述第三中心点和所述第四中心点的第二连线的中点;根据所述第一法向量和所述中点构造十字交叉面。
进一步地,获取所述右房室沟交线的第一中心点,该方法包括:获取所述右房室沟交线上的任意一个第一目标点,到除所述第一目标点以外的所有点的距离的第一平均值;确定最小的所述第一平均值对应的所述第一目标点为所述第一中心点。
进一步地,获取所述右房室沟平面的第三中心点,该方法包括:获取所述右房室沟平面上的任意一个第二目标点,到除所述第二目标点以外的所有点的距离的第二平均值;确定最小的所述第二平均值对应的所述第二目标点为所述第三中心点。
进一步地,根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,该方法还包括:获取左右心室室间沟交线。
进一步地,所述中心线树包括左冠和右冠,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法包括:将所述右冠中穿过所述十字交叉面,并向后延伸最远的中心线命名为主PL;获取所述中心线树中除所述主PL以外的中心线,与所述主PL的分叉点;获取预定中心线的分叉点到末点的线段上的所有点,到所述左右心室室间沟交线的最小距离的平均值;获取所述预定中心线的分叉点到末点的线段上的所有点,到所述左右心室室间沟交线的距离的最小值;所述平均值和所述最小值的乘积最小,对应的所述预定中心线为PD。
进一步地,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法还包括:将所述右冠的起点至所述主PL与所述PD的分叉点部分,命名为RCA。
进一步地,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法还包括:将所述RCA分为第一前段、中段和第一后段;将从所述第一前段发出,并且所述RCA的分叉点和末点构成的向量的第二法向量与所述右房室沟平面的第三法向量的点积小于预定数值的部分,命名为CA;将从所述中段发出,并且所述RCA的分叉点和末点构成的向量的第四法向量与所述右房室沟平面的第五法向量的点积小于所述预定数值的部分,命名为Rv;将从所述第一后段发出,并且所述RCA的分叉点和末点构成的向量的第六法向量与所述右房室沟平面的第七法向量的点积小于所述预定数值的部分,命名为AM。
进一步地,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法还包括:将所述RCA分为第二前段和第二后段;将从所述第二前段发出,并且所述RCA的分叉点和末点构成的向量的第八法向量与所述右房室沟平面的第九法向量的点积小于预定数值的部分,命名为SAN;将从所述第二后段发出,并且所述RCA的分叉点和末点构成的向量的第十法向量与所述右房室沟平面的第十一法向量的点积小于预定数值的部分,命名为RA Branch。
进一步地,所述中心线树包括左冠和右冠,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法包括:将所述左冠中与所述左右心室室间沟交线评分最小的中心线命名为LAD,其中,所述评分是根据中心线上的所有点到所述左右心室室间沟交线的最小距离的平均值,与中心线上的所有点到所述左右心室室间沟交线的距离的最小值之和,再除以中心线的长度得到的;获取所述左冠中的中心线中除所述LAD以外的中心线,与所述LAD的分叉点;确定所述分叉点点序最小的点为LM分叉点,其中,所述分叉点点序指的是所有除所述LAD以外的中心线与所述LAD中心线的分叉点排序得到分叉点点序;根据所述LM分叉点将所述左冠分为LAD冠脉树和LCX冠脉树。
进一步地,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法还包括:将所述LAD上从所述LM分叉点起至所述LAD末点的所有点,拟合为LAD平面;计算所述LAD冠脉树中的中心线与所述LAD的分叉点与中心线末点构成的向量,与所述LAD平面的法向量的点积;将所述点积小于第一预定值的中心线命名为D,将所述点积大于或者等于第一预定值的中心线命名为S;所述中心线若为D,且所述中心线与所述LAD的分叉点与所述LM分叉点之间的距离小于预定距离的情况下,将所述中心线命名为RM。
进一步地,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,该方法还包括:将所述LCX冠脉树中与所述左房室沟交线评分最小的中心线命名为主LCX,评分规则为:所述LM分叉点后中心线上的所有点到所述左房室沟交线的最小距离的平均值,与所述LM分叉点后中心线上的所有点到所述左房室沟交线的距离的最小值之和,再除以中心线的长度得到所述评分;计算所述LCX冠脉树中除所述主LCX以外的中心线,与所述主LCX的分叉点至中心线末点构成的向量,与所述左房室沟平面的法向量的点积;将所述点积小于第二预定值的中心线命名为OM,将所述点积大于第三预定值的中心线命名为LA Branch,否则,命名为LCX,其中,所述第三预定值大于所述第二预定值;在所述中心线为OM的情况下,且所述中心线与所述主LCX的分叉点与所述LM分叉点之间的距离小于预定距离的情况下,将所述中心线命名为RM。
根据本申请的另一方面,提供了一种心脏冠脉血管的命名装置,该装置包括:第一获取单元,用于获取心脏冠脉血管的3D影像,且从所述3D影像中提取出所述心脏冠脉血管的中心线树;第二获取单元,用于获取心脏心房心室分割模型;确定单元,用于根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,所述解剖学结构至少包括房室沟交线和房室沟平面;命名单元,用于根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名。
根据本申请的另一方面,还提供了一种计算机可读存储介质,计算机可读存储介质包括存储的程序,其中,在程序运行时控制计算机可读存储介质所在设备执行上述的心脏冠脉血管的命名方法。
根据本申请的另一方面,还提供了一种处理器,处理器用于运行程序,其中,程序运行时执行上述的心脏冠脉血管的命名方法。
应用本申请的技术方案,首先获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;并且获取心脏心房心室分割模型;之后根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;最后根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本申请通过心脏心房心室分割模型,确定心房心室的解剖学结构的方式,实现了根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,从而达到了在心脏冠脉血管的3D影像的复杂背景中,可以连贯、准确地提取出血管信息并进行命名的技术效果,进而解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是根据本申请实施例的一种心脏冠脉血管的命名方法流程图;
图2是根据本申请实施例的一种可选的中心线树模型示意图;
图3是根据本申请实施例的一种可选的心脏心房心室分割模型示意图;
图4是根据本申请实施例的一种可选的房室沟交线示意图;
图5是根据本申请实施例的一种可选的房室沟平面示意图;
图6是根据本申请实施例的一种可选的中心点说明示意图;
图7是根据本申请实施例的一种可选的左右心室室间沟交线示意图;
图8是根据本申请实施例的一种可选的心脏冠脉血管命名的结果图;
图9是根据本申请实施例的一种可选的心脏冠脉血管命名的RCA示意图;
图10是根据本申请实施例的一种可选的冠脉树中心线LAD命名结果示意图;
图11是根据本申请实施例的一种可选的冠脉树中心线LCX命名结果示意图;
图12是根据本申请实施例的心脏冠脉血管的命名装置示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分的实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本申请的实施例。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
应该理解的是,当元件(诸如层、膜、区域、或衬底)描述为在另一元件“上”时,该元件可直接在该另一元件上,或者也可存在中间元件。而且,在说明书以及权利要求书中,当描述有元件“连接”至另一元件时,该元件可“直接连接”至该另一元件,或者通过第三元件“连接”至该另一元件。
实施例1
根据本申请的实施例,提供了一种心脏冠脉血管的命名方法。需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机***中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本申请实施例的心脏冠脉血管的命名方法流程图。如图1所示,该方法包括以下步骤:
步骤S101,获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树。
上述步骤中的心脏冠脉血管的3D影像可以通过CTA(Computed TomographicAngiography,计算机断层摄影血管造影术)得到,CTA是一种非创伤性的评价血管***的检查方法;心脏冠脉血管的中心线树可以是使用阈值分割方法进行提取,进一步地,在获取中心线树后可以使用骨架法对中心线树中的中心线进行提取,可以得到心脏冠状血管的中心线信息,示例性的,如图2所示为一种可选地中心线树模型示意图,需要说明的是,此处不对阈值分割方法做出具体的限定,同时中心线的提取也不仅限于骨架法,可以根据实际的需求选取合适的方法进行中心线树及中心线的提取。
步骤S102,获取心脏心房心室分割模型。
上述步骤中的心脏心房心室分割模型可以是使用AI分割方法,进行心脏心房心室的提取并建立心脏心房心室分割模型,示例性的,如图3所示为一种可选地心脏心房心室分割模型示意图。
步骤S103,根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面。
上述步骤中的心脏心房心室分割模型包括右心房右心室模型和左心房左心室模型,心房心室的解剖学结构可以通过确定左房室沟交线、左房室沟平面以及右房室沟交线、右房室沟平面确定,示例性的,如图4所示为房室沟交线示意图,图5所示为房室沟平面示意图。
步骤S104,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。
上述步骤中的相对位置关系可以包括延长线,分叉点,距离的最大值或最小值等,但需要说明的是,无论是哪一种位置关系都是基于中心线树和心房心室的解剖学结构进行说明。
应用本申请的技术方案,首先获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;并且获取心脏心房心室分割模型;之后根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;最后根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本申请通过心脏心房心室分割模型,确定心房心室的解剖学结构的方式,实现了根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,从而达到了在心脏冠脉血管的3D影像的复杂背景中,可以连贯、准确地提取出血管信息并进行命名的技术效果,进而解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
进一步地,心脏心房心室分割模型包括右心房右心室模型和左心房左心室模型,根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面,该方法包括:将右心房右心室模型扩大预定倍数,使得右心房和右心室产生第一交集;根据第一交集确定右房室沟交线;根据右房室沟交线确定右房室沟平面;将左心房左心室模型扩大预定倍数,使得左心房和左心室产生第二交集;根据第二交集确定左房室沟交线;根据左房室沟交线确定左房室沟平面。
具体地,右房室沟平面和左房室沟平面的法向朝向心房。
在一个可选的实施例中,首先将右心房右心室模型扩大至1.1倍,使得右心房右心室产生交集,即上述的第一交集,取右心房右心室的第一交集,得到右房室沟交线,基于右房室沟交线上所有点,使用相关算法进行拟合确定右房室沟平面,其中,相关算法可以使用最小二乘法,同样的方法应用于左心房左心室,将左心房右心室模型扩大至1.1倍,根据获取的第二交集确定左房室沟交线,并拟合确定左房室沟平面。
进一步地,根据心脏心房心室分割模型,确定心房心室的解剖学结构,该方法还包括:获取右房室沟交线的第一中心点和左房室沟交线的第二中心点;获取第一中心点和第二中心点的第一连线的第一法向量;获取右房室沟平面的第三中心点和左房室沟平面的第四中心点;获取第三中心点和第四中心点的第二连线的中点;根据第一法向量和中点构造十字交叉面。
在一个可选的实施例中,根据左右房室沟交线中点连线的法向量,即上述的第一法向量,以及左右房室沟平面中心连线的中点,即上述的第二连线的中点可以确定唯一的十字交叉面。
进一步地,获取右房室沟交线的第一中心点,该方法包括:获取右房室沟交线上的任意一个第一目标点,到除第一目标点以外的所有点的距离的第一平均值;确定最小的第一平均值对应的第一目标点为第一中心点。
在一个可选的实施例中,上述步骤中的第一中心点,如图5所示可知右房室沟交线为一个不规则的闭合曲线,为了进一步的说明,如图6所示,假设图中所示轨迹为某房室沟交线,理论上房室沟交线上是由无数个点构成的,也即,有无数个第一目标点,为了方便说明,取A、B、C三点进行描述,也即A、B、C三点为上述的第一目标点,在A、B、C点均可有很多个到除第一目标点以外的点的距离,图6中A、B、C点分别选择三个距离进行说明,图中由A点确定的第一平均值是通过计算三条实线的长度,得到一个平均值,同理,B、C点的第一平均值分别是根据由B、C点出发的三条虚线分别计算平均值,得到A、B、C三点的第一平均值后,比较三点的平均值的大小,最小的平均值对应的为右房室沟交线的第一中心点,同理,第二中心点也是通过同样的方式确定。
进一步地,获取右房室沟平面的第三中心点,该方法包括:获取右房室沟平面上的任意一个第二目标点,到除第二目标点以外的所有点的距离的第二平均值;确定最小的第二平均值对应的第二目标点为第三中心点。
上述步骤中的第三中心点的确定与上述第一中心点和第二中心点的确定有相似之处,不同的地方在于房室沟平面是一个平面,也即,在确定第三中心点时,若房室沟平面是一个规则的圆,则第三中心点可能位于平面的中心,同理,第四中心点和第三中心点的确定方法是一致的。
进一步地,根据心脏心房心室分割模型,确定心房心室的解剖学结构,该方法还包括:获取左右心室室间沟交线。
在一个可选的实施例中,如图7所示为左右心室室间沟交线示意图。
进一步地,中心线树包括左冠和右冠,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法包括:将右冠中穿过十字交叉面,并向后延伸最远的中心线命名为主PL;获取中心线树中除主PL以外的中心线,与主PL的分叉点;获取预定中心线的分叉点到末点的线段上的所有点,到左右心室室间沟交线的最小距离的平均值;获取预定中心线的分叉点到末点的线段上的所有点,到左右心室室间沟交线的距离的最小值;平均值和最小值的乘积最小,对应的预定中心线为PD。
在一个可选的实施例中,如图8所示,图8为心脏冠脉血管命名的结果图,十字交叉面垂直于图8,因为图片的镜像关系,右冠为图示的左半部分,可以根据十字交叉面判断主PL,右冠中穿过十字交叉面并向后延伸最远的中心线就是主PL,除主PL外的中心线与主PL的分叉点,计算被处理中心线分叉点后至末点的部分,与室间沟交线评分最小的为PD,图8中可以看到有两个分叉点,即PD1、PD2,其中,评分规则如下:分叉点后中心线上每一点到室间沟交线的最小距离的均值与中心线上所有点到室间沟交线的距离的最小值之积,即上述的平均值和最小值的乘积最小。
在另一个可选的实施例中,如果中心线分叉点到末点所成向量与右房室沟平面向量点积大于0.3,则为AVN,否则为PL。
进一步地,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法还包括:将右冠的起点至主PL与PD的分叉点部分,命名为RCA。
在一个可选的实施例中,基于图8所示,将右冠的起点至主PL与PD的分叉点部分,命名为RCA,如图9所示,图9为心脏冠脉血管命名的RCA示意图。
进一步地,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法还包括:将RCA分为第一前段、中段和第一后段;将从第一前段发出,并且RCA的分叉点和末点构成的向量的第二法向量与右房室沟平面的第三法向量的点积小于预定数值的部分,命名为CA;将从中段发出,并且RCA的分叉点和末点构成的向量的第四法向量与右房室沟平面的第五法向量的点积小于预定数值的部分,命名为Rv;将从第一后段发出,并且RCA的分叉点和末点构成的向量的第六法向量与右房室沟平面的第七法向量的点积小于预定数值的部分,命名为AM。
在一个可选的实施例中,仍以图9为例进行说明,将RCA再均分为前中后段,即上述的第一前段、中段和第一后段;第一前段发出并且与RCA分叉点至末点所成向量的法向与右房室沟平面法向的点积小于等于0.5,命名为CA;中段发出并且与RCA分叉点至末点所成向量的法向与右房室沟平面法向的点积小于等于0.5,命名为Rv,图中包含Rv1、Rv2;第一后段发出并且与RCA分叉点至末点所成向量的法向与右房室沟平面法向的点积小于等于0.5,命名为AM。
进一步地,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法还包括:将RCA分为第二前段和第二后段;将从第二前段发出,并且RCA的分叉点和末点构成的向量的第八法向量与右房室沟平面的第九法向量的点积小于预定数值的部分,命名为SAN;将从第二后段发出,并且RCA的分叉点和末点构成的向量的第十法向量与右房室沟平面的第十一法向量的点积小于预定数值的部分,命名为RA Branch。
在一个可选的实施例中,再将RCA均分为前后两段,即上述的第二前段和第二后段,第二前段发出并且与RCA分叉点至末点所成向量的法向与右房室沟平面法向的点积大于0.5,命名为SAN;第二后段发出并且与RCA分叉点至末点所成向量的法向与右房室沟平面法向的点积大于0.5,命名为RA Branch。
进一步地,中心线树包括左冠和右冠,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法包括:将左冠中与左右心室室间沟交线评分最小的中心线命名为LAD,其中,评分是根据中心线上的所有点到左右心室室间沟交线的最小距离的平均值,与中心线上的所有点到左右心室室间沟交线的距离的最小值之和,再除以中心线的长度得到的;获取左冠中的中心线中除LAD以外的中心线,与LAD的分叉点;确定分叉点点序最小的点为LM分叉点,其中,分叉点点序指的是所有除所述LAD以外的中心线与所述LAD中心线的分叉点排序得到分叉点点序;根据LM分叉点将左冠分为LAD冠脉树和LCX冠脉树。
在一个可选的实施例中,根据左冠中心线与室间沟交线评分最小的为LAD,评分规则如下:分叉点后中心线上每一点到室间沟交线的最小距离的均值与中心线上所有点到室间沟交线的距离的最小值之和,再除以中心线长度得到评分;再根据LAD与除LAD外其余所有左冠中心线求分叉点,得到分叉点点序最小的点为LM分叉点,由此可得到LAD冠脉树和LCX冠脉树。
进一步地,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法还包括:将LAD上从LM分叉点起至LAD末点的所有点,拟合为LAD平面;计算LAD冠脉树中的中心线与LAD的分叉点与中心线末点构成的向量,与LAD平面的法向量的点积;将点积小于第一预定值的中心线命名为D,将点积大于或者等于第一预定值的中心线命名为S;中心线若为D,且中心线与LAD的分叉点与LM分叉点之间的距离小于预定距离的情况下,将中心线命名为RM。
在一个可选的实施例中,如图10所示,图10为冠脉树中心线LAD命名结果示意图,LAD上从LM分叉点起至LAD末点部分中心线点使用最小二乘法拟合平面,简称LAD平面。计算LAD冠脉树中心线与LAD的分叉点开始,至中心线末点所成向量与LAD平面法向点积,小于0.1则为D,如图10所示D1、D2、D3、D4,否则为S,如图10所示S1、S2,其中,若为D且与LAD分叉点到LM分叉点距离小于0.5mm,则命名为RM。
进一步地,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,该方法还包括:将LCX冠脉树中与左房室沟交线评分最小的中心线命名为主LCX,评分规则为:LM分叉点后中心线上的所有点到左房室沟交线的最小距离的平均值,与LM分叉点后中心线上的所有点到左房室沟交线的距离的最小值之和,再除以中心线的长度得到评分;计算LCX冠脉树中除主LCX以外的中心线,与主LCX的分叉点至中心线末点构成的向量,与左房室沟平面的法向量的点积;将点积小于第二预定值的中心线命名为OM,将点积大于第三预定值的中心线命名为LA Branch,否则,命名为LCX,其中,第三预定值大于第二预定值;在中心线为OM的情况下,且中心线与主LCX的分叉点与LM分叉点之间的距离小于预定距离的情况下,将中心线命名为RM。
在一个可选的实施例中,如图11所示,图11为冠脉树中心线LCX命名结果示意图,根据中心线与左房室沟交线评分最小的为主LCX,评分规则如下:LM分叉点后中心线上每一点到左房室沟交线的最小距离的均值与中心线上所有点到左房室沟交线的距离的最小值之和,再除以中心线长度得到评分,得到主LCX后,判断其余中心线与主LCX分叉点,并得到分叉点至末点所成向量。计算该向量与左房室沟平面法向点积。若点积小于-0.2,则为OM,如图11所示OM1、OM2、OM3、OM4、OM5;若点积大于0.1,则为LA Branch,否则为LCX。其中,若命名为OM,且分叉点与LM分叉点距离小于0.5mm,则命名为RM。
实施例2
本申请实施例还提供了一种心脏冠脉血管的命名装置,需要说明的是,本申请实施例的心脏冠脉血管的命名装置可以用于执行本申请实施例所提供的用于心脏冠脉血管的命名方法。以下对本申请实施例提供的心脏冠脉血管的命名装置进行介绍。
图12是根据本申请实施例的心脏冠脉血管的命名装置示意图。如图12所示,该装置包括:
第一获取单元1201,用于获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;
第二获取单元1202,用于获取心脏心房心室分割模型;
确定单元1203,用于根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;
命名单元1204,用于根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。
应用本申请的技术方案,通过第一获取单元获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;通过第二获取单元获取心脏心房心室分割模型;通过确定单元根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;通过命名单元根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本申请通过心脏心房心室分割模型,确定心房心室的解剖学结构的方式,实现了根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,从而达到了在心脏冠脉血管的3D影像的复杂背景中,可以连贯、准确地提取出血管信息并进行命名的技术效果,进而解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
进一步地,确定单元包括:第一扩大模块,用于将右心房右心室模型扩大预定倍数,使得右心房和右心室产生第一交集;第一确定模块,用于根据第一交集确定右房室沟交线;第二确定模块,用于根据右房室沟交线确定右房室沟平面;第二扩大模块,用于将左心房左心室模型扩大预定倍数,使得左心房和左心室产生第二交集;第三确定模块,用于根据第二交集确定左房室沟交线;第四确定模块,用于根据左房室沟交线确定左房室沟平面。
进一步地,确定单元还包括:第一获取模块,用于获取右房室沟交线的第一中心点和左房室沟交线的第二中心点;第二获取模块,用于获取第一中心点和第二中心点的第一连线的第一法向量;第三获取模块,用于获取右房室沟平面的第三中心点和左房室沟平面的第四中心点;第四获取模块,用于获取第三中心点和第四中心点的第二连线的中点;构造模块,根据第一法向量和中点构造十字交叉面。
进一步地,第一获取模块包括:第一获取子模块,用于获取右房室沟交线上的任意一个第一目标点,到除第一目标点以外的所有点的距离的第一平均值;第一确定子模块,用于确定最小的第一平均值对应的第一目标点为第一中心点。
进一步地,第三获取模块包括:第二获取子模块,用于获取右房室沟平面上的的任意一个第二目标点,到除第二目标点以外的所有点的距离的第二平均值;第二确定子模块,用于确定最小的第二平均值对应的第二目标点为第三中心点。
进一步地,确定单元还包括:第五获取模块,用于获取左右心室室间沟交线。
进一步地,命名单元包括:第一命名模块,用于将右冠中穿过十字交叉面,并向后延伸最远的中心线命名为主PL;第六获取模块,用于获取中心线树中除主PL以外的中心线,与主PL的分叉点;第七获取模块,用于获取预定中心线的分叉点到末点的线段上的所有点,到左右心室室间沟交线的最小距离的平均值;第八获取模块,用于获取预定中心线的分叉点到末点的线段上的所有点,到左右心室室间沟交线的距离的最小值;第二命名模块,用于平均值和最小值的乘积最小,对应的预定中心线为PD。
进一步地,命名单元还包括:第三命名模块,用于将右冠的起点至主PL与PD的分叉点部分,命名为RCA。
进一步地,命名单元还包括:第一分段模块,用于将RCA分为第一前段、中段和第一后段;第四命名模块,用于将从第一前段发出,并且RCA的分叉点和末点构成的向量的第二法向量与右房室沟平面的第三法向量的点积小于预定数值的部分,命名为CA;第五命名模块,用于将从中段发出,并且RCA的分叉点和末点构成的向量的第四法向量与右房室沟平面的第五法向量的点积小于预定数值的部分,命名为Rv;第六命名模块,用于将从第一后段发出,并且RCA的分叉点和末点构成的向量的第六法向量与右房室沟平面的第七法向量的点积小于预定数值的部分,命名为AM。
进一步地,命名单元还包括:第二分段模块,用于将RCA分为第二前段和第二后段;第七命名模块,用于将从第二前段发出,并且RCA的分叉点和末点构成的向量的第八法向量与右房室沟平面的第九法向量的点积小于预定数值的部分,命名为SAN;第八命名模块,用于将从第二后段发出,并且RCA的分叉点和末点构成的向量的第十法向量与右房室沟平面的第十一法向量的点积小于预定数值的部分,命名为RA Branch。
进一步地,命名单元还包括:第九命名模块,用于将左冠中与左右心室室间沟交线评分最小的中心线命名为LAD,其中,评分是根据中心线上的所有点到左右心室室间沟交线的最小距离的平均值,与中心线上的所有点到左右心室室间沟交线的距离的最小值之和,再除以中心线的长度得到的;第九获取模块,用于获取左冠中的中心线中除LAD以外的中心线,与LAD的分叉点;第五确定模块,用于确定分叉点点序最小的点为LM分叉点,其中,分叉点点序指的是所有除所述LAD以外的中心线与所述LAD中心线的分叉点排序得到分叉点点序;分类模块,用于根据LM分叉点将左冠分为LAD冠脉树和LCX冠脉树。
进一步地,命名单元还包括:拟合模块,用于将LAD上从LM分叉点起至LAD末点的所有点,拟合为LAD平面;第一计算模块,用于计算LAD冠脉树中的中心线与LAD的分叉点与中心线末点构成的向量,与LAD平面的法向量的点积;第九命名模块,用于将点积小于第一预定值的中心线命名为D,将点积大于或者等于第一预定值的中心线命名为S;第十命名模块,用于中心线若为D,且中心线与LAD的分叉点与LM分叉点之间的距离小于预定距离的情况下,将中心线命名为RM。
进一步地,命名单元还包括:第十一命名模块,用于将LCX冠脉树中与左房室沟交线评分最小的中心线命名为主LCX,评分规则为:LM分叉点后中心线上的所有点到左房室沟交线的最小距离的平均值,与LM分叉点后中心线上的所有点到左房室沟交线的距离的最小值之和,再除以中心线的长度得到评分;第二计算模块,用于计算LCX冠脉树中除主LCX以外的中心线,与主LCX的分叉点至中心线末点构成的向量,与左房室沟平面的法向量的点积;第十二命名模块,用于将点积小于第二预定值的中心线命名为OM,将点积大于第三预定值的中心线命名为LA Branch,否则,命名为LCX,其中,第三预定值大于第二预定值;第十三命名模块,用于在中心线为OM的情况下,且中心线与主LCX的分叉点与LM分叉点之间的距离小于预定距离的情况下,将中心线命名为RM。
所述心脏冠脉血管的命名装置包括处理器和存储器,上述第一获取单元、第二获取单元、确定单元、命名单元等均作为程序单元存储在存储器中,由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
处理器中包含内核,由内核去存储器中调取相应的程序单元。内核可以设置一个或以上,通过调整内核参数来实现根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,从而达到了在心脏冠脉血管的3D影像的复杂背景中,可以连贯、准确地提取出血管信息并进行命名的技术效果,进而解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM),存储器包括至少一个存储芯片。
本发明实施例提供了一种计算机可读存储介质,所述计算机可读存储介质包括存储的程序,其中,在所述程序运行时控制所述计算机可读存储介质所在设备执行所述心脏冠脉血管的命名方法。
本发明实施例提供了一种处理器,所述处理器用于运行程序,其中,所述程序运行时执行所述心脏冠脉血管的命名方法。
本发明实施例提供了一种设备,设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序,处理器执行程序时实现至少以下步骤:步骤S101,获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;步骤S102,获取心脏心房心室分割模型;步骤S103,根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;步骤S104,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本文中的设备可以是服务器、PC、PAD、手机等。
本申请还提供了一种计算机程序产品,当在数据处理设备上执行时,适于执行初始化有至少如下方法步骤的程序:步骤S101,获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;步骤S102,获取心脏心房心室分割模型;步骤S103,根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;步骤S104,根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
从以上的描述中,可以看出,本申请上述的实施例实现了如下技术效果:
1)、应用本申请的技术方案,首先获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;并且获取心脏心房心室分割模型;之后根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;最后根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本申请通过心脏心房心室分割模型,确定心房心室的解剖学结构的方式,实现了根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,从而达到了在心脏冠脉血管的3D影像的复杂背景中,可以连贯、准确地提取出血管信息并进行命名的技术效果,进而解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
2)、应用本申请的技术方案,通过第一获取单元获取心脏冠脉血管的3D影像,且从3D影像中提取出心脏冠脉血管的中心线树;通过第二获取单元获取心脏心房心室分割模型;通过确定单元根据心脏心房心室分割模型,确定心房心室的解剖学结构,解剖学结构至少包括房室沟交线和房室沟平面;通过命名单元根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名。本申请通过心脏心房心室分割模型,确定心房心室的解剖学结构的方式,实现了根据中心线树和心房心室的解剖学结构的相对位置关系,为心脏冠脉血管进行命名,从而达到了在心脏冠脉血管的3D影像的复杂背景中,可以连贯、准确地提取出血管信息并进行命名的技术效果,进而解决了现有技术中对心脏冠脉血管进行命名的方***度较低的技术问题。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (14)
1.一种心脏冠脉血管的命名方法,其特征在于,包括:
获取心脏冠脉血管的3D影像,且从所述3D影像中提取出所述心脏冠脉血管的中心线树;
获取心脏心房心室分割模型;
根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,所述解剖学结构至少包括房室沟交线和房室沟平面;
根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名;
其中,右房室沟交线为右心房和右心室的交线、左房室沟交线为左心房和左心室的交线、右房室沟平面为右房室沟交线拟合得到的平面以及左房室沟平面为左房室沟交线拟合得到的平面,根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,包括:
获取所述右房室沟交线的第一中心点和所述左房室沟交线的第二中心点;
获取所述第一中心点和所述第二中心点的第一连线的第一法向量;
获取所述右房室沟平面的第三中心点和所述左房室沟平面的第四中心点;
获取所述第三中心点和所述第四中心点的第二连线的中点;
根据所述第一法向量和所述中点构造十字交叉面;
其中,所述中心线树包括左冠和右冠,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,包括:
将所述右冠中穿过所述十字交叉面,并向后延伸最远的中心线命名为主PL;
获取所述中心线树中除所述主PL以外的中心线,与所述主PL的分叉点;
获取预定中心线的分叉点到末点的线段上的所有点,到左右心室室间沟交线的最小距离的平均值;
获取所述预定中心线的分叉点到末点的线段上的所有点,到所述左右心室室间沟交线的距离的最小值;
所述平均值和所述最小值的乘积最小,对应的所述预定中心线为PD。
2.根据权利要求1所述的方法,其特征在于,所述心脏心房心室分割模型包括右心房右心室模型和左心房左心室模型,根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,所述解剖学结构至少包括房室沟交线和房室沟平面,包括:
将所述右心房右心室模型扩大预定倍数,使得右心房和右心室产生第一交集;
根据所述第一交集确定右房室沟交线;
根据所述右房室沟交线确定右房室沟平面;
将所述左心房左心室模型扩大所述预定倍数,使得左心房和左心室产生第二交集;
根据所述第二交集确定左房室沟交线;
根据所述左房室沟交线确定左房室沟平面。
3.根据权利要求1所述的方法,其特征在于,获取所述右房室沟交线的第一中心点,包括:
获取所述右房室沟交线上的任意一个第一目标点,到除所述第一目标点以外的所有点的距离的第一平均值;
确定最小的所述第一平均值对应的所述第一目标点为所述第一中心点。
4.根据权利要求1所述的方法,其特征在于,获取所述右房室沟平面的第三中心点,包括:
获取所述右房室沟平面上的任意一个第二目标点,到除所述第二目标点以外的所有点的距离的第二平均值;
确定最小的所述第二平均值对应的所述第二目标点为所述第三中心点。
5.根据权利要求1所述的方法,其特征在于,根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,还包括:
获取左右心室室间沟交线。
6.根据权利要求1所述的方法,其特征在于,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,还包括:
将所述右冠的起点至所述主PL与所述PD的分叉点部分,命名为RCA。
7.根据权利要求6所述的方法,其特征在于,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,还包括:
将所述RCA分为第一前段、中段和第一后段;
将从所述第一前段发出,并且所述RCA的分叉点和末点构成的向量的第二法向量与所述右房室沟平面的第三法向量的点积小于预定数值的部分,命名为CA;
将从所述中段发出,并且所述RCA的分叉点和末点构成的向量的第四法向量与所述右房室沟平面的第五法向量的点积小于所述预定数值的部分,命名为Rv;
将从所述第一后段发出,并且所述RCA的分叉点和末点构成的向量的第六法向量与所述右房室沟平面的第七法向量的点积小于所述预定数值的部分,命名为AM。
8.根据权利要求6所述的方法,其特征在于,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,还包括:
将所述RCA分为第二前段和第二后段;
将从所述第二前段发出,并且所述RCA的分叉点和末点构成的向量的第八法向量与所述右房室沟平面的第九法向量的点积小于预定数值的部分,命名为SAN;
将从所述第二后段发出,并且所述RCA的分叉点和末点构成的向量的第十法向量与所述右房室沟平面的第十一法向量的点积小于预定数值的部分,命名为RA Branch。
9.根据权利要求5所述的方法,其特征在于,所述中心线树包括左冠和右冠,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,包括:
将所述左冠中与所述左右心室室间沟交线评分最小的中心线命名为LAD,其中,所述评分是根据中心线上的所有点到所述左右心室室间沟交线的最小距离的平均值,与中心线上的所有点到所述左右心室室间沟交线的距离的最小值之和,再除以中心线的长度得到的;
获取所述左冠中的中心线中除所述LAD以外的中心线,与所述LAD的分叉点;
确定分叉点点序最小的点为LM分叉点,其中,所述分叉点点序指的是所有除所述LAD以外的中心线与所述LAD中心线的分叉点排序得到分叉点点序;
根据所述LM分叉点将所述左冠分为LAD冠脉树和LCX冠脉树。
10.根据权利要求9所述的方法,其特征在于,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,还包括:
将所述LAD上从所述LM分叉点起至所述LAD末点的所有点,拟合为LAD平面;
计算所述LAD冠脉树中的中心线与所述LAD的分叉点与中心线末点构成的向量,与所述LAD平面的法向量的点积;
将所述点积小于第一预定值的中心线命名为D,将所述点积大于或者等于第一预定值的中心线命名为S;
所述中心线若为D,且所述中心线与所述LAD的分叉点与所述LM分叉点之间的距离小于预定距离的情况下,将所述中心线命名为RM。
11.根据权利要求9所述的方法,其特征在于,根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名,还包括:
将所述LCX冠脉树中与所述左房室沟交线评分最小的中心线命名为主LCX,评分规则为:所述LM分叉点后中心线上的所有点到所述左房室沟交线的最小距离的平均值,与所述LM分叉点后中心线上的所有点到所述左房室沟交线的距离的最小值之和,再除以中心线的长度得到所述评分;
计算所述LCX冠脉树中除所述主LCX以外的中心线,与所述主LCX的分叉点至中心线末点构成的向量,与所述左房室沟平面的法向量的点积;
将所述点积小于第二预定值的中心线命名为OM,将所述点积大于第三预定值的中心线命名为LA Branch,否则,命名为LCX,其中,所述第三预定值大于所述第二预定值;
在所述中心线为OM的情况下,且所述中心线与所述主LCX的分叉点与所述LM分叉点之间的距离小于预定距离的情况下,将所述中心线命名为RM。
12.一种心脏冠脉血管的命名装置,其特征在于,包括:
第一获取单元,用于获取心脏冠脉血管的3D影像,且从所述3D影像中提取出所述心脏冠脉血管的中心线树;
第二获取单元,用于获取心脏心房心室分割模型;
确定单元,用于根据所述心脏心房心室分割模型,确定心房心室的解剖学结构,所述解剖学结构至少包括房室沟交线和房室沟平面;
命名单元,用于根据所述中心线树和所述心房心室的解剖学结构的相对位置关系,为所述心脏冠脉血管进行命名;
其中,右房室沟交线为右心房和右心室的交线、左房室沟交线为左心房和左心室的交线、右房室沟平面为右房室沟交线拟合得到的平面以及左房室沟平面为左房室沟交线拟合得到的平面,所述确定单元包括:
第一获取模块,用于获取所述右房室沟交线的第一中心点和所述左房室沟交线的第二中心点;
第二获取模块,用于获取所述第一中心点和所述第二中心点的第一连线的第一法向量;
第三获取模块,用于获取所述右房室沟平面的第三中心点和所述左房室沟平面的第四中心点;
第四获取模块,用于获取所述第三中心点和所述第四中心点的第二连线的中点;
构造模块,用于根据所述第一法向量和中点构造十字交叉面;
其中,所述中心线树包括左冠和右冠,所述命名单元包括:
第一命名模块,用于将所述右冠中穿过所述十字交叉面,并向后延伸最远的中心线命名为主PL;
第六获取模块,用于获取所述中心线树中除所述主PL以外的中心线,与所述主PL的分叉点;
第七获取模块,用于获取预定中心线的分叉点到末点的线段上的所有点,到左右心室室间沟交线的最小距离的平均值;
第八获取模块,用于获取所述预定中心线的分叉点到末点的线段上的所有点,到所述左右心室室间沟交线的距离的最小值;
第二命名模块,用于所述平均值和所述最小值的乘积最小,对应的所述预定中心线为PD。
13.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质包括存储的程序,其中,在所述程序运行时控制所述计算机可读存储介质所在设备执行权利要求1至11中任意一项所述的方法。
14.一种处理器,其特征在于,所述处理器用于运行程序,其中,所述程序运行时执行权利要求1至11中任意一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111545244.4A CN114418947B (zh) | 2021-12-16 | 2021-12-16 | 心脏冠脉血管的命名方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111545244.4A CN114418947B (zh) | 2021-12-16 | 2021-12-16 | 心脏冠脉血管的命名方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114418947A CN114418947A (zh) | 2022-04-29 |
CN114418947B true CN114418947B (zh) | 2023-11-03 |
Family
ID=81268491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111545244.4A Active CN114418947B (zh) | 2021-12-16 | 2021-12-16 | 心脏冠脉血管的命名方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114418947B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110706770A (zh) * | 2019-09-30 | 2020-01-17 | 上海杏脉信息科技有限公司 | 心脏数据处理设备及处理方法、计算机可读存储介质 |
CN111292314A (zh) * | 2020-03-03 | 2020-06-16 | 上海联影智能医疗科技有限公司 | 冠脉分段方法、装置、图像处理***和存储介质 |
CN111369525A (zh) * | 2020-03-02 | 2020-07-03 | 联影智能医疗科技(北京)有限公司 | 图像分析方法、设备和存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7813785B2 (en) * | 2003-07-01 | 2010-10-12 | General Electric Company | Cardiac imaging system and method for planning minimally invasive direct coronary artery bypass surgery |
-
2021
- 2021-12-16 CN CN202111545244.4A patent/CN114418947B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110706770A (zh) * | 2019-09-30 | 2020-01-17 | 上海杏脉信息科技有限公司 | 心脏数据处理设备及处理方法、计算机可读存储介质 |
CN111369525A (zh) * | 2020-03-02 | 2020-07-03 | 联影智能医疗科技(北京)有限公司 | 图像分析方法、设备和存储介质 |
CN111292314A (zh) * | 2020-03-03 | 2020-06-16 | 上海联影智能医疗科技有限公司 | 冠脉分段方法、装置、图像处理***和存储介质 |
Non-Patent Citations (1)
Title |
---|
李治安 ; 王新房 ; 杨娅 ; 刘俐 ; 吕清 ; 谢明星 ; .冠脉解剖和冠脉血流显像对照研究.医学研究通讯.2000,(10),第13-15页. * |
Also Published As
Publication number | Publication date |
---|---|
CN114418947A (zh) | 2022-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102016959B1 (ko) | 혈관 모델 생성 방법 및 장치 | |
JP6981981B2 (ja) | 心臓モデルガイド付き冠動脈瘤セグメンテーション | |
US10241968B2 (en) | System and method for real-time simulation of patient-specific cardiac electrophysiology including the effect of the electrical conduction system of the heart | |
US11948695B2 (en) | Systems and methods for identifying and modeling unresolved vessels in image-based patient-specific hemodynamic models | |
US10282846B2 (en) | Model-based segmentation of an anatomical structure | |
CN108564574B (zh) | 确定血流储备分数的方法、计算机设备及计算机可读存储介质 | |
Queirós et al. | Fast left ventricle tracking using localized anatomical affine optical flow | |
CN117115150B (zh) | 用于确定分支血管的方法、计算设备和介质 | |
CN114418947B (zh) | 心脏冠脉血管的命名方法及装置 | |
RU2719951C2 (ru) | Моделирование анатомической структуры методом конечных элементов | |
US20230100255A1 (en) | System and method for interactive contouring of medical images | |
CN109155065B (zh) | 用于诊断图像分析和图像质量评估的***及方法 | |
JP6748368B2 (ja) | モデリング装置、モデリング方法、およびモデリングプログラム | |
CN116168810A (zh) | 用于心内超声心动图的智能图像导航 | |
Goyal et al. | MRI image based patient specific computational model reconstruction of the left ventricle cavity and myocardium | |
US10236083B2 (en) | Visualization apparatus and visualization method | |
Veress et al. | The direct incorporation of perfusion defect information to define ischemia and infarction in a finite element model of the left ventricle | |
JP6748370B2 (ja) | 生体モデル生成装置、生体モデル生成方法、および生体モデル生成プログラム | |
Kubik et al. | Novel 3D heart left ventricle muscle segmentation method for PET-gated protocol and its verification | |
Verhülsdonk et al. | Shape of my heart: Cardiac models through learned signed distance functions | |
Frank | Deep Learning Approaches to Multi-Modal Biomedical Image Segmentation | |
Cheng et al. | Elliptic Shape Prior Dynamic Programming for Accurate Vessel Segmentation in MRI Sequences with Automated Optimal Parameter Selection | |
Mountris et al. | Image-based Cardiac Electrophysiology Simulation through the Meshfree Mixed Collocation Method | |
KR20230147467A (ko) | 심방세동 재발을 예측하는 장치 및 예측 방법 | |
Giachetti et al. | Measurable models of abdominal aortic aneurysm on the web |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |