Segmentation of the Clustered Cells with Optimized Boundary Detection in Negative Phase Contrast Images - Fig 4

<p><b>(A) Variation of the mean light intensity and (B) its time derivative along detected contours during contour expansion.</b> The mean light intensity gradient decreases with time and mostly reaches its steady state value after about 40 times iteration. The time derivative of the mean light intensity was used as the termination condition for the iteration.</p

Summary of the false segmentation rates obtained with the proposed method.

<p>The values in parenthesis are corresponding percentages.</p><p>Among the three cell lines, MCF 10A cells have the lowest overall false tracking rate.</p><p>^<b>£</b> Oversegmentation: the number of detected cells is more than their actual number in a given area;</p><p>^<b>†</b> Undersegmentation: the number of detected cells is less than their actual number in a given area;</p><p>^<b>‡</b> Over detection: the debris or artifacts present in the field of view are falsely detected as cells.</p><p>Summary of the false segmentation rates obtained with the proposed method.</p

Contour expansion operation with different gains.

<p>(A) When the gains <i>α</i> and <i>β</i> are large, the internal energy dominates the evolution of contours and the contours tends to shrink to minimize the total energy. (B) With decreased gains, the final detected boundaries expanded and were closer to the cell boundaries. (C) When <i>α</i> and <i>β</i> further decreases to 0.01, optimized estimations of cell boundaries with more details are obtained. (D) Mesh plot of the selected area showing cells and converged contours with different combination of<i>α</i> and <i>β</i>. The combination <i>α</i> = <i>β</i> = 1 gives poor estimation of cell boundaries compared with the other two combinations.</p

Comparison of positive (A) and negative (B) phase contrast images of the MCF 10A cells.

<p>In the positive phase contrast image, cells with larger height show reversed image contrast, while in the negative phase contrast image, all cells have consistent image contrast.</p

Segmentation of the clustered cells through peak detection.

<p>(A) Raw phase contrast image with detected peaks for the clustered cells. (B) Mask area preliminarily detected with the thresholding method. (C) Segmentation of the mask area based on the distance between the pixels and the detected peaks inside the mask area. Each pixel is associated with the peak which has the shortest distance with it. In the figure, the subareas for individual cells are plotted as different colours. (D) Boundaries of the subareas were extracted as the initial contours for contour expansion operation. (E) After contour expansion, the final contours for each cell were obtained with improved estimation of cell boundaries.</p

Contour expansion method for cell boundary detection.

<p>(A) It is assumed that cells have illustrated height profile with one peak located above cell nucleus. In negative phase contrast images, the light intensity of cells is proportional to cell height. Therefore, the light intensity distribution over cell surface is similar to height profile of cells with one peak located above each cell body. (B) Mesh plot of the light inteisity for a selected cell marked with yellow arrows in Fig <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0130178#pone.0130178.g002" target="_blank">2A</a> and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0130178#pone.0130178.g002" target="_blank">2B</a>. The light intensity decreases towards the cell boundary, which is consistent with the illustration shown in (A). (C) Quiver plot of the gradient of light intensity for the selected cell. Over cell surface, the gradient of light intensity pointing outwards. (D) The procedure of the contour expansion method for cell boundary detection: (a) The raw negative phase contrast image of the selected cell. (b) The thresholding method was used to get a preliminary mask for the selected cell. (c) The boundary of the mask was extracted and taken as the initial contour. (d) With contour expansion method, the initial contour is driven by the field of gradient of light intensity to gradually converge to the cell boundary. (e) The contour is finally converged at the boundary of the cell, where the contour achieves the minimum energy. The comparison of the contour detected with the thresholding method and the proposed method is shown in (B). It is clear that the proposed method provides an improved estimation of the cell boundary.</p

Demonstration of cell localization, boundary detection, and segmentation of the clustered cells.

<p>(A) For the negative phase contrast image, peaks of light intensity are detected for all cells, as indicated by red circles. (B) Preliminary masks are obtained with the thresholding method. Masks in green are areas with multiple peaks indicating clustered cells. Masks in yellow are areas with single cell. (C) The boundaries of the preliminarily detected masks are extracted to serve as initial contours for individual cells. (D) Contour expansion method is applied to detect cell boundaries for all cells in the-field-of-view. Except one oversegmentation (marked by a yellow arrow) and one falsely detected cell from a debris (marked by a blue arrow), all the other cells are successfully segmented.</p

Comparison of boundary detection obtained with different methods for isolated cells.

<p>(A) The boundaries detected with the thresholding method, region based active contour method, and the contour expansion method for all isolated cells in the field of view. (B) Comparison of the areas enclosed by the contours detected with different methods. The proposed method detects much larger areas than the other two methods.</p

Cell image segmentation result for phase contrast images from four different experiments.

<p>(A) Raw negative phase contrast images. (B) Segmentation result with false segmented cells pointed by arrows.</p

The contour expansion operation is robust to the selected threshold value.

<p>The contour expansion was performed with the initial contours obtained with threshold value of 45 (A) and 65 (B). The final converged contours are close to each other, as shown in (C), which implies that the contour expansion operation is robust with respect to the selection of threshold value. (D) The enlarged image for the selected area in (C).</p