专利名称：用超声图像进行快速解剖绘图的制作方法
技术领域：
本发明一般地涉及超声成像，具体地讲，涉及根据超声图像进行解剖结构三维 (3D)重建的方法和系统。
背景技术：
多种用于心脏内超声成像的装置和方法为本领域所公知。例如，Biosense Webster Inc. (Diamond Bar, California)提供了用于实时生成3D超声图像的 CartoSound 系统和SoundStar 导管。可通过脉管系统插入心脏的SoundStar导管包括 位置传感器和相控超声换能器阵列。CartoSound系统可通过处理来自位置传感器和超声换 能器的信号来生成心室的3D图像。本领域已知用心脏内超声成像进行心内膜表面非接触式重建的若干方法。例如， PCT国际专利公布WO 00/19908描述了用于心脏内超声成像的可操纵换能器阵列，该专利 的公开内容以引用方式并入本文中。该阵列可形成超声束，该超声束由有源孔径引导到所 需方向。美国专利6，004，269描述了一种以装入导管中的超声装置为基础的声学成像系 统，该专利的公开内容也以引用方式并入本文中。该超声装置可将超声信号朝心脏内的内 部结构引导，以产生超声图像。作为另一个实例，PCT国际专利公布WO 99/55233描述了一种用于描绘患者心脏 的3-D表面的方法，该专利的公开内容以引用方式并入本文中。用训练数据形成3-D网格 模型，以用作一组患者心脏的原型形状。在不同图像平面上取得患者心脏的多个超声图像。 在每个图像上人工识别出解剖位置。将网格模型相对于预定的解剖位置与图像严格对齐。描绘轮廓时的此类人工辅助在基于超声图像的3D重建方法中是常见的。例如，美 国专利申请公布2006/0M1445描述了解剖结构的建模方法，其中用超声波传感器在超声 波传感器的相应的多个空间位置处获得解剖结构的多个超声图像，该专利的公开内容以引 用的方式并入本文中。在多个空间位置的每一处测量超声波传感器的位置和取向坐标。在 一个或多个超声图像中标出与解剖结构特征有关的所关注的轮廓。根据所关注的轮廓以及 测量的位置和取向坐标，构建解剖结构的三维(3D)模型。
发明内容
下文所述的本发明实施例提供了用2D超声图像的集合快速准确地进行3D表面重 建的方法。因此根据本发明的实施例提供了进行三维(3D)绘图(mapping)的方法，该方法包 括获取活体受试者体内腔体的多个二维OD)超声图像，其中2D图像在3D参照系中具有不 同的、相应的位置。在每个2D超声图像中，可识别出与腔体内部各位置对应的像素。将从 多个2D图像中识别的像素记录在3D参照系中，以确定与腔体内部对应的体积。处理器可 重建该体积的外表面，其代表腔体的内表面。在一些实施例中，获取多个2D图像包括将含有超声成像装置的探头插入体内，用位于体内的探头获取2D图像。在本发明所公开的实施例中，插入探头包括将导管插入受试 者的心脏，并且其中体积对应于心室的内部。插入导管通常包括将导管设置为可使得超声 成像装置位于第一心室内，其中除第一心室外，体积可以对应第二心室的内部，以使得重建 外表面包括生成第二心室内表面的3D图。通常，获取多个2D图像包括用含有位置换能器的超声探头获取2D图像，记录识别 的像素包括接收和处理与位置换能器相关的信号，以找出探头在3D参照系中的坐标，以及 用坐标将识别的像素记录在3D参照系中。在一些实施例中，在2D超声图像中识别像素包括将与具有低反射的体内位置相 对应的像素归类为腔体内部。通常，将像素进行归类包括设定阈值，以及将具有低于指定阀 值的相应灰度值的像素归类为腔体内部。在可供选择的实施例中，获取二维OD)超声图像 包括获取多普勒图像，其中颜色代表流动，并且其中识别像素包括根据像素相应的颜色值 将像素归类为腔体内部。在本发明所公开的实施例中，重建外表面包括对体积中的像素应用滚球算法 (ball-pivoting algorithm)。应用滚球算法可以包括用在体积上变化的重建分辨率处理像素。根据本发明的实施例，还提供了用于三维(3D)绘图的包括超声探头的设备，其能 够获取活体受试者体内腔体的多个二维OD)超声图像，其中2D图像在3D参照系中具有不 同的相应位置。处理器能够识别与腔体内部位置对应的每个2D超声图像像素，将从多个2D 图像识别的像素记录在3D参照系中，以确定与腔体内部对应的体积，以及重建代表腔体内 表面的体积外表面。根据本发明的实施例，另外提供了计算机软件产品，其包括可储存程序指令的计 算机可读介质，当处理器读取了指令时，可使处理器接收活体受试者体内腔体的多个二维 (2D)超声图像(这些2D图像在3D参照系中具有不同的相应位置)，以识别与腔体内部位 置对应的每个2D超声图像像素，并且将从多个2D图像识别的像素记录在3D参照系中，以 确定与腔体内部对应的体积，以及重建代表腔体内表面的体积外表面。通过以下结合附图的实施例的详细说明，将更全面地理解本发明。
图1为根据本发明实施例的用于超声成像的基于导管的系统的示意性说明图；图2为图1系统中所用导管的远端的示意性侧视图；图3为根据本发明实施例的用导管获取的超声图像的示意图；图4为流程图，其示意说明根据本发明的实施例的用超声图像进行快速解剖绘图 的方法；以及图5为根据本发明实施例制作的心室3D图的示意图。
具体实施例方式现在参见图1和图2，图中示意说明根据本发明实施例的基于导管的超声成像系 统20。图1为整个系统的说明图，而图2为系统中所用探头(如导管观)的远端的侧视图。 这里以举例说明的方式示出了该系统和导管，有助于理解下文所述的基于超声波的3D绘图方法。然而，这些方法并不限于基于导管的超声感测，还可以通过适当的变更，在体内或 体外使用其他类型的探头获取的2D或3D超声图像来应用这些方法。此外，这些方法不仅 可用于心脏，还可以用于其他解剖腔体的绘图。如图1所示，操作者22(如医生)将导管观插入患者沈体内，使得导管的远端穿 过脉管系统进入患者的心脏M。将导管的近端连接到控制台34上，此控制台通常包括具有 合适的信号处理电路和用户接口电路的处理器38。该处理器可接收和处理来自导管观的 信号，如下文所述。处理器38可包括一般用途的计算机处理器，其可以通过软件编程来实 现本文所述的功能。该软件可以通过(例如)网络以电子形式下载到处理器中。作为另外 一种选择或除此之外，该软件可以储存在有形的计算机可读储存介质中，例如光学、磁性或 电子存储介质。作为另外一种选择或除此之外，该处理器的至少某些功能可以通过数字信 号处理器(DSP)或专门的或可编程的硬件逻辑电路实现。通常，控制台34还让使用者能够观察并控制导管观的功能，以及查看并编辑利用 导管形成的图像。为了这些目的，该控制台包括显示器40和用户界面42。如图2所示，导管观的远端包括超声成像装置50，该装置用于生成体内的超声图 像。装置50通常包括换能器相控阵列52，如本领域所知，它可用来在扫描超声束平面(本 文中称为“束平面”或“图像平面”)中获取二维OD) “扇形”图像，其包含导管的纵向轴线。 这些换能器可接收从束平面中的对象反射的超声波，并根据反射的超声波输出信号。通常， 这些信号由穿过导管观的导线56传送至控制台34，控制台可以处理信号以形成并显示超 声图像和3D图，如下文所述。导管观的远端还包括位置传感器M，其可以生成表示体内导管位置(位置和取 向)的信号。根据这些位置信号，控制台34可以确定由成像装置50获取的每个扇形图像 的位置和取向。这样，处理器38能够确定在扇形图像中出现的对象的坐标，并记录和合并 在不同导管位置处获取的多个2D图像。在图示的实施例中，系统20使用磁性位置感测来确定导管观远端在心脏M内的 位置坐标。为了确定位置坐标，控制台34中的驱动电路36驱动磁场发生器30，使其在患 者沈体内产生磁场。磁场发生器30通常包括线圈，其放置在患者躯干下面的体外已知位 置处。这些线圈在包含心脏M的预定工作体积内产生磁场。传感器M可以包括(例如) 在导管观远端内的一个或多个线圈，其可以响应这些磁场而产生电信号。处理器38可处 理这些信号，以确定导管28远端的位置(位置和取向)坐标。控制台34可以用这些坐标 来驱动显示器40，使其显示导管的位置和状态。该位置感测和处理方法在由Biosense Webster Inc制造的CART0 系统中实 施。这种磁性位置感测在(例如)美国专利6，沈6，551中有详细描述，其公开内容以引用 方式并入本文中。将超声成像与磁性位置感测相结合的其他系统在美国专利6，690，963、 6，716，166和6，773，402中有所描述，这些专利的公开内容也以引用方式并入本文中。作为另外一种选择或除此之外，系统20可包括用于在患者沈体内操纵和操作导 管观的自动化机构(未示出)。在此类实施例中，处理器38根据位置感测系统提供的信号 生成用于控制导管运动的控制输入。尽管图1示出了具体的系统构型，但也可在本发明的可供选择的实施例中使用其 他系统构型。例如，可以用其他类型的位置换能器实施下文描述的方法，如阻抗型或超声位置传感器。如本文所用，术语“位置换能器”是指安装在导管观之上或之内的元件，该元件 可使控制台34接收指示该元件坐标的信号。因而该位置换能器可包括导管中的接收器，如 传感器M，其可根据换能器接收到的能量产生送给控制单元的位置信号；或者该换能器可 包括发射器，其发射探头外部的接收器可感测的能量。此外，在通过下文所述的方法进行绘 图和成像不仅可以使用导管，还可以使用其他类型的探头，可以在心脏和其他身体器官及 区域中使用，也可以使用体外的超声探头。图3为根据本发明的实施例由导管观获取的超声图像60的示意图。图像为2D扇 形的形式，其顶点在成像装置50处。如上文所述，控制台34可以根据从位置传感器M接 收到的信号，确定3D空间中的顶点位置和扇形取向。图像中的暗区域62、64对应于充满血 液而具有低反射的区域，例如心室。较亮的区域通常代表组织，如心脏内壁和外壁。如前文所述，操作者22可以操纵心脏M内的导管28，以便从不同位置、以不同的 方向获取图像。构成图像的反射不仅可以来自导管远端所在的心室，也可以来自其他的心 室和解剖结构。因此，例如可以将导管观插入右心房(比较容易通过腔静脉插入)，导管可 以从右心房或左心房并且可以从心室获取图像。图4为流程图，其示意说明根据本发明的实施例用超声图像进行快速解剖绘图的 方法。在该方法的每次重复中，在图像获取步骤70，传感器50获取2D超声扇形图像，其具 有图像60(图幻的一般形式。例如，可以在心动周期中的某个标注点(例如心脏收缩或心 脏舒张)选通图像获取，同时同步使用心电图(ECG)监测，或者作为另外一种选择，可以连 续获取图像而无需选通。在腔体识别步骤72，处理器38可识别由超声导管获取的每个2D 图像中所关注的心室内部(血池区域)。这些低反射“暗”区域可以通过(例如)对超声图 像的灰度级应用阀值来识别。阈值可以自动或手动设置。可以用本领域已知的任何合适的 方法自动选择阈值，如Otsu法，其根据像素灰度值的集群选择阈值。作为另外一种选择，步骤70中的图像获取可以使用其他超声成像模式，并可以对 步骤72中用来识别腔体的方法作相应的修改。例如，可以用多普勒成像技术来获取2D图 像，例如本领域已知的彩色多普勒、能量多普勒或组织成像多普勒。此类技术用图像颜色 (通常称为伪彩色)代表流动。在彩色多普勒中，图像中将血流区域着色，而组织区域不着 色。在这种情况下，在步骤72，可以将颜色值大于某个阀值的像素识别为属于血池区域。另 一方面，在组织成像多普勒中，将组织区域着色而血液不着色，以使颜色值低于某个阀值的 像素将被识别为属于血池区域。无论使用哪种成像模式，在二值化步骤74，处理器38都会用阈值将2D灰度或彩色 图像转换成二值图像。在二值图像中，将值为“0”的像素归类为血液区域，而值为“1”的像 素属于组织。可以使用其他图像处理操作，以改善血液与组织之间区分的精确度。例如，可 以通过连续使用形态学腐蚀和膨胀来移除组织区域内的小的暗区域，这些区域会被错误地 识别为血池区域。在图像记录步骤76中，处理器38找出获取的2D图像的顶点位置和取向。如上文 所述，处理器可以根据位置传感器M输出的信号计算位置和取向坐标。在获知顶点位置和 图像取向的条件下，处理器可以计算二值化图像中每个像素在磁场发生器30的固定3D参 照系中的3D坐标，从而将2D图像像素记录在3D体积中。在获取完成步骤78中，获得每个2D图像后，操作者移动心脏中的导管顶端，并重复上述2D图像获取和处理步骤，直到处理器在所关注的体积内收集到足够数量的二值像素值。在表面重建步骤80中，处理器38收集值为“0”的像素(血液)，并找出与这些像 素接界的外表面。与心室中的血池接界的外表面为包围该心室的心脏壁内表面。因此，通 过找出包围血池的外表面，处理器可有效地构建所考虑的心室的3D图。任选地，处理器可 以在重建3D表面的同时继续获取2D图像，从而可以逐渐改进该图。将3D域中的血池作为表面重建的基础，可以为解决图像分割问题提供快速有效 的方法。基于此类内点的集合可以用各种算法来重建某个体积的外表面。例如，处理器38 可以使用美国专利6，968，299中所述的滚球算法，该专利的公开内容以引用方式并入本文 中。该算法可以通过在点云上方“旋转”具有一定半径的球，推算出插入给定点云的三角网 格。用该方法找出的三角形的顶点可限定云的外表面。为了减小计算量，不是所有的“0”像素都必须用于建立模型，处理器38可以在不 同区域使用不同的分辨率。上文提到的滚球算法适用于通过在图的不同区域使用不同的球 半径来进行这样的可变分辨率绘图。通常仅在血液-组织界面区域(沿着血池的边缘)需 要使用高分辨率。该区域可在二维图像中自动或手动加以识别。然后处理器仅在组织边界 附近进行高分辨率重建，其他地方采用低分辨率。换句话讲，处理器使用较少的2D图像中 远离组织的区域中的像素，而使用较多的组织界面(通常在大约Imm的厚度内)区域中的 像素。处理器38以上文所述方式重建和显示3D图之后，使用者可以使用图像编辑工具 通过用户界面42 (例如)切割和/或改变模型，以改正伪像并移除不受关注的特征。图5为心室的3D图90的示意图，这是一种可用根据本发明实施例的上述方法制 作的3D图。在这种情况下，该图示出了心脏的左心室，包括肺静脉的一些部分。这样的图 可以用于各种目的，例如用作电解剖图(用上文所述的CARTO系统制成的那种图)的起点 或用作全3D超声或层析图像的片段的起点。应当理解以上述的实施例仅是举例方式的援引，本发明并不限于上文具体示出和 描述的内容。并且，本发明的范围包括上述各种特征的组合和子组合、以及本领域技术人员 在阅读上述说明书时可能想到的并且现有技术中未公开的变型形式和修改形式。
权利要求
1.一种三维(3D)绘图方法，包括获取活体受试者体内腔体的多个二维OD)超声图像，所述2D图像在3D参照系中具有 不同的相应位置；在每个所述2D超声图像中识别与所述腔体内部中的位置相对应的像素；将从所述多个2D图像识别的所述像素记录在所述3D参照系中，以确定与所述腔体内 部对应的体积；以及用处理器重建所述体积的外表面，所述外表面代表所述腔体的内表面。
2.根据权利要求1所述的方法，其中获取所述多个2D图像包括将含有超声成像装置的 探头插入所述体内，以及用所述体内的所述探头获取所述2D图像。
3.根据权利要求2所述的方法，其中插入所述探头包括将导管插入所述受试者的心 脏，并且其中所述体积对应于所述心脏的心室内部。
4.根据权利要求3所述的方法，其中插入所述导管包括将所述导管设置为可使得所述 超声成像装置位于所述心脏的第一心室中，并且其中除所述第一心室外，所述体积还对应于所述心脏的第二心室内部，使得重建所述 外表面包括生成所述第二心室的内表面的3D图。
5.根据权利要求1所述的方法，其中获取所述多个2D图像包括用含有位置换能器的超 声探头获取所述2D图像，并且其中记录所述识别的像素包括接收和处理与所述位置换能器相关的信号，以找出所述 探头在所述3D参照系中的坐标，以及用所述坐标将所述识别的像素记录在所述3D参照系 中。
6.根据权利要求1所述的方法，其中识别所述2D超声图像中的所述像素包括将与所述 体内具有低反射的位置相对应的像素归类为属于所述腔体的内部。
7.根据权利要求6所述的方法，其中将所述像素进行归类包括设定阈值，以及将具有 低于指定阈值的相应灰度值的像素归类为属于所述腔体的内部。
8.根据权利要求1所述的方法，其中获取所述二维OD)超声图像包括获取多普勒图 像，其中颜色代表流动，并且其中识别所述像素包括根据所述像素的相应颜色值将像素归 类为属于所述腔体的内部。
9.根据权利要求1所述的方法，其中重建所述外表面包括对所述体积内的所述像素应 用滚球算法。
10.根据权利要求9所述的方法，其中应用滚球算法包括用在体积上变化的重建分辨 率处理所述像素。
11.一种用于三维(3D)绘图的设备，包括超声探头，所述超声探头被配置成获取活体受试者体内腔体的多个二维OD)超声图 像，所述2D图像在3D参照系中具有不同的相应位置；以及处理器，所述处理器被配置成可在每个所述2D超声图像中识别与所述腔体内部中的 位置相对应的像素，并可将从所述多个2D图像识别的像素记录在所述3D参照系中，以确 定与所述腔体内部对应的体积并重建所述体积的外表面，所述外表面代表所述腔体的内表
12.根据权利要求11所述的设备，其中所述超声探头配置成插入所述体内并获取所述体内的2D图像。
13.根据权利要求12所述的设备，其中所述探头包括可插入所述受试者心脏的导管， 其中所述导管包括超声成像装置，并且其中所述体积对应于所述心脏的心室内部。
14.根据权利要求13所述的设备，其中所述导管配置成插入所述心脏，以使所述超声 成像装置位于所述心脏的第一心室内，并且其中除所述第一心室外，所述体积还对应于所述心脏的第二心室内部，使得重建所述 外表面包括生成所述第二心室的内表面的3D图。
15.根据权利要求11所述的设备，其中所述超声探头包括位置换能器，并且其中所述处理器被配置成可接收和处理与所述位置换能器相关的信号，以找出所述探 头在所述3D参照系中的坐标，并用所述坐标将所述识别的像素记录在所述3D参照系中。
16.根据权利要求11所述的设备，其中所述处理器被配置成可将与所述体内具有低反 射的位置相对应的像素归类为属于所述腔体的内部。
17.根据权利要求17所述的设备，其中所述处理器被配置成可设定阈值，并将具有低 于指定阈值的相应灰度值的像素归类为属于所述腔体的内部。
18.根据权利要求11所述的设备，其中所述二维OD)超声图像包括多普勒图像，其中 颜色代表流动，并且其中所述处理器被配置成可根据所述像素的相应颜色值将像素归类为 属于所述腔体的内部。
19.根据权利要求11所述的设备，其中所述处理器被配置成可对所述体积中的所述像 素应用滚球算法，以重建所述体积的所述外表面。
20.根据权利要求19所述的设备，其中所述导管配置成可用在体积上变化的重建分辨 率来处理所述像素。
21.一种计算机软件产品，包括储存程序指令的计算机可读介质，所述指令当由处理器 读取时，使所述处理器接收活体受试者体内腔体的多个二维OD)超声图像，所述2D图像在 3D参照系中具有不同的相应位置，以在每个所述2D超声图像中识别与所述腔体内部位置 对应的像素，并且将从所述多个2D图像识别的所述像素记录在所述3D参照系中，以确定与 所述腔体内部对应的体积并重建所述体积的外表面，所述外表面代表所述腔体的内表面。
22.根据权利要求21所述的产品，其中所述体积对应于所述心脏的心室内部。
23.根据权利要求22所述的产品，其中用所述心脏的第一心室中的超声成像装置获取 所述2D超声图像，并且其中除所述第一心室外，所述体积还对应于所述心脏的第二心室的 内部，使得重建所述外表面包括生成所述第二心室的所述内表面的3D图。
24.根据权利要求21所述的产品，其中用含有位置换能器的超声探头获取所述2D图 像，并且其中所述指令使所述处理器处理与所述位置换能器相关的信号，以找出所述探头 在所述3D参照系中的坐标，并用所述坐标将所述识别的像素记录在所述3D参照系中。
25.根据权利要求21所述的产品，其中所述指令使所述处理器将与所述体内具有低反 射的位置相对应的像素归类为属于所述腔体的内部。
26.根据权利要求21所述的产品，其中所述指令使所述处理器对所述体积中的所述像 素应用滚球算法，以重建所述外表面。
27.根据权利要求沈所述的产品，其中所述指令使所述处理器用在体积上变化的重建 分辨率处理所述像素。
全文摘要
本发明名称为“用超声图像进行快速解剖绘图”。本发明公开了一种三维(3D)绘图方法，所述方法包括获取活体受试者体内腔体的多个二维(2D)超声图像，所述2D图像在3D参照系中具有不同的相应位置。在每个所述2D超声图像中，与所述腔体内部的一些位置对应的像素被识别。将从所述2D图像识别的所述像素记录在所述3D参照系中，以确定与所述腔体内部相对应的体积。重建所述体积的外表面，所述外表面代表所述腔体的内表面。
文档编号G06T17/00GK102106743SQ20101062092
公开日2011年6月29日 申请日期2010年12月23日 优先权日2009年12月23日
发明者A·C·阿尔特曼, A·戈德伯格, N·基利姆 申请人:韦伯斯特生物官能(以色列)有限公司