Source: https://patents.google.com/patent/US9414770B2/en
Timestamp: 2018-04-27 07:02:07
Document Index: 329666542

Matched Legal Cases: ['art 38', 'art 38', 'art 100', 'art 100', 'art 100', 'art 38', 'art 100', 'Application No. 11195867']

US9414770B2 - Respiratory effect reduction in catheter position sensing - Google Patents
Respiratory effect reduction in catheter position sensing Download PDF
US9414770B2
US9414770B2 US12980488 US98048810A US9414770B2 US 9414770 B2 US9414770 B2 US 9414770B2 US 12980488 US12980488 US 12980488 US 98048810 A US98048810 A US 98048810A US 9414770 B2 US9414770 B2 US 9414770B2
US12980488
US20120172712A1 (en )
A method, including positioning body-electrodes in galvanic contact with a body of a patient and locating a probe in the body of the patient. The method further includes tracking positions of the probe during respiration of the patient and determining indications related to impedances between the body-electrodes during the respiration. The method also includes calculating a function relating the positions of the probe to the indications, and applying the function to identify end-expirium points of the respiration based on subsequent indications related to the impedances.
The present invention relates generally to sensing the position of an object placed within a living body, and specifically to detection and compensation for respiration effects.
A wide range of medical procedures involve placing objects, such as sensors, tubes, catheters, dispensing devices, and implants, within the body. Real-time imaging methods are often used to assist doctors in visualizing the object and its surroundings during these procedures. In most situations, however, real-time three-dimensional imaging is not possible or desirable. Instead, systems for obtaining real-time spatial coordinates of the internal object are often utilized. In these systems compensation for effects of respiration may be required.
U.S. Patent Application 2009/0030307, to Govari et al., whose disclosure is incorporated herein by reference, describes a method for position tracking that includes placing an internal reference probe in a reference location within a heart of a subject, and collecting and processing location coordinates of the probe during one or more respiratory cycles so as to define a range of the location coordinates corresponding to the reference location.
U.S. Pat. No. 6,711,429, to Gilboa et al., whose disclosure is incorporated herein by reference, describes a system and method of displaying at least one point-of-interest of a body during an intra-body medical procedure. The display may be effected by monitoring and displaying the catheter's location throughout a respiratory cycle and also averaging its location over at least one respiratory cycle.
US Patent Application 2009/0182224, to Shmarak et al., whose disclosure is incorporated herein by reference, describes apparatus for generating an organ timing signal relating to an inspected organ within the body of a patient. The disclosure refers to reconstructing a respiratory trajectory.
locating a probe in the body of the patient;
tracking positions of the probe during respiration of the patient;
determining indications related to impedances between the body-electrodes during the respiration;
calculating a function relating the positions of the probe to the indications; and
applying the function to identify end-expirium points of the respiration based on subsequent indications related to the impedances.
Typically, identifying the end-expirium points includes evaluating that the function is less than a predetermined threshold.
In a disclosed embodiment, tracking the positions of the probe includes determining the positions using an electromagnetic tracking system. Alternatively or additionally, tracking the positions of the probe includes computing the positions in response to currents between the probe and the body-electrodes.
In another disclosed embodiment, determining the indications includes measuring respective currents between the probe and the body-electrodes. Alternatively or additionally, determining the indications includes formulating a matrix relating the impedances between the body-electrodes to the positions of the probe.
Formulating the matrix may include determining a direction of the respiration. The method may further include determining the end-experium points for the respiration in response to the matrix and the direction. The method may yet further include determining the end-experium points for the respiration in response to a rate of change of the indications, the matrix, and the direction.
body-electrodes which are positioned in galvanic contact with a body of a patient;
a probe located in the body of the patient; and
track positions of the probe during respiration of the patient,
determine indications related to impedances between the body-electrodes during the respiration,
calculate a function relating the positions of the probe to the indications, and
apply the function to identify end-expirium points of the respiration based on subsequent indications related to the impedances.
FIG. 1A is a schematic, pictorial illustration of a position sensing system, utilizing a hybrid catheter probe, and FIG. 1B is a schematic detailed view showing the distal end of the hybrid catheter, according to an embodiment of the present invention;
FIG. 2A is a flow chart schematically illustrating a process for operating the position sensing system, and
FIG. 2B is a simplified block diagram of the system, according to an embodiment of the present invention;
FIG. 4 is a schematic illustration of an active current location (ACL) patch circuit, according to an embodiment of the present invention;
FIG. 5 shows exemplary graphs of filters used by a filtration module of the position sensing system, according to an embodiment of the present invention; and
FIG. 6 shows schematic graphs of a projected respiration indicator, and parameters derived from the projected indicator, vs. time, according to an embodiment of the present invention.
An embodiment of the present invention provides a method for determining the location of a probe within the body of a patient, typically a subject undergoing a medical procedure. The method compensates for errors in the measured location that are caused by respiration of the patient.
In order to compensate for the respiration, a probe is located in the body of the patient, and body-electrodes are positioned in galvanic contact with the patient's body. A processor tracks positions of the probe during respiration of the patient. In addition, the processor measures currents between an electrode of the probe and the body-electrodes. From the currents the processor derives indications of impedances between the body-electrodes, as respiration indicators, during the respiration. The processor generates an optimal correlation of a parameterized function of the indicators with the probe positions. The parameterized function may be expressed as a matrix correlating the positions and the respiration indicators.
Once the matrix has been determined, it may be used to determine one or more end-experium times in the respiration cycle.
The following description assumes by way of example, and at least in part, an embodiment having two position-determining coordinate sub-systems. A first, electromagnetic (EM), coordinate sub-system determines a position of a catheter probe in a subject by measuring magnetic fields at the probe. A second, current, coordinate sub-system determines the position of the probe by measuring currents from the probe to different electrodes or patches on the subject. The second sub-system is also herein termed the advanced current location (ACL) sub-system. The two coordinate sub-systems are registered with each other, as described hereinbelow, by using a hybrid catheter probe which is able to measure the magnetic fields for the EM sub-system and to act as a source for the currents of the ACL sub-system.
However, while use of the two coordinate sub-systems may be advantageous, it is not a requirement for embodiments of the present invention. Rather, embodiments of the present invention only require one coordinate system such as the ACL sub-system described herein.
Aspects of a tracking system similar to the ACL sub-system described herein are described in U.S. Patent Application 2010/0079158 to Meir Bar-Tal et al., which is assigned to the assignee of the present application and which is incorporated herein by reference.
FIG. 1A is a schematic, pictorial illustration of a position sensing system 36, utilizing a hybrid catheter probe 20, and FIG. 1B is a schematic detailed view showing the distal end of the hybrid catheter, according to an embodiment of the present invention. The hybrid catheter may also be referred to herein as a mapping-catheter. A medical professional 56 is assumed to operate system 36.
By way of example, except where otherwise stated in the description hereinbelow, mapping-catheter 20 is assumed to be used in an invasive procedure within a chamber of a heart 38 of a patient 40, also referred to herein as subject 40. Alternatively, position sensing system 36 may be used with probes similar to catheter 20 in other body cavities. Subject 40 is placed in a magnetic field generated, for example, by positioning under the subject a location pad 43 containing magnetic field generator coils 42. The magnetic fields generated by coils 42 generate electrical signals in coils 24, 26 and 28 of an electromagnetic (EM) sensor 22 located at the distal end of catheter 20. The electrical signals are conveyed to a control unit 44, which analyzes the signals so as to determine the coordinates of the position and of the orientation of catheter 20. Alternatively, the coils in magnetic field sensor 22 may be driven to generate magnetic fields, which are detected by coils 42.
In addition, in system 36 there may be other non-hybrid catheters comprising one or more electrodes, similar to electrodes 30, 32, and 34, but not comprising a sensor such as sensor 22. Non-hybrid catheters are herein also referred to as investigation-catheters, and the electrodes of the investigation-catheters are also referred to as investigation-catheter-conductive electrodes. System 36 is able to track these investigation-catheters. By way of example, one such non-hybrid catheter 21 is shown in FIG. 1A.
System 36 may also include a reference position sensor, such as an internally-placed catheter, inserted into a moving organ of body 40, herein assumed to be heart 38, and maintained in a substantially fixed position relative to the moving organ. Herein the reference sensor is assumed to comprise a coronary sinus reference catheter (CSRC) 27, and is also referred to herein as reference catheter 27. Catheter 27 may be either a hybrid or a non-hybrid catheter.
Typically, system 36 includes other elements, which are not shown in the figures for the sake of simplicity, and which are referred to as necessary in the following description. For example, system 36 may include an ECG monitor, coupled to receive signals from one or more body surface electrodes, so as to provide an ECG synchronization signal to control unit 44. Typically, system 36 also includes an ablation system operated by control unit 44.
The configuration of FIG. 1A is an example configuration, which is chosen purely for the sake of conceptual clarity. In alternative embodiments, any other suitable configuration can also be used. Typically, processor 46 comprises a general-purpose processor, which is programmed in software to carry out the functions described herein. The software may be downloaded to the processor in electronic form, over a network, for example, or it may, alternatively or additionally, be provided and/or stored on non-transitory tangible media, such as magnetic, optical, or electronic memory.
FIG. 2A is a flow chart 100 schematically illustrating a process for operating system 36 during performance of a procedure on a subject, and FIG. 2B is a simplified block diagram of the system, according to an embodiment of the present invention.
In a reference frame correlation step 102, coordinates measured in an EM reference frame and in an active current location (ACL) reference frame are correlated. An EM tracker sub-system 115 generates measurements in the EM reference frame; an ACL tracker sub-system 117 generates measurements in the ACL frame, also herein termed the Body Coordinate System. The EM tracker sub-system measures locations using the electromagnetic fields generated by coils 24, 26, and 28. The ACL tracker sub-system measures locations using currents through ACL patches 60P.
Except where otherwise indicated, the following steps of the flowchart are performed in intermediate processing modules 119, which comprise body coordinate system module 119A, correlation analysis module 119B, patch current calibration module 119C, current projection module 119D, end-expirium detection module 119E, patch effective area compensation module 119F, and filtration module 119G.
In an ACL patch calibration step 106, processor 46, generates currents in the ACL patches to determine differences in individual ACL patch impedances. The differences in the impedances affect the currents in the ACL patches that are measured by the processor.
Steps 102 and 106 are typically only implemented at the beginning of operation of system 36. The remaining steps of the flow chart are typically performed on a continuing basis during operation of the system.
In a patch compensation step 108 processor 46 compensates for changes in the ACL patches effective area. The changes are typically caused by factors such as change of conductivity of the patch, usually because of sweating, and partial peeling of the patch from the skin of the patient. Processor 46 uses currents similar to those generated in step 106 to determine compensation factors.
In a current projection step 110, the processor measures the currents in the ACL patches that are generated by currents injected into catheters being tracked, and applies the adjustments determined in steps 106 and 108 to the currents. The adjusted currents are used to determine preliminary positions of the catheters being tracked, as well as to determine respiration indicators. The respiration indicators are indications, based on impedance measurements, that are representative of a state of respiration of the subject undergoing the procedure.
In a filtration step 112, the time-varying positions and the time-varying respiration indicators determined in step 110 are filtered to permit only components having the very low frequencies associated with respiration to pass. Higher frequency components, such as those associated with the beating heart, are blocked. In filtration step 112, rates of change of the respiration indicators, also herein termed respiration indicator velocities, are also calculated. The filtered positions, indicators, and velocities are stored for future analysis.
After filtration step 112, the process proceeds in one of two paths. A first, training, path 114 is followed during initialization of the process. A second, operational path 116 is followed once initialization has been completed.
In training path 114, in an analysis step 118, filtered time-varying respiration indicators are correlated with filtered time-varying positions of a “training” catheter, i.e., one that is positioned in a fixed position in the heart for the duration of the training. A combination of the different indicators, which provides the best correlation with the position, is generated in the form of a correlation matrix. In addition, the respiration indicators are analyzed to determine a “most prominent direction” vector for the respiration motion.
The training, as shown by a comparison 120, typically continues until results generated by the correlation matrix are within acceptable preset limits. The comparison checks that positions generated by the correlation matrix are within preset limits of actual positions. Once the positions are within the preset limits, the flow chart transfers to operational path 116.
In a detection step 122 in the operational path, the correlation matrix and the respiration direction vector are used to analyze further incoming filtered respiration indicators (generated in step 112). The analysis estimates a gating point at which end-expirium occurs. The analysis may also include use of a threshold by processor 46 to determine if the estimation of the end-expirium point is to be classified as provisional or as final.
In a position step 124, processor 46 recalls stored training catheter positions, corresponding to the respiratory trajectory of the catheter. The trajectory is gated to the end-expirium point, and the gated trajectory is used to estimate the respiratory movement of a catheter, typically a catheter other than the training catheter. The respiratory position that is estimated in this way is subtracted from actual catheter position coordinates, to give catheter position measurements that are invariant over the respiratory cycle.
The following description explains each of the steps of flow chart 100 in detail.
In body coordinate system module 119A, processor 46 applies the relationship in performing reference frame correlation step 102 of flowchart 100. As stated above, system 36 comprises two tracking sub-systems: EM tracker sub-system 115 using sensors such as sensor 22, and ACL tracker sub-system 117 using currents through patches 60P. Each sub-system operates in a respective frame of reference. The EM tracker sub-system operates in an EM frame of reference that is generally fixed with respect to pad 43. The ACL tracker sub-system operates in an ACL frame of reference, the body coordinate system (BCS), that is assumed to be generally fixed with respect to patches 80R. Patches 80R enable measurements made in one of the sub-systems to be converted to the other sub-system. Reference patches 80R are attached to the back of subject 40, so that any motion of the subject with respect to pad 43 is reflected in signal changes in the EM sensors of the reference patches.
Initially, processor 46 analyzes signals from the EM sensors on reference patches 80R, to determine an initial frame of reference for the BCS. Typically, after the initial frame has been determined, the processor analyzes signals from the EM sensors periodically, to determine changes in the location and orientation of the BCS frame of reference. The processor is able to detect if system parameters have changed beyond those expected, and in this event may return to process followed initially.
Initially, the processor accumulates the location of patches 80R in LP coordinates, i.e., coordinates measured relative to location pad (LP) 43, for a time patchInitTime, typically approximately 1 sec.
P -> ⁢ 1 ~ = 1 N ⁢ ∑ i N ⁢ P -> ⁢ 1 i ⁢ ⁢ P -> ⁢ 1 STD = ∑ i N ⁢ 1 N ⁢ ( P -> ⁢ 1 i - P -> ⁢ 1 ~ ) 2 , ( 1 )
is the mean value of {right arrow over (P)}li for each patch l, and
STD is the standard deviation of
Providing the value of each
STD is less than a preset figure, typically approximately 1 mm, the calibration is accepted, in which case the mean {tilde over ({right arrow over (P)})} of all the means is set as the origin the BCS:
P -> ~ = 1 3 ⁢ ∑ 1 3 ⁢ P -> ⁢ 1 ~ ( 2 )
{right arrow over (P)}linit=
−{tilde over ({right arrow over (P)})} (3)
During operation of system 36, patches 80R may move, as exemplified by patches 80′, 82′, and 84′, and processor 46 measures the new positions of the patches periodically, typically with a period of the order of one second. Embodiments of the present invention assume that the axes defined in the calibration phase move as an approximately rigid body, and processor 46 determines the translation and rotation of the axes from the new patch 80R positions during the tracking phase. Prior to the determinations, the new patch positions are filtered to reduce noise, the filtering typically comprising a low pass filter of the type:
The vector p is calculated in position step 124.
Reference is now made to FIG. 4, which is a schematic illustration of an ACL patch circuit, according to an embodiment of the present invention.
All ACL patches have generally similar circuits. Each ACL patch i comprises a defibrillation protection circuit 152 and an ablation protection circuit 154. The two circuits are connected in series between the patch and ground. In FIG. 4, and for the analysis below, for each patch i,
Iij is a current measured through patch i at frequency fi
The processor applies the two procedures referenced above, i.e., the patch impedance calibration procedure, and the patch effective area compensation procedure, typically while the training path, described above, is being followed. The processor may also apply, as necessary, the two procedures while the operational path is being followed.
q ij = q ^ i ⁡ ( f j ) = [ V ij I ij ] ^ ( 8 )
V ~ j ≡ ∑ i ⁢ V ij , and ⁢ ⁢ E ij ≡ δ ij ⁢ E j ( 9 )
Applying Ohm's law and equation (10) to the complete circuit of FIG. 4 gives, for a particular patch i:
X ij = E ij + ( q ij + z ij ) ⁢ I ij = δ ij ⁢ E j + ( q ij + z ij ) ⁢ I ij = δ ij ( V ~ j - ∑ k ⁢ q kj ⁢ I kj ) + ( q ij + z ij ) ⁢ I ij ( 11 )
where I is a matrix of patch currents, and X is a matrix of patch voltages.
The patch currents may also be written as a vector s. The negative sign in equation (12) assumes a convention that positive current flows into body 40, and positive currents are also measured flowing out of the body. Alternatively, an equation similar to equation (12), using an impedance matrix Im, may be written relating matrices I and X.
Equation (13b) shows that is affected by the patch resistances. A calibrated current, that does not depend on the patch resistances and thus does not depend on frequency fk, may be defined as:
The description in this section explains patch compensation step 108 (FIG. 2A). In step 108 processor 46 uses patch effective area module 119F to perform a process that compensates for changes in the effective area of ACL patches i. In this section, patches 60P are differentiated by being referred to as patch i and patch j. Some causes of changes to the effective area of the ACL patches are partial peeling of a patch from patient body 40, and skin conductivity changes, typically due to sweating. A patch effective area compensation model assumes that
m ^ ij = m ij ∑ k = 1 , k ≠ j k = N ⁢ m kj = V j / I ij ∑ k = 1 , k ≠ j k = N ⁢ V j / I kj = 1 / I ij ∑ k = 1 , k ≠ j k = N ⁢ 1 / I kj ( 16 )
Processor 46 calculates and stores the values of mid using equation (16) during implementation of the area compensation process.
m ^ ij = C i ⁢ d ij ∑ n , n ≠ j ⁢ C n ⁢ d nj ( 19 )
Assuming that processor 46 performs eigenvector
analysis of the matrix ATA, the processor selects the eigenvector with the smallest eigenvalue. Typically, the values of {circumflex over (m)}ij and dij for matrices A and AT are filtered with a filter similar to equation (4), where the filter is adjusted to have a time constant of the order of 10 seconds. The smallest eigenvector corresponds to normalized values of the 6 areas Ci, herein termed Eai. Typically, processor 46 calculates the values of Eai periodically, with a period that may be set by operator 56. In one embodiment, the period is of the order of 1 second.
The estimated current vector s, derived from equation (13c), gives 6 respective values of the currents Ii in the ACL patches. To compensate for the patches effective area, Eai, processor 46 forms a normalized value of each of the currents:
In i = I i · Ea i ∑ i = 1 i = 6 ⁢ I i · Ea i ≡ ( In ⁢ ⁢ 6 ) ( 21 )
In addition to performing the normalization of equation (23), in current projection step 110 processor also normalizes respiration indicators. The respiration indicators are measures of changes in resistance between ACL patches caused by the respiration of the patient. The respiration indicators are generally similar in properties to the patch effective areas (Ci), and an expression for the indicators is derived using equation (20). From equation (20),
Averaging equation (24) over j, and assuming as an approximation that dij=1(i≠j); dij=0(i=j), we define a respiration indicator RIi for patch i as:
RI i = ∑ j ⁢ m ^ ij ⁢ ∑ n ⁢ R n ⁢ d nj ^ = ∑ n ⁢ R n ( ∑ n ⁢ m ij ^ ⁢ d nj ^ ) ( 25 )
Assuming RIi is written as an N-dimensional vector {right arrow over (RI)}i,
may be written as an N×N matrix M, with 0 at the diagonals, and
may be written as a matrix D, given by equation (26):
In this case processor 46 iteratively estimates a new vector RIi,t from a previous estimate RIi,t-1 and a new measurement
according to the two step process:
RI i , t _ = M T · D · RI i , t - 1 , then ( 27 ) RI i , t := RI i , t _ ∑ i ⁢ RI i , t
Vector RIi,t is six-dimensional, but has only five degrees of freedom (as in the currents of equation (21)). As shown in equation (28) below, the values of RIi derived from equation (27) are further converted in the same manner as the currents of equation (21), by multiplying by matrix J, to give a five-dimensional respiration indicator vector.
FIG. 5 shows exemplary graphs of filters used by filtration module 119G, according to an embodiment of the present invention. Filtration step 112, implemented in module 119G, filters both the locations derived from the EM system, i.e., the values of {right arrow over (p)}b as given by equation (5), and the five-dimensional respiration indicator vector (RIi5) given by equation (28). Hereinbelow, for simplicity, the respiration indicator vector (RIi5) is also referred to as (RI), having five elements RI. The filtration typically uses a very low pass filter, with a sharp transition between passing and rejected frequencies, so as to preserve only the respiratory frequencies of the signals, while removing frequencies due to the heart beat. In one embodiment, module 119G comprises a finite element response (FIR) filter which has a 0 dB rejection for frequencies less than 0.2 Hz and a 40 dB or greater rejection for frequencies greater than 0.7 Hz, as shown in a graph 180.
In module 119G values of the respiration indicators produced by the FIR filter are also adjusted so that they are centered around zero (so that there is no “DC” component attached to the RI values). The adjustment may be made using a low pass filter such as is defined by equation (29) to estimate an adjustment value.
y i =ay i-1+(1−α)x i, (29)
where yi, yi-1 are current and previous adjustment values,
xi is the current RI measurement, and
α is a factor between 0 and 1. α is typically selected so that there is an effective time constant of approximately 1 s in determining current adjustment values.
The values of yi from equation (29) are subtracted from the values of RI produced by the FIR filter, to remove the RI DC component. The filtered, adjusted respiration indicator vector is herein termed (RIa).
Module 119G also calculates respiration velocities, i.e., rates of change of the respiration indicators, and filters the velocities. In one embodiment, the differentiation and filtration of the respiration indicators RI is produced by convolving the indicators with a velocity kernel. The velocity kernel may be derived from differentiation of a Gaussian function, and a graph of a typical velocity kernel v is given in a graph 182. The convolution is according to equation (30):
Rv=RI{circle around (x)}v (30)
where Rv is the velocity of respiration indicator RI.
During the initialization mode, i.e., in training path 114, values from filtration module 119G are analyzed in analysis step 118. The analysis uses correlation analysis module 119B. At the start of initialization, typically at the beginning of a medical procedure, professional 56 positions a reference sensor, herein assumed to be CSRC 27, in a relatively fixed position in heart 38. The sensor is typically left in this position for two or more respiration cycles, i.e. for a time period of approximately 30 s, herein termed the training period.
During the training period, processor 46 saves position values generated by CSRC 27, and the filtered, adjusted respiration indicators RIa. Since the correlations performed by module 119B are with regard to motion due to respiration, measurements are performed with respect to a mean of the reference sensor's position, determined according to equation (31):
r → i = p → i - 1 K ⁢ ∑ 1 K ⁢ p → i ( 31 )
{right arrow over (p)}i is the sensor position for sample i, and
{right arrow over (r)}i is the adjusted sensor position for sample i, and
K is the total number of samples in the training period.
To avoid using position and respiration indicator values that may be very distant from a central value, e.g., they may have been obtained during sudden motion of the sensor, the filtered positions are sorted to remove such values. In one embodiment the absolute values of {right arrow over (r)}i are assumed to lie between 0 and 20 mm, and the absolute values are sorted into a histogram with 1 mm width bins. Only positions, and associated respiration indicators, within a user-selected percentile of the histogram are used in further analysis by module 119B. Typically the 90th percentile is selected. The respiration indicator vector produced by the sorting is herein referred to as (RIas).
To prevent numerical instability in the further analysis, the elements of (RIas) are scaled to have values close to unity. The scaling may be accomplished by multiplying each of the elements of (RIas), RIask, by a common factor g, given by equation (32):
g = 1 Max ⁡ ( Abs ⁡ ( RIas → k ) ) ( 32 )
In addition to multiplying the elements of (RIas) by the factor g, a unity element is added to vector (RIas), producing a six-dimensional respiration indicator vector (RIk) having elements given by equation (33):
k ={g{right arrow over (RIas)} k,1} (33)
A correlation between the respiration indicators RIk and the corresponding position vector rk, in the form of a (3 row, 6 column) correlation matrix C, is generated from vector (RIk) and vector rk according to equation (34):
C=({right arrow over (r)} k)·({circumflex over (R)}I k)−1 (34)
where ({right arrow over (r)}k) are stacked values of vectors rk,
({circumflex over (R)}I k) are stacked values of vectors (RIk), and
({circumflex over (R)}I k)−1 is the pseudo inverse of ({circumflex over (R)}I k).
Correlation matrix C effectively parameterizes the different respiration indicators RIk so that the parameterized indicators form an optimal correlation with the position vector rk.
In a final part of analysis step 118, as shown by comparison 120, matrix C is validated by estimating values of rk, {circumflex over (r)}k, using matrix C and corresponding values of RIk, according to equation (35):
{circumflex over (r)}k =C
k (35)
The estimated values, {circumflex over (r)}k, are compared with the measured values of rk stored in memory 47. Providing the comparison gives a difference that is within a preset threshold limit, analysis step 118, and thus training mode 114, is assumed to be complete and correlation matrix C is used in operational mode 116. Typically, a root mean square (RMS) of the difference between the estimated and measured values of the location is less than 1 mm, and a maximum value of the difference is less than 3 mm.
If the difference between the estimated and measured values of rk is not within the preset threshold limit, the training continues.
In addition to determining correlation matrix C, analysis module 119B determines a “most prominent direction” for the respiration motion, as a respiration direction vector {right arrow over (rd)}. In one embodiment vector {right arrow over (rd)} is calculated by stacking the adjusted sensor positions found by equation (31) into a matrix R, as in equation (36):
R = [ x 1 y 1 z 1 x 2 y 2 z 2 ⋮ ⋮ ⋮ x n y n z n ] ( 36 )
Singular value decomposition (SVD) is applied to matrix R, and respiration direction vector {right arrow over (rd)} is assumed to be equal to the eigenvector of the largest eigenvalue.
End Expirium Detection
Once the correlation matrix C and respiration direction vector {right arrow over (rd)} have been determined in analysis step 118, processor 46 transfers the matrix and the respiration direction to end-expirium detection module 119E. After the transfer, processor 46 follows operational path 116 of flow chart 100. In detection step 122, implemented in module 119E, processor 46 determines when end-expirium occurs.
FIG. 6 shows schematic graphs of a projected respiration indicator, and parameters derived from the projected indicator, vs. time, according to an embodiment of the present invention. The graphs illustrate results generated in module 119E.
In step 122, processor 46 calculates, for each set of filtered indicator values derived from filtration module 119G, an indicator projected onto the most prominent respiration direction. The projection is performed according to equation (37):
RIp k ={right arrow over (rd)}·C·
k (37)
where RIpk is the projected indicator at a time k,
and {right arrow over (rd)}, C, and
k are as described above with reference to equations (36), (34), and (33).
A graph 200 shows exemplary values of RIp vs. time.
The processor also calculates a projected indicator velocity according to equation (38):
RIv k ={right arrow over (rd)}·C(3,5)·Rv k (38)
where RIvk is the projected indicator velocity at a time k,
{right arrow over (rd)} is as described above,
C(3,5) is a (3 row by 5 column) matrix derived from matrix C by omitting the last column, and
Rvk is the respiration indicator velocity given by equation (30).
A graph 202 shows exemplary values of RIv vs. time.
Typically, maximum values of RIv, RIvMax, are used to determine boundaries of each respiration cycle.
Typical values of RIvMax are illustrated as points 204 on graph 202. Corresponding points on graph 200 are points 206.
In order to reduce false determinations of RIvMax, processor 46 typically ignores any maximum that is within is of a previous RIvMax value, since such a maximum may be caused by penetration of a heartbeat signal through filtration module 119G.
To further reduce false determinations of RIvMax, processor 46 may calculate a decay value with time for a given RIvMax, and may check that a subsequent RIvMax is greater than that given by the decayed value. The decaying values of RIvMax may be calculated iteratively using an equation of the form:
RIvMaxi=α(RIvMaxi-1−β)+β (39)
In one embodiment α=0.99 and β=0.5. Typically, processor 46 prevents the decaying value from reaching zero, limiting the value to a preset lower limit.
Broken lines 208 illustrate the decaying values of RIvMax applied to each RIvMax. As is illustrated in the graphs, each subsequent point 204 is above the decaying value from the previous point 204.
Once the processor has detected a valid maximum velocity RIvMax the processor analyzes values of RIp (graph 200) (up to those corresponding to the previous RIvMax) to determine a minimum of value of RIp. These minima, shown as points 210, correspond to the end-expirium points.
The end-experium points are points in time, and the values of RIp may be corrected so that all RIp values are greater than or equal to zero. In order to accomplish this correction, the processor typically applies a linear interpolation and a translation (herein, by way of example, an elevation) to the RIp values between adjacent end-experium points, setting both the beginning and end points equal to zero. The segments of a line 212 illustrate the line on which the interpolation is based, and a broken line 214 illustrates the translated, interpolated values of RIp, herein termed RIe values.
The interpolation described above relies on interpolation between two known end-experium points. In embodiments of the present invention, after a latest end-expirium point has been determined, the processor may estimate values of RIe at a time i,
ei, until a next end-experium point is determined. The estimation typically performs a similar type of translation as that described above, from the latest determined end-expirium point.
ei values may be used to predict a provisional next end-expirium point, by comparing the
ei values with a user-defined threshold Th, i.e., by finding times when the condition
ei<Th is true. A typical value of Th is 2. Typically, the provisional next end-expirium point is assumed to be correct until an actual end-expirium point (a minimum value of RIp) is found. The actual end-expirium point then supersedes the provisional point.
1. Apparatus for determining position and orientation coordinates of a probe within a body of a patient while compensating for respiration of the patient, the apparatus comprising:
(i) a plurality of body-electrodes adapted for positioning in galvanic contact with a body of a patient;
(ii) a probe having at least one conductive electrode and adapted for positioning in the body of the patient; and
(iii) a control unit comprising:
(a) an alternating current driver for measuring current between the at least one conductive electrode and the plurality of body-electrodes;
(b) a voltage driver for measuring impedance between the at least one conductive electrode and the plurality of body-electrodes;
(c) a processor, which is configured to:
(i) track position and orientation coordinates of the probe during respiration of the patient using the at least one conductive electrode and the plurality of body-electrodes,
(ii) determine respiration indications related to impedance measured between the body-electrodes during the respiration,
(iii) calculate a function relating the tracked position and orientation coordinates of the probe to the respiration indications,
(iv) apply the function to identify end-expirium points of the respiration of the patient based on subsequent indications related to the impedance and evaluating that the function is less than a predetermined threshold; and
(v) determine respiratory compensated position and orientation coordinates of the probe using an estimate of respiratory movement of the probe based on the end-expirium points.
2. The apparatus according to claim 1, further comprising an electromagnetic coordinate sub-system, and wherein the processor is further configured to track the positions of the probe by determining the position and orientation coordinates of the probe using the electromagnetic coordinate sub-system.
3. The apparatus according to claim 1, wherein tracking the positions of the probe comprises computing the positions in response to currents between the at least one conductive electrode of the probe and the plurality of body-electrodes.
4. The apparatus according to claim 1, wherein the processor is further configured to determine the respiration indications by measuring respective currents between the at least one conductive electrode of the probe and the plurality of body-electrodes.
5. The apparatus according to claim 1, wherein the processor is further configured to calculate a function by formulating a matrix relating the impedances respiration indications related to impedances between the plurality of body-electrodes to the positions and orientation coordinates of the probe.
6. The apparatus according to claim 5, wherein the processor is further configured to formulate the matrix by determining a direction vector of respiration motion.
7. The apparatus according to claim 6, wherein the processor is further configured to determine the end-expirium points for the respiration motion in response to the matrix and the direction vector.
8. The apparatus according to claim 6, wherein the processor is configured to determine the end-experium points for the respiration based on a rate of change of the respiration indications, the matrix, and the direction vector.
9. An apparatus for determining position and orientation coordinates of a probe within a body of a patient while compensating for respiration of the patient, the apparatus comprising:
(i) obtain position and orientation coordinates measurements of the probe using the at least one conductive electrode and the plurality of body-electrodes,
(ii) track position and orientation coordinates of the probe during respiration of the patient based on the position and orientation coordinates measurements using the at least one conductive electrode and the plurality of body-electrodes,
(iii) determine respiration indications related to impedances between the body-electrodes during respiration of the patient, the respiration indications being representative of a state of respiration of the patient,
(iv) calculate a function relating the tracked position and orientation coordinates of the probe to the respiration indications,
(v) apply the function to identify end-expirium points of the respiration based on subsequent respiration indications related to impedance measured between the at least one conductive electrode and the plurality of body-electrodes and evaluating that the function is less than a predetermined threshold,
(vi) determine compensated position and orientation coordinates measurements of the probe using an estimate of the respiratory movement based on the end-expirium points, and
(vii) determine respiratory compensated position and orientation coordinates of the probe based on the compensated position and orientation coordinates measurements of step (vi).
US12980488 2010-12-29 2010-12-29 Respiratory effect reduction in catheter position sensing Active 2031-11-28 US9414770B2 (en)
US12980488 US9414770B2 (en) 2010-12-29 2010-12-29 Respiratory effect reduction in catheter position sensing
CA 2762298 CA2762298A1 (en) 2010-12-29 2011-12-15 Respiratory effect reduction in catheter position sensing
JP2011288915A JP6045790B2 (en) 2010-12-29 2011-12-28 Reduction of influence of respiration in the detection of the catheter position
EP20110195867 EP2471482B1 (en) 2010-12-29 2011-12-28 Respiratory effect reduction in catheter position sensing
CN 201110462048 CN102525472B (en) 2010-12-29 2011-12-29 Effects of reduced respiratory catheter position sensing,
US15097626 US20160262655A1 (en) 2010-12-29 2016-04-13 Respiratory effect reduction in catheter position sensing
US15097626 Division US20160262655A1 (en) 2010-12-29 2016-04-13 Respiratory effect reduction in catheter position sensing
US20120172712A1 true US20120172712A1 (en) 2012-07-05
US9414770B2 true US9414770B2 (en) 2016-08-16
ID=45470372
US12980488 Active 2031-11-28 US9414770B2 (en) 2010-12-29 2010-12-29 Respiratory effect reduction in catheter position sensing
US15097626 Pending US20160262655A1 (en) 2010-12-29 2016-04-13 Respiratory effect reduction in catheter position sensing
US (2) US9414770B2 (en)
EP (1) EP2471482B1 (en)
JP (1) JP6045790B2 (en)
CN (1) CN102525472B (en)
CA (1) CA2762298A1 (en)
US20160235339A1 (en) 2015-02-13 2016-08-18 Biosense Webster (Israel) Ltd. Training of impedance based location system using registered catheter images
JP2005319256A (en) 2004-05-07 2005-11-17 Kanazawa Inst Of Technology Estimate method of parasympathetic activity
WO2007109406A1 (en) 2006-03-22 2007-09-27 Koninklijke Philips Electronics, N.V. Respiration-gated cardiography
JP2008058657A (en) 2006-08-31 2008-03-13 Nippon Telegr & Teleph Corp <Ntt> Optical fiber cutting tool and cutting method
WO2009128062A2 (en) 2008-04-15 2009-10-22 Navotek Medical Ltd. Hybrid medical device localization system
CN101657153A (en) 2007-03-09 2010-02-24 圣朱德医疗有限公司房颤分公司 System and method for correction of inhomogeneous fields
US20140221803A1 (en) 2008-09-30 2014-08-07 Biosense Webster (Israel), Ltd. Current localization tracker
JPS54127192A (en) * 1978-03-24 1979-10-02 Tokyo Shibaura Electric Co Breathing monitor device
US20110190647A1 (en) 2006-03-22 2011-08-04 Koninklijke Philips Electronics N. V. Respiration-gated cardiography
JP2009530047A (en) 2006-03-22 2009-08-27 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Cardiography which is gated by respiration
CN101199416A (en) 2006-09-06 2008-06-18 韦伯斯特生物官能公司 Correlation of cardiac electrical maps with body surface measurements
Australian Examination Report, Jun. 14, 2014,from co-pending U.S. Appl. No. 12/980,488 filed on Dec. 29, 2010-pending.
Chinese Official Action and Search Report, Sep. 23, 2014, from co-pending U.S. Appl. No. 12/980,488 filed on Dec. 29, 2010-pending.
Extended European Search Report dated Apr. 12, 2012 from related European Application No. 11195867.4.
Extended European Search Report EP 11 19 5867 dated Apr. 12, 2012.
JP Office Action for JP2011-288915 dated Nov. 10, 2015.
CN102525472A (en) 2012-07-04 application
JP6045790B2 (en) 2016-12-14 grant
US20120172712A1 (en) 2012-07-05 application
US20160262655A1 (en) 2016-09-15 application
CN102525472B (en) 2015-09-09 grant
JP2012139506A (en) 2012-07-26 application
EP2471482B1 (en) 2016-10-26 grant
CA2762298A1 (en) 2012-06-29 application
EP2471482A1 (en) 2012-07-04 application
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BAR-TAL, MEIR;REEL/FRAME:026742/0701