Quantifying changes in shoulder orientation between the prone and supine positions from magnetic resonance imaging.

BACKGROUND: Predicting breast tissue motion using biomechanical models can provide navigational guidance during breast cancer treatment procedures. These models typically do not account for changes in posture between procedures. Difference in shoulder position can alter the shape of the pectoral muscles and breast. A greater understanding of the differences in the shoulder orientation between prone and supine could improve the accuracy of breast biomechanical models. METHODS: 19 landmarks were placed on the sternum, clavicle, scapula, and humerus of the shoulder girdle in prone and supine breast MRIs (N = 10). These landmarks were used in an optimization framework to fit subject-specific skeletal models and compare joint angles of the shoulder girdle between these positions. FINDINGS: The mean Euclidean distance between joint locations from the fitted skeletal model and the manually identified joint locations was 15.7 mm ± 2.7 mm. Significant differences were observed between prone and supine. Compared to supine position, the shoulder girdle in the prone position had the lateral end of the clavicle in more anterior translation (i.e., scapula more protracted) (P < 0.05), the scapula in more protraction (P < 0.01), the scapula in more upward rotation (associated with humerus elevation) (P < 0.05); and the humerus more elevated (P < 0.05) for both the left and right sides. INTERPRETATION: Shoulder girdle orientation was found to be different between prone and supine. These differences would affect the shape of multiple pectoral muscles, which would affect breast shape and the accuracy of biomechanical models.

Vision for the 12 LABOURS Digital Twin Platform.

Clinical translation of personalised computational physiology workflows and digital twins can revolutionise healthcare by providing a better understanding of an individual's physiological processes and any changes that could lead to serious health consequences. However, the lack of common infrastructure for developing these workflows and digital twins has hampered the realisation of this vision. The Auckland Bioengineering Institute's 12 LABOURS project aims to address these challenges by developing a Digital Twin Platform to enable researchers to develop and personalise computational physiology models to an individual's health data in clinical workflows. This will allow clinical trials to be more efficiently conducted to demonstrate the efficacy of these personalised clinical workflows. We present a demonstration of the platform's capabilities using publicly available data and an existing automated computational physiology workflow developed to assist clinicians with diagnosing and treating breast cancer. We also demonstrate how the platform facilitates the discovery and exploration of data and the presentation of workflow results as part of clinical reports through a web portal. Future developments will involve integrating the platform with health systems and remote-monitoring devices such as wearables and implantables to support home-based healthcare. Integrating outputs from multiple workflows that are applied to the same individual's health data will also enable the generation of their personalised digital twin.Clinical Relevance- The proposed 12 LABOURS Digital Twin Platform will enable researchers to 1) more efficiently conduct clinical trials to assess the efficacy of their computational physiology workflows and support the clinical translation of their research; 2) reuse primary and derived data from these workflows to generate novel workflows; and 3) generate personalised digital twins by integrating the outputs of different computational physiology workflows.

Quinapril treatment curtails decline of global longitudinal strain and mechanical function in hypertensive rats.

BACKGROUND: Left ventricular (LV) global longitudinal strain (GLS) has been proposed as an early imaging biomarker of cardiac mechanical dysfunction. OBJECTIVE: To assess the impact of angiotensin-converting enzyme (ACE) inhibitor treatment of hypertensive heart disease on LV GLS and mechanical function. METHODS: The spontaneously hypertensive rat (SHR) model of hypertensive heart disease ( n  = 38) was studied. A subset of SHRs received quinapril (TSHR, n  = 16) from 3 months (mo). Wistar Kyoto rats (WKY, n  = 13) were used as controls. Tagged cardiac MRI was performed using a 4.7 T Varian preclinical scanner. RESULTS: The SHRs had significantly lower LV ejection fraction (EF) than the WKYs at 3 mo (53.0 ±â€Š1.7% vs. 69.6 ±â€Š2.1%, P  < 0.05), 14 mo (57.0 ±â€Š2.5% vs. 74.4 ±â€Š2.9%, P  < 0.05) and 24 mo (50.1 ±â€Š2.4% vs. 67.0 ±â€Š2.0%, P  < 0.01). At 24 mo, ACE inhibitor treatment was associated with significantly greater LV EF in TSHRs compared to untreated SHRs (64.2 ±â€Š3.4% vs. 50.1 ±â€Š2.4%, P  < 0.01). Peak GLS magnitude was significantly lower in SHRs compared with WKYs at 14 months (7.5% ±â€Š0.4% vs. 9.9 ±â€Š0.8%, P  < 0.05). At 24 months, Peak GLS magnitude was significantly lower in SHRs compared with both WKYs (6.5 ±â€Š0.4% vs. 9.7 ±â€Š1.0%, P  < 0.01) and TSHRs (6.5 ±â€Š0.4% vs. 9.6 ±â€Š0.6%, P  < 0.05). CONCLUSIONS: ACE inhibitor treatment curtails the decline in global longitudinal strain in hypertensive rats, with the treatment group exhibiting significantly greater LV EF and GLS magnitude at 24 mo compared with untreated SHRs.

Non-contact quantification of aortic stenosis and mitral regurgitation using carotid waveforms from skin displacements.

Objective. Early diagnosis of heart problems is essential for improving patient prognosis.Approach. We created a non-contact imaging system that calculates the vessel-induced deformation of the skin to estimate the carotid artery pressure displacement waveforms. We present a clinical study of the system in patients (n= 27) with no underlying condition, aortic stenosis (AS), or mitral regurgitation (MR).Main results. Displacement waveforms were compared to aortic catheter pressures in the same patients. The morphologies of the pressure and displacement waveforms were found to be similar, and pulse wave analysis metrics, such as our modified reflection indices (RI) and waveform duration proportions, showed no significant differences. Compared with the control group, AS patients displayed a greater proportion of time to peak (p= 0.026 andp= 0.047 for catheter and displacement, respectively), whereas augmentation index (AIx)was greater for the displacement waveform only (p= 0.030). The modified RI for MR (p= 0.047 andp= 0.004 for catheter and displacement, respectively) was lower than in the controls. AS and MR were also significantly different for the proportion of time to peak (p= 0.018 for the catheter measurements), RI (p= 0.045 andp= 0.002 for the catheter and displacement, respectively), and AIx (p= 0.005 for the displacement waveform).Significance. These findings demonstrate the ability of our system to provide insights into cardiac conditions and support further development as a diagnostic/telehealth-based screening tool.

Markov modeling of phase singularity interaction effects in human atrial and ventricular fibrillation.

Atrial and ventricular fibrillation (AF/VF) are characterized by the repetitive regeneration of topological defects known as phase singularities (PSs). The effect of PS interactions has not been previously studied in human AF and VF. We hypothesized that PS population size would influence the rate of PS formation and destruction in human AF and VF, due to increased inter-defect interaction. PS population statistics were studied in computational simulations (Aliev-Panfilov), human AF and human VF. The influence of inter-PS interactions was evaluated by comparison between directly modeled discrete-time Markov chain (DTMC) transition matrices of the PS population changes, and M/M/∞ birth-death transition matrices of PS dynamics, which assumes that PS formations and destructions are effectively statistically independent events. Across all systems examined, PS population changes differed from those expected with M/M/∞. In human AF and VF, the formation rates decreased slightly with PS population when modeled with the DTMC, compared with the static formation rate expected through M/M/∞, suggesting new formations were being inhibited. In human AF and VF, the destruction rates increased with PS population for both models, with the DTMC rate increase exceeding the M/M/∞ estimates, indicating that PS were being destroyed faster as the PS population grew. In human AF and VF, the change in PS formation and destruction rates as the population increased differed between the two models. This indicates that the presence of additional PS influenced the likelihood of new PS formation and destruction, consistent with the notion of self-inhibitory inter-PS interactions.

Correcting bias in cardiac geometries derived from multimodal images using spatiotemporal mapping.

Cardiovascular imaging studies provide a multitude of structural and functional data to better understand disease mechanisms. While pooling data across studies enables more powerful and broader applications, performing quantitative comparisons across datasets with varying acquisition or analysis methods is problematic due to inherent measurement biases specific to each protocol. We show how dynamic time warping and partial least squares regression can be applied to effectively map between left ventricular geometries derived from different imaging modalities and analysis protocols to account for such differences. To demonstrate this method, paired real-time 3D echocardiography (3DE) and cardiac magnetic resonance (CMR) sequences from 138 subjects were used to construct a mapping function between the two modalities to correct for biases in left ventricular clinical cardiac indices, as well as regional shape. Leave-one-out cross-validation revealed a significant reduction in mean bias, narrower limits of agreement, and higher intraclass correlation coefficients for all functional indices between CMR and 3DE geometries after spatiotemporal mapping. Meanwhile, average root mean squared errors between surface coordinates of 3DE and CMR geometries across the cardiac cycle decreased from 7 ± 1 to 4 ± 1 mm for the total study population. Our generalised method for mapping between time-varying cardiac geometries obtained using different acquisition and analysis protocols enables the pooling of data between modalities and the potential for smaller studies to leverage large population databases for quantitative comparisons.

Roadmap for an imaging and modelling paediatric study in rural NZ.

Our study methodology is motivated from three disparate needs: one, imaging studies have existed in silo and study organs but not across organ systems; two, there are gaps in our understanding of paediatric structure and function; three, lack of representative data in New Zealand. Our research aims to address these issues in part, through the combination of magnetic resonance imaging, advanced image processing algorithms and computational modelling. Our study demonstrated the need to take an organ-system approach and scan multiple organs on the same child. We have pilot tested an imaging protocol to be minimally disruptive to the children and demonstrated state-of-the-art image processing and personalized computational models using the imaging data. Our imaging protocol spans brain, lungs, heart, muscle, bones, abdominal and vascular systems. Our initial set of results demonstrated child-specific measurements on one dataset. This work is novel and interesting as we have run multiple computational physiology workflows to generate personalized computational models. Our proposed work is the first step towards achieving the integration of imaging and modelling improving our understanding of the human body in paediatric health and disease.

Super-resolution 4D flow MRI to quantify aortic regurgitation using computational fluid dynamics and deep learning.

Changes in cardiovascular hemodynamics are closely related to the development of aortic regurgitation (AR), a type of valvular heart disease. Metrics derived from blood flows are used to indicate AR onset and evaluate its severity. These metrics can be non-invasively obtained using four-dimensional (4D) flow magnetic resonance imaging (MRI), where accuracy is primarily dependent on spatial resolution. However, insufficient resolution often results from limitations in 4D flow MRI and complex aortic regurgitation hemodynamics. To address this, computational fluid dynamics simulations were transformed into synthetic 4D flow MRI data and used to train a variety of neural networks. These networks generated super-resolution, full-field phase images with an upsample factor of 4. Results showed decreased velocity error, high structural similarity scores, and improved learning capabilities from previous work. Further validation was performed on two sets of in vivo 4D flow MRI data and demonstrated success in de-noising flow images. This approach presents an opportunity to comprehensively analyse AR hemodynamics in a non-invasive manner.

Automated model calibration with parallel MCMC: Applications for a cardiovascular system model.

Computational physiological models continue to increase in complexity, however, the task of efficiently calibrating the model to available clinical data remains a significant challenge. One part of this challenge is associated with long calibration times, which present a barrier for the routine application of model-based prediction in clinical practice. Another aspect of this challenge is the limited available data for the unique calibration of complex models. Therefore, to calibrate a patient-specific model, it may be beneficial to verify that task-specific model predictions have acceptable uncertainty, rather than requiring all parameters to be uniquely identified. We have developed a pipeline that reduces the set of fitting parameters to make them structurally identifiable and to improve the efficiency of a subsequent Markov Chain Monte Carlo (MCMC) analysis. MCMC was used to find the optimal parameter values and to determine the confidence interval of a task-specific prediction. This approach was demonstrated on numerical experiments where a lumped parameter model of the cardiovascular system was calibrated to brachial artery cuff pressure, echocardiogram volume measurements, and synthetic cerebral blood flow data that approximates what can be obtained from 4D-flow MRI data. This pipeline provides a cerebral arterial pressure prediction that may be useful for determining the risk of hemorrhagic stroke. For a set of three patients, this pipeline successfully reduced the parameter set of a cardiovascular system model from 12 parameters to 8-10 structurally identifiable parameters. This enabled a significant ( > 4 × ) efficiency improvement in determining confidence intervals on predictions of pressure compared to performing a naive MCMC analysis with the full parameter set. This demonstrates the potential that the proposed pipeline has in helping address one of the key challenges preventing clinical application of such models. Additionally, for each patient, the MCMC approach yielded a 95% confidence interval on systolic blood pressure prediction in the middle cerebral artery smaller than ±10 mmHg (±1.3 kPa). The proposed pipeline exploits available high-performance computing parallelism to allow straightforward automation for general models and arbitrary data sets, enabling automated calibration of a parameter set that is specific to the available clinical data with minimal user interaction.

The inspection paradox: An important consideration in the evaluation of rotor lifetimes in cardiac fibrillation.

Background and Objective: Renewal theory is a statistical approach to model the formation and destruction of phase singularities (PS), which occur at the pivots of spiral waves. A common issue arising during observation of renewal processes is an inspection paradox, due to oversampling of longer events. The objective of this study was to characterise the effect of a potential inspection paradox on the perception of PS lifetimes in cardiac fibrillation. Methods: A multisystem, multi-modality study was performed, examining computational simulations (Aliev-Panfilov (APV) model, Courtmanche-Nattel model), experimentally acquired optical mapping Atrial and Ventricular Fibrillation (AF/VF) data, and clinically acquired human AF and VF. Distributions of all PS lifetimes across full epochs of AF, VF, or computational simulations, were compared with distributions formed from lifetimes of PS existing at 10,000 simulated commencement timepoints. Results: In all systems, an inspection paradox led towards oversampling of PS with longer lifetimes. In APV computational simulations there was a mean PS lifetime shift of +84.9% (95% CI, ± 0.3%) (p < 0.001 for observed vs overall), in Courtmanche-Nattel simulations of AF +692.9% (95% CI, ±57.7%) (p < 0.001), in optically mapped rat AF +374.6% (95% CI, ± 88.5%) (p = 0.052), in human AF mapped with basket catheters +129.2% (95% CI, ±4.1%) (p < 0.05), human AF-HD grid catheters 150.8% (95% CI, ± 9.0%) (p < 0.001), in optically mapped rat VF +171.3% (95% CI, ±15.6%) (p < 0.001), in human epicardial VF 153.5% (95% CI, ±15.7%) (p < 0.001). Conclusion: Visual inspection of phase movies has the potential to systematically oversample longer lasting PS, due to an inspection paradox. An inspection paradox is minimised by consideration of the overall distribution of PS lifetimes.

A governing equation for rotor and wavelet number in human clinical ventricular fibrillation: Implications for sudden cardiac death.

BACKGROUND: Ventricular fibrillation (VF) is characterized by multiple wavelets and rotors. No equation to predict the number of rotors and wavelets observed during fibrillation has been validated in human VF. OBJECTIVE: The purpose of this study was to test the hypothesis that a single equation derived from a Markov M/M/∞ birth-death process could predict the number of rotors and wavelets occurring in human clinical VF. METHODS: Epicardial induced VF (256-electrode) recordings obtained from patients undergoing cardiac surgery were studied (12 patients; 62 epochs). Rate constants for phase singularity (PS) (which occur at the pivot points of rotors) and wavefront (WF) formation and destruction were derived by fitting distributions to PS and WF interformation and lifetimes. These rate constants were combined in an M/M/∞ governing equation to predict the number of PS and WF in VF episodes. Observed distributions were compared to those predicted by the M/M/∞ equation. RESULTS: The M/M/∞ equation accurately predicted average PS and WF number and population distribution, demonstrated in all epochs. Self-terminating episodes of VF were distinguished from VF episodes requiring termination by a trend toward slower PS destruction, slower rates of PS formation, and a slower mixing rate of the VF process, indicated by larger values of the second largest eigenvalue modulus of the M/M/∞ birth-death matrix. The longest-lasting PS (associated with rotors) had shorter interactivation time intervals compared to shorter-lasting PS lasting <150 ms (∼1 PS rotation in human VF). CONCLUSION: The M/M/∞ equation explains the number of wavelets and rotors observed, supporting a paradigm of VF based on statistical fibrillatory dynamics.

MITEA: A dataset for machine learning segmentation of the left ventricle in 3D echocardiography using subject-specific labels from cardiac magnetic resonance imaging.

Segmentation of the left ventricle (LV) in echocardiography is an important task for the quantification of volume and mass in heart disease. Continuing advances in echocardiography have extended imaging capabilities into the 3D domain, subsequently overcoming the geometric assumptions associated with conventional 2D acquisitions. Nevertheless, the analysis of 3D echocardiography (3DE) poses several challenges associated with limited spatial resolution, poor contrast-to-noise ratio, complex noise characteristics, and image anisotropy. To develop automated methods for 3DE analysis, a sufficiently large, labeled dataset is typically required. However, ground truth segmentations have historically been difficult to obtain due to the high inter-observer variability associated with manual analysis. We address this lack of expert consensus by registering labels derived from higher-resolution subject-specific cardiac magnetic resonance (CMR) images, producing 536 annotated 3DE images from 143 human subjects (10 of which were excluded). This heterogeneous population consists of healthy controls and patients with cardiac disease, across a range of demographics. To demonstrate the utility of such a dataset, a state-of-the-art, self-configuring deep learning network for semantic segmentation was employed for automated 3DE analysis. Using the proposed dataset for training, the network produced measurement biases of -9 ± 16 ml, -1 ± 10 ml, -2 ± 5 %, and 5 ± 23 g, for end-diastolic volume, end-systolic volume, ejection fraction, and mass, respectively, outperforming an expert human observer in terms of accuracy as well as scan-rescan reproducibility. As part of the Cardiac Atlas Project, we present here a large, publicly available 3DE dataset with ground truth labels that leverage the higher resolution and contrast of CMR, to provide a new benchmark for automated 3DE analysis. Such an approach not only reduces the effect of observer-specific bias present in manual 3DE annotations, but also enables the development of analysis techniques which exhibit better agreement with CMR compared to conventional methods. This represents an important step for enabling more efficient and accurate diagnostic and prognostic information to be obtained from echocardiography.

Efficient Ventricular Parameter Estimation Using AI-Surrogate Models.

The onset and progression of pathological heart conditions, such as cardiomyopathy or heart failure, affect its mechanical behaviour due to the remodelling of the myocardial tissues to preserve its functional response. Identification of the constitutive properties of heart tissues could provide useful biomarkers to diagnose and assess the progression of disease. We have previously demonstrated the utility of efficient AI-surrogate models to simulate passive cardiac mechanics. Here, we propose the use of this surrogate model for the identification of myocardial mechanical properties and intra-ventricular pressure by solving an inverse problem with two novel AI-based approaches. Our analysis concluded that: (i) both approaches were robust toward Gaussian noise when the ventricle data for multiple loading conditions were combined; and (ii) estimates of one and two parameters could be obtained in less than 9 and 18 s, respectively. The proposed technique yields a viable option for the translation of cardiac mechanics simulations and biophysical parameter identification methods into the clinic to improve the diagnosis and treatment of heart pathologies. In addition, the proposed estimation techniques are general and can be straightforwardly translated to other applications involving different anatomical structures.

Passive myocardial mechanical properties: meaning, measurement, models.

Passive mechanical tissue properties are major determinants of myocardial contraction and relaxation and, thus, shape cardiac function. Tightly regulated, dynamically adapting throughout life, and affecting a host of cellular functions, passive tissue mechanics also contribute to cardiac dysfunction. Development of treatments and early identification of diseases requires better spatio-temporal characterisation of tissue mechanical properties and their underlying mechanisms. With this understanding, key regulators may be identified, providing pathways with potential to control and limit pathological development. Methodologies and models used to assess and mimic tissue mechanical properties are diverse, and available data are in part mutually contradictory. In this review, we define important concepts useful for characterising passive mechanical tissue properties, and compare a variety of in vitro and in vivo techniques that allow one to assess tissue mechanics. We give definitions of key terms, and summarise insight into determinants of myocardial stiffness in situ. We then provide an overview of common experimental models utilised to assess the role of environmental stiffness and composition, and its effects on cardiac cell and tissue function. Finally, promising future directions are outlined.

Systematic Comparison of Left Ventricular Geometry Between 3D-Echocardiography and Cardiac Magnetic Resonance Imaging.

Aims: Left ventricular (LV) volumes estimated using three-dimensional echocardiography (3D-echo) have been reported to be smaller than those measured using cardiac magnetic resonance (CMR) imaging, but the underlying causes are not well-understood. We investigated differences in regional LV anatomy derived from these modalities and related subsequent findings to image characteristics. Methods and Results: Seventy participants (18 patients and 52 healthy participants) were imaged with 3D-echo and CMR (<1 h apart). Three-dimensional left ventricular models were constructed at end-diastole (ED) and end-systole (ES) from both modalities using previously validated software, enabling the fusion of CMR with 3D-echo by rigid registration. Regional differences were evaluated as mean surface distances for each of the 17 American Heart Association segments, and by comparing contours superimposed on images from each modality. In comparison to CMR-derived models, 3D-echo models underestimated LV end-diastolic volume (EDV) by -16 ± 22, -1 ± 25, and -18 ± 24 ml across three independent analysis methods. Average surface distance errors were largest in the basal-anterolateral segment (11-15 mm) and smallest in the mid-inferoseptal segment (6 mm). Larger errors were associated with signal dropout in anterior regions and the appearance of trabeculae at the lateral wall. Conclusions: Fusion of CMR and 3D-echo provides insight into the causes of volume underestimation by 3D-echo. Systematic signal dropout and differences in appearances of trabeculae lead to discrepancies in the delineation of LV geometry at anterior and lateral regions. A better understanding of error sources across modalities may improve correlation of clinical indices between 3D-echo and CMR.

Quantifying optical anisotropy in soft tissue membranes using Mueller matrix imaging.

SIGNIFICANCE: A non-destructive technique for accurately characterizing the spatial distribution of optical properties of soft tissue membranes may give improved outcomes in many tissue engineering applications. AIM: This study aimed to develop a non-destructive macroscopic imaging technique that is sensitive to optical anisotropy, typical of fibrous components in soft tissue membranes, and can address some of the difficulties caused by the complex turbid nature of these tissues. APPROACH: A near-infrared Mueller matrix imaging polarimeter employing logarithm decomposition was developed and used to conduct transmission measurements of all the polarization properties across the full thickness of bovine pericardium tissue. RESULTS: The full Mueller matrix was measured across a 70 mm × 70 mm sample of calf bovine pericardium and revealed significant retardance (linear and circular) and depolarization in this tissue. Regions with a uniform axis of optical anisotropy were identified. Mueller matrix imaging demonstrated that the exhibited circular retardance was sufficient to lead to possible misinterpretation of apparent fiber orientation when using conventional polarization imaging techniques for such tissues. CONCLUSIONS: Mueller matrix imaging can identify regional distributions of optical anisotropy in calf bovine pericardium. This new capability is a promising development in non-destructive imaging for tissue selection.

Insights From Computational Modeling Into the Contribution of Mechano-Calcium Feedback on the Cardiac End-Systolic Force-Length Relationship.

In experimental studies on cardiac tissue, the end-systolic force-length relation (ESFLR) has been shown to depend on the mode of contraction: isometric or isotonic. The isometric ESFLR is derived from isometric contractions spanning a range of muscle lengths while the isotonic ESFLR is derived from shortening contractions across a range of afterloads. The ESFLR of isotonic contractions consistently lies below its isometric counterpart. Despite the passing of over a hundred years since the first insight by Otto Frank, the mechanism(s) underlying this protocol-dependent difference in the ESFLR remain incompletely explained. Here, we investigate the role of mechano-calcium feedback in accounting for the difference between these two ESFLRs. Previous studies have compared the dynamics of isotonic contractions to those of a single isometric contraction at a length that produces maximum force, without considering isometric contractions at shorter muscle lengths. We used a mathematical model of cardiac excitation-contraction to simulate isometric and force-length work-loop contractions (the latter being the 1D equivalent of the whole-heart pressure-volume loop), and compared Ca2+ transients produced under equivalent force conditions. We found that the duration of the simulated Ca2+ transient increases with decreasing sarcomere length for isometric contractions, and increases with decreasing afterload for work-loop contractions. At any given force, the Ca2+ transient for an isometric contraction is wider than that during a work-loop contraction. By driving simulated work-loops with wider Ca2+ transients generated from isometric contractions, we show that the duration of muscle shortening was prolonged, thereby shifting the work-loop ESFLR toward the isometric ESFLR. These observations are explained by an increase in the rate of binding of Ca2+ to troponin-C with increasing force. However, the leftward shift of the work-loop ESFLR does not superimpose on the isometric ESFLR, leading us to conclude that while mechano-calcium feedback does indeed contribute to the difference between the two ESFLRs, it does not completely account for it.

Efficient estimation of load-free left ventricular geometry and passive myocardial properties using principal component analysis.

Models of cardiac mechanics require a well-defined reference geometry from which deformations and hence myocardial strain and stress can be calculated. In the in vivo beating heart, the load-free (LF) geometry generally cannot be measured directly, since, in many cases, there is no stage at which the lumen pressures and contractile state are all zero. Therefore, there is a need for an efficient method to estimate the LF geometry, which is essential for an accurate mechanical simulation of left ventricular (LV) mechanics, and for estimations of passive and contractile constitutive parameters of the heart muscle. In this paper, we present a novel method for estimating both the LF geometry and the passive stiffness of the myocardium. A linear combination of principal components from a population of diastolic displacements is used to construct the LF geometry. For each estimate of the LF geometry and tissue stiffness, LV inflation is simulated, and the model predictions are compared with surface data at multiple stages during passive diastolic filling. The feasibility of this method was demonstrated using synthetically deformation data that were generated using LV models derived from clinical magnetic resonance image data, and the identifiability of the LF geometry and passive stiffness parameters were analysed using Hessian metrics. Applications of this method to clinical data would improve the accuracy of constitutive parameter estimation and allow a better simulation of LV wall strains and stresses.