This article has Open Peer Review reports available.
Computer-assisted counting of retinal cells by automatic segmentation after TV denoising
- Kristian Bredies^{1},
- Marcus Wagner^{2},
- Christian Schubert^{3} and
- Peter Ahnelt^{3}Email author
https://doi.org/10.1186/1471-2415-13-59
© Bredies et al.; licensee BioMed Central Ltd. 2013
Received: 17 December 2012
Accepted: 23 August 2013
Published: 20 October 2013
Abstract
Background
Quantitative evaluation of mosaics of photoreceptors and neurons is essential in studies on development, aging and degeneration of the retina. Manual counting of samples is a time consuming procedure while attempts to automatization are subject to various restrictions from biological and preparation variability leading to both over- and underestimation of cell numbers. Here we present an adaptive algorithm to overcome many of these problems.
Digital micrographs were obtained from cone photoreceptor mosaics visualized by anti-opsin immuno-cytochemistry in retinal wholemounts from a variety of mammalian species including primates. Segmentation of photoreceptors (from background, debris, blood vessels, other cell types) was performed by a procedure based on Rudin-Osher-Fatemi total variation (TV) denoising. Once 3 parameters are manually adjusted based on a sample, similarly structured images can be batch processed. The module is implemented in MATLAB and fully documented online.
Results
The object recognition procedure was tested on samples with a typical range of signal and background variations. We obtained results with error ratios of less than 10% in 16 of 18 samples and a mean error of less than 6% compared to manual counts.
Conclusions
The presented method provides a traceable module for automated acquisition of retinal cell density data. Remaining errors, including addition of background items, splitting or merging of objects might be further reduced by introduction of additional parameters. The module may be integrated into extended environments with features such as 3D-acquisition and recognition.
Keywords
Mammalian photoreceptor cells Automatical counting Adaptive algorithm Continuous optimization Total variation denoisingBackground
Introduction
The vertebrate retina contains two or more subtypes of photoreceptors and dozens of interneuron types, thus being organized for effective operation at different light levels and at different bands of the sunlight’s spectrum. Regional shifts in densities and proportions of the subtypes of photoreceptors and interneurons in the retina have been studied intensively as they are assumed to reflect both the evolution of species and specific adaptations to their lifestyle [1, 2]. Deviations within cell density and mosaic regularity from their normal variability are of specific interest for research involved in developmental control or progressive loss of photoreceptor and other cells due to degenerative diseases and other pathologic processes in the visual system.
Obtaining such data implies the acquisition and analysis of a large number of samples, which is often a time-consuming task requiring persons with appropriate training and experience. Given, moreover, the approximately planar organization of retinal layers, the problem of detection and counting of photoreceptor cells is a promising candidate for at least partial automatization. However, while various papers have developed options for mapping and analysis of retinal cell mosaic data, once they are digitized [3–7], approaches to full automatization of the time-consuming procedure of actual identification and counting have been rare (we mention [8–10]). Besides the time and computing power required for acquiring two- and three-dimensional representations of the tissue of interest, the heterogeneity of the tissue itself is still a major challenge. Many parameters such as tissue thickness/transparency, preparatory and manipulatory distortions change across the retinas, and the cells of interest themselves change in size, shape and spacing. Consequently, for reliable detection of the targets and their differentiation from other items such as debris, other cells, local damage or blood vessels highly adaptive algorithms are required.
The present approach focuses on addressing these problems for the segmentation and counting of photoreceptors. We propose an automatic detection and counting procedure, which is based on Rudin-Osher-Fatemi total variation (TV) denoising [11] with subsequent segmentation. As the comparison with manually collected data shows, this method is able to detect the targets at comparable reliability. In the view of the authors, the main advantage of the method consists in the complete traceability of all data processing steps and its reproducibility independently from a particular software platform. Thus a possibility for standardization and direct comparisons of automatic counts for samples obtained within different environments is provided.
In the present study, the method has been implemented as a MATLAB module and has been applied to single frames (and montages). The method, however, is not limited to the processing of two-dimensional data and can be equally implemented for the analysis of three- or multidimensional data stacks. It could as well be integrated into more comprehensive motorized acquisition setups for stereological sampling or complete mapping of retinal populations.
Retinal image data
Retinal image data
No. | Species | Eye | Retina region | Labeled structures | Staining | Density per mm^{2} | File name | Image size | Remarks about image quality and particularities | Shown in Figure(s) |
---|---|---|---|---|---|---|---|---|---|---|
1 | Orangutan | Left | Sup. periphery, | S-cones, | JH455/DAB, | 190 (cones) | OrangWhM | 1024 ×854 | Sufficient quality, | 1A |
between 10 mm from | blood vessels | Collagen IV/DAB | but vessels present | |||||||
fovea and ora serrata | ||||||||||
2 | Orangutan | Left | Inf. periphery, | S-cones, | JH455/DAB, | 156 (cones) | OrUinf-0001pr | 1026 ×854 | Sufficient quality, | 4I–J |
between 11 mm from | blood vessels | Collagen IV/DAB | but vessels present | |||||||
fovea and ora serrata | ||||||||||
3 | Domestic cat | S-cones | JH455/DAB | 860 | loStack | 1024 ×1024 | Sufficient quality and contrast, | |||
but small cones | ||||||||||
4 | Domestic cat | Inf. retina | S-cones | JH455/DAB | 851 | luStack | 1024 ×1024 | Sufficient quality and contrast, | 1B | |
but small cones | ||||||||||
5 | Domestic cat | S-cones | JH455/DAB | 827 | roStack | 1024 ×1024 | Sufficient quality and contrast, | |||
but small cones | ||||||||||
6 | Domestic cat | S-cones | JH455/DAB | 904 | ruStack | 1024 ×1024 | Sufficient quality and contrast, | |||
but small cones | ||||||||||
7 | Manul | S-cones | JH455/DAB | 195 | bcergcomp | 1024 ×1024 | Sufficient quality and contrast, | |||
but small cones | ||||||||||
8 | Eurasian Lynx | Left | Periphery | S-cones | JH455/DAB | 1076 | Movie-17-1 | 1024 ×1024 | Good quality and contrast | 4A–B |
9 | Cheetah | S-cones | JH455/DAB | 130 | 1+5CSCStack | 1024 ×1024 | Good quality and contrast | |||
10 | Cheetah | S-cones | JH455/DAB | 1055 | Stack SCFoc. | 1024 ×1024 | Good quality and contrast | |||
11 | Jaguar | Right | Periphery | S-cones, | JH455/DAB | 440 (cones) | Combination | 1024 ×1024 | Good quality, but additional | 1C, 4G–H |
horizontal cells | cell type present | |||||||||
12 | Pangolin | Right | Mid periphery | M-/L-cones | JH492/DAB | 2290 | 492_20_09 | 1024 ×1024 | Background very blotchy | 1D, 2A–F, |
4C–D | ||||||||||
13 | Agouti | Left | Sup. temp. | S-cones | JH455/DAB | 705 | A4scal | 1024 ×853 | Low contrast, cell debris present | 1E |
periphery | ||||||||||
14 | Agouti | S-cones | JH455/DAB | 790 | A5scal | 1024 ×853 | Low contrast, cell debris present | |||
15 | Agouti | S-cones | JH455/DAB | 855 | A6scal | 1024 ×853 | Low contrast, cell debris present | |||
16 | Agouti | left | Temp. periphery | S-cones | JH455/DAB | 480 | B1scal | 1024 ×853 | Poor quality, cell debris present | 1F, 4E–F |
17 | Agouti | S-cones | JH455/DAB | 905 | B2scal | 1024 ×853 | Poor quality, cell debris present | |||
18 | Agouti | S-cones | JH455/DAB | 3500 | C4scal | 1024 ×853 | Poor quality, cell debris present |
Most of the eyes were obtained from animals delivered to veterinary pathology from zoos and animal parks; some originate from collaborations for other studies [12, 13]. Post mortem times were between 0.5 and 12 hours. Eyes were enucleated and immersed in 0.01M phosphate buffer saline (PBS, pH 7.4) with 4% paraformaldehyde. Some were treated after being opened with a cut along the corneal limbus for faster penetration of fixative. Retina wholemounts were prepared in PBS and flattened by radial cuts, in order to preserve the horizontal and vertical meridian. Cone photoreceptors sensitive to medium/long wavelengths (M-/L-cones) were labeled in isolated Pangolin retina using JH492 antibody, in all other retinas S-cones were labeled using JH455 (both antibodies provided by J. Nathans [14]). In peripheral Jaguar retina, in addition to S-cones horizontal cells are (unintentionally) co-labeled by JH455. Retinal vessels in Orangutan were labeled with rabbit anti-mouse collagen IV (AbD Serotec, 2150-1470). After incubation in primary antisera overnight for up to 3 days visualization was done using goat anti-rabbit igg-biotin (Sigma, B7389), ExtrAvidin-peroxidase conjugate (Sigma, E2886) and the diaminobenzidine (DAB) reaction. After washing in PBS, retinas were gradually transferred up to 90% glycerol, mounted with photoreceptor-side up on a glass slide, and cover slipped.
Manual counting of labeled cones was done within sampling frames of 150×150μm or 300×300μm using an online-video overlay system consisting of Canvas 5 (ACD Systems, USA) software on a Macintosh computer connected with a Hamamatsu 2400 analog camera attached to a Nikon Eclipse microscope. This system allows dual live view of the specimen: through the microscope’s optics or on the video image overlaid by the sampling frames. Optional change of focus and illumination supports optimized online identification of cells and exclusion of artifacts by position, form, color and other details.
The images (8 or 10 bit grey scale) used for computer-assisted cell counting were obtained by using a Photometrix Camera model CH250/A connected to a Nikon Eclipse E600 microscope (magnification factors 200× to 600×) using QED Imaging Software (QED Imaging Inc., Pittsburgh, PA) on a Macintosh computer. In most cases, a projection image (max density or sum) was composed from a stack of images at relevant focus levels using the public domain NIH Image program (developed at the U.S. National Institutes of Health; available at http://rsb.info.nih.gov/nih-image/ (accessed 11.02. 2013)).
For testing the computer-assisted cell counting method, images were chosen by the following criteria: retinas from different species, differences in quality concerning morphology (different amount of cell debris, broken cones etc.) and contrast between labeled cells and background (see Table 1). Images with additional labeled structures like horizontal cells or blood vessels have been analyzed as well. The respective objects of interest were then manually counted by an experienced coauthor (C.S.), and the results were tabulated for comparison with the outcome from the automatic counting method (see below).
Description of the detection and counting method
as “black enough” to belong to images of photoreceptor cells. Finally, all “black” features consisting of a number of connected adjacent pixels, which is bigger than a given threshold size f, are automatically counted, cf. Figure 2D.
In this method, no more than three parameters remain to be adjusted manually. These are: the window size m for the median filter, the parameter c in (2), which influences the contrast differentiation between photoreceptor cells and the background, and the minimal size f of a connected feature to be recognized as a photoreceptor cell.
Our algorithm can be summarized as follows:
Algorithm 1 Automatic segmentation after TV denoising
Implementation
Implementation as a MATLAB tool
The initializations are ${x}^{0}={\stackrel{\u0304}{x}}^{0}={I}^{\left(1\right)}$, p ^{(1),0} = p ^{(2),0} = 0, and the output will be given by $\left({x}_{\mathit{\text{ij}}}\right)=\left({x}_{\mathit{\text{ij}}}^{N}\right)$. As default values of the regularization parameter and the number of iterations, we use α=0.05 and N=50. Note that, in principle, other minimization algorithms could be implemented for solving the TV denoising problem. For instance, we mention the generalized TV approach from [16] and the optimal control method described in [17], which is based on the interior-point solver IPOPT [18, 19].
The thresholding in Step 3 is realized with the help of the built-in MATLAB functions mean and var. Finally, for Step 4, the labelling and counting procedure, the function bwlabel is utilized, which is again part of the image processing toolbox. It yields a labelled image in which each connected component is identified by a positive integer. With this information, the identification and counting of those connected components, which comprise at least f pixels, can be easily realized.
Usage
Usually, in order to analyze the topography of items in a retina preparation, a considerable number of image files has to be generated, each showing a segment. Our software tool was especially designed to cope with multiple files showing similar structures. In this situation, it is possible to start with a manual count within one or two typical images in order to calibrate the parameters m, c and f. This can be done by starting the program and adding a single image to the file list. In most cases, the default parameters give a good starting point, hence one can perform a segmentation in order to decide whether one is satisfied with the results. Otherwise, adjust the parameters m, c and f using the configuration dialog (cf. Figure 3B) and try again. Once the parameters have been adjusted, they can be utilized for the analysis of the whole image set. The batch processing feature of the software allows to perform this analysis without further user interaction: simply add the remaining files to the list and start the segmentation procedure. Finally, a report which lists, for each file in the batch, the number and mean density of detected cells as well as their positions, can be automatically generated. As default values for the parameters m, c and f, we encoded m=30, c=2.5 and f=5. Eventually, for convenience and reproducibility, the tool also provides saving and loading of the file list as well as of a list of all parameters.
Results and discussion
Overall cell counts by different methods
No. | Image | Manual count # | Direct segm. | TV denoising & segm. | m /c /f |
---|---|---|---|---|---|
# / rel. error | # / rel. error | ||||
Orangutan | |||||
1 | OrangWhM | 2328 | 1784 / 24.9 | 2286 / 1.8 | 30 / 2.0 / 5 |
2 | OrUinf-0001pr | 1192 | 793 / 33.5 | 1168 / 2.0 | 30 / 3.0 / 5 |
Domestic cat | |||||
3 | loStack | 557 | 518 / 7.0 | 543 / 2.5 | 30 / 2.5 / 5 |
4 | luStack | 552 | 584 / 5.8 | 562 / 1.8 | 30 / 2.5 / 5 |
5 | roStack | 536 | 539 / 0.6 | 520 / 3.0 | 30 / 2.5 / 5 |
6 | ruStack | 584 | 573 / 1.9 | 573 / 1.9 | 30 / 2.5 / 5 |
Manul | |||||
7 | bcergcomp | 128 | 110 / 14.1 | 128 / 0.0 | 20 / 2.5 / 25 |
Lynx | |||||
8 | Movie-17-1 | 708 | 675 / 4.7 | 659 / 6.9 | 30 / 2.5 / 5 |
Cheetah | |||||
9 | 1+5CSCStack | 85 | 84 / 1.2 | 91 / 7.1 | 30 / 3.5 / 40 |
10 | Stack SCFoc. | 687 | 656 / 4.5 | 634 / 7.7 | 30 / 2.5 / 5 |
Jaguar | |||||
11 | Combination | 289 (S-cones) | 177 / 38.8 | 275 / 4.8 | 30 / 3.5 / 5 |
Pangolin | |||||
12 | 492_20_09 | 369 | 109 / 70.5 | 367 / 0.5 | 20 / 2.5 / 25 |
Agouti | |||||
13 | A4scal | 43 | 14 / 67.4 | 44 / 2.3 | 40 / 2.5 / 150 |
14 | A5scal | 48 | 20 / 58.3 | 54 / 12.5 | 40 / 2.5 / 150 |
15 | A6scal | 52 | 44 / 15.4 | 53 / 1.9 | 40 / 2.5 / 150 |
16 | B1scal | 31 | 35 / 12.9 | 43 / 38.7 | 40 / 2.5 / 150 |
17 | B2scal | 55 | 51 / 7.3 | 56 / 1.8 | 40 / 2.5 / 150 |
18 | C4scal | 213 (?) | 95 / 55.4 | 234 / 9.8 | 30 / 2.5 / 50 |
Mean error | 23.6 | 5.9 | |||
Standard dev. | 23.6 | 8.6 |
Number of correctly recognized cells within the automatic counts
No. | Image # | Manual count | Direct segm. | TV denoising & segm. | m /c /f |
---|---|---|---|---|---|
# / percent | # / percent | ||||
Orangutan | |||||
1 | OrangWhM | 2328 | 1652 / 71.0 | 2176 / 93.4 | 30 / 2.0 / 5 |
2 | OrUinf-0001pr | 1192 | 735 / 61.7 | 1077 / 90.4 | 30 / 3.0 / 5 |
Domestic cat | |||||
3 | loStack | 557 | 474 / 85.1 | 512 / 91.9 | 30 / 2.5 / 5 |
4 | luStack | 552 | 531 / 96.2 | 527 / 95.5 | 30 / 2.5 / 5 |
5 | roStack | 536 | 502 / 93.6 | 499 / 93.1 | 30 / 2.5 / 5 |
6 | ruStack | 584 | 530 / 90.8 | 545 / 93.3 | 30 / 2.5 / 5 |
Manul | |||||
7 | bcergcomp | 128 | 80 / 62.5 | 126 / 98.4 | 20 / 2.5 / 25 |
Lynx | |||||
8 | Movie-17–1 | 708 | 642 / 90.7 | 640 / 90.4 | 30 / 2.5 / 5 |
Cheetah | |||||
9 | 1+5CSCStack | 85 | 81 / 95.3 | 85 / 100.0 | 30 / 3.5 / 40 |
10 | Stack SCFoc. | 687 | 637 / 92.7 | 624 / 90.8 | 30 / 2.5 / 5 |
Jaguar | |||||
11 | Combination | 289 (cones) | 164 / 56.7 | 253 / 87.5 | 30 / 3.5 / 5 |
Pangolin | |||||
12 | 492_20_09 | 369 | 104 / 28.2 | 297 / 80.5 | 20 / 2.5 / 25 |
Agouti | |||||
13 | A4scal | 43 | 12 / 27.9 | 40 / 93.0 | 40 / 2.5 / 150 |
14 | A5scal | 48 | 16 / 33.3 | 45 / 93.8 | 40 / 2.5 / 150 |
15 | A6scal | 52 | 40 / 76.9 | 50 / 96.2 | 40 / 2.5 / 150 |
16 | B1scal | 31 | 29 / 93.5 | 31 / 100.0 | 40 / 2.5 / 150 |
17 | B2scal | 55 | 50 / 90.9 | 54 / 98.2 | 40 / 2.5 / 150 |
18 | C4scal | 213 (?) | 77 / 36.1 | 185 / 86.8 | 30 / 2.5 / 50 |
Mean percentage | 71.3 | 93.0 | |||
Standard dev. | 24.5 | 4.8 |
Artifacts within the automatic counts produced by the different methods
No. | Image | Direct segm. | TV denoising & segm. | m /c /f |
---|---|---|---|---|
# / percent | # / percent | |||
Orangutan | ||||
1 | OrangWhM | 132 / 7.4 | 110 / 4.8 | 30 / 2.0 / 5 |
2 | OrUinf-0001pr | 58 / 7.3 | 91 / 7.8 | 30 / 3.0 / 5 |
Domestic cat | ||||
3 | loStack | 44 / 8.5 | 31 / 5.7 | 30 / 2.5 / 5 |
4 | luStack | 53 / 9.1 | 35 / 6.2 | 30 / 2.5 / 5 |
5 | roStack | 37 / 6.9 | 21 / 4.0 | 30 / 2.5 / 5 |
6 | ruStack | 43 / 7.5 | 28 / 4.9 | 30 / 2.5 / 5 |
Manul | ||||
7 | bcergcomp | 30 / 23.4 | 2 / 1.6 | 20 / 2.5 / 25 |
Lynx | ||||
8 | Movie-17–1 | 33 / 4.9 | 19 / 2.9 | 30 / 2.5 / 5 |
Cheetah | ||||
9 | 1+5CSCStack | 3 / 3.6 | 6 / 6.6 | 30 / 3.5 / 40 |
10 | Stack SCFoc. | 19 / 2.9 | 10 / 1.6 | 30 / 2.5 / 5 |
Jaguar | ||||
11 | Combination | 13 / 7.3 | 22 / 8.0 | 30 / 3.5 / 5 |
Pangolin | ||||
12 | 492_20_09 | 5 / 4.6 | 70 / 19.1 | 20 / 2.5 / 25 |
Agouti | ||||
13 | A4scal | 2 / 14.3 | 4 / 9.1 | 40 / 2.5 / 150 |
14 | A5scal | 4 / 20.0 | 9 / 16.7 | 40 / 2.5 / 150 |
15 | A6scal | 4 / 9.1 | 3 / 5.7 | 40 / 2.5 / 150 |
16 | B1scal | 6 / 17.1 | 12 / 27.9 | 40 / 2.5 / 150 |
17 | B2scal | 1 / 2.0 | 2 / 3.6 | 40 / 2.5 / 150 |
18 | C4scal | 18 / 18.9 | 49 / 20.9 | 30 / 2.5 / 50 |
Mean percentage | 9.7 | 8.7 | ||
Standard dev. | 6.1 | 7.2 |
The quality requirements for automatic counts depend on the specific goal of interest. When assessing general distribution patterns or gradients for certain types of receptors, a relative error of 10% may be acceptable while for the analysis of degeneration or proliferation phenomena or the detection of initial points of density changes, an error of about 5% and less is desirable. Table 2 shows that the latter goal has been realized by the TV/segmentation method in 12 of 18 cases while the error is below 10% in 16 of 18 cases. The mean error amounts to 5.9%. Moreover, Table 3 shows that the method recognizes correctly more than 90% of the manually marked cells in 15 of 18 cases (93.0% in the mean). The number of artifacts contained in the counts, as listed in Table 4, amounts to 8.7% in the mean. Thus, given that the TV/segmentation method makes no use of additional information about the shape of the cells or the variation of their size, these results are quite satisfactory.
As to be expected, our results show that automatic counting by the TV/segmentation method is superior to direct segmentation in every respect. The mean relative error produced by the latter amounts to as much as 23.6%, and the relative error goes below 10% in 8 of 18 cases only. The loss of precision is mostly caused by the fact that, by direct segmentation, a significantly smaller number of marked cells is recognized than by the TV/segmentation method (cf. Figure 2D–F) while the ratio of artifacts produced by both methods is comparable. The superiority of the TV/segmentation method can be seen as well by comparing the standard deviations of the indicators. Let us remark that, during our experiments, we further observed that the TV/segmentation method is even superior to direct segmentation after subtraction of the median (Steps 1, 3 and 4 of Algorithm 1).
As the patterns within the Orangutan retina closely resemble the organization within the human one, similar results are to be expected for human samples.
Let us briefly compare the proposed method with other approaches pursued in the literature. In [8], the authors perform the image processing steps by use of a commercial software package which, unfortunately, comes as a “black box” without documentation of the internally utilized algorithms. In [9, 10, 20, 21], after certain presmoothing/denoising steps, watershed segmentation is employed. Additionally, before segmenting, in [20] an illumination correction is performed while the authors in [10] interpose a contrast enhancement step. A common feature of all approaches is the necessity to select a number of parameters, including the (expected or minimal) feature size, by the experimenter.
Although preprocessing steps, particularly denoising or smoothing of the raw image data, are crucial for the quality of the results of subsequent segmentation, they have not been thoroughly documented in the cited references (if at all), and their dependence on additional, manually tuned parameters remains unclear. Moreover, any denoising method generates artifacts, thus modifying fine structures within the images in a specific way. In contrast to this situation, the preprocessing steps involved in our method (median filtering and TV denoising) are traceable and reproducible, including the manual setting of the single parameter m. For the TV denoising method, the characteristic artifacts are well-investigated (see e. g. the discussion in [16], pp. 519 ff.). In fact, this method has been deliberately selected in order to take advantage of its well-known “cartooning” effect. As a further difference, we segment by intensity thresholding instead of watershedding. The latter approach is well suited for the analysis of large, clumpy cell aggregates while our method is better suited for the detection of single cells to be differentiated from a more or less blotchy background, which may contain additional structures like vessels or different cell types. The segmentation depends on no more than two further parameters c and f, which have to be selected on the base of the observed contrast as well as of a reasonable guess of the feature size. A possible improvement could be the introduction of an additional upper bound for the size of the recognized features, thus reducing and possibly underestimating the number of cells since larger aggregates formed of merged cell images will then be excluded.
The reliability of our method, when evaluated by the mean relative error of automatic counting in relation to manual counts (5.9% as documented above), fits well within the range of errors documented in the cited references: [9], p. 1969, Figure one: ≤ 10% in 11 of 23 cases; [10], p. 641, Figure two: ≤ 5% in 31 of 40 cases; [20], p. R100.7, Figure two(A): 6% and 17%; [21], p. 589, Figure three: ≤ 10% in roughly half of the cases.
Due to its traceability, the TV denoising method offers the further advantage of easy reimplementation. On the one hand, this is possible with respect to the use of non-proprietary software, on the other hand, the method may be carried over (after optimization of the code as necessary) to the analysis of three- or multidimensional data stacks. Concerning the runtime behaviour, the analysis of an 1024 ×1024-pixel image takes typically less than 40 sec (on a Mini-PC equipped with four processors Intel(R)Core(TM)i3 CPU M380 @ 2.53 GHz) where approximately half of the time is consumed by the median filtering procedure. No particular attempts for tuning have been made.
Conclusions
In conclusion, the presented approach provides a reproducible method for the segmentation of labeled cellular objects, in particular retinal photoreceptors from grey scale micrographs. Relatively simple outlines and little overlap between the targets in the projection plane are clearly important preconditions for the current implementation, but within these limitations the procedure is shown to deliver robust results well comparable to those of manual counts by experienced observers. The method is well documented and can be transferred to other tasks including those with more automatization for both image acquisition and further refinements such as from 2D- or 3D-object recognition criteria.
Availability and requirements
The software, to which applies the GNU General Public License v.2, is stored at the location http://www.meduniwien.ac.at/counttool/ within the archive cacount_tool_2013_02_08.zip. Its execution requires a current version of MATLAB (e.g., v. 7.14.0.739 (R2012a) and higher) together with the MATLAB Image Processing Toolbox (e.g., v. 8.0 (R2012a) and higher). Within this environment, the tool runs platform-independently. The code has been written in MATLAB; its listing will be provided as Additional file 1. To use the tool, deflate the zip-archive, include its location as well as location of the image data into the MATLAB path and type the command main. Further details on the usage can be found in the accompanying online documentation.
Appendix
Rudin-Osher-Fatemi TV denoising: mathematical background
over suitable function spaces, e. g. spaces of Sobolev functions or functions of bounded variation. The first member within F, the data fidelity term, aims for a least-square approximation of the captured data I(s) while the second one, the so-called regularization term, has been purposely introduced in order to ensure existence as well as uniqueness of a minimizing solution x(s). The influence of the second term is weighted by a number α > 0, which is called the regularization parameter. Note that the regularization term depends on the gradient ∇x(s), thus favorizing a certain edge structure within the minimizing solution x(s). A rigorous mathematical development of this idea relies on a closer analysis of the Euler-Lagrange equation, which must be satisfied as a second-order PDE by the minimizers of (9) as a necessary condition, cf. [22], pp. 64–66.
and the minimization is performed over all functions of bounded variation on Ω (for more details, cf. [24], pp. 355 ff., and [25], pp. 175 ff.). Loosely speaking, | x |_{ T V } is a measure for the oscillation of the function x(s) and depends on the (generalized) derivatives of x(s) through a duality expression. A closer mathematical analysis reveals that the solution x(s) of the Rudin-Osher-Fatemi TV denoising problem largely preserves the edge structure of the input data I(s). After a suitable discretization, the particular structure of | x |_{ T V } allows for the application of highly efficient primal-dual numerical solvers like the one utilized in our present study, cf. again [25], pp. 179 ff.
Declarations
Acknowledgements
PA and CS have been supported by the Austrian Science Foundation grant I-433 B13 within the E-Rare joint project Rhorcod. KB and MW have been supported by the Austrian Science Foundation within the Special Research Program F 32 “Mathematical Optimization and Applications in Biomedical Sciences” (Graz).
Authors’ Affiliations
References
- Ahnelt PK, Kolb H: The mammalian photoreceptor mosaic-adaptive design. Prog Retinal Eye Res. 2000, 19: 711-777. 10.1016/S1350-9462(00)00012-4.View ArticleGoogle Scholar
- Peichl L: Diversity of mammalian photoreceptor properties: adaptations to habitat and lifestyle?. Anat Rec A Discov Mol Cell Evol Biol. 2005, 287: 1001-1012.View ArticlePubMedGoogle Scholar
- Curcio CA, Sloan KR: Packing geometry of human cone photoreceptors: variation with eccentricity and evidence for local anisotropy. Visual Neurosci. 1992, 9: 169-180. 10.1017/S0952523800009639.View ArticleGoogle Scholar
- Fernández E, Cuenca N, De Juan: A compiled BASIC program for analysis of spatial point patterns: application to retinal studies. J Neurosci Methods. 1993, 50: 1-15. 10.1016/0165-0270(93)90051-R.View ArticlePubMedGoogle Scholar
- Galli-Resta L, Novelli E, Kryger Z, Jacobs GH, Reese BE: Modelling the mosaic organization of rod and cone photoreceptors with a minimal-spacing rule. Eur J Neurosci. 1999, 11: 1461-1469. 10.1046/j.1460-9568.1999.00555.x.View ArticlePubMedGoogle Scholar
- Martinez Mozos, Bolea JA, Ferrandez JM, Ahnelt PK, Fernandez E: V-Proportion: a method based on the Voronoi diagram to study spatial relations in neuronal mosaics of the retina. Neurocomputing. 2010, 74: 418-427. 10.1016/j.neucom.2010.07.020.View ArticleGoogle Scholar
- Wässle H, Riemann HJ: The mosaic of nerve cells in the mammalian retina. Proc Roy Soc B — Biol Sci. 1978, 200 (1141): 441-461. 10.1098/rspb.1978.0026.View ArticleGoogle Scholar
- Clérin E, Wicker N, Mohand-Saïd S, Poch O, Sahel JA, Léveillard T: e-conome: an automated tissue counting platform of cone photoreceptors for rodent models of retinitis pigmentosa. BMC Ophthalmol. 2011, 11: 38-10.1186/1471-2415-11-38. 10.1186/1471-2415-11-38View ArticlePubMedPubMed CentralGoogle Scholar
- Filippopoulos T, Danias J, Chen B, Podos SM, Mittag TW: Topographic and morphologic analyses of retinal ganglion cell loss in old DBA/2NNia mice. Invest Ophthalmol Visual Sci. 2006, 47: 1968-1974. 10.1167/iovs.05-0955.View ArticleGoogle Scholar
- Salinas-Navarro M, Mayor-Torroglosa S, Jiménez-López M, Avilés-Trigueros M, Holmes TM, Lund RD, Villegas-Pérez MP, Vidal-Sanz M: A computerized analysis of the entire retinal ganglion cell population and its spatial distribution in adult rats. Vis Res. 2009, 49: 115-126. 10.1016/j.visres.2008.09.029.View ArticlePubMedGoogle Scholar
- Rudin LI, Osher S, Fatemi E: Nonlinear total variation based noise removal algorithms. Physica D. 1992, 60: 259-268. 10.1016/0167-2789(92)90242-F.View ArticleGoogle Scholar
- Ahnelt PK, Fernández E, Martinez O, Bolea JA, Kübber-Heiss A: Irregular S-cone mosaics in felid retinas. Spatial interaction with axonless horizontal cells, revealed by cross correlation. J Opt Soc Am A: Optics Image Sci Vis. 2000, 17: 580-588. 10.1364/JOSAA.17.000580.View ArticleGoogle Scholar
- Ahnelt PK, Schubert C, Kübber-Heiss A, Schiviz A, Anger E: Independent variation of retinal S and M cone photoreceptor topographies: a survey of four families of mammals. Vis Neurosci. 2006, 23: 429-435.View ArticlePubMedGoogle Scholar
- Chiu MI, Nathans J: A sequence upstream of the mouse blue visual pigment gene directs blue cone-specific transgene expression in mouse retinas. Vis Neurosci. 1994, 11: 773-780. 10.1017/S0952523800003072.View ArticlePubMedGoogle Scholar
- Chambolle A, Pock T: A first-order primal-dual algorithm for convex problems with applications to imaging. J Math Imaging Vis. 2011, 40: 120-145. 10.1007/s10851-010-0251-1.View ArticleGoogle Scholar
- Bredies K, Kunisch K, Pock T: Total generalized variation. SIAM J Imaging Sci. 2010, 3: 492-526. 10.1137/090769521.View ArticleGoogle Scholar
- Franek L, Franek M, Maurer H, Wagner M: A discretization method for the numerical solution of Dieudonné-Rashevsky type problems with application to edge detection within noisy image data. Opt Control Appl Meth. 2012, 33: 276-301. 10.1002/oca.996.View ArticleGoogle Scholar
- Laird C, Wächter A: Introduction to IPOPT: a tutorial for downloading, installing, and using IPOPT. Revision No. 1863. Electronically published: http://www.coin-or.org/Ipopt/documentation/ (accessed at 11.02.2013)
- Wächter A, Biegler LT: On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Math Program Ser A. 2006, 106: 25-57. 10.1007/s10107-004-0559-y.View ArticleGoogle Scholar
- Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang I, Friman O, Guertin DA, Chang J, Lindquist RA, Moffat J, Golland P, Sabatini DM: CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 2006, 7: R100-10.1186/gb-2006-7-10-r100. 10.1186/gb-2006-7-10-r100.sView ArticlePubMedPubMed CentralGoogle Scholar
- Danias J, Shen F, Goldblum D, Chen B, Ramos-Esteban J, Podos SM, Mittag T: Cytoarchitecture of the retinal ganglion cells in the rat. Invest Ophthalmol Vis Sci. 2002, 43: 587-594.PubMedGoogle Scholar
- Aubert G, Kornprobst P: Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations, 2nd ed. 2006, New York: SpringerGoogle Scholar
- Chambolle A: Mathematical Problems In Image Processing. Inverse Problems In Image Processing And Image Segmentation: Some Mathematical And Numerical Aspects. : , ICTP Lecture Notes, II. Trieste: Abdus Salam International Centre for Theoretical Physics; 2000. (electronic)Google Scholar
- Bredies K, Lorenz D: Mathematische Bildverarbeitung. Einführung in Grundlagen und moderne Theorie. : , Wiesbaden: Vieweg + Teubner Verlag / Springer Fachmedien Wiesbaden GmbH; 2011Google Scholar
- Chan TF, Shen J: Image Processing and Analysis. Variational, PDE, Wavelet, and Stochastic Methods. : , Philadelphia: SIAM; 2005Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2415/13/59/prepub
Pre-publication history
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.