通过磁共振图像的自动分割进行CT替代

摘要

在对患者体内的解剖结构建模用于诊断或治疗规划时，可以访问经分割的共同配准的CT和MRI解剖结构参考图像的图谱集(26)，可以选择一个或多个这种结构的图像并使其与患者的临床图像中对应结构的MRI图像重叠。用户可以点击并拖动参考MRI图像上的界标或分割段边缘以使参考MRI图像变形，从而与患者的MRI图像对齐。患者MRI图像中的界标到参考MRI图像的配准也将患者MRI图像界标与共同配准的参考CT图像中的对应界标配准，并自动将来自参考CT图像界标的电子密度信息归属到对应的经配准的患者MRI界标。

说明书

本申请尤其应用于医疗成像系统中。然而，要认识到，所描述的技术也可以应用于其他类型的成像系统、治疗规划系统和/或其他医疗应用中。

对于辐射处理规划、组合的磁共振(MR)/正电子发射断层摄影(PET)及其他应用而言，MR成像模态由于其更好的软组织分辨能力而优于计算机断层摄影(CT)。与具有绝对尺度，即Hounsfield尺度的CT图像不同的是，MRI图像没有绝对尺度。这使得使用MR图像取代CT图像困难且不精确。

利用CT测量的Hounsfield值为在放射治疗规划(RTP)和PET衰减校正中例如基于图像的剂量计算提供了必需的信息。对于基于MRI的RTP和MR/PET系统而言，CT和这样测量的Hounsfield单位是不可用的。这就需要一种直接从MRI估计吸收系数或停止供能的机制。

现有技术中对便于克服上述缺陷的系统和方法的需求仍未得到满足。

根据一个方面，一种在磁共振图像(MRI)体积中提供电子密度信息的系统包括：存储器，所述存储器存储包括彼此共同配准的参考计算机断层摄影(CT)图像和参考MRI图像的图谱集(26)；显示器，所述显示器呈现患者MRI图像和参考MRI图像的视图；以及处理器，所述处理器执行用于将来自所述参考CT图像中的对应界标的电子密度信息与所述患者MRI图像中选定界标相关联的算法。

根据另一方面，一种为磁共振成像(MRI)图像中的体素提供电子密度信息的方法包括：从感兴趣区域的多个相应的CT和MRI图像产生类属共同配准的计算机断层摄影(CT)和MRI参考图像；接收患者MRI图像；以及配准所述患者MRI图像和所述参考MRI图像。该方法还包括：为经配准的患者MRI图像从所述共同配准的CT图像的对应分割段检索电子密度信息；以及将检索到的电子密度信息与经配准的患者MRI图像相关联。

根据另一方面，一种解剖结构的共同配准的计算机断层摄影(CT)和磁共振成像(MRI)参考图像的图谱集包括：从一个或多个对象的解剖结构的扫描图像产生的解剖结构的多个CT和MRI参考图像，其中参考图像被共同配准以将来自所述CT参考图像中的体素的电子密度信息与相应MRI参考图像中对应体素相关联。所述图谱集还包括机器可读介质，所述机器可读介质存储共同配准的所述多个CT和MRI参考图像，供操作员调用和操纵。

一个优点是向参考MRI图像体积传递了CT电子密度信息。

另一个优点在于向MRI图像上映射了空间值信息。

在阅读并理解了下述详细说明的情况下，本领域普通技术人员将认识到本创新的更多优点。

可以通过各种部件和部件设置，以及通过各种步骤和步骤设置的形式实现本创新。附图的目的在于对各方面进行图示，不应将其解读为对本发明构成限制。

图1图示了用于从由磁共振成像装置摄取的图像产生CT特有信息的系统，以辅助即将进行的PET/MR扫描的衰减校正、基于MR图像的辐射处理的剂量计算等。

图2图示了用于从由MRI装置产生的图像产生CT特有信息的方法。

图3图示了包括诸如CT、MRI等的多个成像装置的示范性医院系统，成像装置产生成像数据，由独立或共享的重建处理器重建成像数据以产生3D图像表示。

图1图示了一种系统10，其用于从由磁共振成像装置摄取的图像产生CT特有信息以辅助：对即将进行的PET/MR扫描的衰减校正、基于MR图像进行辐射处理的剂量计算等。系统10采用算法从MRI数据集导出CT特有信息，并且包括耦合到CT扫描器13和MRI装置14的用户界面12。MRI装置可以是常规MRI装置、集成的线性加速MRI(ILAMRI)装置等。在一个实施例中，使用组合的CT-MRI成像装置。用户界面包括处理器16，处理器16执行存储器18中存储的算法和/或例程以执行这里所述的各种功能。用户界面还包括用于向用户呈现信息的显示器20以及输入装置22(例如，鼠标、触针、键盘、麦克风或其他输入装置)。存储器18包括图像数据24(例如，由CT和/或MRI装置收集的)、一个或多个图谱集26、界标数据库28和参考结构数据库30。图谱集包括来自多个对象的通用共同配准的CT和MRI数据，其变形以收敛到解剖特性的公共通用映射。界标数据库包括解剖结构公共的多个界标，使用界标以允许操作员将来自数据库30的参考结构手动变形(例如点击并拖放等)到患者图像中的感兴趣区域或结构。

在一个实施例中，通过摄取具有不同形状和尺寸的不同患者或对象的若干CT和MRI图像并对每幅图像手动配准界标来填充图谱集26。在利用MRI对新患者成像以进行临床评估(例如，诊断、治疗规划、辐射剂量评价等)时，在新患者的MRI图像上重叠类属(generic)或通用MRI参考图像并在显示器20上加以显示。操作员对齐对应的界标，以可变形地将新患者的MRI图像配准到参考MRI图像上。因为先前MR和CT参考图像被彼此配准，所以新MR图像到参考MR图像的配准也将新MR图像配准到参考CT图像。访问参考CT图像以选择与CT图像中的对应界标、结构或区域相关联的Hounsfield值，并将它们归属于新患者MRI图像中的体素。通过这种方式，向MRI患者图像提供CT强度信息。

处理器16产生患者体内要被检查的区域的图谱集。为了产生图谱集26，由CT扫描器13和MRI装置14产生各患者体内的感兴趣区域的CT图像和MR图像。对这些图像都进行(例如手动地)分割以确定图谱集的CT和MRI部分中类似组织的对应区域。在辐射处理规划期间，产生当前患者体内感兴趣区域的MRI图像。处理器利用弹性变换将患者的MRI图像与图谱集MRI图像配准。患者MRI图像与MRI图谱集图像的配准也将患者的MRI图像与CT图谱集图像配准。基于这种配准，为患者MRI图像的每个体素选择对应的CT或Hounsfield值以在剂量和其他规划计算中使用。或者，可以将单一CT值分配给分割段的每个，并将患者MRI图像的分割段之内的所有体素与同一CT或Hounsfield值相关。

在辐射规划中，选择通过目标区域(例如，肿瘤等)的一系列轨迹。还为每个轨迹选择治疗处理的形状、强度和停留时间。执行选择以向目标区域递送期望的最小辐射剂量，同时不超过周围组织的最大可接受剂量。最大可接受剂量随着组织类型而变化。

在处理射束向目标前进时，其通过组织，组织会使射束衰减。衰减是Hounsfield或CY数的函数，并且在计算递送给目标的剂量时被计入考虑。如果在处理会话期间患者有运动，射束在通往目标的路径上穿过的组织可能会改变。进行运行中计算以调整处理方案，以确保向目标递送最小累积剂量而在其他组织中不超过最大值。

系统10不仅能够用于辐射处理规划中，而且能够在辐射处理期间使用以监测运动或患者与处理射束之间的失准。基于所感测的失准，可以修改处理方案，或者，如果患者移动太远，可以关闭处理射束。

在组合PET/MRI扫描中，知道针对被成像区域的每个体素的Hounsfield数也是有利的。基于Hounsfield数，可以进行PET辐射衰减校正。利用所描述的MRI/CT图谱集，可以从MRI图像确定每个体素的电子密度或Hounsfield数，以用在PET或SPECT成像技术的衰减校正中。

根据示例，用户使解剖学图谱集变形，以匹配个体患者的MRI数据集，其中解剖学图谱集包括3D CT数据集和一组解剖学参考结构30和解剖学界标28。变形过程提供了CT Hounsfield单位从图谱集到患者MRI数据集的映射。基于图谱集的变形能够自动根据MRI搜索确定最佳扫描平面。

在一个实施例中，从不同患者摄取的目标成像模态MR的一组手动分割图像数据集生成解剖结构的图谱集。该分割由感兴趣解剖结构的描绘以及突出的解剖界标构成。例如，借助于点分布模型在统计模型中对患者群内的描绘和界标的空间变化进行编码。

用户或技术人员通过检查图像的灰度值设置突出界标。亦即，周围的图像信息是对于给定界标的特征。因此，周围图像信息的模型伴随着界标的空间分布模型。在训练数据间界标以一对一的关系对应，借助于弹性图像配准，处理器16向公共灰度值图像弯曲所有训练图像。那么，这个图谱集图像包含每个界标附近的特征灰度值图像信息。

处理器16还针对各CT图像执行同一解剖结构的描绘。可以从不同患者获取训练图像(例如MRI、CT等)。由处理器从这个数据产生放射(radio)密度到空间区域/解剖结构的映射。由于CT有绝对Hounsfield尺度，多个数据集不是必要的；虽然如此，仍可以使用多个数据集检验不同患者间所描绘区域的放射密度值的一致性。

在另一实施例中，不是直接使用映射的Hounsfield值，而是处理器使用变形的图谱集来导出MRI数据的粗分割，以导出灰度值组织模型。

根据另一实施例，在组合的MR-PET医疗扫描器(未示出)中采用该系统。在处理器16从组合扫描器接收患者的MR数据时，检索出训练数据上定义的突出界标。例如，处理器可以采用标准图像处理技术并利用图谱集图像中的灰度值特征。为了进一步便于识别患者图像中感兴趣区域的空间位置和/或取向，可以限制基于空间位置的统计变化模型的搜索空间。在定位界标之后，可以从MR参考数据向患者的MR图像传送所描绘的结构。这是利用与用于产生灰度值Hounsfield参考图像相同的基于点对应的弹性方法来实现的。然后基于所存储的放射密度值到解剖区域的映射向患者MR图像中的区域分配放射(例如电子)密度值。给定带有放射密度值标签的MR图像，可以执行对应PET图像的适当衰减校正。

在涉及放射治疗处理规划系统的另一实施例中，在没有伴随的PET图像时采用类似方案。将所得的带放射密度标签的MR图像用于处理规划系统的剂量计算，并省却了用于执行期望任务时对CT图像的需求。

图2示出了用于从由MRI装置产生的图像来产生CT特有信息的方法。在50，从多个对象和/或其中的感兴趣区域的CT和MRI扫描产生参考图像的图谱集。从这一群体构造通用或类属参考CT图像体积和通用或类属参考MRI图像体积。CT参考图像具有与其中的体素、界标、结构和/或区域相关联的Hounsfield值。CT和MRI参考图像彼此共同配准，从而使得新的MRI图像到参考MRI图像的配准也将新的MRI图像自动配准到CT参考图像。

在52，产生患者MRI图像体积。在54，将患者MRI图像界标配准(例如自动或手动)到参考MRI图像中的对应界标。在56，电子密度信息从共同配准的参考CT图像界标和/或区域中被检索出来并被归属到对应的经配准的患者MRI图像界标和/或区域。然后可以将电子密度信息用于对后续的MR/PET或MR/SPECT患者图像执行衰减校正。例如，可以从共同配准的CT参考图像界标检索对于所选的经配准MR参考图像界标的Hounsfield单位值并将其归属于患者MRI图像界标。该方法然后返回到54，进一步重复患者MRI图像界标的选择和配准，并继续进行，直到已经配准预定或充足数量的患者MRI图像界标以执行期望任务(例如，治疗规划、辐射剂量计算、衰减校正等)。

参考图3，示范性医院系统可以包括多个成像装置100，例如CT 13、MRI 14等，成像装置100产生成像数据，由个体或共享的重建处理器102重建成像数据以产生3D图像表示。通过网络104向中央存储器106发送图像表示。

在与网络连接的站110，操作员使用用户界面12将选定的3D患者MRI图像表示从中央存储器移动到本地存储器114。例如，视频处理器116从参考图像存储器30在参考图像的图谱集26中选择用于在显示器20的第一观察窗118₁中显示的参考图像。可以由操作员手动、半自动地或自动地选择参考图像。在半自动选择时，为操作员呈现从参考图像的较大集合中自动选择的密切匹配感兴趣区域中器官(一个或多个)或组织的尺寸、个性特征(identity)和/或形状的参考图像的子集，操作员从该子集中进行选择。在自动选择时，处理器116自动选择符合预定标准(例如，患者图像中感兴趣区域中器官(一个或多个)或组织的尺寸、形状、个性特征等)的最佳匹配参考图像。

在第二观察窗118₂中显示患者MRI图像。第三观察窗118₃能够显示参考图像和患者图像的重叠，以允许用户将参考图像中的界标配准到患者图像中的对应结构或界标。例如，操作员通过界面12(例如，使用鼠标、触针或其他适当的用户输入装置)选择与参考图像中的界标对应的患者图像界标。视频处理器然后将来自参考图像的电子密度体素值归属到患者图像，以改善患者图像的分割。

在一个实施例中，不同类型的组织具有不同的体素值(例如，骨为“1”、脂肪组织为“2”等)，且每个体素值与电子密度相关联。患者图像到参考图像的配准将体素值与患者图像相关联，由此便于对患者图像中的结构进行密度估计。

在另一实施例中，操作员使用用户界面从图谱集26选择参考图像，并且视频处理器在观察窗118₃中将参考图像叠加到患者MRI图像上。为了使参考图像符合患者图像中一个或多个结构的形状，操作员使用输入装置选择患者图像的一个或多个分割段上的界标。操作员然后将参考图像上的第二界标指定为与患者图像上的界标对应。然后将与参考图像界标相关联的CT数据归属到患者MRI图像界标。

一旦已经配准了充足数量的界标，患者图像就可以被存储在中央存储器106中或直接用于另一过程中。例如，治疗规划(例如辐射、烧蚀等)站130能够使用轮廓来规划治疗会话。一旦规划得令操作员满意，就将规划的治疗传输到实施所规划的会话的治疗装置132。其他站可以在各种其他规划过程中使用成型的轮廓。

在另一实施例中，将具有替代CT数据的MRI患者图像用于辐射剂量的计算，例如，在配准期间使用与患者图像中的界标相关联的Hounsfield值。

在另一实施例中，可以调节观察窗118₃中显示的重叠图像以相对于CT图像数据对MRI图像数据加权，或反之。例如，可以调节滑动杆或旋钮(未示出)以改变MRI图像数据或CT图像数据的权重，其中该调节可以是机械的或出现于显示器20上并利用输入装置进行的。在一个示例中，操作员可以将观察窗118₃中的图像，从纯粹CT图像数据(如观察窗118₁中所示)，经过CT和MRI图像数据的多个和/或连续组合，调节成纯粹MRI图像数据(如观察窗118₂中所示)。例如，可以从0∶1到1∶0离散地或连续地调节CT图像数据与MRI图像数据的比例。

已经参考若干实施例描述了本创新。在阅读并理解了前述详细描述的同时，本领域技术人员可以想到修改和变化。旨在将本创新解读为包括所有此类落在权利要求及其等同替代的范围内的修改和变化。