用于分析影像的方法和仪器与流程

本发明涉及一种用于分析脑影像的仪器和方法，更具体地说，本发明涉及一种用于从血管可见的脑影像确定脑移位的计算机实现的技术。

背景技术：
这些年来，已经发展了各种神经外科技术。这些技术涉及将诸如导管或电极之类的神经外科器械插入到脑实质内的标靶。这些神经外科器械可以使用例如在EP1509153中描述类型的立体定向技术以亚毫米精度引导到一位置。尽管可获得能够将神经外科器械精确地引导到预定部位的立体定向仪器，但是无法保证这样的器械能够实际到达由外科医生通过术前MRI扫描等标识的脑内的标靶结构。例如，脑结构能够被定靶的精度可能会受到所谓的“脑移位”的影响，在用来标识脑标靶的MRI扫描之后经常发生这种“脑移位”。这种脑移位可能发生的原因多种多样。例如，可能由于患者头部在影像获取与手术台上的放置之间的取向的任何变化而产生脑移位。手术过程中脑脊髓流体（CSF）的泄漏和空气进入颅骨内也会导致脑恶化。而且，当神经外科器械被强制通过脑组织时，颅骨内的脑可能发生位移或膨胀。迄今为止，外科医生已经注意到了脑移位的存在，并且已经设计并修改外科手术程序以使其影响最小化。例如，已经尝试通过在外科手术过程中将患者头部始终以一定取向放置（不过，出于外科手术接近之目的，这种取向可能是次最优的）以使脑移位最小化。类似地，神经外科医生经常基于从之前的外科手术程序获得的关于可能发生的脑移位的大小的经验通过向MRI影像中标识的标靶位置添加偏移来补偿脑移位。尽管这种已知技术能够补偿大多数情况下的脑移位的影响，但是它们是高度主观的，并且完全依赖于神经外科医生的技巧和判断。脑移位的存在还意味着，通常不可能在外科手术程序本身过程中知道所植入的神经外科器械（例如电极、导管等）是否已经到达期望的标靶。这只能通过后续MRI扫描随后建立。之前已经针对图像引导的神经外科手术提出了许多技术，这些技术试图测量或预测脑移位。例如，WO2005/025404描述了一种皮层表面配准技术。具体地说，WO2005/025404描述了使用线条激光扫描器来术中测量皮层表面的三维轮廓。然后将该表面轮廓与之前获取的MRI数据进行配准。WO2010/037850描述了如何使用生物力学模型来估计脑移位，该生物力学模型分析术前获取的MRI影像和术中获取的多普勒超声影像。在2007年8月22至26日在法国里昂举行的医学和生物学协会的IEEE工程的2007年度国际会议的公报“Frameworkforalow-costintra-operativeimage-guidedneuronavigatorincludingbrainshiftcompenstion（用于包括脑移位补偿的低成本术中影像引导的神经导航仪的框架）”（MBucki等人，872至875页）（ISBN:987-1-4244-0787-3）中也描述了用于估计脑移位的生物力学模型。WO2006/028474描述了使用脑表面的术中影像生成脑位移的统计学图集。然而，这些技术中没有一种技术能够提供可用来评估和/或提高在复杂（例如深脑）神经外科手术程序中的装置放置精度的足够精确和可靠的脑移位测量。

技术实现要素：
根据本发明，提供了一种通过分析脑影像确定脑移位的计算机实现的方法，该方法包括如下步骤：（i）取得对象的脑的第一影像，该第一影像示出了所述脑内的血管相对于参考位置的位置；（ii）取得所述脑的第二影像，该第二影像是在对对象进行介入之后获取的，所述第二影像已经通过计算机断层扫描（CT）成像仪器获取并示出了所述脑内的血管相对于所述参考位置的位置；以及（iii）从所述脑的至少一个感兴趣区域中的血管相对于所述参考位置的的位置在所述第一影像和所述第二影像之间的移位来计算脑移位。因而，本发明包括一种计算机实现的方法，该方法用来从对象的脑的第一和第二影像建立脑移位。这些图像示出了在介入之前和之后对象的脑内的血管相对于一定参考位置的位置。所述脑的所述第二影像使用CT成像仪器（例如术中CA血管造影成像仪器）获取，因而可以被认为提供了脑解剖的空间精确表示。具体地说，CT成像仪器的使用提供了关于深脑区域的精确空间信息。所述影像可以示出整个脑或脑的选定部分。所述方法包括评估在第一和第二影像之间血管相对于固定参考位置（例如，颅骨或公共坐标系中的点）的运动，由此允许对脑移位的量进行量化。血管运动的量可以使用第一和第二影像中存在的全部血管或优选仅仅一些血管来建立。如下所述，步骤（iii）可以包括直接计算脑移位（例如，建立向量或标量脑移位值）或间接确定脑移位（例如，通过建立限定脑内位置的坐标系的修正或调节，或者通过修改这种坐标系内的一个或更多个标靶位置和/或轨迹）。本发明的方法具有的优点在于，该方法可以完全在手术室内进行（例如，使用在术中CT仪器获取的第一和第二CT血管造影影像），以便在神经外科手术过程中在一个或更多个情形下测量脑内的一个或更多个部位的脑移位。脑移位的这种测量可以用来在标定脑内结构时帮助外科医生（例如通过允许调节手动立体定向仪器的设定）或有助于更新外科手术机器人系统的标靶坐标。这样，能够提高标定精度，并因而提高处理效率。另选地，可以将术前获取的第一影像与术中（或术后）获取的第二影像比较。有利地，步骤（iii）包括使用共配准过程来确定所述至少一个感兴趣区域中的血管的位置的移位。换言之，执行共配准程序以发现将（所述至少一个感兴趣区域中的）血管从所述第一影像最佳地映射到所述第二影像或从所述第二影像最佳地映射到所述第一影像的变换。优选地，这种配准程序使用基于强度的配准技术执行。该配准程序可以包括非线性配准过程。这种非线性配置程序可以包括非可变形配准程序。该非可变形配准程序可以包括仿设配准过程。优选地，所述配准程序包括刚性配准程序。这种刚性配准程序可以包括使用刚性变换，其中一个影像被平移、旋转和/或缩放以发现所述至少一个感兴趣区域中存在的血管之间的最佳对应关系。在一个优选实施方式中，所述第二影像可以留下不改变，而所述第一影像可以进行平移、旋转和缩放中的至少一个以发现所述最佳对应关系。当所述第一和第二影像二者均使用CT成像仪器换取时，使用刚性配准程序是特别优选的，因而可以认为所述影像是空间精确的（例如基本没有扭曲）。由于公共参考位置在所述第一和所述第二影像中均有限定，因此所述变换提供了在感兴趣区域中的血管运动（并因而脑移位）的测量。步骤（iii）可以包括确定所述脑中的一个感兴趣区域中的脑移位。另选地，步骤（iii）可以包括为所述脑中的多个感兴趣区域建立脑移位。如果使用了多个感兴趣区域，这些区域可以重叠或相互不包括。可以在所述脑中的标靶部位（或多个标靶部位）附近的感兴趣区域（或多个感兴趣区域）中和/或在沿着穿过所述脑到达这种标靶部位（或多个标靶部位）的轨迹（或多个轨迹）的感兴趣区域（或多个感兴趣区域）中确定脑移位。感兴趣区域可以包括深脑标靶。优选的是，所述脑内的每个感兴趣区域的体积小于整个脑的体积。优选地，每个感兴趣区域的体积不大于整个脑体积的50%，更优选不大于整个脑体积的25%，更优选不大于整个脑体积的10%。在绝对方面，每个感兴趣区域优选小于200cm3，小于100cm3，更优选小于50cm3。使用脑的这种子区域作为感兴趣区域确保在不受到脑的血管运动可以相当小的其他区域（例如在颅骨附近）影响的情况下在该感兴趣区域中确定局部脑移位。当为深标靶（在深标靶处，脑移位的量可能与皮层表面附近出现的脑移位明显不同）建立脑移位时，这是特别有利的。所述第一和第二影像中的血管可以包括动脉和/或静脉。尽管脑中的大部分动脉与它们所穿过的脑组织一起运动，应注意的是，一些静脉可能表现出较少的运动，例如周围静脉可能某种程度地锚固到颅骨的内侧。因此，有利但绝非必要的是，所述第一和第二影像示出由动脉构成的血管。还优选但不是必要的是，任何感兴趣区域都不包括脑的包括倾向于不与脑组织一起运动的血管（诸如相对于颅骨某种程度上固定不动的周围静脉）的区域。方便地，使用在所述第一和第二影像中可见的参考结构来建立所述第一和第二影像二者共用的所述参考位置。在这种情况下，可以执行两个配准过程。可以执行第一配准过程，以建立在每个影像中可见的所述参考结构之间的第一变换。该第一配准过程可以基于所述第一影像和所述第二影像的全部的比较或者这些影像的其中所述参考结构可见的选定部分的比较。然后可以执行第二配准过程，以建立在每个影像的选定感兴趣区域中存在的血管之间的第二变换。于是，所述第一变换和所述第二变换之间的不同可以提供相对于所述参考结构的脑移位的测量（由血管运动证明）。可以使用在所述第一影像和所述第二影像中可见的任何适当参考结构。有利地，所述参考结构包括对象的在体内具有基本不变位置的解剖学特征。例如，所述解剖学特征可以包括对象的颅骨或颅骨的一部分。另选地，所述参考结构可以包括相对于对象的颅骨固定不动的基准标记组。例如，可以将一基准标记组附装（例如拧或粘结）到颅骨。除了提供在所述影像中可见的参考结构之外，脑的所述第一影像和所述第二影像可以通过同一成像仪器在对象的头部相对于所述成像仪器放置在同一固定位置的情况下获取。在这种实施方式中，所述参考位置可以包括所述成像仪器的位置参考点或原点。换言之，所述参考位置可以是由该成像仪器的成像系统限定的内部原点或一些已知参考点。如上所述，用来获取所述第二影像的成像仪器包括计算机化/计算机断层扫描（CT）成像仪器，例如x射线计算机断层扫描（CT）成像仪器。用来获取所述第一影像的成像仪器可以包括任何适当的医学成像仪器。例如，其可以包括断层扫描仪、计算机断层扫描仪（CT）、x射线计算机断层扫描仪、x射线、磁共振成像（MRI）、正电子断层扫描仪（PET）、超声波等。优选地，所述医学成像仪器提供三维或立体影像，或者这种立体影像能够从使用该仪器收集的数据构建。可以使用不同的成像仪器，包括不同类型的成像仪器来收集所述第一和第二影像。优选地，使用同一个成像仪器来收集所述第一影像和所述第二影像。在优选实施方式中，使用手术室内计算机化断层扫描血管造影（CTA）成像仪器来获得所述第一影像和所述第二影像。在这种示例中，所述第一影像和所述第二影像包括第一和第二CTA影像。这种CTA仪器可以方便地包括能够在手术室中使用的CT扫描仪装置，诸如由Medtronics公司出售的O-arm（注册商标）外科手术（CT）成像系统。使用手术室内医学成像仪器的优点在于，可以在患者经受外科手术的同时收集使用本发明的方法分析的所述第一影像和所述第二影像。因而能够快速地（例如在几秒钟或几分钟内）进行脑移位测量，并且用来帮助外科医生的结果使得随后的外科手术动作最优化。有利地，所述方法包括在所述第一影像中限定至少第一标靶部位的步骤。该步骤例如包括外科医生评阅所述第一影像并且标识需要放置一个或更多个神经外科器械（例如导管、电极）的部位。然后，可以直接在所述第一影像上限定由外科医生标识的部位或多个部位，例如可以在所述第一影像上的一定位置（电子地）标记一点。如果必要，可以限定多个这种标靶部位。有利地，可以通过从至少一个另外的影像输入标靶信息而在所述第一影像中限定所述第一标靶部位（和任意附加标靶部位）。所述另外的影像可以使用另外的成像仪器获取。该另外的成像仪器可以是与获得所述第一影像和所述第二影像不同类型的成像仪器。将标靶从所述另外的影像映射到所述第一影像上的步骤优选使用（例如以上描述类型的）配准技术执行。方便地，所述另外的影像包括使用MRI扫描仪获取的影像。有利地，所述另外的影像示出了脑内的血管。可以使用许多不同的MRI获取协议来获取其中血管可见的MR影像。例如，磁共振血管造影（MRA）可将流过血管的血与静止组织区分开。血管还可以在“T1”MRI扫描上看到，这种影像的对比可以可选地通过将对比剂注射到患者血流内（提供所谓的对比MRA影像）来提高。该配准程序因而可以基于在所述第一影像和所述另外影像中可见的血管的对比（例如配准）。优选地，该配准过程使用示出了脑内的至少一个感兴趣区域的数据来执行。这种感兴趣区域可以是围绕一定标靶部位的体积或感兴趣体积，并且可以与比较所述第一影像和所述第二影像时使用的感兴趣区域相同或不同。有利地，以上描述的共配准过程提供了可将位置在所述另外的图像和第一图像之间映射的必要位置关联或变换。所述另外的图像本身可以供外科医生使用以为将插入对象脑内的器械标识标靶部位或轨迹。另选地，可以从相对于所述另外的图像具有已知位置关系的一个或更多个其他图像的分析输入标靶信息。例如，将所述第一影像与另外的影像共配准的步骤有利地允许将也是使用所述另外的成像仪器（即，用来获取在配准过程中使用的所述另外的图像的仪器）的外科手术计划影像（或多个这种外科手术计划影像）与所述第一影像配准。所述外科手术计划影像可以包括由神经外科医生或通过使用适当的计划软件标识的标靶或多个标靶。这种外科手术计划影像例如可以包括术前获取的影像（例如使用在T2获取模式下操作以示出相关脑结构的MRI仪器获取），外科手术团队在几天或几周的时间内有时间更详细地分析这些术前影像。尽管可以使用任何合适的医学成像仪器，但是优选的是，所述另外的影像包括使用包括MRI扫描仪的另外的仪器取得的血管造影照片。有利地，所述外科手术计划影像使用（例如在T2模式下操作的）同一个MR扫描仪术前获取。当获取所述另外的影像和所述外科手术计划影像时，可以将对象保持在MRI扫描仪内的同一位置。这样，所述另外的影像中的位置将对应于所述外科手术计划影像中的位置。如果必要，在MR影像获取过程中可以使用MR可见的基准，以允许将所有获取的MRI影像关联至公共参考位置。这也可以用来说明扫描之间的对象的任何运动。由于MRI数据可能经受扭曲，因此还可以使得产生所述另外的（MRI）影像和所述第一影像之间的变换所用的配准过程是弹性的。然后可以变换（即拉伸/压缩而不仅仅平移/旋转）另外的（MRI）影像，以实现与所述第一影像的最佳拟合。这样，MRI扭曲对标定精度的影响可以降低或基本消除。另选地，可以使用刚性配准程序。使用本发明的方法分析的所述第二影像是在对患者进行介入之后获取的。所发生的介入可以包括可能导致脑在颅骨内移位的任何事件或与对象的相互作用。所述介入可以包括外科手术介入（例如切割对象的脑或将器械插入对象的脑内）或非外科手术介入（例如重定向对象）。应注意的是，外科手术介入不构成这里描述的本发明的第一方面的方法的一部分。有利地，所述第二影像（即如步骤（ii）中所限定的）已经在包括将神经外科器械插入所述第一标靶部位的介入之后获取。所述第二影像因而可以示出所插入的神经外科器械。介入还可以已经包括将另外的神经外科器械插入另外的标靶部位。步骤（iii）于是可以进一步包括评估所插入的神经外科器械相对于所述第一标靶部位的位置的步骤。除了测量脑移位之外，这还允许测量任何定位误差（例如涉及没有预料到的器械偏转等）。步骤（iii）可以包括直接或间接地建立脑移位。例如，步骤（iii）可以包括通过计算限定脑移位的大小的标量来直接确定脑移位。有利地，步骤（iii）包括通过计算既描述脑移位的大小又描述脑移位的方向（例如相对于所述参考位置）的脑移位向量来直接确定脑移位。在步骤（iii）中确定的脑移位可以有利地用于更新例如由外科手术机器人或立体定向系统使用的标靶坐标。步骤（iii）可以另选地包括间接建立脑移位。例如，步骤（iii）可以包括计算考虑到了脑移位的新的或更新的标靶坐标或建立相对于坐标系的偏移以说明脑移位。而且，这种脑移位的间接测量可以有利地用于提供已经考虑到了脑移位的作用的（例如外科手术机器人或立体定向系统的）标靶坐标。应注意的是，在外科手术过程中补偿脑移位的作用的能力是特别有利的。应注意的是，这里使用术语“影像”来指代从成像仪器（例如从立体医学成像仪或断层扫描仪）收集的数据，并不意味着这种数据必须以某种形式可视地显示。类似地，对影像中的“可见”结构的参考是为了说明，关于这种结构的信息在影像中存在，并且术语“可见”不应解读为要求这种影像或结构图形地显示。短语“取得影像”旨在表示搜取已经获取的影像（例如通过将影像加载到计算机内）和/或使用适当的成像仪器获取影像。以上描述的方法使用计算机来实现。因此，本发明还涵盖被编程以实现以上方法的计算机。当在适当的计算机上运行时实现以上方法的计算机软件也可想到。包括被构造成实现以上描述的方法的处理器的仪器也被本发明所涵盖。这里还描述了用于测量脑移位的仪器，该仪器包括：用于取得对象的脑的第一影像的装置，该第一影像示出了所述脑内的血管相对于参考位置的位置；取得所述脑的第二影像的装置，该第二影像是在对对象进行介入之后获取的，所述第二影像示出了脑内的血管相对于所述参考位置的位置；以及从所述脑的至少一个感兴趣区域中的血管相对于所述参考位置的位置在所述第一影像和所述第二影像之间的移位来计算脑移位的装置。所述仪器还包括以上在相关方法方面中描述的任一个特征。用于测量脑移位的仪器可以与外科手术机器人组合地使用，由所述外科手术机器人使用的标靶坐标和/或坐标系被更新以说明由用于确定脑移位的装置确定的脑移位。这里还描述了一种外科手术机器人系统，该外科手术机器人系统包括：机器人，该机器人用于将外科手术器械引导到位于对象的脑实质内的至少一个标靶部位；医学成像仪器（例如CT成像仪器），该医学成像仪器用来对对象的头部的至少一部分进行成像，该医学成像仪器适合于收集头部的示出了脑实质内的血管的一系列影像；和计算机控制器，该计算机控制器用来更新由所述机器人使用的所述至少一个标靶部位的坐标以说明脑移位，所述脑移位从通过所述医学成像仪器获取的多个影像确定。本发明的另一个方面提供了用于对患者执行机器人神经外科手术的仪器。该仪器优选包括机动化的神经外科机器人，该机器人用于将医疗器械移动到在机器人坐标系中限定的位置。该机器人可以是计算机控制的。该仪器还有利地提供了成像仪器，该成像仪器用来获取患者的脑的术中影像。例如，该仪器可以包括CT成像仪器（例如O-arm）。每个术中影像中的至少一个位置在所述机器人坐标系中优选是已知的。例如，优选在所述机器人坐标系和所述术中影像的坐标之间提供映射或关联。还提供了一种标定仪器，该标定仪器用来计算患者脑内的标靶位置，医学器械将前进到该标靶位置。所述标定仪器优选取得患者的至少一个之前（例如术前）获取的影像（在该影像中已经标识了标靶），并通过将所述至少一个之前获取的影像与使用所述成像仪器取得的至少一个术中影像配准来计算所述机器人坐标系中的相对应标靶位置。所述至少一个之前获取的影像可以已经术前获取（例如使用MRI仪器）或者在外科手术程序过程中的较早时间点获取（例如通过用来获取术中影像的成像仪器）。在一个优选实施方式中，所述标定仪器使用所述影像中可见的血管将所述至少一个之前获取的影像与所述至少一个术中影像配准。如这里详细描述的，可以将一定的感兴趣区域中的血管配准以提供所述机器人坐标系中的标靶部位。还描述了一种用于配准CT和MRI影像的方法，该方法包括如下步骤：（a）取得对象的脑的第一影像，该第一影像是使用CT仪器获取的，该第一影像示出了脑中的血管的位置；（b）取得对象的脑的第二影像，该第二影像是使用MRI仪器获取的，该第二影像示出了脑中的血管的位置；和（c）使用在所述第一影像和所述第二影像中可见的血管将所述第一影像和所述第二影像配准。在步骤（c）中可以使用以上描述的配准技术中的任一种。例如，可以使用刚性配准。优选地，将所述第一影像和所述第二影像配准的步骤（c）包括使所述第二影像弹性扭曲，从而使得在所述第二影像中可见的血管重叠于在所述第一影像中可见的相对应血管上。所述第二（MRI）影像的这种弹性扭曲将使用这种MRI仪器获取的影像中固有地存在的任何扭曲消除。方便地，将所述第一影像和所述第二影像配准的步骤（c）包括在所述第一影像和所述第二影像内的至少一个感兴趣区域中执行配准。优选地，每个感兴趣区域的体积不大于整个脑体积的50%，更优选不大于整个脑体积的25%，并且更优选不超过整个脑体积的10%。然后可以将在MRI坐标系中限定的标靶映射到CT坐标系中。这里还描述了一种外科手术仪器，该外科手术仪器包括：用于对患者进行外科手术的外科手术机器人，该外科手术机器人包括在机器人坐标系中具有已知位置的多关节型臂；和医学成像仪器，该医学成像仪器用于获取患者的术中影像，由所述医学成像仪器获取的影像的位置在影像坐标系中限定，其中所述外科手术仪器包括基准标记单元，该基准标记单元由所述多关节型臂保持，使其在所述机器人坐标系中具有已知位置，所述基准标记单元包括在由所述医学成像仪器获取的影像中可见的参考结构。这样，能够发现所述机器人坐标系和所述医学成像仪器之间的位置关系，而无需将基准标记直接附装（例如螺接）至患者的颅骨。这消除了在执行神经外科手术之前辅助医生将标记或框架附装到患者颅骨上的需要。优选地，所述基准标记单元以能够保持在外科手术机器人的多关节型臂的工具支座中的工具的形式提供。优选地，所述基准标记单元可以在可重复的位置（例如使用运动或虚拟运动附装机构）附装至所述多关节型臂。这里还描述了一种用于确定脑移位的方法，该方法包括如下步骤：（i）取得对象的脑的第一影像，该第一影像示出了所述脑内的血管相对于参考位置的位置；（ii）取得所述脑的第二影像，该第二影像是在对对象进行介入之后获取的，所述第二影像示出了脑内的血管相对于所述参考位置的位置；以及（iii）从所述脑的至少一个感兴趣区域中的血管相对于所述参考位置的位置在所述第一影像和所述第二影像之间的移位来计算脑移位。该方法还可以包括这里描述的步骤或特征中的任一个或多个。附图说明现在将参照附图仅以示例方式描述本发明，其中：图1示出了术前获取的MRI（T2）数据；图2示出了术前获取的MRA数据；图3示出了在神经外科手术之前获取的手术室内患者的CT血管造影照片；图4示出了在已经将导管外科地植入之后获取的手术室内患者的CT血管造影照片；图5示出了外科手术机器人和手术室内的CT扫描仪；图6是概述适当的配准过程的流程图；图7是概述用于配准术前CT数据和术后CT数据的过程的流程图；以及图8是概述用于配准MRA数据和CTA数据的过程的流程图。具体实施方式本发明涉及一种测量脑移位的方法，现在将在能够引入这种脑移位的外科手术程序的环境下描述本发明。应该注意的是，这里描述的用于分析影像的计算机实现的方法涉及分析在这种外科手术程序过程中取得的医学影像，并不旨在包括以下概述的任何一个外科手术步骤。下面描述的这些外科手术步骤纯粹用来帮助理解如何使用本发明的计算机实现的分析方法。参照图1，在非常一般的意义下图示了术前换取的（T2）MRI影像。该影像示出了患者的颅骨2和脑4。这样的T2MRI影像还将详细示出脑内的不同结构，不过为了容易理解将这些结构从图1中省略了。如在神经外科手术技术中的传统一样，这样的MRI影像由外科医生使用来计划外科手术程序并且通常是术前（例如，提前几天或几周）获取的。可能使用这样的MRI影像计划的外科手术程序的一个示例是放置导管以便将治疗剂输送到脑内的具体结构。图1还示出了标靶区域或标靶部位6和轨迹8，导管将沿着该轨迹8插入而到达位于标识的脑结构处的标靶。这些通常由神经外科医生手动地计划。图1所示的影像内的位置在MRI坐标系内限定。这样的MRI坐标系例如可以通过参考MR可见的基准标记来限定，这些基准标记在成像之前固定至患者头部（但是在图1中没有示出）。接下来参照图2，示出了在将MRI对比剂注射到患者血流中之后取得的（T1）MRI影像的图示。该影像可以被称为对比MR血管造影照片或“对比MRA”影像。该影像示出了脑内的各种血管10，但是没有示出脑内各种结构的任何细节。因此，该影像通常不适合于标标靶识和详细的外科手术计划目的。应该注意的是，除了对比MRA之外，还可以获取标准MRA（其中血管内的血流与这些血管周围的静止组织区分开）。图2中所示的对比MRA影像是与图1的MRI影像一样的患者的影像，并且所有位置均在同一个MRI坐标系内描述。这例如可以通过参考相对于颅骨具有固定位置的MR可见的基准标记来实现。另选地，患者头部可以在扫描之间保持在同一位置，由此消除基准标记的需要。在与MRI影像相同的坐标系内获取对比MRA影像允许将来自图1所示的MRI影像的标靶部位6和轨迹8映射到图2的其中血管可见的对比MRA影像。参照图3，其示出了患者的第一CT血管造影（CTA）影像。该CTA影像是在患者位于手术台上其头部固定至该手术台的期望位置的同时使用手术室内CT仪器（其优选是可动的）获得的，以用于随后的神经外科手术程序。在向患者血流内注射合适的CT对比剂之后获得的第一CTA影像示出了在CT坐标系中脑内的血管20的位置。在该示例中，该CT坐标系相对于在第一CTA影像中也可见的呈现颅骨30形式的参考结构来限定。然后执行配准程序以将图2的对比MRA影像（其包括标靶部位和轨迹信息）与图3的第一CTA影像配准。该配准过程涉及限定标靶部位6周围的感兴趣区域32，然后执行共配准过程以发现在该感兴趣区域中CTA影像和对比MRA影像之间的最佳位置匹配。该配准过程优选是基于强度的，并利用在CTA影像和对比MRA影像二者中均可见的血管。为了说明MRA扭曲，该配准过程可以是全部或部分弹性的，以允许将对比MRA影像变换（压缩/拉伸）而适合CTA影像的血管结构。下面更详细地概述优选的共配准程序的其他细节，但基本结果是数学变换，该数学变换允许将MRI坐标系中测量的位置映射到CT坐标系中。因此，该配准过程允许将在MR坐标系中限定的标靶部位6和轨迹8映射到CT坐标系中，如图3所示。因而，上述过程允许将基于术前获取的MRI数据的标靶和轨迹信息映射到在外科手术当天在手术室内建立的CT坐标系中。重要的是，这可以在不需要在术前MRI扫描和外科手术程序之间保持附着于患者的基准标记的情况下进行。如下面将更详细地描述的，CT坐标系还可以关联到机器人或立体定向坐标系，该机器人或立体定向坐标系允许诸如立体定向框架或外科手术机器人之类的外部设备将外科手术器械插入到已经在基于CT的坐标中描述的脑内的期望标靶。参照图4，示意性地示出了第二CTA影像，该第二CTA影像是在导管40插入之后取得的。该第二CTA影像是利用与获取图3的第一CTA影像使用的相同手术室内CT仪器取得的，并且患者的头部仍然固定至手术台（注意，在外科手术过程中可能发生一些头部运动）。如从图4可以看到，导管已经沿着图3中所示的轨迹精确地插入，并且使其顶端精确地位于使用图3中所示的第一CTA影像限定的标靶坐标处。然而，在外科手术程序过程中，在一些点处脑已经相对于颅骨移位（注意，这种脑移位在图4中极其夸大地示出）。该脑移位（注意，该脑移位导致血管采取图4中的如20’所示的位置）意味着，导管的顶端并没有实际放置在脑内的标靶部位（治疗剂的输送需要该标靶部位）。换言之，有由图4中的向量B限定的脑移位，该脑移位导致导管的顶端从其在脑实质内的期望标靶偏离。应注意的是，在图4中，假定导管精确位于其原来计划的位置，并且所有误差均源自于脑移位。当然，导管放置时可能还有其他误差，或者导管可能由于脑组织的密度变化而从期望轨迹偏离开。在这种情况下，也可以从第二CTA影像建立实际导管位置从计划位置的任何变化。通常，上述类型的脑移位只能在术后（例如在随后的MRI扫描过程中）或者当治疗失败时标识出。然而，第二CTA影像的获取允许以如下方式在手术前后测量任何脑移位。首先，执行配准程序，以建立描述第一和第二CTA影像中所示的颅骨（或者颅骨的一个或更多个预定区域）之间的空间关系的参考变换。然后在标靶部位周围的脑的感兴趣区域32中将第一CTA影像与第二CTA影像配准，这提供了脑（血管）变换。参考变换和脑（血管）变换之间的差描述了脑内的感兴趣区域32在第一和第二CTA影像之间相对于颅骨已经移动了多少。换言之，使用颅骨作为固定位置参考特征可以建立描述感兴趣区域32中的脑移位的向量。这种测量脑移位的方式可以完全地用作将伴随一定神经外科手术程序的脑移位的量量化。如果脑移位超过一定量，则手术可能撤销（例如通过移出导管）并且再次进行（例如通过将导管植入到移位的标靶部位）。另选地，在外科手术程序过程中可以多次测量脑移位，并且如果测量到了任何明显的脑移位，可以对标靶坐标进行适当修正。例如，可以在颅骨中钻出开口并刺穿脑硬膜之后测量脑移位。所测量的脑移位然后可以用来更新将在随后的立体定向或机器人手术过程中用以插入神经外科器械的标靶坐标。使用上述方法测量脑移位还可以在神经外科器械植入过程中周期性地测量，并且相应地更新标靶坐标或外科手术坐标系的原点。现在将参照图5来描述一种技术，该技术用来将神经外科立体定向机器人的坐标系（例如由RenishawMayfield生产的Neuro|mate机器人系统50）关联到相关的CT仪器52的坐标系。具体地说，图5示出了Neuro|mate机器人系统5，该系统50安排用来对患者54执行外科手术，患者的头部在关联至机器人50的手术台56上刚性地保持在适当位置。一个移动式手术室内CT扫描仪52位于患者头部周围，并被布置成在包围头部的CT成像空间58中获取CT影像。机器人系统50还包括多关节型机械臂60，该机械臂60在其远端处具有支座62，神经外科器械等可以附装至该支座62。在图5中，该多关节型机械臂60携带CT可见标记设备64，该标记设备64相对于支座62具有已知（预校准）的位置。也就是说，该CT可见的标记设备64在该机器人坐标系内具有已知位置。为了建立机器人坐标系和CT坐标系之间的位置关联，多关节型机械臂60将CT可见的标记设备64放置在CT成像空间58中，并且移动式CT扫描仪52执行CT扫描。在这种扫描过程中用于该CT可见的标记设备64的合适位置在图5中示出为64’。因而，所得到的CT影像（其可以是以上描述的CT血管造影影像）示出了患者头部内相对于CT可见的标记设备64的位置（例如血管和标靶部位的位置）。因而，对CT影像的分析允许在CT坐标系和机器人坐标系之间建立关联或变换，由此允许将已经在CT坐标系中限定的标靶位置和轨迹（如上所述，其可能已经从MRI坐标系映射到CT坐标系）映射到机器人坐标系内。然后可以该标记设备64从支座上拆下并替换为外科器械，可以沿着期望轨迹将该外科器械引导到标靶坐标。该标记设备64可以在需要时（在同一位置）再次附装至支座，并再次成像。例如，在CT坐标系和机器人坐标系之间的映射需要重新建立或校验时。应注意的是，也可以在患者出现之前建立CT坐标系和机器人坐标系之间的映射。例如，可以在患者被安置在手术台56上之前通过机器人的机械臂将CT可见的标记设备保持在CT扫描仪的成像空间内的各种不同位置。然而，这是假定当患者进入时CT扫描仪和外科手术机器人之间的相对位置不变时才这样。另选地，患者54可以将该CT可见的标记设备附着于它们的颅骨，例如作为头部框架的一部分或通过直接（例如螺纹）附装。这种附装件可以包括一特征以与由机器人50的机械臂60保持的器械接合以便能够在机器人坐标系内建立CT可见的标记的位置。一旦已经建立CT坐标系和机器人坐标系之间的变换，如果在机器人坐标系内限定的标靶位置改变（例如由于脑移位），则可以更新这些标靶位置。例如，可以以上述方式（即通过查找相对于颅骨的这种脑移位）测量脑移位的量。假定患者的颅骨相对于机器人不移动，则可以将所测量的脑移位直接用于修正机器人坐标系内的标靶坐标。应该注意的是，因为在上述方法中相对于颅骨来测量脑移位，所以CT扫描仪在测量之间可相对于患者移动。而且，这种用于脑移位的修正可以在手术中周期地或在一定时间点来确定。然后可以在手术过程中相应地改变机器人坐标系内的标靶位置以说明任何脑移位。还应该注意的是，机器人坐标系的原点可以调节以说明脑移位。这与改变标靶位置具有相同作用。尽管以上示例提及了评估标靶区域中的脑移位，但应注意的是，可以对脑内的多个感兴趣区域测量脑移位。这些感兴趣区域可以对应于多个标靶部位（例如，如果植入多个神经外科器械的话）和/或沿着插入轨迹的多个点（例如为了确保插入的器械不损坏敏感组织或刺穿任何大血管）。重要的是注意到，尽管上述示例示出了二维影像，但这些二维影像仅仅是三维或立体影像的切片。提供上述二维图示仅仅是为了容易理解，本领域技术人员将充分认识到，在大多数实际情况下，坐标系和任何脑移位都是以三维限定和测量的。以上在一般方面描述了影像的配准，但是现在将参照图6至8描述合适的配准过程的更多细节。可以使用各种已知配准过程来执行表明血管的两个数据组（例如，Angio-MR、Angio-CT、血管模型或血管图集数据组）的共配准。例如，合适的配准方法包括如下三组技术中的任一组。首先，可以使用基于强度的配准方法。如果3D血管数据包括血管模型或图集，优选在配准之前将其转换成3D体积数据。在这种技术中，通过使强度相似性最大化来确定变换。其次，可以使用基于特征的配准方法。这些方法通常从两个数据组提取一定特征。合适的特征可以包括血管连接部和/或血管形状。可以通过使不同组的标识特征之间的距离最小化来估计最佳变换。最后，可以使用基于强度和特征配准方法二者的混合技术。这种混合方法包括所谓的主动外观模型和颗粒系统方法。现在参照图6，示出了基于强度的配准过程的流程图。如从图6可以看出，步骤70包括取得所谓的移动数据（诸如以上提到的第二CTA影像），步骤72包括取得固定数据（诸如以上提到的第一CTA影像）。针对固定和移动数据限定感兴趣区域，例如在脑中限定的标靶部位周围的小区域。执行步骤74，该步骤74利用变换（T移动→固定）将（感兴趣区域中的）移动数据扭曲。在步骤76中通过成本函数的计算比较扭曲的移动数据和固定数据的强度。步骤78然后试图选择用于在变换（T移动→固定）中使用的、反馈给步骤74中的值的改进参数。该过程反复迭加，直到无法进一步提高扭曲的移动数据和固定数据之间的相关性。然后在步骤80中将变换的最终的最佳值返回。该配准过程也可以看作是两个基本部分的叠加。第一个部分是变换模型的选择，第二个部分是变换参数的估计。而且，大多数影像配准问题都可以看作是最优化问题，该最优化问题在一般数学方面可描述如下：其中I移动和I固定是（待扭曲的）移动和固定（标靶）数据，T移动→固定是所选的变换模型中的参数向量，W(I移动|T移动→固定)是具有参数T移动→固定的变换模型下的扭曲移动数据，而f成本(·)是用来测量扭曲移动数据和固定数据之间的强度相似性的成本函数。接下来参照图7，示出了描述使用配准技术测量术前和术后CTA影像之间的脑移位的流程图。该配准技术使用基于两个区域的刚性配准（三个平移和三个旋转角），该刚性配准使用术前和术后CTA数据。因而，该方法通过测量在待测量脑移位的点（例如在探针的顶端处）处的相对于颅骨的空间改变来估计脑移位。这给出了一个向量，该向量标识脑移位相对于颅骨的大小和方向。该配准基于体积CTA数据，并且其使得切片距离和切片厚度问题最小化。由于其是基于CT的，其也避免了任何MRI扭曲。更详细地参照图7，步骤90和92分别包括取得术前影像和术后影像。两个影像示出了参考结构（例如颅骨）和脑内的血管。执行步骤94，该步骤94包括基于参考结构（例如颅骨）的术前和术后CTA影像之间的刚性配准，这产生变换还执行步骤96，该步骤96包括将术前和术后CTA影像中的预定感兴趣区域中的血管共配准。这产生了变换该感兴趣区域的尺寸可以例如是50×40×15mm。然后，按照如下计算在两个探针（左右探针）的顶端在插入之后所位于的两个点处的脑移位。首先，在CT机器的UVW坐标中的术后CTA影像中确定每个探针的顶端位置。这些位置被限定为（a）r固定和l固定。然后使用变换在术前CTA中计算探针顶端的位置，由此给出位置（a）和（b）然后使用在术前CTA中计算探针顶端的位置，这提供了（a）和（b）也使得可以从如下表达式计算每个探针的位置处的脑移位（Sl和Sr）：参照图8，图示了用于将MRA数据与CTA数据配准的过程。步骤100和102分别包括取得（例如找回或获取）CTA数据和MRA数据。然后执行步骤104，该步骤104在血管可见的区域中限定感兴趣区域。该感兴趣区域可以位于将脑内的要插入神经外科器械的标靶结构或部位周围。然后执行步骤106，以将通过确定描述两个（MRA和CTA）数据组之间的对应关系的最佳变换来共配准这两个数据组。然后可以将来自MRA数据（例如关于标靶结构的位置）映射到CTA数据上。以上描述的技术优选使用适当变成的计算机来实现。应该记住，上述示例是非限制性的，并且仅仅是为了帮助理解本发明而提供的。本领域技术人员将完全意识到根据以上概述的一般教导而可能实现的上述示例的许多变型。