专利名称：基于点的自适应弹性图像配准的制作方法
技术领域：
本发明涉及数字成像领域。特别是，本发明涉及一种将第一图像与第 二图像配准的方法，且涉及一种用于将第一图像与第二图像配准的图像处 理设备和软件程序。
背景技术：
图像配准的目的，例如，在医学成像应用中，是补偿图像间的差异， 例如，由于患者移动、不同的扫描器形式、解剖学中的变化等引起的差异。 诸如刚体变换或仿射变换那样的全局配准方法通常无法应付局部差异。一 种已知的相关解决方案称为弹性配准。对医学图像进行鲁棒性的弹性配准 是一个困难的问题，这是目前加强研究的主题。
基于点的弹性配准包括以下步骤定义与第一图像相关的一组控制点；
然后在这些控制点上执行第一图像的弹性变形，以便实现第一图像与第二 图像的最佳空间对应，其中，用相似性度量来量化该对准情况。在参数化 几何变换的情况下，通过计算最佳参数设置来达到该最佳对准，该最佳参 数设置对于弹性配准而言通常意味着控制点的最佳数量和位置以及在这些 控制点上的位移参数(定义了第一图像的弹性变形程度)。
最广泛使用的弹性图像配准的变形类型是B样条，其是在控制点的规 则网格上定义的。 一般来说，当需要对第一图像进行高度弹性变形时，就 需要定义高密度的控制点。在控制点的规则网格的情况下，就需要对于整 个第一图像提供高密度，即使是仅对于第一图像中的较小面积才需要这类 高弹性变形。至少需要确定关于每一个控制点的位移参数，以致于在此情 况下将需要优化大量的参数，这就需要较长的计算时间。
可以利用以控制点的不规则网格为基础的变换来克服上述缺陷。固定 数量的控制点在第一图像上的位置都认为是自由参数(要被优化)，在优化处 理过程中，这些参数以及控制点位移参数都能够改变。这允许根据需要移动控制点，并能够对于第一图像中需要高度弹性变形的区域提供高密度的 控制点而在其他图像区域中，控制点密度可以更小。例如，国际专利申请
号WO2005/057495公开了一种弹性变形方法，其中，在几个控制点上对第 一图像施加力场，并自动查找在被施加了力的控制点的最佳位置，以便最 小化第一个与第二图像之间的差异。
然而，在图像配准处理开始时就固定了控制点的数量，并在整个处理 过程中始终保持该固定数量。由于在配准处理之前无法获知控制点的最佳 数量和初始相对位置，所以将需要比必要数量的控制点更多数量的控制点 来获得可接受的图像配准结果，这进而又意味着执行优化处理所需的计算 量和时间也不必要地较高。
本发明的目的是提供一种图像配准方法，其中，可以优化控制点的数 量，以便最小化执行图像配准所需的计算量和时间，同时不会降低配准结 果的质量。本发明的目的还在于提供一种相应的图像处理设备和软件程序。
发明内容
根据本发明，提供了一种用于将第一图像与第二图像进行配准的方法，
所述方法包括
在所述第一图像内放置至少一个控制点(S6)，为确定关于所述至少一个 所述控制点的、确定表示所述第一图像的弹性变形的、定义了位置和位移
参数的第一参数设置(S7)，其中，所述弹性变形提高了变形后的所述第一图 像与所述第二图像之间的相似性，以便使所述第一图像弹性变形(S8)，从而
提高所述第一图像与所述第二图像之间的相似性，并且重复以下步骤 在所述第一图像内放置至少一个附加控制点(S6);
基于所述第一参数设置，确定关于为所述至少一个附加控制点确定表 示所述第一图像的进一步弹性变形的、的、定义了位置和位移参数的第二
参数设置(S7)，直至满足了预定准则，其中，所述进一步弹性变形以便使所 述第一图像弹性变形(S8)，从而进一步提高了在进一步弹性变形后的所述第 一图像与所述第二图像之间的所述相似性；直至满足了预定准则为止。
因而，本发明的目的是通过以下实现的从一个或多个控制点(优选为
单个控制点)开始，并在每一次弹性变形操作之后迭代地添加新的控制点，直至满足预定准则。照这样，不必提前指定控制点的数量，而且控制点的 数量可以自动适应于变形区域的复杂性。在优选的实施例中，在每一次弹 性变形操作之后迭代地添加新的控制点，直至第一图像与第二图像之间的 相似性至少达到预定级别。
优选地，在每次添加一个或更多附加控制点时，确定关于第一图像内
所有控制点的最佳参数设置。因而， 一般来说，优化一组N个控制点，而 且以所得到的配置充当关于接下来对一组N+M个控制点进行优化的开始 点，其中N和M是整数。在本发明的一个示范性实施例中，M=l，而且在 下一次优化操作之前，优选地在第一图像内随机地放置单个附加控制点。 因此，在一个优选的示范性实施例中，从仅一个(N4)随机放置的控制点的 初始配置开始，逐个添加控制点，直至第一图像与第二图像之间的相似性 再无(显著的)提高为止。
有利地，对每一个控制点的参数设置进行优化，以便优化相似性度量(例 如，其可以是第一图像与第二图像之间的平方差，但是也可以使用许多其 他类型的相似性度量，包括互信息或互相关，而且本发明并非必须受限于 这一点)。在优选实施例中，在每一次弹性变形操作之后获得一个相似性度 量，并且可以确定对第一图像与第二图像之间的相似性进行了提高的量(即， 由最后一次迭代引起的相似性度量的提高)，并且可以将该提高的量和预定 准则进行比较，其中，只有在没有满足所述预定的准则的情况下才添加一 个或多个附加控制点。
此外，根据本发明，提供了一种用于执行第一图像与第二图像的配准 的图像处理设备，所述设备包括用于存储所述第二图像的存储器、用于接 收关于所述第一图像的图像数据的装置以及如下配置的处理装置
一开始，在所述第一图像内放置至少一个控制点，并确定关于所述至 少一个控制点的、定义了位置和位移参数的第一参数设置，以便使所述第 一图像弹性变形，从而提高所述第一图像与所述第二图像之间的相似性； 然后重复以下步骤-
在所述第一图像内放置至少一个附加控制点，确定关于所述至少一个 附加的控制点的、定义了位置和位移参数的第二参数设置，以便使所述第 一图像弹性变形，从而进一步提高所述第一图像与所述第二图像之间的所述相似性；直至满足预定准则为止。
还根据本发明所述提供了一种用于将第一图像与第二图像进行配准的 软件程序，其中，所述软件程序使处理器能够执行-
一开始，在所述第一图像内放置至少一个控制点，并确定关于所述至 少一个控制点的、定义了位置和位移参数的第一参数设置，以便使所述第 一图像弹性变形，从而提高所述第一图像与所述第二图像之间的相似性； 然后重复以下步骤
在所述第一图像内放置至少一个附加控制点，确定关于所述至少一个 附加的控制点的、定义了位置和位移参数的第二参数设置，以便使所述第 一图像弹性变形，从而进一步提高所述第一图像与所述第二图像之间的所 述相似性；直至满足预定准则为止。
本发明的这些及其他特征将参照此处所描述的实施例而显而易见且得 以阐明。
现在将仅借助于实例并参照附图来描述本发明的实施例，其中
图1表示根据本发明的一个示范性实施例的图像处理设备的示意图，
该图像处理设备适合于执行根据本发明的一个示范性实施例的方法；以及 图2表示根据本发明的方法的一个示范性实施例的简化流程图。
具体实施例方式
图1描述了根据本发明的图像处理设备的一个示范性实施例，该图像 处理设备用于执行根据本发明的方法的一个示范性实施例。图1所示的图 像处理设备包括连接到存储器2的中央处理单元(CPU)或图像处理器1，存 储器2用于至少存储第一和第二图像、控制点的参数设置、以及第一和第 二相似性度量。图像处理器1可以连接到多个输入/输出网络或诊断设备， 例如MR设备、CT设备或超声扫描仪。图像处理器1还连接到显示设备4(例 如，计算机监视器)，以显示在图像处理器1内计算或调整的信息或图像。 操作者可以通过键盘5和/或其他在图1中未示的输入/输出设备来与图像处 理器1交互。虽然在以下是参照医学应用来描述该方法的，但是应当注意的是，本 发明可以应用于所有要求配准的多维数据集或图像。例如，本发明可以应 用于产品的质量检验，在该情况中将实际产品的图像与参考产品的图像进 行比较。另外，该方法可以应用于材料测试，例如，用于监测感兴趣的对 象在特定时间段内的变化。
图2表示根据本发明的、用于将第一图像与第二图像进行配准的方法 的一个示范性实施例的流程图。在步骤S1中开始之后，在步骤S2中，在 第一图像区域内随机地放置单个控制点。接下来，在步骤S3中计算关于该 单个控制点的最佳参数设置，这种参数设置至少包括该控制点在第一图 像区域内的最佳位置，以及用于定义将应用于由此定位的该控制点上的弹 性变形的程度的位移参数。由此优化这些参数设置，以利用单个控制点来 实现第一和第二图像的最佳对准。 一旦在步骤S4中在单个控制点上将所需 的弹性变形应用于第一图像，就在步骤S5中计算相似性度量，该相似性度 量表示第一个与第二图像之间的对准程度，这是利用单个控制点来实现的。 一种适当的相似性度量是第一与第二图像之间的平方差，而且本发明的该 示范性实施例的方法旨在优化相似性度量，以实现两个图像之间的最佳对 准，同时最小化执行图像配准所需的计算量和时间。
接下来，在步骤S6中，在第一图像区域内随机地放置一个附加的控制 点，并在步骤S7中计算关于第一图像区域内所有控制点的最佳参数设置， 以实现第一和第二图像的最佳对准。
在步骤S8中， 一旦以恰当定位的控制点实现了第一图像的弹性变形， 就在步骤S9中计算新的相似性度量。在步骤S10中，根据某种预定的停止 准则，来比较该新的相似性度量和先前计算的相似性度量(例如，用一个阈 值来比较该差异)。如果在步骤Sll中不满足该预定的停止准则(例如，当前 相似性度量与先前相似性度量之间的差异至少等于该阈值，该阈值指示已 经对第一与第二图像之间的相似性进行了至少某个预定量的提高)，则方法 返回步骤S6，在该步骤S6中，在此添加另一个控制点并重复以上处理。一 旦满足了该停止准则(例如，当前相似性度量与先前相似性度量之间的差异 低于上述阈值)，则在步骤S12中该方法结束并完成图像配准处理。
一般来说，对两个图像I,、 12的配准包括建立t变换，以使得t(W与l2之间的差异按照一个相似性度量^^是最小的。在图像配准中，通常形成优
化问题，以使得表示理想变换tc的参数矢量c将将会使目标函数f(c)=COrr (tc(iai2)最大化。因此，根据本发明的示范性实施例所述，可以将该优化问 题公式化成对于每一次迭代査找第一图像内的给定控制点集合的最佳位置 以及它们的最佳位移参数。对于本领域技术人员将会显而易见的是，可以 使用许多不同类型的变换，而且可以在以下中找到实例例如在1999年的 Journal of Computer Assisted Tomography 23(5)第800-805页上由D. Rueckert 等人所著的"Comparison and evaluation of rigid, affine and non-rigid registration of breast MR images", 以及待干U出的Physics in Medicine and Biology 2005中由V. Pekar禾口 K. Rohr所著的"An adaptive irregular grid approach for 3-D deformable image registration"。
可以利用标准的数字优选技术来解决公式化后的优化问题，例如，如 1965年的Computer Journal,(7)第308-313页上由J. A. Nelder和R. Mead所 著的"A simplex method for function minimization"所述的下山单纯形 (downhill simplex)方法。
因而，从第一图像区域内随机放置的单个控制点开始，利用局部收敛 优化策略来査找关于该控制点集合的最佳配置，其中所有控制点(包括在以 前步骤中优化的控制点)的位置和位移参数都认为是自由参数。在最初的几 次迭代中，由于将要优化较少数量的参数，所以可以很快地执行关于仅一 个或少量控制点的优化步骤。和采用基于固定数量控制点的局部优化策略 的图像配准的现有技术方法相比，该提议的方法利用数量少的多的控制点 得到了可比拟的甚至更好的结果。因此，该提议的方法可以显著地提高图 像配准处理的速度，而且满足应用指定的质量要求。这在严格要求时间的 应用中是最为重要的，所述严格要求时间的应用例如内-外科配准，在该情 况中必须仅在所关注的应用指定的区域上(例如，临床焦点)实现最佳配准精 确度。此外，控制点数量上的迭代式增加增强了配准算法的鲁棒性。对于 仅能利用大量控制点来实现的、要求高精确度的应用而言，可以按照恰当 的方式来定义所述终止准则。
应当注意的是，本发明可以应用于CT图像、核磁共振图像(MRI)、正 电子发射断层造影图像(PET)、单光子发射计算机断层造影图像(SPECT)、或超声(us)形式。而且，可以使用其他数据集。
应当注意的是，上述实施例是说明性的而非限制本发明，本领域技术 人员能够在不脱离附带的权利要求所定义的本发明范围的情况下设计许多 备选的实施例。在权利要求中，任何置于括号内的参考标记都不应构成对 该权利要求的限制。"包括"和"包含"等词语并不排除除了在任何一个权 利要求或整个说明书中所列的那些元件或步骤之外还存在其他元件或步 骤。对元件的单数称谓并不排除对这类元件的复数称谓，反之亦然。可以 借助于包含了几个不同元件的硬件来执行本发明，并且还可以借助于适当 程序化后的计算机来执行本发明。在罗列了若干个装置的设备权利要求中， 可以将这些若干个装置合并成一个和相同的硬件项。在相互不同的从属权 利要求中记载特定手段的这一事实并不表示不能有利地使用这些手段的组合。
权利要求
1、一种用于将第一图像与第二图像进行配准的方法，所述方法包括在所述第一图像内放置至少一个控制点(S6)，为所述控制点确定表示所述第一图像的弹性变形的、定义了位置和位移参数的第一参数设置(S7)，其中，所述弹性变形提高了变形后的所述第一图像与所述第二图像之间的相似性，并且重复以下步骤在所述第一图像内放置至少一个附加控制点(S6)；基于所述第一参数设置，为所述附加控制点确定表示所述第一图像的进一步弹性变形的、定义了位置和位移参数的第二参数设置(S7)，其中，所述进一步弹性变形进一步提高了在进一步弹性变形后的所述第一图像与所述第二图像之间的所述相似性；直至满足了预定准则为止。
2、 根据权利要求1所述的方法，其中，所述预定准则包括所述第一图像与第二图像之间的相似性至少达到了预定级别。
3、 根据权利要求1所述的方法，其中，在所述方法开始时(S1)，在所 述第一图像内随机地放置单个控制点(S2)，并且确定关于所述单个控制点的 参数设置(S3)。
4、 根据权利要求1所述的方法，其中，在每次添加一个或多个附加控 制点时，确定关于所述第一图像内所有控制点的最佳参数设置(S7)。
5、 根据权利要求1所述的方法，其中，向所述第一图像逐个添加控制 点(S6)，并执行关于各个控制点集合的弹性变形(S8)，直至满足了所述预定 准则为止。
6、 根据权利要求1所述的方法，其中，优化每一个控制点的参数设置， 以便优化相似性度量。
7、 根据权利要求1所述的方法，其中，在每个弹性变形操作之后获得 相似性度量(S9)，并确定对所述第一图像与所述第二图像之间的相似性进行了提高的量(SIO)，并且将所述量与停止准则进行比较(Sll)，其中，仅在没有满足所述停止准则的情况下，添加一个或多个附加控制点。
8、 根据权利要求7所述的方法，其中，将所述量与阈值进行比较，并 且其中，只有在所述量超过了所述阈值的情况下，添加一个或多个附加控 制点。
9、 一种用于执行第一图像与第二图像的配准的图像处理设备，所述设 备包括用于存储所述第二图像的存储器(2)、用于接收关于所述第一图像的 图像数据的装置、以及如下配置的处理装置(l):一开始，在所述第一图像内放置至少一个控制点，并且为所述控制点 确定表示所述第一图像的弹性变形的、定义了位置和位移参数的第一参数 设置，其中，所述弹性变形提高了变形后的所述第一图像与所述第二图像 之间的相似性；然后步骤为在所述第一图像内放置至少一个附加控制点；以及基于所述第一参数设置，为所述附加控制点确定表示所述第一图像的 进一步弹性变形的、定义了位置和位移参数的第二参数设置，其中，所述 进一步弹性变形进一步提高了在进一步弹性变形后的所述第一图像与所述 第二图像之间的所述相似性；直至满足了预定准则为止。
10、 一种用于将第一图像与第二图像进行配准的软件程序，其中，所 述软件程序使处理器(l)能够执行一开始，在所述第一图像内放置至少一个控制点，并且为所述控制点 确定表示所述第一图像的弹性变形的、定义了位置和位移参数的第一参数 设置，其中，所述弹性变形提高了变形后的所述第一图像与所述第二图像 之间的相似性；然后步骤为-在所述第一图像内放置至少一个附加控制点；以及基于所述第一参数设置，为所述附加控制点确定表示所述第一图像的 进一步弹性变形的、定义了位置和位移参数的第二参数设置，其中，所述 进一步弹性变形进一步提高了在进一步弹性变形后的所述第一图像与所述 第二图像之间的所述相似性；直至满足了预定准则为止。
全文摘要
一种基于点的弹性配准方法，用于将第一图像与第二图像进行配准。一开始，在源图像区域内随机地放置单个控制点(S2)，并确定关于该单个控制点的最佳参数设置(S3)，用以执行关于该第一图像的弹性变形(S4)，以便优化相似性度量。然后逐个添加附加控制点(S6)，并且每次都对于新的控制点集合重复该弹性变形处理(S8)，直至满足预定停止准则，例如，所得到的该相似性度量的提高不再超过某个预定阈值。因而，提供了一种高速、高质量的配准方法，而不必在一开始就指定控制点的数量。
文档编号G06T3/00GK101305395SQ200680041753
公开日2008年11月12日 申请日期2006年11月9日 优先权日2005年11月10日
发明者A·弗伦茨, H·维施曼, I-C·卡尔森, S·雷尼施 申请人:皇家飞利浦电子股份有限公司