The instant invention relates generally to medical imaging and navigation and more particularly to a system and method for compensation of motion in a moving organ using an internal reference sensor.
Systems and methods for obtaining and displaying two-dimensional and three-dimensional images are known in the art, for example, as seen by reference to U.S. Pat. No. 7,386,339 entitled “MEDICAL IMAGING AND NAVIGATION SYSTEM” to Strommer et al., hereby incorporated by reference in its entirety. Strommer et al. disclose a medical imaging and navigation system that has a capability for constructing and displaying three-dimensional images of moving organs, synchronously with the actual movement of these organs and synchronously with an invasive surgical tool, such as a catheter. The system includes a medical positioning system (MPS) for ascertaining the location and orientation of multiple MPS sensors, a two-dimensional imaging system having an image detector for obtaining two-dimensional images of the moving organ and a superimposing processor. The MPS system includes a sensor mounted on the surgical tool and a sensor attached to the body of the patient for a positional reference (“Patient Reference Sensor”, or PRS). The system acquires a plurality of two-dimensional images (and respective location/orientation data and organ timing data, e.g., ECG signal) and records the sets of positions and orientation of all sensors. The system reconstructs a three-dimensional image from the combination of 2-D images and sensor data. When a physician inserts the surgical tool into the body of the patient, the system also detects the location and orientation of the MPS sensor that is mounted on the tool. The superimposing processor super-imposes a representation of the surgical tool on the currently displayed two-dimensional and three-dimensional images, which may be played back in accordance with real-time ECG data.
The PRS is provided so that the sensors associated with the surgical tools remain in a co-registered coordinate system to the X-ray imager at all times. The system detects movements of the patient using the PRS (e.g., patient body movements and respiration induced movements). The movements (as sensed by the PRS) are used to shift the coordinate system relative to the coordinate system in which the two-dimensional images were acquired. Therefore, in Strommer et al., the projection of real-time location information on previously recorded 2-D or 3-D images is both ECG synchronized and respiration compensated. However, in some situations, there is little or no correlation between the external motion compensation signals being used (i.e., the ECG signal and the PRS readings) and the internal motion of a region of interest. For example only, in the case of atrial fibrillation, the ECG signal may not effectively serve as a predictor or correlation input for the motion of the atria.
There is therefore a need for a system and method for compensation for the motion of a moving organ that minimizes or eliminates one or more of the problems set forth above.
One advantage of the methods and apparatus described, depicted and claimed herein relates to the ability to accurately compensate for the motion of a moving region of interest in a patient's body (e.g., a moving organ such as the heart), such as may be needed when superimposing a representation of a catheter tip on an image acquired at a time different than a time when the position of the catheter tip was acquired.
This disclosure is directed to a method and apparatus for displaying a moving region of interest (ROI) located within a body. One embodiment of the method involves tracking the motion of the ROI over time and generating a motion compensation function. Next, determining a position and orientation (P&O) of an invasive medical device, such as, for example, a catheter. The next step involves correcting the determined P&O using the motion compensation function to thereby compensate for the motion of the ROI between a first time at which an image of the ROI was acquired and a second time (different than the first time) at which the P&O was determined. The next step involves superimposing a representation of the medical device onto the image in accordance with the corrected P&O.
In a preferred embodiment, tracking the motion of the region of interest involves associating a first localization sensor with the moving region of interest such that the sensor moves with the region of interest. Through this step, the localization system (e.g., a medical position system (MPS) in one embodiment) can acquire a first series of P&O readings, which readings define not only the motion of the sensor but also the motion of the region of interest. The method further involves acquiring a second series of P&O readings from a second localization sensor associated with the medical device. When the correlation between the first and second series of P&O readings exceeds a threshold, the system is enabled to perform motion compensation since the same motion of the region of interest can be assumed to influence the motion of the medical device. In one embodiment, the motion compensation function may comprise a time-varying vector displacement. Thus, for a given spatial position, for a given time, the function defines a vector displacement (and potentially rotation) by which the P&O of the medical device will need to be corrected so as to match its corresponding value at a time when the image was taken.
In still further embodiments, a plurality of localization sensors are deployed, where the compensation function is a weighted summation of the individual displacement vectors respectively attributed to the movements detected by the plurality of localization sensors. A weighting factor associated with each input may correspond to the correlation level observed for that input relative to the motion of the medical device.
These and other benefits, features, and capabilities are provided according to the structures, systems, and methods depicted, described and claimed herein.
Referring now to the drawings wherein like reference numerals are used to identify identical components in the various views,
Before proceeding to a detailed description keyed to the drawings, a general overview concerning motion compensation will be set forth. As a starting point, there is a desire to reduce a patient's exposure to x-rays, such as may be used in fluoroscopy. It is therefore desirable to be able to use, and reuse to the greatest extent possible, an image (or a sequence of images defining a cine loop) of a region of interest acquired in the past. This will reduce the need for continuous exposure or subsequent additional exposures for the purpose of acquiring updated imaging. Navigation of a medical instrument using the previously-acquired image or cine-loop is made possible by ascertaining the position and orientation (P&O) of the instrument and then superimposing a projection of that instrument's P&O onto the image. A problem arises over time, however, because both the patient as well as his or her internal organs can move (e.g., beating heart), changing positions relative to the time at which the image was taken. Absent compensation for these varying types of motion, the P&O readings reflecting the real time position of the medical instrument would be inaccurately represented on the image (i.e., the representation could be superimposed in the “wrong” location on the image). U.S. Pat. No. 7,386,339 referred to in the Background discloses motion compensation for patient movements and respiration-induced movements by providing a patient reference sensor (PRS). By interpreting P&O readings that track the motion of a catheter relative to the P&O readings of the PRS, a certain type of motion compensation can be achieved. In other words, the movements detected by the PRS shift the coordinate system relative to the coordinate system in which the two dimensional images were acquired. However, as also described in the Background, the PRS P&O readings may have little or no correlation to the movements of an internal moving organ.
With continued reference to
The main control 12, in a computer-implemented embodiment, is programmed to perform a plurality of functions, including those shown in block form in
The image database 18 is configured to store image information of relating to the patient's body, including the moving region of interest, and which may comprise (1) one or more two-dimensional still images acquired at respective, individual times in the past; (2) a plurality of related two-dimensional images obtained in real-time from an image acquisition device (e.g., fluoroscopic images from an x-ray imaging apparatus, such as that shown in exemplary fashion in
The MPS system 20 is configured to acquire positioning (localization) data (i.e., position and orientation—P&O) of one or more MPS sensors. The P&O may be expressed as a position (i.e., a coordinate in three axes X, Y and Z) and orientation (i.e., an azimuth and elevation) of the magnetic field sensor in the magnetic field relative to a magnetic field generator(s)/transmitter(s).
The internal MPS position reference sensor 241 is associated with a moving region of interest (ROI) in the body, which may be a moving organ, and more specifically may be the heart and/or chambers or portions thereof (e.g., atria). The internal position reference sensor 241 is associated with the ROI in such a way that it will move together with the moving ROI, and thus fairly indicate the motion of the region of interest. Generally, associating the sensor 241 with the region of interest (ROI) may be done in any one or more ways: (1) placing the sensor 241, or an interventional device like a catheter carrying the sensor 241, in an anatomical area where it is held by the anatomy itself, for example, a catheter that has been maneuvered in a tubular organ like the coronary sinus; (2) fixing the sensor 241, or an interventional device like a catheter carrying the sensor, to the anatomy in the region of interest using a fixation mechanism, active or passive, for the duration of the procedure; (3) holding the sensor 241, or an interventional device like a catheter carrying the sensor, in steady contact with the anatomy in the region of interest; and (4) placing sensor 241 (or interventional device carrying the sensor) in a non-MPS-enabled device that is in turn affixed to the anatomy in the region of interest. As to approach (2), where the region of interest is the heart, an example may include placing the sensor 241 epicardially in the surface of the heart. As to approach (3), an example may include associating the sensor 241 with a catheter that is maneuvered into steady contact with the heart interior. As to approach (4), an example may include placing an MPS-enabled guidewire (having the sensor 241) in the lumen of a pacing lead that is in turn affixed to the tissue of a heart chamber.
One or more additional, optional internal position sensors may be provided, for example, as shown by sensor 242. The additional one or more sensors 242 may be associated with either or both of the (1) the moving region of interest; or (2) the medical device 26. The additional sensors 242, are configured to provide additional data points (P&O readings) with respect to either the moving region of interest or medical device, as the case may be, thereby providing addition information concerning their respective motions over time.
The patient reference sensor (PRS) 243 is configured to provide a stable, positional reference of the patient's body so as to allow motion compensation for gross patient body movements and/or respiration-induced movements, as described above. The PRS 243 may be attached to the patient's manubrium sternum, a stable place on the chest, or other location that is relatively positionally stable.
In a magnetic field-based embodiment, the P&O may be based on capturing and processing the signals received from the magnetic field sensor while in the presence of a controlled low-strength AC magnetic field. Accordingly, the internal sensors may each comprise one or more magnetic field detection coil(s), and it should be understood that variations as to the number of coils, their geometries, spatial relationships, the existence or absence of cores and the like are possible. From an electromagnetic perspective—all sensors are created equal: voltage is induced on a coil residing in a changing magnetic field, as contemplated here. The sensors 24 are thus configured to detect one or more characteristics of the magnetic field(s) in which they are disposed and generate an indicative signal, which is further processed to obtain the P&O thereof. For one example of a sensor, see U.S. Pat. No. 7,197,354 entitled SYSTEM FOR DETERMINING THE POSITION AND ORIENTATION OF A CATHETER issued to Sobe, hereby incorporated by reference in its entirety.
The electro-cardiogram (ECG) monitor 22 is configured to continuously detect an electrical timing signal of the heart organ through the use of a plurality of ECG electrodes (not shown), which may be externally-affixed to the outside of a patient's body. The timing signal generally corresponds to the particular phase of the cardiac cycle, among other things. The ECG signal may be used by the main control 12 for ECG synchronized play-back of a previously captured sequences of images (cine loop). The ECG monitor 22 and ECG-electrodes may comprise conventional components.
The main control 12, as configured by way of super-imposing function block 32, includes the capability of producing (and superimposing) a projection of the real-time location information (P&O) of a medical device on previously recorded x-ray images or in the case of cine-loops (CL), onto each image in the sequence. In addition, with the availability of the ECG signal and a PRS position signal, the main control 12 can replay a cine loop in an ECG synchronized and respiration-induced motion compensated manner. In a specific case of ECG synchronizing playback of a cine loop of the heart, the sequence is replayed in concordance with a real-time ECG signal (cardiac phase) of the patient. The main control 12 may also be configured to include a respiration compensation algorithm configured to learn the motion induced by the patient's respiration, based on P&O readings from the PRS. The main control 12 then calculates a respiration correction factor to apply to P&O measurements that are to be projected onto a sequence of cine-loop images. The PRS position signal allows for motion compensation for any patient's body movements, as the medical device's position (i.e., P&O measurement) may preferably be taken relative to the P&O measurements from the PRS.
However, as noted above, there are situations where there is very little or no correlation between the internal motion of the region of interest and the external signals (i.e., ECG signals and PRS signal) conventionally used for motion compensation. For example, in the case of atrial fibrillation, the ECG signal cannot serve as a predictor or correlation input for the motion of the atria.
Accordingly, one or more of the internal (i.e., inside the body) position reference sensors (e.g., sensor 24i) are located in the vicinity of the region of interest, or are otherwise associated with the region of interest (e.g., affixed) such that the internal MPS reference sensor moves together with the region of interest over time. As the region of interest moves, the MPS system 20 acquire a series of location (i.e., position and orientation) readings from the sensor. The motion compensation function block 28 (
In
In
In
The series of detected first, second and third positions 1041, 1042 and 1043 of the sensor 241 acquired by the MPS system 20 over time defines not only the motion of the sensor itself but also defines the motion of the region of interest. The motion compensation function block 28 may determine the motion of the region of interest 100 directly in accordance with the motion of the sensor 241. This same motion can be assumed to influence the motion of the medical device 26, provided predefined criteria are met.
The criteria include verifying that an adequate level of correlation exists between the motion of the medical device 26 and the motion of the region of interest 100. One approach to verifying correlation is to compare the respective motions relative to a common time-line. For example, over some time interval, the system 10 may track the motion of the device 26, as indicated by the detected P&O's 1061, 1062 and 1063 shown in
The system 10 may additionally verify that a minimum level of correlation exists between the motion of the device 26 and the other compensation signals described above (i.e., the ECG signal(s) as well as the PRS signal). If there is only poor correlation between the motion of the device 26 and these compensation signals then compensation will not be performed at all. When motion correlation has been verified, the assumption that the motion of the region of interest will influence the motion of the device 26 can be relied on. After correlation has been verified, the MPS system 20 is then enabled to provide motion compensation.
Generate Motion Compensation Function. The MPS system 20 will generate data adequate to track the motion of a moving region of interest over time (e.g., via the internal sensor 241) and allow the compensation function block 28 to generate a time varying motion compensation function. Just as the detected movements of the PRS allows shifting of the coordinate system (as described in U.S. Pat. No. 7,386,339), the motion of the internal sensor (e.g., sensor 241) provides data adequate to implement a similar compensation function. For example in
The compensation function produced based on the motion of the sensor 241 is a time-varying spatial function which accounts for the motion of the region of interest between a first time (at which the image was acquired) and a second time (at which the P&O of the device was measured). Assume that a two-dimensional image was acquired at a time when the region interest was in the first activity state 1001 (i.e.,
In sum, for a given spatial position (measured P&O) of the device 26 for a given time, the compensation function will constitute a vector displacement (and potentially rotation) by which the measured P&O of the MPS-enabled device 26 has to be corrected to match a given time in the past (i.e., at which the image was acquired). The displacement vector may be weighted in accordance with a weighting factor, which in turn may be calculated based on the calculated correlation level described above. Motion compensation approaches may be used as disclosed in U.S. Pat. No. 7,386,339 referred to above as well U.S. Pat. No. 7,343,195 (application Ser. No. 09/949,160 filed Sep. 7, 2001) entitled “METHOD AND APPARATUS FOR REAL TIME QUANTITATIVE THREE-DIMENSIONAL IMAGE RECONSTRUCTION OF A MOVING ORGAN AND INTRA-BODY NAVIGATION” to Strommer et al., the entire disclosures of both U.S. patents being hereby incorporated by reference herein.
P&O Correction. The P&O correction block 30 is configured to correct the P&O reading obtained at the given time using the compensation function. P&O correction function 30 adjusts the measured P&O of the medical device in accordance with the calculated displacement vector (and potentially rotation) described above.
Projection. Finally, a projection of the corrected P&O (three-dimensional) is made onto the two-dimensional image, with a representation of the medical device being superimposed on the image (e.g., may be cross-hairs representing the tip of a catheter or other representation). The resulting image may then be displayed on the display 16. One approach for projecting the corrected P&O onto a 2-D image is a direct consequence of the association of the MPS 3D coordinate system with the X-ray 2D coordinate system, as seen by reference to U.S. Pat. Pub. 2006/0058647, application Ser. No. 11/233,420 entitled METHOD AND SYSTEM FOR DELIVERING A MEDICAL DEVICE TO A SELECTED POSITION WITHIN A LUMEN, to Strommer et al., hereby incorporated by reference in its entirety. Once the coordinate systems are co-registered (a process that may be referred to a magnetic-optical calibration, which may be performed at the installation of the MPS system 20, as noted above), the coordinates of any 3D object (e.g., sensor, landmark or other artifacts) which needs to be displayed on a 2D image are multiplied by a coordinate transformation matrix that computes the corresponding 2D coordinates on the displayed image. This approach is exemplary only and not limiting in nature.
Multiple Inputs, Internal Sensors. In embodiments where the external PRS and/or additional internal position reference sensor inputs are used for motion compensation, the compensation function block 28 implements a composite motion compensation function that is formed by the summation of individual motion compensation contributions, i.e., the individual displacement vectors and (potentially rotations) attributable to each motion/sensor input (provided that correlation requirements are met, as described above). For example, additional sensors may be located on the medical device 26 or other medical devices and/or tools, for example, another MPS sensor disposed on a catheter, guide-wire, etc. The P&O readings from additional internal sensors may reveal other movements or other aspects to the movements of the region of interest and/or the medical device.
The composite compensation function, for example, may include a number of terms where each term corresponds to an input, i.e., one term being provided with respect to the PRS, another term being provided with respect to the internal sensor 241, still another term being provided with respect to an additional sensor 242, and so on. In another embodiment, the inputs from the PRS, the ECG signals and the one or more internal MPS reference sensors may be used in combination to provide for robust motion compensation. In these embodiments, the individual inputs are weighted by a respective weighting factor to form a composite motion compensation function. The respective level of correlation is a principal factor according to which each weighting factor is determined. The weighted function can be depicted as weighted vector summation of the compensation function vectors.
MPS system 110 includes a location and orientation processor 150, a transmitter interface 152, a plurality of look-up table units 1541, 1542 and 1543, a plurality of digital to analog converters (DAC) 1561, 1562 and 1563, an amplifier 158, a transmitter 160, a plurality of MPS sensors 1621, 1622, 1623 and 162N, a plurality of analog to digital converters (ADC) 1641, 1642, 1643 and 164N and a sensor interface 166.
Transmitter interface 152 is connected to location and orientation processor 150 and to look-up table units 1541, 1542 and 1543. DAC units 1561, 1562 and 1563 are connected to a respective one of look-up table units 1541, 1542 and 1543 and to amplifier 158. Amplifier 158 is further connected to transmitter 160. Transmitter 160 is also marked TX. MPS sensors 1621, 1622, 1623 and 162N are further marked RX1, RX2, RX3 and RXN, respectively. Analog to digital converters (ADC) 1641, 1642, 1643 and 164N are respectively connected to sensors 1621, 1622, 1623 and 162N and to sensor interface 166. Sensor interface 166 is further connected to location and orientation processor 150.
Each of look-up table units 1541, 1542 and 1543 produces a cyclic sequence of numbers and provides it to the respective DAC unit 1561, 1562 and 1563, which in turn translates it to a respective analog signal. Each of the analog signals is respective of a different spatial axis. In the present example, look-up table 1541 and DAC unit 1561 produce a signal for the X axis, look-up table 1542 and DAC unit 1562 produce a signal for the Y axis and look-up table 1543 and DAC unit 1563 produce a signal for the Z axis.
DAC units 1561, 1562 and 1563 provide their respective analog signals to amplifier 158, which amplifies and provides the amplified signals to transmitter 160. Transmitter 160 provides a multiple axis electromagnetic field, which can be detected by MPS sensors 1621, 1622, 1623 and 162N. Each of MPS sensors 1621, 1622, 1623 and 162N detects an electromagnetic field, produces a respective electrical analog signal and provides it to the respective ADC unit 1641, 1642, 1643 and 164N connected thereto. Each of the ADC units 1641, 1642, 1643 and 164N digitizes the analog signal fed thereto, converts it to a sequence of numbers and provides it to sensor interface 166, which in turn provides it to location and orientation processor 150. Location and orientation processor 150 analyzes the received sequences of numbers, thereby determining the location and orientation of each of the MPS sensors 1621, 1622, 1623 and 162N. Location and orientation processor 150 further determines distortion events and updates look-up tables 1541, 1542 and 1543, accordingly.
It should be understood that variations and other uses other than for the described imaging embodiments are possible. High-fidelity device positioning as provided through the motion compensation embodiments described herein may be used for alternate purposes such as for placing accurate landmarks (i.e., to serve as navigation references of other devices), for co-registration with other modalities (e.g., Ensite NavX, computed tomography (CT)), as well as for determining when or distinguishing between a “real” motion (i.e., like the actual moving of a catheter by the physician) has occurred versus what seems to be motion but is actually an external event, such as patient motion.
It should be understood that the system 10, particularly main control 12, as described above may include conventional processing apparatus known in the art, capable of executing pre-programmed instructions stored in an associated memory, all performing in accordance with the functionality described herein. It is contemplated that the methods described herein, including without limitation the method steps of embodiments of the invention, will be programmed in a preferred embodiment, with the resulting software being stored in an associated memory and where so described, may also constitute the means for performing such methods. Implementation of the invention, in software, in view of the foregoing enabling description, would require no more than routine application of programming skills by one of ordinary skill in the art. Such a system may further be of the type having both ROM, RAM, a combination of non-volatile and volatile (modifiable) memory so that the software can be stored and yet allow storage and processing of dynamically produced data and/or signals.
Although numerous embodiments of this invention have been described above with a certain degree of particularity, those skilled in the art could make numerous alterations to the disclosed embodiments without departing from the spirit or scope of this invention. All directional references (e.g., plus, minus, upper, lower, upward, downward, left, right, leftward, rightward, top, bottom, above, below, vertical, horizontal, clockwise, and counterclockwise) are only used for identification purposes to aid the reader's understanding of the present invention, and do not create limitations, particularly as to the position, orientation, or use of the invention. Joinder references (e.g., attached, coupled, connected, and the like) are to be construed broadly and may include intermediate members between a connection of elements and relative movement between elements. As such, joinder references do not necessarily infer that two elements are directly connected and in fixed relation to each other. It is intended that all matter contained in the above description or shown in the accompanying drawings shall be interpreted as illustrative only and not limiting. Changes in detail or structure may be made without departing from the spirit of the invention as defined in the appended claims.
This application is a divisional of U.S. application Ser. No. 16/059,557, filed 9 Aug. 2018 (the '557 application), now pending, which is a continuation of U.S. application Ser. No. 12/650,932, filed 31 Dec. 2009 (the '932 application), now U.S. Pat. No. 10,069,668. The '557 application and the '932 application are both hereby incorporated by reference as though fully set forth herein.
Number | Date | Country | |
---|---|---|---|
Parent | 16059557 | Aug 2018 | US |
Child | 17141777 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 12650932 | Dec 2009 | US |
Child | 16059557 | US |