systems and methods can be used to determine conduction velocity and generate one or more conduction velocity and/or pattern maps to facilitate identification of arrhythmogenic mechanisms.
|
25. A method, comprising:
defining a region of interest around a given node on a three-dimensional surface of interest, electrophysiological signals being stored, in memory, for each of a plurality of nodes distributed across the surface of interest;
orthogonally projecting the region of interest for the given node from the surface of interest onto a two-dimensional surface tangential to the given node;
estimating a two-dimensional electrical conduction velocity vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the region of interest;
estimating a three-dimensional conduction velocity vector for the given node based on the two-dimensional electrical conduction velocity vector; and
determining a unitary conduction velocity vector field for the given node based on the three-dimensional conduction velocity vector.
1. A method, comprising:
defining a region of interest around a given node on a three-dimensional surface of interest, the defining the region of interest comprising selecting a set of neighboring nodes in the region of interest for the given node according to temporal and spatial criteria with respect to the given node, electrophysiological signals being stored, in memory, for each of a plurality of nodes distributed across the surface of interest;
orthogonally projecting the region of interest for the given node from the surface of interest onto a two-dimensional surface tangential to the given node;
estimating a two-dimensional electrical conduction velocity vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the region of interest; and
estimating a three-dimensional conduction velocity vector for the given node based on the two-dimensional electrical conduction velocity vector.
26. A system comprising:
memory to store machine readable instructions and data, the data comprising electrical data representing electrophysiological signals for a plurality of nodes distributed across a three-dimensional surface of interest; and
one or more processors to access the memory and execute the instructions, the instructions programmed to at least:
determine a region of interest around a given node of the plurality of nodes on the surface of interest;
project the region of interest around the given node from the surface of interest onto a two-dimensional surface tangential to the given node;
estimate a two-dimensional electrical conduction velocity vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the region of interest;
estimate a three-dimensional conduction velocity vector for the given node on the three-dimensional surface of interest based on the two-dimensional electrical conduction velocity vector; and
determine a unitary conduction velocity vector field for the given node from the three-dimensional conduction velocity vector.
14. A system comprising:
memory to store machine readable instructions and data, the data comprising electrical data representing electrophysiological signals for a plurality of nodes distributed across a three-dimensional surface of interest; and
one or more processors to access the memory and execute the instructions, the instructions programmed to at least:
determine a region of interest around a given node of the plurality of nodes on the surface of interest, the region of interest being determined to define a set of neighboring nodes in the region of interest for the given node according to temporal and spatial criteria with respect to the given node;
project the region of interest around the given node from the surface of interest onto a two-dimensional surface tangential to the given node;
estimate a two-dimensional electrical conduction velocity vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the region of interest; and
estimate a three-dimensional conduction velocity vector for the given node on the three-dimensional surface of interest based on the two-dimensional electrical conduction velocity vector.
2. The method of
3. The method of
the delay being compared to a threshold to remove a spatially identified neighboring node from the set of neighboring nodes.
4. The method of
5. The method of
6. The method of
7. The method of
8. The method of
9. The method of
10. The method of
11. The method of
12. The method of
13. One or more non-transitory media comprising instructions, which when executed by a processor, perform the method of
15. The system of
16. The system of
determine a delay between an activation time of the given node the activation time each of the neighboring nodes in the region of interest, wherein the delays are computed independent for the neighboring nodes on the three-dimensional surface of interest; and
remove a spatially identified neighboring node from the set of neighboring nodes based on the delay relative to a threshold,
wherein the threshold is determined for each of the neighboring nodes in the region of interest based on the activation time of the given node, the activation time of the respective neighboring node and a cycle length.
17. The system of
18. The system of
19. The system of
20. The system of
21. The system of
22. The system of
23. The system of
24. The system of
wherein the instructions are further programmed to reconstruct the electrophysiological signals on to each of the plurality of nodes of the surface of interest based on the electrical measurement data and geometry data.
|
This application claims the benefit of priority from U.S. Provisional Patent Application No. 62/503,747, filed May 9, 2017 and entitled CONDUCTION VELOCITY AND PATTERN MAPPING, which is incorporated herein by reference in its entirety.
This disclosure relates to conduction velocity mapping and pattern mapping.
Ventricular (VT) and atrial tachycardias (AT) are some of the most common clinical cardiac tachyarrhythmias. They are characterized by abnormal and fast muscular activity, with an organized electrical activity resulting from altered electrical propagation. Atrial tachycardia (AT) management is challenging due to mechanism complexity and altered electrogram (EGM) voltages, especially after AF ablation. However, post AF ablated ATs are still poorly investigated through electrocardiographic imaging (ECGI) that is currently used in clinics for noninvasive electro-anatomic mapping (EAM) of atrial arrhythmias.
This disclosure relates to conduction velocity mapping and pattern mapping.
As one example, a method includes defining a region of interest around a given node on a three-dimensional surface of interest. Electrophysiological signals are stored, in memory, for each of a plurality of nodes distributed across the surface of interest. The method also includes orthogonally projecting the region of interest for the given node from the surface of interest onto a two-dimensional surface tangential to the given node. The method also includes. The method also includes estimating a two-dimensional conduction velocity vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the region of interest. The method also includes estimating a three-dimensional conduction velocity vector for the given node based on the two-dimensional conduction velocity vector. In another example, one or more non-transitory computer readable media can be programmed with instructions which when executed by one or more processor perform such method.
As yet another example, a system may include memory to store machine readable instructions and data, the data including electrical data representing electrophysiological signals for a plurality of nodes distributed across a three-dimensional surface of interest. One or more processors may access the memory and execute the instructions. The instructions are programmed to determine a region of interest around a given node of the plurality of nodes on the surface of interest. The instructions are further programmed to project the region of interest around the given node from the surface of interest onto a two-dimensional surface tangential to the given node. The instructions are further programmed to estimate a two-dimensional conduction velocity vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the region of interest. The instructions are further programmed to estimate a three-dimensional conduction velocity vector for the given node on the three-dimensional surface of interest based on the two-dimensional conduction velocity vector.
This disclosure provides a noninvasive diagnostic tool that helps to identify arrhythmogenic mechanisms, including atrial tachycardia (AT) and ventricular tachycardia (VT) mechanisms, using conduction velocity and/or associated electrical propagation pattern features. For example, conduction velocity and/or propagation pattern features across a three-dimensional (3D) region of interest (ROI) are determined and displayed in a graphical user interface (GUI), such as an electroanatomic map of a cardiac envelope (e.g., displaying patterns across a hear surface, such as an epicardial surface).
Clinical accuracy of the propagation mapping, the rotation and dispersion features across the 3D ROI using mathematical operators for noninvasive diagnosis of complex ATs and VTs may be improved, such as by displaying one or more graphical maps of 3D conduction velocity (CV) and/or rhythm patterns derived from the CV map. For example, such 3D CV and/or rhythm pattern maps may provide additional information across the 3D ROI that can visualize the effects of treatment (e.g., before and after ablation at one or more sites) and further identify one or more other treatment delivery sites to be targeted subsequently. As a result of the systems and methods disclosed herein, treatment (e.g., ablation) procedural time can be reduced and the overall procedure simplified. The approach herein also affords online use, such as in clinics. Additionally, the results pattern maps can be utilized for detection of specific AT or VT propagation pattern features (e.g., rotor and foci) for automated or semi-automated diagnosis.
In the flow diagrams disclosed herein
The method begins at 12 to include preprocessing of electrical signals. The signals represent electrical activity measured across a surface of a patient's body over one or more time intervals. For example, a plurality of electrodes can be positioned on the patient's body to measure electrical activity (unipolar electrograms) across the body surface. The electrodes thus can measure the electrical signals non-invasively and the signals can stored as electrical data in memory for the preprocessing. As an example, the preprocessing at 12 may include one or more of the following:
At 14, the method includes reconstructing the preprocessed signals to nodes on a cardiac envelope based on electrical data and geometry data (e.g., provided at 12). For example, the cardiac envelope is three-dimensional (3D) surface model that defines a mesh surface of a cardiac envelope (e.g., an epicardial surface). The mesh may be described by geometry data (stored in memory), such as includes a collection of nodes (vertices) that define a plurality of polygons (triangle or other polygons) forming the surface. Each polygon can include three or more vertices, such that adjacent vertices define edges that surround a face of each respective polygon. Therefore, at 14, corresponding electrogram signals are computed for each of the nodes on the surface based on the preprocessed signals derived from the body surface electrical signal measurements and geometry data.
As one example the reconstruction of signals on to the cardiac envelope may include solving the inverse problem (e.g., using the Method of Fundamental Solutions or other approach) to reconstruct unipolar epicardial electrograms on the cardiac envelope (epicardial surface) based on geometry data. For instance the geometry data includes patient atrium and torso 3D meshes (complex atrium geometry, location, orientation and size). The resulting geometry (3D spatial coordinates) and electrical data (reconstructed electrograms) are stored in memory. Examples of methods to compute reconstructed electrograms from measured body surface electrical signals are disclosed in U.S. Pat. Nos. 6,772,004 and 7,983,743, which are incorporated herein by reference. The reconstructed electrical signals across the mesh surface may be stored in memory as electroanatomical data, including electrical data representing the time-varying electrical potential (e.g., unipolar electrograms) estimated for each of the nodes over one or more time intervals and geometry data described nodes across the surface.
At 16, activation time is computed based on the reconstructed electrical signals on the surface envelope. For example, each node n of the mesh is defined by its 3-D Cartesian coordinates (x, y, z) and its activation time is computed based on the reconstructed electrical signals computed at 14. The activation time may be computed for the nodes distributed across one or more regions of interest on the surface, which may be for any sized area up to the entire surface envelope. For example, the activation time is computed from a method using the reconstructed electrograms and delays between neighboring nodes on the mesh, such as disclosed in U.S. Pat. No. 9,610,023, which is incorporated herein by reference.
At 18, a conduction velocity map is generated. The conduction velocity map is generated by computing a velocity vector based on the computed activation time for nodes distributed across one or more regions of interest, such as any sized area up to the entire surface (of the heart). For example, the conduction velocity vector may be estimated for each of the nodes based on spatial and temporal criteria derived from a set of neighboring nodes within a region of interest. As an example, the neighboring nodes can be selected to be within the region of a given node if both spatial and temporal criteria are satisfied, such as disclosed herein.
By way of example, the region of interest for each given node may be orthogonally projecting onto a plane tangential to the given node. With this projection a two-dimensional conduction velocity vector is estimated for the given node in the projection of the region of interest based on analysis of the activation computed at the given node relative to other nodes residing in the region of interest. A 3D conduction velocity vector for the given node is then estimated on the original 3D surface, such as by projecting the 2D velocity vector (in the 2D projection) back to the original 3D space, e.g., corresponding to surface mesh.
In the example where surface meshes are used, the estimated local CV vector may be considered an apparent velocity because, for example, the computed speed and direction in the velocity vector for a given node are a combination of activation wave front velocities from the entire local myocardium mass within the region of such given node. The resulting conduction velocity map provides a local vector for each of the nodes that includes speed and/or propagation direction for the nodes in each region of interest over one or more time intervals. The conduction velocity map thus provides an efficient estimate of propagation pattern.
In some examples, the projection from the 2D domain to the 3D domain may be guided by aligning one or more anatomical landmarks (e.g., on the 3D cardiac surface). Additionally, depending on the relative density of the mesh structures in 2D and 3D domains, interpolation may be applied to project and fit computed 2D conduction velocity data for nodes in 2D space (e.g., a given distribution of nodes) to corresponding velocity data at a distribution of nodes across an anatomical 3D mesh (e.g., a different distribution of nodes). Such interpolation may be from a dense distribution (in 2D space) to a more sparse distribution of nodes (in 3D space) or from a sparse distribution (in 2D space) to more dense node distribution (in 3D space). This allows the calculated velocity values to be visualized on detailed anatomy geometry (e.g., 3D cardiac surface geometry), such as may be derived from an imaging modality (e.g., CT and/or MRI).
At 20, one or more pattern maps are generated, to facilitate locating specific tachycardia propagation patterns. As an example, one or more rhythm features can be derived from applying mathematical operators to the conduction velocity computed at 18. For instance, the maps may include direction maps, speed maps, dispersion maps and rotation maps which may be visualized on a display, such as shown in
As one example, dispersion features are computed by applying the divergent operator to the conduction velocity vector over one or more time intervals. A corresponding dispersion map (or maps) can be output (as an electrical signal via an output port) to a display device to render a corresponding visualization of the dispersion map(s). The resulting map can visualize divergent and collision areas, such as to help diagnose of focal tachycardias and isthmus barriers.
As another example, rotation features are computed by applying the curl operator to the conduction velocity vector. A corresponding rotation map (or maps) can be output (as an electrical signal via an output port) to a display device to render a corresponding visualization of the rotation map(s). The resulting rotation map can visualize clockwise and counter-clockwise rotations, such as to help diagnose macro-reentrant tachycardia and rotation areas.
At the end of the process, three fields are computed: the CV vector field, the divergence, and the curl fields. The aim of the proposed propagation pattern mapping tool is to compute a CV vector field and derive from it the divergence and curl fields on the 3-D mesh.
Conduction Velocity (CV) Mapping
Selection of the Region of Interest (ROI)
At 52, an ROI is selected for a given node ni. As one example of a 2D reference technique, a ROI is selected around the given node ni regarding spatial and temporal axes. First, the neighboring nodes are identified spatially using the triangulation of the mesh. For example, the first and second-ranks (spatially) neighboring nodes are selected for each respective node ni. For instance, first rank nodes are directly connected to the given node ni by one edge and second rank nodes are connected to the given node by two contiguous edges.
Next, a temporal selection is performed with respect to the given node ni, such as to ensure the continuity of the propagation in case of reentrant circuits and thereby enable estimation of a realistic CV vector field. For example, spatially selected neighboring nodes for which the activation is too delayed (e.g., delayed an amount greater than a time threshold) from the activation at n0 are discarded. For example, in contrast to using an arbitrary temporal threshold, in this example, temporal selection is done automatically based on the cycle length of the activation tcycle. The cycle length may be computed as a time between sequential activations or some other identifiable signal feature between sequential activations. As a further example, neighboring node nj, activated at time tj is selected to be part of the region of the given node n0 activated at time t0 if:
For example, the time window tw is set based on an estimated cycle length of the arrhythmia, such as estimated from a body surface ECG (e.g., a standard 12 lead ECG) or from the unipolar electrograms (e.g., on the body surface or reconstructed onto the cardiac envelope of interest). The time window may be automatically determined or it may be set in response to a user input. At the end of this spatio-temporal ROI selection, the region of interest centered on node n0 is composed of N neighboring nodes. This automated ROI selection method further enables selecting the nodes without an arbitrary threshold that could lead to artifactual results if they are badly identified while preserving the local aspect of the estimation. In some examples, if the resulting number of neighboring nodes in region of interest is less than a predetermined number, the conduction velocity vector may not be computed for the given node. The predetermined number can be fixed number or set in response to a user input, such as to ensure a sufficient number of neighboring nodes in each neighborhood of the given node. In other examples, the size of the neighborhood may be expanded.
Estimation of the Delays
At 54, delays are estimated. In some examples, the delays are more relevant than the activation time marker placement. Accordingly, the delays are included in the example estimation for non-invasive reconstruction of unipolar electrograms (EGMs). The quality of the CV vector estimation depends on the number of neighbors in the region of interest around n0, i.e., N. The adding of the activation time delays between the given node and its neighbors increases the redundancy of the temporal information and so the accuracy of the estimation.
To mitigate uncertainty of the placement of the activation time markers on unipolar EGMs, delays are estimated, for example, using Shors' method, such as disclosed in S. M. Shors, et al., “A method for determining high-resolution activation time delays in unipolar cardiac mapping,” IEEE Trans. Biomed. Eng., vol. 43, no. 12, pp. 1192-1196, 1996. By way of example, to calculate the delays, the delay is estimated by applying the Hilbert transformation to the cross-correlation between the unipolar reconstructed EGM of a neighbor node n belonging to the ROI and the given node no. The time of the maximum correlation correspond to the delay between these two nodes. The corresponding delay τ is the timing for which this transform is equal to zero. The delay may be estimated independently of determining the orthogonal projection (i.e., the delay may be computed prior to identifying the plane tangential to the ROI).
Orthogonal Projection of the Region of Interest
At 56, an orthogonal projection of the ROI is computed. For example, a non-destructive local projection of the ROI (selected at 52) is used. From the geometry of the 3D mesh obtained with high-density electroanatomic map and the local aspect of the ROI, for example, the curvature may be assumed low or negligible. Thus, to estimate the 3-D CV vector at n0, an orthogonal projection of the N nodes onto the plane Π0 tangential at n0 to its region of interest. The use of a local tangential plane is more appropriate than a global tangential plane for less smooth surface geometries, such as atria.
As one example, as the region of interest is a part of an ellipsoid centered on the given node n0, the normal vector of Π0 is collinear to one of the three semi-principal axes of the ROI ellipsoid. Hence, an efficient way to estimate the normal vector of the tangential plane is to use singular value decomposition applied on all the 3-D coordinates of the N neighboring nodes for the given node n0 (e.g., wherein N is the remaining neighboring nodes that have not been removed during 52). For example, let P0=[p0,x p0,y p0,z] be the 3×3 orthonormal transfer matrix obtained from this decomposition. Then, locally, it is assumed the curvature is the lowest spatial dispersion of the ROI. As a result, with the conventional ordering for singular value decomposition, the two first vectors of the decomposition span Π0, whereas the third vector is orthogonal to Π0 which is linked to curvature. Then, the matrix defined by [p0,x p0,y] is used to orthogonally project the ROI onto Π0. This enables the 2D spatial coordinates of the N nodes in the region of interest and the studied node n0 to be determined as {(x,y)}1 . . . N and (x0,y0), respectively.
Estimation of the 3D Local Velocity Vector
At 58, a local 3D velocity vector is estimated for the given node. For example, after determining the orthogonal projection of the ROI (at 54), a 2D local velocity vector for the ROI is estimated (e.g., on the 2D projection plane based spatial and temporal electrical information for each ROI. The spatial and temporal information can include activation times and delays determined for each of the neighboring nodes in the ROI. As an example, the activation times and delays are determined (e.g., at 54) independently prior to or otherwise from the local tangential plane for the ROI on which the 2D local velocity vector is estimated.
In some examples, the spatio-temporal information of the N nodes in the ROI is fit to a polynomial surface model t={tilde over (t)}(x, y). This model may be customized with the neighboring activation times and estimated delays, and locally assumes a quadratic activation wave front at n0, such as shown below in Eq. (2).
where (x,y)={tilde over (t)}(x,y)−{tilde over (t)}(x0,y0).
In some examples, if N is greater than 3 (N>3), the model parameters (e.g., a, b, c, e, and f) are estimated by minimizing the least square errors between 1) the recorded activation times of the N nodes in the region of interest {t}1 . . . N and the model estimated {{tilde over (t)}(x,y)}1 . . . N; and 2) the estimated delays {}1 . . . N and the model estimated {
(x,y)}1 . . . N, such as shown in the example of Eq. (3):
The local 2D velocity vector w0 for the ROI then may be calculated by applying the 2D spatial gradient ∇ on {tilde over (t)}, estimated at n0, such as demonstrated in Eq. (4):
Finally, the 3-D velocity vector v0 is directly obtained from the projection of w0 back on to the original space with [p0,x p0,y] (5).
ν0=[p0,xp0,y]·w0 (5)
At 60, a determination is made whether there are additional nodes across the surface for which the conduction velocity is to be determined. In some examples, the set of nodes may include the entire surface or, in other examples, it may be for nodes distributed across the surface corresponding to one or more regions (e.g., an atrial region, both atria, ventricle or ventricles etc.). If there are more nodes, the method proceeds to 62 in which the next node is identified as the given node n0 for which the conduction velocity is to be computed, and the method repeats. This process (e.g., 52-62) is repeated for each node of interest across the 3-D mesh. At the end of this process, in response to determining that there are no additional nodes, the 3-D CV vector field is defined with a velocity, v, at each node and the value is stored in memory and used to generate a CV map at 64 for each time index in the selected one or time intervals. For example, the each node can include a velocity vector for each time index, including both a direction of propagation value and apparent local speed value associated with each respective node. At 64, the corresponding conduction velocity map may thus be generated and presented on a display based on the stored CV data, such as disclosed herein.
Rhythm Feature Mapping
At 102, an orthogonal projection of the ROI and the CV vector (from
At 104, the mathematical divergent (div) of the unitary CV vector is computed. The div of the unitary CV vector (U0) may be computed based on the CV vectors determined for each of nodes in its local 2D projection over one or more time intervals. At 106, the mathematical curl (rot) of the unitary CV vector U0 is computed. The curl of the unitary CV vector U0 may be computed based on the CV vectors determined for each of nodes in its local 2D projection over one or more time intervals.
By way of example let U0={u}1 . . . N be the unitary CV vector field respectively associated with {v}1 . . . N, the CV vectors of the N nodes belonging to the region of interest of a given node n0. It is assumed that the local curvature of the region of interest is low. It further may be assumed that the N vectors of U0 are coplanar within the local plane Π0. In these conditions, estimating the local curl and divergence in the 3-D Euclidian basis amounts to estimating them in the orthonormal basis, such as may be defined by the triplet P0=[p0,x p0,y P0,z] as demonstrated in Eqs. (6) and (7) below with W0={w}1 . . . N, the unitary CV vectors U0 expressed in the basis defined by P0, as follows:
div(U0)=div(W0) (6)
rot(U0)=P0·rot(W0) (7)
In some examples, assuming the 3D mesh is spherical, the curl operator is modified to get the vector (bi
At 108, local operators for the div and curl expressions are estimated. To estimate locally the divergent div(W0) and the curl rot(W0), a linear model linking each direction component of the local direction of the ROI is parameterized. For example, the linear model can be used to link the spatial coordinates of the N nodes in P0, {(x, y, z)}1 . . . N, to W0. As these unitary vectors are assumed to be coplanar and belong to the plane Π0 defined by the normal vector p0,z, the variations along this normal axis are equal to zero. Thus, to estimate div(W0) and rot(W0), two 2D linear propagation models can be used, namely {tilde over (w)}x and {tilde over (w)}y, linking respectively the component of W0 along p0,x (8) and p0,y (9) with the 2D spatial coordinates x and y such as follows:
{tilde over (w)}x(x,y)=αx(x−x0)+βx(y−y0)+γx (8)
{tilde over (w)}y(x,y)=αy(x−x0)+βy(y−y0)+γy (8)
where αx, αy, βx, βy, γx and γy are the model coefficients.
Then, the model components are calculated by minimizing the least square error between respectively the model estimated unitary vector components {tilde over (w)}x and {tilde over (w)}x as in Eq. (10) for the first model, and between {tilde over (w)}y and wy as in Eq. (11) for the second model.
From this model, each of the partial derivatives in (6) and (7) are estimated at n0 to get respectively its local divergence (12) and curl (13) in the original space at the given node n0, such as may be expressed as:
div(U0)=αx+βy (12)
rot(U0)=(αy−βx)p0,z (13)
The systems and methods disclosed herein facilitate non-invasively identification of atrial tachycardia mechanisms. This can be performed as a bedside application to establish the strategy before the procedure. This is accomplished by one or more of the following:
As an example, a catheter, such as a pacing or an ablation catheter, having one or more therapy delivery devices 156 affixed thereto can be inserted into a patient's body 154 (e.g., into or onto) the patient's heart 152, endocardially or epicardially. Various types and configurations of therapy delivery devices 156 that can be utilized, which can vary depending on the type of treatment and the procedure. For instance, the therapy device 156 can be configured to deliver electrical therapy, chemical therapy, sound wave therapy, thermal therapy or any combination thereof. The placement of the device 156 can be guided based on mapping conduction velocity and/or rhythm patterns across the cardiac surface as disclosed herein. The guidance can be automated, semi-automated or be manually implemented based on information provided.
As a further example, the delivery device 156 can include one or more electrodes located at a tip of an ablation catheter configured to generate heat for ablating tissue in response to electrical signals (e.g., radiofrequency energy) supplied by a therapy system 158. In other examples, the therapy delivery device 156 can be configured to deliver cooling to perform ablation (e.g., cryogenic ablation), to deliver chemicals (e.g., drugs), ultrasound ablation, high-frequency ablation, or a combination of these or other therapy mechanisms. In still other examples, the therapy delivery device 156 can include one or more electrodes located at a tip of a pacing catheter to deliver electrical stimulation, such as for pacing the heart, in response to electrical signals (e.g., pacing pulses) supplied by a therapy system 158. Other types of therapy can also be delivered via the therapy system 158 and the invasive therapy delivery device 156 that is positioned within the body.
In some examples, the therapy system 158 can be located external to the patient's body 154 and be configured to control therapy that is being delivered by the device 156. For instance, the therapy system 158 includes controls (e.g., hardware and/or software) 160 that can communicate (e.g., supply) electrical signals via a conductive link electrically connected between the delivery device (e.g., one or more electrodes) 156 and the therapy system 158. The control system 160 can control parameters of the signals supplied to the device 156 (e.g., current, voltage, repetition rate, trigger delay, sensing trigger amplitude) for delivering therapy (e.g., ablation or stimulation) via the electrode(s) 154 to one or more location of the heart 152. The controls 160 can set the therapy parameters and apply stimulation based on automatic, manual (e.g., user input) or a combination of automatic and manual (e.g., semiautomatic) controls. One or more sensors (not shown) can also communicate sensor information (e.g., electrical and/or thermal measurements) back to the therapy system 158. The position of the device 156 relative to the heart 152 can be determined and tracked intraoperatively via an imaging modality (e.g., fluoroscopy, x-ray), a mapping system 162, direct vision or the like. The location of the device 156 and other therapy parameters thus can be combined and used by the controls 160 to determine corresponding therapy parameter data.
Before, during and/or after delivering a therapy via the therapy system 158, another system or subsystem can be utilized to acquire electrophysiology information for the patient. In the example of
One or more sensors may also be located on the device 156 that is inserted into the patient's body. Such sensors can be utilized separately or in conjunction with the non-invasive sensors 164 for mapping electrical activity for an endocardial surface, such as the wall of a heart chamber, as well as for an epicardial surface. Additionally, such electrode can also be utilized to help localize the device 156 within the heart 152, which can be registered into an image or map that is generated by the system 150. Alternatively, such localization can be implemented in the absence of emitting a signal from an electrode within or on the heart 152.
In each of such example approaches for acquiring patient electrical information, including invasively, non-invasively, or a combination of invasive and non-invasive sensing, the sensor array(s) 164 provide the sensed electrical information (e.g., time-varying, unipolar electrical potential measurements distributed across the thorax) to a corresponding measurement system 166. The measurement system 166 can include appropriate controls and signal processing circuitry 168 for providing corresponding measurement data 170 that describes electrical activity detected by the sensors in the sensor array 164. The measurement data 170 can include analog and/or digital information.
The control 168 can also be configured to control the data acquisition process (e.g., sample rate, line filtering) for measuring electrical activity and providing the measurement data 170. In some examples, the control 168 can control acquisition of measurement data 170 separately from the therapy system 158 operation, such as in response to a user input. In other examples, the measurement data 170 can be acquired concurrently with and in synchronization with delivering therapy by the therapy system, such as to detect electrical activity of the heart 152 that occurs in response to applying a given therapy (e.g., according to therapy parameters). For instance, appropriate time stamps can be utilized for indexing the temporal relationship between the respective measurement data 170 and therapy parameters use to deliver therapy as to facilitate the evaluation and analysis thereof.
Since the measurement system 166 can measure electrical activity of a predetermined region or the entire heart concurrently (e.g., where the sensor array 164 covers the entire thorax of the patient's body 154), the resulting output data (e.g., visualizing attributes of identified conduction velocity and associated rhythm patterns, such as propagation, rotational activity and/or other characteristsics) thus can also represent concurrent data for the predetermined region or the entire heart in a temporally and spatially consistent manner. The one or more time interval for which the output data/maps are computed can be selected based on user input (e.g., selecting a timer interval from one or more waveforms). Additionally or alternatively, the selected intervals can be synchronized with the application of therapy by the therapy system 158.
The mapping system 162 includes signal processing methods 178 to preprocess the electrical measurement data 170. For example, such preprocessing implement any of the types of preprocessing disclosed at 12 in
For the example where the electrical measurement data is obtained non-invasively (e.g., via body surface sensor array 164), electrogram reconstruction 180 can be programmed to compute an inverse solution and provide corresponding reconstructed electrograms (e.g., unipolar electrograms) on a cardiac envelope based on the preprocessed electrical signals (from signal preprocessing 178) and the geometry data 172. For example, the reconstructed electrograms thus can be estimated for each of the nodes on across a 3D mesh surface that defines a cardiac envelope over one or more time intervals.
Examples of inverse algorithms that can be implemented by the electrogram reconstruction 180 in the system 150 include those disclosed in U.S. Pat. Nos. 7,983,743 and 6,772,004, which are incorporated herein by reference. The EGM reconstruction 180 thus can reconstruct the body surface electrical activity measured via the sensor array 164 onto a multitude of locations (nodes) on the cardiac envelope (e.g., greater than 1000 locations, such as about 2000 locations or more). In some examples, the EGM reconstruction 180 of the mapping system 162 can compute electrical activity over a sub-region of the heart based on electrical activity measured invasively, such as via a basket catheter or other form of measurement probe.
As disclosed herein, the cardiac envelope can correspond to a three dimensional surface geometry corresponding to a patient's heart, which surface can be epicardial or endocardial. Alternatively or additionally, the cardiac envelope can correspond to a three-dimensional geometric surface that resides between the epicardial surface of a patient's heart and the surface of the patient's body where the sensor array 164 has been positioned. Additionally, the geometry data 172 that is utilized by the electrogram reconstruction 180 can correspond to actual patient anatomical geometry, a preprogrammed generic model or a combination thereof (e.g., a model that is modified based on patient anatomy).
As an example, the geometry data 172 may includes data representing the cardiac envelope as a three dimensional polygonal mesh, including nodes connected by edges, as well as coordinates for the locations on the patient's torso where each of the sensing electrodes in the array 164 reside, such as derived from image data acquired for the patient. For instance, image processing can include extraction and segmentation of anatomical features, including one or more organs and other structures, from a digital image set. Additionally, a location for each of the electrodes in the sensor array 164 can be included in the patient geometry data 172, such as by acquiring the image while the electrodes are disposed on the patient and identifying the electrode locations in a coordinate system through appropriate extraction and segmentation. Other non-imaging based techniques can also be utilized to obtain the position of the electrodes in the sensor array, such as a digitizer or manual measurements.
As mentioned above, the geometry data 172 can correspond to a mathematical model, such as can be a generic model or a model that has been constructed based on image data for the patient. Appropriate anatomical or other landmarks, including locations for the electrodes in the sensor array 164 can be identified in the geometry data 172 to facilitate registration of the electrical measurement data 170 and performing the electrogram reconstruction 180 thereon. The identification of such landmarks can be done manually (e.g., by a person via image editing software) or automatically (e.g., via image processing techniques).
By way of further example, the geometry data 172 can be acquired using nearly any imaging modality based on which a corresponding representation of the 3D geometrical surface can be constructed (e.g., a generic or patient-specific mesh structure), such as a triangular mesh of mesh of other size polygons as described herein. Such imaging may be performed concurrently with recording the electrical activity that is utilized to generate the patient measurement data 170 or the imaging can be performed separately (e.g., before or after the measurement data has been acquired).
The mapping system 162 includes a CV calculator 182 programmed (e.g., executable software instructions corresponding to method 50) to compute a three-dimensional conduction velocity map for nodes across the cardiac envelope. By way of brief example, the CV calculator includes an ROI calculator to determine a region of interest around a given node of the plurality of nodes on the surface of interest. The CV calculator further includes a 3D-to-2D projection component programmed to project the region of interest around the given node from the surface of interest onto 2D surface (e.g., a plane) tangential to the given node. The CV calculator further includes a 2D CV calculator to estimate a two-dimensional CV vector for the given node based on analysis of the electrophysiological signal for the given node relative to other nodes in the ROI. The CV calculator 182 further is programmed to estimate the 3D CV vector for the given node on the 3D surface of interest (e.g., of the heart surface) based on the 2D CV vector.
As a further example, the CV calculator can compute the CV at each of nodes distributed across the mesh surface that represents the cardiac envelope, as provided by the geometry data 172. The CV calculator 182 may operate to compute the CV map data in real-time (e.g., concurrently) with acquisition of body surface electrical potential measurements (e.g., by sensor array 164). Alternatively or additionally, the CV calculator 182 may operate off-line to compute the CV map data following the acquisition of body surface electrical potential measurements and reconstruction of the electrograms by EGM reconstruction 180.
The mapping system 162 also may include a pattern calculator 184 programmed (e.g., executable software instructions corresponding to method 100) to compute one or more pattern feature maps based on the conduction velocity (e.g., computed by CV calculator 182). For example, the pattern calculator 184 may compute one or more pattern feature maps, such as including propagation direction, speed, rotation and/or dispersion, such as disclosed herein.
The mapping system 162 includes an output generator 186 that is programmed to generate output data 174 corresponding to graphical maps (e.g., computed electrophysiological information spatially represented on a cardiac envelope) based on computed output data, such as disclosed herein above. For example, the output generator 186 can convert the computed map data (e.g., CV and/or pattern map data) into output data (e.g., primitives) to drive a graphics pipeline of a computing system on which the mapping system 162 is implemented. The output data 174 can include data provided to the display 192 for visualization of one or more graphical maps 194, such as disclosed herein.
In some examples, the output generator 186 includes a graphical user interface (GUI) 190 to control the output generator and to set one or more parameters associated with the displayed graphical representation, corresponding to an output visualization of the computed map. For example, a user can configure and select various parameters in response to a user input via the GUI 190, such as including selecting a time interval, selecting a type of map (or maps) that is to be presented in the visualization and the like. Additionally, a user can employ the GUI 190 to selectively program one or more parameters (e.g., temporal and spatial thresholds, filter parameters and the like) utilized by the CV calculator 182 and/or pattern calculator 184 and/or to select one or more sample time intervals to set a time duration for the electrical data 170 that is utilized by the mapping system 162.
The output generator 186 thus provides output data 174, which can be used to visualize a corresponding graphical output in a display 192, such as including an electrocardiographic map 194. For example, the system 162 can generate maps and other output visualizations, such as including but not limited to the maps and other output visualizations disclosed herein. The mapping system 162 and its respective functions 178, 180, 182, 184, 186 and 190 can be implemented as machine-readable software code instructions stored in memory and executable by one or more processors of a computer system. The geometry data 172 and measurement data 170 also may be stored in the memory and be accessed by the processor(s) for performing such functions.
Additionally, in some examples, the mapping system 162 can provide the output data 174 to be utilized by the therapy system 158. The control 160 that is implemented can be fully automated control, semi-automated control (partially automated and responsive to a user input) or manual control based on the output data 174. In some examples, the control 160 of the therapy system can utilize the output data 174 to control one or more therapy parameters. As an example, the control 160 can control delivery of ablation therapy to a site of the heart (e.g., epicardial or endocardial wall) based on activation timing or activation wave fronts determined by one or both of the calculators 182 and 184. In other examples, an individual can view the map 194 generated in the display to manually control the therapy system. Other types of therapy and devices can also be controlled based on the output data 174 and corresponding graphical map 194.
The GUI 200 also includes buttons for selectively displaying rhythm pattern features, including rotation patterns and/or dispersion patterns, across the cardiac surface. For example, a user can selectively activate/deactivate a rotation map (via button 220), as well as counter-clockwise flutter in the rotation map (via button 222) and clockwise flutter in the rotation map (via button 224). The rotation map can be computed, for example, from applying a curl operator to the CV map as disclosed herein. In further examples a user can selectively activate/deactivate a divergent map (via button 226), counter-clockwise flutter in the rotation map (via button 228) and clockwise flutter in the rotation map (via button 230). The divergent map can be computed, for example, to determine the dispersion features by applying divergent operators to the CV map as disclosed herein. One or more intervals of reconstructed electrical data can also be selected by a user in response to a user input activating a select interval button 234 to implement functions for selecting the interval(s) of interest.
In view of the foregoing structural and functional description, those skilled in the art will appreciate that portions of the invention may be embodied as a method, data processing system, or computer program product. Accordingly, these portions of the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware.
Furthermore, portions of the invention may be a computer program product on a computer-usable storage medium having computer readable program code on the medium. Any suitable computer-readable medium may be utilized including, but not limited to, static and dynamic storage devices, hard disks, optical storage devices, and magnetic storage devices.
Certain embodiments of the invention have also been described herein with reference to block illustrations of methods, systems, and computer program products. It will be understood that blocks of the illustrations, and combinations of blocks in the illustrations, can be implemented by computer-executable instructions. These computer-executable instructions may be provided to one or more processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus (or a combination of devices and circuits) to produce a machine, such that the instructions, which execute via the processor, implement the functions specified in the block or blocks.
These computer-executable instructions may also be stored in computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory result in an article of manufacture including instructions which implement the function specified in the flowchart block or blocks. The computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks.
What have been described above are examples. It is, of course, not possible to describe every conceivable combination of components or methods, but one of ordinary skill in the art will recognize that many further combinations and permutations are possible. Accordingly, the invention is intended to embrace all such alterations, modifications, and variations that fall within the scope of this application, including the appended claims. Where the disclosure or claims recite “a,” “an,” “a first,” or “another” element, or the equivalent thereof, it should be interpreted to include one or more than one such element, neither requiring nor excluding two or more such elements. As used herein, the term “includes” means includes but not limited to, the term “including” means including but not limited to. The term “based on” means based at least in part on.
Dubois, Rémi, Dallet, Corentin
Patent | Priority | Assignee | Title |
Patent | Priority | Assignee | Title |
10285611, | Aug 07 2015 | Boston Scientific Scimed, Inc | Cardiac mapping using a 3D grid |
20030158483, | |||
20110137171, | |||
20150216438, | |||
20160045133, | |||
20160100770, | |||
20170014086, | |||
WO2015134237, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
May 09 2018 | CARDIOINSIGHT TECHNOLOGIES, INC. | (assignment on the face of the patent) | / | |||
May 09 2018 | DUBOIS, RÉMI | CARDIOINSIGHT TECHNOLOGIES, INC | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 045777 | /0792 | |
May 09 2018 | DALLET, CORENTIN | CARDIOINSIGHT TECHNOLOGIES, INC | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 045777 | /0792 |
Date | Maintenance Fee Events |
May 09 2018 | BIG: Entity status set to Undiscounted (note the period is included in the code). |
Aug 20 2024 | M1551: Payment of Maintenance Fee, 4th Year, Large Entity. |
Date | Maintenance Schedule |
Mar 30 2024 | 4 years fee payment window open |
Sep 30 2024 | 6 months grace period start (w surcharge) |
Mar 30 2025 | patent expiry (for year 4) |
Mar 30 2027 | 2 years to revive unintentionally abandoned end. (for year 4) |
Mar 30 2028 | 8 years fee payment window open |
Sep 30 2028 | 6 months grace period start (w surcharge) |
Mar 30 2029 | patent expiry (for year 8) |
Mar 30 2031 | 2 years to revive unintentionally abandoned end. (for year 8) |
Mar 30 2032 | 12 years fee payment window open |
Sep 30 2032 | 6 months grace period start (w surcharge) |
Mar 30 2033 | patent expiry (for year 12) |
Mar 30 2035 | 2 years to revive unintentionally abandoned end. (for year 12) |