Volume reconstruction of large tissue specimens from serial physical sections using confocal microscopy and correction of cutting deformations by elastic registration

A set of methods leading to volume reconstruction of biological specimens larger than the field of view of a confocal laser scanning microscope (CLSM) is presented. Large tissue specimens are cut into thin physical slices and volume data sets are captured from all studied physical slices by CLSM. Overlapping spatial tiles of the same physical slice are stitched in horizontal direction. Image volumes of successive physical slices are linked in axial direction by applying an elastic registration algorithm to compensate for deformations because of cutting the specimen. We present a method enabling us to keep true object morphology using a priori information about the shape and size of the specimen, available from images of the cutting planes captured by a USB light microscope immediately before cutting the specimen by a microtome. The errors introduced by elastic registration are evaluated using a stereological point counting method and the Procrustes distance. Finally, the images are enhanced to compensate for the effect of the light attenuation with depth and visualized by a hardware accelerated volume rendering. Algorithmic steps of the reconstruction, namely elastic registration, object morphology preservation, image enhancement, and volume visualization, are implemented in a new Rapid3D software package. Because confocal microscopes get more and more frequently used in scientific laboratories, the described volume reconstruction may become an easy‐to‐apply tool to study large biological objects, tissues, and organs in histology, embryology, evolution biology, and developmental biology. In this work, we demonstrate the reconstruction using a postcranial part of a 17‐day‐old laboratory Wistar rat embryo. Microsc. Res. Tech., 2009. © 2008 Wiley‐Liss, Inc.     

Algorithms for accurate 3D registration of neuronal images acquired by confocal scanning laser microscopy

This paper presents automated and accurate algorithms based on high‐order transformation models for registering three‐dimensional (3D) confocal images of dye‐injected neurons. The algorithms improve upon prior methods in several ways, and meet the more stringent image registration needs of applications such as two‐view attenuation correction recently developed by us. First, they achieve high accuracy (≈ 1.2 voxels, equivalent to 0.4 µm) by using landmarks, rather than intensity correlations, and by using a high‐dimensional affine and quadratic transformation model that accounts for 3D translation, rotation, non‐isotropic scaling, modest curvature of field, distortions and mechanical inconsistencies introduced by the imaging system. Second, they use a hierarchy of models and iterative algorithms to eliminate potential instabilities. Third, they incorporate robust statistical methods to achieve accurate registration in the face of inaccurate and missing landmarks. Fourth, they are fully automated, even estimating the initial registration from the extracted landmarks. Finally, they are computationally efficient, taking less than a minute on a 900‐MHz Pentium III computer for registering two images roughly 70 MB in size. The registration errors represent a combination of modelling, estimation, discretization and neuron tracing errors. Accurate 3D montaging is described; the algorithms have broader applicability to images of vasculature, and other structures with distinctive point, line and surface landmarks.     

Attenuation correction in confocal laser microscopes: a novel two-view approach

Confocal microscopy is a three‐dimensional (3D) imaging modality, but the specimen thickness that can be imaged is limited by depth‐dependent signal attenuation. Both software and hardware methods have been used to correct the attenuation in reconstructed images, but previous methods do not increase the image signal‐to‐noise ratio (SNR) using conventional specimen preparation and imaging. We present a practical two‐view method that increases the overall imaging depth, corrects signal attenuation and improves the SNR. This is achieved by a combination of slightly modified but conventional specimen preparation, image registration, montage synthesis and signal reconstruction methods. The specimen is mounted in a symmetrical manner between a pair of cover slips, rather than between a slide and a cover slip. It is imaged sequentially from both sides to generate two 3D image stacks from perspectives separated by approximately 180° with respect to the optical axis. An automated image registration algorithm performs a precise 3D alignment, and a model‐based minimum mean squared algorithm synthesizes a montage, combining the content of both the 3D views. Experiments with images of individual neurones contrasted with a space‐filling fluorescent dye in thick brain tissue slices produced precise 3D montages that are corrected for depth‐dependent signal attenuation. The SNR of the reconstructed image is maximized by the method, and it is significantly higher than in the single views after applying our attenuation model. We also compare our method with simpler two‐view reconstruction methods and quantify the SNR improvement. The reconstructed images are a more faithful qualitative visualization of the specimen's structure and are quantitatively more accurate, providing a more rigorous basis for automated image analysis.     

Three‐dimensional imaging of porous media using confocal laser scanning microscopy

In the last decade, imaging techniques capable of reconstructing three‐dimensional (3‐D) pore‐scale model have played a pivotal role in the study of fluid flow through complex porous media. In this study, we present advances in the application of confocal laser scanning microscopy (CLSM) to image, reconstruct and characterize complex porous geological materials with hydrocarbon reservoir and CO₂ storage potential. CLSM has a unique capability of producing 3‐D thin optical sections of a material, with a wide field of view and submicron resolution in the lateral and axial planes. However, CLSM is limited in the depth (z‐dimension) that can be imaged in porous materials. In this study, we introduce a ‘grind and slice’ technique to overcome this limitation. We discuss the practical and technical aspects of the confocal imaging technique with application to complex rock samples including Mt. Gambier and Ketton carbonates. We then describe the complete workflow of image processing to filtering and segmenting the raw 3‐D confocal volumetric data into pores and grains. Finally, we use the resulting 3‐D pore‐scale binarized confocal data obtained to quantitatively determine petrophysical pore‐scale properties such as total porosity, macro‐ and microporosity and single‐phase permeability using lattice Boltzmann (LB) simulations, validated by experiments.     

Intensity correction of fluorescent confocal laser scanning microscope images by mean-weight filtering

This paper addresses the problem of intensity correction of fluorescent confocal laser scanning microscope images. Confocal laser scanning microscope images are frequently used in medicine for obtaining 3D information about specimen structures by imaging a set of 2D cross sections and performing 3D volume reconstruction afterwards. However, 2D images acquired from fluorescent confocal laser scanning microscope images demonstrate significant intensity heterogeneity, for example, due to photo‐bleaching and fluorescent attenuation in depth. We developed an intensity heterogeneity correction technique that (a) adjusts the intensity heterogeneity of 2D images, (b) preserves fine structural details and (c) enhances image contrast, by performing spatially adaptive mean‐weight filtering. Our solution is obtained by formulating an optimization problem, followed by filter design and automated selection of filtering parameters. The proposed filtering method is experimentally compared with several existing techniques by using four quality metrics: contrast, intensity heterogeneity (entropy) in a low frequency domain, intensity distortion in a high frequency domain and saturation. Based on our experiments and the four quality metrics, the developed mean‐weight filtering outperforms other intensity correction methods by at least a factor of 1.5 when applied to fluorescent confocal laser scanning microscope images.     

The distribution of caveolin-3 immunofluorescence in skeletal muscle fibre membrane defined by dual channel confocal laser scanning microscopy, fast Fourier transform and image modelling

Membrane domains rich in caveolin‐3 overlie sarcomeric actin in skeletal muscle. The membrane exhibits a regular array of caveolin‐3 immunofluorescence using confocal laser scanning microscopy (CLSM). Fourier analysis of tissue imaged by CLSM accurately defines a repeating intensity with a long‐axis spacing of 1.48 µm confirmed by measurement of direct images. Reverse fast Fourier transform (FFT) and image‐modelling allow reconstruction of the pattern. Mathematical modelling has allowed replication of several features of the FFT, including the second order maxima that confirm the relatively high information content of the original images. Measurements of membrane‐pattern primary long‐axis spacings are consistent with our measurements of the I‐band sarcomere repeat in similarly prepared specimens labelled with fluorescent phalloidin or imaged using differential interference contrast microscopy. Dual‐channel CLSM analysis of the sarcomeric banding pattern of actin and the repeating pattern of muscle fibre membrane caveolin showed that caveolae overlie the I‐band. The anti‐caveolin immunofluorescence is deficient over the Z‐disc and maximal toward each of the I‐band extremities. A mechanism of membrane shape change in which membrane–lipid molecules are interposed between more stable anchored rafts associated with caveolae can be envisaged. Thus, increasing girth and reducing length of the sarcolemma in rapid contraction may be explained.     

Compensation of inhomogeneous fluorescence signal distribution in 2D images acquired by confocal microscopy

In images acquired by confocal laser scanning microscopy (CLSM), regions corresponding to the same concentration of fluorophores in the specimen should be mapped to the same grayscale levels. However, in practice, due to multiple distortion effects, CLSM images of even homogeneous specimen regions suffer from irregular brightness variations, e.g., darkening of image edges and lightening of the center. The effects are yet more pronounced in images of real biological specimens. A spatially varying grayscale map complicates image postprocessing, e.g., in alignment of overlapping regions of two images and in 3D reconstructions, since measures of similarity usually assume a spatially independent grayscale map. We present a fast correction method based on estimating a spatially variable illumination gain, and multiplying acquired CLSM images by the inverse of the estimated gain. The method does not require any special calibration of reference images since the gain estimate is extracted from the CLSM image being corrected itself. The proposed approach exploits two types of morphological filters: the median filter and the upper Lipschitz cover. The presented correction method, tested on images of both artificial (homogeneous fluorescent layer) and real biological specimens, namely sections of a rat embryo and a rat brain, proved to be very fast and yielded a significant visual improvement. Microsc. Res. Tech., 2011. © 2010 Wiley‐Liss, Inc.     

MRT letter: A unified accelerated maximum likelihood technique for widefield,confocal, and super‐resolution 4Pi microscopy

We propose an algorithmic technique for accelerating maximum likelihood (ML) algorithm for image reconstruction in fluorescence microscopy. This is made possible by integrating Biggs–Andrews (BA) method with ML approach. The results on widefield, confocal, and super‐resolution 4Pi microscopy reveal substantial improvement in the speed of 3D image reconstruction (the number of iterations has reduced by approximately one‐half). Moreover, the quality of reconstruction obtained using accelerated ML closely resembles with nonaccelerated ML method. The proposed technique is a step closer to realize real‐time reconstruction in 3D fluorescence microscopy. Microsc. Res. Tech. 78:331–335, 2015. © 2015 Wiley Periodicals, Inc.     

Four-dimensional factor analysis of confocal image sequences (4D-FAMIS) to detect and characterize low copy numbers of human papillomavirus DNA by FISH in HeLa and SiHa cells

Visualization and localization of specific DNA sequences were performed by fluorescence in situ hybridization, confocal laser scanning microscopy (CLSM), and four-dimensional factor analysis of biomedical image sequences (4D-FAMIS). HeLa and SiHa cells containing, respectively 20–50 and 1–2 copies per cell of human papillomavirus (HPV) DNA type 18 and 16 integrated in cellular DNA were used as models. HPV-DNA was identified using DNA probes containing the whole genome of HPV-DNA type 18 or 16, and DNA–DNA hybrids were revealed by alkaline phosphatase and Fast Red. Cell nuclei were counterstained with thiazole orange (TO) or TOTO-iodide. 4D image sequences were obtained using successive dynamic or spectral sequences of images on different optical sections from CLSM. The location of fluorescent signals within the preparations was determined by FAMIS. This original method summarizes image sequences into a reduced number of images called factor images, and curves called factors. Factors estimate different individual physical behaviours in the sequence such as extinction velocity, spectral patterns and depth emission profiles. Factor images correspond to spatial distributions of the different factors. We distinguished between Fast Red and nucleus stainings in HPV-DNA hybridization signals by taking into account differences in their extinction velocities (fluorescence decay rate) or spectral patterns, and in their focus (depth emission profiles). In HeLa cells, factor images showed that Fast-Red-stained targets could be distinguished from nucleus stainings, and were located on different focal planes of the nuclei. In SiHa cells, 4D-FAMIS determined as few as 1–2 copies per cell of HPV-DNA type 16 located in continuous focal planes. Therefore, 4D-FAMIS, together with CLSM, made the detection and characterization of low copy numbers of genes in whole cells possible.     

Metrological challenges for reconstruction of 3‐D microstructures by focused ion beam tomography methods

The development of combined focused ion beam and scanning electron microscopes has enabled significant advances in the characterization of the 3‐D structure of materials. The repeated removal of thin layers or slices with an ion beam and imaging or mapping the chemical or crystallographic structure of each slice enables a 3‐D reconstruction from the images or maps. The accuracy of the reconstruction thus depends on the accuracy with which the slice thickness is measured and maintained throughout the process, and the alignment accuracy of the slices achieved during acquisition or by postacquisition corrections. A survey of papers published in this field suggests that the reconstruction accuracy is not often considered or reported. Using examples from examination of the 3‐D structure of hardmetals, issues affecting the accuracy of slice thicknesses and image realignments are examined and illustrated and potential errors quantified by the use of fiducial markers and the expected isotropy of the hardmetal structure itself.     

Confocal microscopic image sequence compression using vector quantization and three-dimensional pyramid

The three-dimensional (3-D) pyramid compressor project at the University of Glasgow has developed a compressor for images obtained from confocal laser scanning microscopy (CLSM) device. The proposed method using a combination of image pyramid coder and vector quantization techniques has good performance at compressing confocal volume image data. An experiment was conducted on several kinds of CLSM data using the presented compressor compared with other well-known volume data compressors, such as MPEG-1. The results showed that the 3-D pyramid compressor gave a higher subjective and objective image quality of reconstructed images at the same compression ratio, and presented more acceptable results when applying image processing filters on reconstructed images.     

Confocal laser scanning microscopy for the study of the morphological changes of the postextraction sites

A better understanding of the remodeling process of postextraction sockets is essential in dental treatment planning. The aim of this study was to evaluate whether confocal laser scanning microscopy (CLSM) can be applied to imaging contour changes of postextraction sites, as well as to its quantification with image analysis of obtained three‐dimensional images. This work describes a new application of the CLSM technique. The system used was the OLS3100‐USS, LEXT model (Olympus®). CLSM was used for the surface analysis of the extraction site. The measurements taken with CLSM were: (1) mesio‐distal distance, (2) alveolar ridge thickness, and (3) vestibular and lingual alveolar ridge height. Results of study cast scanning at baseline, 1 and 3 months after tooth extraction, with CLSM are well‐detailed images of postextraction areas. The CLSM technique used in study casts is a valid method to measure the dimensional changes that happen in the edentulous area after tooth extraction. This technique allows the evaluation of changes in mesio‐distal distance, thickness of the alveolar ridge and alveolar ridge height based on the measurements on the alveolar contours. Microsc. Res. Tech., 2011. © 2011 Wiley‐Liss, Inc.     

Registration of serial sections of mouse liver cell nuclei

Image registration of biological tissue is essential for 3D reconstruction, which is important for visualizing and quantifying the 3D relationships between internal structures of an object. The biological role of DNA organization, which is an extremely complex 3D architecture within the cell nucleus, has come into focus since it has become clear that the chromatin structure in itself functions as a regulator of DNA. Thus, 3D reconstruction of cell nuclei based on consecutive series of high-resolution ultrathin slices may provide new information about the chromatin structure and its organizational changes during carcinogenesis. This work focuses mainly on the problem of registering successive serial transmission electron micrographs of ultrathin sections of mouse liver cell nuclei to analyse the 3D chromatin structure. A five-step semiautomatic interactive registration method is proposed. The first two steps of the procedure correct the rotation and translation components by using the phase correlation. The third, fourth and fifth steps correct the global distortion, employing a point mapping method based on different ways of selecting the control points. In step three, the control points were automatically computed by phase correlating corresponding subimages of the reference and sensed image. A semiautomatic method is used in the fourth step to select the control points, i.e. an automated method for computing the centre of mass of manually identified anatomical structures in neighbouring slices. For the sections which could not be properly corrected by the four steps, a final step is introduced, where control points are manually selected in the reference and sensed images. An algorithm is proposed to examine the spatial distribution of selected control points. Four sets of serial sections of mouse liver cell nuclei, each with approximately 100 sections, are registered by the proposed method and also registered manually for the comparison of registration accuracy. Artificial X–Z and Z–Y sections of registered series were visually compared for the smoothness of the nuclear membrane. To quantify the registration accuracy and the extent of registration, the correlation coefficient (C) and the overlap index (C₀) were computed over the registered structure of interest. In addition to the visual comparison and the comparison of C and C₀, the registered serial sets were compared by 3D GLCM-based texture features in the Z direction. The results demonstrate that the proposed semiautomatic registration technique achieved accurate results comparable to the manual registration. The proposed registration method relies only on the operator for rough pinpointing of cellular structures. Therefore, it should provide better reproducibility, and allow the user to operate the system faster and in a more relaxed manner than in a manual registration.     

A simple method allowing DIC imaging in conjunction with confocal microscopy

Current optical methods to collect Nomarski differential interference contrast (DIC) or phase images with a transmitted light detector (TLD) in conjunction with confocal laser scanning microscopy (CLSM) can be technically challenging and inefficient. We describe for the first time a simple method that combines the use of the commercial product QPm (Iatia, Melbourne Australia) with brightfield images collected with the TLD of a CLSM, generating DIC, phase, Zernike phase, dark-field or Hoffman modulation contrast images. The brightfield images may be collected at the same time as the confocal images. This method also allows the calculation of contrast-enhanced images from archival data. The technique described here allows for the creation of contrast-enhanced images such as DIC or phase, without compromising the intensity or quality of confocal images collected simultaneously. Provided the confocal microscope is equipped with a motorized z-drive and a TLD, no hardware or optical modifications are required. The contrast-enhanced images are calculated with software using the quantitative phase-amplitude microscopy technique ( Barone-Nugent et al., 2002 ). This technique, being far simpler during image collection, allows the microscopist to concentrate on their confocal imaging and experimental procedures. Unlike conventional DIC, this technique may be used to calculate DIC images when cells are imaged through plastic, and without the use of expensive strain-free objective lenses.     

A novel method for three-dimensional observation of the vascular networks in the whole mouse brain

A novel method for acquiring serial images suitable for three-dimensional reconstruction of vascular networks in the whole brain of mouse was developed. The brain infused with a White India ink-gelatin solution was fixed and embedded in paraffin containing Sudan Black B through xylene also containing Sudan Black B. Each sliced surface of the paraffin block was coated with liquid paraffin and its image was serially acquired. Coating with liquid paraffin extremely improved the quality of the image. The series of serial images was free of distortion and a three-dimensional image was reconstructed without the problem of the alignment and registration of adjacent images. The volume-rendered image indicated three-dimensional distribution of blood vessels in a whole brain. No ghost or shadow was observed on a volume-rendered image of the White India ink-gelatin infused brain. The z-axial resolution examined on the orthogonal sections reconstituted from serial images obtained at an interval of 5 mum showed no cross talk, indicating that the z-axial resolution was no larger than 5 mum. A proper understanding of the vascular system in a whole brain is indispensable to reveal the development of the vascular system in the brain of normal and genetically manipulated mouse and vascular alterations in pathological situation, such as stroke and neurodegenerative disease. Although simple and inexpensive, this method will provide fundamental information on the vascular system in a whole brain.     

MRT letter: Extended depth from focus reconstruction method for stretch zone measurement in 15-5PH steel

The stretch zone width (SZW) data for 15‐5PH steel CTOD specimens fractured at ?150°C to + 23°C temperature were measured based on focused images and 3D maps obtained by extended depth‐of‐field reconstruction from light microscopy (LM) image stacks. This LM‐based method, with a larger lateral resolution, seems to be as effective for quantitative analysis of SZW as scanning electron microscopy (SEM) or confocal scanning laser microscopy (CSLM), permitting to clearly identify stretch zone boundaries. Despite the worst sharpness of focused images, a robust linear correlation was established to fracture toughness (K_C) and SZW data for the 15‐5PH steel tested specimens, measured at their center region. The method is an alternative to evaluate the boundaries of stretched zones, at a lower cost of implementation and training, since topographic data from elevation maps can be associated with reconstructed image, which summarizes the original contrast and brightness information. Finally, the extended depth‐of‐field method is presented here as a valuable tool for failure analysis, as a cheaper alternative to investigate rough surfaces or fracture, compared to scanning electron or confocal light microscopes. Microsc. Res. Tech. 75:1155–1158, 2012. © 2012 Wiley Periodicals, Inc.     

Enamel erosion and prevention efficacy characterized by confocal laser scanning microscope

The aim of this study was to evaluate the erosion‐inhibiting effect of two toothpastes on the development of erosion‐like lesions, by a confocal laser scanning microscope (CLSM). Forty human enamel blocks were divided into five groups (n = 8), in accordance to evaluate the GC MI Paste Plus and Oral B with stannous fluoride, applied as slurries and associated with toothbrush. Specimens were submitted to an erosion challenge from citric acid (0.5%, pH = 2.8), for 5 min, six times a day, alternating in artificial saliva immersions. Reference group was not exposed to treatment. Part of specimens (Groups 02 and 03) was exposed twice daily just to slurries, for 2 min, therefore specimens from Groups 04 and 05 were also abraded, for 30 s. The enamel surfaces were morphological characterized using CLSM images, with mineral loss being measured using the resulting 3D images referenced to an un‐challenged portion of the sample. Step values were compared using the one‐way ANOVA test. CLSM was shown to be a viable, noncontact, and simple technique to characterize eroded surfaces. The statistical difference in the step size was significant between the groups (P = 0.001) and using multiple comparisons a statistically significant protective effect of toothpastes was shown when these were applied as slurries. Although groups submitted to tooth brush showed mineral loss similar to reference control group, due to the damages of abrasion associated. Microsc. Res. Tech. 77:439–445, 2014. © 2014 Wiley Periodicals, Inc.     

A dual path programmable array microscope (PAM): simultaneous acquisition of conjugate and non-conjugate images

A programmable array microscope (PAM) incorporates a spatial light modulator (SLM) placed in the primary image plane of a widefield microscope, where it is used to define patterns of illumination and/or detection. We describe the characteristics of a special type of PAM collecting two images simultaneously. The conjugate image (I_c) is formed by light originating from the object plane and returning along the optical path of the illumination light. The non‐conjugate image (I_nc) receives light from only those regions of the SLM that are not used for illuminating the sample. The dual‐signal PAM provides much more time‐efficient excitation than the confocal laser scanning microscope (CLSM) and greater utilization of the available emission light. It has superior noise characteristics in comparison to single‐sided instruments. The axial responses of the system under a variety of conditions were measured and the behaviour of the novel I_nc image characterized. As in systems in which only I_c images are collected (Nipkow‐disc microscopes, and previously characterized PAMs), the axial response to thin fluorescent films showed a sharpening of the axial response as the unit cell of the repetitive patterns decreased in size. The dual‐signal PAM can be adapted to a wide range of data analysis and collection strategies. We investigated systematically the effects of patterns and unit cell dimensions on the axial response. Sufficiently sparse patterns lead to an I_c image formed by the superposition of the many parallel beams, each of which is equivalent to the single scanning spot of a CLSM. The sectioning capabilities of the system, as given by its axial responses, were similar for a given scan pattern and for processed pseudorandom sequence (PRS) scans with the same size of the unit cell. For the PRS scans, optical sectioning was achieved by a subtraction of an I_nc image or, alternatively, a scaled widefield image from the I_c image. Based on the comparative noise levels of the two methods, the non‐conjugate subtraction was significantly superior. A point spread function for I_c and I_nc was simulated and properties of the optical transfer functions (OTFs) were compared. Simulations of the OTF in non‐conjugate imaging did not suffer from the missing cone problem, enabling a high quality deconvolution of the non‐conjugate side alone. We also investigated the properties of images obtained by subjecting the I_c and I_nc data to a combined maximum likelihood deconvolution.     

Automated identification of neurons in 3D confocal datasets from zebrafish brainstem

Many kinds of neuroscience data are being acquired regarding the dynamic behaviour and phenotypic diversity of nerve cells. But as the size, complexity and numbers of 3D neuroanatomical datasets grow ever larger, the need for automated detection and analysis of individual neurons takes on greater importance. We describe here a method that detects and identifies neurons within confocal image stacks acquired from the zebrafish brainstem. The first step is to create a template that incorporates the location of all known neurons within a population – in this case the population of reticulospinal cells. Once created, the template is used in conjunction with a sequence of algorithms to determine the 3D location and identity of all fluorescent neurons in each confocal dataset. After an image registration step, neurons are segmented within the confocal image stack and subsequently localized to specific locations within the brainstem template – in many instances identifying neurons as specific, individual reticulospinal cells. This image-processing sequence is fully automated except for the initial selection of three registration points on a maximum projection image. In analysing confocal image stacks that ranged considerably in image quality, we found that this method correctly identified on average ∼80% of the neurons (if we assume that manual detection by experts constitutes 'ground truth'). Because this identification can be generated approximately 100 times faster than manual identification, it offers a considerable time savings for the investigation of zebrafish reticulospinal neurons. In addition to its cell identification function, this protocol might also be integrated with stereological techniques to enhance quantification of neurons in larger databases. Our focus has been on zebrafish brainstem systems, but the methods described should be applicable to diverse neural architectures including retina, hippocampus and cerebral cortex.     

An automatic segmentation algorithm for 3D cell cluster splitting using volumetric confocal images

With the rapid advance of three-dimensional (3D) confocal imaging technology, more and more 3D cellular images will be available. Segmentation of intact cells is a critical task in automated image analysis and quantification of cellular microscopic images. One of the major complications in the automatic segmentation of cellular images arises due to the fact that cells are often closely clustered. Several algorithms are proposed for segmenting cell clusters but most of them are 2D based. In other words, these algorithms are designed to segment 2D cell clusters from a single image. Given 2D segmentation methods developed, they can certainly be applied to each image slice with the 3D cellular volume to obtain the segmented cell clusters. Apparently, in such case, the 3D depth information with the volumetric images is not really used. Often, 3D reconstruction is conducted after the individualized segmentation to build the 3D cellular models from segmented 2D cellular contours. Such 2D native process is not appropriate as stacking of individually segmented 2D cells or nuclei do not necessarily form the correct and complete 3D cells or nuclei in 3D. This paper proposes a novel and efficient 3D cluster splitting algorithm based on concavity analysis and interslice spatial coherence. We have taken the advantage of using the 3D boundary points detected using higher order statistics as an input contour for performing the 3D cluster splitting algorithm. The idea is to separate the touching or overlapping cells or nuclei in a 3D native way. Experimental results show the efficiency of our algorithm for 3D microscopic cellular images.