CN105574878B - 使用局部加权拟合的图像配准 - Google Patents

使用局部加权拟合的图像配准 Download PDF

Info

Publication number
CN105574878B
CN105574878B CN201510971118.3A CN201510971118A CN105574878B CN 105574878 B CN105574878 B CN 105574878B CN 201510971118 A CN201510971118 A CN 201510971118A CN 105574878 B CN105574878 B CN 105574878B
Authority
CN
China
Prior art keywords
image
map
coordinates
points
weights
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510971118.3A
Other languages
English (en)
Other versions
CN105574878A (zh
Inventor
A.普赖斯
R.塔尔
Z.基特拉罗
Z.帕特
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Biosense Webster Inc
Original Assignee
Biosense Webster Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Biosense Webster Inc filed Critical Biosense Webster Inc
Publication of CN105574878A publication Critical patent/CN105574878A/zh
Application granted granted Critical
Publication of CN105574878B publication Critical patent/CN105574878B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42—Details of probe positioning or probe attachment to the patient
    • A61B8/4245—Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • A61B8/4254—Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient using sensors mounted on the probe
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0044—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/757—Matching configurations of points or features
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G06T2207/10136—3D ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular

Abstract

一种用于图像配准的计算机实施的方法包括接受器官的三维(3‑D)图像。该图像包括具有相应图像坐标的多个第一图像点。接受定义器官上的生理学参数值的分布的图。该图包括具有相应图坐标的多个第二图点。使置信水平与一个或多个图坐标和图像坐标相关联。根据各个图坐标和图像坐标的置信水平，对至少一些图点和图像点分配权重。通过基于权重计算在图坐标和图像坐标之间的几何变换来将3‑D图像与图拟合。

Description

使用局部加权拟合的图像配准
技术领域
本发明一般而言涉及医学成像和绘图，尤其涉及利用三维(3-D)图像对生理学参数的图进行配准的方法和系统。
背景技术
在各种医学应用中，身体部分或器官的计算机生成的重建被显示给用户。例如，美国专利6,456,867描述了一种对体积的图进行重建的方法，该专利的公开被结合于此以作参考。该方法包括：确定在该体积的表面上位置的坐标，并且产生在接近于所确定的位置的3-D空间中限定重建表面的点的网格。对于该网格上的每个点，根据在网格上的一个或多个点与一个或多个位置之间的位移来定义各自的向量。然后根据各自的向量通过充分地移动网格上的每个点来调整该重建表面，使得该重建表面变形以与该表面的构型类似。
在一些应用中，来自多模态或成像源的图像被相互拟合并被呈现。例如，美国专利5,568,384描述了一种用于把3-D多模态图像集合并成单个合成图像的方法，该专利的公开被结合于此以作参考。最初使用半自动分段技术从待匹配的从两个或更多个不同图像中提取表面。这些表面被表示为具有待匹配的共同特征的轮廓。对一个表面图像执行距离变换，并且使用该距离图像来产生用于匹配过程的成本函数。几何变换包括三维平移、旋转和缩放以适应不同位置、方向和大小的图像。该匹配过程包括：搜索多参数空间，并且调整一个表面或多个表面以找到它们之间的最佳拟合，这使成本函数最小化。
发明内容
一些医学方法和系统测量并呈现用于显示在成像的器官上某个生理学参数的分布的图。例如，显示在心内膜(心脏的内表面)上局部电位的分布的电解剖图在各种心脏过程中是有用的。
有时期望显示参数图连同器官的3-D图像，例如基于磁共振成像(MRI)、计算机断层扫描(CT)、X射线或超声的图像。为了共同显示该图和3-D图像，器官的这两种视图应该首先被配准到一个共同的坐标系。在一些情况下，参数图和/或3-D图像的坐标可以具有不同的置信水平。具有低置信水平的图或图像的点或区域会对配准过程的精度具有不利影响。
为了改善配准过程的精度，本发明的实施例提供用于在考虑到图坐标和/或图像坐标的不同置信水平时将参数图与3-D图像进行配准的方法和系统。3-D图像通常包括具有相应图像坐标的多个图像点。在一些实施例中，图像点定义了成像的器官或其部分的3-D实体积模型、表面模型或3-D线框模型。在一些实施例中，根据图和图像坐标的置信水平，对图点与图像点分配权重。然后通过基于权重计算和应用几何变换来将该图和3-D图像拟合在一起。
在此描述的一些示例性实施例中，使用磁位置跟踪系统来获得参数图。在这些实施例中，诸如导管之类的探头装有传感器，例如电极，其测量在成像的器官中不同位置的生理学参数。位置跟踪系统在进行测量时用于跟踪探头的坐标，以便产生用于指示在器官上参数值的分布的图。
因此，根据本发明的实施例，提供一种用于图像配准的计算机实施的方法，该方法包括：
接受器官的三维(3-D)图像，该图像包括具有相应图像坐标的多个第一图像点；
接受定义器官上的生理学参数值的分布的图，该图包括具有相应图坐标的多个第二图点；
使相应置信水平与一个或多个图坐标和图像坐标相关联；
根据各个图坐标和图像坐标的置信水平，对至少一些图点和图像点分配权重；以及
通过基于权重计算在图坐标和图像坐标之间的几何变换来将3-D图像与图拟合。
在一个实施例中，接受图包括使用位置跟踪系统来测量各个图坐标的生理学参数的值。测量这些值和图坐标可以包括：移动探头，该探头包括位置传感器和在器官中的多个位置上测量生理学参数的传感器；以及对在这些位置的生理学参数和位置传感器的位置执行相应的多次测量。
在另一实施例中，该方法包括：使用几何变换将生理学参数的值叠加到3-D图像上，并且将所叠加的值和3-D图像显示给用户。
在又一实施例中，将3-D图像与图拟合包括：定义在至少一些图像点和相应的图点之间的距离，将分配给各个图点和图像点的权重应用于这些距离以产生加权距离，并且基于这些距离来计算累积的加权距离量度，该量度指示在将3-D图像与图拟合中几何变换的精度。
在再一实施例中，应用权重包括将距离的平方乘以各自的权重以产生加权距离，以及计算累积的加权距离量度包括对加权距离求和。
在一个实施例中，将3-D图像与图拟合进一步包括修改几何变换以便最小化累积的加权距离量度。
在另一实施例中，将3-D图像与图拟合包括：识别出现在3-D图像中的器官外部的解剖标志，在该解剖标志的位置获得相应图点，并且配准3-D图像中的解剖标志的图像坐标与所获得的图点的图坐标。
在又一实施例中，使置信水平相关联包括接受来自用户的置信水平。另外或可选择地，使置信水平相关联包括基于图点和图像点来自动地确定置信水平。自动地确定置信水平可以包括：检测在图坐标和图像坐标中的至少一个中的失真，并且根据该失真的水平来分配置信水平。另外或可选择地，自动地确定置信水平包括根据下述中的至少一个来分配置信水平：
在图中相应的一个或多个图坐标的密度；
在相应的图坐标和图像坐标附近的电组织特征；
在相应的图坐标和图像坐标附近的生理学组织特征；
相应的图坐标和图像坐标距在图和3-D图像中都识别出的配准标志的距离；以及
配准标志的生理学特征。
在一个实施例中，接受3-D图像包括接受由至少一些图像点所限定的器官的实体模型、表面模型和线框模型中的至少一个。
根据本发明的实施例，另外提供用于图像配准的装置，包括：
接口，其被布置成接受器官的三维(3-D)图像，该图像包括具有相应图像坐标的多个第一图像点，以及被布置成进一步接受定义器官上的生理学参数值的分布的图，该图包括具有相应图坐标的多个第二图点；以及
处理器，其被布置成使相应置信水平与一个或多个图坐标和图像坐标相关联，根据各个图坐标和图像坐标的置信水平对至少一些图点和图像点分配权重，并且通过基于权重计算在图坐标和图像坐标之间的几何变换来将3-D图像与图拟合。
根据本发明的实施例，还提供一种用于图像配准的计算机软件产品，该产品包括一种其中存储指令的计算机可读介质，所述指令在由计算机读取时使得计算机：接受预先获得的器官的三维(3-D)图像，该图像包括具有相应图像坐标的多个第一图像点；接受定义器官上的生理学参数值的分布的图，该图包括具有相应图坐标的多个第二图点；使相应置信水平与一个或多个图坐标和图像坐标相关联；根据各个图坐标和图像坐标的置信水平，对图点和图像点分配权重；并且通过基于权重计算在图坐标和图像坐标之间的几何变换来将3-D图像与图拟合。
附图说明
根据与附图一起采用的本发明实施例的后面的详细描述，将更充分地理解本发明，在附图中：
图1A是根据本发明的实施例用于磁位置跟踪的系统的示意性图示；
图1B是根据本发明的实施例示意性说明显示处理器的框图；
图2是根据本发明的实施例的导管的示意性图示；
图3A和图3B是根据本发明的实施例示意性说明投射到心脏的三维(3-D)图像上的电解剖图的图；以及
图4是根据本发明的实施例示意性说明用于配准电解剖图与预先获得的3-D图像的方法的流程图。
具体实施方式
系统描述
图1A是根据本发明的实施例用于对患者的心脏24进行成像和绘图的磁位置跟踪系统20的示意性图示。图1A的系统用于说明一个示例性系统，该系统可以用于获得并呈现参数图。尽管在此描述的实施例主要涉及使用磁跟踪系统所获得的电解剖图，但是可以使用现有技术中已知的任何其它手段来获得参数图，例如使用基于超声或阻抗的跟踪系统。在美国专利5,983,126、6,456,864和5,944,022以及在2005年1月7日提交的美国专利申请11/030,934中描述了示例性的基于阻抗的跟踪系统，所述专利和专利申请的公开被结合于此以作参考。
在系统20中，医生通过静脉或动脉将导管28插入心脏24的室中。系统20确定并显示导管28的位置和方向坐标。尽管下文描述的实施例主要涉及心脏的成像，但是在此描述的方法和系统可以用于对任何合适的器官进行成像。因此，在本发明的上下文中以及在权利要求中，术语“器官”用于描述任何成像的身体部分或组织，其可以包括整个器官或其部分。
系统20包括场发生器，例如场发生线圈30，其位于患者外部的固定的已知位置。线圈30在心脏24周围的预定工作体积中产生磁场。安装在导管28的远端的位置传感器检测在其附近的磁场。
响应于所检测的场，位置传感器产生位置信号并将其传送至控制台34。控制台包括跟踪处理器36，该跟踪处理器36基于由位置传感器所发送的位置信号来计算导管28相对于线圈30的位置和方向。使用显示器44将导管28的位置和方向坐标显示给医生。
系统20包括显示处理器42，其接受由位置跟踪系统所测量的信息以及外部提供的预先获得的3-D图像。显示处理器42配准参数图与预先获得的3-D图像，并将配准的图和图像呈现给医生，正如下面将要解释的。
图1B是根据本发明的实施例示意性说明显示处理器42的细节的框图。处理器42包括外部接口45，其接受预先获得的3-D图像和参数图。处理器42还包括配准处理器46，其执行在此所述的配准方法，并将配准的3-D图像和参数图输出给显示器44。使用诸如键盘、鼠标或跟踪球之类的用户输入装置47，处理器42与诸如医生之类的用户进行交互。如上所述，显示处理器42可以接受来自任何合适的系统或源的3-D图像和参数图。
通常，跟踪处理器36和显示处理器42包括通用计算机，其被以软件进行编程以执行在此描述的功能。该软件可以例如经由网络以电子的形式下载到处理器，或者它可以可选择地以有形介质比如CD-ROM被提供给处理器。处理器36和42可以是独立单元。可选择地，处理器36和42的功能可以相互组合，或者与系统20的其它元件集成。
图2是根据本发明的实施例的导管28的远端的示意性图示。导管28包括位置传感器50，其检测由线圈30所产生的磁场。另外，导管28包括电极52，其检测在其附近的生理学参数。
例如，在一些实施例中，电极52检测在其附近的局部电位。在这些实施例中，使导管28的远端非常接近在心内膜(心脏的内表面)上的各种点。然后使用电极52来检测在心内膜表面上的局部电位。在可选择实施例中，导管28可以包括一个或多个电极52，以用于检测任何其它合适的生理学参数，例如温度和压力。在一些实施例中，导管可以装有用于检测在其附近的特定化学成分的化学传感器。
使用穿过导管28的合适电缆54或者使用无线连接(未示出)，将位置传感器50和电极52连接到控制台34。跟踪处理器36使测量的导管28的位置坐标与电极52所测量的参数相关联。由此，通过移动导管通过心脏中的多个位置，医生能够获得参数的三维图。
在图1A和图2的示例性结构中，例如系统20获得并显示电解剖图，该电解剖图显示在心内膜表面上局部电位的分布。这种电解剖图可以用于检测和/或预测各种心脏失常，以及用于设计相关的医学过程，例如消融。在可选择实施例中，可以使用其它类型的电极和传感器来获得感兴趣的其它参数的图。通常，参数图包含多个图点。每个图点包括参数的特定值，并且与相应图坐标相关联。
图1A和图2的系统结构是为了概念清楚而选择的示例性结构。在此描述的方法可以与任何其它合适的位置跟踪系统或方法结合使用。例如在PCT专利公布WO 96/05768、美国专利5,391,199、5,443,489、6,690,963、6,239,724、6,618,612和6,332,089以及美国专利申请公布2002/0065455 A1、2003/0120150 A1和2004/0068178 A1中描述了一些用于磁位置跟踪的系统和方法，其公开都被结合于此以作参考。如上所述，可选择地，可以使用任何其它合适的方法或系统来产生参数图。
参数图与图像的配准
在一些实际情况下，期望显示投射到所绘图的器官的3-D图像上的参数图。在大多数情况下，3-D图像在绘图过程之前被预先获得，并且因此被称为“预先获得的”图像。可选择地，可以与绘图过程同时或者之后获得3-D图像。无论如何，组合了图和图像信息的双模态显示通常提高了医生对参数图的理解。例如，使用这种共同显示可以更容易定位和检查在成像器官中感兴趣的特征。
3-D图像可来源于磁共振成像(MRI)系统、计算机断层扫描(CT)系统、X射线成像系统、超声成像系统或者任何其它合适的成像系统或源。由于来源于不同系统的3-D图像通常具有不同的坐标系，所以在显示它们之前，成像的器官的两个视图应当进行配准，或者相互拟合。
通常，3-D图像作为多个图像点被表示和提供给显示处理器42。每个图像点包括一个值(例如灰度或颜色)，并且与相应的图像坐标相关联。在一些情况下，成像的器官或其部分可以在3-D图像中被表示为一个由至少部分图像点所限定的实体积模型、表面模型或线框模型。在本专利申请的上下文中以及在权利要求中，术语“图像点”用于描述3-D图像的任何合适的表示，例如由图像点限定的各种体积、表面和线框的表示。在此描述的配准方法可以用于配准参数图与具有任何这种表示的3-D图像。3-D图像通常作为数据文件被提供给显示处理器42，并且可以使用任何方便的数字手段例如网络连接或磁存储装置、或者使用模拟手段例如帧接收器来提供。
配准参数图与3-D图像的过程可以被描述为寻找一种将图与图像最佳拟合的几何变换的任务。这种几何变换可以包括图像坐标的平移、旋转和/或缩放，以便最佳拟合该图。可选择地，几何变换同样可以被定义为将参数图拟合到3-D图像上，或者将图和图像都变换到任意的参考坐标系。下面描述的方法可以用于确定任何这种几何变换。
在一些实施例中，使用的两个模态中之一包括实时模态，而另一个模态包括非实时模态。例如，由位置跟踪系统实时测量的参数图可以与预先获得的MRI或CT图像配准。在这种实施例中，尽管并非是强制性的，但是通常优选的是将几何变换限定为将非实时模态的坐标系平移到实时模态的坐标系。
配准过程的质量常常用累积距离量度来评估，该累积距离量度量化在图和3-D图像之间的距离。例如，在一些实施例中，表示为Mi的每个图点与表示为Ii的对应图像点相关联。在图点Mi和其对应图像点Ii之间的距离被表示为di。距离di通常表示在图点和图像点的坐标之间计算的三维欧氏距离。由
Figure GDA0002276880590000071
给出一个示例性累积距离量度，其中i在相应图点和图像点的对上进行求和。因此，图与3-D图像的配准包括确定一个最小化该累积距离量度D的几何变换。
然而，在一些情况下，关于所讨论的点的位置的精度，图点和/或图像点的坐标具有不同的置信水平。例如，参数图的一些区域可能具有高密度图点，从而导致较高的置信水平。该图的其它区域可能具有较少的图点，从而常常导致较低的置信水平。
也可以根据在图或3-D图像中检测到的失真来分配置信水平(例如，在特定图像点中或者在3-D图像的实体、表面或线框表示的区域中)。例如，显示出坐标值高度变化的图像区域可以被假定包含较高水平的测量噪声，并且因此与平稳变化的区域相比具有较低的置信水平。另外或可选择地，已知由高度时间变化表征的成像器官的区域(例如其形状随时间显著和/或快速改变的区域)也可以被认为是具有较低的置信水平。
图中的失真也可以包括图的一个区域的几何形状与3-D图像中的对应区域相比的变形。在上面引用的美国专利6,456,867中也描述了一些示例性方法，这些方法可以用于在所测量的图中确定低置信水平的区域。另外或可选择地，对图点分配置信水平可以基于任何其它合适的标准或标准的组合。
3-D图像中的图像点或区域也可能具有与它们的位置精度相关的不同置信水平。例如，3-D图像的区域可能具有一个错误分段(例如，图像中遗漏的器官的区域或者与靶器官结合的其它器官或器官部分的剩余部分)。3-D图像的一些区域可能具有较低的空间或时间分辨率。
具有低置信水平的图点和/或图像点可能将错误引入几何变换中，从而导致非最佳配准。因此可以认识到，在评估几何变换时考虑置信水平可以显著地提高配准的质量。可以防止低置信水平的区域使变换失真，并且可以加强较高置信的区域。
为了降低具有低置信水平的图和图像点的影响，在此描述的方法根据其置信水平来对图和/或图像点分配权重。然后将权重应用于为提供累积距离量度而求和的距离。所得到的量度被称作加权累积距离量度。通常，该权重与置信水平成比例，使得具有低置信水平的图和图像点被分配较低的权重，并且反之亦然。
例如，在一些实施例中，每个图点Mi基于该图点的置信水平而被分配一个表示为Wi的权重。在这些实施例中，由下式给出一个示例性加权累积距离量度：
Figure GDA0002276880590000081
可选择地，也可以使用其它合适的基于应用于距离的权重的累积距离量度。确定一个最小化加权累积距离量度的几何变换考虑到图和图像点的不同置信水平，提供了参数图与预先获得的3-D图像的高精度配准。
在一些实施例中，配准过程可以利用位于器官外的解剖标志。在这些实施例中，医生识别出现在3-D图像中的器官外部的一个或多个解剖特征，以用作配准标志或锚点。例如，当对心脏绘图时，医生可以选择主动脉上的一个或多个点，例如主动脉弓的主血管分支。作为另一个例子，医生可以选择下或上腔静脉的开口。然后医生操纵导管28以使它的末端接触所讨论的特征，并且在显示器44上指示预先获得的3-D图像的相应位置。
通过配准该标志，医生有效地识别已知相互匹配的一个图点和一个相应图像点。当使用主动脉或腔静脉作为标志时，医生在通过血管将导管插入到心脏中时可以配准该标志。因此，当导管28进入心脏以开始绘图过程时，已经用预先获得的3-D图像预配准了电解剖图。
可以由系统自动地或者由用户手动地分配权重。在一些实施例中，可以在过程期间的任何时候执行手动权重分配。自动权重分配可能依赖于相关区域中的组织特性。当使用外部标志时，权重分配可能依赖于该标志距所讨论的图或图像点的距离。例如，下腔静脉中的标志由于其柔软和松软特性以及其距左心房的距离而可以被自动地分配一个较低的权重。
图3A和3B是根据本发明的实施例示意性说明配准左心房和相关血管的电解剖图56与三维(3-D)图像58的过程的图。图像58包括由CT扫描器所产生的左心房和肺静脉的分段的预先获得的图像。
图3A示出在配准过程之前的图和3-D图像。如在图中可以看到的，图和3-D图像未被配准到一个共同坐标系，并且通常相对于彼此具有平移、旋转和缩放偏差。
如上所述，电解剖图56包括一个使用导管和磁位置跟踪系统所获得的电解剖图。该图包括多个图点60以及相应的图坐标，每个图点相对于参考注释具有一个局部活化值(以毫秒表示)。不同的局部活化值在图中用不同的阴影图案来示出。医生在其中已经执行了射频(RF)消融的位置用圆圈标记。
3-D图像包括多个图像点62，每个图像点具有相应的图像坐标。在电解剖图56和3-D图像58中都识别出表示为LAA、MPV、LIPV和MVA的四个解剖标志。这四个标志分别指示左心耳、医学肺静脉、左下肺静脉和二尖瓣环的位置。注意，所有的标志都位于靶器官即左心房的外部。
图3B示出在根据下面图4的方法执行配准过程之后的电解剖图56和3-D图像58。如在图中可以看到的，现在图和3-D图像用一个共同的坐标系配准并且一起显示。
配准方法描述
图4是根据本发明实施例示意性说明用于配准电解剖图与预先获得的3-D图像的方法的流程图。由显示处理器42所执行的图4的方法通过在评估并最小化加权累积距离量度时给不同图和图像点分配不同权重，来考虑到在参数图和/或3-D图像中不同的置信水平。
该方法从处理器42在图像接受步骤70接受预先获得的器官的3-D图像开始。该3-D图像由多个图像点62和相关坐标表示。如上所述，至少一部分图像点可以在一些情况下限定成像的器官或其部分的3-D实体、表面或线框模型。处理器42在图像接受步骤72进一步接受器官的电解剖图。在一些实施例中，如上所述，该图由磁位置跟踪系统20产生。该图由多个图点60和相关坐标表示。使用接口45来接受3-D图像和参数图并将其提供给配准处理器46。
对于至少一些图点和/或图像点，在置信分配步骤74，配准处理器46确定相应图和/或图像坐标的置信水平。在一些实施例中，由医生或者由另一用户将一些或全部的置信水平指示给处理器46。例如，医生可以检查显示在显示器44上的图和/或3-D图像，并且使用输入装置47来标记低的或高的置信水平的区域或点。
另外或可选择地，与图相关的一些或全部的置信水平可以由在执行位置测量期间或之后产生图的系统(例如位置跟踪系统)自动地确定，并且与电解剖图一起被提供给处理器46。进一步，另外或可选择地，一些或全部的置信水平可以由处理器46使用图和/或3-D图像基于不同标准例如上述的标准自动地计算。
处理器46在变换计算步骤76计算使图与3-D图像最佳拟合的几何变换。如上所述，每个图点与相应的图像点相关，并且在这两点之间基于它们的坐标定义了距离。为了减少具有低置信水平的图和/或图像点的影响并且增加具有高置信水平的图和/或图像点的影响，处理器46基于置信水平给这些距离分配权重。处理器46评估加权的累积距离量度，例如在上面的公式[2]中给出的量度。
在配准步骤78，处理器46确定最小化加权累积距离量度的几何变换，由此使图与3-D图像最佳拟合。所得到的几何变换通常包括应用于图坐标、图像坐标、或二者的平移、旋转和/或缩放操作的组合。然后使用显示器44将配准的图和3-D图像显示给医生，例如在上面的图3B中所示。
尽管在此描述的实施例主要涉及一种使用装备有电位检测电极的导管的心脏成像系统，但是在此描述的方法可以与产生包括与坐标相关的测量参数的图的任何其它合适的系统一起使用。参数可以在本地或远程测量。例如，装有超声传感器和位置传感器的导管可以用于产生生理学参数例如组织密度的图。这些参数图可以利用上文所述的方法与预先获得的3-D图像配准。
另外或可选择地，测量参数的传感器可以被安装在任何其它合适的体内探头例如内窥镜中。进一步，另外或可选择地，参数的测量可以使用外部传感器例如外部超声换能器来执行。在此描述的方法可以用于处理任何其它合适的器官、组织或其部分的图像和参数图，例如主动脉(AAA)、冠状动脉、脊柱、神经学上的或任何其它的腔，其中执行微创过程。
因此将会认识到，上述的实施例作为例子被引用，并且本发明不限于上文已经特别示出和描述的内容。而是，本发明的范围包括上文所述的各种特征的组合和子组合、以及本领域技术人员一阅读以上描述就将想到的并且未在现有技术中公开的其变化和修改。

Claims (26)

1.一种用于图像配准的计算机实施的方法，包括：
接受器官的三维（3-D）图像，所述图像包括具有相应图像坐标的多个第一图像点；
接受定义所述器官上生理学参数的值的分布的图，所述图包括具有相应图坐标的多个第二图点；
使相应的置信水平与所述图坐标和图像坐标的一个或多个相关联；
响应于所述相应图坐标和图像坐标的置信水平，将权重分配给所述图点和图像点中的至少一些；以及
通过基于所述权重计算在所述图坐标和所述图像坐标之间的几何变换来将所述3-D图像与所述图拟合，
其中，所述权重与所述置信水平成比例，使得具有低置信水平的图和图像点被分配较低的权重，而具有高置信水平的图和图像点被分配较高的权重，以及
其中，所述权重的分配进一步考虑到所述器官内相关区域中的组织特性。
2.根据权利要求1所述的方法，其中接受所述图包括使用位置跟踪系统来测量在所述相应图坐标的生理学参数的值。
3.根据权利要求2所述的方法，其中测量在所述相应图坐标的生理学参数的值包括：
在所述器官中的多个位置上移动探头，所述探头包括位置传感器和测量所述生理学参数的传感器；以及
对在所述位置处的所述生理学参数和所述位置传感器的位置执行相应的多次测量。
4.根据权利要求1所述的方法，并且包括使用所述几何变换将所述生理学参数的所述值叠加到所述3-D图像上，并且将所叠加的值和3-D图像显示给用户。
5.根据权利要求1所述的方法，其中将所述3-D图像与所述图拟合包括：
定义在至少一些所述图像点和相应图点之间的距离，
将分配给所述相应图点和图像点的所述权重应用于所述距离以产生加权距离，并且基于所述距离来计算累积的加权距离量度，所述量度指示在将所述3-D图像与所述图拟合中所述几何变换的精度。
6.根据权利要求5所述的方法，其中应用所述权重包括将所述距离的平方乘以相应的权重以产生所述加权距离，以及其中计算所述累积的加权距离量度包括对所述加权距离求和。
7.根据权利要求5所述的方法，其中将所述3-D图像与所述图拟合进一步包括修改所述几何变换以便最小化所述累积的加权距离量度。
8.根据权利要求1所述的方法，其中将所述3-D图像与所述图拟合包括：
识别出现在所述3-D图像中的所述器官外部的解剖标志，
在所述解剖标志的位置处获得相应的图点，并且配准所述3-D图像中所述解剖标志的图像坐标与所获得图点的图坐标。
9.根据权利要求1所述的方法，其中使所述置信水平相关联包括接受来自用户的置信水平。
10.根据权利要求1所述的方法，其中使所述置信水平相关联包括基于所述图点和图像点来自动地确定所述置信水平。
11.根据权利要求10所述的方法，其中自动地确定所述置信水平包括：
检测在所述图坐标和图像坐标中的至少一个中的失真，并且响应于所述失真的水平来分配所述置信水平。
12.根据权利要求10所述的方法，其中自动地确定所述置信水平包括响应于下述中的至少一个来分配所述置信水平：
在所述图中所述图坐标的相应一个或多个的密度；
在所述相应图坐标和图像坐标附近的电组织特征；
在所述相应图坐标和图像坐标附近的生理学组织特征；
所述相应图坐标和图像坐标距在所述图和所述3-D图像两者中所识别的配准标志的距离；以及
所述配准标志的生理学特征。
13.根据权利要求1所述的方法，其中接受所述3-D图像包括接受由所述图像点的至少一些所限定的所述器官的实体模型、表面模型和线框模型中的至少一个。
14.一种用于图像配准的装置，包括：
用于接受器官的三维（3-D）图像的部件，所述图像包括具有相应图像坐标的多个第一图像点；
用于接受定义所述器官上生理学参数的值的分布的图的部件，所述图包括具有相应图坐标的多个第二图点；
用于使相应的置信水平与所述图坐标和图像坐标的一个或多个相关联的部件；
用于响应于所述相应图坐标和图像坐标的置信水平分配权重给所述图点和图像点中的至少一些的部件；以及
用于通过基于所述权重计算在所述图坐标和所述图像坐标之间的几何变换来将所述3-D图像与所述图拟合的部件，
其中，所述权重与所述置信水平成比例，使得具有低置信水平的图和图像点被分配较低的权重，而具有高置信水平的图和图像点被分配较高的权重，以及
其中，用于分配权重的所述部件在执行所述权重分配时进一步考虑到所述器官内相关区域中的组织特性。
15.根据权利要求14所述的装置，其中用于接受所述图的所述部件进一步包括用于使用位置跟踪系统来测量在所述相应图坐标的生理学参数的值的部件。
16.根据权利要求15所述的装置，其中用于使用所述位置跟踪系统来测量所述生理学参数的所述值的所述部件包括用于在所述器官中的多个位置上移动所述位置跟踪系统的探头以便对在所述位置的所述生理学参数和位置传感器的位置执行相应的多次测量的部件，所述探头包括所述位置传感器和测量所述生理学参数的传感器。
17.根据权利要求14所述的装置，进一步包括：用于使用所述几何变换将所述生理学参数的所述值叠加到所述3-D图像上并且将所叠加的值和3-D图像显示给用户的部件。
18.根据权利要求14所述的装置，其中用于将所述3-D图像与所述图拟合的所述部件包括：
用于定义在至少一些所述图像点和相应图点之间的距离的部件，
用于将分配给所述相应图点和图像点的所述权重应用于所述距离以产生加权距离的部件，以及
用于基于所述距离来计算累积的加权距离量度的部件，所述量度指示在将所述3-D图像与所述图拟合中所述几何变换的精度。
19.根据权利要求18所述的装置，其中用于应用所述权重的所述部件包括用于将所述距离的平方乘以相应的权重以产生所述加权距离的部件，以及其中用于计算所述累积的加权距离量度的所述部件包括用于对所述加权距离求和以产生所述累积的加权距离量度的部件。
20.根据权利要求18所述的装置，其中用于将所述3-D图像与所述图拟合的所述部件进一步包括用于通过修改所述几何变换以便最小化所述累积的加权距离量度来将所述3-D图像与所述图拟合的部件。
21.根据权利要求14所述的装置，其中用于将所述3-D图像与所述图拟合的所述部件进一步包括：用于接受出现在所述3-D图像中的所述器官外部的解剖标志的识别以及在所述解剖标志的位置处获得的相应图点的识别，并且配准所述3-D图像中的所述解剖标志的图像坐标与所获得图点的图坐标的部件。
22.根据权利要求14所述的装置，其中用于关联相应的置信水平的所述部件进一步包括用于经由接口接受来自用户的置信水平的部件。
23.根据权利要求14所述的装置，其中用于关联相应的置信水平的所述部件进一步包括用于基于所述图点和图像点来自动地确定所述置信水平的部件。
24.根据权利要求23所述的装置，其中用于自动地确定所述置信水平的所述部件进一步包括：用于检测在所述图点和图像点中的至少一个中的失真并且响应于所述失真的水平来分配所述置信水平的部件。
25.根据权利要求23所述的装置，其中用于自动地确定所述置信水平的所述部件进一步包括用于响应于下述中的至少一个来分配所述置信水平的部件：
在所述图中所述图坐标的相应一个或多个的密度；
在所述相应图坐标和图像坐标附近的电组织特征；
在所述相应图坐标和图像坐标附近的生理学组织特征；
所述相应图坐标和图像坐标距在所述图和所述3-D图像两者中所识别的配准标志的距离；以及
所述配准标志的生理学特征。
26.根据权利要求14所述的装置，其中所述3-D图像包括由所述图像点的至少一些所限定的所述器官的实体模型、表面模型和线框模型中的至少一个。
CN201510971118.3A 2006-03-21 2007-03-21 使用局部加权拟合的图像配准 Active CN105574878B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11/385,060 US7855723B2 (en) 2006-03-21 2006-03-21 Image registration using locally-weighted fitting
US11/385060 2006-03-21
CNA200710101605XA CN101051387A (zh) 2006-03-21 2007-03-21 使用局部加权拟合的图像配准

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CNA200710101605XA Division CN101051387A (zh) 2006-03-21 2007-03-21 使用局部加权拟合的图像配准

Publications (2)

Publication Number Publication Date
CN105574878A CN105574878A (zh) 2016-05-11
CN105574878B true CN105574878B (zh) 2020-04-28

Family

ID=38293155

Family Applications (2)

Application Number Title Priority Date Filing Date
CNA200710101605XA Pending CN101051387A (zh) 2006-03-21 2007-03-21 使用局部加权拟合的图像配准
CN201510971118.3A Active CN105574878B (zh) 2006-03-21 2007-03-21 使用局部加权拟合的图像配准

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CNA200710101605XA Pending CN101051387A (zh) 2006-03-21 2007-03-21 使用局部加权拟合的图像配准

Country Status (11)

Country Link
US (1) US7855723B2 (zh)
EP (1) EP1837828B1 (zh)
JP (1) JP5345290B2 (zh)
KR (1) KR20070095788A (zh)
CN (2) CN101051387A (zh)
AT (1) ATE546798T1 (zh)
AU (1) AU2007201200B2 (zh)
BR (1) BRPI0701778A (zh)
CA (1) CA2582703C (zh)
IL (1) IL182000A (zh)
MX (1) MX2007003312A (zh)

Families Citing this family (65)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7855723B2 (en) * 2006-03-21 2010-12-21 Biosense Webster, Inc. Image registration using locally-weighted fitting
EP2052362B1 (en) * 2006-08-08 2018-10-17 Koninklijke Philips N.V. Registration of electroanatomical mapping points to corresponding image data
US9549689B2 (en) 2007-03-09 2017-01-24 St. Jude Medical, Atrial Fibrillation Division, Inc. System and method for correction of inhomogeneous fields
US10433929B2 (en) * 2007-03-09 2019-10-08 St. Jude Medical, Atrial Fibrillation Division, Inc. System and method for local deformable registration of a catheter navigation system to image data or a model
US7995864B2 (en) * 2007-07-03 2011-08-09 General Electric Company Method and system for performing image registration
US20100316270A1 (en) * 2007-12-20 2010-12-16 Koninklijke Philips Electronics N.V. 3d reconstruction of a body and of a body contour
US8219179B2 (en) * 2008-03-06 2012-07-10 Vida Diagnostics, Inc. Systems and methods for navigation within a branched structure of a body
JP5335280B2 (ja) * 2008-05-13 2013-11-06 キヤノン株式会社 位置合わせ処理装置、位置合わせ方法、プログラム、及び記憶媒体
US8160325B2 (en) * 2008-10-08 2012-04-17 Fujifilm Medical Systems Usa, Inc. Method and system for surgical planning
US8160326B2 (en) * 2008-10-08 2012-04-17 Fujifilm Medical Systems Usa, Inc. Method and system for surgical modeling
US8478379B2 (en) * 2008-11-12 2013-07-02 Biosense Webster, Inc. Probe visualization based on mechanical properties
US8406497B2 (en) * 2009-02-25 2013-03-26 Dongwook YOON Method for population-driven identification of body landmarks
JP2012531996A (ja) * 2009-07-06 2012-12-13 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 生理学的パラメータの視覚化
US8391578B2 (en) * 2009-11-11 2013-03-05 General Electric Company Method and apparatus for automatically registering images
US8311303B2 (en) * 2010-01-12 2012-11-13 Siemens Corporation Method and system for semantics driven image registration
US9129392B2 (en) 2010-02-25 2015-09-08 The Trustees Of The University Of Pennsylvania Automatic quantification of mitral valve dynamics with real-time 3D ultrasound
JP5624350B2 (ja) * 2010-04-02 2014-11-12 株式会社東芝 医用画像処理装置
WO2011127247A2 (en) * 2010-04-07 2011-10-13 Sanjay Krishna Apparatus and techniques of non-invasive analysis
DE102010019421A1 (de) * 2010-05-05 2011-11-10 Siemens Aktiengesellschaft Bildgebendes Verfahren zur Darstellung von Ergebnissen intravaskulärer Bildgebung und CFD-Ergebnissen sowie Medizinsystem zur Durchführung des Verfahrens
US9763587B2 (en) * 2010-06-10 2017-09-19 Biosense Webster (Israel), Ltd. Operator-controlled map point density
US10524765B2 (en) * 2010-12-27 2020-01-07 St. Jude Medical, Atrial Fibrillation Division, Inc. Refinement of an anatomical model using ultrasound
JP6166663B2 (ja) * 2011-03-02 2017-07-19 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． ナビゲーションガイダンスのための視覚化
JP5263995B2 (ja) * 2011-03-18 2013-08-14 富士フイルム株式会社 ネットワーク構築装置および方法ならびにプログラム
US9600856B2 (en) * 2011-06-16 2017-03-21 Koninklijke Philips N.V. Hybrid point-based registration
CN102846337B (zh) * 2011-06-29 2015-09-02 清华大学 三维超声系统及其目标点的定位方法和装置
KR102070427B1 (ko) 2012-08-08 2020-01-28 삼성전자주식회사 종양의 위치를 추적하는 방법 및 장치
CN103857344B (zh) * 2012-09-20 2016-09-28 东芝医疗系统株式会社 图像处理系统、x射线诊断装置以及图像处理方法
JP6310464B2 (ja) 2012-10-05 2018-04-11 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 改善されたｘ線画像を供する医療撮像システム及び方法
JP5987640B2 (ja) * 2012-11-05 2016-09-07 コニカミノルタ株式会社 超音波を用いた被検体の３次元復元方法および装置
US9091628B2 (en) 2012-12-21 2015-07-28 L-3 Communications Security And Detection Systems, Inc. 3D mapping with two orthogonal imaging views
US8995738B2 (en) * 2013-02-18 2015-03-31 Wisconsin Alumni Research Foundation System and method for magnetic resonance imaging parametric mapping using confidence maps
DE102013205278A1 (de) * 2013-03-26 2014-10-02 Siemens Aktiengesellschaft Verfahren zur Darstellung von Signalwerten eines kombinierten Magnetresonanz-Positronenemissionstomographie-Geräts sowie entsprechend ausgestaltetes Magnetresonanz-Positronenemissionstomographie-Gerät
US9576107B2 (en) * 2013-07-09 2017-02-21 Biosense Webster (Israel) Ltd. Model based reconstruction of the heart from sparse samples
US9351657B2 (en) * 2013-07-19 2016-05-31 Biosense Webster (Israel) Ltd. Cardiac activity visualization with frequency discrimination
US9572535B2 (en) 2013-12-05 2017-02-21 Biosense Webster (Israel) Ltd. Dynamic mapping point filtering using a pre-acquired image
WO2016015994A1 (en) * 2014-07-29 2016-02-04 Koninklijke Philips N.V. Ultrasound imaging apparatus
US9754372B2 (en) * 2014-08-15 2017-09-05 Biosense Webster (Israel) Ltd. Marking of fluoroscope field-of-view
EP2989988B1 (en) * 2014-08-29 2017-10-04 Samsung Medison Co., Ltd. Ultrasound image display apparatus and method of displaying ultrasound image
US9955889B2 (en) * 2014-11-03 2018-05-01 Biosense Webster (Israel) Ltd. Registration maps using intra-cardiac signals
US9872628B2 (en) * 2015-07-03 2018-01-23 Elwha Llc Electrocardiogram systems and related methods
CN111887786A (zh) 2015-08-14 2020-11-06 直观外科手术操作公司 用于图像引导外科手术的配准系统和方法
US11202680B2 (en) 2015-08-14 2021-12-21 Intuitive Surgical Operations, Inc. Systems and methods of registration for image-guided surgery
JP6938469B2 (ja) * 2015-09-16 2021-09-22 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 血管を特徴付ける装置
JP6392192B2 (ja) 2015-09-29 2018-09-19 富士フイルム株式会社 画像位置合せ装置、画像位置合せ装置の作動方法およびプログラム
US11413006B2 (en) * 2016-04-26 2022-08-16 Koninklijke Philips N.V. 3D image compounding for ultrasound fetal imaging
EP3454734B1 (en) 2016-05-11 2024-04-17 Affera, Inc. Anatomical model generation
WO2017197294A1 (en) 2016-05-12 2017-11-16 Affera, Inc. Three-dimensional cardiac representation
US10152786B2 (en) * 2016-10-11 2018-12-11 Biosense Webster (Israel) Ltd. Registration of a magnetic tracking system with an imaging device
US10403053B2 (en) 2016-11-15 2019-09-03 Biosense Webster (Israel) Ltd. Marking sparse areas on maps
US11083517B2 (en) 2017-01-19 2021-08-10 Biosense Webster (Israel) Ltd. Enhancing efficiency of repeat ablation by merging current and previous maps
US10314658B2 (en) * 2017-04-21 2019-06-11 Biosense Webster (Israel) Ltd. Registration of an anatomical image with a position-tracking coordinate system based on visual proximity to bone tissue
US10643330B2 (en) * 2017-04-21 2020-05-05 Biosense Webster (Israel) Ltd. Registration of an anatomical image with a position-tracking coordinate system based on proximity to bone tissue
US10682181B2 (en) * 2017-09-06 2020-06-16 Biosense Webster (Israel) Ltd. Methods and systems for modeling and registration of 3-dimensional images of the heart
CN108143489A (zh) * 2017-12-15 2018-06-12 泗洪县正心医疗技术有限公司 一种基于体表静脉特征进行穿刺导航的模具及其制作方法
US10832422B2 (en) * 2018-07-02 2020-11-10 Sony Corporation Alignment system for liver surgery
US20210267694A1 (en) * 2018-07-04 2021-09-02 Navix International Limited Incorporating new location readings to old models
CN113168734A (zh) 2018-11-09 2021-07-23 维得诊断公司 管状结构的切割面显示
US11602280B2 (en) 2019-04-08 2023-03-14 Avent, Inc. In-scale flexible display for medical device position guidance
US11517217B2 (en) * 2019-04-08 2022-12-06 Avent, Inc. In-scale tablet display for medical device position guidance
US11273288B2 (en) 2019-04-08 2022-03-15 Avent, Inc. System and method for medical device position guidance
US10881353B2 (en) * 2019-06-03 2021-01-05 General Electric Company Machine-guided imaging techniques
US11527002B2 (en) 2019-12-05 2022-12-13 Biosense Webster (Israel) Ltd. Registration of an image with a tracking system
EP3866176A1 (en) * 2020-02-17 2021-08-18 Siemens Healthcare GmbH Machine-based risk prediction for peri-procedural myocardial infarction or complication from medical data
US11875459B2 (en) 2020-04-07 2024-01-16 Vida Diagnostics, Inc. Subject specific coordinatization and virtual navigation systems and methods
US11478182B2 (en) 2021-01-07 2022-10-25 Biosense Webster (Israel) Ltd. Incorporating a confidence level into an electrophysiological (EP) map

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6560354B1 (en) * 1999-02-16 2003-05-06 University Of Rochester Apparatus and method for registration of images to physical space using a weighted combination of points and surfaces

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5568384A (en) * 1992-10-13 1996-10-22 Mayo Foundation For Medical Education And Research Biomedical imaging and analysis
US5391199A (en) * 1993-07-20 1995-02-21 Biosense, Inc. Apparatus and method for treating cardiac arrhythmias
ES2210662T3 (es) 1994-08-19 2004-07-01 Biosense, Inc. Sistemas medicos de diagnosis, de tratamiento y de imagen.
US5876336A (en) * 1994-10-11 1999-03-02 Ep Technologies, Inc. Systems and methods for guiding movable electrode elements within multiple-electrode structure
JPH08131403A (ja) * 1994-11-09 1996-05-28 Toshiba Medical Eng Co Ltd 医用画像処理装置
US6690963B2 (en) * 1995-01-24 2004-02-10 Biosense, Inc. System for determining the location and orientation of an invasive medical instrument
US5697377A (en) * 1995-11-22 1997-12-16 Medtronic, Inc. Catheter mapping system and method
IL125761A (en) * 1996-02-15 2005-05-17 Biosense Inc Independently positionable transducers for location system
AU709081B2 (en) * 1996-02-15 1999-08-19 Biosense, Inc. Medical procedures and apparatus using intrabody probes
US5944022A (en) * 1997-04-28 1999-08-31 American Cardiac Ablation Co. Inc. Catheter positioning system
US6239724B1 (en) * 1997-12-30 2001-05-29 Remon Medical Technologies, Ltd. System and method for telemetrically providing intrabody spatial position
US6226542B1 (en) * 1998-07-24 2001-05-01 Biosense, Inc. Three-dimensional reconstruction of intrabody organs
JP2000194863A (ja) * 1998-12-28 2000-07-14 Nippon Telegr & Teleph Corp <Ntt> ３次元構造獲得・復元方法、装置、および３次元構造獲得・復元プログラムを記録した記録媒体
US6650927B1 (en) * 2000-08-18 2003-11-18 Biosense, Inc. Rendering of diagnostic imaging data on a three-dimensional map
US7729742B2 (en) * 2001-12-21 2010-06-01 Biosense, Inc. Wireless position sensor
US20040068178A1 (en) * 2002-09-17 2004-04-08 Assaf Govari High-gradient recursive locating system
JP4109075B2 (ja) * 2002-10-10 2008-06-25 Ｓｒｉスポーツ株式会社 球体の回転特性と飛行特性の測定方法及び球体の回転特性と飛行特性の測定装置
DE10340544B4 (de) * 2003-09-01 2006-08-03 Siemens Ag Vorrichtung zur visuellen Unterstützung einer elektrophysiologischen Katheteranwendung im Herzen
US7409108B2 (en) * 2003-09-22 2008-08-05 Siemens Medical Solutions Usa, Inc. Method and system for hybrid rigid registration of 2D/3D medical images
JP2005189204A (ja) * 2003-12-26 2005-07-14 Fuji Xerox Co Ltd ３次元形状計測装置および方法
US7672705B2 (en) * 2004-07-19 2010-03-02 Resonant Medical, Inc. Weighted surface-to-surface mapping
US20070049817A1 (en) * 2005-08-30 2007-03-01 Assaf Preiss Segmentation and registration of multimodal images using physiological data
US7855723B2 (en) * 2006-03-21 2010-12-21 Biosense Webster, Inc. Image registration using locally-weighted fitting

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6560354B1 (en) * 1999-02-16 2003-05-06 University Of Rochester Apparatus and method for registration of images to physical space using a weighted combination of points and surfaces

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A survey of medical image registration;J. B. Antoine Maintz and Max A. Viergever;《Medical Image Analysis》;19981231;第2卷(第1期);第1页 *

Also Published As

Publication number Publication date
JP5345290B2 (ja) 2013-11-20
AU2007201200B2 (en) 2012-02-09
CN101051387A (zh) 2007-10-10
EP1837828A2 (en) 2007-09-26
IL182000A0 (en) 2007-07-04
AU2007201200A1 (en) 2007-10-11
KR20070095788A (ko) 2007-10-01
IL182000A (en) 2012-02-29
US7855723B2 (en) 2010-12-21
CA2582703C (en) 2012-09-04
CA2582703A1 (en) 2007-09-21
EP1837828B1 (en) 2012-02-22
US20070223794A1 (en) 2007-09-27
MX2007003312A (es) 2008-11-18
EP1837828A3 (en) 2010-07-14
BRPI0701778A (pt) 2007-11-27
ATE546798T1 (de) 2012-03-15
JP2007268259A (ja) 2007-10-18
CN105574878A (zh) 2016-05-11

Similar Documents

Publication Publication Date Title
CN105574878B (zh) 使用局部加权拟合的图像配准
CN111163692B (zh) 依据体内测量重建解剖结构
US20200170623A1 (en) Methods for using radial endobronchial ultrasound probes for three-dimensional reconstruction of images and improved target localization
US7996060B2 (en) Apparatus, method, and computer software product for registration of images of an organ using anatomical features outside the organ
KR100686477B1 (ko) 삼차원적으로 재구성되는 신체내 기관의 벡터 매핑 및 디스플레이 방법
US20110160569A1 (en) system and method for real-time surface and volume mapping of anatomical structures
JP7184626B2 (ja) Ｅｎｔの骨距離が色符号化された顔のマップ
EP3430999A1 (en) Improving impedance-based position tracking performance using scattered interpolant
Timinger et al. Motion compensated coronary interventional navigation by means of diaphragm tracking and elastic motion models
US8948476B2 (en) Determination of cardiac geometry responsive to doppler based imaging of blood flow characteristics
EP3614338B1 (en) Post-mapping automatic identification of pulmonary veins

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant