CN103106653A

CN103106653A - 医学图像数据组的配准质量的可视化方法和装置

Info

Publication number: CN103106653A
Application number: CN2012103799784A
Authority: CN
Inventors: S·劳滕施拉格; M·特里尼
Original assignee: Siemens AG
Current assignee: Siemens AG
Priority date: 2011-08-12
Filing date: 2012-08-13
Publication date: 2013-05-15
Also published as: KR20130018168A; DE102011080905B4; US20130038629A1; DE102011080905A1

Abstract

本发明涉及一种医学图像数据组的配准质量的可视化方法(1)，包含下述方法步骤：S1)采集基准图像(10)和目标图像(20)；S2)通过包括非弹性配准方法部分和弹性配准方法部分的配准方法把基准图像(10)与目标图像(20)配准成融合图像(40)；S3)确定变形区，其中，变形区包括从目标图像(20)的图像点(25)到基准图像(10)的图像点(15)的移动矢量(45)，所述移动矢量是由弹性配准方法部分造成；S4)使融合图像(40)的至少一个区域与叠加掩模(50)叠加，变形区的特征参数进入所述叠加掩模；S5)可视化融合图像(40)的至少一个区域和叠加的叠加掩模(50)。此外，本发明还设计一种用于执行这样的方法(1)的相应的装置(100)。

Description

医学图像数据组的配准质量的可视化方法和装置

技术领域

本发明涉及一种医学图像数据组的配准质量的可视化方法。除此之外，本发明还涉及一种相应的用于可视化医学图像数据组的配准质量的装置。

背景技术

在介入放射学中，为了判断一个病例，常常融合，即，叠加地显示不同模态，例如，二维和三维X射线影像(XA)，计算机X射线断层扫描技术(CT)，磁共振断层成像(MR)，正电子发射断层显影(PET)，单光子发射计算机断层扫描技术(SPECT)或超声波扫描术(US)的图像数据组。因此，医生在遇到特殊病例时可以组合不同方法的各自的优势。在理想的情况下为了叠加两个图片只需执行基础的数据图像处理操作，诸如旋转或移动，以正确覆盖图像。在这里指的是两个照片或图像数据组的严格配准或刚性配准。在实践中，将被叠加的影像经常是在不同的时间点或通过不同的照相机位置实现并且因此可以在不同的位置或以不同的形式显示可移动的或柔性的目标，像血管，器官或肌肉。例如，在两个拍摄期间的呼吸伪影或不同的手臂姿势可能是为什么一个血管或一个器官，例如，肝脏，在第一个图像中所占据的位置与在第二个图像中不同的原因。严格配准在这种情况下可能会提供令人不满意的叠加质量。在这些情况下经常使用一个非严格的，柔性的或弹性的配准，其将一个图像中的目标被转换成另一个图像中的相应的目标。在使用一个非严格配准算法的时候可能遇到的问题是，对于配准结果，即，融合图像的观察者来说，很难判定一个图像的哪一部分“弹性”转换程度较高而哪一部分稍微差一些。如果，例如，一个图像中对于诊断特别重要的区域是通过一个非严格配准算法改变，则可能会导致后果深远的误判。

发明内容

因此，本发明的任务是，提供一种医学图像数据组的配准质量的可视化方法。

本发明的基本思想是一种医学图像数据组的配准质量的可视化方法，包含下述方法步骤：

S1)采集基准图像和目标图像；

S2)通过一种包括非弹性配准方法部分和弹性配准方法部分的配准方法把基准图像与目标图像配准成融合图像；

S3)确定变形区，其中，变形区包括从目标图像的图像点到基准图像的图像点的移动矢量，移动矢量是通过弹性配准方法部分造成；

S4)使融合图像的至少一个区域与叠加掩模叠加，变形区的特征参数进入叠加掩模；

S5)可视化融合图像的至少一个区域和叠加的叠加掩模。

在第一方法步骤中，采集基准图像和目标图像。为此可以使用引言中所述的成像法。可以通过相同的或不同的成像法获取这两个图像。这些图像可能是立体图，即，三维图，二维图或三维图像数据组的二维剖面图。图像和图像数据组以下将作为同义词使用。在一个立体图中，图像点被叫做体素，在二维图中，被叫做像素。基准图像和目标图像通常包含检查对象的，例如，人或动物的不同的图像区，但是，交集不是空的。

在第二方法步骤中，基准图像与目标图像通过配准方法配准成一个融合图像。两个图像的配准可以理解为，数字图像处理方法，通过这种方法可以使至少一个相似场景的两个图片尽可能地相互一致。基准图像未发生改变。对于目标图像，确定进行变换，通过变换使目标图像尽可能地与基准图像匹配。就是说，这里涉及的是一个优化问题。图像配准是医学图像处理方面的一个常见任务，对此已经有大量的解决方案。例如，作为配准方法的可使用的优化方法有梯度下降法，下山单纯形法，爬山算法，和模拟退火法。通常，尤其是可以区分为非弹性和弹性配准方法。非弹性配准方法可以理解为一种图像配准方法，在这种图像配准方法中可以使用刚体变换，诸如平移和旋转，仿射变换，诸如缩放和剪切，和射影变换。弹性配准方法指的是一种图像配准方法，在这种图像配准方法可以使用弹性变换，在英语中也叫做“非刚体变换”，像，以插条或多项式为基础的变换。按照本发明的第二个方法步骤使用的是一种包含非弹性配准方法部分和弹性配准方法部分的配准方法。这意味着，如果在配准方法使用了刚体变换，仿射变换，和射影变换，它们可以被归类到非弹性配准方法部分。也可以考虑，通过操作者来执行尤其是移动，旋转和重新调整尺寸。弹性变换可以归类到弹性配准方法部分。替代一种由非弹性配准方法部分和弹性配准方法部分组成的配准方法，也可以先使用一种非弹性配准方法和随后使用一种弹性配准方法。

在第三个方法步骤中确定变形区，其中，变形区包括从目标图像的图像点到基准图像的图像点的移动矢量，而移动矢量是通过弹性配准方法部分造成。尤其是，变形区可以是一个矢量矩阵，所述矩阵重新定位目标图像的每个图像点，使目标图像与基准图像达到最佳重合。

在第四个方法步骤中将融合图像的至少一个区域与叠加掩模叠加，变形区的特征参数进入叠加掩模。因此，可以从变形区中计算出确定叠加掩模及其元素的特征参数，例如，间距尺寸或导数。叠加掩模，在英语中也被叫做Overlay，然后在正确的地点被分配给融合图像的一个区域或整个的融合图像。

在第五个方法步骤中融合图像的至少一个区域和叠加的叠加掩模被可视化。可视化融合图像的至少一个区域可以理解为，尤其是在监视器上显示融合图像的这个区域。叠加掩模及其元素，如上所述，通过变形区的特征参数，例如，间距尺寸被图形地叠加到融合图像的这个区域上。

优选，叠加掩模包含矩阵，矩阵的元素是从变形区的雅可比矩阵中确定。因此，导数或变形区的梯度进入叠加掩模，从而叠加掩模的元素是目标图像与基准图像的匹配强度或变形强度的度量。

在一个优选的结构形式中，叠加掩模含有矩阵，矩阵的元素是从对应矢量的长度中确定，而该对应矢量描述的是目标图像的图像点与基准图像的图像点的对应情况。因此，尤其是来自变形区的移动矢量和尤其是来自变形区的移动矢量的特征参数长度进入叠加掩模。长度或对应矢量值是一种衡量目标图像与基准图像匹配强度或变形强度的度量。

按照本发明的一个优选的结构形式规定，可视化融合图像的至少一个区域和叠加的叠加掩模包含对叠加掩模的元素用颜色进行编码，其中，尤其是可给定的值域可以分配给一个可给定的配色盘，或者其中特别是叠加掩模的元素的值域与预先给定的颜色盘对应。从叠加掩模出发，其元素说明的是，目标图像的图像点通过弹性配准方法部分被移动到那种程度才能与基准图像的相应的图像点重合，而通过颜色编码可以对人类观察者凭直觉就能理解地描述变形程度。优选使用熟悉的颜色模板，诸如透明代表没有变形，绿色代表代表稍微变形，黄色到红色代表变形程度越来越大。变形程度和颜色的对应可以由使用者规定。这时，这种对应可以是绝对地给定的或可给定的配色盘遍及叠加掩模的值域，就是说，例如，最大变形强度的0～10％的变形强度不配给颜色或者是透明的，10％～20％时赋予绿色，等等直到最大变形强度的90％～100％时赋予红色。

在一个优选的结构形式中，在基准图像内至少有一个感兴趣区域可以尤其是通过操作者给定，所述感兴趣区域在基准图像与目标图像配准方法的非弹性配准方法部分之后由配准方法的弹性配准方法部分考虑。这种做法在诊断时的优点是，只有感兴趣区域(英语：Region of Interest＝ROI)需由配准方法的弹性配准方法部分考虑。因此，这个方法步骤描述如下：在基准图像中例如通过操作者标记一个或多个感兴趣区域。当基准图像和目标图像非弹性配准之后进行弹性配准，但是所述弹性配准局限于一个感兴趣区域或多个感兴趣区域。因此，叠加掩模及其与图像融合的可视化只包含感兴趣区域。这种结构形式可以被称作冲压操作，因为，感兴趣区域仿佛是从基准图像中“冲压”出来的并且只有这些感兴趣区域才进行弹性配准。

可以考虑，基准图像和目标图像含有三维图像数据组并且配准图像通过叠加的叠加掩模被可视化成二维剖面图。

本发明的另一个基本思想涉及的是一种用于可视化医学图像数据组的配准质量的装置。这种装置含有接收装置，用于接收被采集的基准图像和被采集的目标图像。此外，该装置还包含至少计算装置，用于借助由非弹性配准方法部分和弹性配准方法部分组成的配准方法使基准图像和目标图像配准，获得一个融合图像，和用于确定变形区，其中，变形区包括通过弹性配准方法部分造成的从目标图像的图像点到基准图像的图像点的移动矢量，和用于使融合图像的至少一个区域与叠加掩模叠加，变形区的特征参数进入叠加掩模。此外，该装置还含有用于可视化融合图像的至少一个区域和叠加的叠加掩模的显示装置。

本发明的一个优选的结构形式规定，装置的计算装置和显示装置的设计功能是执行前述的方法。

使用按照其中一个按照本发明的方法或按照本发明的装置可以实现下述优势：

-提高诊断质量，因为强变形区被特别地标出；

-为医生节省时间，因为“关键的”区域被自动标记，尤其是用颜色醒目地标记出来；

-计算的简单实现和针对病人的剂量应用

-对于医生而言可以进行弹性配准的可能性监控，因为，例如，对于骨头不预计其会变形。

附图说明

下述详细介绍的实施示例是本发明优选的结构形式。

其它优选的结构形式产生于附图及其说明。附图简介：

图1按照本发明方法可视化医学图像数据组的配准质量的流程图；

图2按照现有技术对两个医学图像进行配准的实施示例；

图3对两个医学图像进行配准的实施示例和配准质量的按照本发明的可视化；

图4从目标图像的图像点到基准图像的图像点的移动矢量；

图5可视化叠加掩模的实施示例；

图6用于可视化医学图像数据组的配准质量的装置。

具体实施方式

图1显示的是用于可视化医学图像数据组的配准质量的按照本发明方法1的流程图。方法1包含方法步骤S1～S5，在方法步骤5结束之后以“END”结束。下面是对方法步骤的详细描述：

S1)采集基准图像和目标图像；

S2)通过包括非弹性配准方法部分和弹性配准方法部分的配准方法把基准图像与目标图像配准成一个融合图像；

S5)可视化融合图像的至少一个区域和叠加的叠加掩模。

图2图示的是一个按照现有技术对两个医学图像10和20进行配准的实施示例；基准图像10例如是利用X射线设备通过数字减影血管造影(DSA)提取并且显示部分血管。可以看见下部13，中部11和上部12。目标图像20显示的是同一个血管，但是是在一个稍晚些的时刻拍摄的并且与基准图像10最直观的区别首先在于另一个图像片段和在拍摄时几何投射设置发生变化。血管总体上是按顺时针转动。血管的下部21与基准图像的中部11相对应，上部22已经超出了基准图像的图像片段。可以间接认出的是，血管部分21和22不同，尽管它们在位置上应必须保持相同。原因可能是，比如，病人在两次照相期间轻微地动了一下。使用按照现有技术的配准方法得到的是一个融合图像30。该图中的部分31，32，33与基准图像10中的11，12，13对应。融合图像30多了图像部分34，这个部分源自目标图像20的部分24。此外，所使用的配准方法的非弹性配准方法部分使目标图像20这样转动和在位置上移动，使得无论是基准图像还是目标图像所包含的图像部分尽可能地完全重合。问题是由不同的图像部分12和22的叠加引起的图像部分32。在该示例中，如果两个图像的差异是通过弹性配准方法造成，则使用基准图像的显示，就是说，目标图像的图像点被映射到基准图像的图像点上。另一个方案是，例如，使用两个图像的平均值。最终，融合图像的观察者不能判断，融合图像中哪部分是可靠的，因为它们在这两个图中看起来都一样，而哪些部分可能是伪影。这样的话在诊断时误诊的几率可能会提高。

图3显示对两个医学图像10和20进行配准的实施示例和配准质量的按照本发明可视化。内容与图2中所述的实施示例类似。出发点是，例如为了获取较大的图像片段，基准图像10和目标图像20被相互组合。为此，基准图像10和目标图像20被配准成一个融合图像40。为此使用了一种包含弹性配准方法部分的配准方法，所述方法部分例如转动和移动目标图像并且其结果是融合图像40的图像部分41，43和44。所使用的配准方法的弹性配准方法部分使目标图像的目标变形，使该目标与基准图像的目标尽可能地完全重合。现有技术中公知有许多可使用的配准方法。按照本发明现在确定变形区。变形区含有从目标图像的图像点到基准图像的相应图像点的移动矢量。此外，融合图像被叠加一个叠加掩模。叠加掩模例如是矩阵，矩阵的元素从对应矢量的长度中确定，对应矢量描述的是目标图像的图像点与基准图像的图像点的对应情况。因此，尤其是来自变形区的移动矢量和尤其是来自变形区的移动矢量的特征参数长度进入叠加掩模。长度或对应矢量值是衡量目标图像与基准图像匹配强度或变形强度的度量。最后，融合图像和叠加的叠加掩模被可视化，即，叠加掩模被叠加给融合图像并且这两个都显示。在融合图像40示意性图示结果。通过叠加的叠加掩模可以清楚地识别出图像部分42。融合图像的观察者能够容易地识别出，图像部分42已经变形并且这个区域的配准质量比图像的其他部分要差。

图4示意性显示的是从目标图像的图像点25到基准图像的图像点15的移动矢量45。弹性配准方法把描述目标图像的曲线22的图像点25移动到描述基准图像的曲线12的图像点15。这些移动矢量45进入变形区。

图5显示的是一个可视化叠加掩模50的实施示例。例如，来自变形区的移动矢量长度进入叠加掩模50。就是说，为了与基准图像的相应图像点重合，目标图像的图像点被移动得越远，变形区内相应的矢量就越大并且叠加掩模内相应的值越大。优选，为了可视化，对叠加掩模的元素值颜色编码。这意味着，可给定的值域与可给定的配色盘对应。例如，微小的值，即微小的移动用绿色表示，而大的值用醒目的红色表示。在叠加掩模50内这个内容是通过实线55表示，这条线标记大的移动或大的变形，虚线56代表微小的变形。

最后，图6显示的是一种用于可视化医学图像数据组的配准质量的装置100的实施示例。这种装置100含有接收装置112，例如，电子接口卡，用于接收被采集的基准图像和被采集的目标图像。此外，接收装置112通过电子连接装置106，例如，电气总线，与C型臂X射线仪101连接。C型臂X射线仪101含有固定在支架105上的C型臂102并且在C型臂的两个相对端安装有X射线源103和X射线探测器104。C型臂X射线仪101的设计功能是，通过C型臂102的转动108从不同的拍摄位置拍摄多个X光图像并从中制作三维图像数据组。在患者卧榻107上躺着患者109，例如，人。此外，装置100还包含计算装置110，例如，计算机，用于借助一种由非弹性配准方法部分和弹性配准方法部分组成的配准方法使基准图像和目标图像配准为融合图像，和用于确定变形区，其中，变形区包括通过弹性配准方法部分造成的从目标图像的图像点到基准图像的图像点的移动矢量，和用于使融合图像的至少一个区域与叠加掩模叠加，变形区的特征参数进入叠加掩模。此外，装置100还含有用于可视化融合图像的至少一个区域和叠加的叠加掩模的显示装置111，例如，监视器。

Claims

1.医学图像数据组的配准质量的可视化方法(1)，包含下述方法步骤：

S1)采集基准图像(10)和目标图像(20)；

S2)通过包括非弹性配准方法部分和弹性配准方法部分的配准方法把基准图像(10)与目标图像(20)配准成融合图像(40)；

S3)确定变形区，其中，变形区包括从目标图像(20)的图像点(25)到基准图像(10)的图像点(15)的移动矢量(45)，所述移动矢量是由弹性配准方法部分造成；

S4)使融合图像(40)的至少一个区域与叠加掩模(50)叠加，变形区的特征参数进入所述叠加掩模；

S5)可视化融合图像(40)的至少一个区域和叠加的叠加掩模(50)。

2.根据权利要求1所述的方法(1)，其特征在于，所述叠加掩码(50)包含矩阵，所述矩阵的元素是从变形区的雅可比矩阵中确定的。

3.根据权利要求1所述的方法(1)，其特征在于，所述叠加掩模(50)含有矩阵，所述矩阵的元素是从对应矢量的长度中确定，而所述对应矢量描述的是目标图像(20)的图像点(25)与基准图像(10)的图像点(15)的对应情况。

4.根据前述权利要求中任一项所述的方法(1)，其特征在于，可视化融合图像(40)的至少一个区域和叠加的叠加掩码(50)包含对叠加掩码(50)的元素用颜色进行编码，其中，尤其是给定的值域被分配给给定的配色盘或尤其是叠加掩模(50)的元素的值域被分配给可给定的配色盘。

5.根据前述权利要求中任一项所述的方法(1)，其特征在于，在基准图像(10)内的至少一个感兴趣区域(12)可以尤其是通过操作者给定，所述感兴趣区域在基准图像与目标图像配准方法的非弹性配准方法部分之后由配准方法的弹性配准方法部分考虑。

6.根据前述权利要求中任一项所述的方法(1)，其特征在于，所述基准图像(10)和目标图像(20)含有三维图像数据组并且配准图像(40)通过叠加的叠加掩码(50)被可视化成二维剖面图。

7.一种用于可视化医学图像数据组的配准质量的装置(100)，含有接收装置(112)，用于接收被采集的基准图像(10)和被采集的目标图像(20)，至少一个计算装置(110)，用于借助由非弹性配准方法部分和弹性配准方法部分组成的配准方法使基准图像(10)和目标图像(20)配准为融合图像(40)，和用于确定变形区，其中，变形区包括通过弹性配准方法部分造成的从目标图像(20)的图像点(25)到基准图像(20)的图像点(25)的移动矢量(45)，和用于使融合图像(40)的至少一个区域与叠加掩码(50)叠加，变形区的特征参数进入所述叠加掩码，和用于可视化融合图像(40)的至少一个区域和叠加的叠加掩码(50)的显示装置(111)。

8.根据权利要求7所述的装置(100)，其特征在于，所述计算装置(110)和显示装置(111)的设计功能是执行权利要求2～6中任一项所述的方法。