Nuclear Segmentation

Revision as of 19:13, 18 April 2009

This page provides a brief overview of the nuclear segmentation algorithm, usage instruction, parameter description, and sample results.

Overview of the Algorithm

Automatic segmentation of cell nuclei is a widely needed essential step several biological applications. Although much progress has been made, this step continues to require improvements in accuracy, speed, level of automation, and ease of adaptability to new applications.

The segmentation process starts by binarizing the image based on the Graph-Cuts algorithm with automatic learning using minimum error thresholding. In the second step, the nuclear channel image is transformed into another one in which nuclei are represented by blobs of which local maxima points are used to represent seed points of nuclei. The transformation is achieved using a multi-scale Laplacian of Gaussian (LoG) filter, with the addition of automatic and adaptive scale selection. Given a set of seed points and the output of the LoG filter, the next step is to delineate each nucleus. For this, we use a fast clustering technique, called local maximum clustering. In the last step, nuclear contours are refined using a multi-label graph-cuts approach with alpha-expansions. In order to reduce memory requirements and processing time, alpha-expansions is preceded by a graph-coloring step that permits processing a large number of nuclei in parallel. Figure 1 shows a flow-chart illustrating the different steps of our work.

Flow chart outlining the main steps in the nuclear segmentation algorithm

Running the Program

In addition to its use in the GUI, a stand-alone nuclear segmentation executable is available.

Usage:

segment_nuclei.exe <InputImageFileName> <OutputImageFileName> <ParametersFileName>

• The input file is any 2-D or 3-D image with an ITk supported type.
• The output file is a 'tif' image holding the labeled image.
• The parameters file is a text file containing the parameters values as described below.

Description of the Parameters

A sample parameters file:

 ! Segmentation parameters File
 ! All parameters are case sensitive
 high_sensitivity	 :	1
 LoG_size	         :	30
 min_scale	         :	4
 max_scale	         :	6
 xy_clustering_res	 :	2
 z_clustering_res	 :	1
 finalize_segmentation	 :	1
 sampling_ratio_XY_to_Z :	2
 Use_Distance_Map	 :	1
 refinement_range	 :	5

The parameters are as follows: high_sensitivity: a boolean parameter used to enable(1)/disable(0) high sensitivity binarization. Default value is (0). LoG_size: Not used anymore and will be removed min_scale: the minimum scale of the laplacian of gaussian filter. max_scale: the maximum scale of the laplacian of gaussian filter. xy_clustering_res: resolution of the local maximum clustering (in the x-y plane) used in the initial binarization step z_clustering_res: resolution of the local maximum clustering (along z) used in the initial binarization step finalize_segmentation: a boolean parameter sed to enable(1)/disable(0) the segmentation refinement step sampling_ratio_XY_to_Z: ratio of image sampling along Z to the sampling in X-Y Use_Distance_Map: a boolean parameter used to enable(1)/disable(0) the use of the distance constraint to select the LoG scales refinement_range: If segmentation finalization (refinement) is enabled, this parameter sets the maximum distance an initial contour can be shifted

Sample Results

References