CN107392872B - 一种微导管塑形器的生成方法及生成系统 - Google Patents

一种微导管塑形器的生成方法及生成系统 Download PDF

Info

Publication number
CN107392872B
CN107392872B CN201710628619.0A CN201710628619A CN107392872B CN 107392872 B CN107392872 B CN 107392872B CN 201710628619 A CN201710628619 A CN 201710628619A CN 107392872 B CN107392872 B CN 107392872B
Authority
CN
China
Prior art keywords
blood vessel
point
path
micro
catheter
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
Application number
CN201710628619.0A
Other languages
English (en)
Other versions
CN107392872A (zh
Inventor
杨光明
李建军
秦岚
刘文哲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qianglian Zhichuang Beijing Technology Co ltd
Original Assignee
Qianglian Zhichuang Beijing Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qianglian Zhichuang Beijing Technology Co ltd filed Critical Qianglian Zhichuang Beijing Technology Co ltd
Priority to CN201710628619.0A priority Critical patent/CN107392872B/zh
Publication of CN107392872A publication Critical patent/CN107392872A/zh
Application granted granted Critical
Publication of CN107392872B publication Critical patent/CN107392872B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M25/00Catheters; Hollow probes
    • A61M25/0009Making of catheters or other medical or surgical tubes
    • G06T5/70

Abstract

本发明涉及一种微导管塑形器的生成方法及生成系统,生成方法包括:根据三维血管图像得到三维血管中心轴;建立连接网络;读取用户选取的原点和末点坐标,计算最短路径;对三维血管中心轴周围的像素进行处理,得到初步截取的血管;根据初步截取的血管,获得血管表面;根据得到的血管表面,计算得到微导管折线路径;对微导管折线路径进行平滑处理,得到微导管路径;根据微导管路径上的折线点,得到塑形器的折线路径;对塑形器的折线路径进行平滑处理,得到表示微导管塑形器的二元三维矩阵;根据微导管塑形器的二元三维矩阵,得到微导管塑形器的表面文件。本发明能够有效地指导设计人员设计微导管塑形器的工作,提高设计人员的工作效率。

Description

一种微导管塑形器的生成方法及生成系统
技术领域
本发明属于医疗器械领域,具体涉及一种微导管塑形器的生成方法及生成系统。
背景技术
微导管是介入手术治疗中常用的一种器械。例如,在颅内动脉瘤的弹簧圈介入栓塞手术中,首先要将相应微导管选择性送达动脉瘤内。该手术过程中很重要的一个环节就是微导管的成功塑形。对微导管前端的形态进行良好的塑形,能够大大提高介入手术中微导管的到位准确性、微导管在栓塞过程中的稳定性和微导管的操控灵活性。成功的微导管塑形是手术成功的基本保证,微导管塑形对于床突旁及微小动脉瘤的介入治疗尤其重要。
然而在实际的手术过程当中,微导管的塑形成功与否与医生的技术水平和临床经验直接相关。虽然颅内动脉瘤栓塞手术中的微导管塑形已有不少国内外临床报告和学术论文进行了深入讨论与总结,脑血管和颅内动脉瘤的三维空间形态可通过数字血管减影技术和3D重建技术获得,血管的精确尺寸也可以通过三维图像测量获得,微导管的塑形还是需要依赖医生丰富的知识和经验,而临床上一直以来都缺少精确有效的辅助设计手段。
发明内容
为了解决现有技术存在的上述问题,本发明提供了一种微导管塑形器的生成方法及生成系统。
本发明所采用的技术方案为:一种微导管塑形器的生成方法包括以下步骤:
根据三维血管图像得到三维血管中心轴;
根据三维血管中心轴上各点与其周围所有点的相邻关系,建立连接网络;
读取用户选取的原点和末点坐标,计算沿三维血管中心轴从与原点距离最近的点到与末点距离最近的点的最短路径;
对三维血管中心轴周围的像素进行处理,得到初步截取的血管;
根据初步截取的血管,获得临界点坐标,将所有临界点坐标相连得到血管表面;
根据得到的血管表面,计算得到微导管折线路径;
对微导管折线路径进行平滑处理,得到微导管路径;
根据微导管路径上的折线点,得到塑形器的折线路径;
对塑形器的折线路径进行平滑处理,得到表示微导管塑形器的二元三维矩阵;
根据微导管塑形器的二元三维矩阵,得到微导管塑形器的表面文件。
进一步地,所述根据三维血管图像得到三维血管中心轴的步骤中,采用图像侵蚀算法对表示三维血管图像的二元三维矩阵重复进行对称的图像侵蚀处理,直到得到三维血管中心轴。
进一步地,所述根据三维血管中心轴上各点与其周围所有点的相邻关系建立连接网络的步骤中,预设三维血管中心轴上两个点之间的最小直线距离为L,如果三维血管中心轴上两个点之间的距离小于或等于
Figure BDA0001362412030000021
则判定这两个点为相邻关系。
进一步地,所述计算沿三维血管中心轴从与原点距离最近的点到与末点距离最近的点的最短路径的具体过程为:
分别计算三维血管中心轴上各个点与读取的原点和末点之间的距离,得到三维血管中心轴上与原点之间距离最近的点以及三维血管中心轴上与末点之间距离最近的点;其中,三维血管中心轴上与原点之间距离最近的点为起点,三维血管中心轴上与末点之间距离最近的点为终点,
根据得到的连接网络,采用迪杰斯特拉算法得到沿三维血管中心轴从起点到终点的最短路径。
进一步地,所述对三维血管中心轴周围的像素进行处理得到初步截取的血管的具体过程为:
根据获取的最短路径,得到最短路径上各个点的坐标;
对于二元三维矩阵的每个非0的血管像素点,定义该像素点与最短路径上各个点沿最短路径的距离为该像素点的坐标(x,y,z)同最短路径上各个点坐标的距离的最小值;
计算二元三维矩阵的每个非0的血管像素点与最短路径上各个点坐标之间的距离;
预设距离临界值,如果计算得到的距离大于距离临界值,则将该距离对应的血管像素值设为背景值0,以此将与沿最短路径较远的血管点设为背景值,仅保留沿最短路径附近的血管,得到初步截取的血管。
进一步地,所述得到血管表面的具体过程为:
以起点为种子,采用区域生长的方法,对与种子不相连但在得到初步截取的血管的过程中保留的像素点进行清除,得到新的二元三维矩阵;
在新的二元三维矩阵的基础上,获得满足条件的临界点坐标,将所有临界点坐标相连得到血管表面;其中,临界点坐标满足的条件包括与种子相连且与最短路径上各个点坐标之间的距离小于距离临界值。
进一步地,所述根据得到的血管表面得到微导管折线路径的具体过程为:
以通过起点的沿最短路径方向为起始方向,以直线的方式延伸路径,方向记为
Figure BDA0001362412030000031
如果路径抵达血管壁,则记录下对应血管壁的坐标W,找到与坐标W距离最小的沿最短路径上的点C,求得方向
Figure BDA0001362412030000033
相对于点C的对称方向
Figure BDA0001362412030000032
并保证从坐标W开始,延方向
Figure BDA0001362412030000034
延伸会使得路径朝向血管内部,直到再次碰到血管壁;
重复上述过程,直到路径与终点之间的距离小于距离临界值,得到微导管折线路径。
进一步地,所述根据微导管路径上的折线点得到塑形器的折线路径的具体过程为:
根据微导管路径上的折线点,计算出每段线段的长度及方向,其中方向由三维单元向量表示,并计算出由两段线段决定的坐标系到下两段线段决定的坐标系的旋转矩阵;
利用旋转矩阵得到相邻线段组之间的角度;每个线段组中包含两条相邻的线段;
假定微导管塑形器的角度与微导管路径的折线角度有固定的倍数关系α,根据折线点的坐标和计算得到的旋转矩阵,从第一条线段开始,将微导管路径旋转矩阵对应的旋转角度乘以给定的倍数α,得到塑形器下一条线段相对于前一条线段的旋转矩阵,并依据该矩阵得到塑形器下一条线段的方向,同时保证塑形器下一条线段的长度同微导管路径对应线段的长度相同,由此确定下一条线段,直到得到塑形器的全部折线路径。
进一步地,所述得到微导管塑形器的表面文件的具体过程为:
假定微导管的半径为r1,塑形器的半径为r2;
对每个背景像素,计算其坐标(x,y,z)与微导管路径上每个点的距离,对得到的距离值进行比较后得到最小距离值d;
如果d<=r1,则将背景像素的像素值设为1;如果r1<d<=r2,则将背景像素的像素值设为2;如果d>r2,则将背景像素的像素值设为0;
选取像素值为0和2的表面作为微导管塑形器的外表面,像素值为1和2的表面作为微导管塑形器的内表面,得到微导管塑形器的表面文件。
一种微导管塑形器的生成系统包括三维血管中心轴生成模块、血管截取模块、血管表面生成模块、微导管路径生成模块、塑形器路径生成模块和塑形器表面文件生成模块;
所述三维血管中心轴生成模块用于根据三维血管数据生成三维血管中心轴;所述血管截取模块用于根据三维血管中心轴上各点之间的连接关系对三维血管中心轴周围的像素进行处理,得到截取的血管;所述血管表面生成模块用于根据截取的血管生成血管表面文件;所述微导管路径生成模块用于根据血管表面文件处理得到微导管路径;所述塑形器路径生成模块用于根据微导管路径处理得到微导管塑形器的路径;所述塑形器表面文件生成模块用于根据微导管塑形器的路径生成微导管塑形器的表面文件。
由于采用以上技术方案,本发明的有益效果为:本发明根据三维血管图像 得到初步截取的血管,采用Marching Cube算法得到血管表面,进而得到微导管路径,根据微导管路径得到微导管塑形器的表面文件,采用本发明能够自动生成微导管塑形器的表面文件,进而为塑形器的制作提供设计参考。本发明能够有效地指导设计人员设计微导管塑形器的工作,提高设计人员的工作效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。显而易见,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一实施例提供的一种微导管塑形器的生成方法的流程图;
图2是本发明一实施例提供的一种微导管塑形器的生成系统的结构框图。
图中:1-三维血管中心轴生成模块;2-血管截取模块;3-血管表面生成模块;4-微导管路径生成模块;5-塑形器路径生成模块;6-塑形器表面文件生成模块。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行详细的描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式,都属于本发明所保护的范围。
微导管在神经介入手术中发挥着相当重要的作用,在治疗血管瘤和血管栓塞中,针对不同病人的病变情况和血管束设计出最合适的微导管路径及微导管塑形器能够大大提高手术的效率和治疗效果。
手术前利用高分辨率的DSA(digital subtraction angiography,数字减影心血管造影术)、CTA(CT血管造影术)或MRA(Magnetic Resonance Angiography,磁共振血管成像),通过图像分割的方法,可以得到完整的三维血管图像。三维血管图像通常是以二元三维矩阵的形式进行表示,其中,0表示背景,1表示 血管及血管内部血液。
基于完整的三维血管图像,如图1所示,本发明提供了一种微导管塑形器的生成方法,其包括以下步骤:
S1、根据三维血管图像得到三维血管中心轴。
以三维血管数据为基准,参照人体坐标建立三维坐标系。其中,X轴对应人体左右方向,Y轴对应人体前后方向,Z轴对应人体头脚方向。三维血管图像中每个像素点对应三维坐标系的整数点,例如第一个像素点对应的三维坐标为(1,1,1)。采用图像侵蚀算法对表示三维血管图像的二元三维矩阵重复进行对称的图像侵蚀处理,直到得到三维血管中心轴。
图像侵蚀算法能够在保证血管连接性不变的前提下,将血管缩小成多条相互连接的线段。由于图像侵蚀算法在缩小的过程中保持对称性,因此得到的线段为血管束的中心线。血管束中心线以中心线上所有点的坐标来表示,由于每个像素点对应整数点三维坐标,血管束中心线上的每个点也对应整数点三维坐标。
S2、根据三维血管中心轴上各点与其周围所有点的相邻关系,建立连接网络。
预设三维血管中心轴上两个点之间的最小直线距离为L,如果三维血管中心轴上两个点之间的距离小于或等于
Figure BDA0001362412030000061
则判定这两个点为相邻关系。L通常取值为1。例如,三维血管中心轴上一个点的坐标为(x,y,z),另一个点的坐标为(x+1,y+1,z+1),这两个点之间的直线距离为
Figure BDA0001362412030000062
则判定这两个点为相邻关系。三维血管中心轴上点(x,y,z)与点(x+2,y,z)之间的直线距离为2,则判定这两个点不相邻。
例如,三维血管中心轴上四个点A、B、C、D之间的关系为:A<->B相邻,B<->C相邻,C<->D相邻,其它均不相邻,则建立的连接网络为点A、B、C、D构成的网络。
S3、读取用户选取的原点和末点坐标,计算得到三维血管中心轴上与原点距离最近的点以及与末点距离最近的点,采用迪杰斯特拉算法得到沿三维血管 中心轴从与原点距离最近的点到与末点距离最近的点的最短路径,其具体过程为:
1)分别计算三维血管中心轴上各个点与读取的原点和末点之间的距离,得到三维血管中心轴上与原点之间距离最近的点,该点记为起点S;得到三维血管中心轴上与末点之间距离最近的点,该点记为终点E。
2)根据步骤S2得到的连接网络,每条从起点S到终点E的可能路径为从起点S到终点E的相邻点,再到起点S的相邻点的相邻点,依次类推,直到相邻点中包含终点E,采用迪杰斯特拉算法得到沿三维血管中心轴从起点S到终点E的最短路径。例如,S->A->B->E,起点S与点A相邻,点A与点B相邻,点B与终点E相邻。
S4、对三维血管中心轴周围的像素进行处理,得到初步截取的血管,其具体过程为:
1)根据步骤S3获取的最短路径,得到最短路径上各个点的坐标。
2)对于二元三维矩阵的每个非0的血管像素点,定义该像素点与最短路径上各个点沿最短路径的距离为该像素点的坐标(x,y,z)同最短路径上各个点坐标的距离的最小值。
3)计算二元三维矩阵的每个非0的血管像素点与最短路径上各个点坐标之间的距离。
4)预设距离临界值,如果步骤3)计算得到的距离大于距离临界值,则将该距离对应的血管像素值设为背景值0,以此将与沿最短路径较远的血管点设为背景值,仅保留沿最短路径附近的血管,得到初步截取的血管。
S5、根据初步截取的血管,采用Marching Cube算法获得满足条件的临界点坐标,将所有临界点坐标相连得到血管表面,其具体过程为:
1)以起点S为种子,采用区域生长的方法,对与种子不相连但在步骤S4中保留的像素点进行清除,即将其像素值设为背景值0,得到新的二元三维矩阵。
该新的二元三维矩阵中,值为1的像素点同时满足步骤S4中与最短路径上 各个点坐标之间的距离小于距离临界值以及与种子相连这两个条件。该新的二元三维矩阵能够代表起点S到终点E之间的血管。
2)在新的二元三维矩阵的基础上,1和0的临界点定义为左右或上下或前后两个相邻点的像素分别为1和0的点。例如,坐标为(x,y,z)的A点对应像素值为1,坐标为(x,y+1,z)的B点对应像素值为0,A点与B点之间的临界点则为坐标为(x,y+0.5,z)的点。
S6、根据得到的血管表面,计算得到微导管折线路径。
为计算微导管路径,以起点S为起始点,以通过该起始点的沿最短路径方向为起始方向,以直线的方式延伸路径,方向记为
Figure BDA0001362412030000081
如果路径抵达血管壁,则记录下对应血管壁的坐标W,找到与坐标W距离最小的沿最短路径上的点C,求得方向
Figure BDA0001362412030000084
相对于点C的对称方向
Figure BDA0001362412030000082
并保证从坐标W开始,延方向
Figure BDA0001362412030000083
延伸会使得路径朝向血管内部,直到再次碰到血管壁。然后重复上述过程,直到路径与终点E之间的距离小于距离临界值,得到微导管折线路径。其中,依照完成截取的二元化三维矩阵,是否为血管壁以0和1的边界点为判断依据。
S7、对微导管折线路径进行平滑处理,用平滑曲线代替折线,得到微导管路径。
S8、根据微导管路径上的折线点,得到塑形器的折线路径,其具体过程为:
1)根据微导管路径上的折线点,计算出每段线段的长度及方向,其中方向由三维单元向量表示,并计算出由两段线段决定的坐标系到下两段线段决定的坐标系的旋转矩阵。利用旋转矩阵得到相邻线段组之间的角度。每个线段组中包含两条相邻的线段。
从两条线段到三维坐标系的计算方法为:假定v1,v2为不平行的两段相邻线段,取与v1同方向的单元向量为坐标系X轴,v1叉乘v2得到的方向为坐标系Y轴,Y轴单元向量和X轴单元向量叉乘得到Z轴单元向量。
2)假定微导管塑形器的角度与微导管路径的折线角度有固定的倍数关系α,根据折线点的坐标和步骤1)中计算得到的旋转矩阵,从第一条线段开始,将微导管路径旋转矩阵对应的旋转角度乘以给定的倍数α,得到塑形器下一条 线段相对于前一条线段的旋转矩阵,并依据该矩阵得到塑形器下一条线段的方向,同时保证塑形器下一条线段的长度同微导管路径对应线段的长度相同,由此确定下一条线段,直到得到塑形器的全部折线路径。其中,倍数α通常取值为2。
S9、对塑形器的折线路径进行平滑处理,得到表示微导管塑形器的二元三维矩阵。
S10、根据微导管塑形器的二元三维矩阵,利用Marching Cube的方法,选取像素值为0和2的表面作为微导管塑形器的外表面,像素值为1和2的表面作为微导管塑形器的内表面,得到微导管塑形器的表面文件,其具体过程为:
1)假定微导管的半径为r1,塑形器的半径为r2。
2)对每个背景像素(像素值为0),计算其坐标(x,y,z)与微导管路径上每个点的距离,对得到的距离值进行比较后得到最小距离值d。
3)如果d<=r1,则将背景像素的像素值设为1。如果r1<d<=r2,则将背景像素的像素值设为2。如果d>r2,则将背景像素的像素值设为0。
4)选取像素值为0和2的表面作为微导管塑形器的外表面,像素值为1和2的表面作为微导管塑形器的内表面,得到微导管塑形器的表面文件。
本发明微导管塑形器的生成方法还包括以下步骤,在生成的微导管路径的基础上,对微导管折线路径上的折线点的位置进行调整,得到调整后的微导管路径。
使用本发明微导管塑形器的生成方法生成微导管塑形器时,用户给定大概的微导管原点和末点,本发明微导管塑形器的生成方法自动计算出微导管的最优路径,并可依据原点和末点对血管束进行简化显示,以便于使用户更清晰的了解手术相关区域的血管构造。基于得到的微导管路径,用户可以通过对关键点的位置进行必要的调整,以便达到最优效果,在优化的微导管路径基础上,自动生成微导管塑形器的表面文件,根据生成的表面文件得到微导管塑形器。
如图2所示,本发明还提供了一种微导管塑形器的生成系统,其包括三维血管中心轴生成模块1、血管截取模块2、血管表面生成模块3、微导管路径生 成模块4、塑形器路径生成模块5和塑形器表面文件生成模块6。其中,三维血管中心轴生成模块1用于根据三维血管数据生成三维血管中心轴。血管截取模块2用于根据三维血管中心轴上各点之间的连接关系对三维血管中心轴周围的像素进行处理,得到截取的血管。血管表面生成模块3用于根据截取的血管生成血管表面文件。微导管路径生成模块4用于根据血管表面文件处理得到微导管路径。塑形器路径生成模块5用于根据微导管路径处理得到微导管塑形器的路径。塑形器表面文件生成模块6用于根据微导管塑形器的路径生成微导管塑形器的表面文件。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明所披露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (8)

1.一种微导管塑形器的生成方法,其特征在于,包括以下步骤:
根据三维血管图像得到三维血管中心轴;
根据三维血管中心轴上各点与其周围所有点的相邻关系,建立连接网络;
读取用户选取的原点和末点坐标,计算沿三维血管中心轴从与原点距离最近的点到与末点距离最近的点的最短路径;
对三维血管中心轴周围的像素进行处理,得到初步截取的血管;
根据初步截取的血管,获得临界点坐标,将所有临界点坐标相连得到血管表面;
根据得到的血管表面,计算得到微导管折线路径;
对微导管折线路径进行平滑处理,得到微导管路径;
根据微导管路径上的折线点,得到塑形器的折线路径;
对塑形器的折线路径进行平滑处理,得到表示微导管塑形器的二元三维矩阵;
根据微导管塑形器的二元三维矩阵,得到微导管塑形器的表面文件;
所述对三维血管中心轴周围的像素进行处理得到初步截取的血管的具体过程为:
根据获取的最短路径,得到最短路径上各个点的坐标;
对于二元三维矩阵的每个非0的血管像素点,定义该像素点与最短路径上各个点沿最短路径的距离为该像素点的坐标(x,y,z)同最短路径上各个点坐标的距离的最小值;
计算二元三维矩阵的每个非0的血管像素点与最短路径上各个点坐标之间的距离;
预设距离临界值,如果计算得到的距离大于距离临界值,则将该距离对应的血管像素值设为背景值0,以此将与沿最短路径较远的血管点设为背景值,仅保留沿最短路径附近的血管,得到初步截取的血管;
所述根据得到的血管表面得到微导管折线路径的具体过程为:
以通过起点的沿最短路径方向为起始方向,以直线的方式延伸路径,方向记为
Figure FDA0002282257940000021
如果路径抵达血管壁,则记录下对应血管壁的坐标W,找到与坐标W距离最小的沿最短路径上的点C,求得方向
Figure FDA0002282257940000022
相对于点C的对称方向
Figure FDA0002282257940000023
并保证从坐标W开始,延方向
Figure FDA0002282257940000024
延伸会使得路径朝向血管内部,直到再次碰到血管壁;
重复上述过程,直到路径与终点之间的距离小于距离临界值,得到微导管折线路径。
2.如权利要求1所述的一种微导管塑形器的生成方法,其特征在于,所述根据三维血管图像得到三维血管中心轴的步骤中,采用图像侵蚀算法对表示三维血管图像的二元三维矩阵重复进行对称的图像侵蚀处理,直到得到三维血管中心轴。
3.如权利要求1所述的一种微导管塑形器的生成方法,其特征在于,所述根据三维血管中心轴上各点与其周围所有点的相邻关系建立连接网络的步骤中,预设三维血管中心轴上两个点之间的最小直线距离为L,如果三维血管中心轴上两个点之间的距离小于或等于
Figure FDA0002282257940000025
则判定这两个点为相邻关系。
4.如权利要求1所述的一种微导管塑形器的生成方法,其特征在于,所述计算沿三维血管中心轴从与原点距离最近的点到与末点距离最近的点的最短路径的具体过程为:
分别计算三维血管中心轴上各个点与读取的原点和末点之间的距离,得到三维血管中心轴上与原点之间距离最近的点以及三维血管中心轴上与末点之间距离最近的点;其中,三维血管中心轴上与原点之间距离最近的点为起点,三维血管中心轴上与末点之间距离最近的点为终点,
根据得到的连接网络,采用迪杰斯特拉算法得到沿三维血管中心轴从起点到终点的最短路径。
5.如权利要求1或2或3或4所述的一种微导管塑形器的生成方法,其特征在于,所述得到血管表面的具体过程为:
以起点为种子,采用区域生长的方法,对与种子不相连但在得到初步截取的血管的过程中保留的像素点进行清除,得到新的二元三维矩阵;
在新的二元三维矩阵的基础上,获得满足条件的临界点坐标,将所有临界点坐标相连得到血管表面;其中,临界点坐标满足的条件包括与种子相连且与最短路径上各个点坐标之间的距离小于距离临界值。
6.如权利要求1或2或3或4所述的一种微导管塑形器的生成方法,其特征在于,所述根据微导管路径上的折线点得到塑形器的折线路径的具体过程为:
根据微导管路径上的折线点,计算出每段线段的长度及方向,其中方向由三维单元向量表示,并计算出由两段线段决定的坐标系到下两段线段决定的坐标系的旋转矩阵;
利用旋转矩阵得到相邻线段组之间的角度;每个线段组中包含两条相邻的线段;
假定微导管塑形器的角度与微导管路径的折线角度有固定的倍数关系a,根据折线点的坐标和计算得到的旋转矩阵,从第一条线段开始,将微导管路径旋转矩阵对应的旋转角度乘以给定的倍数a,得到塑形器下一条线段相对于前一条线段的旋转矩阵,并依据该矩阵得到塑形器下一条线段的方向,同时保证塑形器下一条线段的长度同微导管路径对应线段的长度相同,由此确定下一条线段,直到得到塑形器的全部折线路径。
7.如权利要求1或2或3或4所述的一种微导管塑形器的生成方法,其特征在于,所述得到微导管塑形器的表面文件的具体过程为:
假定微导管的半径为rl,塑形器的半径为r2;
对每个背景像素,计算其坐标(x,y,z)与微导管路径上每个点的距离,对得到的距离值进行比较后得到最小距离值d;
如果d<=rl,则将背景像素的像素值设为1;如果rl<d<=r2,则将背景像素的像素值设为2;如果d>r2,则将背景像素的像素值设为0;
选取像素值为0和2的表面作为微导管塑形器的外表面,像素值为1和2的表面作为微导管塑形器的内表面,得到微导管塑形器的表面文件。
8.如权利要求1-7任一项所述的一种微导管塑形器的生成系统,其特征在于,它包括三维血管中心轴生成模块、血管截取模块、血管表面生成模块、微导管路径生成模块、塑形器路径生成模块和塑形器表面文件生成模块;
所述三维血管中心轴生成模块用于根据三维血管数据生成三维血管中心轴;所述血管截取模块用于根据三维血管中心轴上各点之间的连接关系对三维血管中心轴周围的像素进行处理,得到截取的血管;所述血管表面生成模块用于根据截取的血管生成血管表面文件;所述微导管路径生成模块用于根据血管表面文件处理得到微导管路径;所述塑形器路径生成模块用于根据微导管路径处理得到微导管塑形器的路径;所述塑形器表面文件生成模块用于根据微导管塑形器的路径生成微导管塑形器的表面文件。
CN201710628619.0A 2017-07-27 2017-07-27 一种微导管塑形器的生成方法及生成系统 Active CN107392872B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710628619.0A CN107392872B (zh) 2017-07-27 2017-07-27 一种微导管塑形器的生成方法及生成系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710628619.0A CN107392872B (zh) 2017-07-27 2017-07-27 一种微导管塑形器的生成方法及生成系统

Publications (2)

Publication Number Publication Date
CN107392872A CN107392872A (zh) 2017-11-24
CN107392872B true CN107392872B (zh) 2020-04-28

Family

ID=60341830

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710628619.0A Active CN107392872B (zh) 2017-07-27 2017-07-27 一种微导管塑形器的生成方法及生成系统

Country Status (1)

Country Link
CN (1) CN107392872B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108478918A (zh) * 2018-06-01 2018-09-04 成都真实维度科技有限公司 一种动脉瘤介入手术微导管塑形辅助支具及其制作方法
CN109447967B (zh) * 2018-10-26 2022-03-01 强联智创(北京)科技有限公司 一种颅内动脉瘤图像的分割方法及系统
CN110288693A (zh) * 2019-06-06 2019-09-27 昆明医科大学第二附属医院 一种颅内动脉瘤介入手术智能微导管塑形系统
CN110507415B (zh) * 2019-08-13 2021-10-15 广州新脉科技有限公司 微导管虚拟介入的路径模拟方法及装置
CN114366296B (zh) * 2021-12-31 2023-05-30 杭州脉流科技有限公司 改进的微导管路径生成方法、塑形方法、设备和存储介质
CN114788919B (zh) * 2022-03-09 2023-12-08 北京柯莱文科技咨询有限公司 动脉瘤介入手术微导管塑形辅助模具及其制作方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2843373Y (zh) * 2005-11-24 2006-12-06 张绪平 介入用导管快速塑形器
CN102763135A (zh) * 2009-12-23 2012-10-31 通用电气公司 用于自动分割和时间跟踪的方法
CN102871729A (zh) * 2011-05-20 2013-01-16 西门子公司 使用医疗器械规划血管内操作用的装置和方法
CN204520841U (zh) * 2015-02-15 2015-08-05 李晓青 动脉瘤栓塞术中微导管塑形装置
CN105530887A (zh) * 2013-09-04 2016-04-27 皇家飞利浦有限公司 机器人系统
CN106037853A (zh) * 2016-05-04 2016-10-26 李泽福 一种用于颅内动脉瘤栓塞术的微导管塑形支架及制备方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7794424B2 (en) * 1991-07-15 2010-09-14 Paskar Larry D Catheter with out-of-plane configurations
US8460236B2 (en) * 2010-06-24 2013-06-11 Hansen Medical, Inc. Fiber optic instrument sensing system
US9057600B2 (en) * 2013-03-13 2015-06-16 Hansen Medical, Inc. Reducing incremental measurement sensor error

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2843373Y (zh) * 2005-11-24 2006-12-06 张绪平 介入用导管快速塑形器
CN102763135A (zh) * 2009-12-23 2012-10-31 通用电气公司 用于自动分割和时间跟踪的方法
CN102871729A (zh) * 2011-05-20 2013-01-16 西门子公司 使用医疗器械规划血管内操作用的装置和方法
CN105530887A (zh) * 2013-09-04 2016-04-27 皇家飞利浦有限公司 机器人系统
CN204520841U (zh) * 2015-02-15 2015-08-05 李晓青 动脉瘤栓塞术中微导管塑形装置
CN106037853A (zh) * 2016-05-04 2016-10-26 李泽福 一种用于颅内动脉瘤栓塞术的微导管塑形支架及制备方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"3D打印辅助微导管塑形在颅内动脉瘤栓塞术中应用";徐超 等;《介入放射学杂志》;20170131;第26卷(第1期);1-5 *
"Microcatheter Shaping for Intracranial Aneurysm Coiling Using the 3-Dimensional Printing Rapid Prototyping Technology: Preliminary Result in the First 10 Consecutive Cases";KatsunariNamba et al;《World Neurosurgery》;20150731;第84卷(第1期);178-186 *
"Tailor-made shaping of microcatheters using three-dimensional printed vessel models for endovascular coil embolization";Toshihiro Ishibashi et al;《Computers in Biology and Medicine》;20161031;第77卷(第1期);59-63 *

Also Published As

Publication number Publication date
CN107392872A (zh) 2017-11-24

Similar Documents

Publication Publication Date Title
CN107392872B (zh) 一种微导管塑形器的生成方法及生成系统
CN107468334B (zh) 一种三维微导管塑形辅助设计系统及设计方法
US9472017B2 (en) Fast rendering of curved reformation of a 3D tubular structure
WO2020083374A1 (en) Method and system for measuring morphological parameters of an intracranial aneurysm image
JP5129480B2 (ja) 管状臓器の3次元再構成を行うシステム及び血管撮像装置の作動方法
JP4891541B2 (ja) 血管狭窄率解析システム
CN105659289B (zh) 血管分割方法
US20150257850A1 (en) Image processing apparatus, image processing method and program
JP6230608B2 (ja) 物体の3dモデリング
CN105389815B (zh) 一种乳腺图像配准方法及装置
JP2002504385A (ja) 自動経路立案方法
US20150257655A1 (en) Medical image processing apparatus
CN105913479B (zh) 一种基于心脏ct图像的血管曲面重建方法
US11704819B2 (en) Apparatus and method for aligning 3-dimensional data
JP5194138B2 (ja) 画像診断支援装置およびその動作方法、並びに画像診断支援プログラム
US20230113035A1 (en) 3d pathfinder visualization
KR101001754B1 (ko) 초음파 시스템에서 초음파 영상의 프레임 보간 방법 및장치
CN111784751B (zh) 基于3d/2d配准的导丝3d模拟跟踪方法及装置
KR102347029B1 (ko) 2d xa 영상과 3d cta 영상 간의 강체 정합 방법, 이를 수행하기 위한 기록 매체 및 장치
Mohamed et al. Computer-aided planning for endovascular treatment of intracranial aneurysms (CAPETA)
JP2019150358A (ja) 画像処理装置と画像処理プログラムと画像処理方法
CN111243084B (zh) 一种用于心脏放射治疗中的心脏三维模型构建方法
JP7155670B2 (ja) 医用画像処理装置、医用画像処理方法、プログラム、及びデータ作成方法
Sgouritsa et al. Neck localization and geometry quantification of intracranial aneurysms
JP6937187B2 (ja) 乳房形状変化予測方法、乳房形状変化予測システム、及び乳房形状変化予測プログラム

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
CP02 Change in the address of a patent holder

Address after: 100176 room 901, building 3, yard 2, Ronghua South Road, Daxing Economic and Technological Development Zone, Beijing

Patentee after: QIANGLIAN ZHICHUANG (BEIJING) TECHNOLOGY Co.,Ltd.

Address before: 207 room 3, building 538, Yongfeng Tun, Beijing, Haidian District, 100000

Patentee before: QIANGLIAN ZHICHUANG (BEIJING) TECHNOLOGY Co.,Ltd.

CP02 Change in the address of a patent holder