A novel tool for segmenting 3D medical images based on generalized cylinders and active surfaces

Automatic Aortic Dissection Centerline Extraction Via Morphology-Guided CRN Tracker

Aortic dissection (AD) centerline extraction has important clinical value in the quantitative diagnosis and treatment of AD disease. However, AD centerline extraction is a difficult task and quantitative evaluation is rarely studied. In this work, we propose a fully automatic algorithm to extract AD centerline based on a convolutional regression network (CRN) and the morphological properties of AD. To this end, we first design a topological model to describe the complex topology of AD. With this model, CRNs are trained to estimate the position, tangential vector, and scale of the centerline. The tracking accuracy is further improved by centerline continuity and a gradient-based penalty function. In addition, seed points are extracted on the basis of random regression and line clustering to ensure automated vessel tracking. The proposed method has been evaluated on an AD database and a public aortic database, and achieved high overlapping ratios of 0.9610 and 1.0000, respectively. The tracked centerline is very close to the ground truth and shows good stability, with low average distance errors of 1.4720 mm and 1.8748 mm, respectively.

Accelerating Chan-Vese model with cross-modality guided contrast enhancement for liver segmentation

Accurate and fast liver segmentation remains a challenging and important task for clinicians. Segmentation algorithms are slow and inaccurate due to noise and low quality images in computed tomography (CT) abdominal scans. Chan-Vese is an active contour based powerful and flexible method for image segmentation due to superior noise robustness. However, it is quite slow due to time-consuming partial differential equations, especially for large medical datasets. This can pose a problem for a real-time implementation of liver segmentation and hence, an efficient parallel implementation is highly desirable. Another important aspect is the contrast of CT liver images. Liver slices are sometimes very low in contrast which reduces the overall quality of liver segmentation. Hence, we implement cross-modality guided liver contrast enhancement as a pre-processing step to liver segmentation. GPU implementation of Chan-Vese improves average speedup by 99.811 (± 7.65) times and 14.647 (± 1.155) times with and without enhancement respectively in comparison with the CPU. Average dice, sensitivity and accuracy of liver segmentation are 0.656, 0.816 and 0.822 respectively on the original liver images and 0.877, 0.964 and 0.956 respectively on the enhanced liver images improving the overall quality of liver segmentation.

Fast parallel vessel segmentation

BACKGROUND AND OBJECTIVE

Accurate and fast vessel segmentation from liver slices remain challenging and important tasks for clinicians. The algorithms from the literature are slow and less accurate. We propose fast parallel gradient based seeded region growing for vessel segmentation. Seeded region growing is tedious when the inter connectivity between the elements is unavoidable. Parallelizing region growing algorithms are essential towards achieving real time performance for the overall process of accurate vessel segmentation.

METHODS

The parallel implementation of seeded region growing for vessel segmentation is iterative and hence time consuming process. Seeded region growing is implemented as kernel termination and relaunch on GPU due to its iterative mechanism. The iterative or recursive process in region growing is time consuming due to intermediate memory transfers between CPU and GPU. We propose persistent and grid-stride loop based parallel approach for region growing on GPU. We analyze static region of interest of tiles on GPU for the acceleration of seeded region growing.

RESULTS

We aim fast parallel gradient based seeded region growing for vessel segmentation from CT liver slices. The proposed parallel approach is 1.9x faster compared to the state-of-the-art.

CONCLUSION

We discuss gradient based seeded region growing and its parallel implementation on GPU. The proposed parallel seeded region growing is fast compared to kernel termination and relaunch and accurate in comparison to Chan-Vese and Snake model for vessel segmentation.

GPU acceleration of liver enhancement for tumor segmentation

BACKGROUND AND OBJECTIVE

Medical image segmentation plays a vital role in medical image analysis. There are many algorithms developed for medical image segmentation which are based on edge or region characteristics. These are dependent on the quality of the image. The contrast of a CT or MRI image plays an important role in identifying region of interest i.e. lesion(s). In order to enhance the contrast of image, clinicians generally use manual histogram adjustment technique which is based on 1D histogram specification. This is time consuming and results in poor distribution of pixels over the image. Cross modality based contrast enhancement is 2D histogram specification technique. This is robust and provides a more uniform distribution of pixels over CT image by exploiting the inner structure information from MRI image. This helps in increasing the sensitivity and accuracy of lesion segmentation from enhanced CT image. The sequential implementation of cross modality based contrast enhancement is slow. Hence we propose GPU acceleration of cross modality based contrast enhancement for tumor segmentation.

METHODS

The aim of this study is fast parallel cross modality based contrast enhancement for CT liver images. This includes pairwise 2D histogram, histogram equalization and histogram matching. The sequential implementation of the cross modality based contrast enhancement is computationally expensive and hence time consuming. We propose persistence and grid-stride loop based fast parallel contrast enhancement for CT liver images. We use enhanced CT liver image for the lesion or tumor segmentation. We implement the fast parallel gradient based dynamic seeded region growing for lesion segmentation.

RESULTS

The proposed parallel approach is 104.416 ( ±  5.166) times faster compared to the sequential implementation and increases the sensitivity and specificity of tumor segmentation.

CONCLUSION

The cross modality approach is inspired by 2D histogram specification which incorporates spatial information existing in both guidance and input images for remapping the input image intensity values. The cross modality based liver contrast enhancement improves the quality of tumor segmentation.

Semi-automatic segmentation of femur based on harmonic barrier

BACKGROUND AND OBJECTIVE

Segmentation of the femur from the hip joint in computed tomography (CT) is an important preliminary step in hip surgery planning and simulation. However, this is a time-consuming and challenging task due to the weak boundary, the varying topology of the hip joint, and the extremely narrow or blurred space between the femoral head and the acetabulum. To address these problems, this study proposed a semi-automatic segmentation framework based on harmonic fields for accurate segmentation.

METHODS

The proposed method comprises three steps. First, with high-level information provided by the user, shape information provided by neighboring slices as well as the statistical information in the mask, a region selection method is proposed to effectively locate joint space for the harmonic field. Second, incorporated with an improved gradient, the harmonic field is used to adaptively extract a curve as the barrier that separates the femoral head from the acetabulum accurately. Third, a divide and conquer segmentation strategy based on the harmonic barrier is used to combine the femoral head part and body part as the final segmentation result.

RESULTS

We have tested 40 hips with considerately narrow or disappeared joint spaces. The experimental results are evaluated based on Jaccard, Dice, directional cut discrepancy (DCD) and receiver operating characteristic (ROC), and we achieve the higher Jaccard of 84.02%, Dice of 85.96%, area under curve (AUC) of 89.3%, and the lower error with DCD of 0.52mm. The effective ratio of our method is 79.1% even for cases with severe malformation. The results show that our method performs best in terms of effectiveness and accuracy on the whole data set.

CONCLUSIONS

The proposed method is efficient to segment femurs with narrow joint space. The accurate segmentation results can assist the physicians for osteoarthritis diagnosis in future.

Fully automatic and nonparametric quantification of adipose tissue in fat-water separation MR imaging

Despite increasing demand and research efforts, currently there is no consensus on the protocol for automated and reliable quantification of adipose tissue (AT) and visceral adipose tissue (VAT) using MRI. The purpose of this study was to propose a novel computational method with enhanced objectiveness for the quantification of AT and VAT in fat-water separation MRI. 3T data from IDEAL were acquired for the fat-water separation. Fat tissues were separated from nonfat regions (background air, bone, water, and other nonfat tissues) using K-means clustering (K = 2). From the binary fat mask, arm regions were separated from body based on the relative size of connected component. AT was obtained from the binary body fat mask. With the initial contour as the outer boundary of body fat, the subcutaneous adipose tissue (SAT) and VAT were separated using deformable model driven by a specifically generated deformation field pointing to the inner boundary of SAT. The proposed method was tested on 16 patients with dyslipidemia and evaluated by comparing the correlation with semi-automatic segmentation results. Good robustness was also observed in the proposed method from the Bland-Altman plots. Compared to other established fat segmentation methods, the proposed method is highly objective for fat-water separation MRI with minimal variability induced by subjective parameter settings.

MBIS: multivariate Bayesian image segmentation tool

We present MBIS (Multivariate Bayesian Image Segmentation tool), a clustering tool based on the mixture of multivariate normal distributions model. MBIS supports multichannel bias field correction based on a B-spline model. A second methodological novelty is the inclusion of graph-cuts optimization for the stationary anisotropic hidden Markov random field model. Along with MBIS, we release an evaluation framework that contains three different experiments on multi-site data. We first validate the accuracy of segmentation and the estimated bias field for each channel. MBIS outperforms a widely used segmentation tool in a cross-comparison evaluation. The second experiment demonstrates the robustness of results on atlas-free segmentation of two image sets from scan-rescan protocols on 21 healthy subjects. Multivariate segmentation is more replicable than the monospectral counterpart on T1-weighted images. Finally, we provide a third experiment to illustrate how MBIS can be used in a large-scale study of tissue volume change with increasing age in 584 healthy subjects. This last result is meaningful as multivariate segmentation performs robustly without the need for prior knowledge.

A novel method for pulmonary embolism detection in CTA images

In this paper, we propose a new computer-aided detection (CAD) - based method to detect pulmonary embolism (PE) in computed tomography angiography images (CTAI). Since lung vessel segmentation is the main objective to provide high sensitivity in PE detection, this method performs accurate lung vessel segmentation. To concatenate clogged vessels due to PEs, the starting region of PEs and some reference points (RPs) are determined. These RPs are detected according to the fixed anatomical structures. After lung vessel tree is segmented, the region, intensity, and size of PEs are used to distinguish them. We used the data sets that have heart disease or abnormal tissues because of lung disease except PE in this work. According to the results, 428 of 450 PEs, labeled by the radiologists from 33 patients, have been detected. The sensitivity of the developed system is 95.1% at 14.4 false positive per data set (FP/ds). With this performance, the proposed CAD system is found quite useful to use as a second reader by the radiologists.

A marker-based watershed method for X-ray image segmentation

Digital X-ray images are the most frequent modality for both screening and diagnosis in hospitals. To facilitate subsequent analysis such as quantification and computer aided diagnosis (CAD), it is desirable to exclude image background. A marker-based watershed segmentation method was proposed to segment background of X-ray images. The method consisted of six modules: image preprocessing, gradient computation, marker extraction, watershed segmentation from markers, region merging and background extraction. One hundred clinical direct radiograph X-ray images were used to validate the method. Manual thresholding and multiscale gradient based watershed method were implemented for comparison. The proposed method yielded a dice coefficient of 0.964±0.069, which was better than that of the manual thresholding (0.937±0.119) and that of multiscale gradient based watershed method (0.942±0.098). Special means were adopted to decrease the computational cost, including getting rid of few pixels with highest grayscale via percentile, calculation of gradient magnitude through simple operations, decreasing the number of markers by appropriate thresholding, and merging regions based on simple grayscale statistics. As a result, the processing time was at most 6s even for a 3072×3072 image on a Pentium 4 PC with 2.4GHz CPU (4 cores) and 2G RAM, which was more than one time faster than that of the multiscale gradient based watershed method. The proposed method could be a potential tool for diagnosis and quantification of X-ray images.