Incorporation of a left ventricle finite element model defining infarction into the XCAT imaging phantom

Development of physiologically-informed computational coronary artery plaques for use in virtual imaging trials

BACKGROUND

As a leading cause of death, worldwide, cardiovascular disease is of great clinical importance. Among cardiovascular diseases, coronary artery disease (CAD) is a key contributor, and it is the attributed cause of death for 10% of all deaths annually. The prevalence of CAD is commensurate with the rise in new medical imaging technologies intended to aid in its diagnosis and treatment. The necessary clinical trials required to validate and optimize these technologies require a large cohort of carefully controlled patients, considerable time to complete, and can be prohibitively expensive. A safer, faster, less expensive alternative is using virtual imaging trials (VITs), utilizing virtual patients or phantoms combined with accurate computer models of imaging devices.

PURPOSE

In this work, we develop realistic, physiologically-informed models for coronary plaques for application in cardiac imaging VITs.

METHODS

Histology images of plaques at micron-level resolution were used to train a deep convolutional generative adversarial network (DC-GAN) to create a library of anatomically variable plaque models with clinical anatomical realism. The stability of each plaque was evaluated by finite element analysis (FEA) in which plaque components and vessels were meshed as volumes, modeled as specialized tissues, and subjected to the range of normal coronary blood pressures. To demonstrate the utility of the plaque models, we combined them with the whole-body XCAT computational phantom to perform initial simulations comparing standard energy-integrating detector (EID) CT with photon-counting detector (PCD) CT.

RESULTS

Our results show the network is capable of generating realistic, anatomically variable plaques. Our simulation results provide an initial demonstration of the utility of the generated plaque models as targets to compare different imaging devices.

CONCLUSIONS

Vast, realistic, and variable CAD pathologies can be generated to incorporate into computational phantoms for VITs. There they can serve as a known truth from which to optimize and evaluate cardiac imaging technologies quantitatively.

SPECT and CT misregistration reduction in [99mTc]Tc-MAA SPECT/CT for precision liver radioembolization treatment planning

PURPOSE

Respiration and body movement induce misregistration between static [99mTc]Tc-MAA SPECT and CT, causing lung shunting fraction (LSF) and tumor-to-normal liver ratio (TNR) errors for 90Y radioembolization planning. We aim to alleviate the misregistration between [99mTc]Tc-MAA SPECT and CT using two registration schemes on simulation and clinical data.

METHODS

In the simulation study, 70 XCAT phantoms were modeled. The SIMIND Monte Carlo program and OS-EM algorithm were used for projection generation and reconstruction, respectively. Low-dose CT (LDCT) at end-inspiration was simulated for attenuation correction (AC), lungs and liver segmentation, while contrast-enhanced CT (CECT) was simulated for tumor and perfused liver segmentation. In the clinical study, 16 patient data including [99mTc]Tc-MAA SPECT/LDCT and CECT with observed SPECT and CT mismatch were analyzed. Two liver-based registration schemes were studied: SPECT registered to LDCT/CECT and vice versa. Mean count density (MCD) of different volumes-of-interest (VOIs), normalized mutual information (NMI), LSF, TNR, and maximum injected activity (MIA) based on the partition model before and after registration were compared. Wilcoxon signed-rank test was performed.

RESULTS

In the simulation study, compared to before registration, registrations significantly reduced estimation errors of MCD of all VOIs, LSF (Scheme 1: - 100.28%, Scheme 2: - 101.59%), and TNR (Scheme 1: - 7.00%, Scheme 2: - 5.67%), as well as MIA (Scheme 1: - 3.22%, Scheme 2: - 2.40%). In the clinical study, Scheme 1 reduced 33.68% LSF and increased 14.75% TNR, while Scheme 2 reduced 38.88% LSF and increased 6.28% TNR compared to before registration. One patient may change from 90Y radioembolization untreatable to treatable and other patients may change the MIA up to 25% after registration. NMI between SPECT and CT was significantly increased after registrations in both studies.

CONCLUSION

Registration between static [99mTc]Tc-MAA SPECT and corresponding CTs is feasible to reduce their spatial mismatch and improve dosimetric estimation. The improvement of LSF is larger than TNR. Our method can potentially improve patient selection and personalized treatment planning for liver radioembolization.

Technical Note: Respiratory impacts on static and respiratory gated 99m Tc-MAA SPECT/CT for liver radioembolization- A simulation study

PURPOSE

We aimed to evaluate respiratory impacts on static and respiratory gated (RG) ^99m Tc-MAA SPECT in terms of respiratory motion (RM) blur, attenuation correction (AC) and volume-of-interest (VOI) segmentation on lung shunt faction (LSF) and tumor-to-normal liver ratio (TNR) estimation for liver radioembolization therapy planning.

METHODS

The XCAT phantom was used to simulate a population of 300 phantoms, modelling various anatomical variations, tumor characteristics, respiratory motion amplitudes, LSFs and TNRs. One hundred and twenty noisy projections of average activity maps near end-expiration (End-EX) and whole respiratory cycle were simulated analytically, modeling attenuation and geometric collimator-detector-response (GCDR). The OS-EM algorithm was employed for reconstruction, modeling AC and GCDR. RM effect was evaluated for static SPECT, while AC and VOI mismatch effects were investigated independently and together for static and RG SPECT utilizing one gate, i.e., End-EX. LSF and TNR errors were measured based on the ground truth. Lesions with different characteristics were also investigated for static and RG SPECT.

RESULTS

RM overestimates LSF and underestimates TNR. The VOI mismatch caused the largest errors in both RG and static SPECT for LSF and TNR estimation, reaching 160% and -52% correspondingly with extremely mismatched VOIs for RG SPECT, even larger than those for static SPECT. With matched AC and VOIs, RG SPECT has better performance than static SPECT. Larger TNR errors are associated with tumors of smaller sizes and higher TNR for static SPECT.

CONCLUSIONS

The VOI segmentation mismatch has a stronger impact, followed by RM and AC in static ^99m Tc-MAA SPECT/CT. This effect is more pronounced for RG SPECT. When VOI masks are derived from a matched CT, RG SPECT is generally superior to static SPECT for LSF and TNR estimation. The performance of RG SPECT could be worse than static SPECT when a mismatched CT is used for segmentation. This article is protected by copyright. All rights reserved.

Novel Methodology for Measuring Regional Myocardial Efficiency

Our approach differs from the usual global measure of cardiac efficiency by using PET/MRI to measure efficiency of small pieces of cardiac tissue whose limiting size is equal to the spatial resolution of the PET scanner. We initiated a dynamic cardiac PET study immediately prior to the injection of 15.1 mCi of 11C-acetate acquiring data for 25 minutes while simultaneously acquiring MRI cine data. 1) A 3D finite element (FE) biomechanical model of the imaged heart was constructed by utilizing nonrigid deformable image registration to alter the Dassault Systèmes FE Living Heart Model (LHM) to fit the geometry in the cardiac MRI cine data. The patient specific FE cardiac model with estimates of stress, strain, and work was transformed into PET/MRI format. 2) A 1-tissue compartment model was used to calculate wash-in (K1) and the linear portion of the decay in the PET 11C-acetate time activity curve (TAC) was used to calculate the wash-out k2(mono) rate constant. K1 was used to calculate blood flow and k2(mono) was used to calculate myocardial volume oxygen consumption ( MVO2 ). 3) Estimates of stress and strain were used to calculate Myocardial Equivalent Minute Work ( MEMW ) and Cardiac Efficiency = MEMW/MVO2 was then calculated for 17 tissue segments of the left ventricle. The global MBF was 0.96 ± 0.15 ml/min/gm and MVO2 ranged from 8 to 17 ml/100gm/min. Six central slices of the MRI cine data provided a range of MEMW of 0.1 to 0.4 joules/gm/min and a range of Cardiac Efficiency of 6 to 18%.

Innovations in Computer Technologies Have Impacted Radiation Dosimetry Through Anatomically Realistic Phantoms and Fast Monte Carlo Simulations

Radiological physics principles have not changed in the past 60 y when computer technologies advanced exponentially. The research field of anatomical modeling for the purpose of radiation dose calculations has experienced an explosion in activity in the past two decades. Such an exciting advancement is due to the feasibility of creating three-dimensional geometric details of the human anatomy from tomographic imaging and of performing Monte Carlo radiation transport simulations on increasingly fast and cheap personal computers. The advent of a new type of high-performance computing hardware in recent years-graphics processing units-has made it feasible to carry out time-consuming Monte Carlo calculations at near real-time speeds. This paper introduces the history of three generations of computational human phantoms (the stylized medical internal radiation dosimetry-type phantoms, the voxelized tomographic phantoms, and the boundary representation deformable phantoms) and new development of the graphics processing unit-based Monte Carlo radiation dose calculations. Examples are given for research projects performed by my students in applying computational phantoms and a new Monte Carlo code, ARCHER, to problems in radiation protection, imaging, and radiotherapy. Finally, the paper discusses challenges and future opportunities for research.

Incorporation of the Living Heart Model into the 4D XCAT Phantom for Cardiac Imaging Research

The 4D extended cardiac-torso (XCAT) phantom has provided a valuable tool to study the effects of anatomy and motion on medical images, especially cardiac motion. One limitation of the XCAT was that it did not have a physiological basis which to realistically simulate variations in cardiac function. In this work, we incorporate into the XCAT anatomy the four-chamber FE Living Heart Model (LHM) developed by the Living Heart Project (LHP). The LHM represents the state of the art in cardiac FE simulation because of its ability to accurately replicate the biomechanical motion of the entire heart and its variations. We create a new series of 4D phantoms capable of simulating patients with varying body sizes and shapes; cardiac positions, orientations, and dynamics. While extendable to other imaging modalities and technologies, our goal is to use the FE-enhanced XCAT models to investigate the optimal use of computed tomography (CT) for the evaluation of coronary artery disease (CAD). With the ability to simulate realistic, predictive, patient quality 4D imaging data, the enhanced XCAT models will enable optimization studies to identify the most promising systems or system parameters for further clinical validation.

A finite element model of myocardial infarction using a composite material approach

Computational models are effective tools to study cardiac mechanics under normal and pathological conditions. They can be used to gain insight into the physiology of the heart under these conditions while they are adaptable to computer assisted patient-specific clinical diagnosis and therapeutic procedures. Realistic cardiac mechanics models incorporate tissue active/passive response in conjunction with hyperelasticity and anisotropy. Conventional formulation of such models leads to mathematically-complex problems usually solved by custom-developed non-linear finite element (FE) codes. With a few exceptions, such codes are not available to the research community. This article describes a computational cardiac mechanics model developed such that it can be implemented using off-the-shelf FE solvers while tissue pathologies can be introduced in the model in a straight-forward manner. The model takes into account myocardial hyperelasticity, anisotropy, and active contraction forces. It follows a composite tissue modeling approach where the cardiac tissue is decomposed into two major parts: background and myofibers. The latter is modelled as rebars under initial stresses mimicking the contraction forces. The model was applied in silico to study the mechanics of infarcted left ventricle (LV) of a canine. End-systolic strain components, ejection fraction, and stress distribution attained using this LV model were compared quantitatively and qualitatively to corresponding data obtained from measurements as well as to other corresponding LV mechanics models. This comparison showed very good agreement.

PET Image Reconstruction and Deformable Motion Correction Using Unorganized Point Clouds

Quantitative positron emission tomography imaging often requires correcting the image data for deformable motion. With cyclic motion, this is traditionally achieved by separating the coincidence data into a relatively small number of gates, and incorporating the inter-gate image transformation matrices into the reconstruction algorithm. In the presence of non-cyclic deformable motion, this approach may be impractical due to a large number of required gates. In this paper, we propose an alternative approach to iterative image reconstruction with correction for deformable motion, wherein unorganized point clouds are used to model the imaged objects in the image space, and motion is corrected for explicitly by introducing a time-dependence into the point coordinates. The image function is represented using constant basis functions with finite support determined by the boundaries of the Voronoi cells in the point cloud. We validate the quantitative accuracy and stability of the proposed approach by reconstructing noise-free and noisy projection data from digital and physical phantoms. The point-cloud-based maximum likelihood expectation maximization (MLEM) and one-pass list-mode ordered-subset expectation maximization (OSEM) algorithms are validated. The results demonstrate that images reconstructed using the proposed method are quantitatively stable, with noise and convergence properties comparable to image reconstruction based on the use of rectangular and radially-symmetric basis functions.

Quantitative analysis of hypertrophic myocardium using diffusion tensor magnetic resonance imaging

Systemic hypertension is a causative factor in left ventricular hypertrophy (LVH). This study is motivated by the potential to reverse or manage the dysfunction associated with structural remodeling of the myocardium in this pathology. Using diffusion tensor magnetic resonance imaging, we present an analysis of myocardial fiber and laminar sheet orientation in ex vivo hypertrophic (6 SHR) and normal (5 WKY) rat hearts using the covariance of the diffusion tensor. First, an atlas of normal cardiac microstructure was formed using the WKY b0 images. Then, the SHR and WKY b0 hearts were registered to the atlas. The acquired deformation fields were applied to the SHR and WKY heart tensor fields followed by the preservation of principal direction (PPD) reorientation strategy. A mean tensor field was then formed from the registered WKY tensor images. Calculating the covariance of the registered tensor images about this mean for each heart, the hypertrophic myocardium exhibited significantly increased myocardial fiber derangement ([Formula: see text]) with a mean dispersion of 38.7 deg, and an increased dispersion of the laminar sheet normal ([Formula: see text]) of 54.8 deg compared with 34.8 deg and 51.8 deg, respectively, in the normal hearts. Results demonstrate significantly altered myocardial fiber and laminar sheet structure in rats with hypertensive LVH.

Image reconstruction in higher dimensions: myocardial perfusion imaging of tracer dynamics with cardiac motion due to deformation and respiration

Myocardial perfusion imaging (MPI) using slow rotating large field of view cameras requires spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration. In vivo, MPI contains additional degrees of freedom involving unavoidable motion of the heart due to quasiperiodic beating and the effects of respiration, which can severely degrade the quality of the images. This work develops a technique for a single photon emission computed tomography (SPECT) that reconstructs the distribution of the radiotracer concentration in the myocardium using a tensor product of different sets of basis functions that approximately describe the spatiotemporal variation of the radiotracer concentration and the motion of the heart. In this study the temporal B-spline basis functions are chosen to reflect the dynamics of the radiotracer, while the intrinsic deformation and the extrinsic motion of the heart are described by a product of a discrete set of Gaussian basis functions. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it deforms due to cardiac beating, and is displaced due to respiratory motion. These results are compared with the conventional 4D-spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. The higher dimensional reconstruction method proposed here improves bias, yet the signal-to-noise ratio (SNR) decreases slightly due to redistribution of the counts over the cardiac-respiratory gates. Additionally, there is a trade-off between the number of gates and the number of projections per gate to achieve high contrast images.

The direct incorporation of perfusion defect information to define ischemia and infarction in a finite element model of the left ventricle

This paper describes the process in which complex lesion geometries (specified by computer generated perfusion defects) are incorporated in the description of nonlinear finite element (FE) mechanical models used for specifying the motion of the left ventricle (LV) in the 4D extended cardiac torso (XCAT) phantom to simulate gated cardiac image data. An image interrogation process was developed to define the elements in the LV mesh as ischemic or infarcted based upon the values of sampled intensity levels of the perfusion maps. The intensity values were determined for each of the interior integration points of every element of the FE mesh. The average element intensity levels were then determined. The elements with average intensity values below a user-controlled threshold were defined as ischemic or infarcted depending upon the model being defined. For the infarction model cases, the thresholding and interrogation process were repeated in order to define a border zone (BZ) surrounding the infarction. This methodology was evaluated using perfusion maps created by the perfusion cardiac-torso (PCAT) phantom an extension of the 4D XCAT phantom. The PCAT was used to create 3D perfusion maps representing 90% occlusions at four locations (left anterior descending (LAD) segments 6 and 9, left circumflex (LCX) segment 11, right coronary artery (RCA) segment 1) in the coronary tree. The volumes and shapes of the defects defined in the FE mechanical models were compared with perfusion maps produced by the PCAT. The models were incorporated into the XCAT phantom. The ischemia models had reduced stroke volume (SV) by 18-59 ml. and ejection fraction (EF) values by 14-50% points compared to the normal models. The infarction models, had less reductions in SV and EF, 17-54 ml. and 14-45% points, respectively. The volumes of the ischemic/infarcted regions of the models were nearly identical to those volumes obtained from the perfusion images and were highly correlated (R² = 0.99).

Continuous flow left ventricular pump support and its effect on regional left ventricular wall stress: finite element analysis study

Left ventricular assist device (LVAD) support unloads left ventricular (LV) pressure and volume and decreases wall stress. This study investigated the effect of systematic LVAD unloading on the 3-dimensional myocardial wall stress by employing finite element models containing layered fiber structure, active contractility, and passive stiffness. The HeartMate II(®) (Thoratec, Inc., Pleasanton, CA) was used for LV unloading. The model geometries and hemodynamic conditions for baseline (BL) and LVAD support (LVsupport) were acquired from the Penn State mock circulatory cardiac simulator. Myocardial wall stress of BL was compared with that of LVsupport at 8,000, 9,000, 10,000 RPM, providing mean pump flow (Q(mean)) of 2.6, 3.2, and 3.7 l/min, respectively. LVAD support was more effective at unloading during diastole as compared to systole. Approximately 40, 50, and 60% of end-diastolic wall stress reduction were achieved at Q(mean) of 2.6, 3.2, and 3.7 l/min, respectively, as compared to only a 10% reduction of end-systolic wall stress at Q(mean) of 3.7 l/min. In addition, there was a stress concentration during systole at the apex due to the cannulation and reduced boundary motion. This modeling study can be used to further understand optimal unloading, pump control, patient management, and cannula design.

An exponential growth of computational phantom research in radiation protection, imaging, and radiotherapy: a review of the fifty-year history

Radiation dose calculation using models of the human anatomy has been a subject of great interest to radiation protection, medical imaging, and radiotherapy. However, early pioneers of this field did not foresee the exponential growth of research activity as observed today. This review article walks the reader through the history of the research and development in this field of study which started some 50 years ago. This review identifies a clear progression of computational phantom complexity which can be denoted by three distinct generations. The first generation of stylized phantoms, representing a grouping of less than dozen models, was initially developed in the 1960s at Oak Ridge National Laboratory to calculate internal doses from nuclear medicine procedures. Despite their anatomical simplicity, these computational phantoms were the best tools available at the time for internal/external dosimetry, image evaluation, and treatment dose evaluations. A second generation of a large number of voxelized phantoms arose rapidly in the late 1980s as a result of the increased availability of tomographic medical imaging and computers. Surprisingly, the last decade saw the emergence of the third generation of phantoms which are based on advanced geometries called boundary representation (BREP) in the form of Non-Uniform Rational B-Splines (NURBS) or polygonal meshes. This new class of phantoms now consists of over 287 models including those used for non-ionizing radiation applications. This review article aims to provide the reader with a general understanding of how the field of computational phantoms came about and the technical challenges it faced at different times. This goal is achieved by defining basic geometry modeling techniques and by analyzing selected phantoms in terms of geometrical features and dosimetric problems to be solved. The rich historical information is summarized in four tables that are aided by highlights in the text on how some of the most well-known phantoms were developed and used in practice. Some of the information covered in this review has not been previously reported, for example, the CAM and CAF phantoms developed in 1970s for space radiation applications. The author also clarifies confusion about 'population-average' prospective dosimetry needed for radiological protection under the current ICRP radiation protection system and 'individualized' retrospective dosimetry often performed for medical physics studies. To illustrate the impact of computational phantoms, a section of this article is devoted to examples from the author's own research group. Finally the author explains an unexpected finding during the course of preparing for this article that the phantoms from the past 50 years followed a pattern of exponential growth. The review ends on a brief discussion of future research needs (a supplementary file '3DPhantoms.pdf' to figure 15 is available for download that will allow a reader to interactively visualize the phantoms in 3D).

MRXCAT: Realistic numerical phantoms for cardiovascular magnetic resonance

BACKGROUND

Computer simulations are important for validating novel image acquisition and reconstruction strategies. In cardiovascular magnetic resonance (CMR), numerical simulations need to combine anatomical information and the effects of cardiac and/or respiratory motion. To this end, a framework for realistic CMR simulations is proposed and its use for image reconstruction from undersampled data is demonstrated.

METHODS

The extended Cardiac-Torso (XCAT) anatomical phantom framework with various motion options was used as a basis for the numerical phantoms. Different tissue, dynamic contrast and signal models, multiple receiver coils and noise are simulated. Arbitrary trajectories and undersampled acquisition can be selected. The utility of the framework is demonstrated for accelerated cine and first-pass myocardial perfusion imaging using k-t PCA and k-t SPARSE.

RESULTS

MRXCAT phantoms allow for realistic simulation of CMR including optional cardiac and respiratory motion. Example reconstructions from simulated undersampled k-t parallel imaging demonstrate the feasibility of simulated acquisition and reconstruction using the presented framework. Myocardial blood flow assessment from simulated myocardial perfusion images highlights the suitability of MRXCAT for quantitative post-processing simulation.

CONCLUSION

The proposed MRXCAT phantom framework enables versatile and realistic simulations of CMR including breathhold and free-breathing acquisitions.

Left ventricular finite element model bounded by a systemic circulation model

A series of models were developed in which a circulatory system model was coupled to an existing series of finite element (FE) models of the left ventricle (LV). The circulatory models were used to provide realistic boundary conditions for the LV models. This was developed for the JSim analysis package and was composed of a systemic arterial, capillary, and venous system in a closed loop with a varying elastance LV and left atria to provide the driving pressures and flows matching those of the FE model. Three coupled models were developed, a normal LV under normotensive aortic loading (116/80 mm Hg), a mild hypertension (137/89 mm Hg) model, and a moderate hypertension model (165/100 mm Hg). The initial step in the modeling analysis was that the circulation was optimized to the end-diastolic pressure and volume values of the LV model. The cardiac FE models were then optimized to the systolic pressure/volume characteristics of the steady-state JSim circulatory model solution. Comparison of the stress predictions for the three models indicated that the mild hypertensive case produced a 21% increase in the average fiber stress levels, and the moderate hypertension case had a 36% increase in average stress. The circulatory work increased by 18% and 43% over that of the control for the mild and moderate hypertensive cases, respectively.

A Comparison of Hyperelastic Warping of PET Images with Tagged MRI for the Analysis of Cardiac Deformation

The objectives of the following research were to evaluate the utility of a deformable image registration technique known as hyperelastic warping for the measurement of local strains in the left ventricle through the analysis of clinical, gated PET image datasets. Two normal human male subjects were sequentially imaged with PET and tagged MRI imaging. Strain predictions were made for systolic contraction using warping analyses of the PET images and HARP based strain analyses of the MRI images. Coefficient of determination R (2) values were computed for the comparison of circumferential and radial strain predictions produced by each methodology. There was good correspondence between the methodologies, with R (2) values of 0.78 for the radial strains of both hearts and from an R (2) = 0.81 and R (2) = 0.83 for the circumferential strains. The strain predictions were not statistically different (P ≤ 0.01). A series of sensitivity results indicated that the methodology was relatively insensitive to alterations in image intensity, random image noise, and alterations in fiber structure. This study demonstrated that warping was able to provide strain predictions of systolic contraction of the LV consistent with those provided by tagged MRI Warping.

A magnetic resonance software simulator for the evaluation of myocardial deformation estimation

This paper proposes a methodology to design a physiologically realistic computer simulator of images of the left ventricle myocardium based on a patient-specific biomechanical model. The simulator takes a magnetic resonance image of a given patient at end diastole, uses a manual segmentation of that image to model the geometry of the myocardium and sets the parameters of the constitutive model used for biomechanical simulation according to a regional labeling of the contractility of the myocardium for that patient. The simulated deformations are used to warp the magnetic resonance dataset throughout the cardiac cycle to generate different image modalities. The simulator is validated by quantifying its ability to model actual deformations in a set of patients affected by an acute myocardial infarction. Specifically a high correlation has been encountered between the ejection fraction derived from the simulated end systolic deformation of the myocardium and the myocardium segmented from actual data. Additionally, most of the parameters that describe the simulated deformation compare well with reported values. Overall, the simulator is intended as a testbed for extensive comparisons of myocardial motion tracking methods due to its ability to relate the impaired myocardial function with the associated ventricular remodeling, a novel contribution in the literature of cardiac image simulators.

Generation of synthetic but visually realistic time series of cardiac images combining a biophysical model and clinical images

We propose a new approach for the generation of synthetic but visually realistic time series of cardiac images based on an electromechanical model of the heart and real clinical 4-D image sequences. This is achieved by combining three steps. The first step is the simulation of a cardiac motion using an electromechanical model of the heart and the segmentation of the end diastolic image of a cardiac sequence. We use biophysical parameters related to the desired condition of the simulated subject. The second step extracts the cardiac motion from the real sequence using nonrigid image registration. Finally, a synthetic time series of cardiac images corresponding to the simulated motion is generated in the third step by combining the motion estimated by image registration and the simulated one. With this approach, image processing algorithms can be evaluated as we know the ground-truth motion underlying the image sequence. Moreover, databases of visually realistic images of controls and patients can be generated for which the underlying cardiac motion and some biophysical parameters are known. Such databases can open new avenues for machine learning approaches.

Four-Dimensional Image Reconstruction Strategies in Cardiac-Gated and Respiratory-Gated PET Imaging

Cardiac and respiratory movements pose significant challenges to image quality and quantitative accuracy in PET imaging. Cardiac and/or respiratory gating attempt to address this issue, but instead lead to enhanced noise levels. Direct four-dimensional (4D) PET image reconstruction incorporating motion compensation has the potential to minimize noise amplification while removing considerable motion blur. A wide-ranging choice of such techniques is reviewed in this work. Future opportunities and the challenges facing the adoption of 4D PET reconstruction and its role in basic and clinical research are also discussed.

Weighted least-squares finite element method for cardiac blood flow simulation with echocardiographic data

As both fluid flow measurement techniques and computer simulation methods continue to improve, there is a growing need for numerical simulation approaches that can assimilate experimental data into the simulation in a flexible and mathematically consistent manner. The problem of interest here is the simulation of blood flow in the left ventricle with the assimilation of experimental data provided by ultrasound imaging of microbubbles in the blood. The weighted least-squares finite element method is used because it allows data to be assimilated in a very flexible manner so that accurate measurements are more closely matched with the numerical solution than less accurate data. This approach is applied to two different test problems: a flexible flap that is displaced by a jet of fluid and blood flow in the porcine left ventricle. By adjusting how closely the simulation matches the experimental data, one can observe potential inaccuracies in the model because the simulation without experimental data differs significantly from the simulation with the data. Additionally, the assimilation of experimental data can help the simulation capture certain small effects that are present in the experiment, but not modeled directly in the simulation.

Modeling and Representation of Human Hearts for Volumetric Measurement

This paper investigates automatic construction of a three-dimensional heart model from a set of medical images, represents it in a deformable shape, and uses it to perform volumetric measurements. This not only significantly improves its reliability and accuracy but also makes it possible to derive valuable novel information, like various assessment and dynamic volumetric measurements. The method is based on a flexible model trained from hundreds of patient image sets by a genetic algorithm, which takes advantage of complete segmentation of the heart shape to form a geometrical heart model. For an image set of a new patient, an interpretation scheme is used to obtain its shape and evaluate some important parameters. Apart from automatic evaluation of traditional heart functions, some new information of cardiovascular diseases may be recognized from the volumetric analysis.

Clinical benefits of integrating cardiac and vascular models

INTRODUCTION

Recent advances in computational methods and medical imaging techniques have enabled non-invasive exploration of cardiovascular pathologies, from cardiac level to complex arterial networks. The potential of cardiac and vascular modeling in guiding and monitoring therapies could be further extended through the integration of the two systems.

AREAS COVERED

This review includes advances in methods for cardiac electromechanics and vascular flow simulations. The results of a literature search depicting the state of the art in cardiac and vascular modeling are reviewed. The paper goes on to address the benefits and challenges of combined cardiovascular modeling, highlighting the relevance of specific cardiovascular features and implementation. Various alternative approaches and insights on future directions are presented and analyzed with respect to their applicability to clinical practice.

EXPERT OPINION

The article has emerged from the exploration of currently available cardiac and vascular mathematical tools and their corresponding clinical application. The summarized analysis suggests that future efforts should be aimed at developing more accurate and patient-specific mathematical models integrating cardiac and vascular functions to enhance the knowledge of cardiovascular pathologies.