Three-dimensional cardiovascular imaging-genetics: a mass univariate framework.

Motivation: Left ventricular (LV) hypertrophy is a strong predictor of cardiovascular outcomes, but its genetic regulation remains largely unexplained. Conventional phenotyping relies on manual calculation of LV mass and wall thickness, but advanced cardiac image analysis presents an opportunity for high-throughput mapping of genotype-phenotype associations in three dimensions (3D). Results: High-resolution cardiac magnetic resonance images were automatically segmented in 1124 healthy volunteers to create a 3D shape model of the heart. Mass univariate regression was used to plot a 3D effect-size map for the association between wall thickness and a set of predictors at each vertex in the mesh. The vertices where a significant effect exists were determined by applying threshold-free cluster enhancement to boost areas of signal with spatial contiguity. Experiments on simulated phenotypic signals and SNP replication show that this approach offers a substantial gain in statistical power for cardiac genotype-phenotype associations while providing good control of the false discovery rate. This framework models the effects of genetic variation throughout the heart and can be automatically applied to large population cohorts. Availability and implementation: The proposed approach has been coded in an R package freely available at https://doi.org/10.5281/zenodo.834610 together with the clinical data used in this work. Contact: declan.oregan@imperial.ac.uk. Supplementary information: Supplementary data are available at Bioinformatics online.

Machine Learning of Three-dimensional Right Ventricular Motion Enables Outcome Prediction in Pulmonary Hypertension: A Cardiac MR Imaging Study.

Purpose To determine if patient survival and mechanisms of right ventricular failure in pulmonary hypertension could be predicted by using supervised machine learning of three-dimensional patterns of systolic cardiac motion. Materials and Methods The study was approved by a research ethics committee, and participants gave written informed consent. Two hundred fifty-six patients (143 women; mean age ± standard deviation, 63 years ± 17) with newly diagnosed pulmonary hypertension underwent cardiac magnetic resonance (MR) imaging, right-sided heart catheterization, and 6-minute walk testing with a median follow-up of 4.0 years. Semiautomated segmentation of short-axis cine images was used to create a three-dimensional model of right ventricular motion. Supervised principal components analysis was used to identify patterns of systolic motion that were most strongly predictive of survival. Survival prediction was assessed by using difference in median survival time and area under the curve with time-dependent receiver operating characteristic analysis for 1-year survival. Results At the end of follow-up, 36% of patients (93 of 256) died, and one underwent lung transplantation. Poor outcome was predicted by a loss of effective contraction in the septum and free wall, coupled with reduced basal longitudinal motion. When added to conventional imaging and hemodynamic, functional, and clinical markers, three-dimensional cardiac motion improved survival prediction (area under the receiver operating characteristic curve, 0.73 vs 0.60, respectively; P < .001) and provided greater differentiation according to difference in median survival time between high- and low-risk groups (13.8 vs 10.7 years, respectively; P < .001). Conclusion A machine-learning survival model that uses three-dimensional cardiac motion predicts outcome independent of conventional risk factors in patients with newly diagnosed pulmonary hypertension. Online supplemental material is available for this article.

Relationship between body composition and left ventricular geometry using three dimensional cardiovascular magnetic resonance.

BACKGROUND: Although obesity is associated with alterations in left ventricular (LV) mass and volume which are of prognostic significance, widely differing patterns of remodelling have been attributed to adiposity. Our aim was to define the relationship between body composition and LV geometry using three-dimensional cardiovascular magnetic resonance. METHODS: In an observational study 1530 volunteers (55 % female, mean age 41.3 years) without known cardiovascular disease underwent investigation including breath-hold high spatial resolution 3D cines. Atlas-based segmentation and co-registration was used to create a statistical model of wall thickness (WT) and relative wall thickness (RWT) throughout the LV. The relationship between bio-impedence body composition and LV geometry was assessed using 3D regression models adjusted for age, systolic blood pressure (BP), gender, race and height, with correction to control the false discovery rate. RESULTS: LV mass was positively associated with fat mass in women but not in men (LV mass: women ß = 0.11, p < 0.0001; men ß = -0.01, p = 0.82). The 3D models revealed that in males fat mass was strongly associated with a concentric increase in relative wall thickness (RWT) throughout most of the LV (ß = 0.37, significant area = 96 %) and a reduced mid-ventricular cavity (ß = -0.22, significant area = 91 %). In women the regional concentric hypertrophic association was weaker, and the basal lateral wall showed an inverse relationship between RWT and fat mass (ß = -0.11, significant area = 4.8 %). CONCLUSIONS: In an adult population without known cardiovascular disease increasing body fat is predominately associated with asymmetric concentric hypertrophy independent of systolic BP, with women demonstrating greater cavity dilatation than men. Conventional mass and volume measurements underestimate the impact of body composition on LV structure due to anatomic variation in remodelling.

Population-based studies of myocardial hypertrophy: high resolution cardiovascular magnetic resonance atlases improve statistical power.

BACKGROUND: Cardiac phenotypes, such as left ventricular (LV) mass, demonstrate high heritability although most genes associated with these complex traits remain unidentified. Genome-wide association studies (GWAS) have relied on conventional 2D cardiovascular magnetic resonance (CMR) as the gold-standard for phenotyping. However this technique is insensitive to the regional variations in wall thickness which are often associated with left ventricular hypertrophy and require large cohorts to reach significance. Here we test whether automated cardiac phenotyping using high spatial resolution CMR atlases can achieve improved precision for mapping wall thickness in healthy populations and whether smaller sample sizes are required compared to conventional methods. METHODS: LV short-axis cine images were acquired in 138 healthy volunteers using standard 2D imaging and 3D high spatial resolution CMR. A multi-atlas technique was used to segment and co-register each image. The agreement between methods for end-diastolic volume and mass was made using Bland-Altman analysis in 20 subjects. The 3D and 2D segmentations of the LV were compared to manual labeling by the proportion of concordant voxels (Dice coefficient) and the distances separating corresponding points. Parametric and nonparametric data were analysed with paired t-tests and Wilcoxon signed-rank test respectively. Voxelwise power calculations used the interstudy variances of wall thickness. RESULTS: The 3D volumetric measurements showed no bias compared to 2D imaging. The segmented 3D images were more accurate than 2D images for defining the epicardium (Dice: 0.95 vs 0.93, P<0.001; mean error 1.3 mm vs 2.2 mm, P<0.001) and endocardium (Dice 0.95 vs 0.93, P<0.001; mean error 1.1 mm vs 2.0 mm, P<0.001). The 3D technique resulted in significant differences in wall thickness assessment at the base, septum and apex of the LV compared to 2D (P<0.001). Fewer subjects were required for 3D imaging to detect a 1 mm difference in wall thickness (72 vs 56, P<0.001). CONCLUSIONS: High spatial resolution CMR with automated phenotyping provides greater power for mapping wall thickness than conventional 2D imaging and enables a reduction in the sample size required for studies of environmental and genetic determinants of LV wall thickness.

A prospective evaluation of cardiovascular magnetic resonance measures of dyssynchrony in the prediction of response to cardiac resynchronization therapy.

BACKGROUND: Many patients with electrical dyssynchrony who undergo cardiac resynchronization therapy (CRT) do not obtain substantial benefit. Assessing mechanical dyssynchrony may improve patient selection. Results from studies using echocardiographic imaging to measure dyssynchrony have ultimately proved disappointing. We sought to evaluate cardiac motion in patients with heart failure and electrical dyssynchrony using cardiovascular magnetic resonance (CMR). We developed a framework for comparing measures of myocardial mechanics and evaluated how well they predicted response to CRT. METHODS: CMR was performed at 1.5 Tesla prior to CRT. Steady-state free precession (SSFP) cine images and complementary modulation of magnetization (CSPAMM) tagged cine images were acquired. Images were processed using a novel framework to extract regional ventricular volume-change, thickening and deformation fields (strain). A systolic dyssynchrony index (SDI) for all parameters within a 16-segment model of the ventricle was computed with high SDI denoting more dyssynchrony. Once identified, the optimal measure was applied to a second patient population to determine its utility as a predictor of CRT response compared to current accepted predictors (QRS duration, LBBB morphology and scar burden). RESULTS: Forty-four patients were recruited in the first phase (91% male, 63.3 ± 14.1 years; 80% NYHA class III) with mean QRSd 154 ± 24 ms. Twenty-one out of 44 (48%) patients showed reverse remodelling (RR) with a decrease in end systolic volume (ESV) ≥ 15% at 6 months. Volume-change SDI was the strongest predictor of RR (PR 5.67; 95% CI 1.95-16.5; P = 0.003). SDI derived from myocardial strain was least predictive. Volume-change SDI was applied as a predictor of RR to a second population of 50 patients (70% male, mean age 68.6 ± 12.2 years, 76% NYHA class III) with mean QRSd 146 ± 21 ms. When compared to QRSd, LBBB morphology and scar burden, volume-change SDI was the only statistically significant predictor of RR in this group. CONCLUSION: A systolic dyssynchrony index derived from volume-change is a highly reproducible measurement that can be derived from routinely acquired SSFP cine images and predicts RR following CRT whilst an SDI of regional strain does not.

Evaluation of current algorithms for segmentation of scar tissue from late gadolinium enhancement cardiovascular magnetic resonance of the left atrium: an open-access grand challenge.

BACKGROUND: Late Gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) imaging can be used to visualise regions of fibrosis and scarring in the left atrium (LA) myocardium. This can be important for treatment stratification of patients with atrial fibrillation (AF) and for assessment of treatment after radio frequency catheter ablation (RFCA). In this paper we present a standardised evaluation benchmarking framework for algorithms segmenting fibrosis and scar from LGE CMR images. The algorithms reported are the response to an open challenge that was put to the medical imaging community through an ISBI (IEEE International Symposium on Biomedical Imaging) workshop. METHODS: The image database consisted of 60 multicenter, multivendor LGE CMR image datasets from patients with AF, with 30 images taken before and 30 after RFCA for the treatment of AF. A reference standard for scar and fibrosis was established by merging manual segmentations from three observers. Furthermore, scar was also quantified using 2, 3 and 4 standard deviations (SD) and full-width-at-half-maximum (FWHM) methods. Seven institutions responded to the challenge: Imperial College (IC), Mevis Fraunhofer (MV), Sunnybrook Health Sciences (SY), Harvard/Boston University (HB), Yale School of Medicine (YL), King's College London (KCL) and Utah CARMA (UTA, UTB). There were 8 different algorithms evaluated in this study. RESULTS: Some algorithms were able to perform significantly better than SD and FWHM methods in both pre- and post-ablation imaging. Segmentation in pre-ablation images was challenging and good correlation with the reference standard was found in post-ablation images. Overlap scores (out of 100) with the reference standard were as follows: Pre: IC = 37, MV = 22, SY = 17, YL = 48, KCL = 30, UTA = 42, UTB = 45; Post: IC = 76, MV = 85, SY = 73, HB = 76, YL = 84, KCL = 78, UTA = 78, UTB = 72. CONCLUSIONS: The study concludes that currently no algorithm is deemed clearly better than others. There is scope for further algorithmic developments in LA fibrosis and scar quantification from LGE CMR images. Benchmarking of future scar segmentation algorithms is thus important. The proposed benchmarking framework is made available as open-source and new participants can evaluate their algorithms via a web-based interface.

3WJ RNA Nanoparticles-Aptamer Functionalized Exosomes From M2 Macrophages Target BMSCs to Promote the Healing of Bone Fractures.

Up to now, impaired bone regeneration severely affects the healing of bone fractures, thus bringing tremendous suffering to patients. As a vital mediator between inflammatory response and bone regeneration, M2 macrophage-derived exosomes (M2-Exos) attenuate inflammation and promote tissue repair. However, due to a lack of specific targeting property, M2-Exos will be rapidly eliminated after systematic administration, thus compromising their effectiveness in promoting bone regeneration. To solve this hurdle, we initially harvested and characterized the pro-osteogenic properties of M2-Exos. A bone marrow mesenchymal stem cell (BMSC)-specific aptamer was synthesized and 3-way junction (3WJ) RNA nanoparticles were applied to conjugate the BMSC-specific aptamer and M2-Exos. In vitro assays revealed that M2-Exos bore the representative features of exosomes and significantly promoted the proliferation, migration, and osteogenic differentiation of BMSCs. 3WJ RNA nanoparticles-aptamer functionalized M2-Exos (3WJ-BMSCapt/M2-Exos) maintained the original physical characteristics of M2-Exos, but bore a high specific binding ability to BMSCs. Furthermore, when being systemically administered in the mice model with femoral bone fractures, these functionalized M2-Exos mainly accumulated at the bone fracture site with a slow release of exosomal cargo, thereby significantly accelerating the healing processes compared with the M2-Exos group. Our study indicated that the 3WJ-BMSCapt/M2-Exos with BMSCs targeting ability and controlled release would be a promising strategy to treat bone fractures.

Remodeling after acute myocardial infarction: mapping ventricular dilatation using three dimensional CMR image registration.

BACKGROUND: Progressive heart failure due to remodeling is a major cause of morbidity and mortality following myocardial infarction. Conventional clinical imaging measures global volume changes, and currently there is no means of assessing regional myocardial dilatation in relation to ischemic burden. Here we use 3D co-registration of Cardiovascular Magnetic Resonance (CMR) images to assess the long-term effects of ischemia-reperfusion injury on left ventricular structure after acute ST-elevation myocardial infarction (STEMI). METHODS: Forty six patients (age range 33-77 years) underwent CMR imaging within 7 days following primary percutaneous coronary intervention (PPCI) for acute STEMI with follow-up at one year. Functional cine imaging and Late Gadolinium Enhancement (LGE) were segmented and co-registered. Local left ventricular wall dilatation was assessed by using intensity-based similarities to track the structural changes in the heart between baseline and follow-up. Results are expressed as means, standard errors and 95% confidence interval (CI) of the difference. RESULTS: Local left ventricular remodeling within infarcted myocardium was greater than in non-infarcted myocardium (1.6%±1.0 vs 0.3%±0.9, 95% CI: -2.4% - -0.2%, P=0.02). One-way ANOVA revealed that transmural infarct thickness had a significant effect on the degree of local remodeling at one year (P<0.0001) with greatest wall dilatation observed when infarct transmurality exceeded 50%. Infarct remodeling was more severe when microvascular obstruction (MVO) was present (3.8%±1.3 vs -1.6%±1.4, 95% CI: -9.1% - -1.5%, P=0.007) and when end-diastolic volume had increased by >20% (4.8%±1.4 vs -0.15%±1.2, 95% CI: -8.9% - -0.9%, P=0.017). CONCLUSIONS: The severity of ischemic injury has a significant effect on local ventricular wall remodeling with only modest dilatation observed within non-ischemic myocardium. Limitation of chronic remodeling may therefore depend on therapies directed at modulating ischemia-reperfusion injury. CMR co-registration has potential for assessing dynamic changes in ventricular structure in relation to therapeutic interventions.

Metabolic pathways associated with right ventricular adaptation to pulmonary hypertension: 3D analysis of cardiac magnetic resonance imaging.

AIMS: We sought to identify metabolic pathways associated with right ventricular (RV) adaptation to pulmonary hypertension (PH). We evaluated candidate metabolites, previously associated with survival in pulmonary arterial hypertension, and used automated image segmentation and parametric mapping to model their relationship to adverse patterns of remodelling and wall stress. METHODS AND RESULTS: In 312 PH subjects (47.1% female, mean age 60.8 ± 15.9 years), of which 182 (50.5% female, mean age 58.6 ± 16.8 years) had metabolomics, we modelled the relationship between the RV phenotype, haemodynamic state, and metabolite levels. Atlas-based segmentation and co-registration of cardiac magnetic resonance imaging was used to create a quantitative 3D model of RV geometry and function-including maps of regional wall stress. Increasing mean pulmonary artery pressure was associated with hypertrophy of the basal free wall (ß = 0.29) and reduced relative wall thickness (ß = -0.38), indicative of eccentric remodelling. Wall stress was an independent predictor of all-cause mortality (hazard ratio = 1.27, P = 0.04). Six metabolites were significantly associated with elevated wall stress (ß = 0.28-0.34) including increased levels of tRNA-specific modified nucleosides and fatty acid acylcarnitines, and decreased levels (ß = -0.40) of sulfated androgen. CONCLUSION: Using computational image phenotyping, we identify metabolic profiles, reporting on energy metabolism and cellular stress-response, which are associated with adaptive RV mechanisms to PH.

Multi-atlas segmentation with augmented features for cardiac MR images.

Multi-atlas segmentation infers the target image segmentation by combining prior anatomical knowledge encoded in multiple atlases. It has been quite successfully applied to medical image segmentation in the recent years, resulting in highly accurate and robust segmentation for many anatomical structures. However, to guide the label fusion process, most existing multi-atlas segmentation methods only utilise the intensity information within a small patch during the label fusion process and may neglect other useful information such as gradient and contextual information (the appearance of surrounding regions). This paper proposes to combine the intensity, gradient and contextual information into an augmented feature vector and incorporate it into multi-atlas segmentation. Also, it explores the alternative to the K nearest neighbour (KNN) classifier in performing multi-atlas label fusion, by using the support vector machine (SVM) for label fusion instead. Experimental results on a short-axis cardiac MR data set of 83 subjects have demonstrated that the accuracy of multi-atlas segmentation can be significantly improved by using the augmented feature vector. The mean Dice metric of the proposed segmentation framework is 0.81 for the left ventricular myocardium on this data set, compared to 0.79 given by the conventional multi-atlas patch-based segmentation (Coupé et al., 2011; Rousseau et al., 2011). A major contribution of this paper is that it demonstrates that the performance of non-local patch-based segmentation can be improved by using augmented features.

Multiatlas whole heart segmentation of CT data using conditional entropy for atlas ranking and selection.

PURPOSE: Cardiac computed tomography (CT) is widely used in clinical diagnosis of cardiovascular diseases. Whole heart segmentation (WHS) plays a vital role in developing new clinical applications of cardiac CT. However, the shape and appearance of the heart can vary greatly across different scans, making the automatic segmentation particularly challenging. The objective of this work is to develop and evaluate a multiatlas segmentation (MAS) scheme using a new atlas ranking and selection algorithm for automatic WHS of CT data. Research on different MAS strategies and their influence on WHS performance are limited. This work provides a detailed comparison study evaluating the impacts of label fusion, atlas ranking, and sizes of the atlas database on the segmentation performance. METHODS: Atlases in a database were registered to the target image using a hierarchical registration scheme specifically designed for cardiac images. A subset of the atlases were selected for label fusion, according to the authors' proposed atlas ranking criterion which evaluated the performance of each atlas by computing the conditional entropy of the target image given the propagated atlas labeling. Joint label fusion was used to combine multiple label estimates to obtain the final segmentation. The authors used 30 clinical cardiac CT angiography (CTA) images to evaluate the proposed MAS scheme and to investigate different segmentation strategies. RESULTS: The mean WHS Dice score of the proposed MAS method was 0.918 ± 0.021, and the mean runtime for one case was 13.2 min on a workstation. This MAS scheme using joint label fusion generated significantly better Dice scores than the other label fusion strategies, including majority voting (0.901 ± 0.276, p < 0.01), locally weighted voting (0.905 ± 0.0247, p < 0.01), and probabilistic patch-based fusion (0.909 ± 0.0249, p < 0.01). In the atlas ranking study, the proposed criterion based on conditional entropy yielded a performance curve with higher WHS Dice scores compared to the conventional schemes (p < 0.03). In the atlas database study, the authors showed that the MAS using larger atlas databases generated better performance curves than the MAS using smaller ones, indicating larger atlas databases could produce more accurate segmentation. CONCLUSIONS: The authors have developed a new MAS framework for automatic WHS of CTA and investigated alternative implementations of MAS. With the proposed atlas ranking algorithm and joint label fusion, the MAS scheme is able to generate accurate segmentation within practically acceptable computation time. This method can be useful for the development of new clinical applications of cardiac CT.

A bi-ventricular cardiac atlas built from 1000+ high resolution MR images of healthy subjects and an analysis of shape and motion.

Atlases encode valuable anatomical and functional information from a population. In this work, a bi-ventricular cardiac atlas was built from a unique data set, which consists of high resolution cardiac MR images of 1000+ normal subjects. Based on the atlas, statistical methods were used to study the variation of cardiac shapes and the distribution of cardiac motion across the spatio-temporal domain. We have shown how statistical parametric mapping (SPM) can be combined with a general linear model to study the impact of gender and age on regional myocardial wall thickness. Finally, we have also investigated the influence of the population size on atlas construction and atlas-based analysis. The high resolution atlas, the statistical models and the SPM method will benefit more studies on cardiac anatomy and function analysis in the future.

4D Blood Flow Reconstruction Over the Entire Ventricle From Wall Motion and Blood Velocity Derived From Ultrasound Data.

We demonstrate a new method to recover 4D blood flow over the entire ventricle from partial blood velocity measurements using multiple 3D+t colour Doppler images and ventricular wall motion estimated using 3D+t BMode images. We apply our approach to realistic simulated data to ascertain the ability of the method to deal with incomplete data, as typically happens in clinical practice. Experiments using synthetic data show that the use of wall motion improves velocity reconstruction, shows more accurate flow patterns and improves mean accuracy particularly when coverage of the ventricle is poor. The method was applied to patient data from 6 congenital cases, producing results consistent with the simulations. The use of wall motion produced more plausible flow patterns and reduced the reconstruction error in all patients.

Precursors of Hypertensive Heart Phenotype Develop in Healthy Adults: A High-Resolution 3D MRI Study.

OBJECTIVES: This study used high-resolution 3-dimensional cardiac magnetic resonance to define the anatomical and functional left ventricular (LV) properties associated with increasing systolic blood pressure (SBP) in a drug-naïve cohort. BACKGROUND: LV hypertrophy and remodeling occur in response to hemodynamic stress but little is known about how these phenotypic changes are initiated in the general population. METHODS: In this study, 1,258 volunteers (54% women, mean age 40.6 ± 12.8 years) without self-reported cardiovascular disease underwent 3-dimensional cardiac magnetic resonance combined with computational modeling. The relationship between SBP and wall thickness (WT), relative WT, end-systolic wall stress (WS), and fractional wall thickening were analyzed using 3-dimensional regression models adjusted for body surface area, sex, race, age, and multiple testing. Significantly associated points in the LV model (p < 0.05) were identified and the relationship with SBP reported as mean ß coefficients. RESULTS: There was a continuous relationship between SBP and asymmetric concentric hypertrophic adaptation of the septum and anterior wall that was associated with normalization of wall stress. In the lateral wall an increase in wall stress with rising SBP was not balanced by a commensurate hypertrophic relationship. In normotensives, SBP was positively associated with WT (ß = 0.09) and relative WT (ß = 0.07) in the septal and anterior walls, and this regional hypertrophic relationship was progressively stronger among pre-hypertensives (ß = 0.10) and hypertensives (ß = 0.30). CONCLUSIONS: These findings show that the precursors of the hypertensive heart phenotype can be traced to healthy normotensive adults and that an independent and continuous relationship exists between adverse LV remodeling and SBP in a low-risk population. These adaptations show distinct regional variations with concentric hypertrophy of the septum and eccentric hypertrophy of the lateral wall, which challenge conventional classifications of LV remodeling.

Right ventricle segmentation from cardiac MRI: a collation study.

Magnetic Resonance Imaging (MRI), a reference examination for cardiac morphology and function in humans, allows to image the cardiac right ventricle (RV) with high spatial resolution. The segmentation of the RV is a difficult task due to the variable shape of the RV and its ill-defined borders in these images. The aim of this paper is to evaluate several RV segmentation algorithms on common data. More precisely, we report here the results of the Right Ventricle Segmentation Challenge (RVSC), concretized during the MICCAI'12 Conference with an on-site competition. Seven automated and semi-automated methods have been considered, along them three atlas-based methods, two prior based methods, and two prior-free, image-driven methods that make use of cardiac motion. The obtained contours were compared against a manual tracing by an expert cardiac radiologist, taken as a reference, using Dice metric and Hausdorff distance. We herein describe the cardiac data composed of 48 patients, the evaluation protocol and the results. Best results show that an average 80% Dice accuracy and a 1cm Hausdorff distance can be expected from semi-automated algorithms for this challenging task on the datasets, and that an automated algorithm can reach similar performance, at the expense of a high computational burden. Data are now publicly available and the website remains open for new submissions (http://www.litislab.eu/rvsc/).

Understanding the need of ventricular pressure for the estimation of diastolic biomarkers.

The diastolic function (i.e., blood filling) of the left ventricle (LV) is determined by its capacity for relaxation, or the decay in residual active tension (AT) generated during systole, and its constitutive material properties, or myocardial stiffness. The clinical determination of these two factors (diastolic residual AT and stiffness) is thus essential for assessing LV diastolic function. To quantify these two factors, in our previous work, a novel model-based parameter estimation approach was proposed and successfully applied to multiple cases using clinically acquired motion and invasively measured ventricular pressure data. However, the need to invasively acquire LV pressure limits the wide application of this approach. In this study, we address this issue by analyzing the feasibility of using two kinds of non-invasively available pressure measurements for the purpose of inverse mechanical parameter estimation. The prescription of pressure based on a generic pressure-volume (P-V) relationship reported in literature is first evaluated in a set of 18 clinical cases (10 healthy and 8 diseased), finding reasonable results for stiffness but not for residual active tension. We then investigate the use of non-invasive pressure measures, now available through imaging techniques and limited by unknown or biased offset values. Specifically, three sets of physiologically realistic synthetic data with three levels of diastolic residual active tension (i.e., impaired relaxation capability) are designed to quantify the percentage error in the parameter estimation against the possible pressure offsets within the physiological limits. Maximum errors are quantified as 11 % for the magnitude of stiffness and 22 % for AT, with averaged 0.17 kPa error in pressure measurement offset using the state-of-the-art non-invasive pressure estimation method. The main cause for these errors is the limited temporal resolution of clinical imaging data currently available. These results demonstrate the potential feasibility of the estimation diastolic biomarkers with non-invasive assessment of pressure through medical imaging data.

Multi-atlas spectral PatchMatch: application to cardiac image segmentation.

The automatic segmentation of cardiac magnetic resonance images poses many challenges arising from the large variation between different anatomies, scanners and acquisition protocols. In this paper, we address these challenges with a global graph search method and a novel spectral embedding of the images. Firstly, we propose the use of an approximate graph search approach to initialize patch correspondences between the image to be segmented and a database of labelled atlases, Then, we propose an innovative spectral embedding using a multi-layered graph of the images in order to capture global shape properties. Finally, we estimate the patch correspondences based on a joint spectral representation of the image and atlases. We evaluated the proposed approach using 155 images from the recent MICCAI SATA segmentation challenge and demonstrated that the proposed algorithm significantly outperforms current state-of-the-art methods on both training and test sets.

Model-guided directional minimal path for fully automatic extraction of coronary centerlines from cardiac CTA.

Extracting centerlines of coronary arteries is a challenging but important task in clinical applications of cardiac CTA. In this paper, we propose a model-guided approach, the directional minimal path, for the centerline extraction. The proposed method is based on the minimal path algorithm and a prior coronary model is used. The model is first registered to the unseen image. Then, the start point and end point for the minimal path algorithm are provided by the model to automate the centerline extraction process. Also, the direction information of the coronary model is used to guide the path tracking of the minimal path procedure. This directional tracking improves the robustness and accuracy of the centerline extraction. Finally, the proposed method can automatically recognize the branches of the extracted coronary artery using the prior information in the model. We validated the proposed method by extracting the three main coronary branches. The mean accuracy of the 56 cases was 1.32+/-0.81 mm and the detection ratio was 88.7%.

A probabilistic patch-based label fusion model for multi-atlas segmentation with registration refinement: application to cardiac MR images.

The evaluation of ventricular function is important for the diagnosis of cardiovascular diseases. It typically involves measurement of the left ventricular (LV) mass and LV cavity volume. Manual delineation of the myocardial contours is time-consuming and dependent on the subjective experience of the expert observer. In this paper, a multi-atlas method is proposed for cardiac magnetic resonance (MR) image segmentation. The proposed method is novel in two aspects. First, it formulates a patch-based label fusion model in a Bayesian framework. Second, it improves image registration accuracy by utilizing label information, which leads to improvement of segmentation accuracy. The proposed method was evaluated on a cardiac MR image set of 28 subjects. The average Dice overlap metric of our segmentation is 0.92 for the LV cavity, 0.89 for the right ventricular cavity and 0.82 for the myocardium. The results show that the proposed method is able to provide accurate information for clinical diagnosis.

Temporal sparse free-form deformations.

FFD represent a widely used model for the non-rigid registration of medical images. The balance between robustness to noise and accuracy in modelling localised motion is typically controlled by the control point grid spacing and the amount of regularisation. More recently, TFFD have been proposed which extend the FFD approach in order to recover smooth motion from temporal image sequences. In this paper, we revisit the classic FFD approach and propose a sparse representation using the principles of compressed sensing. The sparse representation can model both global and local motion accurately and robustly. We view the registration as a deformation reconstruction problem. The deformation is reconstructed from a pair of images (or image sequences) with a sparsity constraint applied to the parametric space. Specifically, we introduce sparsity into the deformation via L1 regularisation, and apply a bending energy regularisation between neighbouring control points within each level to encourage a grouped sparse solution. We further extend the sparsity constraint to the temporal domain and propose a TSFFD which can capture fine local details such as motion discontinuities in both space and time without sacrificing robustness. We demonstrate the capabilities of the proposed framework to accurately estimate deformations in dynamic 2D and 3D image sequences. Compared to the classic FFD and TFFD approach, a significant increase in registration accuracy can be observed in natural images as well as in cardiac images.