Method and system for navigating a catheter probe

A system for navigating a catheter probe through a body cavity includes a sensing coil affixed to a distal end of the probe. Magnetic fields are projected into the body cavity to induce voltage signals in the sensing coil that are sufficient to describe the orientation and position of the probe. A set of magnetic coils each generates a substantially uniform field in a single respective dimension. The orientation angles of the sensing coil may be determined from known values of the unidirectional fields and the measured induced voltage signals. Gradient magnetic fields with components in two dimensions are projected into the body cavity to induce another group of voltage signals. The geometrical intersection of constant voltage surfaces developed by certain gradient fields that produce the measured induced voltage signals is a set of lines on which the catheter is located. The point of intersection of such lines yields the positional coordinates.

Skip to: Description  ·  Claims  ·  References Cited  · Patent History  ·  Patent History
Description
CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a reissue of U.S. Pat. No. 5,592,939 issued on Jan. 14, 1997. The disclosures of the above patent is incorporated herein by reference.

U.S. patent application Ser. No. 09/494,213 filed on Jan. 24, 2000 is also a reissue of U.S. Pat. No. 5,592,939 issued on Jan. 14, 1997 and claims benefit under 35 U.S.C. § 120 as a continuation of this U.S. patent application Ser. No. 09/231,854 filed on Jan. 14, 1999.

FIELD OF THE INVENTION

The present invention relates to catheter navigation systems and, more particularly, to a method and system for determining the position and orientation of a catheter probe being used during a surgical procedure.

BACKGROUND OF THE INVENTION

Various configurations have been proposed to guide and detect a catheter probe through the internal spaces of a patient undergoing a surgical procedure. These proposed configurations are characterized by several alternative approaches including, inter alia, procedures for solving equations to determine unknown location parameters, the generation and detection of magnetic fields, and the use of sensing devices affixed to the catheter probe.

U.S. Pat. No. 4,905,698 to Strohl, Jr. et al. discloses a locator device external to a subject for generating an electromagnetic field that projects into the subject. A catheter inserted into the subject is fitted with a sensing coil at its distal end. The phase of the voltage that is induced in the coil in response to the field is compared to the phase of the generated field. When an in-phase condition occurs, this is an indication that the locator is behind the coil; alternatively, an out-of-phase condition indicates that the locator is beyond the coil. Positions intermediate these two rough approximations of the coil position are not determined other than by a beeping indicator that signifies that this intermediate positioning has been reached.

U.S. Pat. No. 4,821,731 to Martinelli et al. discloses an electroacoustical transducer means secured to the distal end of a catheter that is inserted into a subject for generating acoustical pulses that propagate along an imaging axis and reflect from an anatomical area of interest. The acoustic echoes are converted by the transducer means into electrical signals representative of an image of the anatomical area under reflection and the relative position of the transducer means and angular orientation of the sensing/imaging axis.

U.S. Pat. No. 4,642,786 to Hansen discloses a magnetic position and orientation measurement system that determines the location of an object in space with various configurations, each characterized by the attachment of a retransmitter to the object consisting of passive resonant circuits. The retransmitter is in a predetermined position and orientation with respect to the object. A magnetic field is generated at a resonant frequency of the retransmitter which then retransmits a magnetic field for subsequent reception. The position and orientation of the object may be calculated based upon the induced signals as developed by the reception of the retransmitted magnetic field. The original transmission and reception may be implemented with an integrated transceiver, separate transmitter and receiver elements, or a single transmitter and an array of receiver coils.

U.S. Pat. No. 4,317,078 to Weed et al. discloses how the location of a magnetically sensitive element may be determined by moving a magnetic field source along specified reference axes to induce signals in the sensor so as to identify a set of null points representative of certain flux linkage values. The null point locations are used to calculate the sensor position.

U.S. Pat. No. 3,868,565 describes a system where a magnetic field is generated which rotates about a known pointing vector. The generated field is sensed along at least two axes by a sensor attached to the object to be located or tracked. Based upon the relationship between the sensed magnetic field components, the position of the object relative to the pointing vector can be computed.

U.S. Pat. No. 4,173,228 to Van Steenwyk et al. discloses a catheter locating system that includes a sensor attached to the distal end of the catheter. An electromagnetic field is projected into the body cavity with magnetic probe coils. The field is detected by the sensor, which generates an induced signal whose magnitude and phase are representative of field strength, separation of sensor and probe coils, and relative orientation of sensor and probe coils. The probe coil undergoes linear and rotational movement to identify orientations and locations of the probe coil where minima and maxima occur in the measured signal induced in the sensor. This information is representative of the position and orientation of the sensor.

U.S. Pat. No. 5,211,165 to Dumoulin et al. discloses a modified catheter device that includes a small RF transmit coil attached to its distal end. The transmit coil is driven by an RF source to create an electromagnetic field that induces electrical signals in an array of receive coils distributed around a region of interest. Alternatively, the receive coils can be placed on the invasive device and the transmit coils are distributed outside the patient. A minimum of one transmit coil and three receive coils is necessary to precisely determine the location of the invasive device. A series of equations is developed to solve for the unknowns x-y-z-φ-θ.

PCT Application No. WO94/04938 to Bladen et al. describes how the location and orientation of a single sensing coil may be determined from induced signals developed in response to a sequence of applied magnetic fields emanating from three groups of field generators each including three mutually orthogonal coils.

The positioning methodology developed by Bladen et al. involves calculating the distance from the sensing coil to each group of field generators as a function of the induced voltage developed in the sensing coil by the field generator. The distance calculation is used to define the radius of a sphere centered on the respective field generator. The intersection (i.e., overlap) of the spheres is used to calculate an estimate of the sensor position, using the spherical radius extending from the known location of the field generators as the estimate for each generator.

The orientation algorithm of Bladen et al. develops general equations for induced voltage including the entire set of unknown variables (x-y-z location and φ-θ orientation). The algorithm specifically solves for the orientation parameters by substituting the measured induced voltage and the computed x-y-z coordinates into the general induced voltage equation, and then reduces the equations to the unknown variables φ-θ.

In an alternative orientation algorithm described by Bladen et al., the induced voltage is treated as a vector quantity, allowing the angle between the magnetic field at the generator and the radial vector joining the sensor to the generator to be calculated with a dot product computation. The angle between the radial vector and the sensor axis can be determined from the computed field angle using the dipole equations that define the generator fields. This sensor angle and the radial position as determined by the position algorithm together define the sensor position for use in the alternative orientation algorithm. These values are used to compute the angular orientation φ and θ.

OBJECTS OF THE INVENTION

It is a general object of the present invention to obviate the above-noted and other disadvantages of the prior art.

It is a more specific object of the present invention to provide a catheter navigation system capable of determining the location of a catheter probe.

It is a further object of the present invention to develop a catheter navigation system employing a sensing coil affixed to the end of a catheter probe for generating induced voltage signals that are sufficient to describe the position and orientation of the sensing coil.

It is a further object of the present invention to develop a methodology for generating magnetic fields that are sufficient to create a series of soluble mathematical expressions describing the position and orientation of the sensing coil.

SUMMARY OF THE INVENTION

In one aspect of the present invention, an improved method of determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain comprises the steps of:

    • inducing within said sensing coil a set of orientation signal values each representative of an orientation of said sensing coil and independent of a position of said sensing coil;
    • determining the orientation of said sensing coil using said induced orientation signal values;
    • inducing within said sensing coil a set of positional signal values each representative of the position of said sensing coil; and
    • determining the position of said sensing coil using said positional signal values and said determined orientation.

In another aspect of the present invention, an improved system for determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain comprises:

    • first transmit means for projecting into said navigational domain magnetic energy that is sufficient to induce signal values within said sensing coil representative of an orientation of said sensing coil and independent of the position of said sensing coil;
    • second transmit means for projecting into said navigational domain magnetic energy that is sufficient to induce signal values within said sensing coil representative of the position of said sensing coil; and
    • analysis means, coupled to said first transmit means and said second transmit means, for determining the position and orientation of said sensing coil from said induced signal values.

In another aspect of the present invention, an improved method of determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain comprises the steps of:

    • defining the location of said sensing coil with a set of independent location parameters; and
    • sequentially generating within said navigational domain a sequence of magnetic fields for inducing within said sensing coil a corresponding sequence of induced signals each defined by an induced signal expression that functionally relates said induced signal to certain ones of said location parameters, such that said set of location parameters is determinable by sequentially solving individual signal expression groups each including certain ones of said induced signal expressions and sufficient to represent a subset of said location parameters.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 schematically illustrates a perspective view of a patient-supporting examination deck in accordance with a preferred embodiment of the present invention;

FIGS. 2A-C schematically illustrate a series of magnetic coil sets for generating uniform fields in the x-, y-, and z-directions, respectively, in accordance with a preferred embodiment of the present invention, and which are configured within the deck of FIG. 1;

FIGS. 3 schematically illustrates a magnetic coil assembly for determining the positional coordinates of the sensing coil in accordance with a preferred embodiment of the present invention, and which is configured within the deck of FIG. 1;

FIG. 4 is a flow diagram describing the location algorithm in accordance with the present invention;

FIG. 5 schematically depicts the magnetic coil assembly of FIG. 3 to illustrate representative field patterns generated during an excitation period;

FIG. 6 is a trace representatively illustrating surfaces of constant signal from the sensing coil, as generated by the magnetic assembly of FIG. 3;

FIG. 7 shows an upper plan schematic view of the magnetic assembly of FIG. 3;

FIG. 8 schematically illustrates a perspective view of a patient-supporting examination deck in accordance with another embodiment of the present invention; and

FIGS. 9A-D schematically illustrate a series of magnetic coil assemblies configured in the deck and rails of FIG. 8 for determining the orientation and position of the sensing coil in accordance with another embodiment of the present invention.

Throughout the drawings the same or similar elements are identified by the same reference numeral.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The present invention is directed to a method and system for determining the location of a catheter or endoscopic probe inserted into a selected body cavity of a patient undergoing a surgical procedure. The location data is obtained from electrical measurements of voltage signals that are induced within a sensing coil affixed to the distal end of the catheter probe. These induced voltage signals are generated by the sensing coil in response to prespecified electromagnetic fields that project into the anatomical region of interest which contains all prospective locations of the catheters probe. The electrical measurements of the induced signals provide sufficient information to compute the angular orientation and the positional coordinates of the sensing coil, and hence the catheter probe, which collectively define the location of the sensing coil. The present invention is operative as the patient is disposed on a patient-supporting examination deck.

As used herein, “sensing coil” refers to an electrically conductive, magnetically sensitive element that is responsive to time-dependent magnetic fields for generating induced voltage signals as a function of and representative of the applied time-dependent magnetic field. The sensing coil is adaptable for secure engagement to the distal end of a catheter probe.

As used herein, “navigational domain” refers to a fully enclosed spatial region whose internal volume substantially encloses the complete prospective range of movement of the sensing coil. The navigational domain may be defined by any geometrical space but preferably takes the form of a spherical volume. Under surgical operating conditions, the navigational domain will correspond to an anatomical region of the recumbent patient where surgical viewing or investigation is desired (e.g., a diseased area of tissue or an organ).

As used herein, “last navigational point” (hereinafter “the LNP”) refers to the most recently determined location of the sensing coil before another iteration of the location algorithm is performed.

As used herein, “uniform field” refers to a magnetic field having a large magnetic field component in a specified axial dimension and relatively smaller magnetic field components in the other axial dimensions, and characterized by substantially uniform field values throughout the navigational domain. In the x-y-z coordinate system used herein, where the uniform fields of interest are the x-directed, y-directed, and z-directed fields, the induced voltage signals developed by such fields in the sensing coil are designated Vx, Vy and Vz, respectively. The term “unidirectional field” is used interchangeably with “uniform field” when appropriate.

As used herein, “unidirectional coils” refer to a magnetic assembly that is operative to generate a uniform field (as defined above) within the navigational domain. A distinct magnetic assembly is employed for each uniform field. Although the unidirectional coils described herein are preferably implemented with a collection of appropriately designed magnetic coils, this implementation should not be construed as a limitation of the present invention. Rather, the unidirectional coils may be constructed from any magnetic configuration that is sufficient to generate the uniform fields.

As used herein, “gradient field” refers to a time-dependent magnetic field having non-zero field components (i.e., components with a high spatial gradient) in two of the three axial dimensions for the coordinate system of interest (e.g., x-y-z system), and a substantially zero component in the remaining axial dimension. For mathematical purposes, a substantially zero component is generated when its value is small compared to the net vector resulting from the other two field components.

As used herein, “constant signal surface” or “constant voltage surface” refers to a surface contour along which at every possible point of location for the sensing coil the same induced voltage is developed in the sensing coil.

As used herein, “delta coils” refers to a magnetic assembly for generating a gradient field (as defined above) within the navigational domain. As will become more apparent hereinafter, the delta coils will typically be described in the context of delta coil pairs including a long coil set and a short coil set each generating gradient fields with components in the same axial dimensions but whose magnetic field patterns are different. Each of the long and short coil sets may be considered to generate a family of constant signal or constant voltage surfaces from the sensing coil within the navigational domain. Although the delta coils are preferably implemented with an array of appropriately designed magnetic coils (discussed below), this preferred implementation should not serve as a limitation of the present invention as it should be apparent to those skilled in the art that other magnetic configurations may be used to adequately generate the gradient fields.

As used herein, “magnetic look-up-table” (alternatively referenced as “the LUT”) refers to a database including the magnetic field values at every x-y-z coordinate position within the navigational domain for the unidirectional coils and delta coils used by the present invention. Accordingly, input data consisting of an x-y-z coordinate and a magnetic field identifier, which designates a selected magnetic coil assembly, is indexed within the database to a corresponding set of magnetic field values constituting the output data. For the x-y-z coordinate system, the output data is represented by the magnetic field variables Hx Hy Hz where the subscript indicates the axial dimension along which the magnetic field value is being reported. The database is created through a computational analysis of the magnetic field patterns generated by the magnetic coil configurations used herein. The mathematical model to develop the necessary formulae defining the field patterns may be developed, for example, from near field electromagnetic theory. An instructive text for facilitating such an analysis is “Field and Wave Electromagnetics” 2nd edition Addison Wesley (1989) by D. K. Cheng, herein incorporated by reference. The database may be stored in any type of facility including, inter alia, read-only memory, firmware, optical storage, or other types of computer storage. Additionally, the database information may be organized into any type of format such as a spreadsheet. It should be apparent to those skilled in the art that any suitable technique may be used to ascertain or record the magnetic field values for the magnetic coil assemblies used herein.

Separation of Variables Methodology

The mathematical construct underlying the present invention is a methodology termed separation of variables. In accordance with this methodology, appropriate equations are developed to isolate unknown variables in such a manner that renders the equations uniquely soluble. There are five unknown variables (φ-θ-x-y-z) that define the location and orientation of the sensing coil. A typical approach to solving for these variables would be to develop a series of coupled non-linear equations expressing the relationship among the variables. However, these equations are generally not uniquely soluble, i.e., multiple solutions are possible.

The mathematical approach used herein and predicated on the separation of variables concept is directed to the development of a series of signal expression statements functionally relating induced voltage values to certain ones of the unknown location parameters. The relationships defined by these expression statements (i.e., induced voltage equations) are such that the unknown variables are determinable by sequentially solving the expression statements. In accordance with one aspect of the present invention, a class of special magnetic fields is generated with characteristic spatial structuring and shaping that is sufficient to cause the variables (i.e., the location parameters) to separate within these induced voltage equations so as to permit resolution of the parameters x-y-z-φ-θ when the equations are sequentially solved.

In particular, a series of substantially uniform fields is successively generated in the x-, y-, and z-directions with the group of unidirectional coils, thereby developing induced voltage expressions (discussed below) involving only the variables φ and θ independent of the unknown positional variables x-y-z. This uncoupling of the variables φ-θ from the variables x-y-z as accomplished by the substantially uniform fields is specifically evident in the three induced voltage equations expressed in two unknowns (i.e., φ and θ), which are easily soluble.

The separation of variables methodology as applied to the determination of the unknown positional coordinates x-y-z is implemented with a sequence of gradient fields as generated by an appropriate group of delta coils. Since the gradient fields have components in only two of the axial dimensions, these fields induce voltages in the sensing coil that are dependent upon the magnetic field values in only these two dimensions at the sensing coil position. As a result, each gradient field generates a family of constant signal surfaces from the sensing coil, from which a constant signal surface is identified for each gradient field that produces the measured induced voltage in the sensing coil. The intersection of two such constant signal surfaces is a line along which the catheter is located. This intersection line is defined by an expression in two of the unknown positional coordinates, wherein the other equation parameters are known, i.e., the measured induced voltage values, the magnetic field values at every x-y-z coordinate for the coil groups (as supplied by the LUT), and the as-computed φ-θ orientation. If an appropriate delta coil configuration is used (e.g., three delta coil pairs), an appropriate number of such intersection lines (e.g., three) may be produced to sufficiently and uniquely resolve the x-y-z coordinates (e.g., by calculating the intersection of such three intersection lines).

Since there are five unknown variables to completely define the catheter probe location, an equal number of independent equations are needed to sufficiently describe its location. These unknowns may be determined using one coil and five magnetic fields (as described herein), two coils and three magnetic fields, or three coils and two magnetic fields.

Overview of Location Algorithm

In accordance with the present invention, a location algorithm was developed for determining the location of a sensing coil affixed to the distal end of a catheter probe that is navigated through an anatomical region of interest within a recumbent patient. The location of the sensing coil is defined by an angular orientation and positional coordinates. The angular orientation is represented by an angle φ corresponding to the angle of departure from the z-axis and an angle θ corresponding to the angle between the x-axis and the projection onto the x-y plane of the vector coincident with the longitudinal axis of the coil. In the coordinate system for describing the present invention, the z-axis coincides with the longitudinal dimension extending from the patient's head to foot. The x-axis coincides with a lateral dimension across the patient's body, and the y-axis is perpendicular to the planar top of the pallet or examination deck. These dimensions are identified as the patient is disposed in the recumbent position on the pallet. As discussed below, the angular orientation is determined from signals induced in the sensing coil in response to a sequence of substantially uniform, unidirectional fields generated successively within the navigational domain. The positional coordinates are determined from signals induced in the sensing coil in response to the gradient magnetic fields.

By way of background, the time-dependent magnetic fields projected into the navigational domain induce voltages in the coil that are representative of the orientation of the coil axis relative to the lines of magnetic flux. The development of an induced voltage in an electrical conductor in response to a changing magnetic field is defined by Faraday's law. If one considers any closed stationary path in space which is linked by a changing magnetic field, it is found that the induced voltage Vind around this path is equal to the negative time rate of change of the total magnetic flux through the closed path. Denoting a closed path with the variable C, the magnetic flus through C is given by, Ψ = ∫ S ⁢ B ·   ⁢ ⅆ S
where S is any surface bounded by the closed path C. Thus, the mathematical statement of Faraday's law is V ind = - δ δ ⁢   ⁢ t ⁢ ∫ S ⁢ B ·   ⁢ ⅆ S
Basically, the law states that a changing magnetic field will induce an electric field which exists in space regardless of whether a conducting wire is present. If a conducting wire is present in the electric field, an induced voltage will develop in the conductor. For a single-turn coil of wire of radius d located in a uniform magnetic field B=Bo sinωt, where the axis of the sensing coil is displaced at an angle θ with respect to the lines of magnetic flux, the induced voltage measured between the two open ends of the coil is expressed as: V ind = - ∂ ∂ t ⁢ ∫ S ⁢ B ·   ⁢ ⅆ S = - ∂ ∂ t ⁢ ( π ⁢   ⁢ ⅆ 2 ⁢ B o ⁢ cos ⁢   ⁢ θ ⁢   ⁢ sin ⁢   ⁢ ω ⁢   ⁢ t ) V ind = - ω ⁢   ⁢ π ⁢ ⅆ 2 ⁢ B o ⁢ cos ⁢   ⁢ θ ⁢   ⁢ cos ⁢   ⁢ ω ⁢   ⁢ t
This relationship for a single coil may be used to determine the induced voltage within a coil of N turns. Assuming that each turn of the coil is separately and equally linked by the magnetic flux (e.g., in tightly wound coils), the induced voltage within the entire coil assembly may be approximated as the summation of the induced voltages developed in each turn. Accordingly, the total voltage across the entire coil assembly is N times the induced voltage for a single turn; hence, the induced voltage Vind is equivalent to
Vind=−Nωπd2Bo cosθωt
Clearly, the induced voltage in the sensing coil will vary with changes in the angular orientation between the coil axis and the direction of the magnetic field lines.

A useful reference frame for spatially conceptualizing the interaction between the sensing coil and the magnetic fields is the Cartesian coordinate system defined by mutually perpendicular axes x-y-z. For purposes of illustration, a nonzero vector â is selected to coincide with the axis through the sensing coil of the present invention (hereinafter “coil axis”).

The angles α, β, and γ that the vector â makes with the unit coordinate vectors, î, ĵ, and {circumflex over (k)}, respectively, are called the direction angles of â; the trigonometric terms cosα, cosβ, and cosγ represent direction cosine values. Employing vector product notation, the following expressions are developed: â·î=∥â∥cosα; â·ĵ=∥â∥cosβ; and â·k=∥â∥cosγ. Referencing the induced voltage equations set forth above, these angles α, β and γ correspond to the angular displacement of the coil axis with respect to uniform fields generated along the x-axis, y-axis, and z-axis directions, respectively. Thus, the correspondence between direction cosine expressions is as follows:

    • cosα corresponds to sinφcosθ;
    • cosβ corresponds to sinθsinφ; and
    • cosγ correspond to cosφ.
      Accordingly, the following relationships illustrate the dependence of induced voltage on the orientation parameters φ and θ:
      Vx≈sinφcosθ;
    • Vy≈sinθsinφ; and
    • Vz≈cosφ,
      where the subscript indicates the direction of the magnetic field that produced the measured induced voltage.

FIG. 4 is a flowchart detailing the location algorithm according to the present invention and should be referenced in connection with the discussion below.

As noted above, the last navigation point (LNP) refers to the x-y-z positional coordinates of the sensing coil as determined by the immediately previous computation cycle of the algorithm. For the first cycle, the LNP is the center of the viewing field.

In accordance with a preferred embodiment of the present invention for implementing the location algorithm, a magnetic assembly of nine individual coil sets are used to generate the magnetic fields sufficient to develop a corresponding set of nine induced voltage signals that are fully representative of the location of the sensing coil. The nine coil sets correspond to a group of three unidirectional coil sets for generating uniform fields in the x, y, and z-directions; a first delta coil group including a short coil set at 0° and a long coil set at 0°; a second delta coil group including a short coil set at 120° and a long coil set at 120°; and a third delta coil group including a short coil set at 240° and a long coil set at 240°. The angular designations associated with the delta coil groups indicate the angle with respect to the z-axis of the coil dimension that is perpendicular to the direction of elongation of the delta coils. Accordingly, the three delta coil groups are arranged pair-wise in a circular orientation about the y-axis at angles of 0°, 120°, and 240°.

The look-up-table (LUT) consists of a database containing the magnetic field values (Hx Hy Hz) at every x-y-z coordinate location within the navigational domain for five coil sets: the unidirectional coil sets for generating the uniform fields in the x, y, and z-directions; the short coil (SC) set at 0°; and the long coil (LC) set at 0°. The magnetic field value data for the short and long coil sets at 120° and 240° may be obtained from the LUT by rotating the field vectors for the long and short coil sets at 0° by the angle (i.e., ±120°) appropriate for the given coil set. The input data for the LUT consists of the x-y-z coordinates and a designation of which coil set is being used to generate the magnetic fields. In response to this input data, the LUT supplies the magnetic field values Hx Hy Hz at the selected x-y-z coordinates for the designated coil set.

The LUT is present to speed up the operational sequence of the location algorithm. Otherwise, an undesirable computational delay exists if the required magnetic fields from the nine coil sets must be individually calculated during each iteration of the algorithm. By predetermining the magnetic field values and storing them in LUT, the location algorithm need only access the LUT to retrieve the appropriate field value without endeavoring into any complex field analysis. At x-y-z coordinates other than those for which magnetic field values are determined in the LUT, an interpolation procedure is employed to calculate the field value.

Determining Angular Orientation of the Sensing Coil

The location algorithm of the present invention initially undertakes a procedure to determine the angular orientation of the sensing coil. An assumption is first made that the coil orientation does not appreciably change during the period between cycle computations. Accordingly, the magnetic field values corresponding to the uniform field pattern at the LNP are used as an approximation for the magnetic field values at the current but as yet undetermined location.

The unidirectional coils are activated in succession, each generating a substantially uniform field that projects into the navigational domain and induces a corresponding voltage signal in the sensing coil. The induced voltage signals are measured by an appropriate detection unit coupled to a proximal end of the catheter device where an electrical connection to the sensing coil is established via suitable connection means extending along the body of the catheter device.

The LUT is then accessed three times to acquire the magnetic field values at the LNP for each of the three unidirectional coils. These values and the measured voltage signals are then substituted into the appropriate equations set forth below to solve for the unknown variables φ and θ that define the coil orientation.

As a general principle, the voltage induced within the sensing coil may be resolved into components along each of the axial dimensions as determined by the extent to which the magnetic flux density is developed along these axial dimensions. For example, a general formula for the induced voltage produced by the unidirectional coil which generates a substantially uniform field in the x-direction is as follows:
Vx=HxxK sinφcosθ+HyxK sinφsinθ+HzxK cosφ
where magnetic field intensity H is related to magnetic flux density by B=μH and K=μo ωNπd2. The first subscript in the field intensity term indicates the axial dimension along which the magnetic field value was determined by accessing the LUT for the given coil set at the LNP, while the second subscript indicates the field-generating coil set. For an x-directed substantially uniform field, the terms Hyz and Hzx are small compared to Hxx. Similar equations are developed below for the induced voltages produced by the unidirectional coils successively generating a y-directed and z-directed substantially uniform field:
Vy=HxyK sinφcosθ+HyyK sinφsinθ+HzyK cosφ,
and
Vz=HxxK sinφcosθ+HyxK sinφsinθ+HzxK cosφ.
The terms Hxy and Hzy in the equation for Vy and the terms Hxz and Hyz in the equation for Vz are small compared to Hyy and Hzz, respectively. After substituting the measured values for the induced voltage signals, the equations are simultaneously solved to determine the unknown variables φ and θ defining the orientation of the sensing coil.

Determining Positional Coordinates of the Sensing Coil

By way of summary, the procedure for determining the positional coordinates of the sensing coil in accordance with the present invention first involves activating each delta coil in succession and measuring the induced voltage thereby developed in the sensing coil. Next, the LUT is accessed to obtain the magnetic field values at the LNP for each specified delta coil. These magnetic field values and the as-computed values for the orientation angles φ and θ are then substituted into the appropriate induced voltage equations to calculate for each delta coil the expected value of the voltage signal induced in the sensing coil. This expected value of the induced signal corresponds to a specific and unique member of the family of constant signal surfaces of the delta coils.

Based on the difference between the measured and expected values for the induced voltage signals, a gradient is calculated (representative of the rate of change of the induced signal) that permits identification of the specific constant signal surface that is responsible for generating the measured value of the induced signal. This procedure is repeated for each delta coil.

For the activation of each delta coil group (comprised of one long coil set and one short coil set), there is an intersection line defined by the intersection of the two constant signal surfaces (which were identified as developing the measured induced signal) on which the sensing coil is located. The intersection of the three such lines from the three delta coil groups uniquely provides the x-y-z coordinates of the sensing coil. Although two such lines are sufficient to describe the position of the sensing coil, greater accuracy and more reliable performance in determining the catheter position is achieved with three lines.

The following is a more detailed description of the procedure summarized above for determining the positional coordinates.

The magnetic field pattern generated by the entire assembly of short coil and long coil sets is characterized by a family of surfaces of constant signal or constant voltage developed by the sensing coil, each having non-zero components in two of the axis directions and a small component in the remaining axis direction. For example, the magnetic field surfaces generated by the short and long coil sets oriented at 0° relative to the x-axis have a small value in the x-direction. The short coil positioned at 0° (i.e., SC(0°)) and long coil positioned at 0° (i.e., LC (0°)) are each independently activated. The induced voltage in the sensing coil is measured for each coil set. The LUT is then accessed to determine the magnetic field values for the SC(0°) and LC(0°) coil set at the LNP.

These magnetic field values (i.e., Hx=small and non-zero Hy Hz components) are used in conjunction with the as-computed orientation angles φ and θ to calculate the values of the induced catheter signals that would be expected from such magnetic field values. The expected and measured induced voltage values are compared, and the difference is used to identify the constant signal surface from each of the SC(0°) and LC(0°) coil sets that would have produced the measured induced signals. The intersection of these identified magnetic constant signal surfaces is a line parallel to the x-axis (thereby resolving the y-z coordinates).

The aforementioned procedure involving the long and short coils oriented at 0° is iteratively repeated for a long and short coil set oriented at 120° (i.e., SC(120°) and LC(120°)) and 240° (i.e., SC(240°) and LC(240°)).

More specifically, the coil sets SC(120°) and LC(120°) are sequentially activated to induce corresponding catheter signals in the sensing coil. In order to utilize the LUT data on the coil sets oriented at 0° for determining the magnetic field components at the LNP generated by the coil sets SC(120°) and LC(120°), a modified LNP is calculated that is equivalent to the original LNP rotationally displaced by 120°. The LUT is then accessed with the modified LNP to determine the magnetic field values generated by the SC(120°) and LC(120°) coil sets at the modified LNP. The field vectors produced by the LUT for both the long coil and short coil are then rotated (−120°) to go from the modified LNP to the actual LNP. Based upon these field values, a pair of induced catheter signals are calculated that correspond to the expected signal values arising from the magnetic field values for the SC(120°) and LC(120°) coil sets. The difference between the measured and expected induced catheter signals is used to identify the magnetic constant signal surface for each of the SC(120°) and LC(120°) coil sets that could produce the measured catheter signal. The intersection of these identified magnetic constant signal surfaces is a line oriented at 120° to the x-axis.

A similar procedure is used involving a modified LNP that is rotationally displaced 240° to simulate the magnetic field patterns for the SC(240°) and LC(240°) coil sets using the SC(0°) and LC(0°) field data. A line oriented at 240° to the x-axis is then identified along which the catheter is located.

Each of the field lines oriented at 0°, 120° and 240° to the x-axis is weighted according to the strength of the measured catheter signals. For example, a weak measurement indicates a relatively imprecise identification of the intersection line, resulting in a weaker weighting. This weighting reflects the accuracy of the estimation used to determine the location of the catheter with the specified coil set. An averaging technique is used to compute a weighted estimate of the intersection of the lines L(0°), L(120°) and L(240°). The intersection is the new value for x-y-z and will replace the x-y-z of the old LNP to become the next LNP. The algorithm iteratively repeats the aforementioned operations using the updated LNP to arrive at the location of the sensing coil after each computation cycle (e.g., every 0.1 s).

Magnetic Assembly for Determining Angular Orientation of Sensing Coil

FIG. 1 schematically illustrates a perspective view of an examination deck that facilitates implementation of the location algorithm in accordance with a preferred embodiment of the present invention, and which employs a magnetic coil assembly arranged in a flat configuration. The examination deck includes a planar top platform 10 suitable for accommodating a recumbent patient disposed lengthwise on the planar top. The navigational domain is illustratively depicted as the spherical volume 12 enclosing a sensing coil 14 attached via suitable connection means 16 to an external signal detection apparatus (not shown). The coil sets embedded in platform 10 (and described in connection with FIGS. 2A-C and 3) are activated by a signal drive unit (not shown) connected via line 18. The examination deck is preferably constructed from a suitable magnetically-permeable material to facilitate magnetic coupling between the embedded coil sets and the overlying sensing coil.

Coil Set for Generating X-Directed Field

FIG. 2A schematically illustrates the unidirectional coil set for generating a substantially uniform x-directed field throughout the navigational domain 12. The coil set includes a first coil pair with elements 20 and 24 and a second coil pair with elements 22 and 26, where the current flow as supplied by drive unit 28 is indicated by the arrow symbol. Coil elements 20 and 22 are disposed in the major surface of platform 10, while elements 24 and 26 are disposed in the lateral walls of platform 10. Elements 24 and 26 are preferably used as compensation coils to substantially cancel undesirable field components generated by elements 20 and 22 in the y- and z-directions. The coils cumulatively generate a substantially uniform x-directed field as indicated by representative field line 27.

Coil Set for Generating Y-Directed Field

FIG. 2B schematically illustrates the unidirectional coil set for generating a substantially uniform y-directed field throughout the navigational domain 12. The coil set includes a coil pair with elements 30 and 32 disposed in spaced-apart and parallel relationship within platform 10, with the indicated current flow as supplied by drive unit 34. The coils generate a substantially uniform y-directed field as indicated by representative field line 33.

Coil Set for Generating Z-Directed Field

FIG. 2C schematically illustrates the unidirectional coil set for generating a substantially uniform z-directed field throughout the navigational domain 12. The coil set includes a first coil pair with elements 36 and 40 and a second coil pair with elements 38 and 42, with the indicated current flow as supplied by drive unit 44. Coil elements 36 and 38 are disposed in the major surface of platform 10, while elements 40 and 42 are disposed in the lateral walls of platform 10. Elements 40 and 42 are preferably used as compensation coils (e.g., Cunard coils) to substantially cancel undesirable field components generated by elements 36 and 38 in the x- and y-directions. The coils cumulatively generate a substantially uniform z-directed field as indicated by representative field line 43.

The coil configurations shown in the Figures are only illustrative and should not be construed as a limitation of the present invention, as it should be apparent to those skilled in the art that other coil configurations are possible within the scope of the present invention provided such other configurations produce the desired magnetic field patterns. A suitable connection means (not shown) couples the sensing coil 14 to a signal measuring device.

Magnetic Assembly for Determining Positional Coordinates of Sensing Coil

FIGS. 3 and 5 show the coil configuration used to determine the positional coordinates of the sensing coil in accordance with a preferred embodiment of the present invention. The configuration includes six coils grouped into three pairs of long and short delta coils (50-52, 54-56, 58-60). The delta coils are mutually coplanar and are disposed in the planar top of the examination deck immediately beneath the recumbent patient. Interconnection means between a signal drive unit (not shown) and the delta coil groups is shown representatively for only coils 50-52.

The coils are preferably arranged in a circular orientation about the y-axis such that there is an axis perpendicular to the direction of elongation of the coils at 0°, 120° and 240° relative to the z-axis. The magnetic field generated by the first group of long (50) and short delta coils (52) is shown representatively by the field lines extending from the upper region of the coils. The field lines from this delta coil group form the family of constant signal surfaces shown within the navigational domain 12. Superposition of the constant signal surfaces generated by the long and short coils of a delta coil group produces a fishnet pattern as shown in FIG. 6. The intersection of two such constant signal surfaces generated by a short and long coil pair is a single line represented by the dotted line 70.

A constant signal surface (72 and 74) is identified for each short coil and long coil activation of a delta coil pair by determining the surface that matches the induced signals developed in the sensing coil. This procedure is repeated for the other two delta coil pairs to produce two other lines comparable to line 70. The intersection of these three lines determines the position of the catheter.

FIG. 7 shows an upper plan schematic view of the entire delta coil arrangement relative to an inner circular space representing the projection of the navigational domain into the plane of the delta coils. It is an object of the present invention to design coils having high spatial gradience in two of the axis dimensions and a substantially zero field value in the remaining axial dimension. This particular design is accomplished by modifying the termination points of the coils with compensation coils such that the modified coil is effectively operative as an infinitely long coil. The long coil sets are further compensated by a central “sucker” coil 88. Accordingly, each of the long coils and short coils is modified by representative compensation coils 80-82, 84-86, 88 and 90-94, 92-96 respectively, disposed at the indicated endpoints and center of the corresponding delta coil. The long coil and short coil configurations are shown schematically for only sets 50-52, but similar configurations likewise exist for the coil sets 54-56 and 58-60 shown representatively as the indicated lines.

The quality of the coils, as measured by the degree of uniformity of the uniform field coils or how close to zero is the field in the non-gradient direction for the delta coils, determines the size of the navigational domain over which the variable separation technique for navigating the catheter will converge and therefore be capable of initially finding the catheter, and hence be of functional utility.

FIG. 8 schematically depicts an examination deck in accordance with another embodiment of the present invention. The deck includes a first rail member 100 and a second rail member 102 in opposed spaced-apart relationship and attached to the platform along respective supporting edges. The navigational domain is illustratively depicted as the spherical volume 12. The deck includes an apertured opening 104. Each rail member has an inner wall and an outer wall. The railed configuration is characterized by the embedding of coil sets in both the planar top and in the rail members. The examination deck is preferably constructed from a magnetically permeable material.

FIGS. 9A-C schematically illustrate the unidirectional coils for implementing the railed configuration used in conjunction with the examination deck of FIG. 8. The magnet assembly for the x-directed unidirectional coil set is shown in FIG. 9A and includes two coil elements 110 and 112 each embedded in a respective rail member. Each coil pair is designed to project a substantially uniform field in the x-direction throughout the navigational domain. FIG. 9B schematically depicts the y-directed unidirected coils including coil elements 114 and 116 each embedded in respective rail members, and further including coil elements 118 and 120 embedded in the planar top of the examination deck. FIG. 9C schematically depicts the z-directed unidirected coils including coil elements 122-124 in one rail member and elements 126-128 in the other rail member. The current flow through each coil configuration is indicated by the arrows. FIG. 9D shows the delta coil arrangement used in the railed configuration. This arrangement is the same as used in the flat configuration described above.

In accordance with another embodiment of the present invention, a second sensing coil is used for stabilization purposes. Inaccurate readings of the catheter probe location may occur from motion artifacts due to breathing action, heart motion, or patient movement. The stabilized location coordinates may be determined by placing a second sensing coil on the sternum of the patient at a known location within the navigational domain. The incremental movement experienced by the second sensing coil due to motion artifacts is detected and subtracted from the measured location value of the probe to arrive at the actual location coordinates of the probe. Further extensions of the present invention are possible to facilitate multi-catheter applications by attaching an additional sensing coil to the distal end of each additional catheter.

Since certain changes may be made in the above apparatus and method without departing from the scope of the invention herein described, it is intended that all matter contained in the above description or shown in the accompanying drawings shall be interpreted in an illustrative and not in a limiting sense.

Claims

1. A method of determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain, comprising the steps of:

inducing within said sensing coil a set of orientation signal values each representative of an orientation of said sensing coil and independent of a position of said sensing coil;
determining the orientation of said sensing coil using said induced orientation signal values;
inducing within said sensing coil a set of positional signal values each representative of the position of said sensing coil; and
determining the position of said sensing coil using said positional signal values and said determined orientation.

2. The method as recited in claim 1, wherein the step of inducing said set of orientation signal values comprises the steps of:

generating from outside said body a series of magnetic fields each penetrating at least said navigational domain and characterized substantially by a principal magnetic component in one axial dimension and relatively smaller magnetic components in two other axial dimensions.

3. The method as recited in claim 1, wherein the step of inducing said set of positional signal values comprises the steps of:

generating from outside said body a series of magnetic fields each penetrating at least said navigational domain and characterized substantially by two principal gradient magnetic components in respective axial dimensions and a relatively smaller magnetic components in a third axial dimension.

4. The method as recited in claim 3, wherein said generating step further includes the steps of:

generating said fields to provide a plurality of constant signal surfaces for the sensing coil such that an intersection between two such surfaces with components in the same axial dimensions produces a line along which said sensing coil is located;
wherein said two such surfaces are identified from among said plurality of constant signal surfaces by their ability to induce one of said positional signal values.

5. The method as recited in claim 4, further comprises the steps of:

weighting each line in accordance with a signal strength of said corresponding constant signal surface; and
determining an intersection of said weighted lines.

6. The method as recited in claim 5, wherein six constant signal surfaces are generated to produce three intersection lines.

7. A system for determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain, comprising:

first transmit means for projecting into said navigational domain magnetic energy that is sufficient to induce signal values within said sensing coil representative of an orientation of said sensing coil and independent of the position of said sensing coil;
second transmit means for projecting into said navigational domain magnetic energy that is sufficient to induce signal values within said sensing coil representative of the position of said sensing coil; and
analysis means, coupled to said first transmit means and said second transmit means, for determining the position and orientation of said sensing coil from said induced signal values.

8. A system for determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain, comprising:

first signal-inducing means for inducing within said sensing coil orientation signals that are representative of the orientation of said sensing coil;
analysis means, coupled to said first signal-inducing means, for determining the orientation of said sensing coil using said induced orientation signals and independent from a position of said sensing coil;
second signal-inducing means for inducing within said sensing coil position signals that are representative of the position of said sensing coil; and
analysis means, coupled to said second signal-inducing means, for determining the position of said sensing coil using said determined orientation and said induced position signals.

9. The system as recited in claim 8, wherein the first signal-inducing means comprises:

field generation means for successively generating magnetic field patterns projected into said navigational domain, each characterized substantially by a principal magnetic field component in one direction and relatively smaller magnetic components in two other directions.

10. The system as recited in claim 9, wherein said field generation means comprises a set of magnetic coils.

11. The system as recited in claim 10, wherein said magnetic coils are disposed in a planar top of an examination deck upon which a patient is disposed during a surgical procedure.

12. The system as recited in claim 10, wherein said magnetic coils are disposed in a planar top and in rail members edge supported by said planar top for an examination deck upon which a patient is disposed during a surgical procedure.

13. The system as recited in claim 8, wherein the second signal-inducing means comprises:

field generation means for successively generating magnetic field patterns each characterized by a first and second gradient field component in respective directions and a relatively smaller third component in another direction.

14. The system as recited in claim 13, wherein the field generation means comprises a magnetic coil assembly.

15. A method of determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain, comprising the steps of:

defining the location of said sensing coil with a set of independent location parameters; and
sequentially generating within said navigational domain a sequence of magnetic fields for inducing within said sensing coil a corresponding sequence of induced signals each defined by an induced signal expression that functionally relates said induced signal to certain ones of said location parameters, such that said set of location parameters is determinable by sequentially solving individual signal expression groups each including certain ones of said induced signal expressions and sufficient to represent a subset of said location parameters.

16. The method as recited in claim 15, wherein said sequence of magnetic fields comprises:

a series of unidirectional magnetic fields each characterized substantially by a principal magnetic field component in one direction and relatively smaller magnetic components in two other directions; and
a series of gradient magnetic fields each characterized by a first and second gradient field component in respective directions and a relatively smaller third component in another direction.

17. The method as recited in claim 16, wherein said signal expression groups include:

an orientation group including induced signal expressions each functionally related to a respective one of said unidirectional magnetic fields and an orientation of said sensing coil, and independent of a position of said sensing coil; and
a position group including induced signal expressions each functionally related to a respective one of said gradient magnetic fields, the orientation of said sensing coil, and the position of said sensing coil.

18. The method as recited in claim 17, wherein the step of sequentially solving said individual signal expression groups includes the steps of:

initially solving the induced signal expressions of said orientation group; and
next solving the induced signal expressions of said position group.

19. A system for determining the location of a magnetically-sensitive, electrically conductive sensing coil affixed to a distal end of a catheter probe partially inserted into a body cavity within a navigational domain, comprising:

means for defining the location of said sensing coil with a set of independent location parameters; and
field generation means for sequentially generating within said navigational domain a sequence of magnetic fields for inducing within said sensing coil a corresponding sequence of induced signals each defined by an induced signal expression that functionally relates said induced signal to certain ones of said location parameters, such that said set of location parameters is determinable by sequentially solving individual signal expression groups each including certain ones of said induced signal expressions and sufficient to represent a subset of said location parameters.

20. The system as recited in claim 19, wherein said sequence of magnetic fields comprises:

a series of unidirectional magnetic fields each characterized substantially by a principal magnetic field component in one direction and relatively smaller magnetic components in two other directions; and
a series of gradient magnetic fields each characterized by a first and second gradient field component in respective directions and a relatively smaller third component in another direction.

21. The system as recited in claim 20, wherein said signal expression groups include:

an orientation group including induced signal expressions each functionally related to a respective one of said unidirectional magnetic fields and an orientation of said sensing coil, and independent of a position of said sensing coil; and
a position group including induced signal expressions each functionally related to a respective one of said gradient magnetic fields, the orientation of said sensing coil, and the position of said sensing coil.

22. The system as recited in claim 21, wherein said field generation means comprises:

analysis means for solving the induced signal expressions of said orientation group; and
analysis means for solving the induced signal expressions of said position group.

23. A method of determining the location of at least one magnetically-sensitive, electrically conductive sensing coil affixed to a medical device partially inserted into a body cavity within a navigational domain, comprising the steps of:

inducing within the at least one sensing coil a set of induced signal values corresponding to a set of location parameters, by sequentially generating within the navigational domain (i) a series of unidirectional magnetic fields, each characterized by a principal magnetic field component in one direction, and relatively smaller magnetic components in two other directions, and (ii) a series of gradient magnetic fields, each characterized by a first and second gradient field component in respective directions, and a relatively smaller third component in another direction;
wherein each of the induced signals is defined by an induced signal expression that functionally relates the induced signal to certain ones of the location parameters, such that the set of location parameters is determinable by sequentially solving individual signal expression groups each including certain ones of the induced signal expressions and sufficient to represent a subset of the location parameters.

24. A method according to claim 23, wherein the signal expression groups include (i) an orientation group including induced signal expressions each functionally related to a respective one of the unidirectional magnetic fields and an orientation of the sensing coil, and independent of a position of the sensing coil, and (ii) a position group including induced signal expressions each functionally related to a respective one of the gradient magnetic fields, the orientation of the sensing coil and the position of the sensing coil.

Referenced Cited
U.S. Patent Documents
1576781 March 1926 Phillips
1735726 November 1929 Bornhardt
2407845 September 1946 Nemeyer
2650588 September 1953 Drew
2697433 December 1954 Sehnder
3016899 January 1962 Stenvall
3017887 January 1962 Heyer
3061936 November 1962 Dobbeleer
3073310 January 1963 Mocarski
3294083 December 1966 Alderson
3367326 February 1968 Frazier
3439256 April 1969 Kahne
3577160 May 1971 White
3674014 July 1972 Tilander
3702935 November 1972 Carey et al.
3704707 December 1972 Halloran
3868565 February 1975 Kuipers
3941127 March 2, 1976 Froning
4037592 July 26, 1977 Kronner
4052620 October 4, 1977 Brunnett
4054881 October 18, 1977 Raab
4117337 September 26, 1978 Staats
4173228 November 6, 1979 Van Steenwyk et al.
4202349 May 13, 1980 Jones
4262306 April 14, 1981 Renner
4287809 September 8, 1981 Egli et al.
4314251 February 2, 1982 Raab
4317078 February 23, 1982 Weed et al.
4328813 May 11, 1982 Ray
4339953 July 20, 1982 Iwasaki
4358856 November 9, 1982 Stivender et al.
4368536 January 11, 1983 Pfeiler
4396885 August 2, 1983 Constant
4418422 November 29, 1983 Richter et al.
4422041 December 20, 1983 Lienau
4431005 February 14, 1984 McCormick
4485815 December 4, 1984 Amplatz
4543959 October 1, 1985 Sepponen
4548208 October 22, 1985 Niemi
4572198 February 25, 1986 Codrington
4584577 April 22, 1986 Temple
4613866 September 23, 1986 Blood
4618978 October 21, 1986 Cosman
4621628 November 11, 1986 Buldermann
4625718 December 2, 1986 Olerud et al.
4642786 February 10, 1987 Hansen
4645343 February 24, 1987 Stockdale et al.
4649504 March 10, 1987 Krouglicof et al.
4651732 March 24, 1987 Frederick
4653509 March 31, 1987 Oloff et al.
4673352 June 16, 1987 Hansen
4706665 November 17, 1987 Gouda
4719419 January 12, 1988 Dawley
4722056 January 26, 1988 Roberts et al.
4722336 February 2, 1988 Kim et al.
4727565 February 23, 1988 Ericson
4737794 April 12, 1988 Jones
4737921 April 12, 1988 Goldwasser et al.
4750487 June 14, 1988 Zanetti
4771787 September 20, 1988 Wurster et al.
4791934 December 20, 1988 Brunnett
4793355 December 27, 1988 Crum et al.
4797907 January 10, 1989 Anderton
4803976 February 14, 1989 Frigg et al.
4821206 April 11, 1989 Arora
4821731 April 18, 1989 Martinelli et al.
4836778 June 6, 1989 Baumrind et al.
4845771 July 4, 1989 Wislocki et al.
4849692 July 18, 1989 Blood
4862893 September 5, 1989 Martinelli
4889526 December 26, 1989 Rauscher et al.
4905698 March 6, 1990 Strohl et al.
4923459 May 8, 1990 Nambu
4931056 June 5, 1990 Ghajar et al.
4945305 July 31, 1990 Blood
4945914 August 7, 1990 Allen
4951653 August 28, 1990 Fry et al.
4977655 December 18, 1990 Martinelli
4989608 February 5, 1991 Ratner
4991579 February 12, 1991 Allen
5002058 March 26, 1991 Martinelli
5005592 April 9, 1991 Cartmell
5013317 May 7, 1991 Cole et al.
5016639 May 21, 1991 Allen
5027818 July 2, 1991 Bova et al.
5030196 July 9, 1991 Inoue
5030222 July 9, 1991 Calandruccio et al.
5031203 July 9, 1991 Trecha
5042486 August 27, 1991 Pfeiler et al.
5050608 September 24, 1991 Watanabe et al.
5054492 October 8, 1991 Scribner et al.
5057095 October 15, 1991 Fabian
5059789 October 22, 1991 Salcudean
5079699 January 7, 1992 Tuy et al.
5086401 February 4, 1992 Glassman et al.
5094241 March 10, 1992 Allen
5097839 March 24, 1992 Allen
5099845 March 31, 1992 Besz et al.
5105829 April 21, 1992 Fabian et al.
5107839 April 28, 1992 Houdek et al.
5107843 April 28, 1992 Aarnio et al.
5107862 April 28, 1992 Fabian et al.
5109194 April 28, 1992 Cantaloube
5119817 June 9, 1992 Allen
5142930 September 1, 1992 Allen et al.
5152288 October 6, 1992 Hoenig et al.
5160337 November 3, 1992 Cosman
5161536 November 10, 1992 Vilkomerson et al.
5178164 January 12, 1993 Allen
5178621 January 12, 1993 Cook et al.
5186174 February 16, 1993 Schlondorff et al.
5187475 February 16, 1993 Wagener et al.
5188126 February 23, 1993 Fabian et al.
5190059 March 2, 1993 Fabian et al.
5197476 March 30, 1993 Nowacki et al.
5197965 March 30, 1993 Cherry et al.
5198768 March 30, 1993 Keren
5198877 March 30, 1993 Schulz
5211164 May 18, 1993 Allen
5211165 May 18, 1993 Dumoulin et al.
5212720 May 18, 1993 Landi et al.
5214615 May 25, 1993 Bauer
5219351 June 15, 1993 Teubner et al.
5222499 June 29, 1993 Allen et al.
5228442 July 20, 1993 Imran
5233990 August 10, 1993 Barnea
5237996 August 24, 1993 Waldman et al.
5249581 October 5, 1993 Horbal et al.
5251127 October 5, 1993 Raab
5251635 October 12, 1993 Dumoulin et al.
5253647 October 19, 1993 Takahashi et al.
5255680 October 26, 1993 Darrow et al.
5257636 November 2, 1993 White
5265610 November 30, 1993 Darrow et al.
5265611 November 30, 1993 Hoenig et al.
5269759 December 14, 1993 Hernandez et al.
5271400 December 21, 1993 Dumoulin et al.
5273025 December 28, 1993 Sakiyama et al.
5274551 December 28, 1993 Corby, Jr.
5279309 January 18, 1994 Taylor et al.
5291199 March 1, 1994 Overman et al.
5295483 March 22, 1994 Nowacki et al.
5297549 March 29, 1994 Beatty et al.
5299254 March 29, 1994 Dancer et al.
5299288 March 29, 1994 Glassman et al.
5305091 April 19, 1994 Gelbart et al.
5305203 April 19, 1994 Raab
5309913 May 10, 1994 Kormos et al.
5315630 May 24, 1994 Sturm et al.
5316024 May 31, 1994 Hirschi et al.
5318025 June 7, 1994 Dumoulin et al.
5320111 June 14, 1994 Livingston
5325728 July 5, 1994 Zimmerman et al.
5325873 July 5, 1994 Hirschi et al.
5329944 July 19, 1994 Fabian et al.
5333168 July 26, 1994 Fernandes et al.
5353795 October 11, 1994 Souza et al.
5353800 October 11, 1994 Pohndorf et al.
5353807 October 11, 1994 DeMarco
5368030 November 29, 1994 Zinreich et al.
5375596 December 27, 1994 Twiss et al.
5377678 January 3, 1995 Dumoulin et al.
5383454 January 24, 1995 Bucholz
5385146 January 31, 1995 Goldreyer
5385148 January 31, 1995 Lesh et al.
5386828 February 7, 1995 Owens et al.
5389101 February 14, 1995 Heilbrun et al.
5391199 February 21, 1995 Ben-Haim
5394457 February 28, 1995 Leibinger et al.
5397329 March 14, 1995 Allen
5399146 March 21, 1995 Nowacki et al.
5400384 March 21, 1995 Fernandes et al.
5402801 April 4, 1995 Taylor
5403321 April 4, 1995 DiMarco
5408409 April 18, 1995 Glassman et al.
5417210 May 23, 1995 Funda et al.
5419325 May 30, 1995 Dumoulin et al.
5423334 June 13, 1995 Jordan
5425367 June 20, 1995 Shapiro et al.
5425382 June 20, 1995 Golden et al.
5426683 June 20, 1995 O'Farrell, Jr. et al.
5426687 June 20, 1995 Goodall et al.
5427097 June 27, 1995 Depp
5429132 July 4, 1995 Guy et al.
5433198 July 18, 1995 Desai
RE35025 August 22, 1995 Anderton
5437277 August 1, 1995 Dumoulin et al.
5443066 August 22, 1995 Dumoulin et al.
5443489 August 22, 1995 Ben-Haim
5444756 August 22, 1995 Pai et al.
5445144 August 29, 1995 Wodicka et al.
5445150 August 29, 1995 Dumoulin et al.
5445166 August 29, 1995 Taylor
5446548 August 29, 1995 Gerig et al.
5447154 September 5, 1995 Cinquin et al.
5448610 September 5, 1995 Yamamoto et al.
5453686 September 26, 1995 Anderson
5456718 October 10, 1995 Szymaitis
5458718 October 17, 1995 Venkitachalam
5464446 November 7, 1995 Dreessen et al.
5478341 December 26, 1995 Cook et al.
5478343 December 26, 1995 Ritter
5480422 January 2, 1996 Ben-Haim
5483961 January 16, 1996 Kelly et al.
5485849 January 23, 1996 Panescu et al.
5487391 January 30, 1996 Panescu
5487729 January 30, 1996 Avellanet et al.
5487757 January 30, 1996 Truckai et al.
5490196 February 6, 1996 Rudich et al.
5494034 February 27, 1996 Schlondorff et al.
5503416 April 2, 1996 Aoki et al.
5513637 May 7, 1996 Twiss et al.
5515160 May 7, 1996 Schulz et al.
5517990 May 21, 1996 Kalfas et al.
5531227 July 2, 1996 Schneider
5531520 July 2, 1996 Grimson et al.
5542938 August 6, 1996 Avellanet et al.
5543951 August 6, 1996 Moehrmann
5546940 August 20, 1996 Panescu et al.
5546949 August 20, 1996 Frazin et al.
5546951 August 20, 1996 Ben-Haim
5551429 September 3, 1996 Fitzpatrick et al.
5558091 September 24, 1996 Acker et al.
5568809 October 29, 1996 Ben-Haim
5572999 November 12, 1996 Funda et al.
5573533 November 12, 1996 Strul
5575794 November 19, 1996 Walus et al.
5583909 December 10, 1996 Hanover
5588430 December 31, 1996 Bova et al.
5592939 January 14, 1997 Martinelli
5595193 January 21, 1997 Walus et al.
5596228 January 21, 1997 Anderton et al.
5600330 February 4, 1997 Blood
5603318 February 18, 1997 Heilbrun et al.
5617462 April 1, 1997 Spratt
5617857 April 8, 1997 Chader et al.
5619261 April 8, 1997 Anderton
5622169 April 22, 1997 Golden et al.
5622170 April 22, 1997 Schultz
5627873 May 6, 1997 Hanover et al.
5628315 May 13, 1997 Vilsmeier et al.
5630431 May 20, 1997 Taylor
5636644 June 10, 1997 Hart et al.
5638819 June 17, 1997 Manwaring et al.
5640170 June 17, 1997 Anderson
5642395 June 24, 1997 Anderton et al.
5643268 July 1, 1997 Vilsmeier et al.
5645065 July 8, 1997 Shapiro et al.
5647361 July 15, 1997 Damadian
5662111 September 2, 1997 Cosman
5664001 September 2, 1997 Tachibana et al.
5674296 October 7, 1997 Bryan et al.
5676673 October 14, 1997 Ferre et al.
5681260 October 28, 1997 Ueda et al.
5682886 November 4, 1997 Delp et al.
5690108 November 25, 1997 Chakeres
5694945 December 9, 1997 Ben-Haim
5695500 December 9, 1997 Taylor et al.
5695501 December 9, 1997 Carol et al.
5697377 December 16, 1997 Wittkampf
5702406 December 30, 1997 Vilsmeier et al.
5711299 January 27, 1998 Manwaring et al.
5713946 February 3, 1998 Ben-Haim
5715822 February 10, 1998 Watkins
5715836 February 10, 1998 Kliegis et al.
5718241 February 17, 1998 Ben-Haim et al.
5727552 March 17, 1998 Ryan
5727553 March 17, 1998 Saad
5729129 March 17, 1998 Acker
5730129 March 24, 1998 Darrow
5730130 March 24, 1998 Fitzpatrick et al.
5732703 March 31, 1998 Kalfas et al.
5735278 April 7, 1998 Hoult et al.
5738096 April 14, 1998 Ben-Haim
5741214 April 21, 1998 Ouchi et al.
5742394 April 21, 1998 Hansen
5744953 April 28, 1998 Hansen
5748767 May 5, 1998 Raab
5749362 May 12, 1998 Funda et al.
5749835 May 12, 1998 Glantz
5752513 May 19, 1998 Acker et al.
5755725 May 26, 1998 Druais
RE35816 June 2, 1998 Schulz
5758667 June 2, 1998 Slettenmark
5762064 June 9, 1998 Polvani
5767669 June 16, 1998 Hansen et al.
5769789 June 23, 1998 Wang et al.
5769861 June 23, 1998 Vilsmeier
5772594 June 30, 1998 Barrick
5775322 July 7, 1998 Silverstein et al.
5776064 July 7, 1998 Kalfas et al.
5782765 July 21, 1998 Jonkman
5787886 August 4, 1998 Kelly et al.
5792055 August 11, 1998 McKinnon
5795294 August 18, 1998 Luber et al.
5797849 August 25, 1998 Vesely et al.
5799055 August 25, 1998 Peshkin et al.
5799099 August 25, 1998 Wang et al.
5800352 September 1, 1998 Ferre et al.
5800535 September 1, 1998 Howard, III
5802719 September 8, 1998 O'Farrell, Jr. et al.
5803089 September 8, 1998 Ferre et al.
5807252 September 15, 1998 Hassfeld et al.
5810728 September 22, 1998 Kuhn
5810735 September 22, 1998 Halperin et al.
5823192 October 20, 1998 Kalend et al.
5823958 October 20, 1998 Truppe
5828725 October 27, 1998 Levinson
5829444 November 3, 1998 Ferre et al.
5831260 November 3, 1998 Hansen
5833608 November 10, 1998 Acker
5834759 November 10, 1998 Glossop
5836954 November 17, 1998 Heilbrun et al.
5840024 November 24, 1998 Taniguchi et al.
5840025 November 24, 1998 Ben-Haim
5843076 December 1, 1998 Webster, Jr. et al.
5848967 December 15, 1998 Cosman
5851183 December 22, 1998 Bucholz
5865846 February 2, 1999 Bryan et al.
5868674 February 9, 1999 Glowinski et al.
5868675 February 9, 1999 Henrion et al.
5871445 February 16, 1999 Bucholz
5871455 February 16, 1999 Ueno
5871487 February 16, 1999 Warner et al.
5873822 February 23, 1999 Ferre et al.
5884410 March 23, 1999 Prinz
5891034 April 6, 1999 Bucholz
5891157 April 6, 1999 Day et al.
5904691 May 18, 1999 Barnett et al.
5907395 May 25, 1999 Schultz et al.
5913820 June 22, 1999 Bladen et al.
5920395 July 6, 1999 Schulz
5921992 July 13, 1999 Costales et al.
5923727 July 13, 1999 Navab
5928248 July 27, 1999 Acker
5938603 August 17, 1999 Ponzi
5938694 August 17, 1999 Jaraczewski et al.
5947981 September 7, 1999 Cosman
5950629 September 14, 1999 Taylor et al.
5951475 September 14, 1999 Gueziec et al.
5954647 September 21, 1999 Bova et al.
5964796 October 12, 1999 Imran
5967980 October 19, 1999 Ferre et al.
5968047 October 19, 1999 Reed
5971997 October 26, 1999 Guthrie et al.
5976156 November 2, 1999 Taylor et al.
5980535 November 9, 1999 Barnett et al.
5983126 November 9, 1999 Wittkampf
5987349 November 16, 1999 Schulz
5987960 November 23, 1999 Messner et al.
5999837 December 7, 1999 Messner et al.
5999840 December 7, 1999 Grimson et al.
6001130 December 14, 1999 Bryan et al.
6006126 December 21, 1999 Cosman
6016439 January 18, 2000 Acker
6019725 February 1, 2000 Vesely et al.
6024695 February 15, 2000 Taylor et al.
6050724 April 18, 2000 Schmitz et al.
6059718 May 9, 2000 Taniguchi et al.
6063022 May 16, 2000 Ben-Haim
6073043 June 6, 2000 Schneider
6104944 August 15, 2000 Martinelli
6118845 September 12, 2000 Simon et al.
6122538 September 19, 2000 Sliwa, Jr. et al.
6131396 October 17, 2000 Duerr et al.
6139183 October 31, 2000 Graumann
6149592 November 21, 2000 Yanof et al.
6156067 December 5, 2000 Bryan et al.
6175756 January 16, 2001 Ferre et al.
6223067 April 24, 2001 Vilsmeier
6273896 August 14, 2001 Franck et al.
6298262 October 2, 2001 Franck et al.
6332089 December 18, 2001 Acker et al.
6341231 January 22, 2002 Ferre et al.
6351659 February 26, 2002 Vilsmeier
6434415 August 13, 2002 Foley et al.
6445943 September 3, 2002 Ferre et al.
6498944 December 24, 2002 Ben-Haim et al.
6701179 March 2, 2004 Martinelli et al.
Foreign Patent Documents
964149 March 1975 CA
3042343 June 1982 DE
3831278 March 1989 DE
4233978 April 1994 DE
10085137 November 2002 DE
0 319 844 January 1988 EP
0350996 January 1990 EP
0 651 968 August 1990 EP
0 419 729 April 1991 EP
0 581 704 July 1993 EP
0655138 August 1993 EP
0894473 January 1995 EP
2417970 February 1979 FR
2765738 June 1998 JP
WO 88/09151 December 1988 WO
WO 89/05123 June 1989 WO
WO 91/03982 April 1991 WO
WO 91/04711 April 1991 WO
WO 91/07726 May 1991 WO
WO 92/03090 March 1992 WO
WO 92/06645 April 1992 WO
WO 94/04938 March 1994 WO
WO 94/23647 October 1994 WO
WO 94/24933 November 1994 WO
WO 96/11624 April 1996 WO
WO- 9641119 December 1996 WO
WO 98/08554 March 1998 WO
WO 98/38908 September 1998 WO
WO 99/60939 December 1999 WO
WO 01/30437 May 2001 WO
Other references
  • “Prestige Cervical Disc System Surgical Technique”, 12 pgs, Author/Date?.
  • Adams et al., “Orientation Aid for Head and Neck Surgeons,” Innov. Tech. Biol. Med., vol. 13, No. 4, 1992, pp. 409-424.
  • Barrick et al., “Prophylactic Intramedullary Fixation of the Tibia for Stress Fracture in a Professional Athlete,” Journal of Orthopaedic Trauma, vol. 6, No. 2, pp. 241-244 (1992).
  • Barrick et al., “Technical Difficulties with the Brooker-Wills Nail in Acute Fractures of the Femur,” Journal of Orthopaedic Trauma, vol. 6, No. 2, pp. 144-150 (1990).
  • Barrick, “Distal Locking Screw Insertion Using a Cannulated Drill Bit: Technical Note,” Journal of Orthopaedic Trauma, vol. 7, No. 3, 1993, pp. 248-251.
  • Batnitzky et al., “Three-Dimensinal Computer Reconstructions of Brain Lesions from Surface Contours Provided by Computed Tomography: A Prospectus,” Neurosurgery, vol. 11, No. 1, Part 1, 1982, pp. 73-84.
  • Bouazza-Marouf et al.; “Robotic-Assisted Internal Fixation of Femoral Fractures”, IMECHE., pp. 51-58 (1995).
  • Brack et al., “Accurate X-ray Based Navigation in Computer-Assisted Orthopedic Surgery,” CAR '98, pp. 716-722.
  • Bryan, “Bryan Cervical Disc System Single Level Surgical Technique”, Spinal Dynamics, 2002, pp. 1-33.
  • Bucholz et al., “Variables affecting the accuracy of stereotactic localizationusing computerized tomography,” Journal of Neurosurgery, vol. 79, Nov. 1993, pp. 667-673.
  • Champleboux et al., “Accurate Calibration of Cameras and Range Imaging Sensors: the NPBS Method,” IEEE International Conference on Robotics and Automation, Nice, France, May, 1992.
  • Champleboux, “Utilisation de Fonctions Splines pour la Mise au Point D'un Capteur Tridimensionnel sans Contact,” Quelques Applications Medicales, Jul. 1991.
  • Cinquin et al., “Computer Assisted Medical Interventions,” IEEE Engineering in Medicine and Biology, May/Jun. 1995, pp. 254-263.
  • Cinquin et al., “Computer Assisted Medical Interventions,” International Advanced Robotics Programme, Sep. 1989, pp. 63-65.
  • Clarysse et al., “A Computer-Assisted System for 3-D Frameless Localization in Sterotaxic MRI,” IEEE Transactions on Medical Imaging, vol. 10, No. 4, Dec. 1991, pp. 523-529.
  • Feldmar et al., “3D-2D Projective Registration of Free-Form Curves and Surfaces,” Rapport de recherche (Inria Sophia Antipolis), 1994, pp. 1-44.
  • Foley et al., “Fundamentals of Interactive Computer Graphics,” The Systems Programming Series, Chapter 7, Jul. 1984, pp. 245-266.
  • Foley et al., “Image-guided Intraoperative Spinal Localization,” Intraoperative Neuroprotection, Chapter 19, 1996, pp. 325-340.
  • Foley, “The StealthStation: Three-Dimensional Image-Interactive Guidance for the Spine Surgeon,” Spinal Frontiers, Apr. 1996, pp. 7-9.
  • Gildenberg et al., “Calculation of Stereotactic Coordinates from the Computed Tomographic Scan,” Neurosurgery, vol. 10, No. 5, May 1982, pp. 580-586.
  • Gonzalez, “Digital Image Fundamentals,” Digital Image Processing, Second Edition, 1987, pp. 52-54.
  • Gottesfeld Brown et al., “Registration of Planar Film Radiographs with Computer Tomography,” Proceedings of MMBIA, Jun. 1996, pp. 42-51.
  • Gueziec et al., “Registration of Computed Tomography Data to a Surgical Robot Using Fluoroscopy: A Feasibility Study,” Computer Science/Mathematics, Sep. 27, 1996, 6 pages.
  • Hamadeh et al, “Kimematic Study of Lumbar Spine Using Functional Radiographies and 3D/2D Registration,” TIMC UMR 5525-IMAG.
  • Hamadeh et al., “Automated 3-Dimensional Computed Tomographic and Fluoroscopic Image Registration,” Computer Aided Surgery (1998), 3:11-19.
  • Hamadeh et al., “Towards Automatic Registration Between CT and X-ray Images: Cooperation Between 3D/2D Registration and 2D Edge Detection,” MRCAS '95, pp. 39-46.
  • Hatch, “Reference-Display System for the Integration of CT Scanning and the Operating Microscope,” Thesis, Thayer School of Engineering, Oct. 1984, pp. 1-189.
  • Heilbrun et al., “Preliminary experience with Brown-Roberts-Wells (BRW) computerized tomography stereotaxic guidance system,” Journal of Neurosurgery, vol. 59, Aug. 1983, pp. 217-222.
  • Henderson et al., “An Accurate and Ergonomic Method of Registration for Image-guided Neurosurgery,” Computerized Medical Imaging and Graphics, vol. 18, No. 4, Jul.-Aug. 1994, pp. 273-277.
  • Hoerenz, “The Operating Microscope I. Optical Principles, Illumination Systems, and Support Systems, ” Journal of Microsurgery, vol. 1, 1980, pp. 364-369.
  • Hofstetter et al., “Fluoroscopy Based Surgical Navigation—Concept and Clinical Applications,” Computer Assisted Radiology and Surgery, 1997, pp. 956-960.
  • Horner et al., “A Comparison of CT-Stereotaxic Brain Biopsy Techniques,” Investigative Radiology, Sep.-Oct. 1984, pp. 367-373.
  • Hounsfield, “Computerized transverse axial scanning (tomography): Part 1. Description of system,” British Journal of Radiology, vol. 46, No. 552, Dec. 1973, pp. 1016-1022.
  • Jacques et al., “A Computerized Microstereotactic Method to Approach, 3-Dimensionally Reconstruct, Remove and Adjuvantly Treat Small CNS Lesions,” Applied Neurophysiology, vol. 43, 1980, pp. 176-182.
  • Jacques et al., “Computerized three-dimensional stereotaxic removal of small central nervous system lesion in patients,” J. Neurosurg., vol. 53, Dec. 1980, pp. 816-820.
  • Joskowicz et al., “Computer-Aided Image-Guided Bone Fracture Surgery: Concept and Implementation,” CAR '98, pp. 710-715.
  • Kelly et al., “Computer-assisted stereotaxic laser resection of intra-axial brain neoplasms,” Journal of Neurosurgery, vol. 64, Mar. 1986, pp. 427-439.
  • Kelly et al., “Precision Resection of Intra-Axial CNS Lesions by CT-Based Stereotactic Craniotomy and Computer Monitored CO2 Laser,” Acta Neurochirurgica, vol. 68, 1983, pp. 1-9.
  • Laitinen et al., “An Adapter for Computed Tomography-Guided, Stereotaxis,” Surg. Neurol., 1985, pp. 559-566.
  • Laitinen, “Noninvasive multipurpose stereoadapter,” Neurological Research, Jun. 1987, pp. 137-141.
  • Lavallee et al, “Matching 3-D Smooth Surfaces with their 2-D Projections using 3-D Distance Maps,” SPIE, vol. 1570, Geometric Methods in Computer Vision, 1991, pp. 322-336.
  • Lavallee et al., “Computer Assisted Driving of a Needle into the Brain,” Proceedings of the International Symposium CAR '89, Computer Assisted Radiology, 1989, pp. 416-420.
  • Lavallee et al., “Computer Assisted Interventionist Imaging: The Instance of Stereotactic Brain Surgery,” North-Holland MEDINFO 89, Part 1, 1989, pp. 613-617.
  • Lavallee et al., “Computer Assisted Spine Surgery: A Technique For Accurate Transpedicular Screw Fixation Using CT Data and a 3-D Optical Localizer,” TIMC, Faculte de Medecine de Grenoble.
  • Lavallee et al., “Image guided operating robot: a clinical application in stereotactic neurosurgery,” Proceedings of the 1992 IEEE Internation Conference on Robotics and Automation, May 1992, pp. 618-624.
  • Lavallee et al., “Matching of Medical Images for Computed and Robot Assisted Surgery,” IEEE EMBS, Orlando, 1991.
  • Lavallee, “A New System for Computer Assisted Neurosurgery,” IEEE Engineering in Medicine & Biology Society 11th Annual International Conference, 1989, pp. 0926-0927.
  • Lavallee, “VI Adaption de la Methodologie a Quelques Applications Cliniques,” Chapitre VI, pp. 133-148.
  • Leksell et al., “Stereotaxis and Tomography—A Technical Note,” ACTA Neurochirurgica, vol. 52, 1980, pp. 1-7.
  • Lemieux et al., “A Patient-to-Computed-Tomography Image Registration Method Based on Digitally Reconstructed Radiographs,” Med. Phys. 21 (11), Nov. 1994, pp. 1749-1760.
  • Levin et al., “The Brain: Integrated Three-dimensional Display of MR and PET Images,” Radiology, vol. 172, No. 3, Sep. 1989, pp. 783-789.
  • Mazier et al., “Computer-Assisted Interventionist Imaging: Application to the Vertebral Column Surgery,” Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 12, No. 1, 1990, pp. 0430-0431.
  • Mazier et al., Chirurgie de la Colonne Vertebrale Assistee par Ordinateur: Appication au Vissage Pediculaire, Innov. Tech. Biol. Med., vol. 11, No. 5, 1990, pp. 559-566.
  • Pelizzari et al., “Accurate Three-Dimensional Registration of CT, PET, and/or MR Images of the Brain,” Journal of Computer Assisted Tomography, Jan./Feb. 1989, pp. 20-26.
  • Pelizzari et al., “Interactive 3D Patient-Image Registration,” Information Processing in Medical Imaging, 12th International Conference, IPMI '91, Jul. 7-12, 136-141 (A.C.F. Colchester et al. eds. 1991).
  • Pelizzari et al., No. 528—“Three Dimensional Correlation of PET, CT and MRI Images,” The Journal of Nuclear Medicine, vol. 28, No. 4, Apr. 1987, p. 682.
  • Phillips et al., “Image Guided Orthopaedic Surgery Design and Analysis,” Trans Inst. MC, vol. 17, No. 5, 1995, pp. 251-264.
  • Potamianos et al., “Intra-Operative Imaging Guidance for Keyhole Surgery Methodology and Calibration,” First International Symposium on Medical Robotics and Computer Assisted Surgery, Sep. 22-24, 1994, pp. 98-104.
  • Reinhardt et al., “CT-Guided ‘Real Time’ Stereotaxy,” ACTA Neurochirurgica, 1989.
  • Roberts et al., “A frameless stereotaxic integration of computerized tomographic imaging and the operating microscope,” J. Neurosurg., vol. 65, Oct. 1986, pp. 545-549.
  • Rosenbaum et al., “Computerized Tomography Guided Stereotaxis: A New Approach,” Applied Neurophysiology, vol. 43, No. 3-5, 1980, pp. 172-173.
  • Sautot, “Vissage Pediculaire Assiste Par Ordinateur,” Sep. 20, 1994.
  • Schueler et al., “Correction of Image Intensifier Distortion for Three-Dimensional X-Ray Angiography,” SPIE Medical Imaging 1995, vol. 2432, pp. 272-279.
  • Selvik et al., “A Roentgen Stereophotogrammetric System,” Acta Radiologica Diagnosis, 1983, pp. 343-352.
  • Shelden et al., “Development of a computerized microsteroetaxic method for localization and removal of minute CNS lesions under direct 3-D vision,” J. Neurosurg., vol. 52, 1980, pp. 21-27.
  • Smith et al., “Computer Methods for Improved Diagnostic Image Display Applied to Stereotactic Neurosurgery,” Automedical, vol. 14, 1992, pp. 371-382.
  • Smith et al., “The Neurostation™—A Highly Accurate, Minimally Invasive Solution to Frameless Stereotactic Neurosurgery,” Computerized Medical Imaging and Graphics, vol. 18, Jul.-Aug. 1994, pp. 247-256.
  • The Laitinen Stereotactic System, E2-E6 author/date?.
  • Viant et al., “A Computer Assisted Orthopaedic System for Distal Locking of Intramedullary Nails,” Proc. of MediMEC '95, Bristol, 1995, pp. 86-91.
  • Watanabe et al., “Three-Dimensional Digitizer (Neuronavigator): New Equipment for Computed Tomography-Guided Stereotaxic Surgery,” Surgical Neurology, vol. 27, No. 6, Jun. 1987, pp. 543-547.
  • Watanabe, “Neuronavigator,” Igaku-no-Ayumi, vol. 137, No. 6, May 10, 1986, pp. 1-4.
  • Weese et al., “An Approach to 2D/3D Registration of a Vertebra in 2D X-ray Fluoroscopies with 3D CT Images,” pp. 119-128.
  • Germano, “Instrumentation, Technique and Technology”, Neurosurgery, vol. 37, No. 2, Aug. 1995, pp. 348-350.
  • Merloz, et al., “Computer Assisted Spine Surgery”, Clinical Assisted Spine Surgery, No. 337, pp. 86-96.
  • Hatch, et al., “Reference-Display System for the Integration of CT Scanning and the Operating Microscope”, Proceedings of the Eleventh Annual Northeast Bioengineering Conference, May, 1985, pp. 252-254.
  • Edward C. Benzel et al., “Magnetic Source Imaging: a Review of the Magnes System of Biomagnetic Technologies Incorporated,” Nurosurgery, vol. 33, No. 2 (Aug. 1993), p. 252-259.
  • European Search Report for EP0805069 mailed Oct. 7, 2008, which is a divisional of EPSN 96919360.6, filed Jun. 11, 1996; which is a national phase of PCT/US1996/010050, filed Jun. 11, 1996; which claims priority to U.S. Appl. No. 08/490342, filed Jun. 14, 1995.
Patent History
Patent number: RE41066
Type: Grant
Filed: Jan 14, 1999
Date of Patent: Dec 29, 2009
Assignee: Metronic Navigation, Inc. (Louisville, CO)
Inventors: Michael A. Martinelli (Winchester, MA), Wayne C. Haase (Sterling, MA)
Primary Examiner: Ruth S Smith
Attorney: Harness, Dickey & Pierce, P.L.C.
Application Number: 09/231,854