CN109727307B - 一种表面网格切割方法、装置及计算机可读存储介质 - Google Patents

一种表面网格切割方法、装置及计算机可读存储介质 Download PDF

Info

Publication number
CN109727307B
CN109727307B CN201811637377.2A CN201811637377A CN109727307B CN 109727307 B CN109727307 B CN 109727307B CN 201811637377 A CN201811637377 A CN 201811637377A CN 109727307 B CN109727307 B CN 109727307B
Authority
CN
China
Prior art keywords
cutting
plane
point
target tissue
acquiring
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
CN201811637377.2A
Other languages
English (en)
Other versions
CN109727307A (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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare 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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201811637377.2A priority Critical patent/CN109727307B/zh
Publication of CN109727307A publication Critical patent/CN109727307A/zh
Application granted granted Critical
Publication of CN109727307B publication Critical patent/CN109727307B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明实施例提供了一种表面网格切割方法、装置及计算机可读存储介质,涉及医疗技术领域。一方面,本发明实施例根据医学图像,获取所述医学图像中的目标组织的树状结构;从而,基于所述树状结构,构建切割模型;进而,根据所述切割模型,切割所述目标组织的表面网格。因此,本发明实施例提供的技术方案能够解决现有技术切割表面网格花费时间长,切割表面网格效率低的问题。

Description

一种表面网格切割方法、装置及计算机可读存储介质
【技术领域】
本发明涉及医疗技术领域,尤其涉及一种表面网格切割方法、装置及计算机可读存储介质。
【背景技术】
虚拟手术仿真系统是虚拟现实技术在医学领域的一个典型应用,虚拟手术仿真系统建模能力强大,利用各种医学影像数据,采用虚拟现实技术建立一个模拟环境,能够对手术情况和场景进行仿真,模拟病人接收手术的手术效果。虚拟手术是医生进行真实手术前准备的一项重要内容。
在仿真系统中,使用表面网格近似表示物体,也就是说,在虚拟手术仿真系统中用表面网格模拟人体组织。在真实手术中需要对病变组织进行切割和缝合,因此,在仿真系统中对表面网格切割是必不可少,而如何快速地对表面网格进行切割一直是研究的热点与难点。
目前对于切割表面网格的研究相对比较少,而更多的基于可视化界面操作,人为设置切割位置,然后对表面网格进行切割。这种切割方法应用于结构相对复杂的表面网格时,需要花费大量的时间在前期表面网格边界切割处理过程中,切割表面网格效率低。
【发明内容】
有鉴于此,本发明实施例提供了一种表面网格切割方法、装置及计算机可读存储介质,用以解决现有技术切割表面网格花费时间长,切割表面网格效率低的问题。
第一方面,本发明实施例提供了一种表面网格切割方法,包括:
根据医学图像,获取所述医学图像中的目标组织的树状结构;
基于所述树状结构,构建切割模型;
根据所述切割模型,切割所述目标组织的表面网格。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述根据医学图像,获取所述医学图像中的目标组织的树状结构,包括:
根据所述切割模型,切割所述目标组织的表面网格。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述根据医学图像,获取所述医学图像中的目标组织的树状结构,包括:
根据所述医学图像,获取所述目标组织的中心线;
获取所述中心线上每个点的位置参数;
基于每个点的位置参数,构建所述树状结构。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,每个点的位置参数包括该点位置坐标以及该点的切线向量。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述切割模型包括三维模型与二维模型中的至少一种;
其中,所述三维模型包括盒状模型;
其中,所述二维模型包括有界平面和无界平面中至少一种。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述方法还包括:
获取所述目标组织的输出控制参数;
根据所述输出控制参数,输出切割后的表面网格。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述目标组织为血管组织,切割模型为盒状物体,所述根据所述切割模型,切割所述目标组织的表面网格包括:
根据血管管径扩大因子和高度扩大因子,设置所述表面网格与所述盒状物体相交在第一平面;
获取表面网格与所述第一平面的所有交点;
根据所有交点,构建若干个第一多边形,并为每个第一多边形设置边界属性值;
获取表面网格上与第一平面相交的所有三角面片,以及,所有交点,并将表面网格上与第一平面相交的三角面片上的部分点以及所有交点,基于最大化最小角的条件,构建第二多边形。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述目标组织为血管组织,切割模型为无界平面时,所述根据所述切割模型,切割所述目标组织的表面网格包括:
基于表面网格连通性,获取第一三角面片,其中,第一三角面片为距离所述平面最近的三角面片;
获取所述第一三角面片和所述平面的交点;
根据所述交点,构建多边形或者三角化处理后的多边形,并设置边界属性值;
获取表面网格与所述平面相交的三角形以及相交点,并将表面网格与所述平面相交的三角形以及相交点,基于最大化最小角的条件,构建第二三角面片。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述目标组织为血管组织,切割模型为有界平面时,所述根据所述切割模型,切割所述目标组织的表面网格包括:
获取所有与该平面相交的三角面片为第三三角面片;
根据表面网格的连通性,获取多个三角面片集合,其中,所述三角面片集合包括多个相连接的第三三角面片;
计算每个集合到平面设置点的距离,基于距离最小对应的三角面片集合构建第四三角面片;
获取所述第四三角面片和所述平面的交点;
根据所述交点,构建多个多边形或者将三角化处理的多边形,并设置边界属性值;
获取表面网格与所述平面相交的三角形以及相交点,并将表面网格与所述平面相交的三角形部分点以及相交点,基于最大化最小角的条件,构建第五三角面片。
第二方面,本发明实施例提供了一种表面网格切割装置,包括处理器、存储器、输入输出器和显示器,所述输入输出器用于接收或发送数据,所述显示器用于显示数据;所述存储器用于存放程序;所述处理器用于执行所述存储器存储的所述程序,以执行第一方面中任意一项所述的方法
第三方面,本发明实施例提供了一种计算机可读存储介质,包括计算机可读指令,当计算机读取并执行所述计算机可读指令时,使得计算机执行第一方面中任意一项所述的方法。
本发明实施例提供的技术方案,利用目标组织的树状结构,自动构建用于切割表面网格的切割模型,不需要操作者基于可视化界面操作,人为设置切割位置,从而降低对操作者的依赖,提高表面网格切割的灵活性,并且采用本发明提供的自动构建切割表面网格的切割模型,切割表面网格的方法,切割结构相对复杂的表面网格,也不需要花费大量的时间在前期表面网格边界切割处理过程中,进而提高了切割表面网格效率。
【附图说明】
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例所提供的一种表面网格切割方法的流程示意图;
图2是本发明实施例所提供的目标组织的树状结构示意图;
图3是本发明实施例针对图1所述的方法中步骤102的实现方法的流程示意图;
图4是本发明实施例所提供的另一种表面网格切割方法的流程示意图;
图5是本发明实施例所提供的采用盒状物体切割,切割边界为多边形的示意图;
图6是本发明实施例所提供的采用盒状物体切割,切割边界为三角化的示意图;
图7是本发明实施例所提供的无界平面切割冠脉及部分主动脉表面网格,切口封闭成三角形的效果展示图;
图8是本发明实施例所提供的采用无界平面切割冠脉及部分主动脉表面网格,切口不封闭的效果展示图;
图9是本发明实施例所提供的采用多个有界平面切割冠脉及部分主动脉表面网格,切口边界为多边形的效果展示图;
图10是本发明实施例所提供的采用多个有界平面切割冠脉及部分主动脉表面网格,切口边界为三角化的效果展示图;
图11是本发明实施例所提供的一种表面网格切割装置的结构示意图。
【具体实施方式】
为了更好的理解本发明的技术方案,下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。
应当理解,本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,一般表示前后关联对象是一种“或”的关系。
应当理解,尽管在本发明实施例中可能采用术语第一、第二等来描述多边形,但这些多边形不应限于这些术语。这些术语仅用来将多边形彼此区分开。例如,在不脱离本发明实施例范围的情况下,第一多边形也可以被称为第二多边形,类似地,第二多边形也可以被称为第一多边形。
取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。类似地,取决于语境,短语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测(陈述的条件或事件)”。
本发明实施例给出一种表面网格切割方法,请参考图1,其为本发明实施例所提供的方法的流程示意图,如图所示,该方法包括以下步骤:
102,根据医学图像,获取医学图像中的目标组织的树状结构。
其中,医学图像可以是DICOM(Digital Imaging and Communications inMedicine,医学数字成像和通信)图像。
目标组织可以包括但不限于血管组织、肺气管组织等。
104,基于树状结构,构建切割模型。
其中,切割模型包括三维模型与二维模型中的至少一种。三维模型可以是但不限于盒状模型(如,六面体模型、圆柱体模型等);二维模型可以是有界平面和无界平面中至少一种。
请参考图2,其为目标组织的树状结构示意图。具体的,执行主体通过遍历该树状结构的叶子节点,得到叶子节点末端的端点以及指定点,其中,该指定点为在远离该叶子节点靠近该叶子节点上一级节点所在方向上,距离该叶子节点末端的端点一定距离范围内的点,根据得到的叶子节点末端的端点和指定点,构建切割模型。
106,根据切割模型,切割目标组织的表面网格。
本发明实施例提供的技术方案,利用目标组织的树状结构,自动构建用于切割表面网格的切割模型,不需要操作者基于可视化界面操作,人为设置切割位置,降低了对操作者的依赖,提高表面网格切割的灵活性,并且采用本发明提供的自动构建切割表面网格的切割模型,切割表面网格的方法,切割结构相对复杂的表面网格,也不需要花费大量的时间在前期表面网格边界切割处理过程中,进而提高了切割表面网格效率。
进一步地,结合前述流程,本发明针对于步骤102根据医学图像,获取所述医学图像中的目标组织的树状结构的实现,提供了一种实现方式,其流程图如图3所述,步骤102具体包括:
1021,根据医学图像,获取目标组织的中心线。
一张医学图像中可能包括多种组织,为了能够快速地识别目标组织,获取目标组织的中心线,在获取目标组织的中心线之前,可以先将医学图像进行图像分割,得到图像分割后的图像,根据图像分割后的图像,获取目标组织的中心线。
在一个可行的实现方式中,获取目标组织的中心线可以利用骨架算法提取,具体为:首先,将读取到的医学图像DICOM进行图像分割,得到分割后的图像,将分割后的图像转换成二值图像,然后,利用骨架算法去除分割后的图像中目标组织的边界像素,从而获取到目标组织的中心线。
1022,获取中心线上每个点的位置参数。
其中,每个点的位置参数包括但不限于:该点位置坐标以及该点的切线向量。具体的,在根据医学图像,获取到目标组织的中心线后,执行主体可以为为该中心线,构建一个空间坐标系,从而得到该中心线上每个点的位置坐标。
1023,基于每个点的位置参数,构建树状结构。
其中,构建的树状结构主要是用于存储中心线上每个点的位置参数,即用该树状结构存储中心线上每个点的位置坐标和每个点的切线向量。
补充说明的是,在构建树状结构之前,为了使得到的每个点的位置参数准确,在获取到目标组织的中心线后,可以先对该中心线进行优化,去除该中心上偏差较大的数据,得到优化后的中心线,从而基于优化后的中心线,执行步骤1022和步骤1023,从而构建树状结构。
结合上述方法流程,本发明进一步提供了另一种实现方法,如图4所述,该方法还包括:
108,获取目标组织的输出控制参数。
其中,输出控制参数主要是用于控制输出切割后的表面网格。一个原始表面网格经一次切割后至少得到两个切割后的表面网格,执行主体基于输出控制参数,输出相应的切割后的表面网格。
需要说明的是,输出控制参数主要控制执行主体输出切割后的那些表面网格,其对表面网格的切割过程不起主要作用,因此,步骤108可以在步骤106之前执行;或者,步骤108也可以106之后执行;又或者,步骤108和步骤106还可以同时执行,本发明对于步骤108和步骤106执行的先后顺序不做限定。
110,根据所述输出控制参数,输出切割后的表面网格。
在一个具体应用场景中,切割模型为盒状物体,输出控制参数包括0和1,其中,输出控制参数为0,指示执行主体输出盒状物体以外的表面网格,输出控制参数为1,指示执行主体输出盒状物体以内的表面网格。
在另一个具体应用场景中,切割模型为有界平面(或无界平面),输出控制参数包括0和1,其中,输出控制参数为0,指示执行主体输出有界平面(或无界平面)上方的表面网格,输出控制参数为1,指示执行主体输出有界平面(或无界平面)下方的表面网格。
需要说明的是,在本发明中通过设置不同值的输出控制参数,实现控制输出显示给用户的不同的切割后的表面网格,因此,用户可以通过输入不同值的输出控制参数,可以使执行主体有选择性地输出切割后想要的表面网格,提高了用户体验性。
结合前述方法流程中,步骤102的实现方式,在一个可行的实施方式中,切割模型为六面体盒状物体,目标组织为血管,基于树状结构,构建该六面体盒状物体(下面称为盒状物体)的过程为:
(1),从血管的树状结构上,选择叶子节点的末端端点C以及指定一点A,,其中,A点一般设为该树状结构末端与末端端点C距离2到5点时对应的点;
(2),选择点A与末端点C之间的某点,例如点B,点B作为确定构建盒状物体沿着中心线轴向方向的向量一个标记点,点A作为另外一个标记点。点A与点B构成向量
(3),用点A与向量的反向向量,利用点法式确定盒状物体的第一个平面;
(4),选定一个血管高度扩大因子rate,将血管末端点C沿着向量往外延伸到点K,/>点法式确定盒状物体的第二个平面;
(5),用构建盒状物体的第一个平面与表面网格求得一个交点D。这个交点D与点A构成一个向量并求得/>方向上的单位/>
(6),选定一个血管管径扩大因子ratio,沿着向量方向延伸,获得一点坐标F,其中,A,E,F三点共线;
(7),选择点F以及向量点法式确定盒状物体的第三个平面;
(8),选择点F关于点A的对称点G,以及向量点法式确定盒状物体的第四个平面;
(9),用向量与向量/>求这两个向量的外积,求得一个向量/>
(10),用血管管径扩大因子ratio以及向量构建沿着向量/>方向上延伸,得到一个点坐标I;
(11),用点I坐标以及向量点法式确定盒状物体的第五个平面;
(12),求点I关于点A的对称点J。用点J以及向量点法式确定盒状物体的第六个平面。
其中,需要说明的是,经过上述步骤(1)~(12),构成的六个平面构成一个封闭的盒状物体。
下面结合两个具体应用场景,详细说明根据切割模型,切割目标组织的表面网格的具体实现过程:
示例一
假设切割模型为盒状物体(如,前述六面体盒状物体),目标组织为冠状动脉组织,其基于构建的六面体盒状物体切割表面网格的过程为:
步骤1,根据血管管径扩大因子和高度扩大因子,设置表面网格与盒状物体相交在第一平面。
需要说明的是,切割表面网格用的盒状物体可以限定切割范围,其具体为,通过设定血管径扩大因子ratio和高度扩大因子rate,使得表面网格与该盒状物体相交于一个平面上,而该盒状物体的其他几个平面起到限定范围的作用。
步骤2,获取表面网格与第一平面的所有交点。
步骤3,根据所有交点,构建若干个第一多边形,并为每个第一多边形设置边界属性值。
其中,构建第一多边形可以将所有交点首尾相连,从而构建出第一多边形。
需要说明的是,在构建第一多边形后,为每个多边形设置不同的边界属性值,这样可以使血管不同的末端有不同的边界属性值。具体的,请参考图5,其为使用盒状物体切割,切割边界为多边形的示意图。
另外,补充说明的是,执行主体还可以将构建的第一多边形进行三角化处理,使每个被三角化的切面上的三角形都有一个共同的边界属性值。具体的,请参考图6,其为使用盒状物体切割,切割边界为三角化的示意图。
值得说明的是,为每个第一多边形设置边界属性值,可以使得切割后得到的表面网格为带边界属性的表面网格,从而可以为体网格划分提供便利。
步骤4,获取表面网格上与第一平面相交的所有三角面片,以及,所有交点,并将表面网格上与第一平面相交的三角面片上的部分点以及所有交点,基于最大化最小角的条件,构建第二多边形,所述第二多边形可以为血管壁面上新的三角形。
示例二
假设切割模型为无界平面(如,前述六面体盒状物体),目标组织为冠状动脉组织,其基于构建的平面切割表面网格的过程为:
步骤1,基于表面网格连通性,获取第一三角面片。
其中,一个平面与多分支血管表面网格相交时,会与多个三角面片相交,且所有相连的三角面片归为一个三角面片集合里面,因此,产生多个三角面片集合,执行主体将距离该平面设置点最近的三角面片集取出来,得到第一三角面片,所述设置点为采用点法式定义平面时对应的设置点。
步骤2,获取第一三角面片和平面的交点。
步骤3,根据所述交点,构建多边形或者三角化处理后的多边形,并设置边界属性值。
将第一三角面片和平面的交点首尾相连,构成多边形,并赋予边界属性值;或者将构成的多边形进行三角化,形成新的三角面片,并对处于原先同一个多边形内的三角面片赋予同一边界属性值。
步骤4,获取表面网格与所述平面相交的三角形以及交点,并将表面网格与平面相交的三角形以及交点,基于最大化最小角的条件,构建第二三角面片。
实施例三
假设切割模型为有界平面(平面上的有限区域),目标组织为冠状动脉组织,其基于构建的平面切割表面网格的过程为:
步骤1,获取所有与该平面相交的三角面片为第三三角面片;
步骤2,根据表面网格的连通性,获取多个三角面片集合,其中,所述三角面片集合包括多个相连接的第三三角面片;
步骤3,计算每个集合到平面设置点的距离,基于距离最小对应的三角面片集合构建第四三角面片;
步骤4,获取所述第四三角面片和所述平面的交点;
步骤5,根据所述交点,构建多个多边形或者将三角化处理的多边形,并设置边界属性值;
步骤6,获取表面网格与所述平面相交的三角形以及相交点,并将表面网格与所述平面相交的三角形部分点以及相交点,基于最大化最小角的条件,构建第五三角面片。
另外,为了便于本领域技术人员直观地了解本发明技术方案效果,请参考图7至图10,其,图7是采用无界平面切割冠脉及部分主动脉表面网格,切口封闭成三角形的效果展示图;图8是采用无界平面切割冠脉及部分主动脉表面网格,切口不封闭的效果展示图;图9是采用有界平面切割冠脉及部分主动脉表面网格,切口边界为多边形的效果展示图;图10是采用有界平面切割冠脉及部分主动脉表面网格,切口边界为三角化的效果展示图。
需要说明的是,S102~S110的执行主体可以为一种表面网格切割装置,该装置可以位于本地表面网格切割设备的应用,或者还可以为位于本地表面网格切割设备的应用中的插件或软件开发工具包(Software Development Kit,SDK)等功能单元,本发明实施例对此不进行特别限定。
可以理解的是,所述应用可以是安装在终端上的应用程序(nativeApp),或者还可以是终端上的浏览器的一个网页程序(webApp),本发明实施例对此不进行限定。本发明实施例进一步给出实现上述方法实施例中各步骤及方法的装置实施例。
请参考图11,其为本发明实施例所提供的装置的功能方块图。如图所示,该装置包括:包括处理器10、存储器20、输入输出器30和显示器40,输入输出器30用于接收或发送数据,显示器40用于显示数据;存储器20用于存放程序;处理器10用于执行存储器20存储的程序,以执行上述任意一项表面网格切割方法。
由于本实施例中的处理器10能够执行前述表面网格切割方法,本实施例未详细描述的部分,可参考对表面网格切割方法的相关说明。
本发明实施例提供了一种计算机可读存储介质,包括计算机可读指令,当计算机读取并执行所述计算机可读指令时,使得计算机执行上述任意一项表面网格切割方法。
本发明实施例提供的技术方案,利用目标组织的树状结构,自动构建用于切割表面网格的切割模型,不需要操作者基于可视化界面操作,人为设置切割位置,从而降低对操作者的依赖,提高表面网格切割的灵活性,并且采用本发明提供的自动构建切割表面网格的切割模型,切割表面网格的方法,切割结构相对复杂的表面网格,也不需要花费大量的时间在前期表面网格边界切割处理过程中,进而提高了切割表面网格效率。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本发明所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元,可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中,包括若干指令用以使得一台计算机装置(可以是个人计算机,服务器,或者网络装置等)或处理器(Processor)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。

Claims (10)

1.一种表面网络切割方法,其特征在于,所述方法包括:
根据医学图像中目标组织的中心线,获取所述目标组织的树状结构;
通过遍历所述树状结构的叶子节点,得到叶子节点末端的端点和指定点;所述指定点为远离所述叶子节点靠近所述叶子节点上一级节点所在方向上,距离所述叶子节点末端的端点一定距离范围内的点;
根据所述叶子节点末端的端点和所述指定点,构建切割模型;
根据所述切割模型,切割所述目标组织的表面网格。
2.根据权利要求1所述的方法,其特征在于,所述根据医学图像中目标组织的中心线,获取所述目标组织的树状结构,包括:
根据所述医学图像,获取所述目标组织的中心线;
获取所述中心线上每个点的位置参数;
基于每个点的位置参数,构建所述树状结构。
3.根据权利要求2所述的方法,其特征在于,每个点的位置参数包括该点位置坐标以及该点的切线向量。
4.根据权利要求1所述的方法,其特征在于,所述切割模型包括三维模型与二维模型中的至少一种;
其中,所述三维模型包括盒状模型;
其中,所述二维模型包括有界平面和无界平面中至少一种。
5.根据权利要求1所述的方法,其特征在于,所述方法还包括:
获取所述目标组织的输出控制参数;
根据所述输出控制参数,输出切割后的表面网格。
6.根据权利要求1所述的方法,其特征在于,所述目标组织为血管组织,切割模型为盒状物体,所述根据所述切割模型,切割所述目标组织的表面网格包括:
根据血管管径扩大因子和高度扩大因子,设置所述表面网格与所述盒状物体相交在第一平面;
获取表面网格与所述第一平面的所有交点;
根据所有交点,构建若干个第一多边形,并为每个第一多边形设置边界属性值;
获取表面网格上与第一平面相交的所有三角面片,以及,所有交点,并将表面网格上与第一平面相交的三角面片上的部分点以及所有交点,基于最大化最小角的条件,构建第二多边形。
7.根据权利要求1所述的方法,其特征在于,所述目标组织为血管组织,切割模型为无界平面时,所述根据所述切割模型,切割所述目标组织的表面网格包括:
基于表面网格连通性,获取第一三角面片,其中,第一三角面片为距离所述平面最近的三角面片;
获取所述第一三角面片和所述平面的交点;
根据所述交点,构建多边形或者三角化处理后的多边形,并设置边界属性值;
获取表面网格与所述平面相交的三角形以及相交点,并将表面网格与所述平面相交的三角形以及相交点,基于最大化最小角的条件,构建第二三角面片。
8.根据权利要求1所述的方法,其特征在于,所述目标组织为血管组织,切割模型为有界平面时,所述根据所述切割模型,切割所述目标组织的表面网格包括:
获取所有与该平面相交的三角面片为第三三角面片;
根据表面网格的连通性,获取多个三角面片集合,其中,所述三角面片集合包括多个相连接的第三三角面片;
计算每个集合到平面设置点的距离,基于距离最小对应的三角面片集合构建第四三角面片;
获取所述第四三角面片和所述平面的交点;
根据所述交点,构建多个多边形或者将三角化处理的多边形,并设置边界属性值;
获取表面网格与所述平面相交的三角形以及相交点,并将表面网格与所述平面相交的三角形部分点以及相交点,基于最大化最小角的条件,构建第五三角面片。
9.一种表面网格切割装置,其特征在于,包括处理器、存储器、输入输出器和显示器,所述输入输出器用于接收或发送数据,所述显示器用于显示数据;所述存储器用于存放程序;所述处理器用于执行所述存储器存储的所述程序,以执行权利要求1-8任意一项所述的方法。
10.一种计算机可读存储介质,其特征在于,包括计算机可读指令,当计算机读取并执行所述计算机可读指令时,使得计算机执行如权利要求1-8任意一项所述的方法。
CN201811637377.2A 2018-12-29 2018-12-29 一种表面网格切割方法、装置及计算机可读存储介质 Active CN109727307B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811637377.2A CN109727307B (zh) 2018-12-29 2018-12-29 一种表面网格切割方法、装置及计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811637377.2A CN109727307B (zh) 2018-12-29 2018-12-29 一种表面网格切割方法、装置及计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN109727307A CN109727307A (zh) 2019-05-07
CN109727307B true CN109727307B (zh) 2023-08-22

Family

ID=66297951

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811637377.2A Active CN109727307B (zh) 2018-12-29 2018-12-29 一种表面网格切割方法、装置及计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN109727307B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114418836B (zh) * 2022-03-11 2022-06-21 腾讯科技(深圳)有限公司 数据处理方法、装置、设备以及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6608628B1 (en) * 1998-11-06 2003-08-19 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration (Nasa) Method and apparatus for virtual interactive medical imaging by multiple remotely-located users
CN101625766A (zh) * 2009-08-03 2010-01-13 深圳先进技术研究院 医学图像处理方法
CN103345568A (zh) * 2013-05-28 2013-10-09 深圳市旭东数字医学影像技术有限公司 基于三维模型的手术规划的方法及其系统
CN104318056A (zh) * 2014-09-24 2015-01-28 北京航空航天大学 基于位置动力学的软组织变形和切割模拟方法
CN104715475A (zh) * 2015-01-22 2015-06-17 中南大学 一种基于调和场的牙颌三维模型自动分割全部牙冠的方法
CN105184086A (zh) * 2014-06-16 2015-12-23 西门子公司 用于冠状动脉中改进的血液动力学计算的方法及系统
CN108537886A (zh) * 2018-04-09 2018-09-14 哈尔滨理工大学 一种虚拟手术切割中的高质量网格划分和优化方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7040896B2 (en) * 2000-08-16 2006-05-09 Align Technology, Inc. Systems and methods for removing gingiva from computer tooth models
JP2007537815A (ja) * 2004-05-18 2007-12-27 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 三次元変形メッシュモデルを用いてオブジェクトの三次元ツリー状管状表面を自動セグメント化するための画像処理システム
JP4450786B2 (ja) * 2005-11-15 2010-04-14 ザイオソフト株式会社 画像処理方法および画像処理プログラム
US20110071799A1 (en) * 2009-09-21 2011-03-24 Per Arne Slotte Grid models
US9594976B2 (en) * 2014-02-18 2017-03-14 Siemens Healthcare Gmbh Sparse appearance learning-based segmentation
US9824492B2 (en) * 2016-03-24 2017-11-21 Vital Images, Inc. Hollow object model visualization in medical images
CN106875462B (zh) * 2017-01-13 2020-07-07 北京航空航天大学 一种基于元球模型和混合驱动方法的实时数字器官切割方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6608628B1 (en) * 1998-11-06 2003-08-19 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration (Nasa) Method and apparatus for virtual interactive medical imaging by multiple remotely-located users
CN101625766A (zh) * 2009-08-03 2010-01-13 深圳先进技术研究院 医学图像处理方法
CN103345568A (zh) * 2013-05-28 2013-10-09 深圳市旭东数字医学影像技术有限公司 基于三维模型的手术规划的方法及其系统
CN105184086A (zh) * 2014-06-16 2015-12-23 西门子公司 用于冠状动脉中改进的血液动力学计算的方法及系统
CN104318056A (zh) * 2014-09-24 2015-01-28 北京航空航天大学 基于位置动力学的软组织变形和切割模拟方法
CN104715475A (zh) * 2015-01-22 2015-06-17 中南大学 一种基于调和场的牙颌三维模型自动分割全部牙冠的方法
CN108537886A (zh) * 2018-04-09 2018-09-14 哈尔滨理工大学 一种虚拟手术切割中的高质量网格划分和优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种基于标志点的三维表面模型的切割方法;刘奇等;《计算机应用研究》;20070210(第02期);全文 *

Also Published As

Publication number Publication date
CN109727307A (zh) 2019-05-07

Similar Documents

Publication Publication Date Title
KR102016959B1 (ko) 혈관 모델 생성 방법 및 장치
CN109285225B (zh) 一种基于医学影像的虚拟现实辅助手术的建立方法
US10262447B2 (en) Systems and methods for virtual entity animation
EP2688040A2 (en) 3D registration of tomographic models using particularly chosen points
US9984311B2 (en) Method and system for image segmentation using a directed graph
CN111161331B (zh) 一种bim模型和gis模型配准方法
CN110084894B (zh) 三维模型的局部放大展示方法、装置与电子设备
Shabat et al. Design of porous micro-structures using curvature analysis for additive-manufacturing
CN113160189A (zh) 血管中心线提取方法、装置、设备及存储介质
CN109300147A (zh) 血管二维中心线和三维中心线配准的方法及装置
CN114998490B (zh) 一种虚拟对象生成方法、装置、设备及存储介质
CN113724401A (zh) 一种三维模型切割方法、装置、计算机设备和存储介质
CN115222879B (zh) 一种模型减面处理方法、装置、电子设备及存储介质
CN109727307B (zh) 一种表面网格切割方法、装置及计算机可读存储介质
CN112614234A (zh) 对混合现实三维场景进行编辑的方法及混合现实设备
WO2021030995A1 (zh) 基于vrds ai下腔静脉影像的分析方法及产品
JP6569389B2 (ja) ポリゴンモデル構造化装置および立体物造形用データ削減装置
CN113506365A (zh) 图像显示方法及相关装置和电子设备、存储介质
CN113101666A (zh) 游戏角色模型方法、装置、计算机设备及存储介质
JP6565494B2 (ja) 立体物造形用データ削減装置
CN115619986B (zh) 场景漫游方法、装置、设备和介质
US9135697B2 (en) Method and system for determining a boundary surface network
CN108428230B (zh) 三维虚拟器官中处理曲面的方法、装置、存储介质及设备
CN108198213B (zh) 一种标记管状结构的方法及装置、管状结构的处理终端
CN112381825B (zh) 用于病灶区图像几何特征提取的方法和相关产品

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
CB02 Change of applicant information

Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant