CN102982546A - 断层图像的分割方法和装置 - Google Patents

断层图像的分割方法和装置 Download PDF

Info

Publication number
CN102982546A
CN102982546A CN2012104973862A CN201210497386A CN102982546A CN 102982546 A CN102982546 A CN 102982546A CN 2012104973862 A CN2012104973862 A CN 2012104973862A CN 201210497386 A CN201210497386 A CN 201210497386A CN 102982546 A CN102982546 A CN 102982546A
Authority
CN
China
Prior art keywords
segmentation
man
sequence
cut apart
auto
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
Application number
CN2012104973862A
Other languages
English (en)
Other versions
CN102982546B (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.)
Haina Medical Faith (beijing) Software Technology Co Ltd
Original Assignee
Haina Medical Faith (beijing) Software 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 Haina Medical Faith (beijing) Software Technology Co Ltd filed Critical Haina Medical Faith (beijing) Software Technology Co Ltd
Priority to CN201210497386.2A priority Critical patent/CN102982546B/zh
Publication of CN102982546A publication Critical patent/CN102982546A/zh
Application granted granted Critical
Publication of CN102982546B publication Critical patent/CN102982546B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了一种断层图像的分割方法和装置。其中,断层图像的分割方法通过根据用户的选择指令分别确定出两个人机交互分割层面,并对确定出的这两个人机交互分割层面进行人机交互分割,同时还对这两个人机交互分割层面之间的中间层面进行自动分割。通过本发明,解决了现有技术中断层图像的分割方法分割精度较低的问题,进而达到了提高断层图像分割精确度的效果。

Description

断层图像的分割方法和装置
技术领域
本发明涉及断层图像领域,具体而言,涉及一种断层图像的分割方法和装置。
背景技术
在临床诊断中,计算机断层扫描(Computed Tomography,简称CT),核磁共振成像(Magnetic Resonance Imaging,简称MRI)和正电子发射型计算机断层显像(PositronEmission Computed Tomography,简称PET)应用日益广泛。这些设备所产生的图像,统称为断层图像,一般为选择人体某一部位,按照一定方向进行等间距扫描,得到一系列二维断层图像。每张断层图像,对应人体中的一个断层,由灰度像素排列组成。灰度像素的灰度值,根据扫描设备的不同,反应人体内组织的不同信息。在临床医疗和科研工作中,对医学图像中感兴趣的解剖结构或器官进行三维分割,具有重要意义。解剖结构或器官的分割结果,不仅能够用于三维重建可视化,还可以进行多种三维空间测量。并且,三维分割结果,在骨骼校正、整形手术、病灶切除、手术规划、放疗规划中具有重要价值。现有技术中根据断层图像序列进行器官三维分割的方法众多,其中一类常用的分割方法原理为:
1)采用序列图像作为输入图像,以序列图像中感兴趣结构或器官作为分割目标;
2)人工选取分割目标所在的某一层面,作为初始层面;
3)通过人机交互操作,在初始层面中,指定分割算法所需要的初始输入,包括初始种子点、初始轮廓或区域等;
4)采用某种二维分割方法得到分割目标在初始层面中的二维分割结果;
5)采用相邻层面“逐层传递分割”方式,将上一层面的分割结果作为下一层面的初始值,启动下一层面的分割过程;
6)完成分割目标在所有层面的分割过程后,将分割目标在各个层面的二维分割结果,组合构成三维分割结果。
上述现有技术中的分割方法存在的问题是:由于人机交互操作被限制在单一层面,即只允许存在一个人机交互分割层面,在“逐层传递分割”过程中,受各种可能干扰因素的影响,分割算法在某一层面的分割结果有可能发生错误,偏离分割目标实际区域。并且,一旦发生分割偏差,将被传递到下一层面和后续层面,有可能逐层放大。最终三维分割结果,和分割目标相比,存在较大偏差,甚至分割失败。此外,一旦分割算法在某一层面分割失败,将造成后续层面分割无法继续。
即,由于只允许存在一个人机交互分割层面,“逐层传递分割”过程一旦启动,无法再进行人机交互或人工校正,分割结果的稳定性和抗干扰性较低,从而影响分割精度,或存在分割失败可能。
针对相关技术中断层图像的分割方法分割精度较低的问题,目前尚未提出有效的解决方案。
发明内容
本发明的主要目的在于提供一种断层图像的分割方法和装置,以解决现有技术中断层图像的分割方法分割精度较低的问题。
为了实现上述目的,根据本发明的一个方面,提供了一种断层图像的分割方法,包括:接收来自用户的第一选择指令,并从断层图像中获取与第一选择指令相对应的分割目标;接收来自用户的第二选择指令并从断层图像中确定与第二选择指令相对应的层面,得到第一人机交互分割层面;对第一人机交互分割层面进行分割,得到分割目标在第一人机交互分割层面内的分割结果;接收来自用户的第三选择指令并从断层图像中确定与第三选择指令相对应的层面,得到第二人机交互分割层面;对第二人机交互分割层面进行分割,得到分割目标在第二人机交互分割层面内的分割结果;获取第一人机交互分割层面与第二人机交互分割层面之间的中间层面,得到第一自动分割层面;对第一自动分割层面进行分割,得到分割目标在第一自动分割层面内的分割结果;以及将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果和分割目标在第一自动分割层面内的分割结果进行组合,得到分割目标的三维分割结果。
进一步地,对第一自动分割层面进行分割包括:获取第一分割轮廓和第二分割轮廓,其中,第一分割轮廓为分割目标在第一人机交互分割层面内的分割结果对应的分割轮廓,第二分割轮廓为分割目标在第二人机交互分割层面内的分割结果对应的分割轮廓;以及以第一分割轮廓和第二分割轮廓为基础采用层间轮廓自动插值法对第一自动分割层面进行分割。
进一步地,第一自动分割层面包括一个或多个分割层面,以第一分割轮廓和第二分割轮廓为基础采用层间轮廓自动插值法对第一自动分割层面进行分割包括:对第一分割轮廓的曲线长度进行预设份数等分,得到多个第一等分点;分别获取多个第一等分点在第一人机交互分割层面中的图像坐标,得到第一序列;对第二分割轮廓的曲线长度进行预设份数等分,得到多个第二等分点;分别获取多个第二等分点在第二人机交互分割层面中的图像坐标,得到第二序列;确定第一序列中各个序列点与第二序列中各个序列点的对应关系;以及采用预设插值公式按照确定出的对应关系确定分割目标在第一层面内的分割轮廓,得到分割目标在第一层面内的分割结果,其中,第一层面为第一自动分割层面中的任一层面。
进一步地,确定第一序列中各个序列点与第二序列中各个序列点的对应关系包括:计算第一分割轮廓区域的质心,得到第一质心;计算第二分割轮廓区域的质心,得到第二质心;以及通过以下公式确定第一序列中各个序列点与第二序列中各个序列点的对应关系:
D ( i ) = Σ k = 0 N - 1 | P ( mod ( k + i , N ) ) - P C , Qk - Q C | 2 , i = 0,1,2 , . . . N - 1 i ′ = min i ( D ( i ) ) ,
其中,mod(k+i,N)表示取(k+i)被N除的余数,P(mod(k+i,N))=(PX(mod(k+i,N)),PY(mod(k+i,N)))表示第一序列中各个序列点的图像坐标,PC=(PXC,PYC)表示第一质心的图像坐标,P(mod(k+i,N))-PC=(PX(mod(k+i,N))-PXC,PY(mod(k+i,N))-PYC)表示P(mod(k+i,N))相对PC的坐标,Qk=(QX(k),QY(k))表示第二序列中各个序列点的图像坐标,QC=(QCX,QCY)表示第二质心的图像坐标,Qk-QC=(QX(k)-QXC,QY(k)-QYC)表示Qk相对QC的坐标,|A,B|表示坐标点A和B之间的距离,N为预设份数,i′=mini(D(i))表示i=i'时D(i)取值最小,第一序列中的序列点P(mod(k+i′,N))与第二序列中的序列点Qk相对应,
预设插值公式为:
L ( k ) = LX ( k ) = QX ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ PX ( mod ( k + i ′ , N ) ) - QX ( k ) ] LY ( k ) = QY ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ PY ( mod ( k + i ′ , N ) ) - QY ( k ) ] , k = 0,1,2 , . . . , N - 1 ,
其中,L(k)=(LX(k),LY(k))表示分割目标在第一层面内的分割轮廓上的序列点,D1表示第一层面距离第一人机交互分割层面的垂直距离,D2表示第一层面距离第二人机交互分割层面的垂直距离。
进一步地,确定第一序列中各个序列点与第二序列中各个序列点的对应关系包括:计算第一分割轮廓区域的质心,得到第一质心;计算第二分割轮廓区域的质心,得到第二质心;以及通过以下公式确定第一序列中各个序列点与第二序列中各个序列点的对应关系:
D ( i ) = Σ k = 0 N - 1 | P ( mod ( k + i , N ) ) - P C , Qk - Q C | 2 , i = 0,1,2 , . . . N - 1 i ′ = min i ( D ( i ) ) ,
其中,mod(k+i,N)表示取(k+i)被N除的余数,P(mod(k+i,N))=(PX(mod(k+i,N)),PY(mod(k+i,N)))表示第一序列中各个序列点的图像坐标,PC=(PXC,PYC)表示第一质心的图像坐标,P(mod(k+i,N))-PC=(PX(mod(k+i,N))-PXC,PY(mod(k+i,N))-PYC)表示P(mod(k+i,N))相对PC的坐标,Qk=(QX(k),QY(k))表示第二序列中各个序列点的图像坐标,QC=(QCX,QCY)表示第二质心的图像坐标,Qk-QC=(QX(k)-QXC,QY(k)-QYC)表示Qk相对QC的坐标,|A,B|表示坐标点A和B之间的距离,N为预设份数,i′=mini(D(i))表示i=i'时D(i)取值最小,第一序列中的序列点P(mod(k+i′,N))与第二序列中的序列点Qk相对应,在确定第一序列中各个序列点与第二序列中各个序列点的对应关系之后,分割方法还包括:采用公式P′k=P(mod(k+i′,N))对第一序列中的各个序列点进行重新编号,得到第一序列中的各个序列点P′k,
预设插值公式为:
L ( k ) = LX ( k ) = QX ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ P ′ X ( k ) - QX ( k ) ] LY ( k ) = QY ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ P ′ Y ( k ) - QY ( k ) ] , k = 0,1,2 , . . . , N - 1
其中,L(k)=(LX(k),LY(k))表示分割目标在第一层面内的分割轮廓上的序列点,D1表示第一层面距离第一人机交互分割层面的垂直距离,D2表示第一层面距离第二人机交互分割层面的垂直距离。
进一步地,第一自动分割层面包括多个分割层面,分割目标在第一自动分割层面内的分割结果包括多个分割结果,在得到分割目标在第一自动分割层面内的分割结果之后,分割方法还包括:接收来自用户的第四选择指令并从多个分割结果中确定与第四选择指令相对应的分割结果,其中,第四选择指令用于选择多个分割结果中偏差最大的分割结果;获取与第四选择指令相对应的分割结果所对应的分割层面,得到第三人机交互分割层面;对第三人机交互分割层面进行分割,得到分割目标在第三人机交互分割层面内的分割结果;获取第一人机交互分割层面与第三人机交互分割层面之间的中间层面,得到第二自动分割层面;对第二自动分割层面进行分割,得到分割目标在第二自动分割层面内的分割结果;获取第二人机交互分割层面与第三人机交互分割层面之间的中间层面,得到第三自动分割层面;以及对第三自动分割层面进行分割,得到分割目标在第三自动分割层面内的分割结果,其中,将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果、分割目标在第三人机交互分割层面内的分割结果、分割目标在第二自动分割层面内的分割结果和分割目标在第三自动分割层面内的分割结果进行组合,得到三维分割结果。
进一步地,分割目标在第二自动分割层面内的分割结果以及分割目标在第三自动分割层面内的分割结果均包括多个分割结果,在得到分割目标在第二自动分割层面内的分割结果,以及得到分割目标在第三自动分割层面内的分割结果之后,分割方法还包括:重新接收第四选择指令以分别从第二自动分割层面内的多个分割结果和第三自动分割层面内的多个分割结果中确定与第四选择指令相对应的分割结果。
进一步地,在从断层图像中获取与第一选择指令相对应的分割目标之后,分割方法还包括:接收来自用户的第五选择指令并从断层图像中确定与第五选择指令相对应的层面方向,其中,第一人机交互分割层面、第二人机交互分割层面和第一自动分割层面均为沿层面方向上的分割层面。
为了实现上述目的,根据本发明的另一方面,提供了一种断层图像的分割装置,该分割装置用于执行本发明上述内容所提供的任一种分割方法。
为了实现上述目的,根据本发明的另一方面,提供了一种断层图像的分割装置,包括:第一接收单元,用于接收来自用户的第一选择指令,并从断层图像中获取与第一选择指令相对应的分割目标;第二接收单元,用于接收来自用户的第二选择指令并从断层图像中确定与第二选择指令相对应的层面,得到第一人机交互分割层面;第一分割单元,用于对第一人机交互分割层面进行分割,得到分割目标在第一人机交互分割层面内的分割结果;第三接收单元,用于接收来自用户的第三选择指令并从断层图像中确定与第三选择指令相对应的层面,得到第二人机交互分割层面;第二分割单元,用于对第二人机交互分割层面进行分割,得到分割目标在第二人机交互分割层面内的分割结果;第一获取单元,用于获取第一人机交互分割层面与第二人机交互分割层面之间的中间层面,得到第一自动分割层面;第三分割单元,用于对第一自动分割层面进行分割,得到分割目标在第一自动分割层面内的分割结果;以及组合单元,用于将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果和分割目标在第一自动分割层面内的分割结果进行组合,得到分割目标的三维分割结果。
进一步地,第三分割单元包括:获取子单元,用于获取第一分割轮廓和第二分割轮廓,其中,第一分割轮廓为分割目标在第一人机交互分割层面内的分割结果对应的分割轮廓,第二分割轮廓为分割目标在第二人机交互分割层面内的分割结果对应的分割轮廓;以及分割子单元,用于以第一分割轮廓和第二分割轮廓为基础采用层间轮廓自动插值法对第一自动分割层面进行分割。
通过本发明,采用接收来自用户的第一选择指令,并从断层图像中获取与第一选择指令相对应的分割目标;接收来自用户的第二选择指令并从断层图像中确定与第二选择指令相对应的层面,得到第一人机交互分割层面;对第一人机交互分割层面进行分割,得到分割目标在第一人机交互分割层面内的分割结果;接收来自用户的第三选择指令并从断层图像中确定与第三选择指令相对应的层面,得到第二人机交互分割层面;对第二人机交互分割层面进行分割,得到分割目标在第二人机交互分割层面内的分割结果,其中,对第一人机交互分割层面及第二人机交互分割层面的分割可以采用多种现有人机交互方法;获取第一人机交互分割层面与第二人机交互分割层面之间的中间层面,得到第一自动分割层面;对第一自动分割层面进行分割,得到分割目标在第一自动分割层面内的分割结果;以及将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果和分割目标在第一自动分割层面内的分割结果进行组合,得到分割目标的三维分割结果。通过根据用户的选择指令分别确定出两个人机交互分割层面,并对确定出的这两个人机交互分割层面进行人机交互分割,同时还对这两个人机交互分割层面之间的中间层面进行自动分割,实现了允许用户指定一个以上的人机交互分割层面,避免了现有技术中只允许存在一个人机交互分割层面所带来的以下弊端:在“逐层传递分割”过程中,受各种可能干扰因素的影响,分割算法在某一层面的分割结果有可能发生错误,偏离分割目标实际区域;解决了现有技术中断层图像的分割方法分割精度较低的问题,进而达到了提高断层图像分割精确度的效果。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的分割方法的流程图;
图2是根据本发明实施例的分割方法得到的分别目标在各个分割层面上的分割结果示意图;
图3是根据本发明实施例的分割方法得到的各层序列点的示意图;
图4是根据本发明优选实施例的分割方法的流程图;以及
图5是根据本发明实施例的分割装置的示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
本发明实施例提供了一种断层图像的分割方法,以下对本发明实施例所提供的断层图像的分割方法进行具体介绍:
图1是根据本发明实施例的断层图像的分割方法的流程图,如图1所示,该实施例所提供的分割方法包括如下的步骤S101至步骤S108:
S101:接收来自用户的第一选择指令,并从断层图像中获取与第一选择指令相对应的分割目标;即,由用户确定感兴趣的器官或解剖结构,并下发相应的选择指令以从二维断层图像中确定出与选择指令相对应的分割目标。
S102:接收来自用户的第二选择指令并从断层图像中确定与第二选择指令相对应的层面,得到第一人机交互分割层面;具体地,用户可以从步骤S101确定出的分割目标的二维断层图像所构成的三维体数据中逐层浏览层面图像,然后指定分割目标开始出现的初始层面作为第一人机交互分割层面。
S103:对第一人机交互分割层面进行分割,得到分割目标在第一人机交互分割层面内的分割结果;具体地,可以采用多种现有人机交互方法包括但不限于区域生长、水平集(Level Set)和图割(Graph Cut)等方法对第一人机交互分割层面进行分割,类似上述所列出的人机交互方法的共同特点是要求用户指定或输入一定的初始条件,如种子点、起始轮廓等,并且允许用户对分割结果进行手工修正,以确保分割结果准确。
S104:接收来自用户的第三选择指令并从断层图像中确定与第三选择指令相对应的层面,得到第二人机交互分割层面;具体地,类似步骤S102,用户可以从步骤S101确定出的分割目标的二维断层图像所构成的三维体数据中逐层浏览层面图像,然后指定分割目标最后出现的终止层面作为第二人机交互分割层面。
S105:对第二人机交互分割层面进行分割,得到分割目标在第二人机交互分割层面内的分割结果;具体地,可以采用与步骤S103中的人机交互方法对第二人机交互分割层面进行分割。
S106:获取第一人机交互分割层面与第二人机交互分割层面之间的中间层面,得到第一自动分割层面;即,将初始层面和终止层面之间的所有的中间层面作为自动分割层面,该第一自动分割层面可能为一个层面,也可能为多个层面,一般情况下,初始层面和终止层面之间包括多个中间层面的情况居多。
S107:对第一自动分割层面进行分割,得到分割目标在第一自动分割层面内的分割结果;即,对所有的中间层面均进行分割,得到分割目标在各个中间层面上的分割结果。
S108:将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果和分割目标在第一自动分割层面内的分割结果进行组合,得到分割目标的三维分割结果。即,将分割目标在各个层面内的二维分割结果进行组合,生成分割目标的三维分割结果。其中,图2中示意性示出了分割目标在起始层面、多个中间层面和终止层面上的二维分割结果。
本发明实施例所提供的断层图像的分割方法通过根据用户的选择指令分别确定出两个人机交互分割层面,并对确定出的这两个人机交互分割层面进行人机交互分割,同时还对这两个人机交互分割层面之间的中间层面进行自动分割,实现了允许用户指定一个以上的人机交互分割层面,避免了现有技术中只允许存在一个人机交互分割层面所带来的以下弊端:在“逐层传递分割”过程中,受各种可能干扰因素的影响,分割算法在某一层面的分割结果有可能发生错误,偏离分割目标实际区域;解决了现有技术中断层图像的分割方法分割精度较低的问题,进而达到了提高断层图像分割精确度的效果。
优选地,步骤S107中可以采用以分割目标在第一人机交互分割层面内的分割结果对应的分割轮廓以及分割目标在第二人机交互分割层面内的分割结果对应的分割轮廓为基础,利用层间轮廓自动插值法对自动分割层面进行分割,其中,假设第一分割轮廓为分割目标在第一人机交互分割层面内的分割结果对应的分割轮廓,第二分割轮廓为分割目标在第二人机交互分割层面内的分割结果对应的分割轮廓,层间轮廓自动插值法的具体分割步骤为:
a.对第一分割轮廓的曲线长度进行预设份数等分,得到多个第一等分点;即,对层面P中的分割轮廓,按照轮廓的曲线长度进行等分,得到等分点在层面P中的二维轮廓点序列Pk=(PXk,PYk),k=0,…,N-1。其中,N为固定常数,在本发明实施例的分割方法中可以取100,其中,Pk=(PXk,PYk)称作第一序列,对于k=0,…,N-1,从P(0)到P(99)的曲线长度相等。
b.采用与步骤a相同的方法,对第二人机交互分割层面Q中的分割轮廓进行等分划分,并得到第二序列Qk=(QXk,QYk),k=0,…,N-1。N采用与层面P中相同的固定常数。
c.确定第一序列中各个序列点Pk=(PXk,PYk)与第二序列中各个序列点Qk=(QXk,QYk)的对应关系,其中,可以采用总体距离平方和最小的原则来确定这两组序列点的对应关系,具体地,
首先,计算第一分割轮廓区域的质心,得到第一质心,具体地,可以通过在每个层面图像中定义二维坐标系XY,该坐标的坐标原点可以取层面图像区域左下角,向右为X轴正方向,向上为Y轴正方向,则层面中每一个像素点的二维坐标可以写成(x,y),分割目标在第一人机交互分割层面内的分割结果上的各个像素点也均可以通过二维坐标来表示,对第一分割轮廓区域的质心的求取则可以通过如下方法进行:假定第一人机交互分割层面P中分割区域内所有像素点的坐标集合e(i)=(eXi,eYi),i=1,…,M。区域质心PC=(PXC,PYC)定义为PXC=sum(eXi)/M,PYC=sum(eYi)/M,sum(*)表示求和,i=1,…,M。并采用相同的原理计算第二分割轮廓区域的质心,得到第二质心QC=(QCX,QCY)。
其次,通过以下公式确定第一序列中各个序列点与第二序列中各个序列点的对应关系:
D ( i ) = Σ k = 0 N - 1 | P ( mod ( k + i , N ) ) - P C , Qk - Q C | 2 , i = 0,1,2 , . . . N - 1 i ′ = min i ( D ( i ) ) ,
其中,mod(k+i,N)表示取(k+i)被N除的余数,P(mod(k+i,N))=(PX(mod(k+i,N)),PY(mod(k+i,N))表示第一序列中各个序列点在第一人机交互分割层面中的图像坐标,PC=(PXC,PYC)表示第一质心在第一人机交互分割层面中的图像坐标,P(mod(k+i,N))-PC=(PX(mod(k+i,N))-PXC,PY(mod(k+i,N))-PYC)表示P(mod(k+i,N))相对PC的坐标,Qk=(QX(k),QY(k))表示所述第二序列中各个序列点在第二人机交互分割层面中的图像坐标,QC=(QCX,QCY)表示第二质心在第二人机交互分割层面中的图像坐标,Qk-QC=(QX(k)-QXC,QY(k)-QYC)表示Qk相对QC的坐标,|A,B|表示坐标点A和B之间的距离,N为所述预设份数,对于包括CT、MR、PT在内的断层影像N的经验取值为100,i′=mini(D(i))表示当i=i'时,D(i)的取值最小,第一序列中的序列点P(mod(k+i′,N))与第二序列中的序列点Qk相对应,举例说明,当i′=0时,P0和Q0对应、P1和Q1对应、……、P(N-1)和Q(N-1)对应;当i′=1时,P1和Q0对应、P2和Q1对应、……、P(N-1)和Q(N-2)对应、P0和Q(N-1)对应。
d.采用预设插值公式按照步骤c确定出的对应关系确定分割目标在第一层面L内的分割轮廓,得到分割目标在第一层面内的分割结果,其中,假设第一层面L为第一自动分割层面中的任意一个层面,对自动分割层面中每一个层面对应的分割轮廓的确定方式均与分割目标在层面L内的分割轮廓的确定方式相同,具体地,主要通过确定分割目标在第一层面L内的分割轮廓上的序列点来确定在层面L上的分割结果,对分割目标在第一层面L内的分割轮廓上的序列点的确定可以通过以下两种方式进行:
方式一:直接按照下面示出的第一预设差值公式来计算层面L内的分割轮廓上的序列点Lk=(LXk,LYk),第一预设差值公式:
L ( k ) = LX ( k ) = QX ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ PX ( mod ( k + i ′ , N ) ) - QX ( k ) ] LY ( k ) = QY ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ PY ( mod ( k + i ′ , N ) ) - QY ( k ) ] , k = 0,1,2 , . . . , N - 1
方式二:结合步骤a和c对第一序列进行重新编号,令P(mod(k+i,N))的编号为P′k,,此种重新编号的方法不失一般性,举例说明,假定i′=3,N=100,则将P3重新编号为P′0,P4编为P′1……,P99编为P′96,P0编为P′97,P1编为P′98,P2编为P′99经过重新编号后,序列P′k和序列Qk按照标号顺序,逐点对应,按照下面示出的第二预设差值公式来计算层面L内的分割轮廓上的序列点Lk=(LXk,LYk),第二预设差值公式:
L ( k ) = LXk = QXk + [ D 2 / ( D 1 + D 2 ) ] × [ P ′ Xk - QXk ] LYk = QYk + [ D 2 / ( D 1 + D 2 ) ] × [ P ′ Yk - QYk ] , k = 0,1 , . . . , N - 1
其中,第一预设差值公式中的D1与第二预设差值公式中的D1含义相同,第一预设差值公式中的D2与第二预设差值公式中的D2含义相同,D1表示第一层面距离第一人机交互分割层面的垂直距离,D2表示第一层面距离第二人机交互分割层面的垂直距离。采用本发明实施例的层间轮廓自动插值方法得到的各层序列点在图3中示出。
通过采用上述的“层间轮廓自动插值方法”,能够保证自动分割层面中的分割轮廓,从一个人机交互分割层面中的分割轮廓,向另一个人机交互分割层面中分割轮廓平滑过渡,避免了现有“层间传递”方法存在的偏差放大,或分割失败可能。进而达到了进一步地提高断层图像的分割精确度。
图4是根据本发明优选实施例的分割方法的流程图,如图4所示,该优选实施例的分割方法与图1中示出的分割方法相比,图4中示出的优选实施例的分割方法还包括根据用户的选择指令选择分割目标的层面方向,以及对分割目标在各个自动分割层面上的分割结果进行检查以筛选出偏差较大的分割层面进行再次分割,其中,对分割目标在各个自动分割层面上的分割结果进行检查以筛选出偏差较大的分割层面进行再次分割具体地为:
第一、接收来自用户的第四选择指令并从多个分割结果中确定与第四选择指令相对应的分割结果,其中,第四选择执行用于选择多个分割结果中偏差最大的分割结果,也即,对各个自动分割层面上的分割结果进行检查的步骤可以由用户自身进行判断,以准确找出偏差最大的分割结果;
第二、获取与第四选择指令相对应的分割结果所对应的分割层面,得到第三人机交互分割层面,即,获取偏差最大的分割结果对应的分割层面,并将该获取到的分割层面作为第三人机交互分割层面;
第三、对第三人机交互分割层面进行分割,得到分割目标在第三人机交互分割层面内的分割结果;具体地,可以采用与步骤S103中的人机交互方法对第三人机交互分割层面进行分割。
第四、获取第一人机交互分割层面与第三人机交互分割层面之间的中间层面,得到第二自动分割层面;即,将初始层面和第三人机交互分割层面之间的所有的中间层面作为自动分割层面,该第二自动分割层面可能为一个层面,也可能为多个层面,一般情况下,初始层面和第三人机交互分割层面之间包括多个中间层面的情况居多。
第五、对第二自动分割层面进行分割,得到分割目标第二自动分割层面内的分割结果;具体地,可以采用层间轮廓自动插值法对初始层面和第三人机交互分割层面之间的所有中间层面均进行分割,得到分割目标在各个中间层面上的分割结果。
第六、获取第二人机交互分割层面与第三人机交互分割层面之间的中间层面,得到第三自动分割层面;即,将终止层面和第三人机交互分割层面之间的所有的中间层面作为自动分割层面,该第三自动分割层面可能为一个层面,也可能为多个层面,一般情况下,终止层面和第三人机交互分割层面之间包括多个中间层面的情况居多。
第七、对第三自动分割层面进行分割,得到分割目标在第三自动分割层面内的分割结果,具体地,可以采用层间轮廓自动插值法对终止层面和第三人机交互分割层面之间的所有中间层面均进行分割,得到分割目标在各个中间层面上的分割结果。
第八、将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果、分割目标在第三人机交互分割层面内的分割结果、分割目标在第二自动分割层面内的分割结果和分割目标在第三自动分割层面内的分割结果进行组合,得到三维分割结果。也即,当确定出偏差最大的自动分割层面后,将确定出的该自动分割层面作为一个新的人机交互分割层面,以采用人机交互分割方法对该层面进行分割,并依次采用层间轮廓自动插值法对这个新的人机交互分割层面与初始层面之间的中间层面、与终止层面之间的中间层面进行分割,并将最终的分割结果进行合并后替换原有的三维分割结果。其中,为保证三维分割结果的精度,还可以在对原有的三维分割结果进行更新替换后,再返回第一步骤,以对各个自动分割层面的分割结果进行检查验证,保证三维分割结果的精确度。
本发明优选实施例的分割方法通过对自动分割层面中结果偏差最大的分割层面重新采用人机交互分割方法进行分割,实现了允许用户将自动分割的层面指定为新的人机交互分割层面,并基于新增加的人机交互分割层面,更新其余中间层面自动分割结果,以对结果偏差较大的分割层面进行人工修正;同时通过对新产生的人机交互分割层面与位于其上方和下方的原有的人机交互分割层面之间层面进行层间轮廓自动插值法分割,实现了在进行偏差修正后从一个人机交互分割层面中的分割轮廓向另一个人机交互分割层面中分割轮廓平滑过渡,达到了进一步地提高断层图像的分割精确度。同时,通过根据用户的选择指令从断层图像中确定层面方向,达到了选择有利于观察和分割的层面方向。
本发明实施例还提供了一种断层图像的分割装置,以下对本发明实施例所提供的分割装置进行具体介绍:
图5是根据本发明实施例的断层图像分割装置的示意图,如图5所示,该实施例的分割装置包括:第一接收单元10、第二接收单元20、第一分割单元30、第三接收单元40、第二分割单元50、第一获取单元60、第三分割单元70和组合单元80。
具体地,第一接收单元10用于接收来自用户的第一选择指令,并从断层图像中获取与第一选择指令相对应的分割目标;即,由用户确定感兴趣的器官或解剖结构,并通过第一接收单元10上的操作按键或触摸屏等下发相应的选择指令以从二维断层图像中确定出与选择指令相对应的分割目标。
第二接收单元20用于接收来自用户的第二选择指令并从断层图像中确定与第二选择指令相对应的层面,得到第一人机交互分割层面;具体地,用户可以从上述确定出的分割目标的二维断层图像所构成的三维体数据中逐层浏览层面图像,然后通过第二接收单元上的操作按键或触摸屏等下发指令以指定分割目标开始出现的初始层面作为第一人机交互分割层面。
第一分割单元30用于对第一人机交互分割层面进行分割,得到分割目标在第一人机交互分割层面内的分割结果;具体地,可以采用多种现有人机交互方法包括但不限于区域生长、水平集(Level Set)和图割(Graph Cut)等方法对第一人机交互分割层面进行分割,类似上述所列出的人机交互方法的共同特点是要求用户指定或输入一定的初始条件,如种子点、起始轮廓等,并且允许用户对分割结果进行手工修正,以确保分割结果准确。
第三接收单元40用于接收来自用户的第三选择指令并从断层图像中确定与第三选择指令相对应的层面,得到第二人机交互分割层面;具体地,用户可以从上述确定出的分割目标的二维断层图像所构成的三维体数据中逐层浏览层面图像,然后通过第三接收单元上的操作按键或触摸屏等下发指令以指定分割目标开始出现的初始层面作为第二人机交互分割层面。
第二分割单元50用于对第二人机交互分割层面进行分割,得到分割目标在第二人机交互分割层面内的分割结果;具体地,第二分割单元50可以采用与第一分割单元30相同的分割方法。
第一获取单元60用于获取第一人机交互分割层面与第二人机交互分割层面之间的中间层面,得到第一自动分割层面;即,将初始层面和终止层面之间的所有的中间层面作为自动分割层面,该第一自动分割层面可能为一个层面,也可能为多个层面,一般情况下,初始层面和终止层面之间包括多个中间层面的情况居多。
第三分割单元70用于对第一自动分割层面进行分割,得到分割目标在第一自动分割层面内的分割结果;即,对所有的中间层面均进行分割,得到分割目标在各个中间层面上的分割结果。
组合单元80用于将分割目标在第一人机交互分割层面内的分割结果、分割目标在第二人机交互分割层面内的分割结果和分割目标在第一自动分割层面内的分割结果进行组合,得到分割目标的三维分割结果。即,将分割目标在各个层面内的二维分割结果进行组合,生成分割目标的三维分割结果。
本发明实施例所提供的断层图像的分割装置通过根据用户的选择指令分别确定出两个人机交互分割层面,并对确定出的这两个人机交互分割层面进行人机交互分割,同时还对这两个人机交互分割层面之间的中间层面进行分割,实现了允许用户指定一个以上的人机交互分割层面,避免了现有技术中只允许存在一个人机交互分割层面所带来的以下弊端:在“逐层传递分割”过程中,受各种可能干扰因素的影响,分割算法在某一层面的分割结果有可能发生错误,偏离分割目标实际区域;解决了现有技术中断层图像的分割方法分割精度较低的问题,进而达到了提高断层图像分割精确度的效果。
优选地,第三分割单元70包括获取子单元和分割子单元,其中,获取子单元用于获取分割目标在第一人机交互分割层面内的分割结果对应的分割轮廓以及分割目标在第二人机交互分割层面内的分割结果对应的分割轮廓,分割子单元用于以获取子单元获取到的分割轮廓为基础,利用层间轮廓自动插值法对自动分割层面进行分割,其中,假设第一分割轮廓为分割目标在第一人机交互分割层面内的分割结果对应的分割轮廓,第二分割轮廓为分割目标在第二人机交互分割层面内的分割结果对应的分割轮廓,分割子单元可以通过以下各模块采用层间轮廓自动插值法对自动分割层面进行分割:由第一划分模块对第一分割轮廓的曲线长度进行预设份数等分,得到多个第一等分点;由第一获取模块分别获取多个第一等分点在第一人机交互分割层面中的图像坐标,得到第一序列;由第二划分模块对第二分割轮廓的曲线长度进行预设份数等分,得到多个第二等分点;由第二获取模块分别获取多个第二等分点在第二人机交互分割层面中的图像坐标,得到第二序列;由第一确定模块确定第一序列中各个序列点与第二序列中各个序列点的对应关系;以及由第二确定模块采用预设插值公式按照确定出的对应关系确定分割目标在第一层面内的分割轮廓,得到分割目标在第一层面内的分割结果,其中,第一层面为第一自动分割层面中的任一层面。
通过上述各模块采用“层间轮廓自动插值方法”对自动分割层面进行分割,能够保证自动分割层面中的分割轮廓,从一个人机交互分割层面中的分割轮廓,向另一个人机交互分割层面中分割轮廓平滑过渡,避免了现有“层间传递”方法存在的偏差放大,或分割失败可能。进而达到了进一步地提高断层图像的分割精确度。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种断层图像的分割方法,其特征在于,包括:
接收来自用户的第一选择指令,并从断层图像中获取与所述第一选择指令相对应的分割目标;
接收来自所述用户的第二选择指令并从所述断层图像中确定与所述第二选择指令相对应的层面,得到第一人机交互分割层面;
对所述第一人机交互分割层面进行分割,得到所述分割目标在所述第一人机交互分割层面内的分割结果;
接收来自所述用户的第三选择指令并从所述断层图像中确定与所述第三选择指令相对应的层面,得到第二人机交互分割层面;
对所述第二人机交互分割层面进行分割,得到所述分割目标在所述第二人机交互分割层面内的分割结果;
获取所述第一人机交互分割层面与所述第二人机交互分割层面之间的中间层面,得到第一自动分割层面;
对所述第一自动分割层面进行分割,得到所述分割目标在所述第一自动分割层面内的分割结果;以及
将所述分割目标在所述第一人机交互分割层面内的分割结果、所述分割目标在所述第二人机交互分割层面内的分割结果和所述分割目标在所述第一自动分割层面内的分割结果进行组合,得到所述分割目标的三维分割结果。
2.根据权利要求1所述的分割方法,其特征在于,对所述第一自动分割层面进行分割包括:
获取第一分割轮廓和第二分割轮廓,其中,所述第一分割轮廓为所述分割目标在所述第一人机交互分割层面内的分割结果对应的分割轮廓,所述第二分割轮廓为所述分割目标在所述第二人机交互分割层面内的分割结果对应的分割轮廓;以及
以所述第一分割轮廓和所述第二分割轮廓为基础采用层间轮廓自动插值法对所述第一自动分割层面进行分割。
3.根据权利要求2所述的分割方法,其特征在于,所述第一自动分割层面包括一个或多个分割层面,以所述第一分割轮廓和所述第二分割轮廓为基础采用层间轮廓自动插值法对所述第一自动分割层面进行分割包括:
对所述第一分割轮廓的曲线长度进行预设份数等分,得到多个第一等分点;
分别获取所述多个第一等分点在所述第一人机交互分割层面中的图像坐标,得到第一序列;
对所述第二分割轮廓的曲线长度进行预设份数等分,得到多个第二等分点;
分别获取所述多个第二等分点在所述第二人机交互分割层面中的图像坐标,得到第二序列;
确定所述第一序列中各个序列点与所述第二序列中各个序列点的对应关系;以及
采用预设插值公式按照确定出的对应关系确定所述分割目标在第一层面内的分割轮廓,得到所述分割目标在所述第一层面内的分割结果,其中,所述第一层面为所述第一自动分割层面中的任一层面。
4.根据权利要求3所述的分割方法,其特征在于,确定所述第一序列中各个序列点与所述第二序列中各个序列点的对应关系包括:
计算所述第一分割轮廓区域的质心,得到第一质心;
计算所述第二分割轮廓区域的质心,得到第二质心;以及
通过以下公式确定所述第一序列中各个序列点与所述第二序列中各个序列点的对应关系:
D ( i ) = Σ k = 0 N - 1 | P ( mod ( k + i , N ) ) - P C , Qk - Q C | 2 , i = 0,1,2 , . . . N - 1 i ′ = min i ( D ( i ) ) ,
其中,mod(k+i,N)表示取(k+i)被N除的余数,P(mod(k+i,N))=(PX(mod(k+i,N)),PY(mod(k+i,N)))表示所述第一序列中各个序列点的图像坐标,PC=(PXC,PYC)表示所述第一质心的图像坐标,P(mod(k+i,N))-PC=(PX(mod(k+i,N))-PXC,PY(mod(k+i,N))-PYC)表示P(mod(k+i,N))相对PC的坐标,Qk=(QX(k),QY(k))表示所述第二序列中各个序列点的图像坐标,QC=(QCX,QCY)表示所述第二质心的图像坐标,Qk-QC=(QX(k)-QXC,QY(k)-QYC)表示Qk相对QC的坐标,|A,B|表示坐标点A和B之间的距离,N为所述预设份数,i′=mini(D(i))表示i=i'时D(i)取值最小,所述第一序列中的序列点P(mod(k+i′,N))与所述第二序列中的序列点Qk相对应,
所述预设插值公式为:
L ( k ) = LX ( k ) = QX ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ PX ( mod ( k + i ′ , N ) ) - QX ( k ) ] LY ( k ) = QY ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ PY ( mod ( k + i ′ , N ) ) - QY ( k ) ] , k = 0,1,2 , . . . , N - 1
其中,L(k)=(LX(k),LY(k))表示所述分割目标在所述第一层面内的分割轮廓上的序列点,D1表示所述第一层面距离所述第一人机交互分割层面的垂直距离,D2表示所述第一层面距离所述第二人机交互分割层面的垂直距离。
5.根据权利要求3所述的分割方法,其特征在于,确定所述第一序列中各个序列点与所述第二序列中各个序列点的对应关系包括:
计算所述第一分割轮廓区域的质心,得到第一质心;
计算所述第二分割轮廓区域的质心,得到第二质心;以及
通过以下公式确定所述第一序列中各个序列点与所述第二序列中各个序列点的对应关系:
D ( i ) = Σ k = 0 N - 1 | P ( mod ( k + i , N ) ) - P C , Qk - Q C | 2 , i = 0,1,2 , . . . N - 1 i ′ = min i ( D ( i ) ) ,
其中,mod(k+i,N)表示取(k+i)被N除的余数,P(mod(k+i,N))=(PX(mod(k+i,N)),PY(mod(k+i,N)))表示所述第一序列中各个序列点的图像坐标,PC=(PXC,PYC)表示所述第一质心的图像坐标,P(mod(k+i,N))-PC=(PX(mod(k+i,N))-PXC,PY(mod(k+i,N))-PYC)表示P(mod(k+i,N))相对PC的坐标,Qk=(QX(k),QY(k))表示所述第二序列中各个序列点的图像坐标,QC=(QCX,QCY)表示所述第二质心的图像坐标,Qk-QC=(QX(k)-QXC,QY(k)-QYC)表示Qk相对QC的坐标,|A,B|表示坐标点A和B之间的距离,N为所述预设份数,i′=mini(D(i))表示i=i'时D(i)取值最小,所述第一序列中的序列点P(mod(k+i′,N))与所述第二序列中的序列点Qk相对应,
在确定所述第一序列中各个序列点与所述第二序列中各个序列点的对应关系之后,所述分割方法还包括:采用公式P′k=P(mod(k+i′,N))对所述第一序列中的各个序列点进行重新编号,得到所述第一序列中的各个序列点P′k,
所述预设插值公式为:
L ( k ) = LX ( k ) = QX ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ P ′ X ( k ) - QX ( k ) ] LY ( k ) = QY ( k ) + [ D 2 / ( D 1 + D 2 ) ] × [ P ′ Y ( k ) - QY ( k ) ] , k = 0,1,2 , . . . , N - 1
其中,L(k)=(LX(k),LY(k))表示所述分割目标在所述第一层面内的分割轮廓上的序列点,D1表示所述第一层面距离所述第一人机交互分割层面的垂直距离,D2表示所述第一层面距离所述第二人机交互分割层面的垂直距离。
6.根据权利要求1所述的分割方法,其特征在于,所述第一自动分割层面包括多个分割层面,所述分割目标在所述第一自动分割层面内的分割结果包括多个分割结果,在得到所述分割目标在所述第一自动分割层面内的分割结果之后,所述分割方法还包括:
接收来自所述用户的第四选择指令并从所述多个分割结果中确定与所述第四选择指令相对应的分割结果,其中,所述第四选择指令用于选择所述多个分割结果中偏差最大的分割结果;
获取与所述第四选择指令相对应的分割结果所对应的分割层面,得到第三人机交互分割层面;
对所述第三人机交互分割层面进行分割,得到所述分割目标在所述第三人机交互分割层面内的分割结果;
获取所述第一人机交互分割层面与所述第三人机交互分割层面之间的中间层面,得到第二自动分割层面;
对所述第二自动分割层面进行分割,得到所述分割目标在所述第二自动分割层面内的分割结果;
获取所述第二人机交互分割层面与所述第三人机交互分割层面之间的中间层面,得到第三自动分割层面;以及
对所述第三自动分割层面进行分割,得到所述分割目标在所述第三自动分割层面内的分割结果,
其中,将所述分割目标在所述第一人机交互分割层面内的分割结果、所述分割目标在所述第二人机交互分割层面内的分割结果、所述分割目标在所述第三人机交互分割层面内的分割结果、所述分割目标在所述第二自动分割层面内的分割结果和所述分割目标在所述第三自动分割层面内的分割结果进行组合,得到所述三维分割结果。
7.根据权利要求6所述的分割方法,其特征在于,所述分割目标在所述第二自动分割层面内的分割结果以及所述分割目标在所述第三自动分割层面内的分割结果均包括多个分割结果,在得到所述分割目标在所述第二自动分割层面内的分割结果,以及得到所述分割目标在所述第三自动分割层面内的分割结果之后,所述分割方法还包括:
重新接收所述第四选择指令以分别从所述第二自动分割层面内的多个分割结果和所述第三自动分割层面内的多个分割结果中确定与所述第四选择指令相对应的分割结果。
8.根据权利要求1所述的分割方法,其特征在于,在从断层图像中获取与所述第一选择指令相对应的分割目标之后,所述分割方法还包括:
接收来自所述用户的第五选择指令并从所述断层图像中确定与所述第五选择指令相对应的层面方向,
其中,所述第一人机交互分割层面、所述第二人机交互分割层面和所述第一自动分割层面均为沿所述层面方向上的分割层面。
9.一种断层图像的分割装置,其特征在于,包括:
第一接收单元,用于接收来自用户的第一选择指令,并从断层图像中获取与所述第一选择指令相对应的分割目标;
第二接收单元,用于接收来自所述用户的第二选择指令并从所述断层图像中确定与所述第二选择指令相对应的层面,得到第一人机交互分割层面;
第一分割单元,用于对所述第一人机交互分割层面进行分割,得到所述分割目标在所述第一人机交互分割层面内的分割结果;
第三接收单元,用于接收来自所述用户的第三选择指令并从所述断层图像中确定与所述第三选择指令相对应的层面,得到第二人机交互分割层面;
第二分割单元,用于对所述第二人机交互分割层面进行分割,得到所述分割目标在所述第二人机交互分割层面内的分割结果;
第一获取单元,用于获取所述第一人机交互分割层面与所述第二人机交互分割层面之间的中间层面,得到第一自动分割层面;
第三分割单元,用于对所述第一自动分割层面进行分割,得到所述分割目标在所述第一自动分割层面内的分割结果;以及
组合单元,用于将所述分割目标在所述第一人机交互分割层面内的分割结果、所述分割目标在所述第二人机交互分割层面内的分割结果和所述分割目标在所述第一自动分割层面内的分割结果进行组合,得到所述分割目标的三维分割结果。
10.根据权利要求9所述的分割装置,其特征在于,所述第三分割单元包括:
获取子单元,用于获取第一分割轮廓和第二分割轮廓,其中,所述第一分割轮廓为所述分割目标在所述第一人机交互分割层面内的分割结果对应的分割轮廓,所述第二分割轮廓为所述分割目标在所述第二人机交互分割层面内的分割结果对应的分割轮廓;以及
分割子单元,用于以所述第一分割轮廓和所述第二分割轮廓为基础采用层间轮廓自动插值法对所述第一自动分割层面进行分割。
CN201210497386.2A 2012-11-28 2012-11-28 断层图像的分割方法和装置 Active CN102982546B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210497386.2A CN102982546B (zh) 2012-11-28 2012-11-28 断层图像的分割方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210497386.2A CN102982546B (zh) 2012-11-28 2012-11-28 断层图像的分割方法和装置

Publications (2)

Publication Number Publication Date
CN102982546A true CN102982546A (zh) 2013-03-20
CN102982546B CN102982546B (zh) 2015-04-01

Family

ID=47856487

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210497386.2A Active CN102982546B (zh) 2012-11-28 2012-11-28 断层图像的分割方法和装置

Country Status (1)

Country Link
CN (1) CN102982546B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105787978A (zh) * 2016-02-29 2016-07-20 深圳市医诺智能科技发展有限公司 一种医学图像隔层自动勾画方法、装置和系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040202369A1 (en) * 2002-12-06 2004-10-14 Nikolaos Paragios User interactive level set methods for image segmentation
CN102622750A (zh) * 2012-02-24 2012-08-01 西安电子科技大学 基于交互式区域生长的胃部ct序列图像分割方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040202369A1 (en) * 2002-12-06 2004-10-14 Nikolaos Paragios User interactive level set methods for image segmentation
CN102622750A (zh) * 2012-02-24 2012-08-01 西安电子科技大学 基于交互式区域生长的胃部ct序列图像分割方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ANDREA SCHENK ET AL: "《Efficient Semiautomatic Segmentation of 3D Objects in Medical Images》", 《MEDICAL IMAGE COMPUTING AND COMPUTER ASSISTED INTERVENTION-MICCAI 2000》 *
徐伟等: "《一种快速有效的交互式医学序列图像分割方法》", 《现代电子技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105787978A (zh) * 2016-02-29 2016-07-20 深圳市医诺智能科技发展有限公司 一种医学图像隔层自动勾画方法、装置和系统

Also Published As

Publication number Publication date
CN102982546B (zh) 2015-04-01

Similar Documents

Publication Publication Date Title
RU2562898C2 (ru) Определение представляющей интерес области при визуализации сердца
US20210361174A1 (en) System for vascular assessment
US6999811B2 (en) Method and device for the registration of two 3D image data sets
US9554728B2 (en) Visualization for navigation guidance
EP2622571B1 (en) Advanced fiber tracking and medical navigation in a brain
US7986823B2 (en) System and method for consistent detection of mid-sagittal planes for magnetic resonance brain scans
EP2558968B1 (en) System and method for planning a neurosurgical operation
CN105408939A (zh) 用于配准成像设备与跟踪设备的配准系统
CN103608842A (zh) 用于处理医学图像的系统和方法
US10311585B2 (en) Mechanism for advanced structure generation and editing
CN103140855A (zh) 基于知识的自动图像分割
US20140023250A1 (en) Medical imaging systen and method for providing an image representation supporting accurate guidance of an intervention device in a vessel intervention procedure
US20210353361A1 (en) Surgical planning, surgical navigation and imaging system
WO2011019456A1 (en) Ct atlas of musculoskeletal anatomy to guide treatment of sarcoma
EP2584970A2 (en) Mechanism for advanced structure generation and editing
US10325412B2 (en) Cutting three-dimensional image
WO2011041475A1 (en) Medical image analysis system using n-way belief propagation for anatomical images subject to deformation and related methods
US20160300398A1 (en) System for dynamically propagating real-time alterations of medical images
CN102982546B (zh) 断层图像的分割方法和装置
CN103518225B (zh) 心肌聚类与取向
Andrea et al. Validation of stereo vision based liver surface reconstruction for image guided surgery
CN115908479A (zh) 一种脑组织漂变校正方法和手术导航系统
CN108701492A (zh) 医学图像导航系统
US20240164851A1 (en) Systems and methods for a multidimensional tracking system
US20140309476A1 (en) Ct atlas of musculoskeletal anatomy to guide treatment of sarcoma

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant