Vascular Centerline Extraction in 3D MR Angiograms for Phase Contrast MRI Blood Flow Measurement

Automatic Graph-Based Modeling of Brain Microvessels Captured With Two-Photon Microscopy

Graph models of cerebral vasculature derived from two-photon microscopy have shown to be relevant to study brain microphysiology. Automatic graphing of these microvessels remain problematic due to the vascular network complexity and two-photon sensitivity limitations with depth. In this paper, we propose a fully automatic processing pipeline to address this issue. The modeling scheme consists of a fully-convolution neural network to segment microvessels, a three-dimensional surface model generator, and a geometry contraction algorithm to produce graphical models with a single connected component. Based on a quantitative assessment using NetMets metrics, at a tolerance of 60 μm, false negative and false positive geometric error 19 rates are 3.8% and 4.2%, respectively, whereas false nega- 20 tive and false positive topological error rates are 6.1% and 4.5%, respectively. Our qualitative evaluation confirms the efficiency of our scheme in generating useful and accurate graphical models.

Fully Convolutional DenseNets for Segmentation of Microvessels in Two-photon Microscopy

Segmentation of microvessels measured using two-photon microscopy has been studied in the literature with limited success due to uneven intensities associated with optical imaging and shadowing effects. In this work, we address this problem using a customized version of a recently developed fully convolutional neural network, namely, FC-DensNets. To train and validate the network, manual annotations of 8 angiograms from two-photon microscopy was used. Segmentation results are then compared with that of a state-of-the-art scheme that was developed for the same purpose and also based on deep learning. Experimental results show improved performance of used FC-DenseNet in providing accurate and yet end-to-end segmentation of microvessels in two-photon microscopy.

Fully automated calculation of image-derived input function in simultaneous PET/MRI in a sheep model

Background

Obtaining the arterial input function (AIF) from image data in dynamic positron emission tomography (PET) examinations is a non-invasive alternative to arterial blood sampling. In simultaneous PET/magnetic resonance imaging (PET/MRI), high-resolution MRI angiographies can be used to define major arteries for correction of partial-volume effects (PVE) and point spread function (PSF) response in the PET data. The present study describes a fully automated method to obtain the image-derived input function (IDIF) in PET/MRI. Results are compared to those obtained by arterial blood sampling.

Methods

To segment the trunk of the major arteries in the neck, a high-resolution time-of-flight MRI angiography was postprocessed by a vessel-enhancement filter based on the inertia tensor. Together with the measured PSF of the PET subsystem, the arterial mask was used for geometrical deconvolution, yielding the time-resolved activity concentration averaged over a major artery. The method was compared to manual arterial blood sampling at the hind leg of 21 sheep (animal stroke model) during measurement of blood flow with O15-water. Absolute quantification of activity concentration was compared after bolus passage during steady state, i.e., between 2.5- and 5-min post injection. Cerebral blood flow (CBF) values from blood sampling and IDIF were also compared.

Results

The cross-calibration factor obtained by comparing activity concentrations in blood samples and IDIF during steady state is 0.98 ± 0.10. In all examinations, the IDIF provided a much earlier and sharper bolus peak than in the time course of activity concentration obtained by arterial blood sampling. CBF using the IDIF was 22 % higher than CBF obtained by using the AIF yielded by blood sampling.

Conclusions

The small deviation between arterial blood sampling and IDIF during steady state indicates that correction of PVE and PSF is possible with the method presented. The differences in bolus dynamics and, hence, CBF values can be explained by the different sampling locations (hind leg vs. major neck arteries) with differences in delay/dispersion. It will be the topic of further work to test the method on humans with the perspective of replacing invasive blood sampling by an IDIF using simultaneous PET/MRI.

Robust infrarenal aortic aneurysm lumen centerline detection for rupture status classification

The objective of this work is to develop a robust method for human abdominal aortic aneurysm (AAA) centerline detection that can contribute to the accurate computation of features for the prediction of AAA rupture risk. A semiautomatic algorithm is proposed for detecting the lumen centerline in contrast-enhanced abdominal computed tomography images based on online adaboost classifiers, which does not require prior image segmentation. The algorithm was developed and applied to thirty ruptured and thirty unruptured AAA image data sets and the tortuosities of the detected centerline were measured to assess the correlation between AAA tortuosity and the binary ruptured and unruptured labels. The lumen of each data set was segmented manually by a trained radiologist and the resulting centerlines of each data set were defined as the gold standard to evaluate the accuracy of the algorithm and to compare it against two widely used segmentation techniques. The average mean relative accuracy of the offline adaboost classifier is 91.9% with a standard deviation of 1.6%; for the online adaboost classifier it is 93.6% with a standard deviation of 1.9% (p<0.05). The online adaboost classifier outperforms the offline adaboost classifier while their computational costs are similar. Aneurysm tortuosity computed from an accurately derived lumen centerline using online adaboost is statistically higher for ruptured aneurysms compared to unruptured aneurysms, indicating that tortuosity can be used to assess rupture risk in the vascular clinic.

Automated prescription of an optimal imaging plane for measurement of cerebral blood flow by phase contrast magnetic resonance imaging

This study describes and evaluates a semiautomated method for prescribing an optimal imaging plane that is located as close as possible to the skull base, and is simultaneously nearly perpendicular to the four arteries leading blood to the brain [internal carotid arteries (ICAs) and vertebral arteries (VAs)]. Such a method will streamline and improve reliability of the measurement of total cerebral blood flow and intracranial pressure by velocity encoding phase-contrast magnetic resonance imaging. The method first extracts the vessels' centerline from a 2-D time-of-flight magnetic resonance angiogram of the neck by performing distance transformations. An anatomical marker, the V2 segment of the VAs, is then identified to guide the imaging plane to be as close and below the skull base. An imaging plane that is nearly perpendicular to the ICAs and V2 segment of VAs is then identified by minimizing a misalignment value, estimated by a weighted mean of the angles between the plane's normal and the vessel axes at the vessel-plane intersections. The performance of the semiautomated method was evaluated by comparing manually selected planes to those found semiautomatically in nine magnetic resonance angiogram datasets. The semiautomated method consistently outperformed manual prescription with a significantly smaller misalignment value, 8.6° versus 20.7° (P < 0.001), respectively, and significantly improved reproducibility.

Coronary artery segmentation using geometric moments based tracking and snake-driven refinement

Automatic or semi-automatic segmentation and tracking of artery trees from computed tomography angiography (CTA) is an important step to improve the diagnosis and treatment of artery diseases, but it still remains a significant challenging problem. In this paper, we present an artery extraction method to address the challenge. The proposed method consists of two steps: (1) a geometric moments based tracking to secure a rough centerline, and (2) a fully automatic generalized cylinder structure-based snake method to refine the centerlines and estimate the radii of the arteries. In this method, a new line direction based on first and second order geometric moments is adopted while both gradient and intensity information are used in the snake model to improve the accuracy. The approach has been evaluated on synthetic images as well as 8 clinical coronary CTA images with 32 coronary arteries. Our method achieves 94.7% overlap tracking ability within an average distance inside the vessel of 0.36 mm.