Synchrotron-based visualization and segmentation of elastic lamellae in the mouse carotid artery during quasi-static pressure inflation

This dataset contains images that were obtained during quasi-static pressure inflation of mouse carotid arteries. Images were taken with phase propagation imaging at the X02DA TOMCAT beamline of the Swiss Light Source synchrotron at the Paul Scherrer Institute in Villigen, Switzerland. Scans of n=12 left carotid arteries (n-6 Apoe-deficient mice, n=6 wild-type mice, all on a C57Bl6J background) were taken at pressure levels of 0, 10, 20, 30, 40, 50, 70, 90 and 120 mmHg. For analysis we selected 75 images from the center of each stack (starting at the center of the stack, and skipping 2 of every three images in both cranial and caudal axial directions) for each sample and for each pressure level, resulting in a total of 75 x 12 x 9 = 8100 analyzed images from 108 different scans. Segmentation, 3D visualization and geometric analysis is presented in the corresponding manuscript. Files are uploaded in 16bit .tif format and are named: mouseid_pressurelevel_stacknumber, with mouseid consisting of either Apoe (Apoe-deficient) or Bl (wild-type) and the mouse number, pressurelevel varies from P0 to P120 and stacknumber indicates which image from the stack has been uploaded

Semiautomated Skeletonization of the Pulmonary Arterial Tree in Micro-CT Images

We present a simple and robust approach that utilizes planar images at different angular rotations combined with unfiltered back-projection to locate the central axes of the pulmonary arterial tree. Three-dimensional points are selected interactively by the user. The computer calculates a sub- volume unfiltered back-projection orthogonal to the vector connecting the two points and centered on the first point. Because more x-rays are absorbed at the thickest portion of the vessel, in the unfiltered back-projection, the darkest pixel is assumed to be the center of the vessel. The computer replaces this point with the newly computer-calculated point. A second back-projection is calculated around the original point orthogonal to a vector connecting the newly-calculated first point and user-determined second point. The darkest pixel within the reconstruction is determined. The computer then replaces the second point with the XYZ coordinates of the darkest pixel within this second reconstruction. Following a vector based on a moving average of previously determined 3- dimensional points along the vessel\u27s axis, the computer continues this skeletonization process until stopped by the user. The computer estimates the vessel diameter along the set of previously determined points using a method similar to the full width-half max algorithm. On all subsequent vessels, the process works the same way except that at each point, distances between the current point and all previously determined points along different vessels are determined. If the difference is less than the previously estimated diameter, the vessels are assumed to branch. This user/computer interaction continues until the vascular tree has been skeletonized

Contrast-enhanced micro-CT imaging in murine carotid arteries : a new protocol for computing wall shear stress

Background: Wall shear stress (WSS) is involved in the pathophysiology of atherosclerosis. The correlation between WSS and atherosclerosis can be investigated over time using a WSS-manipulated atherosclerotic mouse model. To determine WSS in vivo, detailed 3D geometry of the vessel network is required. However, a protocol to reconstruct 3D murine vasculature using this animal model is lacking. In this project, we evaluated the adequacy of eXIA 160, a small animal contrast agent, for assessing murine vascular network on micro-CT. Also, a protocol was established for vessel geometry segmentation and WSS analysis. Methods: A tapering cast was placed around the right common carotid artery (RCCA) of ApoE(-/-) mice (n = 8). Contrast-enhanced micro-CT was performed using eXIA 160. An innovative local threshold-based segmentation procedure was implemented to reconstruct 3D geometry of the RCCA. The reconstructed RCCA was compared to the vessel geometry using a global threshold-based segmentation method. Computational fluid dynamics was applied to compute the velocity field and WSS distribution along the RCCA. Results: eXIA 160-enhanced micro-CT allowed clear visualization and assessment of the RCCA in all eight animals. No adverse biological effects were observed from the use of eXIA 160. Segmentation using local threshold values generated more accurate RCCA geometry than the global threshold-based approach. Mouse-specific velocity data and the RCCA geometry generated 3D WSS maps with high resolution, enabling quantitative analysis of WSS. In all animals, we observed low WSS upstream of the cast. Downstream of the cast, asymmetric WSS patterns were revealed with variation in size and location between animals. Conclusions: eXIA 160 provided good contrast to reconstruct 3D vessel geometry and determine WSS patterns in the RCCA of the atherosclerotic mouse model. We established a novel local threshold-based segmentation protocol for RCCA reconstruction and WSS computation. The observed differences between animals indicate the necessity to use mouse-specific data for WSS analysis. For our future work, our protocol makes it possible to study in vivo WSS longitudinally over a growing plaque

Initial steps towards automatic segmentation of the wire frame of stent grafts in CT data

For the purpose of obtaining a geometrical model of the wire frame of stent grafts, we propose three tracking methods to segment the stent's wire, and compare them in an experiment. A 2D test image was created by obtaining a projection of a 3D volume containing a stent. The image was modified to connect the parts of the stent's frame and thus create a single path. Ten versions of this image were obtained by adding different noise realizations. Each algorithm was started at the start of each of the ten images, after which the traveled paths were compared to the known correct path to determine the performance. Additionally, the algorithms were applied to 3D clinical data and visually inspected. The method based on the minimum cost path algorithm scored excellent in the experiment and showed good results on the 3D data. Future research will focus on establishing a geometrical model by determining the corner points and the crossings from the results of this method.\u

A Survey on Deep Learning in Medical Image Analysis

Deep learning algorithms, in particular convolutional networks, have rapidly become a methodology of choice for analyzing medical images. This paper reviews the major deep learning concepts pertinent to medical image analysis and summarizes over 300 contributions to the field, most of which appeared in the last year. We survey the use of deep learning for image classification, object detection, segmentation, registration, and other tasks and provide concise overviews of studies per application area. Open challenges and directions for future research are discussed.Comment: Revised survey includes expanded discussion section and reworked introductory section on common deep architectures. Added missed papers from before Feb 1st 201

Automation Process for Morphometric Analysis of Volumetric CT Data from Pulmonary Vasculature in Rats

With advances in medical imaging scanners, it has become commonplace to generate large multidimensional datasets. These datasets require tools for a rapid, thorough analysis. To address this need, we have developed an automated algorithm for morphometric analysis incorporating A Visualization Workshop computational and image processing libraries for three-dimensional segmentation, vascular tree generation and structural hierarchical ordering with a two-stage numeric optimization procedure for estimating vessel diameters. We combine this new technique with our mathematical models of pulmonary vascular morphology to quantify structural and functional attributes of lung arterial trees. Our physiological studies require repeated measurements of vascular structure to determine differences in vessel biomechanical properties between animal models of pulmonary disease. Automation provides many advantages including significantly improved speed and minimized operator interaction and biasing. The results are validated by comparison with previously published rat pulmonary arterial micro-CT data analysis techniques, in which vessels were manually mapped and measured using intense operator intervention

Ultra-high-resolution 3D imaging of atherosclerosis in mice with synchrotron differential phase contrast: a proof of concept study.

The goal of this study was to investigate the performance of 3D synchrotron differential phase contrast (DPC) imaging for the visualization of both macroscopic and microscopic aspects of atherosclerosis in the mouse vasculature ex vivo. The hearts and aortas of 2 atherosclerotic and 2 wild-type control mice were scanned with DPC imaging with an isotropic resolution of 15 μm. The coronary artery vessel walls were segmented in the DPC datasets to assess their thickness, and histological staining was performed at the level of atherosclerotic plaques. The DPC imaging allowed for the visualization of complex structures such as the coronary arteries and their branches, the thin fibrous cap of atherosclerotic plaques as well as the chordae tendineae. The coronary vessel wall thickness ranged from 37.4 ± 5.6 μm in proximal coronary arteries to 13.6 ± 3.3 μm in distal branches. No consistent differences in coronary vessel wall thickness were detected between the wild-type and atherosclerotic hearts in this proof-of-concept study, although the standard deviation in the atherosclerotic mice was higher in most segments, consistent with the observation of occasional focal vessel wall thickening. Overall, DPC imaging of the cardiovascular system of the mice allowed for a simultaneous detailed 3D morphological assessment of both large structures and microscopic details

A framework for computational fluid dynamic analyses of patient-specific stented coronary arteries from optical coherence tomography images

The clinical challenge of percutaneous coronary interventions (PCI) is highly dependent on the recognition of the coronary anatomy of each individual. The classic imaging modality used for PCI is angiography, but advanced imaging techniques that are routinely performed during PCI, like optical coherence tomography (OCT), may provide detailed knowledge of the pre-intervention vessel anatomy as well as the post-procedural assessment of the specific stent-to-vessel interactions. Computational fluid dynamics (CFD) is an emerging investigational tool in the setting of optimization of PCI results. In this study, an OCT-based reconstruction method was developed for the execution of CFD simulations of patient-specific coronary artery models which include the actual geometry of the implanted stent. The method was applied to a rigid phantom resembling a stented segment of the left anterior descending coronary artery. The segmentation algorithm was validated against manual segmentation. A strong correlation was found between automatic and manual segmentation of lumen in terms of area values. Similarity indices resulted >96% for the lumen segmentation and >77% for the stent strut segmentation. The 3D reconstruction achieved for the stented phantom was also assessed with the geometry provided by X-ray computed micro tomography scan, used as ground truth, and showed the incidence of distortion from catheter-based imaging techniques. The 3D reconstruction was successfully used to perform CFD analyses, demonstrating a great potential for patient-specific investigations. In conclusion, OCT may represent a reliable source for patient-specific CFD analyses which may be optimized using dedicated automatic segmentation algorithms

Automation of Hessian-Based Tubularity Measure Response Function in 3D Biomedical Images

The blood vessels and nerve trees consist of tubular objects interconnected into a complex tree- or web-like structure that has a range of structural scale 5 μm diameter capillaries to 3 cm aorta. This large-scale range presents two major problems; one is just making the measurements, and the other is the exponential increase of component numbers with decreasing scale. With the remarkable increase in the volume imaged by, and resolution of, modern day 3D imagers, it is almost impossible to make manual tracking of the complex multiscale parameters from those large image data sets. In addition, the manual tracking is quite subjective and unreliable. We propose a solution for automation of an adaptive nonsupervised system for tracking tubular objects based on multiscale framework and use of Hessian-based object shape detector incorporating National Library of Medicine Insight Segmentation and Registration Toolkit (ITK) image processing libraries