基于CT图像的肺实质分割技术研究

基于CT图像的肺实质分割不仅仅是后续图像处理最基础和最重要的技术,而且是一个典型的亟待解决的问题。基于CT图像的精确的肺实质分割是数学分析、计算机辅助分析和治疗的前提。本文明确了肺实质分割的概念及本质,以及肺实质分割在临床治疗中的重大意义。重点介绍了阈值法、区域生长法、主动轮廓模型和遗传算法,分析了它们的本质以及优缺点,并且列举大量实例说明。最后指出了基于CT图像的肺实质分割方法的发展趋势与面临的挑战。     

Automatic model-based segmentation of the heart in CT images

Automatic image processing methods are a prerequisite to efficiently analyze the large amount of image data produced by computed tomography (CT) scanners during cardiac exams. This paper introduces a model-based approach for the fully automatic segmentation of the whole heart (four chambers, myocardium, and great vessels) from 3-D CT images. Model adaptation is done by progressively increasing the degrees-of-freedom of the allowed deformations. This improves convergence as well as segmentation accuracy. The heart is first localized in the image using a 3-D implementation of the generalized Hough transform. Pose misalignment is corrected by matching the model to the image making use of a global similarity transformation. The complex initialization of the multicompartment mesh is then addressed by assigning an affine transformation to each anatomical region of the model. Finally, a deformable adaptation is performed to accurately match the boundaries of the patient's anatomy. A mean surface-to-surface error of 0.82 mm was measured in a leave-one-out quantitative validation carried out on 28 images. Moreover, the piecewise affine transformation introduced for mesh initialization and adaptation shows better interphase and interpatient shape variability characterization than commonly used principal component analysis.      

Atlas-driven lung lobe segmentation in volumetric X-ray CT images

High-resolution X-ray computed tomography (CT) imaging is routinely used for clinical pulmonary applications. Since lung function varies regionally and because pulmonary disease is usually not uniformly distributed in the lungs, it is useful to study the lungs on a lobe-by-lobe basis. Thus, it is important to segment not only the lungs, but the lobar fissures as well. In this paper, we demonstrate the use of an anatomic pulmonary atlas, encoded with a priori information on the pulmonary anatomy, to automatically segment the oblique lobar fissures. Sixteen volumetric CT scans from 16 subjects are used to construct the pulmonary atlas. A ridgeness measure is applied to the original CT images to enhance the fissure contrast. Fissure detection is accomplished in two stages: an initial fissure search and a final fissure search. A fuzzy reasoning system is used in the fissure search to analyze information from three sources: the image intensity, an anatomic smoothness constraint, and the atlas-based search initialization. Our method has been tested on 22 volumetric thin-slice CT scans from 12 subjects, and the results are compared to manual tracings. Averaged across all 22 data sets, the RMS error between the automatically segmented and manually segmented fissures is 1.96 +/- 0.71 mm and the mean of the similarity indices between the manually defined and computer-defined lobe regions is 0.988. The results indicate a strong agreement between the automatic and manual lobe segmentations.     

3-D segmentation algorithm of small lung nodules in spiral CT images.

Computed tomography (CT) is the most sensitive imaging technique for detecting lung nodules, and is now being evaluated as a screening tool for lung cancer in several large samples studies all over the world. In this report, we describe a semiautomatic method for 3-D segmentation of lung nodules in CT images for subsequent volume assessment. The distinguishing features of our algorithm are the following. 1) The user interaction process. It allows the introduction of the knowledge of the expert in a simple and reproducible manner. 2) The adoption of the geodesic distance in a multithreshold image representation. It allows the definition of a fusion--segregation process based on both gray-level similarity and objects shape. The algorithm was validated on low-dose CT scans of small nodule phantoms (mean diameter 5.3--11 mm) and in vivo lung nodules (mean diameter 5--9.8 mm) detected in the Italung-CT screening program for lung cancer. A further test on small lung nodules of Lung Image Database Consortium (LIDC) first data set was also performed. We observed a RMS error less than 6.6% in phantoms, and the correct outlining of the nodule contour was obtained in 82/95 lung nodules of Italung-CT and in 10/12 lung nodules of LIDC first data set. The achieved results support the use of the proposed algorithm for volume measurements of lung nodules examined with low-dose CT scanning technique.     

Automatic lung segmentation for accurate quantitation of volumetricX-ray CT images

Segmentation of pulmonary X-ray computed tomography (CT) images is a precursor to most pulmonary image analysis applications. This paper presents a fully automatic method for identifying the lungs in three-dimensional (3-D) pulmonary X-ray CT images. The method has three main steps. First, the lung region is extracted from the CT images by gray-level thresholding. Then, the left and right lungs are separated by identifying the anterior and posterior junctions by dynamic programming. Finally, a sequence of morphological operations is used to smooth the irregular boundary along the mediastinum in order to obtain results consistent with those obtained by manual analysis, in which only the most central pulmonary arteries are excluded from the lung region. The method has been tested by processing 3-D CT data sets from eight normal subjects, each imaged three times at biweekly intervals with lungs at 90% vital capacity. We present results by comparing our automatic method to manually traced borders from two image analysts. Averaged over all volumes, the root mean square difference between the computer and human analysis is 0.8 pixels (0.54 mm). The mean intrasubject change in tissue content over the three scans was 2.75% +/- 2.29% (mean +/- standard deviation).     

Precise segmentation of multimodal images.

We propose new techniques for unsupervised segmentation of multimodal grayscale images such that each region-of-interest relates to a single dominant mode of the empirical marginal probability distribution of grey levels. We follow the most conventional approaches in that initial images and desired maps of regions are described by a joint Markov-Gibbs random field (MGRF) model of independent image signals and interdependent region labels. However, our focus is on more accurate model identification. To better specify region borders, each empirical distribution of image signals is precisely approximated by a linear combination of Gaussians (LCG) with positive and negative components. We modify an expectation-maximization (EM) algorithm to deal with the LCGs and also propose a novel EM-based sequential technique to get a close initial LCG approximation with which the modified EM algorithm should start. The proposed technique identifies individual LCG models in a mixed empirical distribution, including the number of positive and negative Gaussians. Initial segmentation based on the LCG models is then iteratively refined by using the MGRF with analytically estimated potentials. The convergence of the overall segmentation algorithm at each stage is discussed. Experiments show that the developed techniques segment different types of complex multimodal medical images more accurately than other known algorithms.     

基于改进C-V模型的肾脏CT图像分割方法

生物组织的自动分割是计算机辅助诊断和病变 检测的关键步骤。在腹腔CT图像中,肾脏组织本身的灰度不均匀性使得传统C-T模型无法准 确实现肾脏的分割。为了解决上述问题,本文结合图像全局和 局部统计信息改进了传统的C-V模型。 基于先验知识,提出了描述肾脏组织皮质 特征的数学表达式。选择感兴趣区域,在预处理阶段获得了CT图像中肾脏的大致初始轮廓。 随后,应用C-V模型进行轮廓演化时引入局域信息,提高了C-V模型的局部适应性。实验结 果表明,与现有方法相比,本文的方法的结果更接近于人工分割结果,其肾脏分割结果的Di ce系数平均值为94.0%。     

Automated segmentation of acetabulum and femoral head from 3-d CT images

This paper describes several new methods and software for automatic segmentation of the pelvis and the femur, based on clinically obtained multislice computed tomography (CT) data. The hip joint is composed of the acetabulum, cavity of the pelvic bone, and the femoral head. In vivo CT data sets of 60 actual patients were used in the study. The 120 (60 /spl times/ 2) hip joints in the data sets were divided into four groups according to several key features for segmentation. Conventional techniques for classification of bony tissues were first employed to distinguish the pelvis and the femur from other CT tissue images in the hip joint. Automatic techniques were developed to extract the boundary between the acetabulum and the femoral head. An automatic method was built up to manage the segmentation task according to image intensity of bone tissues, size, center, shape of the femoral heads, and other characters. The processing scheme consisted of the following five steps: 1) preprocessing, including resampling 3-D CT data by a modified Sine interpolation to create isotropic volume and to avoid Gibbs ringing, and smoothing the resulting images by a 3-D Gaussian filter; 2) detecting bone tissues from CT images by conventional techniques including histogram-based thresholding and binary morphological operations; 3) estimating initial boundary of the femoral head and the joint space between the acetabulum and the femoral head by a new approach utilizing the constraints of the greater trochanter and the shapes of the femoral head; 4) enhancing the joint space by a Hessian filter; and 5) refining the rough boundary obtained in step 3) by a moving disk technique and the filtered images obtained in step 4). The above method was implemented in a Microsoft Windows software package and the resulting software is freely available on the Internet. The feasibility of this method was tested on the data sets of 60 clinical cases (5000 CT images).     

Real-time volume rendering visualization of dual-modality PET/CT images with interactive fuzzy thresholding segmentation.

Three-dimensional (3-D) visualization has become an essential part for imaging applications, including image-guided surgery, radiotherapy planning, and computer-aided diagnosis. In the visualization of dual-modality positron emission tomography and computed tomography (PET/CT), 3-D volume rendering is often limited to rendering of a single image volume and by high computational demand. Furthermore, incorporation of segmentation in volume rendering is usually restricted to visualizing the presegmented volumes of interest. In this paper, we investigated the integration of interactive segmentation into real-time volume rendering of dual-modality PET/CT images. We present and validate a fuzzy thresholding segmentation technique based on fuzzy cluster analysis, which allows interactive and real-time optimization of the segmentation results. This technique is then incorporated into a real-time multi-volume rendering of PET/CT images. Our method allows a real-time fusion and interchangeability of segmentation volume with PET or CT volumes, as well as the usual fusion of PET/CT volumes. Volume manipulations such as window level adjustments and lookup table can be applied to individual volumes, which are then fused together in real time as adjustments are made. We demonstrate the benefit of our method in integrating segmentation with volume rendering in its application to PET/CT images. Responsive frame rates are achieved by utilizing a texture-based volume rendering algorithm and the rapid transfer capability of the high-memory bandwidth available in low-cost graphic hardware.     

Tumor segmentation in lung CT images based on support vector machine and improved level set

In lung CT images, the edge of a tumor is frequently fuzzy because of the complex relationship between tumors and tissues, especially in cases that the tumor adheres to the chest and lung in the pathology area. This makes the tumor segmentation more difficult. In order to segment tumors in lung CT images accurately, a method based on support vector machine (SVM) and improved level set model is proposed. Firstly, the image is divided into several block units; then the texture, gray and shape features of each block are extracted to construct eigenvector and then the SVM classifier is trained to detect suspicious lung lesion areas; finally, the suspicious edge is extracted as the initial contour after optimizing lesion areas, and the complete tumor segmentation can be obtained by level set model modified with morphological gradient. Experimental results show that this method can efficiently and fast segment the tumors from complex lung CT images with higher accuracy.     

Three-dimensional segmentation and growth-rate estimation of small pulmonary nodules in helical CT images

Small pulmonary nodules are a common radiographic finding that presents an important diagnostic challenge in contemporary medicine. While pulmonary nodules are the major radiographic indicator of lung cancer, they may also be signs of a variety of benign conditions. Measurement of nodule growth rate over time has been shown to be the most promising tool in distinguishing malignant from nonmalignant pulmonary nodules. In this paper, we describe three-dimensional (3-D) methods for the segmentation, analysis, and characterization of small pulmonary nodules imaged using computed tomography (CT). Methods for the isotropic resampling of anisotropic CT data are discussed. 3-D intensity and morphology-based segmentation algorithms are discussed for several classes of nodules. New models and methods for volumetric growth characterization based on longitudinal CT studies are developed. The results of segmentation and growth characterization methods based on in vivo studies are described. The methods presented are promising in their ability to distinguish malignant from nonmalignant pulmonary nodules and represent the first such system in clinical use.     

Automated segmentation of the lumbar pedicle in CT images for spinal fusion surgery

Exact information about the shape of a lumbar pedicle can increase operation accuracy and safety during computer-aided spinal fusion surgery, which requires extreme caution on the part of the surgeon, due to the complexity and delicacy of the procedure. In this paper, a robust framework for segmenting the lumbar pedicle in computed tomography (CT) images is presented. The framework that has been designed takes a CT image, which includes the lumbar pedicle as input, and provides the segmented lumbar pedicle in the form of 3-D voxel sets. This multistep approach begins with 2-D dynamic thresholding using local optimal thresholds, followed by procedures to recover the spine geometry in a high curvature environment. A subsequent canal reference determination using proposed thinning-based integrated cost is then performed. Based on the obtained segmented vertebra and canal reference, the edge of the spinal pedicle is segmented. This framework has been tested on 84 lumbar vertebrae of 19 patients requiring spinal fusion. It was successfully applied, resulting in an average success rate of 93.22 % and a final mean error of 0.14 ± 0.05 mm. Precision errors were smaller than 1 % for spine pedicle volumes. Intra- and interoperator precision errors were not significantly different.     

A new model-based framework for lung tissue segmentation in three-dimensional thoracic CT images

A new framework for model-based lung tissue segmentation in three-dimensional thoracic CT images is proposed. In the first stage, a parametric model for lung segmenting surface is created using shape representation based on level sets method. This model is constituted by the sum of a mean distance function and a number of weighted eigenshapes. Consequently, unlike the other model-based segmentation methods, there is no need to specify any marker point in this model. In the second stage, the segmenting surface is varied so as to be matched with the binarized input image. For this purpose, a region-based energy function is minimized with respect to the parameters including the weights of eigenshapes and coefficients of a three-dimensional similarity transform. Finally, the resulted segmenting surface is post-processed in order to improve its fitness with the lung borders of the input image. The experimental results demonstrated the outperformance of the proposed framework over its model-based counterparts in model matching stage. Moreover, it performed slightly better in terms of final segmentation results.     

Automatic segmentation of high-throughput RNAi fluorescent cellular images.

High-throughput genome-wide RNA interference (RNAi) screening is emerging as an essential tool to assist biologists in understanding complex cellular processes. The large number of images produced in each study make manual analysis intractable; hence, automatic cellular image analysis becomes an urgent need, where segmentation is the first and one of the most important steps. In this paper, a fully automatic method for segmentation of cells from genome-wide RNAi screening images is proposed. Nuclei are first extracted from the DNA channel by using a modified watershed algorithm. Cells are then extracted by modeling the interaction between them as well as combining both gradient and region information in the Actin and Rac channels. A new energy functional is formulated based on a novel interaction model for segmenting tightly clustered cells with significant intensity variance and specific phenotypes. The energy functional is minimized by using a multiphase level set method, which leads to a highly effective cell segmentation method. Promising experimental results demonstrate that automatic segmentation of high-throughput genome-wide multichannel screening can be achieved by using the proposed method, which may also be extended to other multichannel image segmentation problems.     

Knowledge-based segmentation of Landsat images

A knowledge-based approach for Landsat image segmentation is proposed. The image segmentation problem is solved by extracting kernel information from the input image to provide an initial interpretation of the image and by using a knowledge-based hierarchical classifier to discriminate between major land-cover types in the study area. The proposed method is designed in such a way that a Landsat image can be segmented and interpreted without any prior image-dependent information. The general spectral land-cover knowledge is constructed from the training land-cover data, and the road information of an image is obtained through a road-detection program     

Intensity-based segmentation of microarray images

The underlying principle in microarray image analysis is that the spot intensity is a measure of the gene expression. This implicitly assumes the gene expression of a spot to be governed entirely by the distribution of the pixel intensities. Thus, a segmentation technique based on the distribution of the pixel intensities is appropriate for the current problem. In this paper, clustering-based segmentation is described to extract the target intensity of the spots. The approximate boundaries of the spots in the microarray are determined by manual adjustment of rectilinear grids. The distribution of the pixel intensity in a grid containing a spot is assumed to be the superposition of the foreground and the local background. The k-means clustering technique and the partitioning around medoids (PAM) were used to generate a binary partition of the pixel intensity distribution. The median (k-means) and the medoid (PAM) of the cluster members are chosen as the cluster representatives. The effectiveness of the clustering-based segmentation techniques was tested on publicly available arrays generated in a lipid metabolism experiment (Callow et al., 2000). The results are compared against those obtained using the region-growing approach (SPOT) (Yang et al., 2001). The effect of additive white Gaussian noise is also investigated.     

Optimal segmentation of cell images

An optimal segmentation algorithm for light microscopic cell images is presented. The image segmentation is performed by thresholding a parametric image approximating the original image. Using the mean squared error between the original and the constructed image as the cost function, the segmentation problem is transformed into an optimisation process where parametric parameters are determined that minimise the defined cost function. The cost function is iteratively minimised using an unsupervised learning rule to adjust the parameters, and a parametric image is constructed at each iteration, based on the obtained parameters. The cell region is extracted by thresholding the final parametric image, where the threshold is one of the image parameters. Application results to real cervical images are provided to show the performance of the proposed segmentation approach. Experimental segmentation results are presented for the proposed optimal algorithm for synthetic cell images corrupted by variant levels of noise; these results are compared with the K-means clustering method and Bayes classifier in terms of classification errors     

Computer-aided characterization of mammographic masses: accuracy of mass segmentation and its effects on characterization.

Mass segmentation is used as the first step in many computer-aided diagnosis (CAD) systems for classification of breast masses as malignant or benign. The goal of this paper was to study the accuracy of an automated mass segmentation method developed in our laboratory, and to investigate the effect of the segmentation stage on the overall classification accuracy. The automated segmentation method was quantitatively compared with manual segmentation by two expert radiologists (R1 and R2) using three similarity or distance measures on a data set of 100 masses. The area overlap measures between R1 and R2, the computer and R1, and the computer and R2 were 0.76 +/- 0.13, 0.74 +/- 0.11, and 0.74 +/- 0.13, respectively. The interobserver difference in these measures between the two radiologists was compared with the corresponding differences between the computer and the radiologists. Using three similarity measures and data from two radiologists, a total of six statistical tests were performed. The difference between the computer and the radiologist segmentation was significantly larger than the interobserver variability in only one test. Two sets of texture, morphological, and spiculation features, one based on the computer segmentation, and the other based on radiologist segmentation, were extracted from a data set of 249 films from 102 patients. A classifier based on stepwise feature selection and linear discriminant analysis was trained and tested using the two feature sets. The leave-one-case-out method was used for data sampling. For case-based classification, the area Az under the receiver operating characteristic (ROC) curve was 0.89 and 0.88 for the feature sets based on the radiologist segmentation and computer segmentation, respectively. The difference between the two ROC curves was not statistically significant.     

Multiscale morphological segmentation of gray-scale images

In this paper, the authors have proposed a method of segmenting gray level images using multiscale morphology. The approach resembles the watershed algorithm in the sense that the dark (respectively bright) features which are basically canyons (respectively mountains) on the surface topography of the gray level image are gradually filled (respectively clipped) using multiscale morphological closing (respectively opening) by reconstruction with isotropic structuring element. The algorithm detects valid segments at each scale using three criteria namely growing, merging and saturation. Segments extracted at various scales are integrated in the final result. The algorithm is composed of two passes preceded by a preprocessing step for simplifying small scale details of the image that might cause over-segmentation. In the first pass feature images at various scales are extracted and kept in respective level of morphological towers. In the second pass, potential features contributing to the formation of segments at various scales are detected. Finally the algorithm traces the contours of all such contributing features at various scales. The scheme after its implementation is executed on a set of test images (synthetic as well as real) and the results are compared with those of few other standard methods. A quantitative measure of performance is also formulated for comparing the methods.     

Adaptive river segmentation in SAR images

There is difficulty for distinguishing of river and shadow in Synthetic Aperture Radar (SAR) images. A method of river segmentation in SAR images based on wavelet energy and gradient is proposed in this paper. It mainly includes two algorithms: coarse segmentation and refined segmentation. Firstly, The river regions are coarsely segmented by the wavelet energy feature,and then refined segmented accurately by the gradient threshold which is got adaptively. The experimental results show the validity of the method, which provides a good foundation for targets detection above the river.