电解剖标测点到对应的图像数据的配准的制作方法

专利名称：电解剖标测点到对应的图像数据的配准的制作方法
技术领域：
本发明大体关于医学成像领域。更具体地，本发明涉及将电解剖标测 点配准到预采集的成像数据集。
背景技术：
用于心血管应用的图像引导疗法包括例如从MRI、 CT、超声或者荧光 透视成像获取的诸如3D成像数据的预采集的体积成像数据与过程中电解 剖标测(EAM)信息的整合，所述过程中电解剖标测(EAM)信息可以由 磁场、电场、或超声技术进行定位。该技术依赖于两个数据集的正确对准， 也就是所公知的配准过程。
现有用于对EAM和成像数据进行配准的方法基于点到表面距离最小 化算法(例如，迭代最近点算法ICP)。这种方法要求在标测过程中对心脏 内表面的点进行特异性选择，然后对EA—M点与从成像数据集分割的腔室表 面边界进行点到表面的配准。
常规用于点采集和后续点到表面配准的方案依赖于技术且耗时费力， 这是由于必须首先操作标测导管尖端以使其到达心内膜上的路标，之后， 必须在EAM数据记录中对尖端位置进行明显标注。这一单独的路标识别和 标注过程必须重复50到100次，以便为基于ICP的配准定义出具有足够细 节的心内膜表面从而使函数准确。为此，在患者研究的临床或诊断或治疗 分量开始之前，仅仅为了配准过程就必须进行15-45分钟的标测。应该认识 到，未与腔室表面接触的任何所采集的EAM点位置反过来影响用点到表面 距离最小化的配准的质量。
因此，改进的点采集方法能够有利地提高灵活性，并能够节约成本和 时间
发明内容
因此，本发明优选地通过提供根据所附权利要求的方法、装置、系统 和计算机可读介质，单独地或以任意组合的方式寻求缓解、减轻或消除一 个或多个上述现有技术中的不足和缺点，并解决至少上面所提到的问题。
在本发明的一方面，提供了一种对感兴趣体积中的结构进行识别的方 法。所述方法包括以连续模式采集与所述结构相关的多个点，随后，将至 少一个所述点配准到所述结构的先前所采集的成像数据集。
在本发明的另一方面，提供了一种对感兴趣体积中的结构进行识别的 装置。所述装置包括以连续模式采集与所述结构相关的多个点的采集硬件， 以及随后将数据点中的每一个配准到所述结构的成像数据集的配准硬件。
在本发明的另一方面，提供了一种包括根据本发明的上述方面的装置 的医学成像系统。
在本发明的又一方面，提供了一种其上嵌有计算机程序的计算机可读 介质，所述计算机程序由计算机进行处理，以对感兴趣体积中的结构进行 识别。所述计算机程序包括以连续模式采集与所述结构相关的多个点的采 集代码段，以及将至少一个点配准到结构的先前所采集的成像数据集的配 准代码段。
在本发明的再一方面，提供了一种对感兴趣体积中的解剖结构进行识 别的方法。所述方法包括以连续模式采集与所述解剖结构相关的多个点， 随后，将至少一个点配准到所述结构的过程前或过程后所采集的成像数据 集。
本发明通过对标测系统进行修改，提供了更快的EAM点的采集，使得 对导管尖端位置进行自动且连续的记录，而不需要显式导航(explicit navigation)到心内膜，也不需要在心内膜上标注基准路标。


本发明能够具有的这些和其他方面、特征和优势将从下述参照附图对 本发明的实施例的描述中变得明显并参照其进行阐述，在附图中-图1为示出EAM点到预采集的成像数据的常规配准的图示； 图2示出了根据实施例的连续点云(簇)采集的方法； 图3为示出根据实施例的方法的图示；图4为示出根据实施例的装置的图示；
图5为示出根据实施例的装置的图示；以及
图6为示出根据实施例的计算机可读介质的图示。
具体实施例
图1示出了 EAM点到预采集的成像数据的常规配准，其需要繁琐的导 管导航和靶标点选择，以识别心内膜上的表面点。在来自术前成像(示出 为带阴影的表面模型)的经分割的心脏表面模型上对这些EAM点进行配准 和覆盖。
下列描述集中于本发明的应用于医学成像的实施例，具体而言为将连 续采集的电解剖标测点配准到成像数据集。然而，应该理解，本发明并不 局限于这一应用，而是可以应用于需要将分散的3D点云(多个点)与多态 性体积成像数据集进行快速融合的用于例如心脏恢复室、EP实验室或图像 引导外科手术室的介入性引导或诊断的任何术中过程。
本发明提供了一种点采集方法，其中，没有对EAM数据集中的心血管 边界进行明确定义，所以不能使用传统的点到表面距离最小化(例如，ICP) 方法。即远离心脏内表面的内部点位置会反过来影响配准质量。
根据图2a，在本发明的实施例中，提供了一种连续点云(非组织和非 靶标的点簇)采集的方法，该方法不同于图像引导标测的常规点采集。所 述方法不需要(一个或多个)导管尖端到心内膜的特定导航。所述方法提 1共了以连续方式以非门控(im-gated)或门控(其中，可以将诸如ECG信 号或者呼吸信号的时间选通信号用作触发以标记在配准过程中应该包括哪 个EAM点)模式进行记录的EAM位置(点)21的自动采集，而电生理学 家/用户用操纵具有一个或多个位置传感器元件的导管贯穿诸如心室的感 兴趣的心血管结构。
通过将所述方法用于连续点云采集，在非常短的时间段中采集大量的 点。例如，如果以非门控模式采集所述点，在5分钟内可以有效地采集大 约10，000个点。由于生成的体积点云包括分散轨迹宽泛的导管尖端位置， 并对所述体积点云进行密集采样，这能够得到关于感兴趣的心血管结构的 几何形状的良好定义，从而使其足够与术前的成像数据进行配准。
7在本发明的实施例中，提供了将连续采集的EAM点云映射到至少一个 成像数据集的方法。所述方法可以概括为两个分开的步骤，其中，第一步 骤包括至少一个成像数据集的图像处理。第二步骤包括将连续采集的点云 对准到所处理的至少一个成像数据集，或者执行连续采集的点云到所处理 的至少一个成像数据集的配准。
图2b中示出了使用连续点云采集的方法的结果，所述方法将连续采集 的点映射到成像数据集，其中，将所述连续采集的点对准到过程前或过程 后所采集的成像数据集。
在一个实施例中，提供了一种识别感兴趣体积中的结构的方法。所述 方法包括以连续模式采集与结构相关的多个点，并随后将所述至少一个点 配准到结构的先前采集的成像数据集。
根据图3，在一个实施例中，配准包括通过分割对成像数据集进行处理 33，以在成像数据集中对感兴趣体积中的结构进行定位，从而生成构成所 述结构的经分割的表面。另外，配准包括基于经分割的表面，通过将位于 结构表面内的体素定义为非零掩模值(例如值1)，同时将结构表面外的体 素赋予零掩模值，来对3D二元掩模(binary mask)进行计算34。
另外，配准包括计算35 3D距离函数或者反映从3D 二元掩模到感兴趣 的结构表面的距离的其他计算度量。步骤33-35的总体目标为从成像体积提 取3D距离测量，可以将所述测量用于下述步骤，以确定将多个点配准到成 像数据集的感兴趣体积中的结构的接近程度。
配准还包括通过使用从例如过程前MR/CT和电解剖标测系统(即 ARTO)的已知坐标系定向导出的表变换而启动36对连续采集的多个点与 ik像数据集的对准。另外，配准包括对未知的配准参数进行迭代估计37， 以相对于成像数据集对多个点进行更新。例如，给定一组估计的参数值， 可以通过将变换矩阵乘以3D点坐标而执行与成像数据集相关的多个点的 空间变换。
配准还包括将每个点与来自距离函数的对应值进行比较38，从而根据 所估计的配准参数对目标函数值进行计算39。在一个实施例中，对于所有 的点，将对在步骤38中于每个点处计算的距离函数值进行求和，从而生成 对应于当前配准变换估计的目标函数值。另外，配准包括例如通过基于共轭梯度的函数最小化，对所估计的配准参数30进行迭代更新，直到所估计 的配准参数的变化小于预定的容差，而将目标函数最小化。 在一个实施例中，多个点为电解剖标测点。
在一个实施例中，至少一个成像数据集为过程前采集的成像数据集。 在一个实施例中，至少一个成像数据集为在连续采集多个点之后进行 的过程后采集。
在一个实施例中，至少一个成像数据集为在连续的多个点采集之前采 集的。
在一个实施例中，成像数据集的感兴趣体积中的结构为心脏结构。
在一个实施例中，至少一个成像数据集的图像处理包括分割和进行掩 模，即，通过使用至少一个分割工具计算3D 二元掩模，所述分割工具可以 为通常可用于成像工作站(例如Siemens Syngo或Philips Brilliance工作站) 上或者图像分析软件(例如 Analyze 、 MATLAB 、 Insight Registration/Segmentation Toolkit)内的诸如用于可变形表面的算法、水平集 分割方法、连通分量区域生长等的各种分割工具中的任一种。处理还包括 基于掩模计算3D距离函数，使得距离值对于经分割的腔室内的所有像素都 为零，而对于腔室外的每个像素，距离值为从成像数据集到经分割的表面 上的最近边界点的距离。
在一个实施例中，分割包括向成像数据集使用在对形状约束表面模型 进行基于标准化区域的演化的基础上的可变形表面。
在一个实施例中，配准步骤包括通过将点云的质心移位以匹配3D 二元 掩模的质心，而将连续采集的点云与经处理的成像数据集(例如，来自3D 旋转血管造影(3DRA)的预采集CT/MR体积数据或术中采集的类似CT 的重建)参照系进行粗略对准。接下来，例如利用将EAM坐标系与MR或 CT成像系统的坐标系进行映射的近似表变换开始，将未知的刚体变换参数 的初始估计用于相对于成像体积更新EAM点云位置。之后，使用EAM点 云内的每个点的当前位置查找距离函数中的对应值，对云内的所有点的距 离值进行求和以计算目标函数值，所述目标函数值依赖于未知的配准参数， 在3D中，所述目标函数值对应于X/Y/Z中的三个未知平移参数并且对应于 围绕X/Y/Z轴的三个未知旋转角。接下来，通过利用例如基于共轭梯度的函数最小化方法，通过对EAM点云和所处理的成像数据集之间的未知的基 于刚性的变换的参数估计进行迭代地更新，而对目标函数进行最小化。预 定阈值(例如运动O.lmm)定义了何时迭代被认为已经收敛。最终的配准 参数表示EAM坐标系和成像系统坐标系之间经估计的变换。
在一个实施例中，至少一个成像数据集为术中采集的成像数据集。 , 在一个实施例中，处理步骤和计算步骤不包括分割和3D 二元掩模生 成。在这一情况下，如在图3中所示，方法从步骤31开始，而非从步骤32 开始。这可以是如下这样的情况对其而言心血管腔室具有相对一致强度 值的成像数据集的情况，例如，对于对比增强的X射线/CT成像数据集而 言或者对于暗血/亮血MRI脉冲序列而言，心血管腔室具有相对一致强度值 的成像数据集。如果不需要进行分割，降低了处理步骤的总时间，在一些 实施例中，这有很大的优势。然而，在一些情况下，当不执行分割时，降 低了多个点的对准准确性。因此，在一些实施例中优选进行分割，而在一 些实施例中，优选不进行分割。不进行分割的实施例是有优势的，这是由 于不需要分立的分割步骤，并且其使得临床工作流程更加快速。在过程前 成像数据的体积可视化足够用于电解剖标测的过程中，可以采用不进行分 割的方法。
在使用经分割的表面模型的其他实施例中，分割方法为优选的实施例。 在一个实施例中，配准步骤包括基于已知的表变换将连续采集的点云进行 粗略对准，将EAM坐标系映射到成像坐标系。接下来，使用点位置在体积 成像数据集中的对应位置处对图像强度值进行采样。所有点位置处的强度 值的分布提供了可以用于执行变换参数估计的强度直方图。在一致增强的 心血管腔室的情况下，强度直方图会显示出围绕平均强度值的相对紧密分 布。对于不良配准估计，强度直方图会显示出源自从位于感兴趣的心血管 结构之外的体素采样的点的更大分布。
根据一个实施例，使用诸如从直方图得出的标准偏差或熵的测量，以 能够对变换参数估计进行迭代，从而将对3D点云中的所有点采样的体素强 度值的相似性最大化。因此，可以在不需要进行分割的情况下得出最后配 准变换估计。
10根据图4，在一个实施例中，提供了识别感兴趣体积中的结构的设备 40。所述设备包括以连续的模式采集与结构相关的多个点的采集硬件41， 和随后将每个数据点配准到结构的成像数据集的配准硬件42。
根据图5，在本发明的一个实施例中，配准硬件还包括通过分割对成像 数据集进行处理的处理硬件51。装置还包括计算3D 二元掩模的第一计算 硬件52，和计算距离函数或者反映到感兴趣的结构表面的距离的其他计算 度量的第二计算硬件53。另外，装置包括对准硬件54,其用于将连续采集 的多个点与成像数据集进行对准；以及估计硬件55，其用于估计配准参数， 以相对于成像数据集对多个点进行更新。另外，装置包括比较硬件56，其 用于将每个点与来自距离函数的对应值相比较；以及第三计算硬件57，其
用于根据所估计的配准参数计算目标函数的值；以及最小化硬件58，其通 过利用基于共轭梯度的函数最小化对所估计的配准参数进行迭代更新，直 到所估计的配准参数的变化小于预定的容差，而将目标函数最小化。
在本发明的一个实施例中，配准硬件能够执行根据一些实施例的方法。 在一些实施例中，装置的采集硬件包括处理器和存储器。处理器能够 启动对多个点的采集，并将多个点存储到存储器中。可以将不同的处理器 和存储器用于采集过程。
在一些实施例中，装置的配准硬件包括处理器和存储器。处理器能够 执行配准过程，而存储器用于存储。可以使用用于配准过程的不同步骤的 不同处理器和存储器。可以将存储器集成在处理器中，诸如处理器内部存 储器。
处理器可以为诸如Intel或AMD处理器、CPU、微处理器、可编程智 能计算机(PIC)微控制器、数字信号处理器(DSP)等的各种存储器的任 一种。然而，本发明的范围不局限于这些特定处理器。
存储器可以为任一种能够存储点信息的存储器，例如双密度RAM (DDR、 DDR2)、单密度RAM (SDRAM)、静态RAM (SRAM)、动态 RAM (DRAM)、图像RAM (VRAM)等的随机存取存储器(RAM)。存 储器也可以为诸如USB、 Compact Flash、 SmartMedia、 MMC存储器、 MemoryStick、SD Card、MiniSD、MicroSD、xD Card、TransFlash和MicroDrive
ii存储器等的FLASH存储器。然而，本发明的范围不局限于这些特定的存储 器。
存储器和任选的处理器内部存储器可以包括执行多个点的采集和配准 的软件。
在一个实施例中，采集硬件还包括通过使用导管上的适当感测元件能 够进行感兴趣体积的连续点采集的导管。所述导管可以具有多个分离的感 测元件以进一步促进点采集过程。
若干种当前可获得的进行点采集的导管可以用于通过以连续的方式进 行问询(poll)而执行连续的点采集，并且将测量数据连续地传送回控制器 软件。
在一个实施例中，导管利用电磁进行点位置的连续电磁感测。 在一个实施例中，导管利用超声迸行点位置的连续超声感测。 在一个实施例中，装置为医学工作站。 在一个实施例中，提供了包括所述装置的系统。
在一个实施例中，所述系统为诸如CT系统、MRI系统或超声系统的 医学成像系统。
根据图6，在一个实施例中，提供了其上嵌有由计算机处理的计算机程 序60的计算机可读介质。计算机程序包括执行连续EAM点云采集的代码 段61。计算机程序60还包括将所采集的连续点云对准到至少一个成像数据 集的配准代码段62。
计算机可以为诸如处理器或存储器(参见上述描述)的具有处理和存 储能力的任何设备。
在一个实施例中，计算机程序具有于执行根据本发明的实施例的方法 的代码段。
在另一实施例中，计算机程序包括控制例如由导管执行的连续点采集 的代码段(控制器软件)。
在本发明的一个实施例中，计算机程序包括在根据一些实施例的装置 的硬件中。
在一个实施例中，感兴趣体积为能够用导管进行点云的连续采样的任 何解剖学空腔。在一个实施例中，感兴趣体积为能够用导管进行点云的连续采样的任 何非解剖学空腔。
在本发明的具体实施方式
中，首先，将导管引入诸如患者心脏的器官。
使用导管对描述器官空腔的数据点进行连续采集。之后，对器官的CT或 MRI成像数据集进行处理，并将所采集的器官空腔的点对准到经处理的成 像数据集。对准之后，可以开始心脏的图像引导疗法。
本发明可以以包括硬件、软件、固件或者这些的任意的任何适合的形 式实现。然而，优选地，将本发明实现为在一个或多个数据处理器和/或数 字信号处理器上运行的计算机软件。本发明的实施例的元件和组件可以以 任何适合的方式进行物理地、功能性地或者逻辑性地实现。实际上，可以 在单个单元上、多个单元上或者作为其他功能单元的一部分实现所述功能
性。这样，本发明可以在单个单元上实现，或者可以物理上地或功能上地 分散到不同单元和处理器。
尽管上面参考具体实施例对本发明进行了描述，但本发明并不局限于 上述实施例中的具体形式。实际上，本发明仅由所附权利要求限定，上述 特定实施例之外的实施例同样落入所附权利要求的范围。
在权利要求中，词语"包括"不排除其他的元件或步骤。另外，虽然 是单独列出的，但是多个硬件、元件或方法步骤也可以例如由单个单元或 咎理器来实现。另夕卜，尽管在不同的权利要求中包括单独的特征，可以有 利地对这些特征进行合并，不同权利要求的内容并不意味着特征的组合不 可行和/或无优势。另夕卜，单数参考不排除复数。词语"一"、"一个""第一" "第二"等不排除多个。仅将权利要求中的参考数字提供作为说明示例， 不应该将其理解为以任何方式限制权利要求的范围。
权利要求
1、一种对感兴趣体积中的结构进行识别的方法，包括-以连续模式采集与所述结构相关的多个点，以及随后-将至少一个所述点配准到所述结构的先前采集的成像数据集。
2、 如权利要求l所述的方法，其中，所述配准包括 -通过分割对所述成像数据集进行处理，-计算3D二元掩模，-计算距离函数或者反映到所述感兴趣结构的表面的距离的其他计算度量， -将所述连续采集的多个点与所述成像数据集进行对准， -估计配准参数以相对于所述成像数据集对所述多个点进行更新， -将每个点与来自所述距离函数的对应值相比较， -根据所述估计的配准参数计算目标函数值，-通过利用基于共轭梯度的函数最小化对所述配准参数进行迭代更新，直到 所述配准参数的变化小于预定的容差，而将所述目标函数最小化。
3、 如权利要求2所述的方法，其中，所述分割包括通常可用于成像工作站 或者在图像分析软件中提供的用于可变形表面的算法、水平集分割方法、 或连通分量区域生长。
4、 如权利要求1所述的方法，包括-基于将所述EAM坐标系映射到所述成像坐标系的已知表变换，将所述连 续采集的多个点进行对准，-基于所述连续采集的点位置对从所述成像数据集内的所述对应位置处的 图像强度导出的值进行采样，以消除对分割的需要， -提供所有点位置处的强度或强度导出值的分布的强度直方图，以及 -基于所述强度直方图执行变换参数估计。
5、 如权利要求4所述的方法，包括使用从所述直方图导出的标准差和/或 熵在变换参数估计上进行迭代，以便使在所述多个点中的所有点上采样的 体素强度值的相似性最大化。
6、 一种对感兴趣体积中的结构进行识别的装置(40)，包括 -采集硬件(41)，其用于以连续模式采集与所述结构相关的多个点，以及 -配准硬件(42)，其用于随后将所述数据点中的每一个配准到所述结构的 成像数据集。
7、 如权利要求6所述的装置(50)，其中，所述配准硬件还包括 -处理硬件(51)，其用于通过分割对所述成像数据集进行处理，-第一计算硬件(52)，其用于计算3D二元掩模，-第二计算硬件(53)，其用于计算距离函数或者反映到所述感兴趣的结构 ，面的距离的其他计算度量，」对准硬件(54)，其用于将所述连续所采集的多个点与所述成像数据集进 行对准，-估计硬件(55)，其用于估计配准参数，以相对于所述成像数据集对所述 多个点进行更新，-比较硬件(56)，其用于将每个点与来自所述距离函数的对应值相比较， -第三计算硬件(57)，其用于根据所述估计的配准参数计算目标函数的值， 以及-最小化硬件(58)，其通过利用基于共轭梯度的函数最小化对所述估计的 配准参数进行迭代更新，直到所述估计的配准参数的变化小于预定的容差， 而将所述目标函数最小化。
8、 如权利要求6所述的装置，其中，所述多个点为电解剖标测点，所述成 像数据集为过程前或过程后采集的成像数据集。
9、 如权利要求8所述的装置，其中，所述采集硬件还包括配置为使用超声 或电磁进行点位置采集的导管。
10、 如权利要求6所述的装置，其中，所述装置为医学工作站。
11、 如权利要求6所述的装置，其中，将所述装置配置为执行如权利要求 1-5的任一个所述的方法。
12、 一种医学成像系统包括如权利要求6^f述的装置。
13、 一种其上嵌有计算机程序(60)的计算机可读介质，所述计算机程序 由计算机处理以对感兴趣体积中的结构进行识别，所述计算机程序包括 -采集代码段(61)，其用于以连续模式采集与所述结构相关的多个点，以 及-配准代码段(62)，其用于将至少一个所述点配准到所述结构的先前采集 的成像数据集。
14、如权利要求13所述的计算机可读介质，包括执行如权利要求1-5的任 一个所述的方法的至少一个代码段。
15、 如权利要求13所述的计算机可读介质，还包括用于控制进行连续点采 集的导管的控制代码段。
16、 一种对感兴趣体积中的解剖结构进行识别的方法，包括 -以连续模式采集与所述解剖结构相关的多个点，以及随后， -将至少一个所述点配准到过程前或过程后采集的所述结构的成像数据集。
全文摘要
本发明提供了一种对感兴趣体积中的结构进行识别的方法。所述方法包括以连续模式采集与所述结构相关的多个点，随后将所述点中的至少一个配准到所述结构的先前采集的成像数据集。本发明还提供了一种装置、系统以及计算机可读介质。本发明通过对标测系统进行修改，提供了对EAM点的更快采集，以便自动且连续地记录导管尖端位置，而不需要显式导航到心内膜，也不需要在心内膜上标注基准路标。
文档编号G06T7/00GK101501727SQ200780029311
公开日2009年8月5日 申请日期2007年8月6日 优先权日2006年8月8日
发明者R·曼茨克, R·陈, V·雷迪, Z·J·马尔察诺 申请人:皇家飞利浦电子股份有限公司