Source: http://patents.com/us-9589379.html
Timestamp: 2017-03-24 22:12:50
Document Index: 379779891

Matched Legal Cases: ['Application No. 2010', 'Application No. 2012', 'Application No. 2013', 'Application No. 2015', 'Application No. 2013', 'Application No. 2015', 'Application No. 2013']

US Patent # 9,589,379. System and method for visualization of cardiac changes under various
pacing conditions - Patents.com
United States Patent 9,589,379
System and method for visualization of cardiac changes under various
pacing conditions
A system and method for visualization of cardiac changes under various
pacing conditions for intervention planning and guidance is disclosed. A
patient-specific anatomical heart model is generated based on medical
image data of a patient. A patient-specific computational model of heart
function is generated based on patient-specific anatomical heart model. A
virtual intervention is performed at each of a plurality of positions on
the patient-specific anatomical heart model using the patient-specific
computational model of heart function to calculate one or more cardiac
parameters resulting from the virtual intervention performed at each of
the plurality of positions. One or more outcome maps are generated
visualizing, at each of the plurality of positions on the
patient-specific anatomical heart model, optimal values for the one or
more cardiac parameters resulting from the virtual intervention performed
at the that position on the patient-specific anatomical heart model.
Inventors: Mansi; Tommaso (Plainsboro, NJ), Passerini; Tiziano (Plainsboro, NJ), Kamen; Ali (Skillman, NJ), Georgescu; Bogdan (Plainsboro, NJ), Comaniciu; Dorin (Princeton Junction, NJ) Applicant: Name City State Country Type SIEMENS HEALTHCARE GMBH Erlangen N/A DE Assignee:
1000002444251
14/741,540
Prior Publication Data Document IdentifierPublication Date US 20150371437 A1Dec 24, 2015 Related U.S. Patent Documents Application NumberFiling DatePatent NumberIssue Date 62016220Jun 24, 2014 Current U.S. Class: 1/1 Current CPC Class: G06T 11/60 (20130101); A61N 1/3627 (20130101); G06T 2207/10081 (20130101); G06T 2207/10088 (20130101); G06T 2207/10132 (20130101); G06T 2207/30048 (20130101)
Current International Class: A61N 1/362 (20060101); G06T 11/60 (20060101)
6236878
7668354
7916919
8295913
2008/0300588
2009/0124896
2009/0142740
2009/0221999
2010/0040272
2011/0015628
2011/0251607
Kruecker et al.
2012/0022843
Ionasec et al.
2012/0041301
2012/0072190
2012/0121151
2012/0203530
2013/0158577
2013/0197881
Mansi et al.
2013/0197884
2013/0216110
2013/0226542
Rapaka et al.
2014/0022250
2014/0088943
Trayanova et al.
2014/0122048
Vadakkumpadan et al.
2014/0136174
Audigier et al.
2015/0042646
2014113746
2014133924
Other References Dikici et al., "Quantification of Delayed Enhancement MR Images", In Proc. MICCAI 2004, LNCS 3216, pp. 250-257, 2004. cited by applicant
.Peyrat et al., "A Computational Framework for the Statistical Analysis of Cardiac Diffusion Tensors: Application to a Small Database of Canine Hearts", IEEE TMI, 26(11): 1500-1514, 2007. cited by applicant
.Mitchell et al., "A Two-Current Model for the Dynamics of Cardiac Membrane", Bulletin of Mathematical Biology, 65(5):767-793, 2003. cited by applicant
.Powell, "Developments of NEWUOA for Minimization Without Derivatives", J. Num. Analysis, 2008. cited by applicant
.Rapaka, et al., "LBM-EP: Lattice-boltzmann method for fast cardiac electrophysiology simulation from 3d images", In: Medical Image Computing and Computer-Assisted Intervention MICCAI, 2012, pp. 33-40, vol. 7511, Springer Berlin Heidelberg. cited by
applicant. Primary Examiner: Johns; Andrew W
62/016,220, filed Jun. 24, 2014, the disclosure of which is herein
Claims The invention claimed is: 1. A method for guiding a cardiac intervention, comprising: generating a patient-specific anatomical heart model based on medical image data of a patient; generating
a patient-specific computational model of heart function based on the patient-specific anatomical heart model; performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the
patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions; and generating one or more outcome maps visualizing, at each of
the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model.
2. The method of claim 1, further comprising: displaying the one or more outcome maps on a display device.
3. The method of claim 1, further comprising: overlaying at least one of the one or more outcome maps on an interventional image acquired during an interventional procedure.
4. The method of claim 1, further comprising: fusing at least one of the one or more outcome maps with a patient-specific model of coronary arteries.
5. The method of claim 1, wherein performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart function to calculate one or more
cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions comprises: performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more
electrical pacing leads using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions.
6. The method of claim 5, wherein performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more electrical pacing leads using the patient-specific computational model of heart
function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions comprises: for each of the plurality of lead positions, performing the virtual pacing using a plurality of pacing
7. The method of claim 6, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the
virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating a respective outcome map visualizing an optimal value for the cardiac parameter at each of the
plurality of lead positions resulting from virtual pacing at each of the plurality of lead positions using the plurality of pacing protocols.
8. The method of claim 7, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the
virtual intervention performed at the that position on the patient-specific anatomical heart model further comprises: for each of the cardiac parameters, generating a second respective outcome map visualizing, for each of the plurality of lead positions,
a pacing protocol resulting in the optimal value for the cardiac parameter at that lead position.
9. The method of claim 5, wherein the virtual pacing is virtual cardiac resynchronization therapy and performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more electrical
pacing leads using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions comprises: performing virtual pacing for a
plurality of left ventricle lead positions with a fixed right ventricle lead position using one or more pacing protocols at each of the plurality of left ventricle positons.
10. The method of claim 9, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the
virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating a respective left ventricle-only outcome map visualizing an optimal value for the cardiac
parameter at each of the plurality of left ventricle lead positions resulting from virtual pacing at each of the plurality of left ventricle lead positions using the one or more pacing protocols.
11. The method of claim 5, wherein the virtual pacing is virtual cardiac resynchronization therapy and performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more electrical
plurality of right ventricle lead positions with a fixed left ventricle lead position using one or more pacing protocols at each of the plurality of right ventricle positons.
12. The method of claim 11, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the
virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating a respective right ventricle-only outcome map showing an optimal value for the cardiac parameter
at each of the plurality of right ventricle lead positions resulting from virtual pacing at each of the plurality of right ventricle lead positions using the one or more pacing protocols.
13. The method of claim 5, wherein the virtual pacing is virtual cardiac resynchronization therapy and performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more electrical
pacing leads using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions comprises: for each of a plurality of left
ventricle lead positions, performing virtual pacing using one or more pacing protocols at each of a plurality of right ventricle lead positons using the patient-specific computational model of heart function to calculate the one or more cardiac
parameters resulting from the virtual pacing performed using the one or more pacing protocols at each of the plurality of right ventricle lead positions.
14. The method of claim 13, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the
virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating bi-lead outcome maps visualizing, for each of the plurality of left ventricle lead positions, an
optimal value for the cardiac parameter at that left ventricle lead position, a corresponding one of the plurality of right ventricle lead positions resulting in the optimal value for the cardiac parameter at that left ventricle lead position, and a
pacing protocol resulting in the optimal value for the cardiac parameter at that left ventricle lead position.
15. The method of claim 1, wherein the one or more cardiac parameters comprise one or more of QRS duration, electrical axis, QT duration, ejection fraction, stroke volume, or myocardium strain.
16. The method of claim 1, wherein performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart function to calculate one or more
cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions comprises: performing a virtual ablation at each of a plurality of ablation positions on the patient-specific anatomical heart model using the
patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual ablation performed at each of the plurality of positions.
17. The method of claim 1, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the
virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: generating, for each of the one or more cardiac parameters, a respective 3D color map representing the optimal values at the plurality of
positions on the patient-specific anatomical heart model for that cardiac parameter.
18. The method of claim 1, wherein the patient-specific computational model of heart function comprises a cardiac electrophysiology model, a cardiac hemodynamics model, and a cardiac biomechanics model.
19. A method for guiding an intervention procedure in a target organ of a patient, comprising: generating a patient-specific anatomical model of the target organ based on medical image data of the patient; generating a patient-specific
computational model of organ function based on the patient-specific anatomical model of the target organ; performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical model of the target organ using the
patient-specific computational model of organ function to calculate one or more outcome parameters resulting from the virtual intervention performed at each of the plurality of positions; and generating one or more outcome maps visualizing, at each of
the plurality of positions on the patient-specific anatomical model of the target organ, optimal values for the one or more outcome parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical model
of the target organ.
20. An apparatus for guiding a cardiac intervention, comprising: means for generating a patient-specific anatomical heart model based on medical image data of a patient; means for generating a patient-specific computational model of heart
function based on the patient-specific anatomical heart model; means for performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart
function to calculate one or more cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions; and means for generating one or more outcome maps visualizing, at each of the plurality of positions on the
patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model.
21. The apparatus of claim 20, further comprising: means for overlaying at least one of the one or more outcome maps on an interventional image acquired during an interventional procedure.
22. The apparatus of claim 20, wherein the means for performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart function to
calculate one or more cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions comprises: means for performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical
heart model for one or more electrical pacing leads using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions.
23. The apparatus of claim 22, wherein the means for performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more electrical pacing leads using the patient-specific
computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions comprises: means for performing the virtual pacing using a plurality of pacing
protocols for each of the plurality of lead positions.
24. The apparatus of claim 23, wherein the means for generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters
resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: means for generating, for each of the cardiac parameters, a respective outcome map visualizing an optimal value for the
cardiac parameter at each of the plurality of lead positions resulting from virtual pacing at each of the plurality of lead positions using the plurality of pacing protocols.
25. The apparatus of claim 24, wherein the means for generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters
resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model further comprises: means for generating, for each of the cardiac parameters, a second respective outcome map visualizing, for each of
the plurality of lead positions, a pacing protocol resulting in the optimal value for the cardiac parameter at that lead position.
26. The apparatus of claim 20, wherein the one or more cardiac parameters comprise one or more of QRS duration, electrical axis, QT duration, ejection fraction, stroke volume, or myocardium strain.
27. The apparatus of claim 20, wherein the means for performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart function to
calculate one or more cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions comprises: means for performing a virtual ablation at each of a plurality of ablation positions on the patient-specific
anatomical heart model using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual ablation performed at each of the plurality of positions.
28. The apparatus of claim 20, wherein the means for generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters
resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: means for generating, for each of the one or more cardiac parameters, a respective 3D color map representing the optimal
values at the plurality of positions on the patient-specific anatomical heart model for that cardiac parameter.
29. A non-transitory computer readable medium storing computer program instructions for guiding a cardiac intervention, the computer program instructions when executed by a processor cause the processor to perform operations comprising:
generating a patient-specific anatomical heart model based on medical image data of a patient; generating a patient-specific computational model of heart function based on the patient-specific anatomical heart model; performing a virtual intervention
at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual intervention performed at each of
the plurality of positions; and generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more cardiac parameters resulting from the virtual
intervention performed at the that position on the patient-specific anatomical heart model.
30. The non-transitory computer readable medium of claim 29, wherein the operations further comprise: overlaying at least one of the one or more outcome maps on an interventional image acquired during an interventional procedure.
31. The non-transitory computer readable medium of claim 29, wherein performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart
function to calculate one or more cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions comprises: performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical
32. The non-transitory computer readable medium of claim 31, wherein performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical heart model for one or more electrical pacing leads using the
patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions comprises: for each of the plurality of lead positions, performing
the virtual pacing using a plurality of pacing protocols.
33. The non-transitory computer readable medium of claim 32, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more
cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating a respective outcome map visualizing an optimal value for
the cardiac parameter at each of the plurality of lead positions resulting from virtual pacing at each of the plurality of lead positions using the plurality of pacing protocols.
34. The non-transitory computer readable medium of claim 33, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more
cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model further comprises: for each of the cardiac parameters, generating a second respective outcome map visualizing, for
each of the plurality of lead positions, a pacing protocol resulting in the optimal value for the cardiac parameter at that lead position.
35. The non-transitory computer readable medium of claim 31, wherein the virtual pacing is virtual cardiac resynchronization therapy and performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical
heart model for one or more electrical pacing leads using the patient-specific computational model of heart function to calculate one or more cardiac parameters resulting from the virtual pacing performed at each of the plurality of lead positions
comprises: performing virtual pacing for a plurality of left ventricle lead positions with a fixed right ventricle lead position using one or more pacing protocols at each of the plurality of left ventricle positons.
36. The non-transitory computer readable medium of claim 35, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more
cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating a respective left ventricle-only outcome map visualizing an
optimal value for the cardiac parameter at each of the plurality of left ventricle lead positions resulting from virtual pacing at each of the plurality of left ventricle lead positions using the one or more pacing protocols.
37. The non-transitory computer readable medium of claim 31, wherein the virtual pacing is virtual cardiac resynchronization therapy and performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical
comprises: performing virtual pacing for a plurality of right ventricle lead positions with a fixed left ventricle lead position using one or more pacing protocols at each of the plurality of right ventricle positons.
38. The non-transitory computer readable medium of claim 37, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more
cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating a respective right ventricle-only outcome map showing an
optimal value for the cardiac parameter at each of the plurality of right ventricle lead positions resulting from virtual pacing at each of the plurality of right ventricle lead positions using the one or more pacing protocols.
39. The non-transitory computer readable medium of claim 31, wherein the virtual pacing is virtual cardiac resynchronization therapy and performing virtual pacing for each of a plurality of lead positions on the patient-specific anatomical
comprises: for each of a plurality of left ventricle lead positions, performing virtual pacing using one or more pacing protocols at each of a plurality of right ventricle lead positons using the patient-specific computational model of heart function to
calculate the one or more cardiac parameters resulting from the virtual pacing performed using the one or more pacing protocols at each of the plurality of right ventricle lead positions.
40. The non-transitory computer readable medium of claim 39, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more
cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: for each of the cardiac parameters, generating bi-lead outcome maps visualizing, for each of the plurality
of left ventricle lead positions, an optimal value for the cardiac parameter at that left ventricle lead position, a corresponding one of the plurality of right ventricle lead positions resulting in the optimal value for the cardiac parameter at that
left ventricle lead position, and a pacing protocol resulting in the optimal value for the cardiac parameter at that left ventricle lead position.
41. The non-transitory computer readable medium of claim 29, wherein the one or more cardiac parameters comprise one or more of QRS duration, electrical axis, QT duration, ejection fraction, stroke volume, or myocardium strain.
42. The non-transitory computer readable medium of claim 29, wherein performing a virtual intervention at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart
function to calculate one or more cardiac parameters resulting from the virtual intervention performed at each of the plurality of positions comprises: performing a virtual ablation at each of a plurality of ablation positions on the patient-specific
43. The non-transitory computer readable medium of claim 29, wherein generating one or more outcome maps visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one or more
cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model comprises: generating, for each of the one or more cardiac parameters, a respective 3D color map representing the
optimal values at the plurality of positions on the patient-specific anatomical heart model for that cardiac parameter. Description BACKGROUND OF THE INVENTION
The present invention relates to visualizing the impact of electrical pacing on cardiac function of a patient, and more particularly to automatically calculating heart parameters for various pacing locations and protocols and visualizing the
cardiac changes resulting from the various pacing locations and protocols to guide lead placement for interventional therapies.
Heart failure is a major cause of death in the western world. Due to insufficient heart function, heart failure causes dyspnea and fatigue, and can also lead to cardiac arrest. For patients with a prolonged QRS-complex (e.g., QRS.gtoreq.120
ms) and low left-ventricular ejection fraction, cardiac resynchronization therapy (CRT) is a well-established treatment. CRT consists of implanting electrodes in the heart to pace the muscle artificially and "resynchronize" cardiac contraction. The
idea is to pace the myocardium at the right ventricle (RV) endocardium and the left ventricle (LV) epicardium, with given pacing intervals, to resynchronize the cardiac motion and hence re-establish the pumping efficiency of the heart. However, despite
strict guidelines, 30% of patients do not respond to CRT. Hence, better patient selection for CRT and improved guidance towards the optimal location of the pacing leads are desirable.
The present invention provides a method and system for visualization of cardiac changes under various pacing conditions for guiding lead placement in interventional therapies. Embodiments of the present invention utilize a patient-specific
computational model of cardiac electro-mechanics to calculate heart parameters for various pacing locations and various pacing protocols, and visualize the results on a patient-specific heart mesh that can be overlaid on interventional images for lead
placement guidance. Embodiments of the present invention generate outcome maps which visualize cardiac parameters for various pacing locations and protocols as colored meshes.
In one embodiment of the present invention, a patient-specific anatomical heart model is generated based on medical image data of a patient. A patient-specific computational model of heart function is generated based on patient-specific
anatomical heart model. A virtual intervention is performed at each of a plurality of positions on the patient-specific anatomical heart model using the patient-specific computational model of heart function to calculate one or more cardiac parameters
resulting from the virtual intervention performed at each of the plurality of positions. One or more outcome maps are generated visualizing, at each of the plurality of positions on the patient-specific anatomical heart model, optimal values for the one
or more cardiac parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical heart model.
In another embodiment of the present invention, a patient-specific anatomical model of a target organ is generated based on medical image data of a patient. A patient-specific computational model of organ function is generated based on
patient-specific anatomical model of the target organ. A virtual intervention is performed at each of a plurality of positions on the patient-specific anatomical model of the target organ using the patient-specific computational model of organ function
to calculate one or more outcome parameters resulting from the virtual intervention performed at each of the plurality of positions. One or more outcome maps are generated visualizing, at each of the plurality of positions on the patient-specific
anatomical model of the target organ, optimal values for the one or more outcome parameters resulting from the virtual intervention performed at the that position on the patient-specific anatomical model of the target organ.
FIG. 1 illustrates an exemplary lead configuration for Cardiac Resynchronization Therapy (CRT);
FIG. 2 illustrates a method for guiding a cardiac interventional therapy according to an embodiment of the present invention;
FIG. 3 illustrates a detailed framework of the computational heart model according to an embodiment of the present invention;
FIG. 4 illustrates exemplary results of generating a patient-specific anatomical heart model;
FIG. 5 illustrates an example of interactive virtual CRT;
FIG. 6 illustrates an exemplary left ventricle (LV)-only outcome maps for QRS duration and electrical axis cardiac parameters;
FIG. 7 illustrates an exemplary 2D unfolded outcome map for the LV;
FIG. 8 illustrates an exemplary outcome map overlaid on an interventional image; and
The present invention relates to visualization of cardiac changes under various pacing conditions for guiding lead placement in interventional therapies. Embodiments of the present invention are described herein to give a visual understanding
of the methods for visualization of cardiac changes and guiding interventional therapies. A digital image is often composed of digital representations of one or more objects (or shapes). The digital representation of an object is often described herein
in terms of identifying and manipulating the objects. Such manipulations are virtual manipulations accomplished in the memory or other circuitry/hardware of a computer system. Accordingly, is to be understood that embodiments of the present invention
may be performed within a computer system using data stored within the computer system or available through a network system.
Cardiac Resynchronization Therapy (CRT) involves placement of an advanced pacemaker to pace the myocardium at the right ventricle (RV) endocardium and the left ventricle (LV) epicardium, with given pacing intervals, to resynchronize the cardiac
motion and hence re-establish the pumping efficiency of the heart. CRT is a well-established therapy for heart failure. FIG. 1 illustrates an exemplary lead configuration for CRT including the LV lead, the RV lead, and a left atrium (LA) lead.
FIG. 2 illustrates a method for guiding a cardiac interventional therapy according to an embodiment of the present invention. The method of FIG. 2 is described herein for guiding CRT, but the present invention is not limited thereto and the
method of FIG. 2 may be similarly applied to ablative treatments for cardiac arrhythmias or any other cardiac intervention procedure that is applied at one or more specific locations of the heart. Furthermore, the method of FIG. 2 may be similarly
applied to interventional therapies performed at specific locations in other organs (e.g., liver tumor ablation therapy) using computational models for other organs (e.g., liver heat capacity), as well. As illustrated in FIG. 2, at step 202, clinical
data and medical image data of a patient are received. The clinical data can include electrophysiology data such as ECG and/or pressure data such as invasive catheter measurements or pressure cuff measurements. The medical image data may be a sequence
of 3D medical images acquired over at least one complete heart cycle. In advantageous embodiments of the present invention, the medical images can be MRI images, CT images and/or ultrasound images, but the present invention is not necessarily limited to
these imaging modalities. The medical images may be received directly from a medical imaging device, such as an MR, CT or ultrasound scanner, or the medical images may be received by loading stored medical images of a patient.
At step 204, a patient-specific computational heart model is generated based on the medical images and the clinical data of the patient. The patient-specific computational heart model is generated by solving an inverse problem to adjust
parameters of the computational heart model such that simulated parameters, such as heart motion, ejection fraction, etc., output by the patient-specific computational heart model match the clinical data and medical images observed for the patient.
FIG. 3 illustrates a detailed framework of the computational heart model according to an embodiment of the present invention. As illustrated in FIG. 3, the patient-specific heart model includes four main models representing cardiac anatomy,
cardiac electrophysiology, cardiac biomechanics, and cardiac hemodynamics. These models are connected as follows. Starting from a medical image 302 of the heart (MRI, CT, Ultrasound), at step 310, a detailed patient-specific anatomical model 312 of the
patient's heart is generated. In one exemplary embodiment, a detailed anatomical model is generated of the two ventricles only, but the present invention is not limited thereto. Arteries and atria are modeled as boundary conditions of the system. The
anatomical model 312 comprises the bi-ventricular geometry, the orientation of myocyte fibers, and any information that varies spatially such as action potential duration, scars, etc. The anatomical model 312 is then given as input to the Heart Model
subsystem (step 320), which will compute myocardium motion over time according to three main sub-parts: the cardiac electrophysiology model 322, the cardiac biomechanics model 324, and the cardiac hemodynamics model 326. Each of these models of cardiac
function can be fitted to patient data using inverse problem approaches and all available preoperative data such as dynamic images 306, ECG 304, invasive pressure measurements 308, etc. As a result, the simulated parameters 330 match the clinical
observations. Once the model is personalized, cardiac therapy, such as CRT, can be simulated by pacing the model at various locations for instance. Clinical parameters are computed from the simulation to quantify the predicted patient response to the
therapy. Depending on the application, the model can include of one, several or all of the above-mentioned sub-parts.
In an exemplary embodiment, in order to generate the patient-specific anatomical model 312, anatomical models of the left ventricle (LV) and right ventricle (RV) are extracted from the medical images. Although it is possible to extract models
for all of the heart chambers, in an advantageous embodiment, only the geometry of LV and RV are explicitly modeled. If the application targets atrial anatomy, at least one of the atria can be segmented and used as computational domain. For each of the
LV and the RV, the heart chamber segmentation can be formulated as a two-step learning problem: anatomical structure localization and boundary delineation. Marginal space learning (MSL) can be used to apply machine learning to 3D object detection. The
idea of MSL is not to learn a monolithic classifier directly in the full parameter space of similarity transformations but to incrementally learn classifiers on marginal spaces. In particular, the detection of each heart chamber can be split into three
problems: position estimation, position-orientation estimation, and position-orientation-scale estimation. A separate classifier is trained based on annotated training data for each of these estimation problems. The classifiers in the lower dimensional
marginal spaces are used to prune the searching space efficiently. This object localization stage results in an estimated transformation (position, orientation, and scale) of the object (e.g., heart chamber). After automatic object localization, the
mean shape model of the object is aligned with the estimated transformation to get a rough estimate of the object shape. The shape is then deformed locally to fit the object boundary using an active shape model (ASM) and a machine learning based
boundary detector. Additional details regarding MSL-based heart chamber segmentation are described in U.S. Pat. No. 7,916,919, issued Mar. 29, 2011, and entitled "System and Method for Segmenting Chambers of a Heart in a Three Dimensional Image",
United States Published Patent Application No. 2010/0040272, and United States Published Patent Application No. 2012/0022843, which are incorporated herein by reference.
The patient-specific LV and RV models are then fused into a single anatomical model of the bi-ventricular myocardium. In particular, the LV and RV anatomies are fused into a single volumetric mesh representation, on which vertices are tagged
into surface zones according to the underlying anatomy. According to an advantageous implementation, tetrahedral elements can be used to accurately represent the details of the bi-ventricular anatomy. Spatial information can then be mapped onto the
anatomical model of the bi-ventricular myocardium. Spatial information, such as scars, border zones, and fibrosis can be identified in images, such as late enhancement MRI. For example, the spatial information may be automatically identified using
trained classifiers or may be manually identified by a clinician. The spatial information is mapped onto the tetrahedral mesh representing the bi-ventricular myocardium. This information is important to simulate the electrical wave around scars, in
particular for wave-reentry assessment and to correctly capture impaired cardiac mechanics due to ill-functioning or dead cells.
A model of fiber orientation can be automatically calculated based on the patient-specific geometry. In an advantageous implementation, the model of fiber orientation can be automatically calculated using a rule-based approach. A generic model
of myocardium fiber architecture that includes fiber and fiber sheets is computed. A rule-based strategy is followed to generate the fiber architecture to cover the entire bi-ventricular myocardium from apex to valves. Below the basal plane, which is
identified automatically using point correspondences of the initial triangulations of the anatomical model, the fiber elevation angle .alpha., i.e. the angle with respect to the short axis plane, varies linearly across the myocardium, e.g., from -70 on
the epicardium to +70 on the endocardium (values that can be defined by the user). Similarly, the sheet direction, which is defined by the angle .beta. with respect to the outward transmural axis, varies transmurally, e.g., from +45 on the epicardium
to -45 on the endocardium (values that can be defined by the user). .alpha. and .beta. are computed for each point of the volumetric bi-ventricular myocardium mesh between the apex and basal plane based on the geodesic distance to the endocardia and
epicardia identified by the facet tags: .alpha.=(d.sub.epi.alpha..sub.endo+d.sub.endo.alpha..sub.epi)/(d.sub.endo- +d.sub.epi), where d.sub.epi, d.sub.endo, .alpha..sub.epi, and .alpha..sub.endo are the distances and angles at the endocardium and
epicardium, respectively. The fiber and sheet orientations are then fixed around each valve. In particular, fibers are longitudinal around the aortic valve and tangential around the mitral, tricuspid, and pulmonary valves, and sheet normals are
oriented towards the barycenter of the valves. The local orthonormal basis is then interpolated from the basal plane to the valve, first by following the myocardium surface, then throughout the myocardium thickness. Additional details regarding
generating the patient-specific anatomical heart model are described in U.S. Published Patent Application No. 2013/0197881 and U.S. Published Patent Application No. 2015/0042464, which are incorporated herein in their entirety by reference. In another
embodiment, the fibers are directly measured in-vivo, in the patient, using MRI diffusion tensor imaging. In another embodiment, an atlas of cardiac fibers could be employed. FIG. 4 illustrates exemplary results of generating a patient-specific
anatomical heart model. As illustrated in FIG. 4, image 400 shows segmentation results for segmenting the left and right ventricles in a medical image, image 410 shows the LV and RV fused into a single volumetric mesh, and image 420 shows the model of
fiber orientation added to the volumetric mesh.
Once the patient-specific anatomical model 312 is generated from the medical image data of the patient, the patient-specific computational heart model 320 is generated by determining patient-specific parameters for the cardiac electrophysiology
model 322, cardiac biomechanics model 324, and cardiac hemodynamics model 326, using the patient-specific anatomical model 312 as input.
The first biophysical phenomenon that needs to be simulated is the cardiac electrophysiology (322), as it commands the cardiac contraction. According to an exemplary embodiment, the cardiac electrophysiology model 322 can compute action
potentials over a time period at each of the nodes of the patient-specific anatomical model using a Lattice-Boltzmann method for Electrophysiology (LBM-EP). In an exemplary embodiment, a Mitchell-Schaeffer cellular model can be employed as a tradeoff
between model fidelity and complexity, but other cellular models can be employed as well. Patient-specific parameters of the cardiac electrophysiological model, such as an electrical diffusivity and an action potential duration at each of the nodes of
the patient-specific anatomical model, can be estimated by solving an inverse problem based on fitting the measured electrophysiology data of the patient. For example, the parameters of the cardiac electrophysiology model can be personalized using ECG
measurements, invasive endocardial mappings, or body surface mappings of the patient. Additional details regarding the cardiac electrophysiology model and estimating patient-specific parameters of the cardiac electrophysiology model are described in
U.S. Published Patent Application No. 2013/0226542 and U.S. Published Patent Application No. 2015/0042646, which are incorporated herein in their entirety by reference.
The cardiac hemodynamics model 326 can compute cardiac blood flow and pressure over a period of time using a lumped model of intra-ventricular pressure that varies according to models of arterial compliance and atrial contraction. A 3-element
Windkessel (WK) model can be used as a model of arterial compliance. The cardiac hemodynamics model can be personalized by estimating patient-specific arterial compliance parameters (e.g., patient-specific parameters of the WK model of arterial
compliance) from flow and pressure data of the patient measured from medical images (e.g., Doppler Ultrasound or 2D/4D PC-MRI), cuff pressure measurements, or invasive pressure catheter measurements. Additional details regarding the cardiac hemodynamics
model and estimating patient-specific parameters of the cardiac hemodynamics model are described in U.S. Published Patent Application No. 2013/0197881, which is incorporated herein in its entirety by reference.
The cardiac biomechanics model 324 simulates movement of the heart chambers over a time period based on the simulated electrophysiology by the cardiac electrophysiology model 322 and the simulated blood flow and pressure by the cardiac
hemodynamics model 326. The cardiac biomechanics model 324 can compute the movement of the heart using the finite element method to solve the dynamics equation Mu+Cu+Ku=f.sub.a+f.sub.p+f.sub.b. u, u, and u denote nodal accelerations, velocities, and
displacements, respectively. M, K, and C are the mass, stiffness, and Rayleigh damping matrices, respectively. Pressure force f.sub.p is calculated from the cardiac hemodynamics model. The active stress f.sub.a is generated by the myocytes and is
computed using a phenomenological model of myocyte contraction, which is mainly governed by the maximum active stress parameter .sigma..sub.0. Boundary conditions f.sub.b capture heart attachment at the valve plane through springs and pericardium
constraints. The orthotropic Holzapfel-Ogden (HO) constitutive law can be used, among others, defined by the stress-strain energy:
.PSI..times..times.e.function..times..times..function.e.function..times..- times..function.e.times..times. ##EQU00001##
Subscripts f and s denote fiber and sheet directions, I.sub.1, I.sub.4{f,s}, and I.sub.8fs are invariants of the deformation tensor, and a, b, a.sub.{f,s,fs}, and b.sub.{f,s,fs} are patient-specific tissue parameters. Equations 1 can be solved
according to multiplicative Jacobian energy decomposition (MJED) formulation or the Total Lagrangian Explicit Dynamics (TLED) method for efficient computations. The cardiac biomechanics model 324 can be personalized by estimating patient-specific
biomechanical tissue parameters (e.g., stiffness and maximum active stress) based on observed heart movement in the medical images of the patient using inverse problem techniques.
As described herein and illustrated in FIG. 3, the patient-specific computational model of heart function (320) includes a personalized cardiac electrophysiology model 322, a personalized cardiac hemodynamics model 326, and a personalized
cardiac biomechanics model 324. It is to be understood that these models are modular and in various embodiments of the present invention, different combinations of these models may be used. For example, the method of FIG. 2 may be implemented using
only the electrophysiology model or the electrophysiology model coupled with the biomechanics model, instead of the full computational heart model 320 shown in FIG. 3. The simulated cardiac parameters (e.g., ejection fraction, QRS duration, QT duration,
etc.) used to determine the optimal lead placement and pacing protocol may vary based on which components of the computations heart model are used.
As described above, the components (cardiac electrophysiology model 322, cardiac hemodynamics model 326, and the cardiac biomechanics model 324) of the computational heart model 320 are personalized using inverse modeling based on the medical
image data and the clinical data of the patient. According to an alternative embodiment, machine learning techniques can be used to directly estimate the personalized model parameters from the medical image data and the clinical data. In addition, the
components of the computational heart model can be updated using interventional information (e.g., EP data, pressure measurements, an ablation region, etc.) obtained during the intervention procedure. In this case, outcome maps generated using the
computational heart model in steps 206 and 208 can be re-generated using the updated computational heart model for finer guidance.
Returning to FIG. 2, at step 206, virtual CRT is performed using a plurality of lead locations and pacing protocols. Based on the patient-specific computational model of cardiac function, virtual CRT is performed by placing LV and RV leads on
the mesh model (i.e., patient-specific anatomical model) and applying electrical stimuli at these locations following a given pacing protocol. The pacing protocol is a set of parameters that controls the pacing at the various leads. In an exemplary
implementation, the pacing protocol includes the frequency of each pacing, the interval between LV and RV pacing, and the interval between the first pacing lead and spontaneous depolarization. Pacing protocols can be input following the standard
nomenclature of CRT devices (DDD, DDT, etc.) or can be freely set by a user. A library of devices can also be made available with their available protocols as listed in the manufacturer description for ease of use. The electrical stimuli are applied as
an additional localized electrical current in the cardiac electrophysiology model. The overall, patient-specific cardiac computational model (e.g., cardiac electrophysiology, cardiac hemodynamics, and cardiac biomechanics) is then recalculated in order
to simulate the cardiac function over a plurality of time frames, and parameters of the cardiac function, such as QRS duration, electrical axis, QT duration, ejection fraction, stroke volume, and myocardium strain, are derived from the computational
Virtual CRT can be performed interactively. In this case, the user selects on the 3D anatomical heart model the positions of the LV lead and the RV lead and the pacing protocol (e.g., using a user input device, such as a mouse, touch-screen,
etc.), and the computational heart model is re-calculated to simulate the cardiac function under the input pacing conditions. Changes in the cardiac electrophysiology, cardiac electrocardiogram (ECG), cardiac motion, and cardiac hemodynamics are
computed and displayed to the user. FIG. 5 illustrates an example of interactive virtual CRT. As illustrated in FIG. 5, image 500 shows an LV lead position 502 and an RV lead position 504 input interactively by a user in the patient-specific anatomical
heart model 506, and image 510 shows an ECG calculated by the patient-specific computational heart model in the virtual CRT performed using the input LV and RV lead positons 502 and 504.
Although virtual CRT can be performed interactively, in an advantageous embodiment of the present invention, virtual CRT is performed automatically at a plurality of pacing locations for the pacing leads using one or more pacing protocols at
each of the pacing locations in order to automatically generate outcome maps that visualize the optimal outcomes for the CRT in order to guide lead placement in the CRT procedure. In an advantageous implementation, the method automatically iterates the
virtual CRT simulation over a series of LV and RV position configurations, and for each lead position configuration automatically iterates the virtual CRT simulation over a series of pacing protocols in order to calculate acute cardiac parameters (e.g.,
QRS duration, QT duration, ejection fraction, stroke volume, etc.) resulting from each of the pacing protocols at each of the lead position configurations. For example, the virtual CRT can be iterated over all points on the LV and RV in the
patient-specific anatomical model. In a possible implementation, one of the LV lead position and the RV lead position can be varied while the other one is fixed in order to independently find optimal LV and RV pacing locations. In another possible
implementation, both the LV lead position and the RV lead position can be varied, such that virtual CRT is performed for all possible combinations of LV and RV lead positions using each of the pacing protocols. Although this method is described herein
as using multiple pacing protocols, it is also possible that a single pacing protocol be used and the lead positions varied for the single pacing protocol.
Returning to FIG. 2, at step 208, outcome maps are generated based on the virtual CRT. The outcome maps visualize optimal outcomes for lead placement based on one or more cardiac parameters resulting from the virtual CRT simulations. The
outcome maps can be used for planning the CRT procedure or guiding the clinician towards the optimal lead location in the CRT procedure.
In an exemplary embodiment, LV-only outcome maps and RV-only outcome maps can be generated and used to respectively guide the lead placement for the LV and RV leads. LV-only outcome maps can be generated as follows: Let X.sub.RV be a fixed RV
lead position. For each possible position of the LV lead, X.sub.LV, at the LV epicardium and LV endocardium, the cardiac function (e.g., electrophysiology, hemodynamics, and biomechanics) and cardiac parameters .theta.=[.theta..sub.1, .theta..sub.2, . . . , .theta..sub.n] (e.g., .theta..sub.1=QRS duration, .theta..sub.2=QT duration, .theta..sub.3=ejection fraction) are computed using each of pacing protocol P.sub.j. Although the cardiac parameters of QRS duration, QT duration, and ejection fraction
are used in the above example, the present invention is not limited thereto and other cardiac parameters such electrical axis, stroke volume, myocardium strain, etc. can be used as well. At each position X.sub.LV on the 3D mesh (i.e., patient-specific
anatomical model), the optimal .theta..sub.i and the associated pacing protocol P.sub.j resulting in the optimal .theta..sub.i are reported for each of the cardiac parameters. The optimal value for each cardiac parameter can be defined based on clinical
knowledge. For example, the optimal value for the QRS duration can be the smallest QRS value and the optimal value for the ejection fraction can be the maximum ejection fraction. The optimal value for each cardiac parameter can also be set by a user. In an exemplary embodiment, a respective outcome map is generated for each cardiac parameter .theta..sub.i by assigning a color to each LV point on the 3D mesh to represent the optimal outcome value for the cardiac parameter .theta..sub.i in that point,
given a pacing protocol. This results in a respective map that visualizes the optimal LV lead placement to achieve the best possible value for each cardiac parameter. For each cardiac parameter .theta..sub.i, a respective outcome map can also be
generated to visualize the pacing protocol P.sub.j associated with the optimal cardiac parameter value at each LV position on the 3D mesh. For example, various colors can be assigned to represent various possible pacing protocols, such that a color
assigned to each LV point on the 3D mesh represents the pacing protocol that results in the optimal value for the cardiac parameter. Accordingly 2.times.n outcome maps can be generated to visualize the optimal cardiac parameter value and optimal pacing
protocol for each of the n cardiac parameters. Each map can be displayed (e.g., on a display device of a computer system) individually on demand or simultaneously on several views. RV-only outcome maps can be generated similarly to the LV-only outcome
maps by fixing the LV lead position X.sub.LV, varying the RV lead location X.sub.RV, and computing the cardiac function and cardiac parameters at each RV lead location X.sub.RV using each of the pacing protocols.
In another embodiment, bi-lead outcome maps can be generated by varying both the LV lead position and the RV lead position. Bi-lead outcome maps can be generated as follows. For every X.sub.LV lead position on the LV epicardium and
endocardium, the RV lead position X.sub.RV can be varied and for each X.sub.RV the parameters .theta. are computed for all available pacing protocols P. At each LV point X.sub.LV, the optimal parameter value (for each of the cardiac parameters) and the
associated RV lead location X.sub.RV and pacing protocol P.sub.j that result in the optimal parameter value are reported. All maps can be displayed simultaneously. Alternatively, the user can display a map related to one parameter .theta..sub.i in
particular and obtain the optimal LV lead position X.sub.LV in terms of optimal outcome for that parameter, the associated RV lead position, and the associated pacing protocol. The bi-lead protocol can be extended to a multi-lead procedure in order to
guide lead placement with devices with additional leads.
FIG. 6 illustrates an exemplary LV-only outcome maps for QRS duration and electrical axis cardiac parameters. The outcome maps in FIG. 6 were generated using a fixed pacing protocol. As illustrated in FIG. 6, image 600 shows an outcome map
representing how optimal QRS varies based on LV lead position for a given pacing protocol. Image 610 shows an outcome map representing how optimal electrical axis varies based on LV lead position for a given pacing protocol. Similar maps can be
obtained for other cardiac parameters calculated by the patient-specific computational model of heart function.
As described above, the outcome maps can be generated using 3D color maps. In another possible embodiment, the outcome maps can be displayed as flat 2D bull's-eye plots. In this case, a polar map is calculated from the 3D mesh for a particular
outcome map and used to draw the corresponding bull's-eye plot. In another possible embodiment, different protocols can be tested over a sampled number of points on the LV and the results can be mapped on a square as if the LV were "unfolded". FIG. 7
illustrates an exemplary 2D unfolded outcome map 700 for the LV. The outcome map 700 arranges the sampled points of the LV from an inferior (front) to posterior (back) direction on the horizontal axis and a direction from the LV apex to the LV base on
the vertical axis. The color of each point of the outcome map 700 represents the electrical axis value at that lead position for LV pacing at 70 ms. A similar map can be displayed for each pacing protocol and/or each cardiac parameter.
As described above virtual CRT is performed at each possible LV lead location and/or each possible RV lead location on the patient-specific anatomical heart model. In alternative embodiment, only a subset of nodes can be tested and
interpolation methods can be used to generate the outcome maps based on the optimal cardiac parameters computed at each of the subset of nodes. Linear, non-linear, or sparse interpolation methods can be employed to perform the interpolation.
Returning to FIG. 2, at step 210, the outcome maps are overlaid on interventional images in an intervention procedure. The 3D outcome maps generated in step 208 can be calculated automatically offline prior to an intervention procedure using
pre-operative medical image data or can be calculated online if interventional 3D imaging of the heart (e.g., MRI, ultrasound, DynaCT) is available. The 3D outcome maps can then be registered to interventional images acquired using an interventional
imaging system (e.g., angiography system, interventional MRI, EP mapping system, or ultrasound) using standard registration techniques. This registration can be performed in real time or near real time during the intervention procedure. By overlaying
one or more 3D outcome maps on the interventional images, the leads and catheters are visualized along the with the outcome map, which enables the user to target the region with optimal outcome as identified by the virtual interventions performed using
patient-specific computational model. FIG. 8 illustrates an exemplary outcome map overlaid on an interventional image. As shown in FIG. 8, a QRS outcome map 800 is overlaid on an angiogram showing the CRT leads 802 and 804. A patient-specific model of
the coronary arteries can also be overlaid. The model of the coronary arteries can be obtained from pre-operative CT scans or intra-operative images, such as DynaCT.
The method of FIG. 2 can be similarly applied to cardiac ablation therapy. At each position of the 3D heart mesh, a virtual ablation is performed and cardiac parameters are recalculated using the patient-specific computational model of heart
function. Virtual ablations using various ablation protocols can be performed at each position. The outcome maps are then generated as described above and can be used to guide cardiac ablation therapy.
Embodiments of the present invention can be similarly applied to other organs (e.g., liver), for example in the case of ablative therapy of abdominal tumors. A patient-specific anatomical model (mesh) of the organ is generated and for each
point of the mesh, a virtual ablation is applied and outcome parameters are calculated. A patient-specific computational model of heat diffusion for the organ may be generated and used to perform the virtual ablations. For example, a computational
model of liver heat diffusion is described in U.S. Publication No. 2014/0136174 and International Publication No. WO/2014/133924 A1, the disclosures of which are incorporated herein by reference in their entirety. The outcome parameters can include
parameters such as recurrence risk, lesion coverage, temperature, etc. The resulting outcome map can then be overlaid on interventional images to guide the ablation procedure.
Although a bi-ventricular anatomical heart model is described above, embodiments of the present invention can also be implemented by modeling the whole heart including all four chambers. For example, this may allow for more precise modeling of
In a possible embodiment, the method of FIG. 2 can be implemented using a cloud-based system to generate the outcome maps. In this embodiment, computations needed to perform the virtual interventions and calculate the outcome maps can be
distributed on a cloud system. In a particular implementation, only the mesh, computational model, and outcome map parameters need to be distributed, thus ensuring privacy and data security.
The above-described methods for visualizing cardiac changes under various pacing conditions and guiding an intervention procedure can be implemented on a computer using well-known computer processors, memory units, storage devices, computer
software, and other components. A high-level block diagram of such a computer is illustrated in FIG. 9. Computer 902 contains a processor 904, which controls the overall operation of the computer 902 by executing computer program instructions which
define such operation. The computer program instructions may be stored in a storage device 912 (e.g., magnetic disk) and loaded into memory 910 when execution of the computer program instructions is desired. Thus, the steps of the methods of FIGS. 2
and 3 may be defined by the computer program instructions stored in the memory 910 and/or storage 912 and controlled by the processor 904 executing the computer program instructions. An image acquisition device 920, such as a CT scanning device, C-arm
image acquisition device, MR scanning device, Ultrasound device, etc., can be connected to the computer 902 to input image data to the computer 902. It is possible to implement the image acquisition device 920 and the computer 902 as one device. It is
also possible that the image acquisition device 920 and the computer 902 communicate wirelessly through a network. In a possible embodiment, the computer 902 may be located remotely with respect to the image acquisition device 920 and may perform the
method steps as part of a server or cloud based service. The computer 902 also includes one or more network interfaces 906 for communicating with other devices via a network. The computer 902 also includes other input/output devices 908 that enable
user interaction with the computer 902 (e.g., display, keyboard, mouse, speakers, buttons, etc.). Such input/output devices 908 may be used in conjunction with a set of computer programs as an annotation tool to annotate volumes received from the image
acquisition device 920. One skilled in the art will recognize that an implementation of an actual computer could contain other components as well, and that FIG. 9 is a high level representation of some of the components of such a computer for
The above-described methods for medical image synthesis may be implemented using computers operating in a client-server relationship. Typically, in such a system, the client computers are located remotely from the server computer and interact
via a network. The client-server relationship may be defined and controlled by computer programs running on the respective client and server computers.
Previous Patent US 9,589,378 | Next Patent US 9,589,380 File A Patent Application Protect your idea -- Don't let