Improved centerline tree detection of diseased peripheral arteries with a cascading algorithm for vascular segmentation.

Vascular segmentation plays an important role in the assessment of peripheral arterial disease. The segmentation is very challenging especially for arteries with severe stenosis or complete occlusion. We present a cascading algorithm for vascular centerline tree detection specializing in detecting centerlines in diseased peripheral arteries. It takes a three-dimensional computed tomography angiography (CTA) volume and returns a vascular centerline tree, which can be used for accelerating and facilitating the vascular segmentation. The algorithm consists of four levels, two of which detect healthy arteries of varying sizes and two that specialize in different types of vascular pathology: severe calcification and occlusion. We perform four main steps at each level: appropriate parameters for each level are selected automatically, a set of centrally located voxels is detected, these voxels are connected together based on the connection criteria, and the resulting centerline tree is corrected from spurious branches. The proposed method was tested on 25 CTA scans of the lower limbs, achieving an average overlap rate of 89% and an average detection rate of 82%. The average execution time using four CPU cores was 70 s, and the technique was successful also in detecting very distal artery branches, e.g., in the foot.

Fast level-set based image segmentation using coherent propagation.

PURPOSE: The level-set method is known to require long computation time for 3D image segmentation, which limits its usage in clinical workflow. The goal of this study was to develop a fast level-set algorithm based on the coherent propagation method and explore its character using clinical datasets. METHODS: The coherent propagation algorithm allows level set functions to converge faster by forcing the contour to move monotonically according to a predicted developing trend. Repeated temporary backwards propagation, caused by noise or numerical errors, is then avoided. It also makes it possible to detect local convergence, so that the parts of the boundary that have reached their final position can be excluded in subsequent iterations, thus reducing computation time. To compensate for the overshoot error, forward and backward coherent propagation is repeated periodically. This can result in fluctuations of great magnitude in parts of the contour. In this paper, a new gradual convergence scheme using a damping factor is proposed to address this problem. The new algorithm is also generalized to non-narrow band cases. Finally, the coherent propagation approach is combined with a new distance-regularized level set, which eliminates the needs of reinitialization of the distance. RESULTS: Compared with the sparse field method implemented in the widely available ITKSnap software, the proposed algorithm is about 10 times faster when used for brain segmentation and about 100 times faster for aorta segmentation. Using a multiresolution approach, the new method achieved 50 times speed-up in liver segmentation. The Dice coefficient between the proposed method and the sparse field method is above 99% in most cases. CONCLUSIONS: A generalized coherent propagation algorithm for level set evolution yielded substantial improvement in processing time with both synthetic datasets and medical images.

A modular simulation framework for colonoscopy using a new haptic device.

We have developed a multi-threaded framework for colonoscopy simulation utilising OpenGL with an interface to a real-time prototype colonoscopy haptic device. A modular framework has enabled us to support multiple haptic devices and efficiently integrate new research into physically based modelling of the colonoscope, colon and surrounding organs. The framework supports GPU accelerated algorithms as runtime modules, allowing the real-time calculations required for haptic feedback.

Identification of the main arterial branches by whole-body contrast-enhanced MRA in elderly subjects using limited user interaction and fast marching.

PURPOSE: To extract a graph model corresponding to a predefined set of arterial branches from whole-body contrast-enhanced magnetic resonance angiography (CE-MRA) data sets in elderly asymptomatic subjects, a high-incidence group. MATERIALS AND METHODS: Maximum intensity projections (MIPs) were used as an interface to place landmarks in the three-dimensional (3D) data sets. These landmarks were linked together using fast marching to form a graph model of the arterial tree. Only vessels of interest were identified. RESULTS: We tested our method on 10 subjects. We were able to build a graph model of the main arterial branches that performed well in the presence of vascular pathologies, such as stenosis and aneurysm. The results were rated by an experienced radiologist, with an overall success rate of 80%. CONCLUSION: We were able to extract chosen arterial branches in 3D whole-body CE-MRA images with a moderate amount of interaction using a single MIP projection.

Whole-body T1 mapping improves the definition of adipose tissue: consequences for automated image analysis.

PURPOSE: To determine whether a whole-body T1-mapping acquisition method improves the definition of adipose tissue (AT) and simplifies automated AT segmentation compared to an image-based method. MATERIALS AND METHODS: The study included 10 subjects. Two whole-body volumes were acquired from each subject using two different flip angles. Whole-body T1 maps were calculated from each pair of whole-body volumes. AT was automatically segmented from the T1 maps and from the original image slices. The results were evaluated using manually segmented slices as reference. RESULTS: The T1-mapping method segmented more of the reference AT than the image-based method, with mean values (standard deviations (SDs)) of 87.7(5.1)% and 81.1(5.2)%, respectively. Compared to the image-based method, the T1-mapping method gives better histogram separation of AT in whole-body volumes. The suggested method also provides an output with smaller in-slice AT intensity variations. CONCLUSION: The T1-mapping method improves the definition of AT. T1-based analysis is superior to analysis based on the original images, and allows fully automated and accurate whole-body AT segmentation.

Centerline-based colon segmentation for CT colonography.

We have developed a fully automated algorithm for colon segmentation, centerline-based segmentation (CBS), which is faster than any of the previously presented segmentation algorithms, but also has high sensitivity as well as high specificity. The algorithm first thresholds a set of unprocessed CT slices. Outer air is removed, after which a bounding box is computed. A centerline is computed for all remaining regions in the thresholded volume, disregarding segments related to extracolonic structures. Centerline segments are connected, after which the anatomy-based removal of segments representing extracolonic structures occurs. Segments related to the remaining centerline are locally region grown, and the colonic wall is found by dilation. Shape-based interpolation provides an isotropic mask. For 38 CT datasets, CBS was compared with the knowledge-guided segmentation (KGS) algorithm for sensitivity and specificity. With use of a 1.5 GHz AMD Athlon-based PC, the average computation time for the segmentation was 14.8 s. The sensitivity was, on average, 96%, and the specificity was 99%. A total of 21% of the voxels segmented by KGS, of which 96% represented extracolonic structures and 4% represented the colon, were removed.

Region-based supine-prone correspondence for the reduction of false-positive CAD polyp candidates in CT colonography.

RATIONALE AND OBJECTIVES: Radiologists often compare the supine and prone data sets of a patient to confirm potential polyp findings in computed tomographic (CT) colonography (CTC). We developed a new automated method that uses region-based supine-prone correspondence for the reduction of false-positive (FP) polyp candidates in computer-aided detection (CAD) for CTC. MATERIALS AND METHODS: Up to six anatomic landmarks are established by use of the extracted region of the colonic lumen. A region-growing scheme with distance calculations is used to divide the colonic lumen into overlapping segments that match in the supine and prone data sets. Polyp candidates detected by means of a CAD scheme are eliminated in colonic segments that have sufficient diagnostic quality and contain polyp candidates in only one of the data sets of a patient. The method was evaluated with 121 CTC cases, including 42 polyps of 5 mm or greater in 28 patients, obtained by use of single- and multidetector CT scanners with standard pre-colonoscopy cleansing. RESULTS: Complete or partial correspondence was established in 71% of cases. Based on a leave-one-patient-out evaluation, application of the method reduced 19% of FP results reported by our CAD scheme at a 90.5% by-polyp detection sensitivity, without loss of any true-positive results. The resulting CAD scheme yielded 2.4 FP results per patient, on average, with the use of the correspondence method, whereas it yielded 3.0 FP results per patient without the use of the method. CONCLUSION: The correspondence method is potentially useful for improving the specificity of CAD in CTC.

Virtual endoscopic visualization of the colon by shape-scale signatures.

We developed a new visualization method for virtual endoscopic examination of computed tomographic (CT) colonographic data by use of shape-scale analysis. The method provides each colonic structure of interest with a unique color, thereby facilitating rapid diagnosis of the colon. Two shape features, called the local shape index and curvedness, are used for defining the shape-scale spectrum. When we map the shape index and curvedness values within CT colonographic data to the shape-scale spectrum, specific types of colonic structures are represented by unique characteristic signatures in the spectrum. The characteristic signatures of specific types of lesions can be determined by use of computer-simulated lesions or by use of clinical data sets subjected to a computerized detection scheme. The signatures are used for defining a two-dimensional color map by assignment of a unique color to each signature region. The method was evaluated visually by use of computer-simulated lesions and clinical CT colonographic data sets, as well as by an evaluation of the human observer performance in the detection of polyps without and with the use of the color maps. The results indicate that the coloring of the colon yielded by the shape-scale color maps can be used for differentiating among the chosen colonic structures. Moreover, the results indicate that the use of the shape-scale color maps can improve the performance of radiologists in the detection of polyps in CT colonography.

Fast and robust computation of colon centerline in CT colonography.

Although several methods for generating the centerline of a colon from CT colonographic scans have been proposed, in general they are time-consuming and do not take into account that the images of the colon may be of nonoptimal quality, with collapsed regions, and stool within the colon. Furthermore, the colonic lumen or wall, which is often used as a basis for computation of a centerline, is not always precisely segmented. In this study, we have developed an algorithm for computation of a colon centerline that is fast compared to the centerline algorithms presented in the reviewed literature, and that relies little on a complete colon segments identification. The proposed algorithm first extracts local maxima in a distance map of a segmented colonic lumen. The maxima are considered to be nodes in a set of graphs, and are iteratively linked together, based on a set of connection criteria, giving a minimum distance spanning tree. The connection criteria are computed from the distance from object boundary, the Euclidean distance between nodes and the voxel values on the pathway between pairs of nodes. After the last iteration, redundant branches are removed and end segments are recovered for each remaining graph. A subset of the initial maxima is used for distinguishing between the colon and noncolonic centerline segments among the set of graphs, giving the final centerline representation. A phantom study showed that, with respect to phantom variations, the algorithm achieved nearly constant computation time (2.3-2.9 s) except for the most extreme setting (20.2 s). The algorithm successfully found all, or most of, the centerline (93% - 100%). Displacement from optimum varied with colon diameter (1.2-6.6 mm). By use of 40 CT colonographic scans, the computer-generated centerlines were compared with the centerlines generated by three radiologists. The similarity was measured based on percent coverage and average displacement. The computer-generated centerlines, when compared with human-generated centerlines, had approximately the same displacement as when the human-generated centerlines were compared among each other (3.8 mm versus 4.0 mm). The coverage of the computer-generated centerlines was slightly less than that of the human-generated centerlines (92% versus 94%). The 40 centerlines were, on average, computed in 10.5 seconds, including computation time for the distance transform, with an Intel Pentium-based 800 MHz computer, as compared with 12-17 seconds or more (excluding computation time for the distance transform needed) per centerline as reported in other studies.

Computerized detection of colorectal masses in CT colonography based on fuzzy merging and wall-thickening analysis.

In recent years, several computer-aided detection (CAD) schemes have been developed for the detection of polyps in CT colonography (CTC). However, few studies have addressed the problem of computerized detection of colorectal masses in CTC. This is mostly because masses are considered to be well visualized by a radiologist because of their size and invasiveness. Nevertheless, the automated detection of masses would naturally complement the automated detection of polyps in CTC and would produce a more comprehensive computer aid to radiologists. Therefore, in this study, we identified some of the problems involved with the computerized detection of masses, and we developed a scheme for the computerized detection of masses that can be integrated into a CAD scheme for the detection of polyps. The performance of the mass detection scheme was evaluated by application to clinical CTC data sets. CTC was performed on 82 patients with helical CT scanners and reconstruction intervals of 1.0-5.0 mm in the supine and prone positions. Fourteen patients (17%) had a total of 14 masses of 30-50 mm, and sixteen patients (20%) had a total of 30 polyps 5-25 mm in diameter. Four patients had both polyps and masses. Fifty-six of the patients (68%) were normal. The CTC data were interpolated linearly to yield isotropic data sets, and the colon was extracted by use of a knowledge-guided segmentation technique. Two methods, fuzzy merging and wall-thickening analysis, were developed for the detection of masses. The fuzzy merging method detected masses with a significant intraluminal component by separating the initial CAD detections of locally cap-like shapes within the colonic wall into mass candidates and polyp candidates. The wall-thickening analysis detected nonintraluminal masses by searching the colonic wall for abnormal thickening. The final regions of the mass candidates were extracted by use of a level set method based on a fast marching algorithm. False-positive (FP) detections were reduced by a quadratic discriminant classifier. The performance of the scheme was evaluated by use of a leave-one-out (round-robin) method with by-patient elimination. All but one of the 14 masses, which was partially cut off from the CTC data set in both supine and prone positions, were detected. The fuzzy merging method detected 11 of the masses, and the wall-thickening analysis detected 3 of the masses including all nonintraluminal masses. In combination, the two methods detected 13 of the 14 masses with 0.21 FPs per patient on average based on the leave-one-out evaluation. Most FPs were generated by extrinsic compression of the colonic wall that would be recognized easily and quickly by a radiologist. The mass detection methods did not affect the result of the polyp detection. The results indicate that the scheme is potentially useful in providing a high-performance CAD scheme for the detection of colorectal neoplasms in CTC.