CN109636915A - 一种基于边增长的医学图像三维重建等值面网格提取方法 - Google Patents
一种基于边增长的医学图像三维重建等值面网格提取方法 Download PDFInfo
- Publication number
- CN109636915A CN109636915A CN201811560982.4A CN201811560982A CN109636915A CN 109636915 A CN109636915 A CN 109636915A CN 201811560982 A CN201811560982 A CN 201811560982A CN 109636915 A CN109636915 A CN 109636915A
- Authority
- CN
- China
- Prior art keywords
- point
- cube
- contour surface
- intersection
- queue
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Image Generation (AREA)
Abstract
本发明属于医学图像技术领域,提供一种医学图像三维重建等值面网格提取方法。本发明,先找到种子三角形,再以三角形的边进行增长。增长过程分为等值面与立方体的交点求解和交点组成三角形网格两个过程。整个重建过程拥有一个队列,存储增长边的两个端点坐标和对应相邻要增长立方体的坐标,除种子三角形外,其它三角形网格都是通过队列中取出数据生成的,当队列为空时,生成指定阈值的等值面三角形网格。
Description
技术领域
本发明属于医学图像领域,具体涉及一种医学图像三维重建等值面网格提取方法。
背景技术
随着医学影像技术的发展,医学影像已成为医生判断病症的重要手段。利用医学的二维断层图像重建出人体组织的三维图像并显示,极大提高了医生对病情的诊断率,且减小了医院对胶片的开销。从医学二维断层图像中提取三维模型的过程称为医学图像的三维重建,三维模型一般用网格表示,这个过程也称为等值面网格提取。当前,医学图像网格生成技术最为成熟的是MC算法,及其MC的衍生算法,因其原理简单,实现容易,得到了广泛的运用。
MC算法通过遍历所有的立方体生成网格,其中大部分时间浪费在空体素(和等值面没有相交的立方体)计算上,且等值面是在各个立方体中单独生成的,没有考虑各个立方体之间的联系,生成的等值面存在一定的二义性。本发明针对这些问题,提出一种新的算法,消除二义性问题,通过各个立方体中的等值面之间的联系进行重建。
发明内容
本发明的目的在于提供一种三维重建的网格提取方法,通过各个立方体中的等值面之间的联系进行重建,消除二义性,形成一个连续的等值面网格。
本发明采取的技术方案分为以下步骤:
A、读取医学图像数据,构造立方体,选取种子三角形;
B、以前景边进行网格增长,求解以前景边增长时等值面与立方体边的交点;
C、将所有交点按求解的顺序3个点一起形成三角形,并将三角形放入三角链表中;
D、判断所有两个相邻交点组成的边是否为增长边,如果是增长边,把增长边的信息放入队列中;
E、从队列中取出数据进行增长,直到队列为空。
根据本发明的一种具体实施方式,在步骤A中,读取医学图像数据后,通过相邻两层医学断层图像各4个相邻像素点构成立方体。遍历中间一层的所有立方体,将立方体8个顶点的灰度值分别与等值面的灰度值(等值面的灰度值简称阈值)进行比较,并将顶点分类为大于阈值和不大于阈值两部分,只有当两部分中顶点数量较小的那个部分的点的数量刚好为1时,才认为该立方体中有种子三角形,并称数量为1时对应的那个点为基本点。以基本点为基准,分别插值求出和基本点相连三条立方体的边与等值面的交点(一条边若和等值面相交,说明边上两个顶点的灰度值,一个大于阈值,一个不大于阈值,假设c为阈值,v1、v2、c1、c2分别为两个顶点的空间坐标和灰度值,则等值面与边的交点坐标为然后将所求的三个交点形成三角形,该三角形就为种子三角形。最后在标记数组中标记该立方体已处理(整个方案拥有一个全局数组,用于记录所有立方体的处理状态,0表示当前立方体没有处理,为1表示当前立方体的信息已经放入了增长队列,为2表示当前立方体已经处理),把种子三角形压入队列中,之后将所求交点进行步骤D的操作,将种子三角形边的信息放入队列中。若中间层没有找到种子三角形,则去相邻其它层寻找,直到找到种子为止。
根据本发明的一种具体实施方式,在步骤B中,以前景边进行网格增长时,先从增长队列中取出数据,称数据中的边为前景边,根据数据中增长立方体的坐标构造立方体。将前景边的两个顶点按放入队列的循序分别命名为I点和II点,以I点→该等值面与立方体的新交点→II点的顺序依次求解该等值面与立方体的所有的新交点。按前景边在立方体中的位置,具体求解过程分为两种情况:
情况一:前景边的两个端点在相连的两条立方体的边上,交点的求解过程为:先找到前景边端点所在两条立方体边相连的立方体顶点,称该顶点为基准点;将基准点的灰度值与阈值进行比较,如果其大于阈值,则称灰度值大于阈值的立方体顶点满足要求的顶点,否者称灰度值不大于阈值的立方体顶点为满足要求的顶点。令前景边所在立方体的面为面SS,面SS的对面为面S_S,面SS上与基准点和I点共线的立方体顶点为SS_I点,面SS上与基准点和II点共线的立方体顶点为SS_II点,面SS上基准点的对角点为SS对交点,面S_S上与面SS上基准点、SS_I点、SS_II点、SS对角点相连的点分别命名为S_S基准点、S_S_I点、S_S_II点、S_S对角点。接下来根据判断条件(立方体上各个顶点灰度值与阈值进行比较,判断该顶点是否是满足要求的点)求等值面与立方体的交点:
a、S_S基准点不是满足要求的点。此时该等值面与立方体只有1个新的相交点,新交点在基准点和S_S基准点的边上插值求得。
b、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点也是满足要求的点。交点求解再细分为三种:①S_S对角点不是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S对角点、S_S对角点与S_S_II点、S_S_II点和SS_II点之间的边上与等值面插值求得;②S_S对角点也是满足要求的点,SS对角点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在SS_I点和S_S_I点、SS对角点点和S_S对角点、S_S_II点和SS_II点之间的边上与等值面插值求得;③S_S对角点和SS对角点也都是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、SS_I点和SS对角点、SS对角点与SS_II点、S_S_II点和SS_II点之间的边上与等值面插值求得。
c、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点中只有S_S_I点是满足要求的点。交点求解再细分为三种:①S_S对角点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S对角点、S_S_II点和S_S基准点之间的边上与等值面插值求得;②S_S对角点也是满足要求的点,但SS对角点不是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、SS对角点和S_S对角点、S_S对角点与S_S_II点、S_S_II点和S_S基准点之间的边上与等值面插值求得;③S_S对角点和SS对角点也都是满足要求的点,此时该等值面与立方体有5个新的相交点,他们分别在SS_I点和S_S_I点、SS_I点和SS对角点、SS对角点与SS_II点、S_S对角点与S_S_II点、S_S_II点和S_S基准点之间的边上与等值面插值求得。
d、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点中只有S_S_II点是满足要求的点。交点求解再细分为三种:①S_S对角点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在S_S基准点和S_S_I点、S_S_II点和S_S对角点、S_S_II点和SS_II点之间的边上与等值面插值求得;②S_S对角点也是满足要求的点,但SS对角点不是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在S_S基准点和S_S_I点、S_S_I点和S_S对角点、S_S对角点和SS对角点、S_S_II点和SS_II点之间的边上与等值面插值求得;③S_S对角点和SS对角点也都是满足要求的点,此时该等值面与立方体有5个新的相交点,他们分别在S_S基准点和S_S_I点、S_S_I点和S_S对角点、SS_I点和SS对角点、SS对角点和SS_II点、S_S_II点和SS_II点之间的边上与等值面插值求得。
e、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点都不是满足要求的点,此时该等值面与立方体有2个新的相交点,他们分别在S_S基准点和S_S_I点、S_S基准点和S_S_II之间的边上与等值面插值求得。
情况二:前景边的两个端点在一个面上相对的两条立方体的边上,交点的求解过程为:先将立方体所有的顶点灰度值与阈值进行比较,将顶点分为大于阈值和不大于阈值的两部分,取数量较小的那部分的顶点为满足要求的点,若两部分点的数量相同,取大于阈值的那部分点为满足要求的点。令前景边所在立方体的面为面SS,面SS的对面为面S_S,面SS上与I点同在一条边的满足要求的点为SS_I点、同在一条边的不是满足要求的点为SS_UP_I点,面SS上与II点同在一条边的满足要求的点为SS_II点、同在一条边的不是满足要求的点为SS_UP_II点,面S_S上与面SS上SS_I点、SS_II点、SS_UP_I点、SS_UP_II点相连的点分别命名为S_S_I点、S_S_II点、S_S_UP_I点、S_S_UP_II点。然后根据判断条件求等值面与立方体的交点:
a、S_S_I点和S_S_II点都不是满足要求的点,此时该等值面与立方体有2个新的相交点,他们分别在SS_I点和S_S_I点、SS_II点和S_S_II之间的边上与等值面插值求得。
b、S_S_I点和S_S_II点都是满足要求的点,此时该等值面与立方体有2个新的相交点,他们分别在S_S_UP_I点和S_S_I点、S_S_UP_II点和S_S_II点之间的边上与等值面插值求得。
c、S_S_I点和S_S_II点中只有S_S_I点是满足要求的点,此时需要再分两种情况:①S_S_UP_I点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在S_S_UP_I点和S_S_I点、S_S_I点和S_S_II点、S_S_II和SS_II点之间的边上与等值面插值求得;②S_S_UP_I点也是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_UP_I点和S_S_UP_I点、S_S_UP_I点和S_S_UP_II点、S_S_I点和S_S_II点、S_S_II和SS_II点之间的边上与等值面插值求得。
d、S_S_I点和S_S_II点中只有S_S_II点是满足要求的点,此时需要再分两种情况:①S_S_UP_II点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S_II点、S_S_II和S_S_UP_II点之间的边上与等值面插值求得;②S_S_UP_II点也是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S_II点、S_S_UP_I点和S_S_UP_II点、S_S_UP_II和SS_UP_II点之间的边上与等值面插值求得。
最后将前景边的两个端点按照II点和I点的顺序排列在求解的新交点后面,以便后续操作,同时标记该立方体已经被处理。
根据本发明的一种具体实施方式,在步骤C中,把交点按步骤B中的求解的顺序从数字0依次递增设定相应的标号,按交点数量分类:交点数为3时,直接三个点形成三角形;交点数为4时,标号为(0,1,2)和(2,3,0)的交点分别组成两个三角形;交点数为5时,标号为(0,1,2)、(2,3,4)和(4,0,2)的交点分别组成三个三角形;交点数为6时,标号为(0,1,2)、(2,3,4)、(4,5,0)、(0,2,4)的交点分别组成四个三角形;交点数为6时,标号为(0,1,2)、(2,3,4)、(4,5,0)、(0,2,4)的交点分别组成四个三角形;交点数为7时,标号为(0,1,2)、(2,3,4)、(4,5,6)、(6,0,2)、(2,4,6)的交点分别组成五个三角形。最后把所有三角形放入三角链表中。
根据本发明的一种具体实施方式,在步骤D中,增长边为等值面与立方体面新交线中那部分能增长的交线,它是相邻两交点组成的边,一条交线能不能增长取决于共用该交线的相邻立方体是否已经处理,只有相邻立方体没有被处理时,该交线才能够增长;要存储增长边的信息为增长边的两个端点坐标和对应相邻要增长的立方体的坐标。
根据本发明的一种具体实施方式,在步骤E中,只要队列中有数据,就从队列中取出数据进行步骤B、C和D的操作。在整个重建过程中,拥有一个全局的增长队列,它存储着增长边的两个端点坐标和对应要增长立方体的坐标,确保所有相连的网格都能够生长出来。
本发明的有益效果:通过以边增长的方式进行医学图像三维重建等值面网格提取,减少了对空体素的计算时间,避免了网格连接时的二义性问题;在使用上,由于是通过网格增长的原因,可以重建指定的组织。
附图说明
图1为从二维图像构成立方体的示意图;
图2为等值面与立方相交时,等值面一侧最少立方体顶点数量为一,在种子三角形选择时,以这种立方体中的三角形为种子;
图3为前景边在要增长的立方体面上的两种情况;
图4为前景边的两个端点在相连的两条边上时,立方体上各个顶点的命名情况;
图5-8为前景边的两个端点在相连的两条边上时,等值面与立方体边的相交情况;
图9为前景边的两个端点在相对的两条边上时,立方体上各个顶点的命名情况;
图10-13为前景边的两个端点在相对的两条边上时,等值面与立方体边的相交情况;
图14为将相交点组合三角形的示意图;
图15为将相交点两个相邻点一组,判断判断共享该两个相邻交点的相邻立方体是否被处理;
图16为已知立方体的一个顶点的下标情况下,其它顶点的小标对应情况。
具体实施方式
为使本发明的特点更加清晰,结合具体实例进行详细说明。
本发明为一种医学图像三维重建等值面提取方法,整体主要基于边增长来实现,主要过程可分为种子三角形的选择、以前景边出发求解等值面与增长立方体的交点、将交点组合成三角形、将相邻两个交点组成的前景边信息放入队列、从队列取出前景边进行增长,具体步骤如下:
A、读取医学影像数据,构造立方体,选取种子三角形;
B、以前景边进行网格增长,求解以前景边增长时等值面与立方体边的交点;
C、将所有交点按求解的顺序3个点一起形成三角形,并将三角形放入三角链表中;
D、判断所有两个相邻交点组合的边是否为增长边,如果是增长边,把增长边的信息放入队列中;
E、从队列中取出数据进行增长,直到队列为空。
在重建算法的实现上,整个重建过程主要分为两个模块,一个模块用于种子三角形的选取,另一模块用于从增长队列中取出数据进行网格增长。
种子选取对应的模块具体实现的内容为:从输入的医学断层图像构成的立方体中,通过相邻两层医学断层图像各4个相邻像素点构造中间层所有立方体(如图1所示)和遍历中间层所有的立方体。对于中间每个立方体,将它的8个顶点对应的灰度值和阈值进行比大小,并记录大于阈值顶点数量size,只有当size为1或者7时,才认为该立方体中存在种子三角形,记size为1时对应大于阈值顶点在该立方体的下标为(i,j,k),或记size为7时对应不大于阈值顶点在该立方体的下标为(i,j,k),其中i、j、k取值都为0或者1,并称被记录坐标的点为基准点。当找到基准点后,其他点的下标可以通过(i+1)%2、(j+1)%2、(k+1)%2求得,如图16。接下来通过找到的基准点,分别插值求出和基准点相连三条立方体的边与等值面的交点,将三个交点相连组成三角形,即为所求的种子三角形,如图2所示。最后,将种子三角形的三条边的信息放入增长队列,记录对应立方体的状态为2(表示已处理),相邻要增长的立方体的状态为1(表示信息已放入增长队列中)。
网格增长对应的模块具体实现内容为:先判断增长队列是否为空,不为空取出数据,根据中数据中的相邻要增长立方体的坐标构造出立方体,如果立方体没有被处理,就进行网格增长。从增长队列中取出前景边的两个端点坐标,以前景边的两个顶点按放入队列的循序分别命名为I点和II点,以I点→该等值面与立方体的新交点→II点的顺序依次求解该等值面与立方体的所有的新交点。按前景边在立方体中的位置,具体求解过程分为两种情况:
情况一:前景边的两个端点在相连的两条立方体的边上,如图3中的a所示,交点的求解过程为:先找到前景边端点所在两条立方体边相连的立方体顶点,如图3中的a中的4号点,称该顶点为基准点,代码实现上主要是找出通过前景边的两个端点坐标和立方体的各个顶点坐标找出基准点在立方体中的下标(i,j,k)。将基准点的灰度值与阈值进行比较,如果其大于阈值,则称灰度值大于阈值的立方体顶点满足要求的顶点,否者称灰度值不大于阈值的立方体顶点为满足要求的顶点。令前景边所在立方体的面为面SS,面SS的对面为面S_S,面SS上与基准点和I点共线的立方体顶点为SS_I点,面SS上与基准点和II点共线的立方体顶点为SS_II点,面SS上基准点的对角点为SS对交点,面S_S上与面SS上基准点、SS_I点、SS_II点、SS对角点相连的点分别命名为S_S基准点、S_S_I点、S_S_II点、S_S对角点,如图4所示,代码实现时,用一个数组记录各个点对应在立方体中的下标,以便后续环节代码的实现,其它顶点的下标由基准点的下标(i,j,k)计算得到,计算原理对应图16。接下来根据判断条件(立方体上各个顶点灰度值与阈值进行比较,判断该顶点是否是满足要求的点)求等值面与立方体的交点:
a、S_S基准点不是满足要求的点。此时该等值面与立方体只有1个新的相交点,新交点在基准点和S_S基准点的边上插值求得,如图2所示。
b、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点也是满足要求的点。交点求解再细分为三种:①S_S对角点不是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S对角点、S_S对角点与S_S_II点、S_S_II点和SS_II点之间的边上与等值面插值求得,如图5中a所示;②S_S对角点也是满足要求的点,SS对角点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在SS_I点和S_S_I点、SS对角点点和S_S对角点、S_S_II点和SS_II点之间的边上与等值面插值求得,如图5中b所示;③S_S对角点和SS对角点也都是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、SS_I点和SS对角点、SS对角点与SS_II点、S_S_II点和SS_II点之间的边上与等值面插值求得,如图5中c所示。
C、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点中只有S_S_I点是满足要求的点。交点求解再细分为三种:①S_S对角点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S对角点、S_S_II点和S_S基准点之间的边上与等值面插值求得,如图6中a所示;②S_S对角点也是满足要求的点,但SS对角点不是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、SS对角点和S_S对角点、S_S对角点与S_S_II点、S_S_II点和S_S基准点之间的边上与等值面插值求得,如图6中b所示;③S_S对角点和SS对角点也都是满足要求的点,此时该等值面与立方体有5个新的相交点,他们分别在SS_I点和S_S_I点、SS_I点和SS对角点、SS对角点与SS_II点、S_S对角点与S_S_II点、S_S_II点和S_S基准点之间的边上与等值面插值求得,如图6中c所示。
d、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点中只有S_S_II点是满足要求的点。交点求解再细分为三种:①S_S对角点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在S_S基准点和S_S_I点、S_S_II点和S_S对角点、S_S_II点和SS_II点之间的边上与等值面插值求得,如图7中a所示;②S_S对角点也是满足要求的点,但SS对角点不是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在S_S基准点和S_S_I点、S_S_I点和S_S对角点、S_S对角点和SS对角点、S_S_II点和SS_II点之间的边上与等值面插值求得,如图7中b所示;③S_S对角点和SS对角点也都是满足要求的点,此时该等值面与立方体有5个新的相交点,他们分别在S_S基准点和S_S_I点、S_S_I点和S_S对角点、SS_I点和SS对角点、SS对角点和SS_II点、S_S_II点和SS_II点之间的边上与等值面插值求得,如图7中c所示。
e、S_S面上S_S基准点是满足要求的点,且S_S_I点和S_S_II点都不是满足要求的点,此时该等值面与立方体有2个新的相交点,他们分别在S_S基准点和S_S_I点、S_S基准点和S_S_II之间的边上与等值面插值求得,如图8所示。
情况二:前景边的两个端点在一个面上相对的两条立方体的边上,如图3中的b所示,交点的求解过程为:先将立方体所有的顶点灰度值与阈值进行比较,将顶点分为大于阈值和不大于阈值的两部分,取数量较小的那部分的顶点为满足要求的点,若两部分点的数量相同,取大于阈值的那部分点为满足要求的点。令前景边所在立方体的面为面SS,面SS的对面为面S_S,面SS上与I点同在一条边的满足要求的点为SS_I点、同在一条边的不是满足要求的点为SS_UP_I点,面SS上与II点同在一条边的满足要求的点为SS_II点、同在一条边的不是满足要求的点为SS_UP_II点,面S_S上与面SS上SS_I点、SS_II点、SS_UP_I点、SS_UP_II点相连的点分别命名为S_S_I点、S_S_II点、S_S_UP_I点、S_S_UP_II点,如图9所示,代码实现时同情况一一样,用一个数组记录各个点对应在立方体中的下标,以便后续环节代码的实现,各个点的下标求解过程为:先根据前景边的端点坐标和立方体各个顶点的坐标求出SS_I点和SS_II点的下标,其它顶点的下标通过与它俩的空间相互关系(如图16)得到。接下来根据判断条件(立方体上各个顶点灰度值与阈值进行比较,判断该顶点是否是满足要求的点)求等值面与立方体的交点:
a、S_S_I点和S_S_II点都不是满足要求的点,此时该等值面与立方体有2个新的相交点,他们分别在SS_I点和S_S_I点、SS_II点和S_S_II之间的边上与等值面插值求得,如图10所示。
b、S_S_I点和S_S_II点都是满足要求的点,此时该等值面与立方体有2个新的相交点,他们分别在S_S_UP_I点和S_S_I点、S_S_UP_II点和S_S_II点之间的边上与等值面插值求得,如图11所示。
c、S_S_I点和S_S_II点中只有S_S_I点是满足要求的点,此时需要再分两种情况:①S_S_UP_I点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在S_S_UP_I点和S_S_I点、S_S_I点和S_S_II点、S_S_II和SS_II点之间的边上与等值面插值求得,如图12中a所示;②S_S_UP_I点也是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_UP_I点和S_S_UP_I点、S_S_UP_I点和S_S_UP_II点、S_S_I点和S_S_II点、S_S_II和SS_II点之间的边上与等值面插值求得,如图12中b所示。
d、S_S_I点和S_S_II点中只有S_S_II点是满足要求的点,此时需要再分两种情况:①S_S_UP_II点不是满足要求的点,此时该等值面与立方体有3个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S_II点、S_S_II和S_S_UP_II点之间的边上与等值面插值求得,如图13中a所示;②S_S_UP_II点也是满足要求的点,此时该等值面与立方体有4个新的相交点,他们分别在SS_I点和S_S_I点、S_S_I点和S_S_II点、S_S_UP_I点和S_S_UP_II点、S_S_UP_II和SS_UP_II点之间的边上与等值面插值求得,如图13中b所示。
最后将前景边的两个端点按照II点和I点的顺序排列在求解的新交点后面,以便后续操作,同时标记该立方体已经被处理。
交点求解完后,把各个交点组成三角形,然后放在三角链表中。把交点按求解的顺序从数字0依次递增设定相应的标号,按交点数量分类:交点数为3时,直接三个点形成三角形;交点数为4时,标号为(0,1,2)和(2,3,0)的交点分别组成两个三角形;交点数为5时,标号为(0,1,2)、(2,3,4)和(4,0,2)的交点分别组成三个三角形;交点数为6时,标号为(0,1,2)、(2,3,4)、(4,5,0)、(0,2,4)的交点分别组成四个三角形;交点数为6时,标号为(0,1,2)、(2,3,4)、(4,5,0)、(0,2,4)的交点分别组成四个三角形;交点数为7时,标号为(0,1,2)、(2,3,4)、(4,5,6)、(6,0,2)、(2,4,6)的交点分别组成五个三角形,如图14所示。
三角形组成完成后,对所有相邻两点组成的边进行处理,判断边是否为增长边,如果是增长边,把两邻两点的坐标和相邻立方体的下标放入队列中。一条边是不是增长边取决于共用该边的相邻立方体是否已经处理。判断它是否已处理,只需找出相邻立方体在所有立方体中对应的坐标,在根据坐标找出对应处理情况。相邻立方体的坐标求解过程为:首先,先找到两个相邻点所在的平面,令垂直该平面的坐标方向为方向i(相邻两点坐标在i方向的数值相等,点和立方体的坐标都有三个方向,对应笛卡尔坐标系的x,y,z),然后和相邻第三点比较各自坐标在i方向上的数值大小,若两相邻点在i方向的坐标数值小于第三点在i方向的坐标数值,则所有相邻立方体的坐标在i方向上的数值为本立方体坐标在i方向数字减一,否则所求相邻立方体的坐标在i方向上的数值为本立方体坐标在i方向数字加一,相邻立方体坐标在其他两个方向数值等于本立方体在这两个方向的数值,如图15,以相邻点b和c为例,b和c的坐标在x方向相等,但b和c的x坐标小于a点x的坐标,说明相邻立方体坐标x方向等于本立方体x坐标减一,其他两个方向相等,若本立方体坐标为(aa,bb,cc),则相邻立方体坐标为(aa-1,bb,cc)。接下来根据所求坐标在标记数组中找出相邻立方体对应的数值,如果对应标记数组中的值为0,说明该边为增长边,先将相应标记数组的数值置为1,再将所求相邻立方体的下标和相邻两交点的坐标按交点求解顺序放入队列中;如果对应标记数组中的值不为0,该边不为增长边,则不再进行处理。
最后,不断从队列中求出数据进行增长,直到队列为空。
本领域技术人员可以对本发明有多种变形方案实现本发明,以上只是本发明的一种可行的实例而已,并非因此限制本发明的权力保护范围。
Claims (6)
1.一种基于边增长的医学图像三维重建等值面网格提取方法,其特征在于,包括以下步骤:
A、读取医学图像数据,构造立方体,选取种子三角形;
B、以前景边进行网格增长,求解以前景边增长时等值面与立方体边的交点;
C、将所有交点按求解的顺序3个点一起形成三角形,并将三角形放入三角链表中;
D、判断所有两个相邻交点组成的边是否为增长边,如果是增长边,把增长边的信息放入队列中;
E、从队列中取出数据进行增长,直到队列为空。
2.根据权利要求1所述的一种基于边增长的医学图像三维重建等值面网格提取方法,其特征在于,在所述步骤A中,种子三角形的选取过程为:从二维断层图像构成的所有立方体中,遍历中间一层的所有立方体,将立方体8个顶点的灰度值分别与等值面的灰度值(等值面的灰度值简称阈值)进行比较,并将顶点分类为大于阈值和不大于阈值两部分,只有当两部分中顶点数量较小的那个部分的点的数量刚好为1时,才认为该立方体中有种子三角形,并称数量为1对应的那个点为基本点,以基本点为基准,分别插值求出和基本点相连三条立方体的边与等值面的交点,然后将所求的三个交点形成三角形,该三角形就为种子三角形;若中间层没有找到种子三角形,则去相邻其它层寻找,直到找到种子三角形为止。
3.根据权利要求1所述的一种基于边增长的医学图像三维重建等值面网格提取方法,其特征在于,在所述步骤B中,以前景边进行增长时,前景边为从队列取出数据中对应的边:将前景边的两个顶点按放入队列的循序分别命名为I点和II点,以I点→该等值面与立方体的新交点→II点的顺序依次求解该等值面与立方体的所有的新交点;具体求解等值面与立方体的交点过程分为两种情况:情况一为前景边的两个端点在相连的两条立方体的边上,情况二为前景边的两个端点在一个面上相对的两条立方体边上。
4.根据权利要求1所述的一种基于边增长的医学图像三维重建等值面网格提取方法,其特征在于,在所述步骤D中,增长边为等值面与立方体面新交线中那部分能增长的交线,一条交线能不能增长取决于共用该交线的相邻立方体是否已经处理,只有相邻立方体没有被处理时,该交线才能够增长。
5.根据权利要求1所述的一种基于边增长的医学图像三维重建等值面网格提取方法,其特征在于,在所述步骤D中,要放入队列的增长边信息为:当确定相邻两个交点组成的边为增长边后,需要将边的两个端点坐标按求解顺序放入队列,同时一同放入队列的还有对应相邻要增长的立方体的坐标。
6.根据权利要求1所述的一种基于边增长的医学图像三维重建等值面网格提取方法,其特征在于,在所述步骤D中,在整个重建过程中,拥有一个队列,它存储着增长边的两个端点坐标和对应相邻要增长立方体的坐标,确保所有相连的网格都能够生长出来。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811560982.4A CN109636915B (zh) | 2018-12-19 | 2018-12-19 | 一种基于边增长的医学图像三维重建等值面网格提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811560982.4A CN109636915B (zh) | 2018-12-19 | 2018-12-19 | 一种基于边增长的医学图像三维重建等值面网格提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109636915A true CN109636915A (zh) | 2019-04-16 |
CN109636915B CN109636915B (zh) | 2020-11-20 |
Family
ID=66075804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811560982.4A Expired - Fee Related CN109636915B (zh) | 2018-12-19 | 2018-12-19 | 一种基于边增长的医学图像三维重建等值面网格提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109636915B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110610478A (zh) * | 2019-07-24 | 2019-12-24 | 浙江大学 | 一种基于邻域拓扑的医学图像三维重建方法 |
CN113888700A (zh) * | 2021-10-20 | 2022-01-04 | 哈尔滨理工大学 | 一种基于体素增长的医学图像三维重建方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101692286A (zh) * | 2009-09-15 | 2010-04-07 | 上海交通大学 | 医学图像三视图获取方法 |
CN101814195A (zh) * | 2010-04-19 | 2010-08-25 | 李楚雅 | 三维建模方法 |
CN101877147A (zh) * | 2010-06-29 | 2010-11-03 | 浙江大学 | 三维三角形网格模型的简化算法 |
CN103679814A (zh) * | 2012-08-31 | 2014-03-26 | 西门子公司 | 一种三角形网格模型的处理方法及装置 |
CN104268934A (zh) * | 2014-09-18 | 2015-01-07 | 中国科学技术大学 | 一种由点云直接重建三维曲面的方法 |
US20150154796A1 (en) * | 2012-08-17 | 2015-06-04 | Google Inc. | Accelerated three-dimensional intersection testing using compressed volumes |
CN108510443A (zh) * | 2018-03-30 | 2018-09-07 | 河北北方学院 | 一种医学图像离线重建定位方法 |
US20180276316A1 (en) * | 2017-03-23 | 2018-09-27 | Autodesk, Inc. | Creating gradients of different materials for three-dimensional models in computer aided design applications |
CN108694726A (zh) * | 2017-04-05 | 2018-10-23 | 贵州利体数字医学科技有限公司 | 医学影像三维重建系统 |
US20180315242A1 (en) * | 2017-04-28 | 2018-11-01 | Baker Hughes Incorporated | Three dimensional surface representation for modeling geological surfaces |
-
2018
- 2018-12-19 CN CN201811560982.4A patent/CN109636915B/zh not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101692286A (zh) * | 2009-09-15 | 2010-04-07 | 上海交通大学 | 医学图像三视图获取方法 |
CN101814195A (zh) * | 2010-04-19 | 2010-08-25 | 李楚雅 | 三维建模方法 |
CN101877147A (zh) * | 2010-06-29 | 2010-11-03 | 浙江大学 | 三维三角形网格模型的简化算法 |
US20150154796A1 (en) * | 2012-08-17 | 2015-06-04 | Google Inc. | Accelerated three-dimensional intersection testing using compressed volumes |
CN103679814A (zh) * | 2012-08-31 | 2014-03-26 | 西门子公司 | 一种三角形网格模型的处理方法及装置 |
CN104268934A (zh) * | 2014-09-18 | 2015-01-07 | 中国科学技术大学 | 一种由点云直接重建三维曲面的方法 |
US20180276316A1 (en) * | 2017-03-23 | 2018-09-27 | Autodesk, Inc. | Creating gradients of different materials for three-dimensional models in computer aided design applications |
CN108694726A (zh) * | 2017-04-05 | 2018-10-23 | 贵州利体数字医学科技有限公司 | 医学影像三维重建系统 |
US20180315242A1 (en) * | 2017-04-28 | 2018-11-01 | Baker Hughes Incorporated | Three dimensional surface representation for modeling geological surfaces |
CN108510443A (zh) * | 2018-03-30 | 2018-09-07 | 河北北方学院 | 一种医学图像离线重建定位方法 |
Non-Patent Citations (3)
Title |
---|
GUO L: "《Three dimension reconstruction of medical images based on an improved marching cubes algorithm》", 《2013 6TH INTERNATIONAL CONFERENCE ON BIOMEDICAL ENGINEERING AND INFORMATICS》 * |
WANG M: "《Three-dimensional computational model simulating the fracture healing process with both biphasic poroelastic finite element analysis and fuzzy logic control》", 《SCIENTIFIC REPORTS》 * |
谭小俊: "《大规模点云的处理及曲面重建技术》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110610478A (zh) * | 2019-07-24 | 2019-12-24 | 浙江大学 | 一种基于邻域拓扑的医学图像三维重建方法 |
CN110610478B (zh) * | 2019-07-24 | 2021-05-25 | 浙江大学 | 一种基于邻域拓扑的医学图像三维重建方法 |
CN113888700A (zh) * | 2021-10-20 | 2022-01-04 | 哈尔滨理工大学 | 一种基于体素增长的医学图像三维重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109636915B (zh) | 2020-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112836618B (zh) | 一种三维人体姿态估计方法及计算机可读存储介质 | |
CN104268934B (zh) | 一种由点云直接重建三维曲面的方法 | |
CN104835202A (zh) | 一种三维虚拟场景快速构建方法 | |
CN101339670A (zh) | 一种计算机辅助的三维颅面复原方法 | |
CN109147025B (zh) | 一种面向rgbd三维重建的纹理生成方法 | |
CN109636915A (zh) | 一种基于边增长的医学图像三维重建等值面网格提取方法 | |
CN108280858A (zh) | 多视图重建中的一种线性全局相机运动参数估计方法 | |
CN112529950B (zh) | 一种用于辐射剂量评估的生物体体素模型实现方法 | |
CN106327479A (zh) | 血管造影中介下先心病术中血管辨识的装置及方法 | |
CN107203988A (zh) | 一种由二维x光图像重建三维体图像的方法及其应用 | |
CN110610478A (zh) | 一种基于邻域拓扑的医学图像三维重建方法 | |
CN110176079A (zh) | 一种基于准共形映射的三维模型变形算法 | |
CN109509259B (zh) | 一种医学图像三维重建等值面网格提取方法 | |
CN109919885A (zh) | 基于b样条和互信息的cvh图像与影像配准融合方法 | |
CN101996415B (zh) | 眼球的三维建模方法 | |
CN108038904A (zh) | 医学影像三维重建系统 | |
CN103345774A (zh) | 一种三维多尺度矢量化的建模方法 | |
CN106097446A (zh) | 一种基于大规模地震数据的三维地质曲面重建方法 | |
CN114494637B (zh) | 一种基于结构体矩阵的砂岩三维真实模型重构方法 | |
CN110322548A (zh) | 一种基于几何图像参数化的三维网格模型生成方法 | |
CN111369662A (zh) | Ct图像中血管的三维模型重建方法及系统 | |
CN109325977A (zh) | 目标区域最优影像选择方法及系统、存储介质、电子设备 | |
CN108538370A (zh) | 一种光照体绘制输出方法及装置 | |
CN113421339A (zh) | 一种基于体素增长的医学图像三维重建等值面网格加速提取方法 | |
CN106097447A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201120 Termination date: 20211219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |