US 6,546,271 B1Grant

Vascular reconstruction

Issue Date:2003-04-08

•26 Claims

•9 Drawing Sheets

Abstract

A method for vascular reconstruction comprises the steps of advancing a catheter having a position sensor into the vessel, acquiring position information from the sensor at a plurality of points in the vessel, calculating a center-line of the vessel based on said position information, and calculating the inner surface of the vessel. The method of the invention preferably further comprises displaying the vessel reconstruction. The invention is also directed to apparatus for reconstructing vessels which comprises a catheter having a position sensor contained therein, means for acquiring position information from said sensor at a plurality of points in said vessel, means for calculating a centerline of the vessel based on said position information, and means for calculating the inner surface of the vessel.

Metadata

Assignee

Bioscience, Inc.

Inventor

Daniel Reisfeld

Application Information

Application Number:US 09/411,154

Filing Date:1999-10-01

Priority Date:1999-10-01

Art Unit:7

Classifications

IPC:

A61B 505

Field of Search:

600407600462600463600467600459600443600424128916128899

Patent Drawings (9 sheets)

Description

Field of the Invention

[0002] The invention relates generally to methods and apparatus for reconstructing blood vessels, and particularly to methods and apparatus for three-dimensional reconstructions.

Background of the Invention

[0003] The present invention relates to methods for reconstructing blood vessels based on position information acquired from a position sensor contained in a catheter that traverses the vessel.

[0004] Blood vessels are visualized today principally by angiography, in which a contrast medium is injected into the bloodstream and the vasculature is imaged using ionizing radiation imaging modalities. Such imaging must be limited however, due to the adverse effects of cumulative radiation on patients. Furthermore, some patients react adversely to the contrast media used in angiography. Accordingly, it would be desirable to have a visualization method that does not depend on ionizing radiation imaging modalities or on the use of contrast agents.

[0005] U.S. Pat. No. 5,546,951 and U.S. patent application No. 08/793,371, which are incorporated herein in their entirety by reference, disclose methods for sensing an electrical property of the heart tissue, for example, local activation time, as a function of the precise location within the heart. The data are acquired with one or more catheters that are advanced into the heart, the catheters having electrical and location sensors in their distal tips. Methods of creating a map of the electrical activity of the heart based on these data are disclosed in commonly assigned U.S. patent applications No. 09/122,137 and 09/357,559 filed on Jul. 24, 1998 and Jul. 22, 1999, respectively, which are also incorporated herein in their entirety by reference. As indicated in these applications, location and electrical activity is preferably initially measured on about 10 to about 20 points on the interior surface of the heart. These data points are then generally sufficient to generate a preliminary reconstruction or map of the cardiac surface to a satisfactory quality. The preliminary map is formed by defining an initial, closed 3-dimensional curved surface, preferably of an ellipsoid, in a reconstruction space in the volume of the sample points. The closed curve is roughly adjusted to a shape which resembles a reconstruction of the sample points. Thereafter, a flexible matching stage is repeatedly performed one or more times to bring the closed curve to resemble the shape of the actual volume being reconstructed. While the above-described methods provide satisfactory reconstructions of organs, the algorithms employed therein do not provide accurate reconstructions of blood vessels, particularly when the vessels are convoluted or tortuous in shape.

[0006] Reconstruction of blood vessels in images or reconstructions of the heart are beneficial because the blood vessels provide confirmatory landmarks that assist in navigation to specific regions in the heart. Accordingly, it would be desirable to have a method of realistically reconstructing blood vessels that complements existing methods for reconstructing chambers of the heart.

[0007] Cardiologists are increasingly attributing atrial fibrillation to defects in the electrical pathways of the heart that originate in the pulmonary vein. Diagnosis and treatment of this condition requires assessing the electrical activity in the pulmonary vein and subsequently ablating defects in or around the vein. Methods for reconstructing the veins with data that may be accumulated with catheters bearing diagnostic or therapeutic components would facilitate the use and outcomes of these procedures.

[0008] Certain interventional and diagnostic procedures such as catheterization of the brain are preceded by the generation of images of the vasculature by modalities such as magnetic resonance imaging (MRI). A reconstruction of the vessels effected during the catheterization that could be registered with a previously acquired image would allow the physician to identify the catheter tip location during the procedure with respect to the previously acquired image.

Summary of the Invention

[0009] The present invention is directed to a method for reconstructing the inner surface of a blood vessel using a catheter having a position sensor contained therein. The method of the invention comprises the following steps:

[0010] a) advancing the catheter into the vessel;

[0011] b) acquiring position information from the sensor at a plurality of points in the vessel;

[0012] c) calculating a center-line of the vessel based on said position information; and

[0013] d) calculating the inner surface of the vessel.

[0014] In a preferred embodiment, the method of the invention further comprises the step of displaying the reconstructed vessel surface.

[0015] The invention is also directed to apparatus for reconstructing the inner surface of vessels, which comprises:

[0016] a) a catheter having a position sensor contained therein;

[0017] b) means for acquiring position information from said sensor at a plurality of points in said vessel;

[0018] c) means for calculating a centerline of the vessel based on said position information; and

[0019] d) means for calculating the inner surface of the vessel.

[0020] In a preferred embodiment, the apparatus of the invention further comprises means for displaying the vessel reconstruction.

[0021] It is an object of the invention to provide a method and apparatus for reconstructing blood vessels without the use of imaging modalities or contrast agents.

[0022] It is another object of the invention to provide a method of reconstructing blood vessels that is complementary to methods and apparatus for reconstructing the heart.

[0023] It is another object of the invention to provide a method of vessel reconstruction that could provide reconstructions that could be registered with previously acquired images of the vasculature.

[0024] These and other objects, features and advantages of the present invention will be more readily apparent for the detailed description set forth below, taken in conjunction with the accompanying drawings.

Brief Description of the Drawings

[0025] FIG. 1 is a schematic drawing of apparatus used to practice the method of the invention;

[0026] FIG. 2A shows a two-dimensional projection of positions that are acquired in a blood vessel using a catheter having a position sensor contained therein;

[0027] FIG. 2B shows the projection of the positions of FIG. 2A rotated by 90°;

[0028] FIG. 3A shows the projection of the positions of FIG. 2A together with a computed vessel centerline;

[0029] FIG. 3B shows the projection of the positions and centerline of FIG. 3A rotated by 90°;

[0030] FIG. 4A shows a projection of the points and centerline of FIG. 3A with normal vectors drawn at increments perpendicular to the centerline;

[0031] FIG. 4B shows a projection of the points, centerline and vectors of FIG. 4A rotated 90° around the vertical axis;

[0032] FIG. 5A shows a wire-frame reconstruction of the vessel positions depicted in FIG. 2A;

[0033] FIG. 5B shows the wire frame reconstruction of vessel of FIG. 5A rotated by 90°;

[0034] FIG. 5C shows an expanded view of a portion of the wireframe reconstruction of FIG. 5B;

[0035] FIG. 6A shows a shaded reconstruction of the vessel positions depicted in FIG. 2A;

[0036] FIG. 6B shows a shaded reconstruction of the vessel of FIG. 6A rotated by 90°;

[0037] FIG. 7 shows a wire frame reconstruction of the vessel positions of FIG. 2A in which the vessel was reconstructed with a variable radius algorithm;

[0038] FIG. 8 shows a shaded reconstruction of the vessel positions of FIG. 2A in which the vessel was reconstructed with a variable radius algorithm;

[0039] FIG. 9A shows a reconstruction of the points of FIG. 2A and FIG. 2B using an algorithm designed for reconstruction of an organ such as the heart;

[0040] FIG. 9B shows the reconstruction of FIG. 9A rotated by 90°.

Detailed Description of Preferred Embodiments

[0041] The method of the invention for reconstructing the inner surface of a blood vessel using a catheter having a position sensor contained therein comprises the following steps:

[0042] a) advancing the catheter into the vessel;

[0043] b) acquiring position information from the sensor at a plurality of points in the vessel;

[0044] c) calculating a center-line of the vessel based on said position information;

[0045] d) calculating the inner surface of the vessel.

[0046] The position sensor may, for example, be an electromagnetic, magnetic or acoustic sensor. Electromagnetic sensors are preferred for use in the process of the invention. Exemplary sensors, catheters containing said sensors and systems incorporating said sensors and catheters are described, for example in U.S. Pat. Nos. 4,173,228; 4,697,595; 4,821,731; 5,042,486; 5,081,993; 5,391,199; 5,558,091; 5,729,129; 5,752,513; 5,913,820; and 5,954,665; and in PCT Applications WO 96/05768 and WO 97/24983, which are incorporated herein in their entirety by reference.

[0047] The position sensor acts as a receiving antenna that receives signals from field generators external to the body. In an alternative embodiment, the position sensor may act as a generator to transmit fields that are received by antennas external to the body.

[0048] The position sensor is preferably contained in the catheter at or proximate to the catheter distal tip.

[0049] Position information acquired by the position sensor contained in the catheter can be used to provide a reconstruction or a map of the physical dimensions of the vessel or vasculature. If additional functional components are contained in the catheter, condition information describing a condition of the blood vessel, may be collected together with position information. As used herein, the term “condition” refers to either a scalar or a vector quantity, and may comprise, for example, an electrical property, a temperature, a pressure, a pH, a measure of local vessel movement or any other condition or combination thereof. For example, the use of a catheter having an electrophysiology electrode in addition to a position sensor permits the simultaneous collection of position and electrical condition information in the vessel. Exemplary electrical condition information that may be collected includes but is not limited to voltage, impedance, conduction velocity and local activation time (LAT). The combined position and condition information may be used to generate a map of the condition information as a function of the spatial coordinates within the vessel.

[0050] The catheters used in the method of the invention may also contain other functional components for effecting therapeutic functions within the body, such as electrodes for selectively ablating tissue within the vessel or other organ within the body.

[0051] FIG. 1 shows a preferred apparatus for carrying out the method of the invention. The apparatus comprises catheter 21 for insertion into the human body. Distal end 24 of catheter 21 includes a sensor 28 that generates signals used to determine the position, and, optionally, the orientation of the catheter within the body. Sensor 28 is preferably an electromagnetic sensor comprising a plurality of coils as described in PCT application WO96/05768, which is incorporated herein in its entirety by reference. This sensor enables continuous generation of up to six dimensions of position and orientation information with respect to externally applied magnetic fields. Alternatively, sensor 28 may comprise only a single sensing coil as described in U.S. Pat. No. 5,913,820. Sensor 28 may comprise other position and/or coordinate sensors as described in U.S. Pat. No. 5,391,199, U.S. Pat. No. 5,443,489 and PCT application WO94/04938 which are incorporated herein by reference. Further, tip 22 may be coated with an opaque marking material to visualize the tip under an imaging apparatus such as a fluoroscope.

[0052] Distal end 24 of catheter 21 optionally includes a functional portion 23 for performing diagnostic and/or therapeutic functions, adjacent to distal tip 22. Functional portion 23 preferably comprises electrodes or sensors for performing electrophysiological measurements, as described, for example, in U.S. Pat. No. 5,391,199 or in PCT application WO97/24983, which are incorporated herein by reference. Alternatively or additionally, functional portion 23 may include other diagnostic apparatus for acquiring parameter values at points within the body. Functional portion 23 may also include therapeutic apparatus as known in the art.

[0053] Sensor 28 is preferably adjacent to and in a fixed relationship relative to functional portion 23 and to tip 22 of catheter 21.

[0054] Catheter 21 preferably includes a handle 30, having controls 32 which are used to steer distal end 24 of catheter 21 in a desired direction. Catheter 21 preferably comprises a steering mechanism in distal end 24 as is known in the art to facilitate repositioning of tip 22.

[0055] Catheter 21 is coupled via an extension cable 25 to a console 34 which enables the user to observe and regulate the function of catheter 21. Console 34 preferably includes a computer 36, keyboard 38, signal processing circuitry 40, which are typically inside computer 36, and display 42. Signal processing circuits 40 typically receive, amplify, filter and digitize signals from catheter 21, including signals from sensor 28 and functional portion 23, whereupon these digitized signals are used by computer 36 to compute the position and/or orientation of catheter tip 22 and to record condition information measured from functional portion 23. Alternatively, appropriate circuitry may be associated with catheter 21 itself so that circuits 40 receive signals that are already amplified, filtered and/or digitized. Preferably, computer 36 includes a memory for storing position and condition information. Computer 36 also comprises means for capturing images from an imaging modality either using a video or a DICOM protocol interface. Computer 36 preferably further comprises dedicated graphics hardware for rapidly calculating the vessel reconstructions and for displaying them on display 42. Preferably, the computer is equipped to receive body surface ECG signals from ECG monitor 73 which is connected to a plurality of ECG body surface leads 52. Alternatively, ECG monitoring may also be conducted directly by circuits 40. Since blood vessels undergo movement as a function of the cardiac cycle, the use of the body surface ECG permits the gating of the position information and the resultant reconstruction to a single point in the cardiac cycle.

[0056] An algorithm is used in the reconstruction of the vessel from the position information acquired at a plurality of acquisition points by the position sensor. The input to the algorithm is a series of three-dimensional sensor locations P₁, P₂, . . . , P_n, that are acquired within the vessel such that P_i=(x_i, y_i, z_i) wherein x_i, y_iand z_iare the particular coordinates of the i^thacquisition point.

[0057] FIG. 2A shows a two dimensional projection of the positions of points 26 at which position information is acquired using the above-described catheter and sensor. FIG. 2B shows the plot of the acquisition points of FIG. 2A rotated 90° about the vertical axis of FIG. 2A.

[0058] The algorithm used in the method of the invention approximates the inner surface of the vessel such that cross-sections of the vessel perpendicular to the vessel centerline are circular or nearly circular in shape. One of ordinary skill in the art will appreciate that the method of the invention may be applied to generate reconstructions having other cross-sectional shapes, such as an elliptical or polygonal shape, for example.

[0059] The vessel is reconstructed in two steps: first, the centerline of the vessel is calculated. Then the vessel wall of appropriate radius or cross-section is calculated around the centerline.

Calculation of the Vessel Centerline

[0060] We use d_if=∥P_i-P_j∥, the distance between sample point P_iand P_j, to find an extreme point, P_m, (28 in FIG. 2A) which is the sample point that is the most remote from all other sample points, i.e.

m = \underset{i}{\max \arg} \sum_{j} d_{ij}^{2} .

[0061] In other words, for each point, we calculate the distance from that point to all other points. The extreme point, P_m, is that point having the longest distance from itself to the other points. The distance between the points from which to base the selection of P_mis assessed, as in the above equation, as the sum of squares of the distances between points.

[0062] We denote the distance between a sample point, P_i, and the extreme point, P_m, by d_i, i.e., d_i=d_mi.

[0063] Let d =max d_i, the maximum distance between the extreme point and any of the other points. We may then define t_ias:

t_{i} = \frac{d_{i}}{d} \in [0, 1] .

[0064] The vessel centerline is a one-dimensional line embedded in three-dimensional space. The centerline represents the locus of all geometric centers of all possible cross-sections of the vessel.

[0065] We seek a parametric representation of the vessel centerline that is of the form:

[0066] F(t)=(X(t), Y(t), Z(t)) tε[0,1]

[0067] wherein X(t), Y(t) and Z(t) are functions of t that represent the coordinates of the centerline.

[0068] Each location P_ihas a corresponding point F(t_i), the point on the centerline closest to P_i. The representation of the centerline F(t) preferably minimizes the distance between the centerline and the sample points in a least squares sense, i.e., Σ∥(t₁)-P_i∥²is minimal over the function space.

[0069] The following series of polynomials of degree k are used for calculating the centerline representation:

X (t) = \sum_{j = 0}^{k} a_{j} t^{j}

Y (t) = \sum_{j = 0}^{k} b_{j} t^{j}

Z (t) = \sum_{j = 0}^{k} c_{j} t^{j} .

[0070] We wish to use a polynomial of high enough degree to be able to fit the data but of low enough degree to avoid over-fitting the data. In practice, a third degree polynomial fit, i.e., in which k=3, is preferred.

[0071] For each of the above polynomial equations (one for each coordinate dimension), we solve a system of linear equations whose unknowns are the coefficients a_j, b_jand c_j=0, . . . ,k). The system of equations equates the polynomials X(t), Y(t) and Z(t) with the respective coordinates of each of the acquired data points according to the following equations:

[0072] X(t_i)=x_ii=1, . . . ,n

[0073] Y(t_i)=y_ii=1, . . . ,n

[0074] Z(t_i)=z_ii=1, . . . ,n

[0075] Singular value decomposition is a robust, preferred method for solving these systems of equations (see for example, “Numerical Recipes in C: The Art of Scientific Computing”, William T. Vetterling (Editor), Saul A. Teukolsky, William H. Press (Editor), and Brian P. Flannery, Cambridge University Press, pp 59-70, 1997).

[0076] FIG. 3A shows a two-dimensional projection of the original acquisition points as well as the computed centerline 29. FIG. 3B shows the points and centerline of FIG. 3A rotated by 90°.

Calculation of the Vessel Wall About the Centerline

[0077] As indicated above, the algorithm used in the method of the invention approximates the inner surface of the vessel wall such that cross-sections of the vessel perpendicular to the vessel centerline are circular or polygonal in shape. Thus, the reconstruction models the vessel as a whole as being tubular in shape.

[0078] A three-dimensional reconstruction about the centerline may be performed to generate a tube of either fixed or variable radius (fixed or variable cross-section).

[0079] Reconstruction of the vessel as a tube is performed by calculating circular or polygonal slices or cross-sections around, and preferably perpendicular to the centerline, and connecting the slices together to form a tube.

[0080] The tangent to the centerline may be expressed by the following series of equations:

D (t) = \frac{ⅆ}{ⅆ t} F (t) = (\frac{ⅆ}{ⅆ t} X (t), \frac{ⅆ}{ⅆ t} Y (t), \frac{ⅆ}{ⅆ t} Z (t))

[0081] where

\frac{ⅆ}{ⅆ t} X (t) = \sum_{j = 1}^{k} a_{j} {jt}^{j - 1}

\frac{ⅆ}{ⅆ t} Y (t) = \sum_{j = 1}^{k} b_{j} {jt}^{j - 1}

\frac{ⅆ}{ⅆ t} Z (t) = \sum_{j = 1}^{k} c_{j} {jt}^{j - 1}

[0082] Denoting

η (V) = \frac{V}{‖ V ‖},

[0083] one of the unit vectors normal to the centerline may be expressed as:

N_{0} (t) = η (- \frac{ⅆ}{ⅆ t} Y (t), \frac{ⅆ}{ⅆ t} X (t), 0)

[0084] FIG. 4A shows the points and centerline of FIG. 3A with normal vectors 41 drawn at increments perpendicular to the centerline. FIG. 4B shows the points, centerline and vectors of FIG. 4A rotated 90° around the vertical axis.

[0085] Another unit vector normal to the centerline that is perpendicular to the previous vector may be expressed by the equation:

[0086] N₁(t)=D(t)×N₀(t)

[0087] Two additional unit vectors are −N₀(t) and −N₁(t). Thus, a set of four vectors N⁰(t) normal to centerline F(t) are ordered counterclockwise in 90° increments about the centerline as follows:

[0088] N⁰(t)=(N₀(t),N_t(t),−N₀(t),−N_t(t))

[0089] N⁰(t) is a very crude sampling of the vectors emanating from the centerline and projecting onto a circle around the centerline F(t). Given a set of vectors Nⁱ(t) that sample the circle around the centerline, additional vectors are in the direction of the sum of previously calculated vectors. Thus, given a set of four vectors, N⁰(t), the set may be expanded to a new set, N¹(t) containing 8 vectors as shown below:

[0090] N_2jⁱ⁺¹(t)=N_jⁱ(t) N_2j+1ⁱ⁺¹(t)=η(N_jⁱ(t)+N_j+1ⁱ(t))

[0091] Another set of vectors, N²(t), containing 16 vectors, may be similarly generated from N¹(t), and so on.

[0092] Finally, let r be the radius of the tube about any point on the centerline. For a tube of fixed radius, a point on the surface of the tube around centerline F(t) which corresponds to the vector N₁(t) is

[0093] S_i(t)=F(t)+r N_i(t)

[0094] The radius r of the tube may be chosen based on the position information of the points and their distances from the centerline. For example, radius r may be chosen as the average or the median of the distances of the points from the centerline. Alternatively, the reconstruction may be performed using a value of r selected by the user that may be typical of the dimensions of the vessel under reconstruction.

[0095] A wire-frame reconstruction of the tubular-shaped vessel is built from small rectangular patches whose vertices are S_i(t), S_i(t+Δ), S_i+1, (t+Δ), S_i+, (t), i.e., points corresponding to two contiguous vectors in a slice, each of said points connected to corresponding points on the next adjacent slice or cross-section of the tube. FIG. 5A shows the wire frame reconstruction of the points, centerline and first vectors of FIG. 4A. FIG. 5B shows the wire-frame reconstruction of FIG. 5A rotated 90° around its vertical axis. FIG. 5C shows a greatly magnified view of a portion of the wire frame reconstruction of FIG. 5B. In the reconstruction as seen in FIG. 5C, each slice of the reconstruction is composed of a 16-sided polygon. It will be appreciated that as the number of sides in the polygon is increased, the vessel reconstruction will appear more circular in cross-section. Each vertex of the polygon represents the radius (not shown) emanating from the centerline (not shown) in the direction of one of the aforesaid unit vectors. The sides of the polygon are made by connecting contiguous vertices within a slice with straight lines. The squares connecting the slices are made by connecting the corresponding vertices on adjacent slices with straight lines.

[0096] The reconstruction is completed by shading the squares making up the wire-frame reconstruction. FIG. 6A and FIG. 6B show the shaded reconstructions of FIG. 5A and FIG. 5B, respectively. Using standard graphic methods known in the art (see for example “OpenGL(r) 1.2 Programming Guide, Third Edition: The Official Guide to Learning OpenGL, Version 1.2”, Mason Woo, et. al., Addison-Wesley Publishing Company, New York, N.Y., 1999), the individual rectangles comprising the wire-frame reconstruction may be shaded using a gray scale or a color scale to lend perspective to the otherwise flat two-dimensional representation of the three dimensional vessel structure. Alternatively, if position information is collected together with condition information using a catheter having both a position and condition sensor, the wire frame reconstruction may be shaded or colored such that the colors or shades of the individual squares represent different values of the condition information as a function of vessel coordinates.

[0097] As indicated hereinabove, the method of the invention may be used to provide a reconstruction of the vessel having either a fixed or a variable cross-section. Reconstruction of a tube with a variable radius or cross-section requires a modified slice formula:

[0098] S_i(t)=F(t)+r(t)N_i(t)

[0099] wherein r(t), the slice radius around centerline F(t), is itself a function of (t).

[0100] As indicated above, each sample point, P_i, has a corresponding point on the centerline, F(t_i), to which it is closest. The distance from P_ito its corresponding point F(t_i) on the centerline, ∥P_i−F(t_i)∥, is an indication of the vessel radius at F(t_i). The radius of a slice may be determined by averaging the distances between sample points and their corresponding centerline points in vicinity of the slice center. For each slice, points closer to the slice may be accorded greater weight in the computation of the average. In other words, the radius of a slice may be computed as a weighted average of the distance between the points to the centerline, with greater weight being accorded to points closer to the slice.

[0101] FIG. 7 and FIG. 8 show the wire-frame and shaded reconstructions, respectively, of the position information of FIG. 2A using the above-described variable radius reconstruction algorithm.

[0102] The decision to reconstruct the vessel with a fixed or variable radius depends on a number of factors. A high quality reconstruction of a tube with variable radius requires more data points taken around the vessel cross-section. This requirement translates into longer procedure time to acquire the points. Alternatively, a fixed radius tube may be reconstructed from fewer data points around the vessel cross-section, which may be acquired, for example by merely traversing the catheter through the vessel. A fixed radius reconstruction would be silent on variations in cross-section of the vessel, but would, nevertheless, be expected to accurately portray the three-dimensional vessel shape.

[0103] As previously indicated, copending commonly assigned applications 09/122,137 and 09/357,559 disclose methods of mapping the electrical activity of the heart. The reconstruction disclosed in these applications is effected by generating a preliminary map by defining an initial, closed 3-dimensional curved surface, preferably of an ellipsoid, in a reconstruction space in the volume of the sample points. The closed curve is roughly adjusted to a shape which resembles a reconstruction of the sample points. Thereafter, a flexible matching stage is repeatedly performed one or more times to bring the closed curve to resemble the shape of the actual volume being reconstructed. The data of FIG. 2A and FIG. 2B were reconstructed using the algorithm of these copending applications, and shaded reconstructions of the data using those methods are shown in FIG. 9A and FIG. 9B. It will be appreciated that the method of the invention successfully in reconstructs vessels more realistically than the above-disclosed methods designed for heart reconstructions.

[0104] Although this invention has been described in connection with its most preferred embodiments, it will become readily apparent to those reviewing this detailed specification that numerous additional embodiments fall well within the scope and spirit of the claimed invention as set forth in the claims which appear below.

Claims

What is claimed is:

1. A method of reconstructing the inner surface of a blood vessel using a catheter having a position sensor contained therein, said sensor generating signals responsive to the position and/or the orientation of the catheter within the body, said method comprising:

a) advancing the catheter into the vessel;

b) acquiring signals from the sensor and calculating position information from said acquired signals at a plurality of acquisition points in the vessel;

c) calculating a centerline of the vessel based on said position information;

d) calculating the inner surface of the vessel.

2. The method of claim 1 wherein the centerline is described as a parametric function.

3. The method of claim 2 wherein the parametric function is of a polynomial form.

4. The method of claim 1 wherein the distances between acquisition points and respective points on the centerline closest to said acquisition points are minimized.

5. The method of claim 4 wherein said distances are minimized in a least squares sense.

6. The method of claim 1 wherein said inner vessel surface is reconstructed in approximately circular cross-section.

7. The method of claim 1 wherein said reconstruction comprises calculating slices of the vessel about said centerline.

8. The method of claim 7 wherein said slices are perpendicular to said centerline.

9. The method of claim 8 wherein said slices perpendicular to said centerline have a fixed cross-section.

10. The method of claim 8 wherein said slices perpendicular to said centerline have a variable cross-section.

11. The method of claim 1 wherein the distance from each point on the reconstruction to the centerline is a function of the distances of the acquisition points to the centerline.

12. The method of claim 11 wherein the distance from each point on the reconstruction to the centerline is the average distance of the acquisition points to the centerline.

13. The method of claim 11 wherein the distance from each point on the reconstruction to the centerline is the median distance of the acquisition points to the centerline.

14. The method of claim 11 wherein said function is a weighted average of the distances of the acquisition points to the centerline.

15. The method of claim 14 wherein said weighted average accords greater weight to proximate acquisition points.

16. The method of claim 11 wherein said distance is user selected.

17. The method of claim 1 which further comprises acquiring condition information at said acquisition points.

18. The method of claim 17 wherein said reconstruction is color-coded to depict values of said condition information.

19. The method of claim 17 wherein values of condition information on the vessel surface intermediate said acquisition points are interpolated.

20. The method of claim 1 which further comprises the step of displaying the reconstructed vessel surface.

21. The method of claim 1 wherein said position information is obtained by traversing the catheter in the vessel.

22. The method of claim 1 wherein said position information is obtained by sampling points on the vessel wall around the vessel cross-section.

23. The method of claim 1 wherein said position sensor is an electromagnetic sensor.

24. Apparatus for reconstructing the inner surface of a blood vessel comprising:

a) catheter having a position sensor contained therein, said sensor generating signals responsive to the position and/or the orientation of the catheter within the body;

b) means for acquiring signals from said sensor at a plurality of points in said vessel and for calculating position information from said acquired signals;

c) means for calculating a centerline of the vessel based on said position information; and

d) means for calculating the inner surface of the vessel.

25. The apparatus of claim 24 which further comprises means for displaying the reconstruction.

26. The apparatus of claim 24 wherein said sensor is an electromagnetic sensor.

Patent Citations (104)

Patent	Date	Inventor	Cited By
US3316896(A)	1967-05-01	Thomasset	Applicant
US4157572(A)	1979-06-01	Kennedy et al.	Applicant
US4173228(A)	1979-11-01	Van Steenwyk et al.	Applicant
US4459990(A)	1984-07-01	Barnea	Applicant
US4522212(A)	1985-06-01	Gelinas et al.	Applicant
US4628937(A)	1986-12-01	Hess et al.	Applicant
US4630203(A)	1986-12-01	Szirtes	Applicant
US4660571(A)	1987-04-01	Hess et al.	Applicant
US4682603(A)	1987-07-01	Franz	Applicant
US4697595(A)	1987-10-01	Breyer et al.	Applicant
US4699147(A)	1987-10-01	Chilson et al.	Applicant
US4821731(A)	1989-04-01	Martinelli et al.	Examiner
US4875165(A)	1989-10-01	Fencil et al.	Applicant
US4878115(A)	1989-10-01	Elion	Applicant
US4892104(A)	1990-01-01	Ito et al.	Applicant
US4898181(A)	1990-02-01	Kessier	Applicant
US4905705(A)	1990-03-01	Kizakevich et al.	Applicant
US4911174(A)	1990-03-01	Pederson et al.	Applicant
US4922912(A)	1990-05-01	Watanabe	Applicant
US4940064(A)	1990-07-01	Desai	Applicant
US4955382(A)	1990-09-01	Franz et al.	Applicant
US4962767(A)	1990-10-01	Brownlee	Applicant
US4979510(A)	1990-12-01	Franz et al.	Applicant
US5022396(A)	1991-06-01	Watanabe	Applicant
US5038791(A)	1991-08-01	Collins et al.	Applicant
US5042486(A)	1991-08-01	Pfeiler et al.	Applicant
US5056524(A)	1991-10-01	Oe	Applicant
US5127403(A)	1992-07-01	Brownlee	Applicant
US5146926(A)	1992-09-01	Cohen	Applicant
US5156151(A)	1992-10-01	Imran	Applicant
US5175773(A)	1992-12-01	Garreau et al.	Applicant
US5215103(A)	1993-06-01	Desai	Applicant
US5227969(A)	1993-07-01	Waggener et al.	Applicant
US5228442(A)	1993-07-01	Imran	Applicant
US5231995(A)	1993-08-01	Desai	Applicant
US5239999(A)	1993-08-01	Imran	Applicant
US5243981(A)	1993-09-01	Hudrlik	Applicant
US5255678(A)	1993-10-01	Deslauriers et al.	Applicant
US5255679(A)	1993-10-01	Imran	Applicant
US5279299(A)	1994-01-01	Imran	Applicant
US5293869(A)	1994-03-01	Edwards et al.	Applicant
US5297549(A)	1994-03-01	Beatty et al.	Applicant
US5311866(A)	1994-05-01	Kagan et al.	Applicant
US5313943(A)	1994-05-01	Houser et al.	Applicant
US5324284(A)	1994-06-01	Imran	Applicant
US5341807(A)	1994-08-01	Nardella	Applicant
US5345936(A)	1994-09-01	Pomeranz et al.	Applicant
US5383917(A)	1995-01-01	Desai et al.	Applicant
US5391199(A)	1995-02-01	Ben-Haim	Applicant
US5409000(A)	1995-04-01	Imran	Applicant
US5433198(A)	1995-07-01	Desai	Applicant
US5443489(A)	1995-08-01	Ben-Haim	Applicant
US5454370(A)	1995-10-01	Avitall	Applicant
US5458116(A)	1995-10-01	Egler	Applicant
US5485849(A)	1996-01-01	Panescu et al.	Applicant
US5487391(A)	1996-01-01	Panescu	Applicant
US5531227(A)	1996-07-01	Schneider	Applicant
US5546951(A)	1996-08-01	Ben-Haim	Applicant
US5549109(A)	1996-08-01	Samson et al.	Applicant
US5558191(A)	1996-09-01	Lai	Applicant
US5588432(A)	1996-12-01	Crowley	Applicant
US5595183(A)	1997-01-01	Swanson et al.	Applicant
US5637090(A)	1997-06-01	McGee et al.	Applicant
US5640967(A)	1997-06-01	Fine et al.	Applicant
US5657755(A)	1997-08-01	Desai	Applicant
US5687737(A)	1997-11-01	Branham et al.	Applicant
US5697377(A)	1997-12-01	Wittkampf	Applicant
US5718241(A)	1998-02-01	Ben-Haim et al.	Applicant
US5729129(A)	1998-03-01	Acker	Applicant
US5730704(A)	1998-03-01	Avitall	Applicant
US5738096(A)	1998-04-01	Ben-Haim	Applicant
US5752513(A)	1998-05-01	Acker et al.	Applicant
US5755664(A)	1998-05-01	Rubenstein	Applicant
US5771895(A)	1998-06-01	Slager	Examiner
US5782773(A)	1998-07-01	Kuo et al.	Applicant
US5797849(A)	1998-08-01	Vesely et al.	Applicant
US5803084(A)	1998-09-01	Olson	Applicant
US5820568(A)	1998-10-01	Willis	Applicant
US5830145(A)	1998-11-01	Tenhoff	Examiner
US5830150(A)	1998-11-01	Palmer et al.	Applicant
US5840031(A)	1998-11-01	Crowley	Applicant
US5842984(A)	1998-12-01	Avitall	Applicant
US5879305(A)	1999-03-01	Yock et al.	Examiner
US5889524(A)	1999-03-01	Sheehan et al.	Applicant
US5913820(A)	1999-06-01	Bladen et al.	Examiner
US5921924(A)	1999-07-01	Avitall	Applicant
US5931835(A)	1999-08-01	Mackey	Applicant
US5931863(A)	1999-08-01	Griffin, III et al.	Applicant
US5951485(A)	1999-09-01	Cyrus et al.	Applicant
US5954665(A)	1999-09-01	Ben-Haim	Applicant
US5999587(A)	1999-12-01	Ning et al.	Applicant
US6047080(A)	2000-04-01	Chen et al.	Applicant
US6052618(A)	2000-04-01	Dahlke et al.	Applicant
US6066094(A)	2000-05-01	Ben-Haim	Applicant
US6083162(A)	2000-07-01	Vining	Examiner
EP501993	1997-06-01		Applicant
EP974936	2000-01-01		Applicant
WO94/04938	1994-03-01		Applicant
WO96/05768	1996-02-01		Applicant
WO97/24983	1997-06-01		Applicant
WO98/12663	1998-03-01		Applicant
WO98/35720	1998-08-01		Applicant
WO98/35720	1998-08-01		Applicant
WO99/05971	1999-02-01		Applicant

Non-Patent Literature (13)

Castleman, K.R.; Digital Image Processing (1996); “Curve and Surface Fitting”; pp. 501-507.Applicant
Jain, A.K.; Fundamentals of Digital Image Processing (1989); “The Back-Projection Operator”; pp. 445.Applicant
Foley J.D., van Dam A., Feiner S.K., Hughes J.F.; 2nd Edition in C Computer Graphics Principles and Practice (1996); “Filling Algorithms”; pp. 979-986.Applicant
Gerstenfeld E., Sahakin A., Swiryn S.; Evidence for Transient Linking of Atrial Excitation During Atrial Fibrillation in Humans (1992); Circulation vol. 86, No. 2, pp 375-382.Applicant
Gerstenfeld E., Sahakian A., Baerman J., Ropella K., Swiryn S.; Detection of Changes in Atrial Endocardial Activation With Use of an Orthogonal Catheter (1991); JACC vol. 18, No. 4, pp 1034-1042.Applicant
Kadish A., Spear J., Levine J., Hanich R., Prood C., Moore E.; Vector Mapping of Myocardial Activation (1986); Laboratory Investigation Arrhythmia vol. 74, No. 3, pp 603-615.Applicant
Kass et al.; Proceedings of First International Conference Vision (1987); Snakes: Active Contour Models; pp 259-268.Applicant
Terzopoulos D.; Transactions on Pattern Analysis and Machine Intelligence (1986), vol. PAMI-8, No. 4; Regularization of Inverse Visual Problems Involving Discontinuities; pp 413-424.Applicant
Lai et al.; IEEE Transactions on Pattern Analysis and Machine Intelligence (1995) vol. 17, No. 11; Deformable Contours: Modeling and Extraction; pp 1084-1090.Applicant
Onnasch et al.; Computers in Cardiology, Long Beach, CA, IEEE Computer Society (1975); A Versatile Program for the Documentation and Comparison of Traced Heart Contours; pp 257-262.Applicant
Duda et al.; Communications of the ACM (1972) vol. 15, No. 1; Use of the Hough Transformation to Detect Lines and Curves in Pictures; pp 11-15.Applicant
Alperin, N. et al.; Automated Analysis of Coronary Lesions From Cineangiograms Using Vessel Tracking and Iterative Deconvolution Techniques; Proceedings. Computers in Cardiology (Cat. No. 89CH2932-2), Jerusalem, Israel; 9/89; pp. 153-156, XP-002156471.Applicant
Kitamura, K. et al.; Estimating the 3-D Skeletons and Transverse Areas of Coronary Arteries from Biplane Angiograms; IEEE Transactions on Medical Imaging, vol. 7, No. 3, 9/88; pp. 173-187; XP-002156470.Applicant

Abstract