目标对象的修改过的图像数据的产生的制作方法
【专利摘要】本发明描述了一种用于产生目标对象（O）的修改过的三维图像数据的方法。为此，首先提供目标对象（O）的至少一个通过当前的、规定的三维表面轮廓（S，RS）限界的体积（V）的体素表示并且采集当前的轮廓输入线（IE）。为此然后基于当前的轮廓输入线（IE）与当前的表面轮廓（S，RS）的交点（i0，i1）确定局部的变形区域（D），以用于修改当前的表面轮廓（S，RS）。然后基于所述轮廓输入线（IE）进行所述变形区域（D）中修改过的三维表面轮廓（RS'）的确定。此外还描述了一种可以用来执行该方法的图像处理装置（10）。
【专利说明】目标对象的修改过的图像数据的产生
【技术领域】
[0001]本发明涉及一种基于成像系统的测量数据产生身体、优选器官内部的目标对象的修改过的三维图像数据的方法。此外，本发明涉及一种用于执行该方法的图像处理装置。
【背景技术】
[0002]基于以现代的成像医学技术系统、例如磁共振断层成像装置、计算机断层成像装置、超声波设备等获得的测量数据不仅可以产生穿过身体内部的任意层的二维图像数据，而且也可以产生整个身体的或单独的单个器官的三维图像数据并且将其按照任意视图示出。典型的应用情况是穿过血管、食道、胃或肠的虚拟漫游。此外，目前也可以产生器官、例如心脏的运动的三维图像并且以任意方式将其显示。体积数据的产生通过或者从首先重建后的层图像数据中(例如通过插值等)重建或者直接从测量数据中重建来进行。如果检查的目标对象是特定的结构，例如特定的器官，特定的骨骼等，则将涉及的目标对象分割为体积数据。
[0003]其中可以辅助性地采用目标对象的这些三维图像数据的、目前越来越重要的应用情况是手术规划。这特别符合在切除或部分切除器官的情况。对此的典型例子是肝脏切除术的规划，其首先可以虚拟地借助计算机执行，以便然后能够尽可能精确估计成功概率或风险。在肝脏切除术中在捐献者或患者中保留的、剩余的肝脏体积的大小起主要作用。通常需要保留健康的肝脏体积的至少百分之40至60，从而肝脏可以完整地再生。在肝脏切除术中的主要问题是延伸穿过肝脏组织的不同的血管结构。例如如果必须随之切除在肿瘤附近的血管，则必须同样切除通过该血管供血的肝脏组织。这导致非常复杂的情景。事先虚拟地检验不同的切除情景在此可以有助于确保保留足够的肝脏体积。为此需要可以尽可能精确和容易使用的工具。迄今为止的可用来确定切除体积的方法原则上可以划分为三个不同类别。第一类别中，策略基于利用一种格栅模型修改体积的轮廓。第二策略基于，在穿过器官的不同的层中分别标记保留的体积然后对三维轮廓进行插值。第三种方法对区域基于其通过血管的供血情况而进行分类，以便然后从中确定切除体积。但是所有这些迄今为止的方法要么不是特别精确要么应用起来太复杂。
【发明内容】
[0004]本发明要解决的技术问题是，提供一种用于产生目标对象的三维图像数据的特别简单和可靠的方法以及一种相应的图像处理装置，其中可以非常简单地以高的精度修改目标对象的体积的表面轮廓(和由此可以修改体积)。
[0005]按照本发明的用于产生目标对象的修改过的三维图像数据的方法优选用于产生对于器官切除术的规划图像数据。在此目标对象或器官优选是肝脏，从而以下通常以此为例。但是方法不限于用于规划器官切除术的应用并且特别地不限于用于规划肝脏切除术的应用。方法按照本发明具有以下方法步骤:
[0006]在第一步骤中首先提供目标对象的至少一个通过当前的、规定的三维表面轮廓限界的体积的体素表示。这样的体素表示可以是所谓的体素堆(Voxel-Haufen)。在此对于该目的，优选地二元的体素掩模(Voxel-Maske)就足够，S卩，无需每个体素有特定的强度值，而是仅需要清楚，体素是否是目标对象的当前体积的部分。体积可以是目标对象完整的体积，但也可以是目标对象的重要的、保留的或要切除的部分体积。例如可以是在第一(虚拟或在前面的介入中实际进行的)切除之后保留的部分体积或反过来是在此切除的部分体积。
[0007]体素表示和由此还有所产生的三维图像数据如上所述大多直接或间接基于成像医学技术系统的测量数据，即，图像数据或体素是基于测量数据重建的。体素表示的提供相应地可以例如通过直接从医学技术系统或还从数据库等的其中事先存储了数据的大容量存储器中加载重建的体积数据、剖视图像数据或测量数据来进行。然后可以基于这些数据产生体素表示。替换地，也可以从存储器或从另一个功能单元接收期望的体积的已经产生的体素表示。
[0008]在另一个步骤中然后采集至少一个当前的轮廓输入线。称谓“当前的”在本发明的范围内理解为其总是涉及相应的修改步骤。
[0009]S卩，当从一开始例如在多个步骤中虚拟地进行器官的切除时，一种起始器官掩模包括分割的器官的、在开始时存在的完全的体积，其中该体积具有通过涉及的器官的外部尺寸给出的三维的当前表面轮廓。为此例如可以在初始化步骤中首先提供目标对象的、通过起始表面轮廓限界的起始体积，其在使用医学技术成像系统的测量数据或从中构造的图像数据的条件下基于模型地(即通过基于模型的分割)被建立。
[0010]第一轮廓输入线在该第一步骤中是“当前的”轮廓输入线。如果然后基于该第一修改步骤已经虚拟地切除了体积的部分(即已经基于第一轮廓输入线进行了表面轮廓的第一修改)，则修改过的表面轮廓或修改过的体积在下一个步骤中是当前的表面轮廓或当前的体积并且在该下一个步骤中输入的轮廓输入线于是是新的“当前的”轮廓输入线。
[0011]在输入轮廓输入线之后基于当前轮廓输入线与体积的当前表面轮廓的交点进行局部的变形区域的确定以修改体积的当前表面轮廓。这是与迄今为止的方法的一个主要区另IJ，因为在此当前表面轮廓的和由此当前体积的变形或者说修改仅仅是局部地，即，如果该当前轮廓输入线通过当前体积延伸，则在当前轮廓输入线的特定规定的、基本上局部限界的周围区域中进行变形或修改。
[0012]最后基于轮廓输入线进行变形区域中修改过的三维表面轮廓的确定。修改过的三维表面轮廓然后例如根据在器官切除时的不同观察方式而反映保留的部分或者反映待切除的部分。利用表面轮廓相应地也修改通过当前表面轮廓限界的当前体积。
[0013]修改过的三维图像数据然后可以按照任意方式输出，例如首先在存储器中存储和/或在图像输出装置，诸如屏幕、打印机等上输出。就此而言，该方法也用于控制图像输出装置，这通过将在变形区域中具有修改过的三维表面轮廓的修改过的体积的图像数据转发给用于图像输出装置的相应驱动器并且由此控制该图像输出装置来实现。
[0014]局部限界的变形的特别的优点在于，修改仅具有有限的作用。这特别地意味着，在一定的范围中的表面的修改绝不会“无意地”改变远处部位的表面。特别地在逐步处理切除建议的情况下该特征是重要的，因为否则的话会无意识地改变已经仔细检查过的和良好诊断的区域。
[0015]另一个优点在于，在本发明的范围内对于操作者利用直观方法可以修改存在的(优选简单的、二元的)体素掩模。特别地不需要，操作者在通过体积的多个平行的剖面中输入切割线，这随着体积的几何特征的复杂性增加而几乎不可能是一目了然的。
[0016]用于产生目标对象的这些修改过的三维图像数据的按照本发明的图像处理装置具有以下部件:
[0017]-图像数据接口，用于接收目标对象通过当前的、规定的三维表面轮廓限界的体积，
[0018]-轮廓线采集接口，用于采集至少一个当前轮廓输入线，
[0019]-变形区域确定单元，用于基于当前轮廓输入线与体积的当前表面轮廓的交点确定用于修改当前表面轮廓的局部变形区域，和
[0020]-表面轮廓确定单元，用于基于轮廓输入线确定体积的修改过的三维表面轮廓。
[0021]此外，图像处理装置优选具有用于将图像例如输出到存储器或图像输出装置的图像输出装置，所述图像输出装置可以由图像处理装置控制为输出所产生的三维图像数据。图像输出装置优选是一个或多个显示器、打印机等。
[0022]按照本发明的图像处理装置可以作为独立的装置存在并且例如经由网络等与其中存储了图像数据或体素掩模的合适存储器连接，或者也可以直接与成像系统连接，在成像系统上产生用于图像数据的测量数据。同样地，图像处理装置也可以例如集成在成像系统的图像重建装置或甚至控制装置中。
[0023]按照本发明的图像处理装置特别地也可以按照软件的形式在具有相应的存储可能性的合适的图像处理计算机上实现。这特别地符合变形区域确定单元和表面轮廓确定单元，其例如可以按照软件模块显示实现。同样这些单元也可以构造为硬件部件，例如按照合适构造的ASIC形式。接口同样可以构造为纯硬件、软件或通过软件支持的硬件接口。尽可能按照软件的实现具有优点，即，可以以简单的方式通过软件更新改型迄今为止已经使用的图像处理装置、例如常规的Viewer (观察站)，以便按照本发明的方式工作。本发明由此也包括一种计算机程序产品，其可以直接加载到图像处理计算机的存储器中，具有程序代码段，用于当程序在图像处理计算机上运行时执行按照本发明的所有步骤。
[0024]实施例以及下面的描述包含了本发明的特别有利的改进方案和扩展方案，其中特别地，一类技术方案也可以类似于另一类技术方案的实施例来扩展。同样地实施例的单个技术特征或技术特征组也可以组合为新的实施例。
[0025]当前轮廓输入线的采集原则上可以按照任意方式进行。例如当前轮廓输入线可以自动地通过另外的装置、例如上层布置的用于手术规划的辅助程序等来规定。但是优选地，可以根据操作者命令来采集当前轮廓输入线。为此特别优选地在图形用户界面上显示体素表示的、位于优选可自由选择的输入剖面中的剖视图。输入剖面在此是当前在监视器上显示的表面，在所述表面中然后可以标出新的输入轮廓线。剖视图可以仅仅是感兴趣区域的如下部分剖面，在所述部分剖面中应当改变轮廓或进行通过目标对象的切割，或者是所涉及的输入剖面中体素表示的完整的剖视图。于是借助图形用户界面采集当前的轮廓输入线，方法是，例如随之规定，操作者在剖视图中例如借助图形用户界面的指针沿着哪条线标出线条。相应地，图像处理装置的轮廓线采集接口优选具有图形用户界面。
[0026]然后优选对当前的轮廓输入线进行预检查或似然性检查，例如，线条本身是否相交或其是否主要位于目标对象的体积中。因为当前的轮廓输入线也可以用于校正例如之前进行的虚拟切除切割和例如又添加到体积中，所以在此完全可能的是，当前的轮廓输入线位于目标对象当前保留的体积外部。但是在这种情况下其有意义地应当位于目标对象的初始的总体积中。原则上可以在每个步骤中采集多于一个当前的轮廓输入线，其然后在以下可以相似地先后在修改表面轮廓的情况下被考虑。但是此外为简单起见假定，在本发明不限于此的情况下，每个步骤仅采集一个轮廓输入线。
[0027]然后确定当前轮廓输入线与当前体积的当前表面轮廓的交点。如果存在多于两个这样的交点，则优选将当前轮廓输入线划分为轮廓输入线片段。当前轮廓输入线与体积的当前表面轮廓的交点由此分别形成在随后修改时的待考察的轮廓输入线的端点或在划分轮廓输入线情况下相应的轮廓输入线片段的端点。在该情况中然后优选对于每个轮廓输入线片段根据当前轮廓输入线与体积的当前表面轮廓的交点确定局部的变形区域。如果仅存在两个这样的交点，则清楚的是，仅确定一个围绕相应的轮廓输入线的局部变形区域。
[0028]此外优选总是对于当前轮廓输入线(或在存在多个轮廓输入线片段的情况下分别对于各个轮廓输入线片段)确定当前的变形中央点，其至少大约位于变形区域的中间。特别优选对于变形区域确定在当前轮廓输入线或当前轮廓输入线片段与体积的当前表面轮廓的两个交点之间的中点，作为当前变形中央点。对于每个轮廓输入线，或在轮廓输入线的划分情况下对于每个轮廓输入线片段，相关的变形中央点必定总是位于当前输入剖面中。
[0029]在此可以简单地将沿着在两个属于轮廓输入线或轮廓输入线片段的交点之间的最短距离的中点用作变形中央点。替换地，也可以这样选择变形中央点，使得其例如位于输入剖面中如下路径的中间，该路径沿着体积在两个交点之间的当前表面轮廓。也就是在该情况中当前变形中央点一方面位于输入剖面中，另一方面也位于当前待修改的表面轮廓上。
[0030]为了精确确定变形区域，在特别优选的方法变形中为了确定用于当前轮廓输入线的(或者在相应划分的情况下分别用于当前轮廓输入线片段的)变形区域的边缘而使用搜索算法，优选A*算法。利用这样的合适的搜索算法可以在使用所限定的成本函数的条件下确定路径，该路径优选近似沿着球面从当前轮廓输入线或当前轮廓输入线片段与当前表面轮廓的第一交点延伸至当前轮廓输入线或当前轮廓输入线片段与当前表面轮廓的第二交点以及反向延伸。该近似的球面或最后还有通过边缘限定的圆形区域的中点然后优选是变形中央点。路径在此在当前的、即体积的待修改的表面轮廓上延伸。
[0031]通过这样的工作方式，利用搜索算法找到封闭的边缘轮廓，该边缘轮廓在体积的待修改的表面轮廓上通过两个交点延伸。由此给出明确的变形边界，该变形边界根据轮廓输入线与待修改表面轮廓的交点的距离来将通过当前轮廓输入线引起的修改限制为局部的，特别是还在垂直于当前输入剖面延伸的方向上限制为局部的。
[0032]特别优选地，为了确定在变形区域中的体积的修改过的三维表面轮廓而对对该变形区域内部的表面轮廓的表面点进行插值，并且具体来说基于当前轮廓输入线或当前轮廓输入线片段和优选也基于沿着表面轮廓的变形区域的边缘的表面点来进行插值。在此可以涉及在边缘上、例如在如上所述的路径上的点，但替换地或附加地也可以是与该路径相邻地位于变形区域外部或内部的点。
[0033]在此原则上可以采用具有合适的边界条件的任何插值方法，其中特别地将当前轮廓输入线上或轮廓输入线片段上的点以及提到的沿着边缘的点规定为固定的边界条件或“支点(StUtzstellen)”。这保证在边缘上或边缘外部平滑地过渡到存在的表面中的新插值过的表面。例如可以这样进行插值，使得找到如下表面函数，所述表面函数这样描述修改过的体积的表面轮廓，使得其正好在表面上是零并且例如在体积内部取负值并且在体积外部取正值，或相反。在此已经已知的支点可以立即置零，即，在这些位置处表面函数具有确定的值零。
[0034]特别优选的插值方法基于在期望的支点处定位的径向基函数。利用这些径向基函数可以相对快和简单地实现非常好的插值结果。利用径向基函数插值表面的原理性方法例如在 Greg Turk 等所著的 “Modelling with implicit surfacesthat interpolate，，，ACMTransactions on Graphics (TOG) , 2002,第 855 至 873 页，以及在 Greg Turk 等所著的“Implicit surfaces that interpolate，，，SMI International Conference on ShapeModeling and Applications, 2002,第 62 至 71 页中描述。
[0035]特别优选地，在插值时附加地考虑在事先采集的、较旧轮廓输入线上的位于当前表面轮廓上的点。在此通常涉及如下轮廓输入线，其位于与当前输入剖面成一个角度(例如与之垂直)并且在前面的修改过程中被输入以限定在该另一平面中的切割。为了能够对较旧轮廓输入线进行考虑，优选在每个步骤中存储当前轮廓输入线数据，例如通过对体素体积的、属于这种轮廓输入线的单个体素进行标记。在后面的方法步骤中在输入新的当前轮廓输入线时于是仅还需要对修改过的当前表面轮廓上的确定区域，即，在当前变形区域中的所有体素进行仔细搜索，看其是否标记为(较早的)轮廓输入线并且由此不应当改变或最多在明确的冲突检查之后改变其。
[0036]但是如果轮廓输入线由操作者在不同的输入剖面中输入，则在最不常见的情况下由此要考虑，这些轮廓输入线直接交叉，即，具有共同的交点。在大多数情况下，轮廓输入线或多或少在彼此旁边延伸经过。这意味着，当前轮廓输入线在特定的区域中反对基于较较旧轮廓输入线进行的对表面轮廓的修改。尽管如此，还应当相应考虑较旧轮廓输入线。
[0037]为了解决该冲突，特别优选地确定当前采集的轮廓输入线(或者在轮廓输入线划分的情况下相应地是轮廓输入线片段)与事先采集的输入线的距离。
[0038]完全特别优选地，为此使用在分别相关的输入剖面的交点中两个轮廓输入线(或轮廓输入线片段)之间的欧氏距离。由此可以首先以简单方式采集偏差，即，冲突的强度。
[0039]特别优选地，可以根据在当前轮廓输入线(或轮廓输入线片段)和事先采集的较旧轮廓输入线之间确定的距离而保持不考虑较较旧轮廓输入线的所限定的区段。即，完全地在插值方法中也考虑较较旧轮廓输入线上的点，如上所述例如作为支点来考虑。但不考虑低于与新的轮廓输入线的确定距离、或优选与其中绘入新轮廓输入线的新输入剖面的确定距离的点，其中该距离例如通过在当前轮廓输入线(或当前轮廓输入线片段)和较较旧轮廓输入线之间的距离来规定。
[0040]插值此外优选这样进行，使得在变形区域内部的表面轮廓的表面点和在变形区域外部的表面轮廓的表面点之间的过渡边是平滑的。通过调整梯度，这例如是可能的。为此也可以在实际的变形区域外部的过渡区域中进行稍微的调整。例如可以对于插值方法合适地选择边界条件，其中可以相应地调整梯度，方法是，将在边缘附近、即在边缘外部的体素也如上面描述那样考虑为支点，并且将其值关于描述修改过的表面轮廓的表面函数选择为使得梯度合适。[0041]在找到了变形区域中体积的修改过的三维表面轮廓之后，在考虑该修改过的表面轮廓的条件下最后修改目标对象的体积的体素表示。图像处理装置为此优选具有体积修改单元。
[0042]为了保持计算开销尽可能小，体积的体素表示的修改在此优选可以利用窄带方法沿着所搜索的、所插值的表面轮廓进行。此外为此特别优选地附近使用包括变形区域的范围盒，以便于是还限制计算开销。
【专利附图】
【附图说明】
[0043]以下参考附图借助实施例来详细解释本发明。其中，
[0044]图1是虚拟的肝脏切除术的示意图(左边:肝脏的剖面中的两个规划的切口 ；右边:肝脏的三个合适的立体3D图，从上开始，中间标记出待切除的第一部分，以及第一部分已切除和标记出待切除的第二部分)，
[0045]图2示意性示出了在如图1中那样的肝脏剖面中的所计划的切口，其中第一待切除部分上的切割线改变，
[0046]图3示出了按照本发明的方法的一个实施例的流程图，
[0047]图4示出了在图形用户界面上显示的、带有所绘入的轮廓输入线的肝脏剖视图，
[0048]图5示出了与按照图4的轮廓输入线匹配的变形区域中的表面轮廓的立体三维图，
[0049]图6是用于阐述对于轮廓输入线找出当前表面轮廓上的变形区域的方法的示意图，
[0050]图7是用于阐述找出和解决在当前表面轮廓上新轮廓输入线与的较旧轮廓输入线的冲突的方法的示意图，
[0051]图8立体地三维示出了基于按照图4的轮廓输入线而在按照图5的变形区域中修改的表面轮廓，和
[0052]图9示出了按照本发明的图像处理装置的简化的方框图。
【具体实施方式】
[0053]与大多数常规的方案不同，在按照本发明的方法中，为了对切除术建模利用特别优选是二元的体素掩模进行工作。在这样的二元体素掩模或体素表示中，切除掩模或部分体积总是与切除表面同义。切除表面隐含地通过切除掩模的边界来限定。这特别以图1为例很好地看出。在此，在左图中在纵向通过肝脏延伸的剖面中示意性示出作为目标对象O的肝脏的轮廓。在该剖面中绘入两条轮廓输入线Ie，其将中间保留的体积Vv与肝脏的两个端部待切除的体积Ve分离。在右边与之匹配地三维地示出了该体积，其中上面示出具有肝脏的完整的表面轮廓S的总体积V，其下面示出了具有标记出的、左边待切除部分的体积Ve的总体积V，并且其下面示出了在切除左边部分之后新剩余的并且具有标记出的待切除右边部分的体积Vv。
[0054]在此，在该当前剩余体积Vv的左面示出了，在切除左边待切除的部分之后如何必然在切除面积上产生保留的体积Vv (剩余体积)的改变的、于是“当前的”表面轮廓RS。就此而言，在本发明的方法中可以的是，在整个过程期间将内部的切除状态分别简单地通过掩丰旲来表不，在该掩|旲中将体素标记为切除的或者保留的。这冋时意味着，规划系统能够从开始就进行任意的切除术，但是在规划的切除术中间也可以改变规划并且又改变存在的体积。
[0055]这例如可以根据图2看出。在此应当通过(如在按照图1的左侧上)剖面中新的输入轮廓线Ie现在确定新的表面轮廓RS’，从而将另一个部分体积Va添加到保留的肝脏体积Vv，以便扩大保留的体积Vv。
[0056]因为原则上也可以将每个体积表示划分为体素，所以原则上目前可供使用的所有的切除术规划系统也可以作为起始点使用，以便在规划时借助按照本发明的方法或按照本发明的装置继续工作。此外在体素示图中工作的方法具有优点，即，不再需要将切除表面明确地与器官匹配。此外无需更大的附加开销就可以实现如在图1和2中所示那样限定多重切除。
[0057]以下根据图3详细解释本方法的优选变形。在此明确指出，在该方法的范围内进行的各个方法步骤也可以不取决于具体描述的其它方法步骤，而因此可以有利地在其它方法步骤的范围内使用。
[0058]该方法首先如下开始，S卩，提供待处理的相关体积的(也就是器官的完整的体积或例如应当进行另一切除术的、仍相关的部分体积的)体素表示、优选为二元的体素掩模。在图3中示出的流程规划在此在开始时是肝脏的完整体积。属于肝脏体积的体素利用通常的分割方法而从借助成像系统生成的图像数据中产生。
[0059]该体素表示然后例如可以在图形用户界面上由操作者在任意方向上旋转，以便这样确定特定的输入剖面P”这在步骤I中进行。然后在图形用户界面上显示该输入剖面。这种示图的例子在图4中示出。在此示出正好位于所选输入剖面中的肝脏剖视图。然而在按照图4的该剖视图中已经虚拟切除了肝脏体积的一部分。在下一个步骤的范围中此时应当修改保留的体积的表面轮廓RS并且具体来说正好在相对于所切除的体积的切除表面上进行修改。为此操作者可以例如利用操作终端的鼠标或以类似方式画出新的轮廓输入线IE。这在步骤II中进行。
[0060]在步骤III中然后首先预检查该轮廓输入线Ie的似然性，例如线本身是否相交或者其到底是否涉及肝脏的体积。如果该预检查完成，则在步骤IV中确定轮廓输入线Ie与当前表面轮廓RS、在此也就是与切除表面RS的交点Ici和i10这些交点在图5中再次示出，其中示出对切除面的表面轮廓RS的观察。在此还绘入了初始轮廓线Ik (即当前的输入剖面P1中的表面轮廓RS上的路径)的走向以及新的轮廓输入线Ie的走向。
[0061]如果当前的轮廓输入线Ie与初始的、现在要修改过的表面轮廓多次相交，则存在多于两个交点。于是将轮廓输入线划Ie分为多个轮廓输入线片段。这可以可选地在步骤V中进行。
[0062]然后对于该轮廓输入线Ie的所有各个片段执行接下来的步骤VI至XIII，其中在其它附图中这仅对于单个轮廓输入线Ie示出。但是原则上清楚的是，该方法可以在任意多个轮廓输入线片段或轮廓输入线Ie的情况下工作。为此在步骤VI中首先选择轮廓输入线片段中的一个。如果仅存在一个简单的轮廓输入线Ie，则取消该步骤。以下如果没有另外指出则为简单起见假定，仅单次输入了轮廓输入线IE。但是在步骤VI至XIII的循环内部的方法步骤对于各个输入线片段相应地如对于一个单个的轮廓输入线Ie那样进行。[0063]在步骤VII中然后标记体素掩模的、位于当前待观察的轮廓输入线Ie上的所有体素。在此考虑，每个输入轮廓线或每个输入轮廓线片段或任意其它线在体素表示内部可以简单地示出为排列在彼此上的一定数量的体素并且由此例如也可以按照各个体素的列表的形式存储。同样，在体积内部的体素可以简单地按照合适的方式标记为属于特定的线。
[0064]因为修改在其它方法中无论如何应当这样进行，S卩，修改过的表面轮廓通过选择的输入线Ie延伸，所以首先将在当前轮廓输入线上的体素在存储器中存储为稍后对新表面轮廓进行插值的支点。
[0065]所有支点的集合在以下利用附图标记C表示，因为在此涉及插值的边界条件(约束(Constraint))。
[0066]如图5所示，当前变形区域D在两边都通过交点L、I1限界，在该变形区域中应局部地对于当前轮廓输入线Ie修改当前表面轮廓RS。变形区域D基本上几乎圆形地围绕变形中央点Cs延伸。为了确定变形区域D或变形中央点Cs，存在不同的可能性。
[0067]在图5中示出的实施例中，作为变形中央点Cs，选择在两个交点L、I1之间在当前输入剖面中在起始轮廓线Ik上的中点。从该点出发然后同心地向外在表面上添加越来越多的属于变形区域的体素，直到达到变形区域D的交点h。这通过径向地向外延伸的箭头表不。
[0068]然而在图3中示出的方法流程中使用另外的方法来确定变形中央点Cs以及变形区域D。以下参考附图6解释该方法。该方法也首先基于其中输入了当前轮廓输入线Ie的当前输入剖面P1,以及轮廓输入线Ie与当前表面轮廓RS相交的两个交点L和当前输入剖面P1以及当前表面轮廓RS都在图6中示出。同样在输入剖面P1中示出当前轮廓输入线Ie以及其与当前表面轮廓RS的交点L、h。在此简单地将当前输入剖面P1中在两个交点Ld1之间的几何中点限定为当前变形中央点。这在步骤VIII中(见图3)进行。
[0069]此外在步骤IX中限定两个点(V、QT。该方法在此如下:首先选择在两个交点L、I1之间的输入轮廓线Ie上的中点。这可以通过沿着在交点ItlU1之间的轮廓输入线Ie简单地计数体素来进行。从该中点出发然后分别垂直于当前输入剖面匕中的轮廓输入线Ie地，在距该轮廓输入线Ie间隔很短处、例如每一边上间隔两个或三个体素处确定点C:、。；。这些点仅被需要用来在后面确定，哪个区域位于插值的表面内部并且哪个区域位于插值的表面外部。点(V，C0-同样可以在稍后对修改过的表面轮廓进行插值时用作支点，但是其中同时将与当前轮廓输入线的距离确定为这两个点的函数值。点(V，C0-由此与其距离值一起作为其它支点或约束C而存储。
[0070]在步骤X中然后确定变形区域D，方法是找到封闭的轮廓或通过交点ItlU1延伸的、围绕变形中央点Cs的边缘Mci, M1O该边缘或者该边缘上或在该边缘Md、M1的区域中的这些点(体素)在后面同样用作用于对所修改的表面轮廓进行插值的支点或边界条件C。
[0071]在如结合图6说明的方法中，使用限于表面的、基于图形的A*搜索算法，以确定变形区域D的边缘。在此确定正好两个半轮廓Mtl和M1，其中一个半轮廓Mtl从交点L到交点I1延伸并且另一个半轮廓M1相反从交点I1返回到交点L。在此，在沿着球形轮廓、即沿着示出的圆周的表面上，从相应的起始点1(|或I1开始分别一步一步地或一个体素一个体素地搜索到相应目标点I1或io的最有希望的方向。在仔细搜索期间将表面体素(即，与表面相交的体素)判读为图形的节点并且分别检查限于这些节点或体素的26邻域(B卩，相应待考察的体素的所有直接的邻居体素，包括对角线上的邻居体素)。为了在每个单个搜索步骤中一个体素一个体素地估计最有希望的搜索方向，使用至相应目标点的估计的距离度量fh。
[0072]为了尽可能快速达到封闭的边缘轮廓的期望条件，使用特殊的成本函数
[0073]
【权利要求】
1.一种用于产生目标对象(O)的修改过的三维图像数据的方法，包括以下方法步骤: -提供所述目标对象(O)的通过当前的三维表面轮廓(S，RS)限界的体积(V)的体素表示， -米集当前的轮廓输入线(IE)， -基于所述当前的轮廓输入线(Ie)与所述当前的表面轮廓(S，RS)的交点(i0, I1)确定局部的变形区域(D)，以便修改所述当前的表面轮廓(S，RS),以及 -基于所述轮廓输入线(Ie)在所述变形区域(D)中确定修改过的三维表面轮廓(RS’)。
2.根据权利要求1所述的方法，其特征在于，在图形用户界面(GUI)上示出所述体素表示(V)的位于输入剖面(P1)中的剖视图(SB)，并且借助所述图形用户界面(⑶I)采集所述当前的轮廓输入线(IE)。
3.根据权利要求1或2所述的方法，其特征在于，基于所述当前的轮廓输入线(Ie)与所述当前的表面轮廓(S，RS)的交点将所述当前的轮廓线(Ie)划分为轮廓输入线片段。
4.根据权利要求3所述的方法，其特征在于，相应地根据所述当前的轮廓输入线与所述当前的表面轮廓(S，RS)的交点分别为轮廓输入线片段确定局部的变形区域(D)。
5.根据权利要求1至4中任一项所述的方法，其特征在于，对于变形区域(D)，确定在所述当前的轮廓输入线(Ie)或所述当前的轮廓输入线片段与所述当前的表面轮廓(S，RS)的两个交点(i。，I1)之间的中点，作为当前的变形中央点(Cs)。
6.根据权利要求1至5中任一项所述的方法，其特征在于，为了确定所述变形区域(D)的边缘(Mtl, M1),相应地对于当前的轮廓输入线(Ie)或者当前的轮廓输入线片段借助搜索算法在使用所限定的成本函数的条件下确定路径，该路径优选近似沿着球面从所述当前的轮廓输入线(Ie)或所述当前的轮廓输入线片段与所述当前的表面轮廓(S，RS)的第一交点(io)延伸至所述当前的轮廓输入线(Ie)或所述当前的轮廓输入线片段与所述当前的表面轮廓(S，RS)的第二交点(^。
7.根据权利要求1至6中任一项所述的方法，其特征在于，首先提供所述目标对象(O)的、通过起始表面轮廓(S)限界的起始体积(V)，所述起始体积在使用成像系统(I)的测量数据(MBD)的条件下基于模型地建立。
8.根据权利要求1至7中任一项所述的方法，其特征在于，为了确定在所述变形区域(D)中的修改过的三维表面轮廓(RS’)，基于 -所述当前的轮廓输入线(Ie),以及 -沿着所述表面轮廓(S，RS)的变形区域(D)的边缘(Mtl, M1)的表面点， 对所述变形区域(D)中的表面轮廓(RS’ )的表面点进行插值。
9.根据权利要求8所述的方法，其特征在于，在插值中考虑其它的、事先采集的轮廓输入线(R1)上的点。
10.根据权利要求9所述的方法，其特征在于，确定当前采集的轮廓输入线(Ie)或轮廓输入线片段与事先采集的轮廓输入线(R1)的距离(d)。
11.根据权利要求10所述的方法，其特征在于，根据在所述当前的轮廓输入线(Ie)或轮廓输入线片段与所述事先采集的轮廓输入线(R1)之间所确定的距离(d)，保持不考虑所述事先采集的轮廓线(R1)的所限定的区段。
12.根据权利要求8至11中任一项所述的方法，其特征在于，将插值进行为使得在所述变形区域(D)内部的表面轮廓(RS’ )的表面点和在所述变形区域(D)外部的表面轮廓(S，RS)的表面点之间的过渡边是平滑的。
13.根据权利要求1至12中任一项所述的方法，其特征在于，在考虑所述变形区域(D)中修改过的三维表面轮廓(RS’ )的条件下修改所述目标对象的体积(V)的体素表示。
14.一种用于产生目标对象(O)的修改过的三维图像数据的图像处理装置(10)，具有以下部件: -图像数据接口(11)，用于接收所述目标对象(O)的通过当前的、规定的三维表面轮廓(S，RS)限界的体积(V)， -轮廓线采集接口( 13)，用于采 集当前的轮廓输入线(Ie), -变形区域确定单元(14)，用于基于所述当前的轮廓输入线(Ie)与所述当前的表面轮廓(S，RS)的交点(U I1)确定局部的变形区域(D)，以便修改所述当前的表面轮廓(S，RS),和 -表面轮廓确定单元(15)，用于基于所述轮廓输入线(Ie)确定所述变形区域(D)中修改过的三维表面轮廓(RS’)。
15.一种可以直接加载到图像处理计算机(10)的存储器中的计算机程序产品，所述计算机程序产品具有程序代码段，用于当程序在所述图像处理计算机(10)上运行时执行按照权利要求1至13中任一项所述方法的所有步骤。
【文档编号】G06T1/00GK103578076SQ201310323763
【公开日】2014年2月12日 申请日期:2013年7月30日 优先权日:2012年7月31日
【发明者】J.克雷特什莫, G.索扎, C.蒂特詹 申请人:西门子公司