A thresholding method for automatic cell image segmentation.

An algorithm for automatic segmentation of PAP-stained cell images and its digital implementation is described. First, the image is filtered in order to eliminate the granularily and small objects in the image which may upset the segmentation procedure. In a second step, information on gradient and compactness is extracted from the filtered image and stored in three histograms as functions of the extinction. From these histograms, two extinction thresholds are computed. These thresholds are suitable to separate the nucleus from the cytoplasm, and the cytoplasm from the background in the filtered image. Masks are determined in this way, and finally used to analyse the nucleus and the cytoplasm in the original image.     

Improving gene quantification by adjustable spot-image restoration

MOTIVATION: One of the major factors that complicate the task of microarray image analysis is that microarray images are distorted by various types of noise. In this study a robust framework is proposed, designed to take into account the effect of noise in microarray images in order to assist the demanding task of microarray image analysis. The proposed framework, incorporates in the microarray image processing pipeline a novel combination of spot adjustable image analysis and processing techniques and consists of the following stages: (1) gridding for facilitating spot identification, (2) clustering (unsupervised discrimination between spot and background pixels) applied to spot image for automatic local noise assessment, (3) modeling of local image restoration process for spot image conditioning (adjustable wiener restoration using an empirically determined degradation function), (4) automatic spot segmentation employing seeded-region-growing, (5) intensity extraction and (6) assessment of the reproducibility (real data) and the validity (simulated data) of the extracted gene expression levels. RESULTS: Both simulated and real microarray images were employed in order to assess the performance of the proposed framework against well-established methods implemented in publicly available software packages (Scanalyze and SPOT). Regarding simulated images, the novel combination of techniques, introduced in the proposed framework, rendered the detection of spot areas and the extraction of spot intensities more accurate. Furthermore, on real images the proposed framework proved of better stability across replicates. Results indicate that the proposed framework improves spots' segmentation and, consequently, quantification of gene expression levels. AVAILABILITY: All algorithms were implemented in Matlab (The Mathworks, Inc., Natick, MA, USA) environment. The codes that implement microarray gridding, adaptive spot restoration and segmentation/intensity extraction are available upon request. Supplementary results and the simulated microarray images used in this study are available for download from: ftp://users:bioinformatics@mipa.med.upatras.gr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

An efficient algorithm for retinal blood vessel segmentation using h-maxima transform and multilevel thresholding

Retinal blood vessel detection and analysis play vital roles in early diagnosis and prevention of several diseases, such as hypertension, diabetes, arteriosclerosis, cardiovascular disease and stroke. This paper presents an automated algorithm for retinal blood vessel segmentation. The proposed algorithm takes advantage of powerful image processing techniques such as contrast enhancement, filtration and thresholding for more efficient segmentation. To evaluate the performance of the proposed algorithm, experiments were conducted on 40 images collected from DRIVE database. The results show that the proposed algorithm yields an accuracy rate of 96.5%, which is higher than the results achieved by other known algorithms.     

An efficient algorithm for retinal blood vessel segmentation using h-maxima transform and multilevel thresholding

Retinal blood vessel detection and analysis play vital roles in early diagnosis and prevention of several diseases, such as hypertension, diabetes, arteriosclerosis, cardiovascular disease and stroke. This paper presents an automated algorithm for retinal blood vessel segmentation. The proposed algorithm takes advantage of powerful image processing techniques such as contrast enhancement, filtration and thresholding for more efficient segmentation. To evaluate the performance of the proposed algorithm, experiments were conducted on 40 images collected from DRIVE database. The results show that the proposed algorithm yields an accuracy rate of 96.5%, which is higher than the results achieved by other known algorithms.     

Segmentation of cervical cell images.

A major problem in the automation of cervical cytology screening is the segmentation of cell images. This paper describes various standard segmentation methods plus one which determines a segmentation threshold based on the stability of the perimeter of the cell as the threshold is varied. As well as contour, certain structural information is used to decide upon the threshold which separates cytoplasm from the background. Once the cytoplasm threshold is found, cytoplasm and nucleus are separated by simple clustering into three groups, cytoplasm, folded cytoplasm and nucleus. These techniques have been tested on 1500 cervical cells that belong to one of eight normal classes and five abnormal classes. A minimum Mahalanobis distance classifier was used to compare results. Manually thresholded cells were classified correctly 66.0% of the time for the 13 class problem and 95.2% of the time on the two (normal-abnormal) class problem. The contour tracing technique was 52.9% and 90.0% correct, respectively.     

High resolution segmentation of cervical cells.

A major problem in the automation of cervical cytology screening is the segmentation of cell images. This paper presents the present status of the work on that problem at the University of Uppsala. A dual resolution system is used. Suspect malignant cells are located at 4 mu resolution. Each such cell is rescanned at 0.5 mu resolution at two different wavelengths, 530 and 570 nm. The nucleus and the cytoplasm are isolated each by two independent methods. For the nucleus adaptive thresholding in the histogram of the 570 nm image and a contouring in a radially transformed version of that image is used. For the cytoplasm a two dimensional thresholding in the 2D histogram and a contouring in a radially transformed version of the 530 nm image is used. If the two nuclear masks agree the surrounding area is checked for disturbing objects. If also the cytoplasm masks agree and are without disturbing objects the whole cell is accepted. The result of the cytoplasm masks agree and are without disturbing objects the whole cell is accepted. The result of the segmentation is thus three categories; free cells, free nuclei and rejected objects. The shape of the objects belonging to the former two categories is checked and irregularly shaped ones are rejected as probably consisting of several overlapping nuclei. Cells passing also this test are classified as normal or malignant. The experience from using this algorithm is discussed and areas for further research are pointed out.     

Automated image detection and segmentation in blood smears.

A simple technique which automatically detects and then segments nucleated cells in Wright's giemsa-stained blood smears is presented. Our method differs from others in 1) the simplicity of our algorithms; 2) inclusion of touching (as well as nontouching) cells; and 3) use of these algorithms to segment as well as to detect nucleated cells employing conventionally prepared smears. Our method involves: 1) acquisition of spectral images; 2) preprocessing the acquired images; 3) detection of single and touching cells in the scene; 4) segmentation of the cells into nuclear and cytoplasmic regions; and 5) postprocessing of the segmented regions. The first two steps of this algorithm are employed to obtain high-quality images, to remove random noise, and to correct aberration and shading effects. Spectral information of the image is used in step 3 to segment the nucleated cells from the rest of the scene. Using the initial cell masks, nucleated cells which are just touching are detected and separated. Simple features are then extracted and conditions applied such that single nucleated cells are finally selected. In step 4, the intensity variations of the cells are then used to segment the nucleus from the cytoplasm. The success rate in segmenting the nucleated cells is between 81 and 93%. The major errors in segmentation of the nucleus and the cytoplasm in the recognized nucleated cells are 3.5% and 2.2%, respectively.     

Segmentation of stained blood cell images measured at high scanning density with high magnification and high numerical aperture optics

In hematological morphology, it is necessary to resolve and analyze the smallest possible cellular details appearing in the light microscope. A prerequisite for computer-aided analysis of subtle morphological features is measuring the cells at a high scanning density with high magnification and high numerical aperture optics. Contrary to visual observations, the information content in a measured picture can be increased by setting the condensor's numerical aperture (NA) greater than the objective's NA. The complexity and heterogeneity of such cell images necessitate a new segmentation method that conserves the morphological information required in the subsequent image analysis, feature extraction, and cell classification. In our segmentation strategy, characteristic color difference thresholds for each nucleus and cytoplasm are combined with geometric operations, probability functions, and a cell model. All thresholds are repeatedly recalculated during the successive improvements of the image masks. None of the thresholds are fixed. This strategy segments blood cell images containing touching cells and large variations in staining, texture, size, and shape. Biological inconsistencies in the calculated cell masks are eliminated by comparing each mask with the cell model criteria integrated into the entire segmentation process. All 20,000 leukocyte images from 120 smears in our leukemia project were segmented with this method.     

Semi-Supervised Segmentation of Ultrasound Images Based on Patch Representation and Continuous Min Cut

Ultrasound segmentation is a challenging problem due to the inherent speckle and some artifacts like shadows, attenuation and signal dropout. Existing methods need to include strong priors like shape priors or analytical intensity models to succeed in the segmentation. However, such priors tend to limit these methods to a specific target or imaging settings, and they are not always applicable to pathological cases. This work introduces a semi-supervised segmentation framework for ultrasound imaging that alleviates the limitation of fully automatic segmentation, that is, it is applicable to any kind of target and imaging settings. Our methodology uses a graph of image patches to represent the ultrasound image and user-assisted initialization with labels, which acts as soft priors. The segmentation problem is formulated as a continuous minimum cut problem and solved with an efficient optimization algorithm. We validate our segmentation framework on clinical ultrasound imaging (prostate, fetus, and tumors of the liver and eye). We obtain high similarity agreement with the ground truth provided by medical expert delineations in all applications (94% DICE values in average) and the proposed algorithm performs favorably with the literature.     

Automatic segmentation of cell nuclei in Feulgen-stained histological sections of prostate cancer and quantitative evaluation of segmentation results

Digital image analysis of cell nuclei is useful to obtain quantitative information for the diagnosis and prognosis of cancer. However, the lack of a reliable automatic nuclear segmentation is a limiting factor for high-throughput nuclear image analysis. We have developed a method for automatic segmentation of nuclei in Feulgen-stained histological sections of prostate cancer. A local adaptive thresholding with an object perimeter gradient verification step detected the nuclei and was combined with an active contour model that featured an optimized initialization and worked within a restricted region to improve convergence of the segmentation of each nucleus. The method was tested on 30 randomly selected image frames from three cases, comparing the results from the automatic algorithm to a manual delineation of 924 nuclei. The automatic method segmented a few more nuclei compared to the manual method, and about 73% of the manually segmented nuclei were also segmented by the automatic method. For each nucleus segmented both manually and automatically, the accuracy (i.e., agreement with manual delineation) was estimated. The mean segmentation sensitivity/specificity were 95%/96%. The results from the automatic method were not significantly different from the ground truth provided by manual segmentation. This opens the possibility for large-scale nuclear analysis based on automatic segmentation of nuclei in Feulgen-stained histological sections.     

超高灵敏度荧光显微成像技术对光敏剂细胞内分布的初步研究

目的：探讨应用基于ICCD的超高灵敏度荧光显微成像系统研究光敏剂细胞内分布的可行性。方法：传代培养内皮细胞、食管癌细胞和肺癌细胞,将不同浓度血卟啉单甲醚(HMME)与细胞共同孵育不同时间。采用荧光显微镜及ICCD组成的荧光显微成像系统采集不同浓度及不同孵育时间条件下HMME的荧光图像,并采用计算机图像处理技术进行图像增强、滤波后计算其细胞浆与细胞核的平均荧光强度比值。同时应用激光共聚焦显微镜图像采集进行对比。结果：HMME浓度为5μg／ml时,荧光显微镜采集到HMME的荧光图像;HMME浓度升高到160μg／ml,激光共聚焦显微镜获得HMME的荧光图像。两组图像的特点都为胞浆中荧光强度较高,细胞核区荧光较弱;细胞浆与细胞核的比值约为2～3：1。结论：荧光显微镜和ICCD采集细胞内光敏剂的荧光图像灵敏度高,方法可靠、实用。HMME较多分布在细胞质中,细胞核吸收较少。     

Segmentation of color images from serous cytology for automated cell classification

OBJECTIVE: To design an automated system for the classification of cells based on analysis of serous cytology, with the aim of segmenting both cytoplasm and nucleus using color information from the images as the main characteristic of the cells. STUDY DESIGN: The segmentation strategy uses color information coupled with mathematical morphology tools, such as watersheds. Cytoplasm and nuclei of all diagnostic cells are retained; erythrocytes and debris are eliminated. Special techniques are used for the separation of clustered cells. RESULTS: A large set of cells was assessed by experts to score the segmentation success rate. All cells were segmented whatever their spatial configurations. The average success rate was 92.5% for nuclei and 91.1% for cytoplasm. CONCLUSION: This color information-based segmentation of images of serous cells is accurate and provides a useful tool. This segmentation strategy will improve the automated classification of cells.     

Computer Simulation of Magnetic Resonance Angiography Imaging: Model Description and Validation

With the development of medical imaging modalities and image processing algorithms, there arises a need for methods of their comprehensive quantitative evaluation. In particular, this concerns the algorithms for vessel tracking and segmentation in magnetic resonance angiography images. The problem can be approached by using synthetic images, where true geometry of vessels is known. This paper presents a framework for computer modeling of MRA imaging and the results of its validation. A new model incorporates blood flow simulation within MR signal computation kernel. The proposed solution is unique, especially with respect to the interface between flow and image formation processes. Furthermore it utilizes the concept of particle tracing. The particles reflect the flow of fluid they are immersed in and they are assigned magnetization vectors with temporal evolution controlled by MR physics. Such an approach ensures flexibility as the designed simulator is able to reconstruct flow profiles of any type. The proposed model is validated in a series of experiments with physical and digital flow phantoms. The synthesized 3D images contain various features (including artifacts) characteristic for the time-of-flight protocol and exhibit remarkable correlation with the data acquired in a real MR scanner. The obtained results support the primary goal of the conducted research, i.e. establishing a reference technique for a quantified validation of MR angiography image processing algorithms.     

A Context-Aware Delayed Agglomeration Framework for Electron Microscopy Segmentation

Electron Microscopy (EM) image (or volume) segmentation has become significantly important in recent years as an instrument for connectomics. This paper proposes a novel agglomerative framework for EM segmentation. In particular, given an over-segmented image or volume, we propose a novel framework for accurately clustering regions of the same neuron. Unlike existing agglomerative methods, the proposed context-aware algorithm divides superpixels (over-segmented regions) of different biological entities into different subsets and agglomerates them separately. In addition, this paper describes a “delayed” scheme for agglomerative clustering that postpones some of the merge decisions, pertaining to newly formed bodies, in order to generate a more confident boundary prediction. We report significant improvements attained by the proposed approach in segmentation accuracy over existing standard methods on 2D and 3D datasets.     

Automatic Region-Based Brain Classification of MRI-T1 Data

Image segmentation of medical images is a challenging problem with several still not totally solved issues, such as noise interference and image artifacts. Region-based and histogram-based segmentation methods have been widely used in image segmentation. Problems arise when we use these methods, such as the selection of a suitable threshold value for the histogram-based method and the over-segmentation followed by the time-consuming merge processing in the region-based algorithm. To provide an efficient approach that not only produce better results, but also maintain low computational complexity, a new region dividing based technique is developed for image segmentation, which combines the advantages of both regions-based and histogram-based methods. The proposed method is applied to the challenging applications: Gray matter (GM), White matter (WM) and cerebro-spinal fluid (CSF) segmentation in brain MR Images. The method is evaluated on both simulated and real data, and compared with other segmentation techniques. The obtained results have demonstrated its improved performance and robustness.     

Automatic Synthesis of Anthropomorphic Pulmonary CT Phantoms

The great density and structural complexity of pulmonary vessels and airways impose limitations on the generation of accurate reference standards, which are critical in training and in the validation of image processing methods for features such as pulmonary vessel segmentation or artery–vein (AV) separations. The design of synthetic computed tomography (CT) images of the lung could overcome these difficulties by providing a database of pseudorealistic cases in a constrained and controlled scenario where each part of the image is differentiated unequivocally. This work demonstrates a complete framework to generate computational anthropomorphic CT phantoms of the human lung automatically. Starting from biological and image-based knowledge about the topology and relationships between structures, the system is able to generate synthetic pulmonary arteries, veins, and airways using iterative growth methods that can be merged into a final simulated lung with realistic features. A dataset of 24 labeled anthropomorphic pulmonary CT phantoms were synthesized with the proposed system. Visual examination and quantitative measurements of intensity distributions, dispersion of structures and relationships between pulmonary air and blood flow systems show good correspondence between real and synthetic lungs (p > 0.05 with low Cohen’s d effect size and AUC values), supporting the potentiality of the tool and the usefulness of the generated phantoms in the biomedical image processing field.     

AXIN-1 protein expression and localization in glioblastoma

The etiology and pathogenesis of tumors of the central nervous system are still inadequately explained. In the present study the expression patterns of a critical molecular component of wnt signaling pathway - axin I was investigated in 42 patients with glioblastoma, the most aggressive form of glial tumors. Immunostaining and image analysis revealed the quantity and localization of the protein. Downregulation of this tumor suppressor expression was observed in 31% of tumors when compared to the levels of axin in healthy brain tissues. Axin was observed in the cytoplasm in 69% of glioblastoma samples, in 21.4% in both the cytoplasm and nucleus and 9.5% had expression solely in the nucleus. Mean values of relative axin's expression obtained by image analysis showed that the highest relative quantity of axin was measured when the protein was in the nucleus and the lowest relative quantity of axin when the protein was localized in the cytoplasm. Investigation on axin's existence at the subcellular level in glioblastomas suggests that axin's expression and spatial regulation is a dynamic process. Despite increasing knowledge on glioma biology and genetics, the prognostic tools for glioblastoma still need improvement. Our findings on expression of axin 1 may contribute to better understanding of glioblastoma molecular profile.