Alejandro F. Frangi

A non-linear gray-level appearance model improves active shape model segmentation

Active Shape Models (ASMs), a knowledge-based segmentation algorithm developed by Cootes and Taylor [1995, 1999], have become a standard and popular method for detecting structures in medical images. In ASMs-and various comparable approaches-the model of the objects shape and of its gray-level variations is based the assumption of linear distributions. In this work, we explore a new way to model the gray-level appearance of the objects, using a k-nearest-neighbors (kNN) classifier and a set of selected features for each location and resolution of the Active Shape Model. The construction of the kNN classifier and the selection of features from training images is fully automatic. We compare our approach with the standard ASMs on synthetic data and in four medical segmentation tasks. In all cases, the new method produces significantly better results (p<0.001).

Simultaneous Super-Resolution and Cross-Modality Synthesis of 3D Medical Images Using Weakly-Supervised Joint Convolutional Sparse Coding

Magnetic Resonance Imaging (MRI) offers high-resolution in vivo imaging and rich functional and anatomical multimodality tissue contrast. In practice, however, there are challenges associated with considerations of scanning costs, patient comfort, and scanning time that constrain how much data can be acquired in clinical or research studies. In this paper, we explore the possibility of generating high-resolution and multimodal images from low-resolution single-modality imagery. We propose the weakly-supervised joint convolutional sparse coding to simultaneously solve the problems of super-resolution (SR) and cross-modality image synthesis. The learning process requires only a few registered multimodal image pairs as the training set. Additionally, the quality of the joint dictionary learning can be improved using a larger set of unpaired images. To combine unpaired data from different image resolutions/modalities, a hetero-domain image alignment term is proposed. Local image neighborhoods are naturally preserved by operating on the whole image domain (as opposed to image patches) and using joint convolutional sparse coding. The paired images are enhanced in the joint learning process with unpaired data and an additional maximum mean discrepancy term, which minimizes the dissimilarity between their feature distributions. Experiments show that the proposed method outperforms state-of-the-art techniques on both SR reconstruction and simultaneous SR and cross-modality synthesis.

Segmentation and Quantification for Angle-Closure Glaucoma Assessment in Anterior Segment OCT

Angle-closure glaucoma is a major cause of irreversible visual impairment and can be identified by measuring the anterior chamber angle (ACA) of the eye. The ACA can be viewed clearly through anterior segment optical coherence tomography (AS-OCT), but the imaging characteristics and the shapes and locations of major ocular structures can vary significantly among different AS-OCT modalities, thus complicating image analysis. To address this problem, we propose a data-driven approach for automatic AS-OCT structure segmentation, measurement, and screening. Our technique first estimates initial markers in the eye through label transfer from a hand-labeled exemplar data set, whose images are collected over different patients and AS-OCT modalities. These initial markers are then refined by using a graph-based smoothing method that is guided by AS-OCT structural information. These markers facilitate segmentation of major clinical structures, which are used to recover standard clinical parameters. These parameters can be used not only to support clinicians in making anatomical assessments, but also to serve as features for detecting anterior angle closure in automatic glaucoma screening algorithms. Experiments on Visante AS-OCT and Cirrus high-definition-OCT data sets demonstrate the effectiveness of our approach.

Semi-supervised Assessment of Incomplete LV Coverage in Cardiac MRI Using Generative Adversarial Nets

Cardiac magnetic resonance (CMR) images play a growing role in diagnostic imaging of cardiovascular diseases. Ensuring full coverage of the Left Ventricle (LV) is a basic criteria of CMR image quality. Complete LV coverage, from base to apex, precedes accurate cardiac volume and functional assessment. Incomplete coverage of the LV is identified through visual inspection, which is time-consuming and usually done retrospectively in large imaging cohorts. In this paper, we propose a novel semi-supervised method to check the coverage of LV from CMR images by using generative adversarial networks (GAN), we call it Semi-Coupled-GANs (SCGANs). To identify missing basal and apical slices in a CMR volume, a two-stage framework is proposed. First, the SCGANs generate adversarial examples and extract high-level features from the CMR images; then these image attributes are used to detect missing basal and apical slices. We constructed extensive experiments to validate the proposed method on UK Biobank with more than 6000 independent volumetric MR scans, which achieved high accuracy and robust results for missing slice detection, comparable with those of state of the art deep learning methods. The proposed method, in principle, can be adapted to other CMR image data for LV coverage assessment.

Machine-learning Support to Individual Diagnosis of Mild Cognitive Impairment Using Multimodal MRI and Cognitive Assessments

Background: Understanding whether the cognitive profile of a patient indicates mild cognitive impairment (MCI) or performance levels within normality is often a clinical challenge. The use of resting-state functional magnetic resonance imaging (RS-fMRI) and machine learning may represent valid aids in clinical settings for the identification of MCI patients. Methods: Machine-learning models were computed to test the classificatory accuracy of cognitive, volumetric [structural magnetic resonance imaging (sMRI)] and blood oxygen level dependent-connectivity (extracted from RS-fMRI) features, in single-modality and mixed classifiers. Results: The best and most significant classifier was the RS-fMRI+Cognitive mixed classifier (94% accuracy), whereas the worst performing was the sMRI classifier (∼80%). The mixed global (sMRI+RS-fMRI+Cognitive) had a slightly lower accuracy (∼90%), although not statistically different from the mixed RS-fMRI+Cognitive classifier. The most important cognitive features were indices of declarative memory and semantic processing. The crucial volumetric feature was the hippocampus. The RS-fMRI features selected by the algorithms were heavily based on the connectivity of mediotemporal, left temporal, and other neocortical regions. Conclusion: Feature selection was profoundly driven by statistical independence. Some features showed no between-group differences, or showed a trend in either direction. This indicates that clinically relevant brain alterations typical of MCI might be subtle and not inferable from group analysis.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre‐operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient‐specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error‐free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post‐treatment intra‐aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random‐effects meta‐analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta‐analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.1385

Tracking and diameter estimation of retinal vessels using Gaussian process and Radon transform

Abstract. Extraction of blood vessels in retinal images is an important step for computer-aided diagnosis of ophthalmic pathologies. We propose an approach for blood vessel tracking and diameter estimation. We hypothesize that the curvature and the diameter of blood vessels are Gaussian processes (GPs). Local Radon transform, which is robust against noise, is subsequently used to compute the features and train the GPs. By learning the kernelized covariance matrix from training data, vessel direction and its diameter are estimated. In order to detect bifurcations, multiple GPs are used and the difference between their corresponding predicted directions is quantified. The combination of Radon features and GP results in a good performance in the presence of noise. The proposed method successfully deals with typically difficult cases such as bifurcations and central arterial reflex, and also tracks thin vessels with high accuracy. Experiments are conducted on the publicly available DRIVE, STARE, CHASEDB1, and high-resolution fundus databases evaluating sensitivity, specificity, and Matthew’s correlation coefficient (MCC). Experimental results on these datasets show that the proposed method reaches an average sensitivity of 75.67%, specificity of 97.46%, and MCC of 72.18% which is comparable to the state-of-the-art.

ApoE ε4 Allele Related Alterations in Hippocampal Connectivity in Early Alzheimer’s Disease Support Memory Performance

BACKGROUNDnWhether the presence of the Apolipoprotein E ε4 allele modulates hippocampal connectivity networks in abnormal ageing has yet to be fully clarified.nnnOBJECTIVEnAllele-dependent differences in this pattern of functional connectivity were investigated in patients with very mild neurodegeneration of the Alzheimers type, carriers and non-carriers of the ε4 allele.nnnMETHODnA seed-based connectivity approach was used. The two groups were similar in demographics, volumetric measures of brain structure, and cognitive profiles.nnnRESULTSnε4-carriers had increased connectivity between the seed area in the left hippocampus and 1) a left insular/lateral prefrontal region and 2) the contralateral right parietal cortex. Moreover, hippocampus- to-parietal connectivity in the group of ε4 carriers was positively associated with memory performance, indicating that the between-group difference reflects compensatory processes. Retrospective analyses of functional connectivity based on patients from the ADNI initiative confirmed this pattern.nnnCONCLUSIONnWe suggest that increased connectivity with areas external to the Default Mode Network (DMN) reflects both compensatory recruitment of additional areas, and pathological interwining between the DMN and the salience network as part of a global ε4-dependent circuital disruption. These differences indicate that the ε4 allele is associated with a more profound degree of DMN network breakdown even in the prodromal stage of neurodegeneration.

Support for Taverna workflows in the VPH-Share cloud platform

BACKGROUND AND OBJECTIVEnTo address the increasing need for collaborative endeavours within the Virtual Physiological Human (VPH) community, the VPH-Share collaborative cloud platform allows researchers to expose and share sequences of complex biomedical processing tasks in the form of computational workflows. The Taverna Workflow System is a very popular tool for orchestrating complex biomedical & bioinformatics processing tasks in the VPH community. This paper describes the VPH-Share components that support the building and execution of Taverna workflows, and explains how they interact with other VPH-Share components to improve the capabilities of the VPH-Share platform.nnnMETHODSnTaverna workflow support is delivered by the Atmosphere cloud management platform and the VPH-Share Taverna plugin. These components are explained in detail, along with the two main procedures that were developed to enable this seamless integration: workflow composition and execution.nnnRESULTSn1) Seamless integration of VPH-Share with other components and systems. 2) Extended range of different tools for workflows. 3) Successful integration of scientific workflows from other VPH projects. 4) Execution speed improvement for medical applications.nnnCONCLUSIONnThe presented workflow integration provides VPH-Share users with a wide range of different possibilities to compose and execute workflows, such as desktop or online composition, online batch execution, multithreading, remote execution, etc. The specific advantages of each supported tool are presented, as are the roles of Atmosphere and the VPH-Share plugin within the VPH-Share project. The combination of the VPH-Share plugin and Atmosphere engenders the VPH-Share infrastructure with far more flexible, powerful and usable capabilities for the VPH-Share community. As both components can continue to evolve and improve independently, we acknowledge that further improvements are still to be developed and will be described.

Three-dimensional reconstruction and NURBS-based structured meshing of coronary arteries from the conventional X-ray angiography projection images

Despite its two-dimensional nature, X-ray angiography (XRA) has served as the gold standard imaging technique in the interventional cardiology for over five decades. Accordingly, demands for tools that could increase efficiency of the XRA procedure for the quantitative analysis of coronary arteries (CA) are constantly increasing. The aim of this study was to propose a novel procedure for three-dimensional modeling of CA from uncalibrated XRA projections. A comprehensive mathematical model of the image formation was developed and used with a robust genetic algorithm optimizer to determine the calibration parameters across XRA views. The frames correspondences between XRA acquisitions were found using a partial-matching approach. Using the same matching method, an efficient procedure for vessel centerline reconstruction was developed. Finally, the problem of meshing complex CA trees was simplified to independent reconstruction and meshing of connected branches using the proposed nonuniform rational B-spline (NURBS)-based method. Because it enables structured quadrilateral and hexahedral meshing, our method is suitable for the subsequent computational modelling of CA physiology (i.e. coronary blood flow, fractional flow reverse, virtual stenting and plaque progression). Extensive validations using digital, physical, and clinical datasets showed competitive performances and potential for further application on a wider scale.