Solid state nuclear track detectors (SSNTD) are materials such as inorganic crystals, plastics and glasses, known to record the path of charged particles. There are several applications for SSNTD in nuclear science; for instance, measurements of radon gas [1], boron neutron capture therapy [2], and age determination by fission track dating [3, 4, 5, 6, 7, 8].
When a charged particle collide in a SSNTD, the ionization and/or the collision with atoms modify the path by which the particle go through. This damage in the SSNTD structure is called latent track. This track becomes visible under an optical microscope after a convenient etching process [9], and these tracks can be counted. Photomicrographs can also be used, which could improve the counting accuracy [10].
Procedures for measuring and counting tracks are time-consuming and involve practical problems, e.g. variation in observer efficiency [11]. An automatic method based on image processing techniques could increase the track counting rate and improve counting reproducibility. However, separating elements in nontrivial images is one of the hardest tasks in image processing [12].
Automatic systems for separating, counting or measuring tracks have been studied for a while, and several solutions were presented (e.g. [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 11, 28, 29]. Still, the precision of automatic methods is not satisfactory yet; an automatic analysis could need to be manually adjusted by the operator, being more time consuming than the usual measure [26, 30]. The major challenges to automatic track counting are detecting overlapping tracks, distinguishing tracks and material defects (e.g. surface scratches due to polishing), and identifying small tracks and defects of comparable size in the background of photomicrographs [11].
To address the problem of identifying overlapping ion tracks in photomicrographs, we propose an algorithm based on the watershed transform [31] using morphological erosions [32] as markers. A similar method was used to separate packings of ellipsoidal particles represented using X-Ray tomography [33]. We tested this method in two datasets of diallyl phthalate (DAP) photomicrographs, and used the results to relate the incident track energy with the mean gray levels and mean diameter products for each sample from the first dataset.
We employed the ISODATA threshold to obtain the binary images used in our tests. The WUSEM algorithm is based on the following techniques: 1. morphological erosion; 2. watershed transform, and 3. labeling. We describe these algorithms in this Section.
2.1. The ISODATA threshold
The Iterative Self-Organizing Data Analysis Technique (A) (ISODATA) threshold [34, 35] is an histogram-based method. It returns a threshold that separates the image into two pixel classes, where the threshold intensity is halfway between their mean intensities.
When applied for two classes, ISODATA is always convergent [36]; in our case, these classes are the tracks (our regions of interest, ROI) and the background. We used the algorithm filters.threshold isodata, implemented in scikit-image [37].
2.2. Structuring elements and morphological erosion
In morphological image processing, a structuring element is a matrix representing a mask, or a shape, that is used to retrieve information about the shapes in an input image [32]. Erosion, a basic operation in image processing, uses a chosen structuring element to shrink the border of all ROI in a binary image; the shrinkage factor is correspondent to the structuring element size [12].
In our tests (Section 3) we used disks as structuring elements, since tracks in DAP are mostly round in shape. The algorithms used were morphology.disk and morphology.erosion, contained in scikit-image.
2.3. Watershed transform
The watershed algorithm is a non-parametric method which defines the contours as the watershed of the gradient modulus of the gray levels of the input image, considered as a relief surface detection method [31].
In this algorithm, the input image is seen in a three-dimensional perspective. Two dimensions correspond to spatial coordinates, and the third represents the gray levels. In this interpretation, we consider three kinds of points [12]:
1. Points in regional minima.
2. Points where a drop of water would flow to a common minimum.
3. Points where a drop of water would flow to different minima. The set of
The aim of watershed algorithms is to find the watershed lines. The method used in this paper is implemented in the function morphology.watershed, from scikit-image.
2.4. Labeling algorithm
In image processing, the labeling algorithm labels connected regions of a binary input image, according to the 2-connectivity sense: all eight pixels surrounding the reference pixel. Pixels receive the same label when they are connected and have the same value. We used the algorithm measure.label from scikit-image, implemented as described in [38].
2.5. Watershed Using Successive Erosions as Markers (WUSEM) algorithm
Here we present the WUSEM (Watershed Using Successive Erosions as Markers) algorithm, which combines morphological erosions, the watershed transform and labeling algorithms to separate regions of interest (ROI) in binary images. The WUSEM algorithm and its steps follow (Figure 1):
1. The user define an initial radius () and an iterative radius (∆r) to a
Figure 1: Representing WUSEM and the functions employed in its implementation.
2. The input binary image is eroded using the structuring element with radius
3. A new structuring element is defined. Its radius is + ∆r. The input
4. The process continues until the eroded image does not have objects. Then,
5. The result is labeled again, to reorder the ROI, and labels with area smaller
WUSEM is implemented in the function segmentation wusem(), available in the Supplementary Material. It receives the arguments str el, initial radius and delta radius, representing the structuring elements, and ∆r, respectively. To exemplify WUSEM’s capabilities, we used it to separate overlapping tracks within the test photomicrographs (Figure 3).
2.6. DAP photomicrographs
We used two sets of diallyl phtalate (DAP, C) photomicrographs to test the WUSEM algorithm. One of them was obtained from detectors irradiated with
tracks, and the other has induced fission tracks.
Figure 2: Representation of disks as structuring elements defined for each segmentation when employing WUSEM in an interest region. The radius of the first, second and third structuring elements are , respectively.
2.6.1. 78Kr tracks
The first dataset contains 362 photomicrographs of tracks from nine different DAP plaques irradiated with ions at a nominal fluence of 1
105
, from a beam perpendicular to the detector surfaces at GSI, Darmstadt, Germany
.
During the irradiation, detectors were covered with aluminum foils forming a moderation layers of thicknesses varying from zero (no cover) to 90 . Detectors were named after their cover thicknesses (K0, K20, K30, ..., K90). Ions arrived at the setup with initial energy of 865 MeV and were slowed down in the aluminum cover before hitting the detector surface. Incidence energies were
Figure 3: WUSEM algorithm application in an input image. First, the image is binarized using the ISODATA threshold. Then, the image is eroded using initial radius () and iterative radius (∆r) equal to 1 and 1 respectively, to ease visualization. The process continues until the eroded image has regions in it. All erosions are summed, and the result is labeled; then, regions with area smaller than 64 px are excluded. Each labeled track receives a color from the nipy spectral colormap. Finally, the function enumerate objects() is used to number the found tracks. Final results are shown to
equal to 10 and 4, respectively. Animation
calculated using the software SRIM [40], and varied from 18 MeV (90 cover) to 865 MeV (no cover). Detectors were etched in a PEW solution (7.5 g KOH, 32.5 g ethanol, 10 g water) solution for 4
2 min at 65
3
.
Images were captured with a CCD camera coupled with a Zeiss microscope, in reflected light mode, under 1250 nominal magnification. Then, the detectors were further etched for 4 minutes, total of 8
3 min, and new images were captured. This way, we obtained 18 subsets of images. Tracks within these photomicrographs are almost isotropic and have the same orientation, resembling circles (Figure 4).
Figure 4: Photomicrograph from the first dataset, presenting tracks in a DAP sample.
2.6.2. Induced fission tracks
The second dataset contains 19 photomicrographs with two different mag-nifications from a DAP plaque used as external detector, coupled to an apatite sample and irradiated with thermal neutrons to induce fission in the atoms inside the mineral
. During the fission process, two fragments are released and, eventually, are detected by the DAP plaque. Fragments arrive at the detector surface with different energies and incidence angles, resulting in a variety of track formats. Hence, counting tracks in these photomicrographs is more complex than in the previous case (Figure 5).
2.7. License and reusability
The WUSEM algorithm and several functions for its implementation lie within the packages Numpy [41], Scipy [42], Matplotlib [43], scikit-image [37], among others. All code published with this paper is written in Python 3 [44]
Figure 5: Input test photomicrograph from the second dataset, presenting tracks in a DAP sample.
and is available under the GNU GPL v3.0 license, and all photomicrographs and figures distributed with this paper are available under the CC-BY 2.0 license.
3.1. Processing photomicrographs of 78Kr tracks
3.1.1. Exemplifying the methodology
Here we use a photomicrograph from the first datasetto exemplify WUSEM (Figure 4). We binarized this photomicrograph using the ISODATA threshold [34, 35] (Figure 6). Different gray levels in some tracks may not be properly separated, complicating the extraction of track features. To address this issue, we
from scipy.
Figure 6: Input photomicrograph (Figure 4) binarized using the ISODATA threshold (threshold = 128) and region filling. Colormap: gray.
After binarizing the input image, the WUSEM algorithm separates the overlapping tracks (Figure 7). For this example, we chose initial radius = 10 and delta radius = 4 as parameters.
Tracks were counted using the “lower right corner” method, i.e., border tracks are only counted if they lie on the right or bottom edges of the image. We used the function clear rd border(), given in the Supplementary Material, to remove the tracks in these edges. WUSEM returns a labeled region, which can be used as parameter to the function enumerate objects() (Figure 8).
3.1.2. Comparison between manual and automatic counting
An experienced observer can easily distinguish tracks in the photomicrograph, even when several tracks are superimposed. For this reason, manual counting is considered the control in the comparison with the automatic counting. In the following, the WUSEM algorithm is applied to the photomicrograph set and the processing parameters are studied.
We established an arbitrary value of up to two tracks less than the mean of the manual counting as a tolerance. WUSEM’s parameters become candidates if the automatic counting lies within the tolerance interval, i.e. 0 < 2, where n is the track number obtained by each approach (Figure 9).
WUSEM’s best parameters for this study are the ones within the tolerance interval for most samples. According to the stated comparison, the best param-
eters are initial radius = 10, delta radius = 20 for 4.5 min samples, and
Using the best parameters defined, we compared manual, classic (flooding) watershed and WUSEM counting for each sample (Figure 10, Table 1). The classic watershed algorithm used is implemented in ImageJ [45], and the full method consisted in binarizing the input image, applying the watershed and removing small objects. The source code for these operations and instructions on how to use it are given in the Supplementary Material.
Since tracks in this dataset have the same shapes, we can attribute an ef-ficiency of 100 % for manual counting. WUSEM counting was initially set to obtain a smaller number of tracks when compared to the manual counting. However, WUSEM counting returns false positives, i.e., incorrectly labels background regions as ROI (points above the 1:1 line in Figure 10 (b) and (d)). To avoid false positives, one could use a more restrictive criterion, such as eccentricity (Section 3.1.3).
Counting reproducibility is important for considering the reliability of track counting. Despite the variations in track characteristics, the efficiencies of WUSEM’s automatic counting remained constant within uncertainties, when compared to manual counting.
Table 1: Mean, standard deviation, and automatic/manual efficiency ratio for each sample of the first dataset.
3.1.3. Relating ion energies with diameter product and mean gray levels
In this application, we use WUSEM to relate the track energies to the product of major and minor diameters (and
, respectively) and the mean gray level of each track. We considered only approximately circular tracks in our analysis, based on an eccentricity (
criterion:
should be equal to or less than 0.3 (Figure 11). This additional criterion can be used also for counting tracks, to ensure that false positives (spurious objects counted as tracks) are avoided.
After separating each track, we can obtain its features such as gray levels and diameters. Once the mean gray level and diameters of each track in a photomicrograph are obtained, we can relate them with the incident energy for each sample. The mean gray level of each sample is obtained getting the mean of all gray levels of the tracks in the images of that sample. We adopted a similar process to obtain the mean diameters. Then, the results are related to the incident energy (Figure 12).
The diameter products roughly reflect the electronic energy loss (dE/dx) curve for in DAP (Figure 12 (a) and (c)), calculated with the software SRIM [40]. The Bragg peak appears around 100 MeV. Further scatter of points can be attributed to poor control of etching conditions. The uncertainty of 3
in temperature may cause a large variation in etching results [46]. Variations in gray level means are impaired probably because this set was acquired in reflected illumination mode, which privileges surface details over depth effects.
3.2. Processing photomicrographs of fission tracks in DAP
In photomicrographs from the first dataset, our main concern was track superposition. However, all tracks were similar, created by a collimated beam of tracks. Here we go one step ahead, applying WUSEM to images where tracks present a variety of shapes. This image dataset was obtained from DAP plaques irradiated with thermal neutrons, coupled with apatite mounts. Fission fragments are born in the interior of the mineral, then emitted at different directions towards the detector. For instance, round tracks were created by perpendicular incident fragments, while the elliptical ones were created by particles hitting DAP surface at shallower angles (Figure 5).
As in the previous analysis, the manual counting performed by an experienced observer is taken as reference because we expect the observer to be able to recognize tracks efficiently. However, in this case, we do not expect the observer efficiency to 100 %. Fragments hitting the detector at lower energies originate tracks that are very difficult to distinguish from detector surfaces. An experienced observer would avoid counting those tracks to keep hers/his counting efficiency constant.
Repeating the previous processes for photomicrographs in the second dataset, we first binarize a test photomicrograph(Figure 5) using the ISODATA threshold. The binarized image is generated for two scenarios: considering and ignoring border tracks. Here, regions in the binary image are also filled using the func-
tion ndimage.morphology.binary fill holes() from scipy. Then we apply the WUSEM algorithm. We chose initial radius = 5 and delta radius = 2
as parameters, and the WUSEM result as parameter in the function enumerate objects()
(Figure 13).
For this dataset, we established an arbitrary tolerance of five tracks less than the mean of the manual counting. In this case, WUSEM’s parameters become candidates if 0 5, where n is the track number obtained by each approach. According to the stated comparison, the best param-
and initial radius = 10, delta radius = 14 for the second one. Using the
best parameters defined, we compared manual, classic watershed and WUSEM
counting for each sample (Figure 14, Table 2).
Table 2: Mean, standard deviation, and automatic/manual efficiency ratio for each magnifi- cation of the second dataset.
WUSEM succeeded in avoiding false positives in this application, when compared to the classic watershed (black line above the 1:1 line in Figure 14 (b)). Still, the user could apply more restrictive criteria. Also, it is worth noting the efficiency variation between the two image sets. Bigger objects are easier to be treated, thus automatic track counting in greater magnification images resulted in a higher counting efficiency (Table 2).
4.1. Structuring elements
In this study, we used disks as structuring elements for processing images in both datasets. Since tracks in photomicrographs from dataset 1 are almost isotropic (as seen in Figure 4), disks are suitable structuring elements to be used in their segmentation. However, tracks within images in dataset 2 do not have a defined format (Figure 5). Employing different structuring elements, e.g. rotated cones or ellipses, could improve their segmentation and the automatic counting result.
4.2. False positives, false negatives and counting efficiency
In track counting, reproducibility is not about counting every track in the image, but counting the same types of tracks every time. It is the primary concern in FTD: for instance, the efficiencies for counting tracks in the standard sample should not vary when using zeta age calibrations [47, 48]. For absolute methods of determining neutron fluence [49], the efficiency should be constant when counting tracks in unknown age samples. The fast-growing areas of FTD using the Laser Ablation Inductively Couple Plasma Mass Spectrometer (LA-ICP-MS, [50, 51, 52, 53]) or the Electron Microprobe ([54, 55]) have the same efficiency issues, and could also benefit from automatic counting.
The major challenge for reproducibility is avoiding false positives, spurious objects such as scratches on the detector or mineral surface and other etching figures which automatic counting algorithms could misrepresent as tracks. In most situations, it is preferable to restrict the criteria, thus increasing the number of false negatives (not counted tracks), even implying in efficiency reduction.
Even more experienced observers expect some decreasing in efficiency due to superposition when counting tracks in high track density samples. This loss tends to be more severe in automatic counting. When applying algorithms as WUSEM, the separation of tracks in objects formed of several tracks is not always possible.
Counting less tracks than the actual number in a cluster is acceptable; however, when processing a large number of clusters in higher track density samples, we expect a lower efficiency when compared with lower track density samples. This effect can be assessed by calibrating the efficiency as a function of track density. Therefore, efficiencies presented for WUSEM (Tables 1 and 2) only hold for the track densities of the used samples.
4.3. Perspective of future development
The WUSEM algorithm represents an advance in automated track counting, opening several possibilities. Differently from the direct application of classic watershed, WUSEM allows adaptation (Figure 15): the user can determine optimal structuring elements and efficiencies using a training image subset and apply these parameters to hers/his sample photomicrographs. This is very convenient in fission track dating, where the tracks present the same shapes regardless of the track source, and efficiency may vary mainly due to superposition of tracks in higher track density samples.
Another possibility is applying WUSEM in separate parts of the input image, allowing the determination of specific parameters for each track configuration. This could bring even better results when dealing with the superposition of two, three or more tracks.
In this paper we present a watershed algorithm using successive erosions as markers, which we call WUSEM. We employ WUSEM to separate overlapping tracks in photomicrographs. WUSEM performs well in images containing overlapping circular tracks (from a collimated beam) and in photomicrographs with fission tracks at various orientations, both in DAP. The results are encouraging: the mean automatic/manual counting efficiency ratio is 0
11 when using WUSEM in the test datasets. We show also that diameter and eccentricity criteria may be used to increase the reliability of this method.
Since WUSEM using circles as structuring elements is aimed to isotropically shaped regions, this technique is suitable for separating etched tracks in DAP. Etching velocity in mineral surfaces depends on the crystal orientation, yielding more complex etching figures. Also, natural minerals are richer in scratches and other etching figures that can be mistaken with fission tracks, especially when using image processing techniques. WUSEM could be studied to separate tracks in mineral surfaces; for that, it would need to use different structuring elements, which have to consider the orientation and shape of each track.
The authors would like to thank Raymond Jonckheere for the second photomicrograph dataset and Matthias Schr¨oter for his insights, and also Christina Trautmann for the sample irradiation at GSI in Darmstadt. This work is supported by the S˜ao Paulo Research Foundation (FAPESP), grants # 2014/22922-0 and 2015/24582-4.
Supplementary Material
The supplementary material contains the complementary results and the code published with this study, and is available at https://github.com/alexandrejaguar/
[10] K. Moriwaki, Y. Shimpuku, M. Makimura, Improvement of Track Count-
[11] A. J. W. Gleadow, S. J. Gleadow, D. X. Belton, B. P. Kohn, M. S.
[12] R. C. Gonzalez, R. E. Woods, Digital Image Processing, 3rd Edition, Pear-
[13] D. Azimi-Garakani, J. G. Williams, Automatic fission track counting using
[14] J. Beer, J. Pipper, W. Heinrich, Automatic measurements of nuclear tracks
[15] P. B. Price, W. Krischer, Semi-automated, three-dimensional measurement
[16] W. Trakowski, B. Sch¨ofer, J. Dreute, S. Sonntag, C. Brechtmann, J. Beer,
[17] S. Masumoto, K. Wadatsumi, Computerized Image-Processing System for
[18] K. Wadatsumi, S. Masumoto, Three-dimensional measurement of fission-
[19] A. P. Fews, Fully automated image analysis of etched tracks in CR-
[20] N. Petford, J. A. Miller, Three-dimensional imaging of fission tracks using
[21] N. Petford, J. A. Miller, J. Briggs, The automated counting of fission
[22] G. Espinosa, R. B. Gammage, K. E. Meyer, C. S. Dudney, Nuclear track
[23] A. Boukhair, A. Haessler, J. C. Adlo, A. Nourreddine, New code for digital
[24] M. Dolleiser, S. R. Hashemi-Nezhad, A fully automated optical microscope
[25] R. Bedogni, Development of a reader for the routine analysis of CR-39 fast
[26] N. Yasuda, K. Namiki, Y. Honma, Y. Umeshima, Y. Marumo, H. Ishii,
[27] H. Tawara, K. Eda, K. Takahashi, T. Doke, N. Hasebe, S. Kodaira,
[28] A. F. de Siqueira, W. M. Nakasuga, A. Pagamisse, C. A. Tello Saenz, A. E.
[29] J. Ruan, H. Li, J. Zhang, Z. Zhang, L. Chen, J. Song, J. Liu, J. Liu,
[30] E. Enkelmann, T. A. Ehlers, G. Buck, A. K. Schatz, Advantages and chal-
[31] S. Beucher, C. Lantuejoul, Use of Watersheds in Contour Detection (1979).
[32] J. Serra, Image analysis and mathematical morphology, Vol. 1, Academic
[33] F. M. Schaller, M. Neudecker, M. Saadatfar, G. Delaney, K. Mecke, G. E.
[34] G. H. Ball, D. J. Hall, A clustering technique for summarizing multivari-
[35] T. W. Ridler, S. Calvard, Picture Thresholding Using an Iterative Selection
[36] F. R. Dias Velasco, Thresholding Using the ISODATA Clustering Algo-
[37] S. van der Walt, J. L. Sch¨onberger, J. Nunez-Iglesias, F. Boulogne, J. D.
[38] C. Fiorio, J. Gustedt, Two linear time Union-Find strategies for image
[39] J. C. Russ, The Image Processing Handbook, 6th Edition, CRC Press, 2011.
[40] J. F. Ziegler, M. D. Ziegler, J. P. Biersack, SRIM - The stopping and range
[41] S. van der Walt, S. C. Colbert, G. Varoquaux, The NumPy Array: A
[42] E. Jones, T. Oliphant, P. Peterson, SciPy: Open source scientific tools for
[43] J. D. Hunter, Matplotlib: A 2D Graphics Environment, Computing in
[44] G. van Rossum, Python Reference Manual, Tech. rep., Centrum voor
[45] C. T. Rueden, J. Schindelin, M. C. Hiner, B. E. DeZonia, A. E. Wal-
[46] S. Guedes, J. C. Hadler, P. J. Iunes, S. R. Paulo, C. A. Tello, On the
[47] A. J. Hurford, P. F. Green, The zeta age calibration of fission-track dating,
[48] A. J. Hurford, Standardization of fission track dating calibration: Recom-
[49] C. J. Soares, I. Alencar, S. Guedes, R. H. Takizawa, B. Smilgys, J. C.
[50] J. C. Hadler, P. J. Iunes, C. A. Tello, F. Chemale, K. Kawashita, E. A. C.
[51] N. Hasebe, J. Barbarand, K. Jarvis, A. Carter, A. J. Hurford, Apatite
[52] C. J. Soares, S. Guedes, J. C. Hadler, R. Mertz-Kraus, T. Zack, P. J.
[53] C. J. Soares, R. Mertz-Kraus, S. Guedes, D. F. Stockli, T. Zack, Character-
[54] D. J. Gombosi, J. I. Garver, S. L. Baldwin, On the development of elec-
[55] A. N. C. Dias, F. Chemale, C. J. Soares, S. Guedes, A new approach
Figure 7: Processing an input image using WUSEM. Each line presents a segmentation step. Left column: original binary image (red) and erosion obtained according to the variation of delta radius (white). Center column: erosion results labeled, used as markers. Right column: watershed results when using generated markers. Parameters for WUSEM algorithm:
Figure 8: (a) Labels generated from tracks in Figure 4 using the WUSEM algorithm. (b) Tracks in (a) enumerated using enumerate objects(). Tracks in the lower or right corners are not counted, according to the “lower right corner” method. Parameters for WUSEM
algorithm: initial radius = 10, delta radius = 4. Colormaps: (a) nipy spectral, (b)
Figure 9: Manual counting mean (top of the blue bar; values on the right) for each sample and automatic counting results with mean within (orange points) and outside (gray points) the tolerance interval (blue bar) for the first dataset.
Figure 10: Comparison between manual and automatic counting for (a, b) 4.5 min etching samples and (c, d) 8.5 min etching samples. (a, c) white: manual counting. Gray: flooding watershed counting. Red line: distribution median. White signal: distribution mean. (b, d) dashed: 1:1 line. Red line: regression for the WUSEM counting data. Black line: regression for the flooding watershed counting data.
Figure 11: Regions from Figure 4 complying with labeled regions. (b) tracks correspondent to (a) in their original gray levels. Colormaps: (a) nipy spectral. (b) magma.
Figure 12: Relation between incident energy versus mean diameter product ((a) 4.5 min; (c) 8.5 min samples, left Y axis) and incident energy versus mean gray levels ((b) 4.5 min; (d) 8.5 min samples, left Y axis). Cyan dashed line: electronic energy loss calculated with SRIM (right Y axis).
Figure 13: (a) Input photomicrograph (Figure 5) binarized using the ISODATA threshold (threshold = 0.59) and region filling. (b) Tracks separated in Figure 5 using the WUSEM algorithm, and then enumerated using the function enumerate objects(). Parameters for
Figure 14: Comparison between manual and automatic counting for photomicrographs in dataset 2. (a) white: manual counting. Gray: flooding watershed counting. Red line: distribution median. White signal: distribution mean. Blue dots: outliers. (b) dashed: 1:1 line. Red line: regression for the WUSEM counting data. Black line: regression for the flooding watershed counting data.
Figure 15: When using suitable input parameters, WUSEM may perform better in certain regions where the classic watershed does not return reliable results. For instance, the highlighted region in (a) presents three tracks. WUSEM separates them correctly, but the region is oversegmented by the classic watershed. The highlighted region in (b), by its turn, is undersegmented by the classic watershed, which returns two tracks. WUSEM returns three tracks, being closer to the real number (four tracks). Left: input photomicrographs with highlighted regions. Center: tracks separated using WUSEM. Right: tracks separated using classic watershed. Parameters for WUSEM algorithm: initial radius = 15, delta radius