CN101268493A - 用于手动调整几何模型的先验信息编码 - Google Patents

用于手动调整几何模型的先验信息编码 Download PDF

Info

Publication number
CN101268493A
CN101268493A CNA200680034615XA CN200680034615A CN101268493A CN 101268493 A CN101268493 A CN 101268493A CN A200680034615X A CNA200680034615X A CN A200680034615XA CN 200680034615 A CN200680034615 A CN 200680034615A CN 101268493 A CN101268493 A CN 101268493A
Authority
CN
China
Prior art keywords
geometric model
adjustment
zone
instrument
view data
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.)
Pending
Application number
CNA200680034615XA
Other languages
English (en)
Inventor
J·威斯
O·艾卡伯特
J·彼得斯
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN101268493A publication Critical patent/CN101268493A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/149Segmentation; Edge detection involving deformable models, e.g. active contour models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20Indexing scheme for editing of 3D models
    • G06T2219/2021Shape modification

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Architecture (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明涉及调整方法(200),用于调整几何模型以适应图像数据,所述调整方法包括:区域选择步骤(230),用于选择所述几何模型的区域;以及手动调整步骤(235),用于利用所述几何模型中所包含的所述区域的一组特征,来手动调整所述几何模型以适应所述图像数据。使用来自选定区域的特征组中的区域特有特征,例如确定变形范围的区域特有参数,来手动调整所述几何模型以适应所述图像数据,需要相对较少的用户交互。

Description

用于手动调整几何模型的先验信息编码
本发明涉及一种调整方法,用于调整几何模型以适应图像数据。
本发明还涉及一种调整系统,用于调整几何模型以适应图像数据。
本发明还涉及一种采集系统,用于获取图像数据,其包括所述调整系统。
本发明还涉及一种工作站,其包括所述调整系统。
本发明还涉及一种由计算机设备加载的计算机程序产品,其包括指令,用于调整几何模型以适应图像数据。
在WO2005/038711中,描述了开始段落中所述的这种调整方法的实施例,以下将其称为参考文献1。该文献描述了几种用于修改几何模型的手动工具,例如高斯拖拉工具和球体推动工具,其目的是为了改善调整几何模型以适应图像数据的结果。例如,在该调整方法的一个实施例中,用户可以选择表示几何模型的网格的一个顶点,并利用高斯拖拉工具来将它拖拉至期望的位置。还可以移动周围的顶点。用平滑函数来控制它们的位移,该平滑函数例如为以该选定顶点为中心的高斯函数。
该方法的缺陷在于几何模型的变形由工具的几何参数来控制。例如,利用诸如高斯函数的半径之类的参数来控制由高斯拖拉工具引起的变形,例如,该高斯函数对被定义为高斯函数的方差的平方根的变形进行建模。将几何模型的一个顶点拖拉向图像数据内的边界会导致几何模型的变形,其没有将几何模型的边界与图像数据内的目标边界对准。最后,可以通过操纵多个顶点直至完成该手动调整为止,来校正在几何模型的边界对准中的所有误差。可替换地,可由用户交互地设置工具的参数,来优化该工具。然而,任一种解决方案都需要大量的用户交互且因此会非常耗费时间。
本发明的目的是提供开始段落中所描述的那种需要相对较少用户交互的调整方法。
本发明的目的是通过以下实现的:用于调整几何模型以适应图像数据的调整方法包括:
区域选择步骤,用于选择所述几何模型的区域;以及
手动调整步骤,用于利用所述几何模型内所包含的所述区域的一组特征,来手动调整所述几何模型以适应所述图像数据。
使用一种利用鼠标操作的标准选择工具来选择所述几何模型中将要被手动调整以适应图像数据的几何模型的区域。同样的工具用于手动调整几何模型,以便将其与图像数据内所包含的察觉到的对象对准。根据本发明所述,由选定区域的一组特征来控制几何模型的变形。该组特征包含在该几何模型内。这允许将优化的工具参数应用于几何模型的选定区域。例如,在由包含多个顶点的自适应网格所表示的几何模型情况下,网格的每一个顶点都是几何模型的一个区域。每一个顶点都可以用其自己的权重函数来表征,权重函数用于控制由应用于该区域的拖拉工具所引起的变形。所述权重函数包括网格中所有顶点的权重,它们中的大部分通常为零。当将拖拉工具应用于选定顶点时,网格中每一个顶点的位移都与选定顶点的权重成正比。选定顶点的权重为1。例如,在某个特定顶点上的权重函数的值可以按指数规律减少。该指数可以与从该特定顶点到选定顶点之间的距离成正比。这样的权重函数用随着与选定区域之间的距离而下降的速率来表征。与选定顶点相关的权重函数的下降速率确定了通过将拖拉工具应用于该顶点而引起的变形范围。例如,如果在选定顶点周围区域中的顶点彼此强耦合,例如表示骨头平坦区域的顶点,则拉工具的下降速率相对较小。利用拖拉工具从该区域中拖拉一个顶点将导致附近顶点的类似位移。另一方面,如果在选定顶点周围区域中的顶点彼此弱耦合,例如表示弯曲软组织的顶点,则权重函数的下降速率应当相对较大。利用拖拉工具从该区域中拖拉一个顶点将对附近顶点的位移起相对有限的影响。仅仅靠近所拖拉顶点的顶点发生位移。顶点的位移将随着与所选顶点之间的距离而快速减少。因此,通过使用选定区域的特征(例如下降速率)来调整几何模型以适应图像数据,本发明的调整方法就需要相对较少的用户交互。
在根据本发明的调整方法的一个实施例中,所述几何模型包括用于手动调整几何模型的工具,并且所述方法还包括用于选择该工具的工具选择步骤。本发明的方法利用诸如参考文献1中的拖拉工具和推动工具之类的工具,来手动调整所述区域以适应图像数据。所述工具可以包含在该区域的特征组内。每一个区域都可以存储其自己的、对所述区域进行调整最为有用的工具。可替换地,该工具可以是综合工具(global tool)。该特征组可以包括用于控制对所述区域进行调整最为有用的工具的工作方式的参数。
在根据本发明的调整方法的另一实施例中,所述几何模型包括用于手动调整几何模型的所述工具的配置,并且所述方法还包括用于选择所述配置的配置选择步骤。所述工具配置可以包含在该区域的特征组内。每一个区域都可以存储其自己的、对所述区域进行调整最有用的工具配置。可选地,所述几何模型可以包括多个配置。拖拉工具配置的实例是用于相对于脉管的脉管调整的工具定向。在相对于脉管中心线的第一预定定向上,该工具充当应用于脉管边界的拖拉工具。在相对于脉管中心线的第二预定定向上,该工具充当应用于脉管中心线的拖拉工具。
在根据本发明的调整方法的另一个实施例中,所述调整方法还包括自动调整步骤,用于自动调整几何模型以适应图像数据,其中,所述手动调整步骤还包括边界条件步骤,用于手动地为所述自动调整步骤设置边界条件。某些用于自动调整的几何模型已经包含了用于调整所述几何模型的一组特征和一个先验信息,其中上述用于自动调整的几何模型例如在J.Weese,V.Pekar,M.Kaus,C.Lorenz,S.Lobregt和R.Truyen所著的、在Proc.IPMI,380-387,Springer Verlag,2001中公开的文献“Shape constrained deformablemodels for 3D medical image segmentation”中所描述的三角形网格,以下将该文献称为参考文献2。在此情况下,用户可以选择几何模型的特定部分,并在手动调整步骤中将这些部分拖拉至所期望的位置。将选定部分以及几何模型中在包围该选定部分的预定形状和大小的体积之外的部分固定到它们的新位置,该新位置定义了边界条件。然后在自动调整步骤中调整几何模型。经调整的几何模型满足强制边界条件。仅调整几何模型中这些未被边界条件所固定的部分。该实施例的优势在于不需要其他的用于手动调整几何模型的区域特征组。在自动调整步骤中使用在这些几何模型中已经可用的区域特征组。可选地,所述几何模型可以包括用于手动设置边界条件的区域特征组。
在根据本发明的调整方法的另一个实施例中,所述调整方法还包括分割步骤,用于分割图像数据。对包含在图像数据中的多个物体使用所述调整方法允许医生描绘所述多个物体。这有助于将图像数据更好的可视化,并使医生能够提取定量的数据,例如包含在图像数据中的物体的几何参数。
本发明的另一个目的是提供开始段落中所述类型的需要相对较少用户交互的调整系统。其通过以下实现的:用于调整几何模型以适应图像数据的适应系统包括:
区域选择单元,用于选择所述几何模型的区域;以及
手动调整单元,用于在所述几何模型内所包含的所述区域的一组特征的基础上,手动调整所述区域以适应图像数据。
本发明的另一个目的是提供开始段落中所述类型的需要相对较少用户交互的图像采集系统。其通过以下实现的:所述图像采集系统包括调整系统,用于调整几何模型以适应图像数据,所述调整系统包括:
区域选择单元,用于选择所述几何模型的区域;以及
手动调整单元,用于在所述几何模型内所包含的所述区域的一组特征的基础上,手动调整所述区域以适应图像数据。
本发明的另一个目的是提供开始段落中所述类型的需要相对较少用户交互的工作站。其通过以下实现的:所述工作站包括调整系统,用于调整几何模型以适应图像数据,所述调整系统包括:
区域选择单元,用于选择所述几何模型的区域;以及
手动调整单元,用于在所述几何模型内所包含的所述区域的一组特征的基础上,手动调整所述区域以适应图像数据。
本发明的另一个目的是提供开始段落中所述类型的需要相对较少用户交互的计算机程序产品。其是通过以下实现的:要由计算机设备加载的计算机程序产品包括指令,用于调整几何模型以适应图像数据,所述计算机设备包括处理单元和存储器,所述计算机程序产品在被加载后提供给所述处理单元执行以下任务的能力:
选择所述几何模型的区域;以及
在所述几何模型内所包含的所述区域的一组特征的基础上,手动调整所述区域以适应图像数据。
基于本说明,本领域技术人员可以实施所述适应系统、图像采集系统、工作站和/或计算机程序产品的修改及其变化,其与所述适应方法的修改及其变化相对应。
本发明的调整方法用于调整几何模型以适应2D、3D和/或4D图像数据。尽管实施例主要描述了调整几何模型以适应3D图像数据的情况,但是在本发明说明书的基础上,将所述调整方法、调整系统、图像采集系统、工作站和/或计算机程序产品延伸到其他维度对本领域技术人员而言是显而易见的。现在可以依惯例通过各种数据采集形式来产生图像数据,例如核磁共振(MRI)、计算机断层造影术(CT)、超声波(US)、正电子发射断层造影术(PET)、单光子发射计算机断层造影术(SPECT)。
将参考附图,相对于以下所述的执行和实施例来阐明根据本发明的适应方法、适应系统、图像采集系统、工作站和计算机程序产品的这些及其它方面,它们由此将变得显而易见,其中:
图1示意性地示出了对股骨的手动调整结果;
图2示出了调整方法的示范性实施例的简化流程图;
图3示意性地示出了对脉管的手动调整结果;
图4示意性地示出了对股骨顶部的手动调整结果;
图5示出了调整系统的实施例的方框图;
图6示意性地示出了图像采集系统的实施例;以及
图7示意性地示出了工作站的实施例。
在全部附图中,相同的参考标记数字用于表明相应的部分。
图1示意性地示出了对股骨100的手动调整效果。将诸如参考文献1中所述的高斯拖拉工具之类的拖拉工具应用于股骨100。在第一个视图101中,将高斯拖拉工具应用于第一个箭头111的起点所指示的区域,并按照第一个箭头111所指示的方向拖拉此区域。由第一个轮廓121来描绘所产生的变形。在第二个视图102中,将高斯拖拉工具应用于第二个箭头112的起点所指示的区域,并按照第二个箭头112所指示的方向拖拉此区域。由第二个轮廓122来描绘所产生的变形。由第一个轮廓121和第二个轮廓122所描绘的变形都是无效的。不应允许这样的变形。不可能使股骨100按照这两个视图所示的方式进行变形。这是因为在第一个视图101和第二个视图102中所示的高斯拖拉工具使用纯粹的几何信息和绝对的本地信息。
在图1中的第三个视图103中,将高斯拖拉工具应用于第三个箭头113的起点所指示的区域,并按照第三个箭头113所指示的方向拉此区域。由第三个轮廓123来描绘所产生的变形。在此,由第三个轮廓123所标记的变形在股骨100的很长范围上扩展。参考文献1的高斯拖拉工具具有一种新特性:其使用了有关股骨100的应用了高斯拖拉工具的区域的先验和更全面的信息。该新特性示出了使用本发明所述的方法超越现有技术的方法的优势。
图2示出了用于调整几何模型以适应图像数据的调整方法200的示范性实施例的简化流程图,调整方法200包括:
初始化步骤210,用于初始化所述几何模型;
区域选择步骤220,用于选择所述几何模型的区域;以及
手动调整步骤235,用于使用所述几何模型内所包含的所述区域的一组特征,手动调整所述区域以适应图像数据。
可选地,调整方法200还包括:
分割步骤205,用于分割图像数据;
自动调整步骤215,用于自动调整所述几何模型以适应图像数据;
工具选择步骤225,用于选择用来手动调整所述区域以适应图像数据的工具;以及
配置选择步骤230,用于设置所述的用来手动调整所述区域以适应图像数据的工具的配置。
进一步参照图2,在开始步骤201之后,方法200前进到GUI步骤202,用于显示图形用户接口,来与用户交流信息。然后方法200前进到分割步骤205。如果在图像数据中没有其他必须对其调整几何模型的对象,则调整方法200就前进到结束步骤299,在此,由调整方法200操纵的处理结束。反之,则选择要被调整来适应图像数据的几何模型。接下来,在初始化步骤210中,初始化选定的几何模型。例如,可以在图像中将几何模型放置在图像数据中的要为其调整几何模型的物体附近。可选的,几何模型还可以进行旋转、平移和/或缩放,以匹配要为其调整几何模型的物体。在初始化步骤210之后,执行自动调整步骤215,用于自动调整初始化后的几何模型以适应图像数据。如果不需要手动调整几何模型,则调整方法200就前进到分割步骤205。如果几何模型中仍然存在需要进行调整以适应图像数据的区域,则在区域选择步骤220中选择该区域。接下来,在工具选择步骤225中,从一组工具中选择一个工具。接下来,在步骤230中,设置所述的与图像数据相关的工具的配置。一旦选择了该区域并且准备了用于手动调整该区域以适应图像数据的工具,则在手动调整步骤235中,手动调整几何模型以适应图像数据。
在手动调整步骤235之后,方法200可以前进到自动调整步骤215,用以在手动调整步骤中定义的新的边界条件的基础上,自动调整几何模型以适应图像数据,方法200前进到区域选择步骤220来选择另一个用于手动调整的区域,方法200前进到工具选择步骤225来选择用于手动调整几何模型以适应图像数据的工具,或方法200前进到配置选择步骤230来选择所述的用于手动调整该几何模型以适应图像数据的工具的配置。如果不需要调整几何模型,方法200就前进到分割步骤205。
在根据本发明的调整方法200的一个实施例中,几何模型基于包括多个顶点的网格。在根据本发明的调整方法200的另一实施例中,网格是多边形网格,例如在参考文献2中所用的三角形网格。多边形网格表示被建模物体的表面。多边形网格相对易于实现。多边形网格的调整很少要求过多的计算时间。可替换地,几何模型可以基于四面体网格。本领域技术人员会理解存在许多可以用本发明的调整方法200来调整的几何模型。在本发明实施例的说明中所用的网格仅是为了说明的目的,并不限制权利要求的范围。
在根据本发明的调整方法200的另一实施例中,调整方法200可选地包括分割步骤205。分割步骤205允许使用调整方法200来调整多个几何模型以适应图像数据集,用以描绘在图像数据内所包含的所关注对象的轮廓。分割步骤205还用作对循环进行控制的控制步骤,该循环包括一个用于调整几何模型以适应图像数据的循环周期。可替换地,调整方法200还可用于调整预定的几何模型以适应图像数据中的对象,以便计算对象的特征。在该情况下,不需要分割步骤205。例如,可以使用调整方法200来计算图像中所包含的心脏对象的体积和其他特征。
在初始化步骤210中,初始化几何模型。例如,可以在图像中将几何模型放置在图像数据中的要为其调整几何模型的物体附近。可选地,几何模型还可以进行旋转、平移和/或缩放,以匹配要为其调整几何模型的物体。初始化步骤可以手动地或自动地进行。
在根据本发明的调整方法200的另一个实施例中,可以在自动调整步骤215中自动调整几何模型以适应图像数据。例如,对于表示几何模型的三角形网格,可以使用参考文献2中所描述的调整方法200。本领域技术人员将了解以下事实:存在许多本领域公知的自动调整方法,并且可以在自动调整步骤215中使用它们的任意组合。
区域选择步骤220是非常重要的步骤,因为选定区域定义了该选定区域的一组特征,其在手动调整步骤中,用来在该手动调整步骤230中将几何模型的所述区域与图像数据对准。可选地,可以使用选定区域来定义用来手动调整几何模型以适应图像数据的一组工具,其中,可以如本发明的另一实施例所述,在工具选择步骤225中从该组工具中选择一种工具。另一种选择是使用在区域选择步骤220中选定的区域和在工具选择步骤中选定的工具,来如本发明的另一实施例所述的,在配置选择步骤230中确定一组工具配置并且确定用于该组工具配置中每个工具配置的所允许参数值。
在根据本发明的调整方法200的另一实施例中,用于在图像数据中选择一个区域的区域选择步骤220被配置为使用GUI工具来旋转和平移所显示的视图(例如投影视图或横截面视图),其包括与根据图像数据渲染的视图交叠的几何模型,以便渲染一个有用的视图来选择要被手动调整以适应图像数据的区域。优选地,为3D图像数据定义旋转和平移。
在根据本发明的调整方法200的另一实施例中,几何模型包括用于手动调整几何模型的工具,并且调整方法200还包括用于选择工具的工具选择步骤225。本发明的调整方法200使用诸如拖拉工具和推动工具之类的工具来手动调整所述区域以适应图像数据。该工具可以包含在该区域的特征组中。每一个区域都可以存储其自己的、对调整所述区域最为有用的工具。可替换地,所述工具可以是综合工具。该特征组可以包括用于控制对所述区域进行调整最为有用的工具的工作方式的参数。可选地,几何模型可以包括一组用于手动调整几何模型的工具。
在基于多边形网格的几何模型的情况下,可以为每一个多边形分配一组可用于手动调整几何模型的区域的工具,所述区域包括该多边形。例如,可以为每一个多边形分配一个拖拉工具和一个推动工具。每一个工具都可以由该多边形的一组特征中所包含的参数来控制。例如,拖拉工具可由用于确定网格的多变形相对于所拖拉的多边形的相对位移的函数来控制,例如高斯函数。可选地,拖拉工具可由多个参数来控制。这对于调整几何模型以适应3D图像数据尤其有用。例如,可以用参考文献1中所述的两个半径来代替高斯拖拉工具的半径。在工具选择步骤225中,当把鼠标指针放在多边形上时,可以由GUI来显示分配给多边形的工具。
在根据本发明的调整方法200的另一实施例中,几何模型包括用于手动调整几何模型的工具配置,并且调整方法200还包括用于选择该配置的配置选择步骤230。工具的配置可以包含在该区域的特征组中。每一个区域可以存储其自己的、对调整所述区域最为有用的工具配置。可选地,几何模型可以包括多个配置。用户可以利用诸如鼠标之类的用户输入设备来操纵工具配置。
预先定义的配置的一个实例是用于与脉管有关的脉管调整的拖拉工具的定向。图3示意性地示出了对表示脉管的几何模型310的手动调整效果。可以将用于调整由等腰三角形330所表示的脉管边界310的工具330定向为平行于或垂直于脉管的中心线320。可以在该工具与脉管中心线320之间的角度的基础上,将工具330的每一个定向分类为平行配置或者垂直配置。在第一个视图301中,工具基本上垂直于中心线。在该配置中,可以将该工具应用于脉管边界310,来使脉管边界310变形。这种调整的结果是如第一个轮廓311所描绘的脉管边界310的变形。在第二个视图302和第三个视图303中,工具330基本上平行于脉管的中心线320。在该配置中,可以将工具330应用于脉管中心线320,来弯曲或伸展脉管边界320。在视图302中用第二个轮廓312描绘了对脉管边界310进行弯曲的结果。在视图303中用第三个轮廓313描绘了对边界310进行伸展的结果。
可替换地,可以在由工具330的应用点所确定的几何模型区域的基础上定义一种类似的效果。如果将工具应用于脉管边界上或其附近,则工具在脉管边界上起作用。如果将工具应用于脉管中心线上或其附近,则工具在脉管中心线上起作用。
可替换地,存在有两种不同的工具,第一种工具用于调整脉管边界,而第二种工具用于调整脉管中心线。可以为几何模型的各个区域分配工具。如果鼠标指针位于脉管边界上或其附近,则第一种工具生效。如果鼠标指针位于脉管中心线上或其附近,则第二种工具生效。
在手动调整步骤235,在几何模型中所包含的区域特征的基础上,手动调整几何模型的所选区域以适应图像数据。所述区域特征包括关于几何模型的该区域的所允许的变形的信息。变形可以取决于许多条件,例如被建模组织的刚度和在该区域中的几何模型表面曲率。因此,用于控制工具的参数取决于该区域。这在图4中进行了图示,图4示意性地示出了对股骨头部400进行手动调整的结果。例如,该工具是高斯拖拉工具。用于控制该工具的参数是对股骨头部的变形进行建模的高斯函数的半径,所述变形被定义为例如高斯函数的方差的平方根。该半径取决于高斯拖拉工具的应用区域。
在图4的第一个视图401中,如箭头411所示地将高斯拖拉工具应用于股骨头部400的半球顶点。用第一个轮廓421来指示所产生的变形。在图4的第二个视图402中,如箭头412所示地将高斯拖拉工具应用于股骨头部400的半球基部。用第二个轮廓422来指示所产生的变形。第一个轮廓421的半径大于第二个轮廓422的半径。导致该差异的原因是以下事实:在股骨头部400的半球顶点的曲率小于在股骨头部400的半球基部的曲率。
本领域技术人员将会理解,可以用多种方式来实现各个部分的位移。例如,在将拖拉工具应用于由网格顶点来表示的几何模型的情况下,特定顶点的位移矢量可以用所拖拉顶点的位移矢量和分配给该特定顶点的权重来定义。可选地,该特定顶点的位移矢量的方向可以不同于所拖拉顶点的位移方向。可替换地,可以使用该特定顶点在以拖拉顶点为原点的本地坐标系统中的本地坐标的矢量值函数,来计算该特定顶点的位移矢量。可替换地,正如本发明另一实施例中所述的,可以优化具有基于所拖拉顶点的位移的边界条件的内部能量表达式,来计算该特定顶点的位移。本领域技术人员将了解此事实:存在其他方法来在由工具引起的网格变形的基础上计算特定顶点的位移。这些方法中的一些可以与拖拉工具相关联;其他的可以与其他工具相关联。
在根据本发明的调整方法200的另一实施例中,手动调整步骤235还包括用于为自动调整步骤215手动设置边界条件的边界条件步骤。在手动调整步骤235中调整几何模型的区域。在手动调整过程中,用户选择几何模型的多个特定部分,例如表示几何模型的网格顶点,并且将这些部分拖拉至所期望的位置。将这些选定部分以及几何模型中在包围这些选定部分的预定形状和大小的体积之外的部分固定到它们的新位置,所述新位置定义了边界条件。然后,在自动调整步骤215中调整几何模型。调整后的几何模型满足强制边界条件。仅调整几何模型中这些未被该边界条件所固定的部分。该实施例的优势在于不需要额外的用于手动调整几何模型的区域特征。基于包含在几何模型中且由自动调整步骤215所使用的特征,来执行整体调整。
在根据本发明的调整方法200的另一实施例中,使用以三角形网格为基础的几何模型和以组合后的内部和外部几何模型势能的最小化为基础的调整方法。在参考文献2中描述了该模型和方法。在该实施例中,用手动选定的点来代替通常用于外部能量的图像驱动边界点。可以通过将非常大的权重分配给外部能量项来实现边界条件,该外部能量项描述了将所调整的部分吸引到其手动选定的位置,所述的非常大的权重比分配给几何模型势能中的其他任何项的权重都大很多。由内部能量参数和边界条件来确定网格的三角形的位移。本领域技术人员将理解,可以通过找到对于力场中运动方程式的固定解,来概括和执行对几何模型势能的最小值的查找。特别是,可以利用模拟来求解运动方程式,所述模拟采用了对能量损耗起作用的项,因此促进了模拟的快速收敛。
可替换地,可以使用以现有技术中已知的另一种成本函数的优化为基础的方法。网格的手动移动的三角形或顶点的位置定义了边界条件,并且所述位置没有被优化。
调整方法200可以有利地用于对自动调整结果进行校正。可替换地,调整方法200可以用作调整几何模型以适应图像数据的可选方法。
在本发明的方法200的所述实施例中的顺序不是强制性的,在不脱离本发明所希望的概念的情况下,本领域技术人员可以改变步骤的顺序或使用线程化模型、多处理器系统或多个进程并行地执行这些步骤。可选地,本发明的方法200的两个步骤可以合并为一个步骤。可选的,可以将本发明的调整方法100的一个步骤分割为多个步骤。
图5示意性地示出了调整系统500的一个实施例,其用于调整几何模型以适应图像数据,调整系统500包括:
初始化单元510,用于初始化几何模型;
区域选择单元520,用于选择几何模型的区域;以及
手动调整单元535,用于在几何模型中所包含的所述区域的区域特征的基础上,手动调整所述区域以适应图像数据。
可选地,调整系统500还包括:
分割单元505,用于分割图像数据;
自动调整单元515,用于自动调整几何模型以适应图像数据;
工具选择单元525,用于选择用来手动调整所述区域以适应图像数据的工具;
配置选择单元530,用于设置工具的配置,例如相对于图像数据来设置;以及
用户接口565,用于与检测系统500进行通信。
在图5所示的调整系统500的实施例中,有三个输入连接器581、582和583,用于输入数据。第一输入连接器581被配置为接收从诸如硬盘、磁带、闪存或光盘之类的数据存储器输入的数据。第二输入连接器582被配置为接收从诸如鼠标或触摸屏之类的用户输入设备输入的数据。第三输入连接器583被配置为接收从诸如键盘之类的用户输入设备输入的数据。输入连接器581、582和583连接到输入控制单元580。
在图5所示的调整系统500的实施例中,有两个输出连接器591和592,用于输出数据。第一输出连接器591被配置为向诸如硬盘、磁带、闪存或光盘之类的数据存储器输出数据。第二输出连接器592被配置为向显示设备输出数据。输出连接器591和592经由输出控制单元590接收各自的数据。
本领域技术人员会理解存在许多方式将输入设备连接到调整系统500的输入连接器581、582和583以及将输出设备连接到调整系统500的输出连接器591和592。这些方式包括但不限于:有线连接与无线连接、诸如局域网(LAN)和广域网(WAN)之类的数字网络、互联网、数字电话网和模拟电话网。
在根据本发明的调整系统500的实施例中,调整系统500包括存储器单元570。存储器单元570被配置为经由任何一个输入连接器581、582和583从外部设备接收输入数据,并在存储器单元570中存储所接收的数据。将数据加载到存储器单元570中允许调整系统500的各个单元对相关数据部分的快速访问。输入数据可以包括但不限于图像数据。可以借助于诸如以下的设备来实现存储器单元570:例如,随机存取存储器(RAM)芯片、只读存储器(ROM)芯片和/或硬盘。优选地,存储器单元570包括RAM,用于存储图像数据集。存储器单元570还被配置为经由存储器总线575,从调整系统500的各个单元接收数据并向其传递数据,调整系统500的各个单元包括:分割单元505、初始化单元510、自动调整单元515、区域选择单元520、工具选择单元525、配置选择单元530、手动调整单元535、用户接口565。存储器单元570还被配置为使得数据经由输出连接器591和592中任意一个对外部设备可用。将来自调整系统500的各个单元的数据存储在存储器单元570中有利地改善了调整系统500的各个单元的性能,以及数据从调整系统500的各个单元到外部设备的传递速率。
可替换地,调整系统500不包括存储器单元570和存储器总线575。调整系统500所使用的输入数据由至少一个外部设备来提供,例如连接到调整系统500的单元的外部存储器或处理器。类似地,由调整系统500产生的输出数据被提供给至少一个外部设备,例如连接到调整系统500的单元的外部存储器或处理器。调整系统500的单元被配置为经由内部连接或经由数据总线在彼此之间接收数据。
在根据本发明的调整系统500的另一实施例中,调整系统500包括用户接口565,用于与调整系统500通信。用户接口565包括:显示单元,用于向用户显示数据;以及选择单元,用于做出选择。将调整系统500与用户接口565相结合允许用户与调整系统500通信。用户接口565还被配置为显示几何模型。用户接口565还被配置为显示用于示出由为适应图像数据而对几何模型进行的调整所造成的几何模型变形的轮廓。用户接口565还被配置为显示用于调整几何模型的工具和工具的配置。用户接口565还被配置为帮助选择工具和配置。可选的,用户接口可以包括调整系统500的多种操作模式,例如手工操作模式和自动操作模式。本领域技术人员会理解,可以在调整系统500的用户接口565中有利地实现更多的功能。
可替换地,调整系统可以使用经由输入连接器582和/或583及输出连接器592连接到调整系统500的外部输入设备和/或外部显示器。本领域技术人员还会理解存在许多用户接口,其可以有利地包括在本发明的调整系统500中。
本发明的调整系统500,例如图5中所示的,可以实现为计算机程序产品,并可以存储在任何适宜的介质上,诸如:例如磁带、磁盘或光盘。该计算机程序可以被加载到包括处理单元和存储器的计算机设备中。计算机程序产品在被加载后,提供给处理单元执行渲染,各种任务的能力。
图6示意性地示出了图像采集系统600的实施例,其使用了本发明的调整系统500,所述图像采集系统600包括:图像采集系统单元610,其经由内部连接与调整系统500相连;输入连接器601;以及输出连接器602。该配置有利地增加了图像采集系统600的性能,为所述图像采集系统600提供了调整系统500的有优势的分割能力。图像采集系统的实例包括但不限于:CT系统、X射线系统、MRI系统、超声波系统、正电子发射断层造影(PET)系统和单光子发射计算机断层造影(SPECT)系统。
图7示意性地示出了工作站700的实施例。所述系统包括系统总线701。处理器710、存储器720、磁盘输入/输出(I/O)适配器730和用户接口(UI)740可操作地连接到系统总线701。磁盘存储设备731可操作地耦合到磁盘I/O适配器730。键盘741、鼠标742和显示器743可操作地耦合到UI740。本发明的调整系统500实现为计算机程序,其存储在盘片存储设备731中。工作站700被配置为将程序和输入数据加载到存储器720中,并在处理器710上执行程序。用户可以用键盘741和/或鼠标742向工作站700输入信息。工作站被配置为向显示设备743和/或向盘片731输出信息。本领域技术人员会理解有许多本领域已知的其它工作站的实施例,并且本实施例起到说明本发明的目的,不能解释为将本发明限于该特定实施例。
应指出,上述实施例是说明性的而不是限制本发明,在不脱离所附权利要求范围的情况下,本领域技术人员将能够设计可替换的实施例。在权利要求中,在括号中的任何参考标记都不应解释为限制权利要求。词语“包括”并不排除未在权利要求中列出的元件或步骤的存在。在元件之前的词语“一”不排除多个此类元件的存在。本发明可以借助于包括几个不同元件的硬件及借助于适当编程的计算机来实现。在系统中,权利要求列举了几个单元,这些单元中的几个可以由同一项硬件或软件来体现。词语第一、第二和第三等的使用并不表示任何排序。这些词语可以被解释为名称。

Claims (9)

1、一种调整方法(200),用于调整几何模型以适应图像数据,所述调整方法包括:
区域选择步骤(230),用于选择所述几何模型的区域;以及
手动调整步骤(235),用于利用所述几何模型中所包含的所述区域的一组特征,来手动调整所述几何模型以适应所述图像数据。
2、如权利要求1所述的调整方法(200),其中,所述几何模型包括用于手动调整所述几何模型的工具,并且所述调整方法(200)还包括工具选择步骤(225),用于选择所述工具。
3、如权利要求1所述的调整方法(200),其中,所述几何模型包括用于手动调整所述几何模型的工具的配置,并且所述调整方法(200)还包括配置选择步骤(230),用于选择所述配置。
4、如权利要求1所述的调整方法(200),还包括自动调整步骤(215),用于自动调整所述几何模型以适应所述图像数据,其中,所述手动调整步骤还包括边界条件步骤,用于为所述自动调整步骤设定边界条件。
5、如权利要求1所述的调整方法(200),还包括分割步骤,用于分割所述图像数据。
6、一种调整系统(200),用于调整几何模型以适应图像数据,所述调整系统(200)包括:
区域选择单元,用于选择所述几何模型的区域;以及
手动调整单元,用于在所述几何模型中所包含的所述区域的区域特征的基础上,手动调整所述区域以适应所述图像数据。
7、一种图像采集系统(600),用于采集图像数据,其包括如权利要求6所述的调整系统(200)。
8、一种工作站(700),其包括如权利要求6所述的调整系统(200)。
9、一种要由计算机设备加载的计算机程序产品,其包括指令,用于调整几何模型以适应图像数据,所述计算机设备包括处理单元和存储器,所述计算机程序产品在被加载后,为所述处理单元提供执行以下任务的能力:
选择所述几何模型的区域;以及
在所述几何模型中所包含的所述区域的区域特征的基础上,手动调整所述区域以适应所述图像数据。
CNA200680034615XA 2005-09-23 2006-09-20 用于手动调整几何模型的先验信息编码 Pending CN101268493A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP05108818.5 2005-09-23
EP05108818 2005-09-23

Publications (1)

Publication Number Publication Date
CN101268493A true CN101268493A (zh) 2008-09-17

Family

ID=37671366

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA200680034615XA Pending CN101268493A (zh) 2005-09-23 2006-09-20 用于手动调整几何模型的先验信息编码

Country Status (4)

Country Link
US (1) US20090115796A1 (zh)
EP (1) EP1929443A1 (zh)
CN (1) CN101268493A (zh)
WO (1) WO2007034424A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8526697B2 (en) * 2008-02-15 2013-09-03 Koninklijke Philips N.V. Apparatus for segmenting an object comprising sub-objects
US20120026168A1 (en) * 2009-04-03 2012-02-02 Koninklijke Philips Electronics N.V. System and method for interactive live-mesh segmentation
US10242430B2 (en) * 2012-03-08 2019-03-26 Brigham Young University Graphical interface for collaborative editing of design space models
US9639992B2 (en) * 2014-10-31 2017-05-02 IntegrityWare, Inc. Methods and systems for multilevel editing of subdivided polygonal data
US10762701B2 (en) 2015-01-28 2020-09-01 Koninklijke Philips N.V. Finite element modeling of anatomical structure

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6804416B1 (en) * 2001-03-16 2004-10-12 Cognex Corporation Method and system for aligning geometric object models with images
AU2003302788A1 (en) * 2002-12-12 2004-06-30 Koninklijke Philips Electronics N.V. Method of segmenting a three-dimensional data set allowing user corrections
WO2005038711A1 (en) * 2003-10-17 2005-04-28 Koninklijke Philips Electronics, N.V. Manual tools for model based image segmentation
US7804997B2 (en) * 2004-06-10 2010-09-28 Technest Holdings, Inc. Method and system for a three dimensional facial recognition system
WO2007034425A2 (en) * 2005-09-23 2007-03-29 Philips Intellectual Property & Standards Gmbh A method of and a system for adapting a geometric model using multiple partial transformations
US20070273688A1 (en) * 2006-05-24 2007-11-29 Yifan Chen System and method of applying geometric surface features to volumetric CAE mesh models

Also Published As

Publication number Publication date
EP1929443A1 (en) 2008-06-11
WO2007034424A1 (en) 2007-03-29
US20090115796A1 (en) 2009-05-07

Similar Documents

Publication Publication Date Title
JP6993334B2 (ja) 自動化された心臓ボリュームセグメンテーション
CN101517526B (zh) 具有3d连通阴影的鼠标指针
Heckel et al. Interactive 3D medical image segmentation with energy-minimizing implicit functions
US8571278B2 (en) System and methods for multi-object multi-surface segmentation
US7773786B2 (en) Method and apparatus for three-dimensional interactive tools for semi-automatic segmentation and editing of image objects
JP4152648B2 (ja) 対象物に含まれる三次元画像をセグメント化する方法
US20060291724A1 (en) Region extraction system, region extraction method and program
US20210132687A1 (en) Device and system for volume visualization and interaction in a virtual reality or augmented reality environment
CN101268492B (zh) 用于使用多个局部变换调整几何模型的方法及系统
US20090202150A1 (en) Variable resolution model based image segmentation
WO2006088429A1 (en) Method and apparatus for editing three-dimensional images
JPWO2012063653A1 (ja) 医用画像表示装置及び医用画像表示方法
CN101203883B (zh) 渐进的基于模型的适配
CN102138161A (zh) 长轴后期增强心脏mri的分割
CN101268493A (zh) 用于手动调整几何模型的先验信息编码
JP2007098028A (ja) モデリング装置、モデリング方法、領域抽出装置、およびプログラム
US20220375099A1 (en) Segmentating a medical image
WO2009101577A2 (en) Interactive selection of a region of interest and segmentation of image data
CN102099833B (zh) 网格冲突避免
US10102347B2 (en) Patient specific anatiomical sketches for medical reports
CN102138159B (zh) 满足预定义空间关系的对象的基于模型的同步分割
JP4736755B2 (ja) モデリング装置、領域抽出装置、モデリング方法及びプログラム
JP4720478B2 (ja) モデリング装置、領域抽出装置およびプログラム
JP2014206931A (ja) 表示処理プログラム、表示処理方法及び表示処理装置
CN102236895A (zh) 用于计算结构的球形度的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20080917