专利名称：基于一致性度量的图像配准的利记博彩app
基于一致性度量的图像配准本发明大体涉及基于一致性度量来配准医疗成像图像，并且本发明尤其应用于计 算机断层摄影(CT)。然而，其也适于其他医疗成像应用和非医疗成像应用。诸如计算机断层摄影(CT)、磁共振(MR)、超声(US)、单光子发射计算机断层摄影 (SPECT)、正电子发射断层摄影(PET)和χ射线的医疗成像形式在诸如癌症的疾病诊断中扮 演重要角色。例如，可以使用它们非侵入地获得指示身体中生理组织的信息，并且可以使用 这些信息帮助确定肿瘤是良性的还是恶性的。典型地，这种非侵入技术比诸如活组织检查 的侵入技术风险更小并且成本更低。另外，对于相对较小的肿瘤，诸如10毫米或者更小的 肿瘤，确定活组织检查针是否击中肿瘤可能相对困难。更具体地，可以使用诸如CT图像的图像执行不同的诊断。通过举例的方式，可以 使用均包括指示同一个肿瘤、但是从在诸如间隔一个月至六个月的不同时刻所采集的数据 生成的信息的两幅CT图像，通过将第一图像中肿瘤的尺寸与第二图像中同一个肿瘤的尺 寸进行比较获取肿瘤随着时间的生长。概括地说，肿瘤尺寸的增大大于预置阈值(例如， 20% )指示该肿瘤是恶性的，而不生长或者小于阈值的生长指示该肿瘤是良性的。不幸的是，当对这些图像进行比较时，由于患者姿态的不同，诸如肺的一些器官可 能不在两幅图像中的相同位置。因而，两幅图像之间的空间配准可能有问题。例如，为了找 到显示用于与来自第一次扫描的第一图像进行比较的肿瘤的图像，临床医生可能必须手动 查看采用在第二次扫描中所采集的数据生成的、第二个图像集合中的许多图像(例如，200 幅或者更多)。即使在空间配准之后，在肺内部诸如肿瘤的结构也可能由于呼吸状态的不同 不在同一个位置。采用自动肿瘤生长评估，将肿瘤在第一图像中的位置自动匹配到第二图像中的位 置。使用弹性配准技术，第一图像到第二图像的完全配准是可能的，即使可能存在患者姿态 差异、解剖区间变化、以及呼吸状态差异，弹性配准技术也允许成像扭曲。不幸的是，弹性配 准可能扭曲图像中的肿瘤，并且因此可能改变肿瘤的尺寸，当执行差异诊断时这是有问题 的。本申请的各方面解决了上面所提出的问题以及其他问题。根据一个方面，一种方法包括将第一图像的第一子部分与第二图像的对应的第二 子部分配准，将第二图像的第二子部分与第一图像的对应的第三子部分配准。第一子部分 包含第一感兴趣对象，并且第三子部分包含第三感兴趣对象。该方法还包括当第一和第三 感兴趣对象充分相似时减小第一子部分的尺寸。该方法还包括重复配准第一子部分、配准 第二子部分、以及减小第一子部分的尺寸直到第一和第三对象不充分相似为止的步骤。根据另一个方面，一种图像配准系统包括将第一图像的第一感兴趣体积与第二图 像的对应的第二感兴趣体积配准和将第二图像的第二感兴趣体积与第一图像的对应的第 三感兴趣体积反向配准的配准部件。第一感兴趣体积包括第一感兴趣对象，并且第三感兴 趣体积包括第三感兴趣对象。该系统还包括确定第一感兴趣对象和第三感兴趣对象之间一 致性值的一致性确定部件。一致性值指示第一感兴趣对象和第三感兴趣对象之间的相似 性。该系统还包括基于一致性值确定第一感兴趣体积的尺寸是否减小的判决部件。
根据另一个方面，一种方法包括迭代地减小第一图像中感兴趣体积的尺寸直到第 一和第二图像之间感兴趣体积的前向和反向配准之间的一致性误差变得大于一致性误差 阈值为止。该方法还包括呈现第一图像的第一感兴趣体积和第二图像的第二感兴趣体积， 其中，第一感兴趣体积和第二感兴趣体积的尺寸与在前向配准和反向配准之间的一致性误 差小于一致性误差阈值时的感兴趣体积的尺寸对应。本发明可以以各种部件和部件的布置，和各种步骤和步骤的布置而变得明显。附 图仅仅是为了图示说明优选实施例的目的，而不应被解读为限制本发明。
图1图示说明了一种医疗成像装置；图2图示说明了示例性图像配准部件；图3图示说明了一种方法；图4图示说明了在两个图像中VOI的初始配准；图5图示说明了在经配准的VOI之间的一致性检查；图6图示说明了在两幅图像中VOI的精确配准。首先参考图1，计算机断层摄影(CT)扫描仪100包括静止台架102，从扫描期间其 一般是静止的意义上来讲其是静止的。然而，可以配置该静止台架102倾斜以及/或者以 其它方式移动。计算机断层摄影(CT)系统100还包括可旋转耦合到静止台架102的旋转台架 104。该旋转台架104围绕检查区域106关于纵轴或者ζ轴108旋转。诸如χ射线管的辐射源110由旋转台架104支撑并与旋转台架104—起围绕检查 区域106旋转。辐射源110 —般发射横贯检查区域106的扇形、楔形、或者锥形辐射。也可 以考虑使用第四代系统。辐射敏感探测器阵列112探测由辐射源110发射的横贯检查区域106的光子，并 且生成指示所探测辐射的投影数据。所图示说明的辐射敏感探测器阵列112包括在ζ轴方 向上延伸的多行辐射敏感光传感器、以及在横向方向上延伸的多列辐射敏感光传感器。也 预期单行探测器阵列结构。重建器114重建来自探测器的投影数据，从而生成指示患者内部解剖体的体积图 像数据。图像处理器116对由重建器114所生成的体积图像数据进行处理，以便以人类可 读形式显示。诸如卧榻的患者支撑物118在检查区域106中支撑患者。患者支撑物118可与旋 转台架104的旋转协同沿着ζ轴108移动，从而促成螺旋状、轴向、或者其他期望扫描轨道。通用计算系统120作为操作员控制台。操作员控制台120包括诸如显示器和/或 打印机的人类可读输出设备以及诸如键盘和/或鼠标的输入设备。驻留在控制台120上的 软件允许操作者控制系统100的操作，例如，通过允许操作员选择扫描协议、发起并且终止 扫描、观看和/或操控体积图像数据、和/或以其它方式与系统100交互。可以使用存储部件122存储由重建器114生成的体积图像数据和/或由图像处理 器116生成的一个或多个图像。使用配准系统124基于一致性度量对在不同时刻采集的数据进行配准。一致性确 定部件126确定一致性度量。如下面更详细描述的，配准系统124执行迭代配准，其中，基于一致性度量使所配准图像中诸如感兴趣区域(ROI)或者感兴趣体积(VOI)的、包含感兴 趣对象的图像子部分的尺寸最优化。在一个实例中，最优化确保图像之间感兴趣对象的一 致明确匹配，同时减少了周围的外部解剖结构。图2进一步图示说明了配准部件124。对于下列讨论，注意到，使用在第一时间采 集的数据生成第一图像，并且使用在第二时间采集的数据生成第二图像，该第二时间可能 是第一时间之后的数周、数月等。感兴趣对象识别器202识别第一图像中的第一感兴趣对象。在一个实例中，感兴 趣对象识别器202基于用户输入识别第一感兴趣对象。例如，用户可以使用鼠标、键盘和/ 或其他输入设备在第一图像中选择第一感兴趣对象。感兴趣体积生成器204在第一图像中在所识别的第一感兴趣对象周围生成第一 感兴趣体积(VOI)。在一个实例中，第一 VOI的初始形状和尺寸是预配置的。合适的形状可 能取决于感兴趣对象和/或感兴趣对象在身体内中的位置。将初始尺寸设置得充分大以包 括对于将VOI明确匹配到另一幅图像足够的背景信息。用户定义的形状和/或尺寸也是预 期的。另外，VOI可以限制于一个或多个特定器官。前向配准部件206将第一图像中的第一感兴趣对象和第一 VOI与第二图像中的对 应的第二感兴趣对象和第二 VOI配准。第二感兴趣对象表示第二图像中被认为是对应于第 一图像中第一感兴趣对象的结构。第一 VOI提供了有助于匹配第一和第二感兴趣对象的背
足fn息。在所说明的示例中使用诸如刚性配准的体积保持配准。通过最优化第一 VOI中的 第一图像值和第二 VOI中的第二图像值之间的相似性量度来计算刚性配准。这可以包括对 于每个体素计算第一图像和第二图像之间的差异。为了最优化，改变第二VOI的位置、取向 和/或比例，直到达到第一感兴趣对象和第二感兴趣对象之间的最佳相似性为止。在所说明的实施例中，仅计算第一和第二 VOI区域之间的相似性量度。相似性量 度可以基于相关性、均方根偏差、交互信息等。为了平滑，可以根据(例如，高斯)加权贡 献对相似性量度求和，其中权重从VOI中心到外围递减。最优化技术可以是穷举、随机、高 斯-牛顿等。反向配准部件208将第二图像中的第二感兴趣对象与第一图像中对应的第三感 兴趣对象配准，第二 VOI提供了有助于这种匹配的背景信息。第三感兴趣对象表示第一图 像中被认为是对应于第二图像中第二感兴趣对象的结构，并且因此，表示第一图像中被认 为是第一感兴趣对象的结构。采用明确配准，所得到的第三感兴趣对象应该与第一感兴趣 对象充分相符。第二VOI提供了有助于匹配第二和第三感兴趣对象的背景信息。再次，使 用诸如刚性配准的体积保持配准。一致性确定部件126测量第一图像中第一感兴趣对象和第三感兴趣对象之间的 一致性。在所说明的实施例中，通过确定第一感兴趣对象和第三感兴趣对象之间诸如欧几 里德距离的距离来计算一致性度量。等式1说明了用于计算第一和第三感兴趣对象之间、 通过体素总数目归一化的欧几里德距离的示例算法。 其中，D代表欧几里德距离，i代表体素总数目，并且X、y和ζ代表体素的坐标。
注意到，在严格意义上，一致性度量不保证能够识别错误匹配。理论上，可能出现 错误的、但反向/前向一致的匹配。然而，一致性度量可以提供对错误匹配的充足指示。判决部件210基于一致性度量确定VOI的尺寸是否是最优化的，该最优化是在该 VOI 一般是可以将第三感兴趣对象与第一感兴趣对象明确匹配的最小尺寸的意义上的。在 一个实例中，判决部件210基于关于一致性度量是大于还是小于预置一致性阈值来进行判 决。如果一致性度量小于阈值，判决部件210就调用感兴趣区域生成器204缩减或者 减小第一 VOI的尺寸。在一个实例中，按照诸如30%的预置百分比减少VOI的体积。在另 一个实例中，按照预置体积减少VOI的体积。也预期其他减少方式。然而，当一致性度量变 得大于阈值时，判决部件210将其一致性度量小于阈值的最后一个VOI识别为最优VOI。存储部件212存储V0I。例如，当一致性度量小于阈值时，将VOI存储在存储部件 212中。同样，如果用于下一个VOI的一致性度量大于一致性阈值，当前VOI是可利用的。讨论各种变化。在所说明的实施例中，用户识别第一感兴趣对象。在可选实施例中，感兴趣对象识 别器202基于期望的组织类型的一个或多个特征识别第一图像中的候选感兴趣对象。这种 识别可以基于第一图像中的灰度级值和/或其他特征。在所说明的实施例中，在感兴趣对象周围生成感兴趣VOI。在可选实施例中，使用 感兴趣区域(ROI)。同样，初始ROI尺寸充足大，并且最终ROI尺寸使得它提供一致明确匹 配、同时减少外部解剖结构。VOI可以是不同形状的。例如，合适的形状包括但是不限于立方体、圆柱体、球体、 椭圆体、和/或其他形状。当使用ROI而不是VOI时，可以使用相对应的二维形状。在所说明的实施例中，将初始VOI设置为相对较大，并且随后减少尺寸直到一致 性度量变得大于一致性阈值为止。在另一个实施例中，可以将初始VOI (或R0I)设置为相 对小，并且随后增大其尺寸直到一致性度量变得小于一致性阈值为止。在所说明的实施例中，使用刚性配准来保持在VOI中的结构的尺寸。在另一个实 施例中，使用弹性配准。在一个实例中，与弹性配准一起使用一系列约束，使得不扭曲感兴 趣对象。应该理解，配准系统124可以是医疗成像系统(如所示)的部分或者是与医疗成 像系统分离的工作站的部分。现在结合图3-6描述操作。在302处，在第一图像中识别第一感兴趣对象。如上所述，这可以经手动或者自动 方法实现。在304处，在第一感兴趣对象周围生成第一 VOI (或者R0I)。对于所选择的感兴趣 对象，为了为第一 VOI与另一幅图像中的VOI的明确配准捕获足够的解剖背景，算法预见以 相对大的固定尺寸的VOI开始。图4示出了具有第一感兴趣对象402和第一 VOI 404的示例性第一图像400。在 本示例中，第一 VOI 404是立方体形的。返回图3，在306处，在第二图像中匹配对应的第二 VOI和感兴趣对象。返回图4，平移、旋转和/或缩放对应的VOI 406，直到基于相似性量度第二图像408中的VOI 406与第一图像400中的VOI 404最相似为止。匹配算法对第二图像408中 最可能的对应感兴趣对象410进行识别。如上所述，第二图像408中VOI 402的配准可以 是刚性的以保持解剖体的形状并且不影响肿瘤生长评估。返回图3，在308处执行反向配准。对于反向配准，基于第二图像中的第二 VOI和 感兴趣对象在第一图像中识别第三VOI和感兴趣对象。如上所述，第三感兴趣对象应该与 第一感兴趣对象匹配。在310处，确定第一感兴趣对象和第三感兴趣对象之间的一致性度量或者误差。 在所说明的实施例中，一致性度量是第一和第三感兴趣对象之间的欧几里德距离，其提供 了对第一和第三感兴趣对象之间匹配的客观量度。图5示出了第一感兴趣对象402 ( “？1”)、第二感兴趣对象410( “P2”)和第三感 兴趣对象(“P3”)、以及Pl和P3之间的距离(“D”)。返回图3，在312处，将一致性度量与预置的阈值进行比较。如果一致性度量小于 阈值，那么在314处减小第一 VOI的尺寸，并且再次执行动作306至312。当一致性度量变 得大于阈值时，则在316处，选择最后的在其处一致性度量小于阈值的VOI尺寸。这样，VOI尺寸逐渐或者迭代地减小，直到丢失了如此多背景信息以至于匹配不再 互逆并且因此一致性误差变得过大为止。换言之，如果VOI尺寸变得如此小以至于VOI的 匹配开始变得不明确，那么前向和后向配准之间的匹配将变得不一致。图6示出了一个示例，其中，基于一致性度量迭代减小初始VOI尺寸602经过VOI 尺寸604和606，直到获得最优尺寸608为止。注意到，在图6中，采用圆柱形V0I，而图4 和图5中VOI的形状是立方体。返回图3，在316处，将产生最后一个仍然可接受的一致性误差的VOI呈现给用户， 该用户可以例如是放射线学者或者自动生长评估算法。这样，当匹配一致性误差增大到可 容忍水平之上时，该过程中止，并且呈现对应于仍产生一致匹配的最佳VOI的匹配解决方案。总之，期望最终VOI的尺寸尽可能小，这是因为最重要的任务是感兴趣对象的匹 配，该感兴趣对象不应该受诸如肺壁、肋骨、血管或者器官的其他结构影响，其中该其他结 构可能在两个图像体积之间已经相对于感兴趣对象改变了位置。另一方面，如果VOI尺寸 选择得过小，那么因为丢失了过多背景信息，所以不再可能得到两个VOI的明确匹配。通常 的根本想法是VOI尺寸应该足够大以建立一致明确匹配，但是足够小以尽可能精确地将所 选结构的位置匹配到另一幅图像中的对应位置，而不会被远处的解剖结构扰乱。在这里结合计算机断层摄影医疗成像应用对实施例进行描述。然而，应该理解，可 以附加地或者备选地与期望在其中捕获辐射频谱方面的其他医疗成像应用和/或非医疗 成像应用一起使用本发明。参考优选实施例对本发明进行了描述。在阅读和理解了前述详细说明的基础上， 可以对其进行修改和改变。旨在将本发明解释为包括所有这些修改和改变只要他们落在所 附权利要求或者其等同替代范围内。
权利要求
一种方法，包括将第一图像的第一子部分与第二图像的对应的第二子部分配准，其中，所述第一子部分包含第一感兴趣对象；将所述第二图像的所述第二子部分与所述第一图像的对应的第三子部分配准，其中，所述第三子部分包含第三感兴趣对象；当所述第一感兴趣对象和所述第三感兴趣对象充分相似时减小所述第一子部分的尺寸；以及重复配准所述第一子部分、配准所述第二子部分、以及减小所述第一子部分的尺寸的步骤，直到所述第一对象和所述第三对象不充分相似为止。
2.如权利要求1所述的方法，还包括确定指示所述第一感兴趣对象和所述第三感兴趣 对象之间的相似性的一致性度量，并且基于所述一致性度量重复配准所述第一子部分、配 准所述第二子部分、以及减小所述第一子部分的尺寸的步骤。
3.如权利要求2所述的方法，还包括通过计算所述第一感兴趣对象和所述第三感兴趣 对象之间的欧几里德距离来计算所述一致性度量。
4.如权利要求3所述的方法，还包括将所述一致性度量与预置的一致性阈值进行比 较，以确定所述第一感兴趣对象和所述第三感兴趣对象是否充分相似。
5.如权利要求1所述的方法，其中，所述第二子部分包含第二感兴趣对象，并且所述方 法还包括呈现对应于所述第一感兴趣对象和所述第三感兴趣对象于其处充分相似的所述 第一子部分和所述第二子部分的所述第一子部分和所述第二子部分。
6.如权利要求5所述的方法，其中，所呈现的第一子部分和第二子部分具有包括足够 解剖内容以将所述第一子部分与所述第二子部分明确配准的尺寸。
7.如权利要求1所述的方法，其中，所述第一子部分和所述第三子部分是感兴趣体积。
8.如权利要求1所述的方法，其中，按照固定百分比减小所述第一子部分的尺寸。
9.如权利要求1所述的方法，还包括使用体积保持配准，其中，仅采用平移、取向和缩 放中的一个或多个来配准各子部分。
10.如权利要求1所述的方法，其中，所述第二子部分包含第二感兴趣对象，并且所述 方法还包括最优化所述第一子部分和所述第二子部分中的图像值与所述第二子部分和所 述第三子部分中的图像值之间的相似性量度。
11.如权利要求1所述的方法，其中，所述感兴趣对象是肿瘤。
12.—种图像配准系统，包括配准部件(124、206、208)，其将第一图像的第一感兴趣体积与第二图像的对应的第二 感兴趣体积配准，并且将所述第二图像的所述第二感兴趣体积与所述第一图像的对应的第 三感兴趣体积反向配准，其中，所述第一感兴趣体积包括第一感兴趣对象，并且所述第三感 兴趣体积包括第三感兴趣对象；一致性确定部件(126)，其确定所述第一感兴趣对象和所述第三感兴趣对象之间的一 致性值，其中，所述一致性值指示所述第一感兴趣对象和所述第三感兴趣对象之间的相似 性；以及判决部件(210)，其基于所述一致性值确定是否减小所述第一感兴趣体积的尺寸。
13.如权利要求12所述的系统，其中，所述判决部件(210)将所述一致性值与一致性阈值进行比较，以确定是否减小所述第一感兴趣体积的尺寸。
14.如权利要求12所述的系统，其中，所述一致性值是所述第一感兴趣对象和所述第 三感兴趣对象之间的欧几里德距离，并且通过用于确定所述一致性值的体素总数目对所述 欧几里德距离进行归一化。
15.如权利要求1所述的方法，其中，按照固定体积减小所述第一感兴趣体积的尺寸。
16.如权利要求1所述的方法，其中，所述一致性值基于相关性、均方根偏差、或者交互 信息中的至少一个以配准所述感兴趣体积。
17.如权利要求1所述的方法，还包括来自所述感兴趣体积的加权贡献，所述感兴趣体 积具有从所述感兴趣体积的中心到所述感兴趣体积的外围递减的权重。
18.如权利要求1所述的方法，其中，所述配准部件(124、206、208)采用穷举配准、随机 配准、或高斯-牛顿配准之一。
19.一种方法，包括迭代地减小第一图像中感兴趣体积的尺寸直到在第一图像和第二图像之间的所述感 兴趣体积的前向配准和反向配准之间的一致性误差变得大于一致性误差阈值为止；以及呈现所述第一图像的第一感兴趣体积和所述第二图像的第二感兴趣体积，其中，所述 第一感兴趣体积和所述第二感兴趣体积的尺寸与在所述前向配准和所述反向配准之间的 所述一致性误差小于所述一致性误差阈值处的所述感兴趣体积的尺寸对应。
20.如权利要求19所述的方法，其中，所述第一感兴趣体积和所述第二感兴趣体积包 括用于将所述第一感兴趣体积与所述第二感兴趣体积明确配准以使得所述第一感兴趣体 积中的第一感兴趣对象充分映射到所述第二感兴趣体积中的第二感兴趣对象的足够解剖 内容，其中，所述第一感兴趣对象和所述第二感兴趣对象对应于同一解剖结构。全文摘要
一种方法，包括将第一图像的第一子部分与第二图像的对应的第二子部分配准，并且将第二图像的第二子部分与第一图像的对应的第三子部分配准。第一子部分包含第一感兴趣对象，并且第三子部分包含第三感兴趣对象。该方法还包括当第一和第三感兴趣对象充分相似时减小第一子部分的尺寸。该方法还包括重复配准第一子部分、配准第二子部分、以及减小第一子部分的尺寸的步骤直到第一和第三对象不充分相似为止。
文档编号G06F19/00GK101903885SQ200880121106
公开日2010年12月1日 申请日期2008年12月12日 优先权日2007年12月18日
发明者R·奥普弗, R·维姆克, S·卡布斯, T·比洛 申请人:皇家飞利浦电子股份有限公司