专利名称：用于分析两张成像图像的方法以及医学成像系统的制作方法
技术领域：
本发明涉及一种用于分析一个器官的两张相对应的、在不同时刻拍摄的成 像图像的方法，以及一种医学成像系统。
背景技术：
为了澄清很多医学问题需要进行纵向成像的检查。这意味着用相同的成像 方法在不同时刻多次对待检查器官、器官系统或身体部位进行检查。通过比较 该器官的成像图像，例如可以确定特定疾病的发展情况。通常只能通过随着时 间出现在某个器官上的变化来识别或至少考虑到特定疾病，尤其是退化的疾病。
这种疾病的例子是阿兹海默症。用常规的检查方法很难诊断出这种疾病， 尤其是当这种疾病处于早期时。这种疾病属于神经退化性疾病，公知在阿兹海 默症发病期间由于神经细胞的死亡而在特定的脑部区域中显示出萎缩，此外在 顶叶和颈叶的脑灰质(也称为脑皮层)也是如此。
虽然这种变化的出现是一种公知事实，但成像方法、尤其是MRT(磁共振 断层造影)只能非常有限地用于诊断阿兹海默症，即使采用体积测定方法来测定 特定的大脑皮层区域。这部分是由于待测量和检测的萎缩的程度在数量上位于 不同个体的特定脑部区域大小的浮动范围之下。
克服该问题的 一种方法是在不同时刻检查患者大脑，并通过比较成像图像 确定萎缩的发展。但是该方法的缺点是待检测的差异一部分特别小，从而比较 该成像图像的用户很容易就忽略了该差异。
此外，如果像在神经学中那样通常借助MRT方法来获得成像图像，则该 成像图像具有这样的特点，即在不同时刻完成的成像图像可能具有不同的扭曲， 尤其是失真。这是由于为了成像而以公知方式采用不同的、最精确地相互协调 的磁场的MRT方法对干扰影响反应灵敏。因此一般不能对这些成像图像始终 保证刚好相同的拍摄条件，从而不同的成像图像始终具有相同的并因此是相类 似的几何失真。仅当借助在拍摄前完成的模型测量检测并平衡这些几何失真之 后，才能更好地平衡这些几何失真。但这由于高昂的费用以及由此带来的实践 中的成本而没有采用。
因此用户通常无法由不同成像图像中的很小的、但是与诊断有关的差异推 断出是几何失真，还是器官中的解剖事件的实际变化。
由于这一事实，MRT检查目前还不属于普遍承认的和稳定的阿兹海默症诊
断方法，而是用于诊断其它疾病。
目前存在各种克服上述问题的方案。在文献Freeborough PA, FoxNC， "The boundary shift integral: an accurate and robust measure of cerebral volume changes from registered repeat MRI", IEEE Trans.Med.Imaging 1997;16:623-629中，公开了
一种用于在重复的、三维MRT拍摄中测量体积变化的方法。其中将不同时刻 完成的MRT照片刚性地相互配准，并将它们的强度值之差积分。利用这种方 法虽然可以检测整个大脑的体积差异，但是该方法没有考虑不同的几何失真。 另外，体积变化不能与特定的大脑区域对应，从而该方法通常在查找或进行特 定的诊断时没有向用户提供足够的帮助。
虽然存在利用该方法只检查特定的脑部区域以获得形态变化的至少粗略 定位，但是还是没有考虑几何失真。
即使该问题是以阿兹海默症以及借助MRT进行诊断为例来描述的，在其 它医学课题中也存在类似的问题。其例子是肿瘤疾病或骨质疏松的过程。同样 在MRT检查也会出现上述几何失真的问题，但是其它检查模态如计算机断层 造影也可能存在类似的问题，例如在有错的校准中。
发明内容
因此本发明要解决的技术问题在于提供一种方法，其在分析一个器官的两 张在不同时刻拍摄的成像图像时以有效的方式向用户提供支持。此外本发明还 要解决的技术问题是提供一种医学成像系统，该系统在分析一个器官的两张在 不同时刻拍摄的成像图像时以有效的方式向用户提供支持。
本发明用于分析至少一个器官系统的两张相对应的、在不同时刻拍摄的成 像图像以确定医学疾病图像中的病理变化的方法，其中所述至少一个器官系统 具有第一区域和在成像图像中与该第一区域不同的第二区域，其中第二区域比 第一区域更强烈地被该医学疾病图像中的病理变化覆盖，该方法包括以下步骤
-确定针对第一区域的配准，从而通过该配准将第一成像图像的第 一区域
与第二成像图像的第一区域相互对准，
-将针对第一区域的配准扩大到扩展的配准，从而通过该扩展的配准一起 采集第二区域，
-借助所述扩展的配准变换两张成像图像之一 ，
-显示经过变换的成像图像和另一张成像图像，和/或通过比较经过变换的 成像图像的第二区域与另 一 张成像图像的第二区域来确定两张成像图像在第二 区域上的差异。
在此，本发明的方法利用了病理变化覆盖第二区域强于第一区域的事实。 通过对两张成像图像中被病理变化较弱涉及的第 一 区域的比较，可以推断出由 成像系统引起的失真。尤其是当第一区域中的病理变化小于通常在成像系统中 出现并用成像图像绘制出的失真时，可以通过将两张成像图像之一 中的第 一 区 域与另 一张成像图像中的第一区域进行比较来以较大的精度检测出由成像系统 引起的失真。
在本发明的方法中这通过确定针对第 一 区域的配准来进4亍，从而通过该配 准将第 一成像图像的第 一 区域与第二成像图像的第 一 区域相互对准。这意味着 确定使得第一成像图像的第一区域和第二成像图像的第一区域相互尽可能重叠 的变换。
在通过确定配准将一张成像图像中的第 一 区域与另 一张成像图像中的第 一区域相互对准并且由此在第 一 区域中识别出由成像系统引起的失真之后，将 针对第一区域的配准扩大为扩张配准，从而通过该扩展的配准一起采集第二区 域。
通过这种方式基本上去掉了第二区域中的由成像系统引起的扭曲，尤其是 失真，从而可以将剩余的差异推断为病理变化。仅当由成像系统引起的扭曲在 成像图像中局部强烈变化时，该方法才会遇到其限制，因为在这种情况下由成 像系统引起的第二区域中的扭曲不能从第 一 区域的扭曲中得到。但是强烈局部 变化的扭曲在目前采用的成像系统通常不会出现。
现在将扩展的配准用于变换两张成像图像之一 。
在对两张成像图像之一进行了变换之后，可以按照不同方式显示和/或确定 两张成像图像之间目前主要存在于第二区域中的剩余差异。按照特别简单的方 式向用户并排显示两张成像图像，即经过变换的成像图像和另一张成像图像， 从而用户可以直观地通过两张成像图像之间的视觉比较而识别出该差异。但是
通过自动或半自动地确定该差异，并且显示所确定的差异，用户也可以在确定 该差异时得到支持。例如可以将两张成像图像相减，从而由此获得两张成像图 像之间的差异。然后所获得的差异例如以减影图像的形式显示出来，或者可以 在所显示的成像图像中标记出所获得的差异。
通常接着将所获得的差异或者两张成像图像一其中有一张经过了变换一 存储在存储介质中。
优选地，所述方法用于第 一 区域基本上没有被医学疾病图像中的病理变化 覆盖的情况。如果第一区域没有被病理变化覆盖，则该区域基本上是恒定的， 由此如果在两张于不同时刻拍摄的成像图像中在第一区域的显示中出现差异， 则该差异可以推断是生成该成像图像的成像系统的原因。现在可以特别准确地 通过将第 一成像图像中的第 一 区域与另 一张成像图像中的第 一 区域对准来确定 由成像系统引起的差异。
优选地，在确定针对第一区域的配准时只确定针对第一区域的配准。通过 这种方式可以更简单、快速和有效地确定该配准，因为只将一个成像图像的第 一区域与另 一 张成像图像的第 一 区域相关。
优选地，针对第一区域的配准是非线性的配准。非线性配准的采用可以更 为准确地采集由成像系统引起的差异，从而提高该方法的精度。
在本方法的优选实施方式中，在扩展所述配准时扩大成像在一张成像图像 中的第一区域，并且将针对第一区域的配准扩展到扩大后的第一区域。由此以 两个分步骤完成了配准的扩展。在此第一分步骤只包括第一区域的扩大，在第 二分步骤中才在扩大后的第 一 区域中确定针对图像点的配准。分为两个子步骤 具有以下优点，即可以更为简单地实施各个子步骤，因为对此存在公知的算法。
第一区域的扩大例如可以通过例如以简单方式用于扩张的形态学图像处 理算子、如所谓的扩张滤波器使第一区域膨胀来实现。尤其是对于凸起的第一 区域来说，还可以通过放大比例来扩大。
在将该配准扩大到经过扩大的第一区域时，通常可以采用不同的方法。例 如可以将经过扩大的第一区域中的一个图像点对应于该第一区域中最靠近该图 像点的 一个图像点，然后将已知针对该第 一 区域的最近图像点的配准转用于经 过扩大的区域中的图像点。在另一种方案中，第一区域的配准可以外推到经过 扩大的第一区域上。
优选地，在扩展所述配准时利用随后的配准扩展来迭代地执行第一区域的
扩大。通过配准的迭代扩展，可以用随后的配准扩展更简单地实施扩大第一区 域的子步骤，因为不必在唯一的一个步骤中保证也一起采集第二区域的配准扩 展。第二区域可以通过该简单方式由连续和重复采用所述扩大来采集。
在本方法的优选实施方式中，在确定针对第 一 区域的配准之前分割第 一 区 域。通过这种方式可以借助事先进行的分割帮助找到将两张成像图像中的第一 区域相互对准的配准，即确定将两张成像图像中的第一区域尽可能置为重叠的 变换。
优选地，两张成像图像是利用磁共振设备或利用计算机断层造影设备完成 的器官系统的三维成像图像。利用三维成像图像，可以更为精确和有效地执行
第一区域的配准和该配准的扩展，因为可以考虑全部三维并且可以;险测和平衡 所有三个空间方向上的失真。尤其是可能由于磁场不均匀性而出现在磁共振图 像中的失真，可以通过这种有效的方式检测出来，并在分析图像时加以考虑。在本方法的优选实施方式中，所述器官系统是患者的大脑。因此该方法对 大脑成像图像的应用特别有利，因为在大脑中待检测的变化通常程度很小，但 是如果存在则具有很高的诊断说服力。
优选地，第一区域是脑白质，第二区域是脑灰质。通过这种方式可以检测 出神经退化过程中出现的微小变化，例如在阿兹海默症(Morbus Alzheimer)中出 现在脑灰质中的变化。
优选地，拍摄所述两张相对应的成像图像的时刻具有1到12个月、优选3 到6个月的间隔。这种时间间隔对于诊断大脑的很多疾病、尤其是对于神经退 化疾病来说是很有利的。
本发明的医学成像系统与计算单元连接，该机算单元实施为用于执行上述 方法。
下面借助附图详细解释本发明以及优选实施方式，但是不限于此。 其中
图la和图lb分别示出大脑的在不同时刻拍摄的横断面， 图2a和2b分别示出脑标记属于该横断面的区域以及两个区域的重叠， 图3示出图2a的一个片段，其中除了脑标记的区域之外还示出两个扩张 区域，
图4a和4b示出借助所确定的配准对第二横断面的变换， 图5示出经过变换的第一横断面和第二横断面的重叠显示，
具体实施例方式
借助图1至图5的人脑的成像图像来说明和解释本发明的方法。图6示出 各个方法步骤的示意概貌。
图la和图lb示出大脑在同一位置的第一横断面l和第二横断面3。两个 横断面在此分别来自 一个利用MRT设备绘制的3D立体数据组。在该3D立体 数据组中明显可相互分隔地对脑皮质5—或脑灰质一以及脑标记7—或脑白质 成像。
两个3D立体数据组在此是在不同时刻完成的。为此而绘制这种3D立体数
问题。在这种情况下优选以3到6个月的时间间隔来绘制这两个3D立体数据 组。
第二横断面3与第一横断面1相比具有不规则的几何失真。在图la和图 lb中通过所示出的矩形9来表示该失真。这种矩形9在此只用于表示否则很难 指明和看见的失真，而且在通常给出的成像图像中看不见。
第二横断面3的失真可以由不同原因引起。在磁共振设备中，该失真是由 于为了完成成像图像而施加的磁场在运行期间出现的微小变化而引起。在这个 过程中出现的失真在此可以位于毫米的范围内。但是还有其它原因引起该失真， 例如对医学成像系统的有错或更改的校准。
如果用户想要分析两个横断面，则他面临可以猜测颈叶区域11中脑皮质 的变化的问题，因为在此在第二横断面3中的脑皮质看起来比第一横断面1中 的脑皮质更窄。但是最后用户不能确定，因为该变化至少部分要归于是失真。
但是由于阿兹海默痴呆症是一种该退化过程会导致主要在脑皮质5中的体 积变化，而脑标记7在形态学上基本上保持不变的疾病，因此通过脑标记7的 比较来确定该失真。下面借助图2a和图2b解释该步骤。
首先分别从第一和第二横断面1、 3中提取出脑标记区域，下面称为第一 脑标记区域13和第二脑标记区域15。由于在3D立体数据组中脑标记7明显地 以脑皮质5为界，因此这例如可以借助基于各个体素的特殊信号强度而识别出
脑标记7的方法来完成。但是还可以釆用其他自动或半自动的常用的分割算法。 例如，用户可以点击脑标记7的一个中心区域，由此触发分割算法，该算法从 选择的起始点开始找寻关联的脑标记7的轮廓。
第一脑标记区域13和第二脑标记区域15具有由失真引起的差异。这在图 2b中示出在第二脑标记区域15中，因为在此第一脑标记区域13用点划线重叠 地示出，从而该失真清楚可见。
通过将第一脑标记区域13与第二脑标记区域15相比较确定将第一和第二 脑标记区域13、 15相互转化的配准。这意味着借助针对第一脑标记区域13的 每个像素的配准可以知道将该像素与第二脑标记区域15的对应像素相对应以 及反过来对应的变换向量。
在此，优选所选择的配准类型与医学成像系统匹配，并且表示配准精度和 获得该配准所需要的计算时间之间的折衷。对于例如通常只出现仿射失真的成 像系统，只确定将两个区域尽可能相关的仿射配准就够了。对于其它成像系统， 如由于磁场非均匀性而导致图像区域主要是非线性失真的MRT设备，将这一 点用于获得匹配的非线性配准。匹配的配准还可以连续获得，为此例如首先获 得刚性或仿射配准，然后以此为基础获得非线性配准。
通过配准两个脑标记区域13、 15,现在知道第一脑标记区域13的图像点 如何与第二脑标记区域15关联，或者换句话说对于一个脑标记区域的每个图像 点来说都知道一个变换向量，该向量将该图像点映射为另 一个脑标记区域的对 应图像点。
在找到将第一和第二脑标记区域13、 15相互转换的配准之后，将该配准 扩展到其它区域。为此在第一步骤中扩张第一(或第二)脑标记区域13。该扩展 例如可以通过简单方式借助用于扩张的形态学的图像处理算子来进行。由于该 配准只是针对第一和第二脑标记区域13、 15而导致扩展后区域的还没有通过配 准采集的像素，可以在第二步骤中配设变换向量，如下面借助图3所解释的。
图3示出第一脑标记区域13中对应于图2的矩形III的一个片段。除了第 一脑标记区域13之外还示出第一扩张区域17。第一扩张区域17在此由扩张第 一脑标记区域13而得到。
第一扩张区域17的像素全部位于第一脑标记区域13附近。现在对位于第 一扩张区域17中以及还没有通过配准采集的的像素21赋予一个变换向量。例 如可以对该像素21赋予在第一脑标记区域13中靠得最近的像素23的变换向 量。替代式地，还可以赋予在脑标记13中多个靠近像素21的像素的变换向量
的平均值(必要时加权)。另外，用于位于第一扩张区域17中的像素21的变换 向量也可以通过外推用于第一脑标记区域13的像素的变换向量来获得。在任何 情况下都将配准从第一脑标记区域13扩大或外推到第一扩张区域17。
在将该配准扩大到第一扩张区域17之后，可以重新迭代地执行该扩大， 从而该扩大还覆盖第二扩张区域19和其它后续区域，直到也覆盖对脑皮层5 成像的区域为止。
经过扩大的配准具有以下优点，即虽然该配准将用于获得该配准的区域一 在这种情况下即脑标记7—相互转换，从而尽可能覆盖这些区域，但是其余区 域一在这种情况下即脑皮层5—通过该扩大的配准只能发生使得由失真引起的 扭曲得到平衡这样的改变。这种由其余区域的解剖特性发生变化而引起的改变，
不能通过该配准得到平衡，从而在采用该配准时还要获得该信息。
在此所示的配准通过第一和第二横断面1、 3的第一和第二脑标记区域13、 15之间的比较来获得，即通过比较两张二维成像图像来获得。但是配准的获得 不限于二维截面图像。脑标记7的提取还可以对整个3D立体数据组进行，从 而分别获知脑标记7的三维结构。然后该配准还通过比较脑标记7的两个三维 结构来进行。
采用三维结构来获得配准具有以下优点，即通过这种方式还可以采集在各 个横断面的排列方向上出现的失真。由此产生更为精确的、分别将一个3D立 体数据组的脑标记和其它3D立体数据组的脑标记相互转换的配准。
在扩大了配准之后，两个横断面之一一在这种情况下是第二横断面3 —用 该配准变换，如图4所示。如果现在将经过变换的第二横断面25与第一横断面 1比较，则作为该配准的起点的第一和第二脑标记区域13、 15基本上重合。相 反，对于脑皮层5来说通过该配准只能基本上校正由失真引起的差异，而由于 病理变化引起的差异大部分还仍然获得。
通过特别简单的方式例如可以向用户并排显示经过变换的第二横断面25 和第一横断面l。通过平衡由成像系统引起的失真，用户可以将其视线集中在 两张成像图像之间可以推断是病理变化的差异上。
病理变化通过经过变换的第二横断面25与第一横断面1的重叠显示29而 清楚可见，如图5所示。两个横断面之间的差异27现在可以简单地可视化，例 如通过将第一横断面l减去经过变换的第二横断面25或者通过有颜色地标识出
差异27。同样在获得该差异后接下来采用例如通过体积测量量化地检测该变化 程度的分析方法。
但是，在此在图1至图5中借助脑部MR图像解释的方法步骤不限于在图 1至图5中示出的特殊器官系统、特殊成像模态或者特殊医学课题。专业人员 可以在其它成像模态中以及对其它医学课题采用该方法，其中器官系统具有在 一张成像图像中有差异显示的第一和第二区域，而且病理变化主要表现在这两 个区域之一 中，而另 一个区域基本未受到病理变化的影响并保持恒定。
该方法的起点是器官系统的第 一成像图像31和与该第 一成像图像31对应 的第二成像图像33，它们分别是用医学成像系统35在不同时刻拍摄的。在此 该器官系统包括分别在两张成像图像31、 33有差异显示的第一区域和第二区 域。
通过本方法在这两张成像图像31、 33上检测可能出现的病理变化，这种 变化属于医学疾病图像。由该病理变化知道该病历变化仅在该器官系统的第二 区域中引起变化，而第 一 区域没有被该疾病图像中的病理变化覆盖。
在第一方法步骤37中确定针对第一区域的配准，从而通过该配准将第一 成像图像31的第一区域和第二成像图像33的第一区域相互对准。
在第二方法步骤39中扩大只用于将两张成像图像中的第一区域相互对准 的配准，从而通过该配准的扩展也一起采集了两张成像图像的第二区域。
在第三方法步骤41中，借助经过扩展的配准变换两张成像图像之一，从 而在第四方法步骤43中通过比较经过变换的成像图像的第二区域与另一张成 像图像的第二区域来自动或半自动获得两张成像图像在第二区域上的差异。
可替换地，在第五方法步骤45中向用户并排显示经过变换的成像图像和 另一张成像图像，从而用户可以通过肉眼观察两张图像而获得两张成像图像之 间的差异。
在第六方法步骤47中，将该方法的结果、即所获得的差异和/或经过变换 的成像图像和另 一张成像图像存储在存储介质中。
优选的，借助与所述医学成像系统连接的计算单元49来执行本方法，因 为这种计算单元49通常具有优选可以用于在此所执行方法的图像分析功能和 图像显示功能。
权利要求
1.一种用于分析至少一个器官系统的两张相对应的、在不同时刻拍摄的成像图像(1，3；31，33)以确定医学疾病图像中的病理变化的方法，其中所述至少一个器官系统具有第一区域(7，13，15)和在成像图像(1，3；31，33)中与该第一区域(7，13，15)不同的第二区域(5)，其中第二区域(5)比第一区域(7，13，15)更为强烈地被该医学疾病图像中的病理变化覆盖，该方法包括以下步骤-确定(37)针对第一区域(7，13，15)的配准，从而通过该配准将第一成像图像(1；31)的第一区域(7，13，15)与第二成像图像(3；33)的第一区域(7，13，15)相互对准，-将针对第一区域(7，13，15)的配准扩大(39)为扩展的配准，从而通过该扩展的配准一起采集第二区域(5)，-借助所述扩展的配准变换(41)两张成像图像之一，-显示(45)经过变换的成像图像和另一张成像图像，和/或通过比较经过变换的成像图像的第二区域(5)与另一张成像图像的第二区域(5)来确定(43)两张成像图像(1，3；31，33)在第二区域上的差异(27)。
9. 根据权利要求1至8中任一项所述的方法，其特征在于，所述器官系 统是患者的大脑。
10. 根据权利要求9所述的方法，其特征在于，所述第一区域(7， 13, 15) 是脑白质，第二区域(5)是脑灰质。
11. 根据权利要求9或10所述的方法，其特征在于，拍摄所述两张相对 应的成像图像(l， 3; 31, 33)的时刻具有3到6个月的间隔。
12. —种具有计算单元(49)的医学成像系统(35)，该机算单元实施为用于实 施根据权利要求1至11中任一项所述的方法。
全文摘要
本发明涉及一种用于分析至少一个器官系统的两张相对应的、在不同时刻拍摄的成像图像(1，3)以确定医学疾病图像中的病理变化的方法。所述器官系统具有在成像图像中不同的第一区域(7，13，15)和第二区域(5)，其中第二区域比第一区域更为强烈地被该医学疾病图像中的病理变化覆盖。通过确定针对第一区域失真。在扩展了配准从而通过该扩展的配准一起采集第二区域之后，通过借助该扩展的配准变换两张成像图像之一来获得、显示和/或存储两张成像图像(1，3)在第二区域上的差异(27)。此外本发明还涉及一种具有执行该方法的计算单元(49)的医学成像系统(35)。
文档编号A61B5/055GK101097628SQ20071012707
公开日2008年1月2日 申请日期2007年6月28日 优先权日2006年6月28日
发明者阿恩·利特曼 申请人:西门子公司