Search

Background and Purpose: A recent paradigm shift in proarrhythmic risk assessment suggests that the integration of clinical, non-clinical, and computational evidence can be used to reach a comprehensive understanding of the proarrhythmic potential of drug candidates. While current computational methodologies focus on predicting the incidence of proarrhythmic events after drug administration, the objective of this study is to predict concentration-response relationships of QTc as a clinical endpoint., Experimental Approach: Full heart computational models reproducing human cardiac populations were created to predict the concentration-response relationship of changes in the QT interval as recommended for clinical trials. The concentration-response relationship of the QT-interval prolongation obtained from the computational cardiac population was compared against the relationship from clinical trial data for a set of well-characterized compounds: moxifloxacin, dofetilide, verapamil, and ondansetron., Key Results: Computationally derived concentration-response relationships of QT interval changes for three of the four drugs had slopes within the confidence interval of clinical trials (dofetilide, moxifloxacin and verapamil) when compared to placebo-corrected concentration-ΔQT and concentration-ΔQT regressions. Moxifloxacin showed a higher intercept, outside the confidence interval of the clinical data, demonstrating that in this example, the standard linear regression does not appropriately capture the concentration-response results at very low concentrations. The concentrations corresponding to a mean QTc prolongation of 10 ms were consistently lower in the computational model than in clinical data. The critical concentration varied within an approximate ratio of 0.5 (moxifloxacin and ondansetron) and 1 times (dofetilide, verapamil) the critical concentration observed in human clinical trials. Notably, no other in silico methodology can approximate the human critical concentration values for a QT interval prolongation of 10 ms., Conclusion and Implications: Computational concentration-response modelling of a virtual population of high-resolution, 3-dimensional cardiac models can provide comparable information to clinical data and could be used to complement pre-clinical and clinical safety packages. It provides access to an unlimited exposure range to support trial design and can improve the understanding of pre-clinical-clinical translation., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Elem Biotech owns the commercial rights to Alya, the computational finite element solver employed in this study. However, any other commercial or academic finite element solver could be employed to reproduce this work. Elem Biotech owns the commercial rights to Alya, the solver employed in this study., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Following the 3 R's principles of animal research-replacement, reduction, and refinement-a high-performance computational framework was produced to generate a platform to perform human cardiac in-silico clinical trials as means to assess the pro-arrhythmic risk after the administrations of one or combination of two potentially arrhythmic drugs. The drugs assessed in this study were hydroxychloroquine and azithromycin. The framework employs electrophysiology simulations on high-resolution three-dimensional, biventricular human heart anatomies including phenotypic variabilities, so as to determine if differential QT-prolongation responds to drugs as observed clinically. These simulations also reproduce sex-specific ionic channel characteristics. The derived changes in the pseudo-electrocardiograms, calcium concentrations, as well as activation patterns within 3D geometries were evaluated for signs of induced arrhythmia. The virtual subjects could be evaluated at two different cycle lengths: at a normal heart rate and at a heart rate associated with stress as means to analyze the proarrhythmic risks after the administrations of hydroxychloroquine and azithromycin. Additionally, a series of experiments performed on reanimated swine hearts utilizing Visible Heart ® methodologies in a four-chamber working heart model were performed to verify the arrhythmic behaviors observed in the in silico trials.The obtained results indicated similar pro-arrhythmic risk assessments within the virtual population as compared to published clinical trials (21% clinical risk vs 21.8% in silico trial risk). Evidence of transmurally heterogeneous action potential prolongations after providing a large dose of hydroxychloroquine was found as the observed mechanisms for elicited arrhythmias, both in the in vitro and the in silico models. The proposed workflow for in silico clinical drug cardiotoxicity trials allows for reproducing the complex behavior of cardiac electrophysiology in a varied population, in a matter of a few days as compared to the months or years it requires for most in vivo human clinical trials. Importantly, our results provided evidence of the common phenotype variants that produce distinct drug-induced arrhythmogenic outcomes., (© 2024. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Stroke is the second leading cause of death worldwide. Nearly two-thirds of strokes are produced by cardioembolisms, and half of cardioembolic strokes are triggered by Atrial Fibrillation (AF), the most common type of arrhythmia. A more recent cause of cardioembolisms is Transcatheter Aortic Valve Replacements (TAVRs), which may onset post-procedural adverse events such as stroke and Silent Brain Infarcts (SBIs), for which no definitive treatment exists, and which will only get worse as TAVRs are implanted in younger and lower risk patients. It is well known that some specific characteristics of elderly patients may lower the safety and efficacy of anticoagulation therapy, making it a real urgency to find alternative therapies. We propose a device consisting of a strut structure placed at the base of the treated artery to model the potential risk of cerebral embolisms caused by dislodged debris of varying sizes. This work analyzes a design based on a patented medical device, intended to block cardioembolisms from entering the cerebrovascular system, with a particular focus on AF, and potentially TAVR patients. The study has been carried out in two stages. Both of them based on computational fluid dynamics (CFD) coupled with Lagrangian particle tracking method. The first stage of the work evaluates a variety of strut thicknesses and inter-strut spacings, contrasting with the device-free baseline geometry. The analysis is carried out by imposing flowrate waveforms characteristic of both healthy and AF patients. Boundary conditions are calibrated to reproduce physiological flowrates and pressures in a patient's aortic arch. In the second stage, the optimal geometric design from the first stage was employed, with the addition of lateral struts to prevent the filtration of particles and electronegatively charged strut surfaces, studying the effect of electrical forces on the clots if they are considered charged. Flowrate boundary conditions were used to emulate both healthy and AF conditions. Results from numerical simulations coming form the first stage indicate that the device blocks particles of sizes larger than the inter-strut spacing. It was found that lateral strut space had the highest impact on efficacy. Based on the results of the second stage, deploying the electronegatively charged device in all three aortic arch arteries, the number of particles entering these arteries was reduced on average by 62.6% and 51.2%, for the healthy and diseased models respectively, matching or surpassing current oral anticoagulant efficacy. In conclusion, the device demonstrated a two-fold mechanism for filtering emboli: while the smallest particles are deflected by electrostatic repulsion, avoiding microembolisms, which could lead to cognitive impairment, the largest ones are mechanically filtered since they cannot fit in between the struts, effectively blocking the full range of particle sizes analyzed in this study. The device presented in this manuscript offers an anticoagulant-free method to prevent stroke and SBIs, imperative given the growing population of AF and elderly patients., Competing Interests: DO, MF, ACM and MV are employed by ELEM Biotech. MV is an investor of ELEM Biotech. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2023 Eguzkitza, Oks, Navia, Houzeaux, Butakoff, Fisa, Campoy Millán and Vázquez.)

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O'Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studies., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Gonzalez-Martin et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

This work intends to study the effect of aortic annulus eccentricity and leaflet rigidity on the performance, thrombogenic risk and calcification risk in bioprosthetic aortic valve replacements (BAVRs). To address these questions, a two-way immersed fluid-structure interaction (FSI) computational model was implemented in a high-performance computing (HPC) multi-physics simulation software, and validated against a well-known FSI benchmark. The aortic valve bioprosthesis model is qualitatively contrasted against experimental data, showing good agreement in closed and open states. Regarding the performance of BAVRs, the model predicts that increasing eccentricities yield lower geometric orifice areas (GOAs) and higher normalized transvalvular pressure gradients (TPGs) for healthy cardiac outputs during systole, agreeing with in vitro experiments. Regions with peak values of residence time are observed to grow with eccentricity in the sinus of Valsalva, indicating an elevated risk of thrombus formation for eccentric configurations. In addition, the computational model is used to analyze the effect of varying leaflet rigidity on both performance, thrombogenic and calcification risks with applications to tissue-engineered prostheses. For more rigid leaflets it predicts an increase in systolic and diastolic TPGs, and decrease in systolic GOA, which translates to decreased valve performance. The peak shear rate and residence time regions increase with leaflet rigidity, but their volume-averaged values were not significantly affected. Peak solid stresses are also analyzed, and observed to increase with rigidity, elevating risk of valve calcification and structural failure. To the authors' knowledge this is the first computational FSI model to study the effect of eccentricity or leaflet rigidity on thrombogenic biomarkers, providing a novel tool to aid device manufacturers and clinical practitioners., (© 2022 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd.)

Background: Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions., Methods: The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC's in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow [Formula: see text], maximum aortic root flow [Formula: see text], average LVAD flow [Formula: see text], and maximum LVAD flow [Formula: see text]); (d) an uncertainty quantification using six validation experiments that include extreme operating conditions., Results: Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer., Conclusions: This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: AS, MV, BE and CB have acted as consultants for Medtronic PLC related to Medtronic’s HVAD. KMN and VV have an ongoing research study for Abbott, Inc for work on Abbott’s HeartMate III. RG and PP need to include the following statement: “The mention of commercial products, their sources, or their use in connection with material reported herein is not to be construed as either an actual or implied endorsement of such products by the U.S. Department of Health and Human Services”.

Venous-arterial extracorporeal membrane oxygenation (VA-ECMO) treatment for acute cardiogenic shock in patients who also have acute lung injury predisposes development of a serious complication called "north-south syndrome" (NSS) which causes cerebral hypoxia. NSS is poorly characterized and hemodynamic studies have focused on cerebral perfusion ignoring the heart. We hypothesized in NSS the heart would be more likely to receive hypoxemic blood than the brain due to the proximity of the coronary arteries to the aortic annulus. To test this, we conducted a computational fluid dynamics simulation of blood flow in a human supported by VA-ECMO. Simulations quantified the fraction of blood at each aortic branching vessel originating from residual native cardiac output versus VA-ECMO. As residual cardiac function was increased, simulations demonstrated myocardial hypoxia would develop prior to cerebral hypoxia. These results illustrate the conditions where NSS will develop and the relative cardiac function that will lead to organ-specific hypoxia. Illustration of the impact of north-south syndrome on organ-specific oxygen delivery. Patients on VA-ECMO have two sources of blood flow, one from the VA-ECMO circuit and one from the residual cardiac function. When there is no residual cardiac function, all organs are perfused with oxygenated blood. As myocardial recovery progresses, blood supply from the two sources will begin to mix resulting in non-homogeneous mixing and differential oxygenation based upon the anatomical site of branching vessels., (© 2021. The Author(s).)

Importance: Being born small for gestational age (SGA), approximately 10% of all births, is associated with increased risk of cardiovascular mortality in adulthood, but mechanistic pathways are unclear. Cardiac remodeling and dysfunction occur in fetuses SGA and children born SGA, but it is uncertain whether and how these changes persist into adulthood., Objective: To evaluate baseline cardiac function and structure and exercise capacity in young adults born SGA., Design, Setting, and Participants: This cohort study conducted from January 2015 to January 2018 assessed a perinatal cohort born at a tertiary university hospital in Spain between 1975 and 1995. Participants included 158 randomly selected young adults aged 20 to 40 years born SGA (birth weight below the 10th centile) or with intrauterine growth within standard reference ranges (controls). Participants provided their medical history, filled out questionnaires regarding smoking and physical activity habits, and underwent incremental cardiopulmonary exercise stress testing, cardiac magnetic resonance imaging, and a physical examination, with blood pressure, glucose level, and lipid profile data collected., Exposure: Being born SGA., Main Outcomes and Measures: Cardiac structure and function assessed by cardiac magnetic resonance imaging, including biventricular end-diastolic shape analysis. Exercise capacity assessed by incremental exercise stress testing., Results: This cohort study included 81 adults born SGA (median age at study, 34.4 years [IQR, 30.8-36.7 years]; 43 women [53%]) and 77 control participants (median age at study, 33.7 years [interquartile range (IQR), 31.0-37.1 years]; 33 women [43%]). All participants were of White race/ethnicity and underwent imaging, whereas 127 participants (80% of the cohort; 66 control participants and 61 adults born SGA) completed the exercise test. Cardiac shape analysis showed minor changes at rest in right ventricular geometry (DeLong test z, 2.2098; P = .02) with preserved cardiac function in individuals born SGA. However, compared with controls, adults born SGA had lower exercise capacity, with decreased maximal workload (mean [SD], 180 [62] W vs 214 [60] W; P = .006) and oxygen consumption (median, 26.0 mL/min/kg [IQR, 21.5-33.5 mL/min/kg vs 29.5 mL/min/kg [IQR, 24.0-36.0 mL/min/kg]; P = .02). Exercise capacity was significantly correlated with left ventricular mass (ρ = 0.7934; P < .001)., Conclusions and Relevance: This cohort of young adults born SGA had markedly reduced exercise capacity. These results support further research to clarify the causes of impaired exercise capacity and the potential association with increased cardiovascular mortality among adults born SGA.

Statistical shape analysis is a powerful tool to assess organ morphologies and find shape changes associated to a particular disease. However, imbalance in confounding factors, such as demographics might invalidate the analysis if not taken into consideration. Despite the methodological advances in the field, providing new methods that are able to capture complex and regional shape differences, the relationship between non-imaging information and shape variability has been overlooked. We present a linear statistical shape analysis framework that finds shape differences unassociated to a controlled set of confounding variables. It includes two confounding correction methods: confounding deflation and adjustment. We applied our framework to a cardiac magnetic resonance imaging dataset, consisting of the cardiac ventricles of 89 triathletes and 77 controls, to identify cardiac remodelling due to the practice of endurance exercise. To test robustness to confounders, subsets of this dataset were generated by randomly removing controls with low body mass index, thus introducing imbalance. The analysis of the whole dataset indicates an increase of ventricular volumes and myocardial mass in athletes, which is consistent with the clinical literature. However, when confounders are not taken into consideration no increase of myocardial mass is found. Using the downsampled datasets, we find that confounder adjustment methods are needed to find the real remodelling patterns in imbalanced datasets., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests., (Copyright © 2020. Published by Elsevier B.V.)

Two-dimensional representation of 3D anatomical structures is a simple and intuitive way for analysing patient information across populations and image modalities. While cardiac ventricles, especially the left ventricle, have an established standard representation (bull's eye plot), the 2D depiction of the left atrium (LA) remains challenging due to its sub-structural complexity including the pulmonary veins (PV) and the left atrial appendage (LAA). Quasi-conformal flattening techniques, successfully applied to cardiac ventricles, require additional constraints in the case of the LA to place the PV and LAA in the same geometrical 2D location for different cases. Some registration-based methods have been proposed but surface registration is time-consuming and prone to errors when the geometries are very different. We propose a novel atrial flattening methodology where a 2D standardised map of the LA is obtained quickly and without errors related to registration. The LA is divided into five regions which are then mapped to their analogue two-dimensional regions. 67 human left atria from magnetic resonance images (MRI) were studied to derive a population-based template representing the averaged relative locations of the PVs and LAA. The clinical application of our methodology is illustrated on different use cases including the integration of MRI and electroanatomical data.

Aims: Late gadolinium enhancement cardiac magnetic resonance (LGE-CMR) may define left atrial (LA) anatomy and structural remodelling, and facilitate atrial fibrillation (AF) ablation. We aimed to assess the intra- and inter-observer reproducibility and agreement of LGE-CMR parameters with direct application to AF ablation techniques., Methods and Results: One experienced and one non-experienced observer performed complete LGE-CMR data analysis twice, on different days, in 40 randomly selected LGE-CMR examinations [20 performed before ablation (pre-ablation) and 20 performed 3 months after ablation (post-ablation)]. Four additional observers (two experienced and two non-experienced) performed complete LGE-CMR data analysis in a subgroup of 30 patients (15 pre-ablation and 15 post-ablation). All LGE-CMR were performed in sinus rhythm. Intra- and inter-observer reproducibility of LA volume, LA area, and sphericity index (SI) was high: coefficient of variation <10% and intraclass correlation coefficient >0.71. Geometric congruency of repeated reconstruction of LA shape was high: maximal error <5 mm for intra-observer and <8 mm for inter-observer. The precision of scar location increased with extent of scar, and was high (Dice coefficient >0.75) when the scar area was >5 cm2 for a single observer and >15 cm2 for multiple observers. Non-experienced observers performed equally well to experienced observers., Conclusion: Late gadolinium enhancement cardiac magnetic resonance measurements of LA area, volume, and SI were reproducible, and geometric congruency of LA shape was high. Location of scar was precise for scar areas >5 cm2 for single observers and >15 cm2 for multiple observers, regardless of the observers' experience. These results may serve as a reference for future studies on the role for substrate-based AF ablation procedures., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author(s) 2019. For permissions, please email: journals.permissions@oup.com.)

Pulmonary vein isolation (PVI) is a common procedure for the treatment of atrial fibrillation (AF) since the initial trigger for AF frequently originates in the pulmonary veins. A successful isolation produces a continuous lesion (scar) completely encircling the veins that stops activation waves from propagating to the atrial body. Unfortunately, the encircling lesion is often incomplete, becoming a combination of scar and gaps of healthy tissue. These gaps are potential causes of AF recurrence, which requires a redo of the isolation procedure. Late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) is a non-invasive method that may also be used to detect gaps, but it is currently a time-consuming process, prone to high inter-observer variability. In this paper, we present a method to semi-automatically identify and quantify ablation gaps. Gap quantification is performed through minimum path search in a graph where every node is a scar patch and the edges are the geodesic distances between patches. We propose the Relative Gap Measure (RGM) to estimate the percentage of gap around a vein, which is defined as the ratio of the overall gap length and the total length of the path that encircles the vein. Additionally, an advanced version of the RGM has been developed to integrate gap quantification estimates from different scar segmentation techniques into a single figure-of-merit. Population-based statistical and regional analysis of gap distribution was performed using a standardised parcellation of the left atrium. We have evaluated our method on synthetic and clinical data from 50 AF patients who underwent PVI with radiofrequency ablation. The population-based analysis concluded that the left superior PV is more prone to lesion gaps while the left inferior PV tends to have less gaps (p < .05 in both cases), in the processed data. This type of information can be very useful for the optimization and objective assessment of PVI interventions., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Aims: We tested the hypothesis that a machine learning (ML) algorithm utilizing both complex echocardiographic data and clinical parameters could be used to phenogroup a heart failure (HF) cohort and identify patients with beneficial response to cardiac resynchronization therapy (CRT)., Methods and Results: We studied 1106 HF patients from the Multicenter Automatic Defibrillator Implantation Trial with Cardiac Resynchronization Therapy (MADIT-CRT) (left ventricular ejection fraction ≤ 30%, QRS ≥ 130 ms, New York Heart Association class ≤ II) randomized to CRT with a defibrillator (CRT-D, n = 677) or an implantable cardioverter defibrillator (ICD, n = 429). An unsupervised ML algorithm (Multiple Kernel Learning and K-means clustering) was used to categorize subjects by similarities in clinical parameters, and left ventricular volume and deformation traces at baseline into mutually exclusive groups. The treatment effect of CRT-D on the primary outcome (all-cause death or HF event) and on volume response was compared among these groups. Our analysis identified four phenogroups, significantly different in the majority of baseline clinical characteristics, biomarker values, measures of left and right ventricular structure and function and the primary outcome occurrence. Two phenogroups included a higher proportion of known clinical characteristics predictive of CRT response, and were associated with a substantially better treatment effect of CRT-D on the primary outcome [hazard ratio (HR) 0.35; 95% confidence interval (CI) 0.19-0.64; P = 0.0005 and HR 0.36; 95% CI 0.19-0.68; P = 0.001] than observed in the other groups (interaction P = 0.02)., Conclusions: Our results serve as a proof-of-concept that, by integrating clinical parameters and full heart cycle imaging data, unsupervised ML can provide a clinically meaningful classification of a phenotypically heterogeneous HF cohort and might aid in optimizing the rate of responders to specific therapies., (© 2018 The Authors. European Journal of Heart Failure © 2018 European Society of Cardiology.)

Aims: Left atrial (LA) fibrosis can be identified by late gadolinium enhancement cardiac magnetic resonance (LGE-CMR) in patients with atrial fibrillation (AF). However, there is limited information about anatomical fibrosis distribution in the left atrium. The aim is to determine whether there is a preferential spatial distribution of fibrosis in the left atrium in patients with AF., Methods and Results: A 3-Tesla LGE-CMR was performed in 113 consecutive patients referred for AF ablation. Images were post-processed and analysed using ADAS-AF software (Galgo Medical), which allows fibrosis identification in 3D colour-coded shells. A regional semiautomatic LA parcellation software was used to divide the atrial wall into 12 segments: 1-4, posterior wall; 5-6, floor; 7, septal wall; 8-11, anterior wall; 12, lateral wall. The presence and amount of fibrosis in each segment was obtained for analysis. After exclusions for artefacts and insufficient image quality, 76 LGE-MRI images (68%) were suitable for fibrosis analysis. Segments 3 and 5, closest to the left inferior pulmonary vein, had significantly higher fibrosis (40.42% ± 23.96 and 25.82% ± 21.24, respectively; P < 0.001), compared with other segments. Segments 8 and 10 in the anterior wall contained the lowest fibrosis (2.54% ± 5.78 and 3.82% ± 11.59, respectively; P < 0.001). Age >60 years was significantly associated with increased LA fibrosis [95% confidence interval (CI) 0.19-8.39, P = 0.04] and persistent AF approached significance (95% CI -0.19% to 7.83%, P = 0.08)., Conclusion: In patients with AF, the fibrotic area is preferentially located at the posterior wall and floor around the antrum of the left inferior pulmonary vein. Age >60 years was associated with increased fibrosis.

Background: In the era of increasingly successful corrective interventions in patients with congenital heart disease (CHD), global and regional myocardial remodeling are emerging as important sources of long-term morbidity/mortality. Changes in organization of the myocardium in CHD, and in its mechanical properties, conduction, and blood supply, result in altered myocardial function both before and after surgery. To gain a better understanding and develop appropriate and individualized treatment strategies, the microscopic organization of cardiomyocytes, and their integration at a macroscopic level, needs to be completely understood. The aim of this study is to describe, for the first time, in 3 dimensions and nondestructively the detailed remodeling of cardiac microstructure present in a human fetal heart with complex CHD., Methods and Results: Synchrotron X-ray phase-contrast imaging was used to image an archival midgestation formalin-fixed fetal heart with right isomerism and complex CHD and compare with a control fetal heart. Analysis of myocyte aggregates, at detail not accessible with other techniques, was performed. Macroanatomic and conduction system changes specific to the disease were clearly observable, together with disordered myocyte organization in the morphologically right ventricle myocardium. Electrical activation simulations suggested altered synchronicity of the morphologically right ventricle., Conclusions: We have shown the potential of X-ray phase-contrast imaging for studying cardiac microstructure in the developing human fetal heart at high resolution providing novel insight while preserving valuable archival material for future study. This is the first study to show myocardial alterations occur in complex CHD as early as midgestation.

During embryogenesis, a mammalian heart develops from a simple tubular shape into a complex 4-chamber organ, going through four distinct phases: early primitive tubular heart, emergence of trabeculations, trabecular remodeling and development of the compact myocardium. In this paper we propose a framework for standardized and subject-independent 3D regional myocardial complexity analysis, applied to analysis of the development of the mouse left ventricle. We propose a standardized subdivision of the myocardium into 3D overlapping regions (in our case 361) and a novel visualization of myocardial complexity, whereupon we: 1) extend the fractal dimension, commonly applied to image slices, to 3D and 2) use volume occupied by the trabeculations in each region together with their surface area, in order to quantify myocardial complexity. The latter provides an intuitive characterization of the complexity, given that compact myocardium will tend to occupy a larger volume with little surface area while high surface area with low volume will correspond to highly trabeculated areas. Using 50 mouse embryo images at 5 different gestational ages (10 subjects per gestational age), we demonstrate how the proposed representation and complexity measures describe the development of LV myocardial complexity. The mouse embryo data was acquired using high resolution episcopic microscopy. The complexity analysis per region was carried out using: 3D fractal dimension, myocardial volume, myocardial surface area and ratio between the two. The analysis of gestational ages was performed on embryos of 14.5, 15.5, 16.5, 17.5 and 18.5 embryonic days, and demonstrated that the regional complexity of the trabeculations increases longitudinally from the base to the apex, with a maximum around the middle. The overall complexity decreases with gestational age, being most complex at 14.5. Circumferentially, at ages 14.5, 15.5 and 16.5, the trabeculations show similar complexity everywhere except for the anteroseptal and inferolateral area of the wall, where it is smaller. At 17.5 days, the regions of high complexity become more localized towards the inferoseptal and anterolateral parts of the wall. At 18.5 days, the high complexity area exhibits further localization at the inferoseptal and anterior part of the wall., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Computational modelling plays an important role in right ventricular (RV) haemodynamic analysis. However, current approaches use smoothed ventricular anatomies. The aim of this study is to characterise RV haemodynamics including detailed endocardial structures like trabeculae, moderator band, and papillary muscles. Four paired detailed and smoothed RV endocardium models (2 male and 2 female) were reconstructed from ex vivo human hearts high-resolution magnetic resonance images. Detailed models include structures with ≥1 mm 2 cross-sectional area. Haemodynamic characterisation was done by computational fluid dynamics simulations with steady and transient inflows, using high-performance computing. The differences between the flows in smoothed and detailed models were assessed using Q-criterion for vorticity quantification, the pressure drop between inlet and outlet, and the wall shear stress. Results demonstrated that detailed endocardial structures increase the degree of intra-ventricular pressure drop, decrease the wall shear stress, and disrupt the dominant vortex creating secondary small vortices. Increasingly turbulent blood flow was observed in the detailed RVs. Female RVs were less trabeculated and presented lower pressure drops than the males. In conclusion, neglecting endocardial structures in RV haemodynamic models may lead to inaccurate conclusions about the pressures, stresses, and blood flow behaviour in the cavity., (© 2018 John Wiley & Sons, Ltd.)

The left atrial appendage (LAA) is a complex and heterogeneous protruding structure of the left atrium (LA). In atrial fibrillation patients, it is the location where 90% of the thrombi are formed. However, the role of the LAA in thrombus formation is not fully known yet. The main goal of this work is to perform a sensitivity analysis to identify the most relevant LA and LAA morphological parameters in atrial blood flow dynamics. Simulations were run on synthetic ellipsoidal left atria models where different parameters were individually studied: pulmonary veins and mitral valve dimensions; LAA shape; and LA volume. Our computational analysis confirmed the relation between large LAA ostia, low blood flow velocities and thrombus formation. Additionally, we found that pulmonary vein configuration exerted a critical influence on LAA blood flow patterns. These findings contribute to a better understanding of the LAA and to support clinical decisions for atrial fibrillation patients., (Copyright © 2018 John Wiley & Sons, Ltd.)

The aim of the present study is to characterize the hemodynamics of left ventricular (LV) geometries to examine the impact of trabeculae and papillary muscles (PMs) on blood flow using high performance computing (HPC). Five pairs of detailed and smoothed LV endocardium models were reconstructed from high-resolution magnetic resonance images (MRI) of ex-vivo human hearts. The detailed model of one LV pair is characterized only by the PMs and few big trabeculae, to represent state of art level of endocardial detail. The other four detailed models obtained include instead endocardial structures measuring ≥1 mm 2 in cross-sectional area. The geometrical characterizations were done using computational fluid dynamics (CFD) simulations with rigid walls and both constant and transient flow inputs on the detailed and smoothed models for comparison. These simulations do not represent a clinical or physiological scenario, but a characterization of the interaction of endocardial structures with blood flow. Steady flow simulations were employed to quantify the pressure drop between the inlet and the outlet of the LVs and the wall shear stress (WSS). Coherent structures were analyzed using the Q-criterion for both constant and transient flow inputs. Our results show that trabeculae and PMs increase the intra-ventricular pressure drop, reduce the WSS and disrupt the dominant single vortex, usually present in the smoothed-endocardium models, generating secondary small vortices. Given that obtaining high resolution anatomical detail is challenging in-vivo , we propose that the effect of trabeculations can be incorporated into smoothed ventricular geometries by adding a porous layer along the LV endocardial wall. Results show that a porous layer of a thickness of 1.2·10 -2 m with a porosity of 20 kg/m 2 on the smoothed-endocardium ventricle models approximates the pressure drops, vorticities and WSS observed in the detailed models.

Background: Cardiovascular magnetic resonance (CMR) imaging has been used to visualise post-ablation atrial scar (PAAS), generally employing a three-dimensional (3D) late gadolinium enhancement (LGE) technique. However the reproducibility of PAAS imaging has not been determined. This cross-over study is the first to investigate the reproducibility of the technique, crucial for both future research design and clinical implementation., Methods: Forty subjects undergoing first time ablation for atrial fibrillation (AF) had detailed CMR assessment of PAAS. Following baseline pre-ablation scan, two scans (separated by 48 h) were performed at three months post-ablation. Each scan session included 3D LGE acquisition at 10, 20 and 30 min post administration of gadolinium-based contrast agent (GBCA). Subjects were allocated at second scan post-ablation to identical imaging parameters ('Repro', n = 10), 3 T scanner ('3 T', n = 10), half-slice thickness ('Half-slice', n = 10) or half GBCA dose ('Half-gad', n = 10). PAAS was compared to baseline scar and then reproducibility was assessed for two measures of thresholded scar (% left atrial (LA) occupied by PAAS (%LA PAAS) and Pulmonary Vein Encirclement (PVE)), and then four measures of non-thresholded scar (point-by-point assessment of PAAS, four normalisation methods). Thresholded measures of PAAS were evaluated against procedural outcome (AF recurrence)., Results: A total of 271 3D acquisitions (out of maximum 280, 96.7%) were acquired. At 20 and 30 min, inter-scan reproducibility was good to excellent (coefficient of variation at 20 min and 30 min: %LA PAAS 0.41 and 0.20; PVE 0.13 and 0.04 respectively for 'Repro' group). Changes in imaging parameters, especially reduced GBCA dose, reduced inter-scan reproducibility, but for most measures remained good to excellent (ICC for %LA PAAS 0.454-0.825, PVE 0.618-0.809 at 30 min). For non-thresholded scar, highest reproducibility was observed using blood pool z-score normalisation technique: inter-scan ICC 0.759 (absolute agreement, 'Repro' group). There was no significant relationship between indices of PAAS and AF recurrence., Conclusion: PAAS imaging is a reproducible finding. Imaging should be performed at least 20 min post-GBCA injection, and a blood pool z-score should be considered for normalisation of signal intensities. The clinical implications of these findings remain to be established in the absence of a simple correlation with arrhythmia outcome., Trial Registration: United Kingdom National Research Ethics Service 08/H0802/68 - 30th September 2008.

Contemporary imaging modalities offer noninvasive quantification of myocardial deformation; however, they make gross assumptions about internal structure of the cardiac walls. Our aim is to study the possible impact of the trabeculations on the stroke volume, strain, and capacity of differently sized ventricles. The cardiac left ventricle is represented by an ellipsoid and the trabeculations by a tissue occupying a fixed volume. The ventricular contraction is modeled by scaling the ellipsoid whereupon the measurements of longitudinal strain, end-diastolic, end-systolic, and stroke volumes are derived and compared. When the trabeculated and nontrabeculated ventricles, having the same geometry and deformation pattern, contain the same amount of blood and contract with the same strain, we observed an increased stroke volume in our model of the trabeculated ventricle. When these ventricles contain and eject the same amount of blood, we observed a reduced strain in the trabeculated case. We identified that a trade-off between the strain and the amount of trabeculations could be reached with a 0.35- to 0.41-cm dense trabeculated layer, without blood filled recesses (for a ventricle with end-diastolic volume of about 150 mL). A trabeculated ventricle can work at lower strains compared to a nontrabeculated ventricle to produce the same stroke volume, which could be a possible explanation why athletes and pregnant women develop reversible signs of left ventricular noncompaction, since the trabeculations could help generating extra cardiac output. This knowledge might help to assess heart failure patients with dilated cardiomyopathies who often show signs of noncompaction., (Copyright © 2017 John Wiley & Sons, Ltd.)

Purpose: Left atrial arrhythmia substrate assessment can involve multiple imaging and electrical modalities, but visual analysis of data on 3D surfaces is time-consuming and suffers from limited reproducibility. Unfold maps (e.g., the left ventricular bull's eye plot) allow 2D visualization, facilitate multimodal data representation, and provide a common reference space for inter-subject comparison. The aim of this work is to develop a method for automatic representation of multimodal information on a left atrial standardized unfold map (LA-SUM)., Methods: The LA-SUM technique was developed and validated using 18 electroanatomic mapping (EAM) LA geometries before being applied to ten cardiac magnetic resonance/EAM paired geometries. The LA-SUM was defined as an unfold template of an average LA mesh, and registration of clinical data to this mesh facilitated creation of new LA-SUMs by surface parameterization., Results: The LA-SUM represents 24 LA regions on a flattened surface. Intra-observer variability of LA-SUMs for both EAM and CMR datasets was minimal; root-mean square difference of 0.008 ± 0.010 and 0.007 ± 0.005 ms (local activation time maps), 0.068 ± 0.063 gs (force-time integral maps), and 0.031 ± 0.026 (CMR LGE signal intensity maps). Following validation, LA-SUMs were used for automatic quantification of post-ablation scar formation using CMR imaging, demonstrating a weak but significant relationship between ablation force-time integral and scar coverage (R 2  = 0.18, P < 0.0001)., Conclusions: The proposed LA-SUM displays an integrated unfold map for multimodal information. The method is applicable to any LA surface, including those derived from imaging and EAM systems. The LA-SUM would facilitate standardization of future research studies involving segmental analysis of the LA.

Background: While individual cardiac myocytes only have a limited ability to shorten, the heart efficiently pumps a large volume-fraction thanks to a cell organization in a complex 3D fibre structure. Subclinical subtle cardiac structural remodelling is often present before symptoms arise. Understanding and early detection of these subtle changes is crucial for diagnosis and prevention. Additionally, personalized computational modelling requires knowledge on the multi-scale structure of the whole heart and vessels., Methods and Results: We developed a rapid acquisition together with visualization and quantification methods of the integrated microstructure of whole in-vitro rodents hearts using synchrotron based X-ray phase-contrast tomography. These images are formed not only by X-ray absorption by the tissue but also by wave propagation phenomena, enhancing structural information, thus allowing to raise tissue contrast to an unprecedented level. We used a (ex-vivo) normal rat heart and fetal rabbit hearts suffering intrauterine growth restriction as a model of subclinical cardiac remodelling to illustrate the strengths and potential of the technique. For comparison, histology and diffusion tensor magnetic resonance imaging was performed., Conclusions: We have developed a novel, high resolution, image acquisition, and quantification approach to study a whole in-vitro heart at myofibre resolution, providing integrated 3D structural information at microscopic level without any need of tissue slicing and processing. This superior imaging approach opens up new possibilities for a systems approach towards analysing cardiac structure and function, providing rapid acquisition of quantitative microstructure of the heart in a near native state., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author 2017. For permissions, please email: journals.permissions@oup.com.)

Reparameterization of surfaces is a widely used tool in computer graphics known mostly from the remeshing algorithms. Recently, the surface reparameterization techniques started to gain popularity in the field of medical imaging, but mostly for convenient 2D visualization of the information initially represented on 3D surfaces (e.g. continuous bulls-eye plot). However, by consistently mapping the 3D information to the same 2D domain, surface reparameterization techniques allow us to put into correspondence anatomical shapes of inherently different geometry. In this paper, we propose a method for anatomical parameterization of cardiac ventricular anatomies that include myocardium, trabeculations, tendons and papillary muscles. The proposed method utilizes a quasi-conformal flattening of the myocardial surfaces of the left and right cardiac ventricles and extending it to cover the interior of the cavities using the local coordinates given by the solution of the Laplace's equation. Subsequently, we define a geometry independent representation for the detailed cardiac left and right ventricular anatomies that can be used for convenient visualization and statistical analysis of the trabeculations in a population. Lastly we show how it can be used for mapping the detailed cardiac anatomy between different hearts, which is of considerable interest for detailed cardiac computational models or shape atlases., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Background: Fetal growth restriction (FGR) affects 5% to 10% of newborns and is associated with increased cardiovascular mortality in adulthood. We evaluated whether prenatal cardiovascular changes previously demonstrated in FGR persist into preadolescence., Methods and Results: A cohort study of 58 FGR (defined as birth weight below 10th centile) and 94 normally grown fetuses identified in utero and followed-up into preadolescence (8-12 years of age) by echocardiography and 3-dimensional shape computational analysis. Compared with controls, FGR preadolescents had a different cardiac shape, with more spherical and smaller hearts. Left ventricular ejection fraction was similar among groups, whereas FGR had decreased longitudinal motion (decreased mitral annular systolic peak velocities: control median, 0.11 m/s [interquartile range, 0.09-0.12] versus FGR median 0.09 m/s [interquartile range, 0.09-0.10]; P<0.01) and impaired relaxation (isovolumic relaxation time: control, 0.21 ms [interquartile range, 0.12-0.35] versus FGR, 0.35 ms [interquartile range, 0.20-0.46]; P=0.04). Global longitudinal strain was decreased (control mean, -22.4% [SD, 1.37] versus FGR mean, -21.5% [SD, 1.16]; P<0.001) compensated by an increased circumferential strain and with a higher prevalence of postsystolic shortening in FGR as compared with controls. These differences persisted after adjustment for parental ethnicity and smoking, prenatal glucocorticoid administration, preeclampsia, gestational age at delivery, days in intensive care unit, sex, age, and body surface area at evaluation., Conclusions: This study provides evidence that cardiac remodeling induced by FGR persists until preadolescence with findings similar to those reported in their prenatal life and childhood. The findings support the hypothesis of primary cardiac programming in FGR for explaining the association between low birth weight and cardiovascular risk in adulthood., (© 2017 American Heart Association, Inc.)

Integration of electrical and structural information for scar characterization in the left ventricle (LV) is a crucial step to better guide radio-frequency ablation therapies, which are usually performed in complex ventricular tachycardia (VT) cases. This integration requires finding a common representation where to map the electrical information from the electro-anatomical map (EAM) surfaces and tissue viability information from delay-enhancement magnetic resonance images (DE-MRI). However, the development of a consistent integration method is still an open problem due to the lack of a proper evaluation framework to assess its accuracy. In this paper we present both: (i) an evaluation framework to assess the accuracy of EAM and imaging integration strategies with simulated EAM data and a set of global and local measures; and (ii) a new integration methodology based on a planar disk representation where the LV surface meshes are quasi-conformally mapped (QCM) by flattening, allowing for simultaneous visualization and joint analysis of the multi-modal data. The developed evaluation framework was applied to estimate the accuracy of the QCM-based integration strategy on a benchmark dataset of 128 synthetically generated ground-truth cases presenting different scar configurations and EAM characteristics. The obtained results demonstrate a significant reduction in global overlap errors (50-100%) with respect to state-of-the-art integration techniques, also better preserving the local topology of small structures such as conduction channels in scars. Data from seventeen VT patients were also used to study the feasibility of the QCM technique in a clinical setting, consistently outperforming the alternative integration techniques in the presence of sparse and noisy clinical data. The proposed evaluation framework has allowed a rigorous comparison of different EAM and imaging data integration strategies, providing useful information to better guide clinical practice in complex cardiac interventions., (Copyright © 2016 Elsevier B.V. All rights reserved.)

The construction of statistical shape models (SSMs) that are rich, i.e., that represent well the natural and complex variability of anatomical structures, is an important research topic in medical imaging. To this end, existing works have addressed the limited availability of training data by decomposing the shape variability hierarchically or by combining statistical and synthetic models built using artificially created modes of variation. In this paper, we present instead a method that merges multiple statistical models of 3D shapes into a single integrated model, thus effectively encoding extra variability that is anatomically meaningful, without the need for the original or new real datasets. The proposed framework has great flexibility due to its ability to merge multiple statistical models with unknown point correspondences. The approach is beneficial in order to re-use and complement pre-existing SSMs when the original raw data cannot be exchanged due to ethical, legal, or practical reasons. To this end, this paper describes two main stages, i.e., (1) statistical model normalization and (2) statistical model integration. The normalization algorithm uses surface-based registration to bring the input models into a common shape parameterization with point correspondence established across eigenspaces. This allows the model fusion algorithm to be applied in a coherent manner across models, with the aim to obtain a single unified statistical model of shape with improved generalization ability. The framework is validated with statistical models of the left and right cardiac ventricles, the L1 vertebra, and the caudate nucleus, constructed at distinct research centers based on different imaging modalities (CT and MRI) and point correspondences. The results demonstrate that the model integration is statistically and anatomically meaningful, with potential value for merging pre-existing multi-modality statistical models of 3D shapes., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Purpose: Highly trained athletes have an increased risk of atrial arrhythmias. Atrial geometrical and functional remodeling may be the underlying substrate. We analyze and relate atrial size, deformation and performance in professional handball players compared with non-sportive subjects., Methods: 24 Professional handball players and 20 non-sportive males were compared. All subjects underwent an echocardiographic study with evaluation of left (LA), right atrial (RA) dimensions and deformation by strain (Sa) and strain rate (SRa). Atrial performance was assessed from the atrial stroke volume (SV). With computational geometrical models, we studied the relation between atrial volumes, strains and SV and compared atrial working conditions. We estimated the functional reserve and a resulting average wall stress., Results: LA and RA volumes were larger in athletes than in controls (35.2 ± 8.8 vs. 24.8 ± 4.3 ml/m(2), p < 0.01 and 29.0 ± 8.4 vs. 19.0 ± 5.1 ml/m(2), p < 0.01 respectively). LASa and RASa during active atrial contraction were decreased in athletes (-12.2 ± 2.0 vs. -14.5 ± 2.1%, p < 0.01 and -12.1 ± 1.8 vs. -14.2 ± 1.5%, p < 0.01 respectively). LASV was similar between groups (6.6 ± 1.4 vs. 7.3 ± 1.1 ml, p = 0.19) and RASV was lower in athletes (6.2 ± 1.3 vs. 7.2 ± 1.1 ml, p < 0.01). Computational models showed that this different operational mode potentially increases performance reserve, but at the cost of higher atrial wall stress., Conclusion: A proportion of athletes with enlarged LA and RA showed different atrial contractile performance, likely resulting in atria working at higher wall stress.

This paper presents a new registration framework for quantifying myocardial motion and strain from the combination of multiple 3D ultrasound (US) sequences. The originality of our approach lies in the estimation of the transformation directly from the input multiple views rather than from a single view or a reconstructed compounded sequence. This allows us to exploit all spatiotemporal information available in the input views avoiding occlusions and image fusion errors that could lead to some inconsistencies in the motion quantification result. We propose a multiview diffeomorphic registration strategy that enforces smoothness and consistency in the spatiotemporal domain by modeling the 4D velocity field continuously in space and time. This 4D continuous representation considers 3D US sequences as a whole, therefore allowing to robustly cope with variations in heart rate resulting in different number of images acquired per cardiac cycle for different views. This contributes to the robustness gained by solving for a single transformation from all input sequences. The similarity metric takes into account the physics of US images and uses a weighting scheme to balance the contribution of the different views. It includes a comparison both between consecutive images and between a reference and each of the following images. The strain tensor is computed locally using the spatial derivatives of the reconstructed displacement fields. Registration and strain accuracy were evaluated on synthetic 3D US sequences with known ground truth. Experiments were also conducted on multiview 3D datasets of 8 volunteers and 1 patient treated by cardiac resynchronization therapy. Strain curves obtained from our multiview approach were compared to the single-view case, as well as with other multiview approaches. For healthy cases, the inclusion of several views improved the consistency of the strain curves and reduced the number of segments where a non-physiological strain pattern was observed. For the patient, the improvement (pacing ON vs. OFF) in synchrony of regional strain correlated with clinician blind assessment and could be seen more clearly when using the multiview approach., (Copyright © 2013 Elsevier B.V. All rights reserved.)

This paper presents a framework for the fusion of multiple point distribution models (PDMs) with unknown point correspondences. With this work, models built from distinct patient groups and imaging modalities can be merged, with the aim to obtain a PDM that encodes a wider range of anatomical variability. To achieve this, two technical challenges are addressed in this work. Firstly, the model fusion must be carried out directly on the corresponding means and eigenvectors as the original data is not always available and cannot be freely exchanged across centers for various legal and practical reasons. Secondly, the PDMs need to be normalized before fusion as the point correspondence is unknown. The proposed framework is validated by integrating statistical models of the left and right ventricles of the heart constructed from different imaging modalities (MRI and CT) and with different landmark representations of the data. The results show that the integration is statistically and anatomically meaningful and that the quality of the resulting model is significantly improved.

Training active shape models requires collecting manual ground-truth meshes in a large image database. While shape information can be reused across multiple imaging modalities, intensity information needs to be imaging modality and protocol specific. In this context, this study has two main purposes: (1) to test the potential of using intensity models learned from MRI simulated datasets and (2) to test the potential of including a measure of reliability during the matching process to increase robustness. We used a population of 400 virtual subjects (XCAT phantom), and two clinical populations of 40 and 45 subjects. Virtual subjects were used to generate simulated datasets (MRISIM simulator). Intensity models were trained both on simulated and real datasets. The trained models were used to segment the left ventricle (LV) and right ventricle (RV) from real datasets. Segmentations were also obtained with and without reliability information. Performance was evaluated with point-to-surface and volume errors. Simulated intensity models obtained average accuracy comparable to inter-observer variability for LV segmentation. The inclusion of reliability information reduced volume errors in hypertrophic patients (EF errors from 17 ± 57% to 10 ± 18%; LV MASS errors from -27 ± 22 g to -14 ± 25 g), and in heart failure patients (EF errors from -8 ± 42% to -5 ± 14%). The RV model of the simulated images needs further improvement to better resemble image intensities around the myocardial edges. Both for real and simulated models, reliability information increased segmentation robustness without penalizing accuracy.

The construction of realistic subject-specific models of the myocardial fiber architecture is relevant to the understanding and simulation of the electromechanical behavior of the heart. This paper presents a statistical approach for the prediction of fiber orientation from myocardial morphology based on the Knutsson mapping. In this space, the orientation of each fiber is represented in a continuous and distance preserving manner, thus allowing for consistent statistical analysis of the data. Furthermore, the directions in the shape space which correlate most with the myocardial fiber orientations are extracted and used for subsequent prediction. With this approach and unlike existing models, all shape information is taken into account in the analysis and the obtained latent variables are statistically optimal to predict fiber orientation in new datasets. The proposed technique is validated based on a sample of canine Diffusion Tensor Imaging (DTI) datasets and the results demonstrate marked improvement in cardiac fiber orientation modeling and prediction.

Active shape models bear a great promise for model-based medical image analysis. Their practical use, though, is undermined due to the need to train such models on large image databases. Automatic building of point distribution models (PDMs) has been successfully addressed and a number of autolandmarking techniques are currently available. However, the need for strategies to automatically build intensity models around each landmark has been largely overlooked in the literature. This work demonstrates the potential of creating intensity models automatically by simulating image generation. We show that it is possible to reuse a 3D PDM built from computed tomography (CT) to segment gated single photon emission computed tomography (gSPECT) studies. Training is performed on a realistic virtual population where image acquisition and formation have been modeled using the SIMIND Monte Carlo simulator and ASPIRE image reconstruction software, respectively. The dataset comprised 208 digital phantoms (4D-NCAT) and 20 clinical studies. The evaluation is accomplished by comparing point-to-surface and volume errors against a proper gold standard. Results show that gSPECT studies can be successfully segmented by models trained under this scheme with subvoxel accuracy. The accuracy in estimated LV function parameters, such as end diastolic volume, end systolic volume, and ejection fraction, ranged from 90.0% to 94.5% for the virtual population and from 87.0% to 89.5% for the clinical population.

In this paper, an unsupervised algorithm, called the Independent Histogram Pursuit (IHP), for segmenting dermatological lesions is proposed. The algorithm estimates a set of linear combinations of image bands that enhance different structures embedded in the image. In particular, the first estimated combination enhances the contrast of the lesion to facilitate its segmentation. Given an N-band image, this first combination corresponds to a line in N dimensions, such that the separation between the two main modes of the histogram obtained by projecting the pixels onto this line, is maximized. The remaining combinations are estimated in a similar way under the constraint of being orthogonal to those already computed. The performance of the algorithm is tested on five different dermatological datasets. The results obtained on these datasets indicate the robustness of the algorithm and its suitability to deal with different types of dermatological lesions. The boundary detection precision using k-means segmentation was close to 97%. The proposed algorithm can be easily combined with the majority of classification algorithms.

This work is framed in the field of statistical face analysis. In particular, the problem of accurate segmentation of prominent features of the face in frontal view images is addressed. We propose a method that generalizes linear Active Shape Models (ASMs), which have already been used for this task. The technique is built upon the development of a nonlinear intensity model, incorporating a reduced set of differential invariant features as local image descriptors. These features are invariant to rigid transformations, and a subset of them is chosen by Sequential Feature Selection for each landmark and resolution level. The new approach overcomes the unimodality and Gaussianity assumptions of classical ASMs regarding the distribution of the intensity values across the training set. Our methodology has demonstrated a significant improvement in segmentation precision as compared to the linear ASM and Optimal Features ASM (a nonlinear extension of the pioneer algorithm) in the tests performed on AR, XM2VTS, and EQUINOX databases.

This paper presents a framework for weighted fusion of several Active Shape and Active Appearance Models. The approach is based on the eigenspace fusion method proposed by Hall et al., which has been extended to fuse more than two weighted eigenspaces using unbiased mean and covariance matrix estimates. To evaluate the performance of fusion, a comparative assessment on segmentation precision as well as facial verification tests are performed using the AR, EQUINOX, and XM2VTS databases. Based on the results, it is concluded that the fusion is useful when the model needs to be updated online or when the original observations are absent.