PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcbioiBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Bioinformatics
 
BMC Bioinformatics. 2010; 11: 373.
Published online 2010 July 8. doi:  10.1186/1471-2105-11-373
PMCID: PMC2914730

Extended morphological processing: a practical method for automatic spot detection of biological markers from microscopic images

Abstract

Background

A reliable extraction technique for resolving multiple spots in light or electron microscopic images is essential in investigations of the spatial distribution and dynamics of specific proteins inside cells and tissues. Currently, automatic spot extraction and characterization in complex microscopic images poses many challenges to conventional image processing methods.

Results

A new method to extract closely located, small target spots from biological images is proposed. This method starts with a simple but practical operation based on the extended morphological top-hat transformation to subtract an uneven background. The core of our novel approach is the following: first, the original image is rotated in an arbitrary direction and each rotated image is opened with a single straight line-segment structuring element. Second, the opened images are unified and then subtracted from the original image. To evaluate these procedures, model images of simulated spots with closely located targets were created and the efficacy of our method was compared to that of conventional morphological filtering methods. The results showed the better performance of our method. The spots of real microscope images can be quantified to confirm that the method is applicable in a given practice.

Conclusions

Our method achieved effective spot extraction under various image conditions, including aggregated target spots, poor signal-to-noise ratio, and large variations in the background intensity. Furthermore, it has no restrictions with respect to the shape of the extracted spots. The features of our method allow its broad application in biological and biomedical image information analysis.

Background

Biological imaging such as confocal fluorescence microscopy and electron microscopy require the use of protein-labeling techniques to localize individual proteins within cells. Biological markers such as green fluorescence protein [1] and a variety of fluorescent dyes [2,3] for fluorescence microscopy, and colloidal gold [4,5] for electron microscopy are widely used. Molecules labeled with biological markers are generally observed as small specific spots against a background of high brightness. Quantitative comprehension of the localization and statistical distribution of the spots are essential for deciphering biological information. In general, cellular microscopic images have a low signal-to-noise ratio (SNR) and the differences in intensity between signal spot and background are not always clear. Moreover, the texture of those backgrounds is complicated. For these reasons, microscopy images are often difficult to manage computationally. Currently, there are several automatic processing and recognition systems for biological images and they have been applied in the quantitative analysis of biological objects ranging from molecules to cells to whole organisms [6-10].

The purpose of this study was to extract and characterize biological spots of intricate morphology and low contrast in an automatic manner. Current standard techniques for spot extraction consist of edge enhancement for image morphology, including discrete convolution by a high-pass mask and the use of first- or second-order differential operators, based on the magnitude of the spatial differences of the spots [11]. One major problem with this approach, however, results from the blurring and degradation of the image contrast during image acquisition. For some spots with weak contrast, edge extraction is not sufficient. In real-world applications, most biological images contain object boundaries, artifacts, and noise. Therefore, edge enhancement filters may cause difficulties in distinguishing the exact edge of the object's structure from artifacts such as trivial geometric features. Additionally, these techniques can amplify background noise in the image while enhancing the object edge [12,13].

In other methods based on conventional frequency-selective filters [14-18], the precise localization of low-contrast spots may not be possible. High-density areas resulting from the integration of many spots may not allow the isolation of individual spots through frequency-selective filters. In addition, the parameter settings are often so complex as to require their modification whenever the target spot images are changed [19,20]. Furthermore, these methods cannot deal with the varied morphology of the spots.

Spot extraction methods based on conventional mathematical morphology [21] effectively capture the spots' location and their shape information [22-26]. These methods employ a morphological algorithm for background subtraction known as the top-hat transformation [27] or rolling-ball transformation [28]. It is well recognized that the principle of these methods is very effective for extracting a target object from a wide variety of image types [29-34].

Morphological operations use small synthetic images called structuring elements (SEs), which are a fundamental tool in mathematical morphology. The SE used as a probe moves along each pixel of the image. To apply morphological filtering for spot extraction from various types of biological images, the procedure to determine the shape and size of the SE is very important. A commonly used SE shape is the square or disk. In the rolling-ball transformation, a ball-shaped SE (such as a disk SE with weights arranged in order to describe a hemisphere in gray scales) is used. In the above-described methods for spot extraction, these SEs were also used. However, most small contiguous spots cannot be individually distinguished, such that several spots are extracted as one connected region because the size (width) of the SEs is wider than the minimum distance between the peaks of adjacent spots. A suitable SE shape for spot extraction includes a straight-line segment (a fuller description of which is given in the Methods and Results sections); however, since processing by common morphological operations with a single line-segment SE is not isotropic, it cannot consider the geometrical details of an intricate image. Thus, for spot extraction, conventional morphological processing is not effective.

Advanced morphological processing with multiple SEs has been reported [35,36]. In this approach, multiple sets of line-segment SEs generated by rotation of the single line-segment SE in different directions are applied. However, in the discrete space of the images, it is difficult to generate a straight line-segment as the SE that can be rotated in an arbitrary direction. This restriction in the rotational direction of the SE prevents adequate spot detection in complicated biomedical images.

In this study, we solved these problems by introducing a simple and practical approach, an extended mathematical morphology, into the automatic detection of spots in biological images. This technique is based on top-hat transformation with a single SE as the straight line-segment. In our algorithm, an original image can be rotated in arbitrary directions with respect to the single SE. This novel method, which we named rotational morphological processing (RMP), can homogeneously treat with geometrical features in an image under various orientations. Top-hat transformation based on RMP has been applied to spot extraction, in the absence of any hypotheses to fit the spots by 2-D Gaussian shape or minimal intensity. Finally, by isotropic processing with the line-segment SE, contiguous spots can be segmented into individual parts. Our novel method was developed in order to automatically extract spots, such as biological markers consisting of antibodies conjugated with fluorescent molecules, from a biological image of intricate morphology and low contrast.

Smal et al. evaluated the spot detection methods most frequently used in fluorescence microscopy [37], including wavelet-based multiscale detecting [16,18], morphological based methods [23,24,38], and the machine learning method [39]. In this study, we compared the performance of our proposed method with other morphological based methods, such as conventional top-hat transformation and h-dome transformation, by using synthetic-noise images.

This report is organized as follows. A brief introduction describing the basics of conventional mathematical morphology is followed, in the Methods section, by a detailed presentation of our spot extraction technique. In the Results section, the application of the detection method to synthetic images as well as to real image data from electron and fluorescence micrographs is discussed. In the final section, the effectiveness of our novel method is summarized and evaluated.

Conventional mathematical morphology

Mathematical morphology is based on set-theory concepts of the shape of an objective image [21]. An image can be represented by a set of pixels. Morphological operations always deal with a set of two images: an objective image and a SE. Each SE has shape and size characteristics as parameters of the operation. Let f denote a gray-scale image function from Z2 into [0, I-1], where I is a positive integer. Let B denote a binary SE. The fundamental operators of mathematical morphology are dilation and erosion.

dilation:

equation image
(1)

erosion:

equation image
(2)

where Df and Db are the domains of the functions f and B, respectively. The opening and closing operations are delivered from dilation and erosion.

opening:

equation image
(3)

closing:

equation image
(4)

The top-hat transformation is one of the commonly used morphological operations for extracting local bright objects from a low contrast image in gray-scale [27]. It is obtained by subtracting from the original image f the opening image γB using the SE B.

Top-hat:

equation image
(5)

It yields an image in which all the residual features (peaks and ridges) are subtracted by the opening operation. Adding these residual features to the original images has the effect of accentuating objective structures with high intensity [27]. If the difference in intensity between the target objects and the background of the image is markedly small, it is difficult to detect these differences with the human eye. However, these low-contrast objects can be extracted and enhanced by the top-hat transformation.

Another method to extract the local bright object in biological images, based on mathematical morphology, is the h-dome transformation [38].

h-dome:

equation image
(6)

where Dh(f) is the h-dome image of a gray-scale image f, (f-h) represents the result of subtracting a constant value h from the gray-scale image, and ρf(f-h) the morphological reconstruction of the gray-scale image from f-h. The gray-level reconstruction is obtained by iterative geodesic dilation of f-h under f until stability is reached [40].

Method

Spot extraction filter: top-hat transformation by RMP

The essential elements required for the spot extraction filter are processing of the biomedical image isotropically and isolation of the adjacent spots from the image background. In order to fulfill these requirements, our method proposes that the objective image can be rotated in arbitrary directions with respect to a single straight line-segment SE whose width has only 1 pixel. The width of this SE ensures determination of the minimum distance between two different spots for isolation. By using this SE, spots separated by distances of only 1 pixel can be distinguished individually. The length of the SE should be adjusted so that it is longer than the size of the target spot. A spot that is smaller than the length of this SE is extracted by the top-hat transformation.

This top-hat transformation by RMP with the straight line-segment SE consists of the following steps:

Algorithm 1 (Top-hat transformation by RMP with line-segment SE)

1. Original image rotation. The original image f (Figure (Figure1a)1a) is rotated in a clockwise direction with respect to the center of the image frame. Assume that dividing a half of the circle (π [rad]) into N equiangles gives us each direction at an angle of π/N [rad] (Figure (Figure1c),1c), which is an increment angle. Namely, fi (Figure (Figure1d,1d, top row) denotes the rotated image of f with the angle of π i/N [rad], where i = 0, 1,..., N-1.

Figure 1
Schematic procedure of rotational mathematical morphology. (a) Original image. (b) Straight line-segment SE B. (c) Setting of the rotation angles. (d) Process of RMP opening. Original image is rotated and each rotated image is processed by opening with ...

2. Opening. All rotated images are subjected to an opening operation with the straight line-segment SE B (Figure (Figure1b).1b). The opening operations of the rotated image fi are represented as γB(fi) (Figure (Figure1d,1d, middle row).

3. Opened image rotation. The opened images (γB(fi)) are rotated π·i/N [rad] in an anticlockwise direction. The rotation at i times of the opened image is denoted by hi (Figure (Figure1d,1d, bottom row).

4. Union of the rotated and opened images (opening by RMP). The processed images (hi) are unified. In union processing, the maximum intensity value, which corresponds to the same pixel coordinate among all opened images, is taken to generate the whole image.

5. Top-hat transformation by RMP. The unified opened image (γ'B(f)) is subtracted from the original image (f).

Figures Figures1a1a and and1b1b show the original image f and the straight line-segment SE B, respectively. Figure Figure1c1c illustrates a setting of the rotation angles. The process in the RMP opening with the straight line-segment SE is represented in Figure Figure1d.1d. RMP opening is defined as follows:

equation image
(7)

The image γ'B(f)(x, y) is defined as the maximum value of h0 (x, y), h1 (x, y),..., hN-1(x, y). The top-hat transformation by RMP is also given by the following equation:

equation image
(8)

This operation is used as the spot extraction filter in our proposed method.

In this study, we utilized the top-hat transform by RMP for spot extraction from biological images obtained with electron and fluorescence microscopy. The entire practical process consists of the following steps:

Algorithm 2 (Spot extraction for practical biological images)

1. Noise reduction: Noise, which is less than the resolution limit of the micrograph or target spot, is removed via opening by RMP (equation (7)) with the straight line-segment SE. The length of the SE is set to be smaller than the diameter of the target spots.

2. Spot extraction: The spots are extracted by the top-hat transformation by RMP (equation (8)) with the straight line-segment SE. The length of SE is set to be larger than the diameter of the target spots.

3. Binarization: The extracted spots are binarized by equation (9) for recognition and measurement of the spots computationally.

Binarization is performed according to the thresholding approach:

equation image
(9)

Namely, the pixels of residual regions by the top-hat transformation by RMP are assigned an intensity of 255 in an 8-bit gray-scale value.

In several cases, the subtraction process and the binarization process leave the small isolated pixels on the image such that they represent residual background noise. The conventional opening operation (equation (3)) can also be applied to remove the remaining noise. This post-processing should be adapted to the particular application.

Results

Isolation of overlapped spots

The performance of our proposed method in the isolation of overlapped spots was compared with that of conventional top-hat transformation. Figure Figure2a2a shows the model image of adjacent spots (f). On the left of Figure Figure2a2a is the original gray-scale image f and on the right is the 3-D topographic map of f. The vertical height of the map represents the intensity. In top-hat transformation, the original image f is first opened and then this opened surface is subtracted from the original surface. In this example, the straight line-segment SE BL (11 × 1 pixels) was used in our method, and the flat disk SE BD (the diameter was 11 pixels) was used in the conventional top-hat transformation. The difference between f and its opening by RMP with BL (i.e., top-hat transformation by RMP) is shown in Figure Figure2b.2b. The result of conventional top-hat transformation with BD of f is shown in Figure Figure2c.2c. In Figure Figure2d,2d, the black line is the profile of the spots image f (the position of the profile line is marked by the arrowhead in Figure Figure2a),2a), which denotes the surface of the spots. The red line is the result of the opening by RMP with BL of f (where the value of the rotational direction N is 36). The green line is the result of the conventional opening with BD. Figure Figure2b2b shows that the two overlapped spots were segmented clearly by our method.

Figure 2
Isolation of overlapped spots. (a) Original image of two overlapped spots. Left: original gray-scale image, right: the 3-D map. (b) The result of the top-hat transformation by RMP with the straight line-segment SE BL (11 × 1 pixels). (c) The result ...

The opening operation can be geometrically processed by pressing the SE up against the surface of the original image f and sliding it underneath the entire surface. The surface of the opened image is constructed from the highest point in the region reached by the SE. Since the line-segment SE BL has a width of 1 pixel, in the process of RMP opening the SE fits the narrow intervening space between the spots. Thus, line-segment SE reaches the baseline of the individual spots (i.e., the upper level of the overlapped region). In contrast, the disk SE BD has a width of 11 pixels as diameter. Since it is larger than the intervening space of the spots, the SE cannot fit the space. Thus, the disk SE cannot reach the level at which the two spots are distinguished. Accordingly, the spots were not isolated by the conventional top-hat transformation.

Verification of optimal number of rotational direction

The optimal number of the image rotational directions (N) was determined. For this experiment, we used a model image containing several artifacts. The image was reconstructed by RMP opening with the straight line-segment SE set at various values for the rotational direction. Peak signal-to-noise ratio (PSNR) was used to provide a quantitative evaluation of this performance. Assuming pixel values in the range of [0: 255], the PSNR is calculated using the following formula:

equation image
(10)

where F(i, j) denotes the original image, F'(i, j) the filtered image by RMP opening, and m1 × m2 the total number of pixels.

Figures Figures3a3a and and3b3b show the original 8-bit gray-scale image (256 × 256 pixels) and the model image overlapped by artifacts, respectively. Some of the round objects (diameter: 9 pixels) are regarded as artifacts to be removed. The RMP opening with the straight line-segment SE (31 × 1 pixels) of N = 8 allowed removal of the artifacts but also of some of the elongated target objects (Figure (Figure3c).3c). At N = 36, artifacts were removed while the other objects were well preserved (Figure (Figure3d).3d). The graph in Figure Figure3e3e shows the variation of PSNR with N, which ranged from 0 to 90 directions. The range of small values of N (from 0 to 8) is enlarged in the insertion. The black line denotes the result of this operation. The highest value of PSNR was obtained at N = 36 and it subsequently remained at almost the same level. In terms of the performance of reconstruction and the computational cost of processing, it is reasonable to assume that the optimal number of image rotational directions is 36.

Figure 3
Verification of the optimal number of rotational directions. (a) Original image (256 × 256 pixels). (b) Artifact-contaminated image. Some round objects (diameter: 9 pixels) are denoted as artifacts. The image was restored by RMP opening with the ...

In addition to this experiment, performance was also tested using the multiple SEs method, shown as a red line in Figure Figure3e.3e. The straight line-segment SE was used for these tests. In the case of N = 1, the single SE (orientation θ = 0 [rad], horizontal direction) was applied; at N = 2, two SEs (θ = 0 and π/2 [rad]) were used, and at N = 4, four SEs (θ = 0, π/2, π/4 and 3π/4 [rad]) were used. These results showed that reconstruction by the multiple SEs was insufficient.

Noise reduction

The noise reduction technique in our method (in algorithm 2), i.e., RMP opening with line-segment SE, was used. Figure Figure4a4a shows the spot model with background noises. These noises were removed by RMP opening (Figure (Figure4b).4b). The length of the line-segment SE was determined from the spatial resolution of the microscopic image, with objects having a shape size smaller than the resolution limit regarded as noise. If the spot size is known in advance, the SE length should be set to remove objects smaller than the spot size.

Figure 4
Noise reduction through RMP opening. (a) Spot model image with background noises. (b) Noises removed from the image by RMP opening with the straight line-segment SE.

Comparison with conventional morphological methods by using synthetic-noise images

The ability of the proposed method was compared with that of conventional morphological methods, such as the conventional top-hat and h-dome transformations. Synthetic-noise images were generated to evaluate the quantitative performance of spot detection. Ten source images of size 512 × 512 pixels were generated, with one such example shown in Figure Figure5a.5a. Each image contained 15 spots, which were modeled as a 2-D Gaussian distribution in random locations. This distribution is a good approximation of the theoretical shape of the spot [41]. Two types of spot models with different intensity levels (255 or 127 in an 8-bit gray-scale) were generated (Figure (Figure5b).5b). First, a 2-D Gaussian distribution with kernel width = 17 (pixels) and standard deviation σ = 2 (pixels) was generated. This distribution was then expanded to the normalized gray-scale range of 0-255. For the second spot model, the distribution was expanded to the range of 0-127.

Figure 5
Experimental results using synthetic images. (a) Example of source image with an 8-bit gray-scale. (b) 2-D Gaussian distributions as spot models. Two spot models with different intensity levels were used in this experiment. The actual image of the spot ...

The Poisson-distributed noise was added to each source image with a uniform (type-A) and a gradient (type-B) background (Figure (Figure5c).5c). This is one of the main sources of noise in fluorescence microscopy imaging [42].

The PSNR between each source image and the corresponding noise-added image was calculated. The average value of the PSNR of the type-A image set was 13.848 ± 0.242 (mean ± SD) dB, while for the type-B image set it was 11.245 ± 0.169 dB.

Since the spots in the synthetic image are prospectively known, the recall, precision, and F-measure rates can be calculated. Let TP denote the number of spots detected correctly, FN the number of spots missing detection, and FP the number of spots false alarmed. If adjacent spots are not separated, those spots are regarded as FN. The recall and precision rates are defined by the equations TP/(TP+FN) and TP/(TP+FP), respectively. The F-measure rate is the harmonic mean recall and precision rate and takes into account both measures:

equation image
(11)

In each image type, the synthetic images had a total of 150 spots; thus, TP + FN = 150. The results detected by our proposed method were compared with those of the conventional top-hat and h-dome transformations.

In this experiment, noise-added images were first smoothed with the Gaussian filter (3 × 3 kernel).

Processing by our proposed method followed algorithm 2. Since the width of the spot domain was 17 pixels (Figure (Figure5b),5b), any structure smaller than this width was regarded as noise or artifact. In step 1, the straight line-segment SE of 13 × 1 pixels was used, and in step 2, the straight line-segment SE of 21 × 1 pixels. In both steps, the rotational direction number (N) of RMP processing was 36. The subtracted image was binarized by the method in step 3.

For the method based on the conventional top-hat transformation, the processing in step 2 differs from that of the proposed method. Top-hat transformation was applied with disk SE (diameter: 21 pixels) and the subtracted image was binarized by the method in step 3.

For the method based on h-dome transformation, this transformation was applied to the smoothed image obtained in step1 of algorithm 2. We set the value of parameter h to 50, and the subtracted image (h-dome) was binarized by the method in step 3.

Finally, binary images obtained from these methods were cleaned by the conventional opening with disk SE (diameter: 13 pixels) as post-processing.

The results of the spot extraction are seen in Figure Figure5c.5c. The type-A and type-B images are shown in the top and bottom row, respectively. All spots were correctly detected by our proposed method but not by the other methods. Figure Figure5d5d shows the actual values of recall, precision, and F-measure for each type image set. The performance in terms of the F-measure for the proposed method was consistently 100% in all type images. Thus, the proposed method is more tolerant of Poisson noise images than other methods.

In addition, the performance of our method was tested using synthetic images under various noise levels. For the source image (Figure (Figure5a),5a), nine degrees of Poisson-distributed noise images with uniform background were generated, with the PSNR decreasing from 17.732 to 7.797 dB. The upper part of Figure Figure5e5e shows the synthetic images, captured only in the rectangular region in Figure Figure5a.5a. The proposed method was applied to these synthetic images with the same procedure and SEs as in the previous experiment. The bottom of Figure Figure5e5e shows that the scores for the recall, precision, and F-measure rates were consistently 100% in the PSNR range of 17.732 to 10.057 dB. Subsequently, the recall rate decreased with decreasing PSNR although the precision rate remained at 100%. This indicates that FP was zero and ensures the accuracy of our method for spot detection.

Spot extraction of colloidal gold particles

The proposed method was applied to electron microscopic images to evaluate its performance in spot extraction. Electron micrographs containing colloidal gold particles were used as test images, with two different sizes of particles (10 nm and 1.8 nm in diameter). Figure Figure6a6a (left) shows a part of the original micrograph with 10-nm-diameter gold particles (British BioCell). The spatial resolution was 0.90 nm/pixel. Figure Figure6b6b (left) shows a part of the original micrograph with 1.8-nm-diameter gold particles (Nickel (II)-Nitrilotriaceticacid-Nanogold, Nano Probes). The spatial resolution of this image was 0.32 nm/pixel. These images were on an 8-bit gray-scale, the intensity values of the routine image quality were inverted. The purpose of this test was to verify whether the diameter of the extracted gold-particle spot as determined by our method was consistent with the nominal diameter.

Figure 6
Spot extraction of colloidal gold particles from electron microscopy images. (a) Original micrograph of 10-nm-diameter gold particles (left). Spot-extracted image (right). The extracted spots are denoted in red regions. Bar: 200 nm. (b) Original micrograph ...

Processing of the spot extraction was carried out with algorithm 2. In the case of the micrograph containing 10-nm gold particles, line-segment SE with a size of 5 × 1 pixels (4.5 × 0.9 nm) and 13 × 1 pixels (11.7 × 0.9 nm) was used in noise reduction and spot extraction, respectively. The extracted spots were binarized and overlaid as seen in the red region on the original image. The result is shown on the right side of Figure Figure6a.6a. For the electron micrograph containing the 1.8-nm gold particles, line-segment SE with a size of 3 × 1 pixels (0.96 × 0.32 nm) and 7 × 1 pixels (2.24 × 0.32 nm) was used for noise reduction and spot extraction, respectively. The result is shown on the right side of Figure Figure6b.6b. For all extracted particles, the averaged Feret's diameter was calculated. The mean ± SD of the diameter of the extracted spots was 10.16 ± 0.77 nm (164 spots) for the 10-nm gold particles, and 1.81 ± 0.15 nm (813 spots) for the 1.8-nm gold particles. This result shows that the diameter of the extracted spots was consistent with the nominal diameter.

Spot extraction of fluorescent antibodies

Our proposed spot extraction method was applied to fluorescence microscopy images in which caveolin-1 molecules in fibroblasts were stained with secondary fluorescent antibodies (Figure (Figure7a).7a). This image had an 8-bit gray-scale and a spatial resolution of 60 nm/pixel. Line-segment SE with a size of 3 × 1 pixels (180 × 60 nm) and 7 × 1 pixels (420 × 60 nm) was used for noise reduction and spot extraction, respectively.

Figure 7
Quantitative spot extraction of fluorescent antibodies from a fluorescence microscopy image. (a) Original image. Fluorescence micrograph of caveolin-1 molecules stained with the primary antibody and a secondary fluorescent antibody in normal fibroblasts. ...

By visual observation, the diameter of each spot was found to be about 5 pixels (ca. 300 nm). Therefore, a SE length longer than the diameter of the target spots was selected for the spot extraction process. Binarization was carried out using equation (9). The regions of the extracted spots were superimposed on the original image as red-colored regions (Figure (Figure7b).7b). From this micrograph, the 627 spots were extracted. The quantitative estimation of the bright spots corresponding to the caveolae is provided in Figure Figure7c.7c. The sum of the intensity of the pixel values in each extracted spot region was calculated and the distribution depicted in a histogram. The total number of extracted spots was 6701 from five micrographs (which included approximately 5 cell regions). The median of the histogram was 1319. The largest cluster can be seen centered at the histogram's median value.

As shown in Figure Figure8,8, the proposed method was applied to the various shapes of the fluorescent spots (top row, original images): (a) Gaussian-like shape (low peak height), (b) Gaussian-like shape (high peak height), (c) irregular shape, and (d) volcano shape. The 3-D maps of the original images are shown in the middle row, which illustrates the relief of the pixel surface. The 3-D maps of the extracted spots, as determined by our method, are shown in the bottom row. All spots were well extracted even though they had intricate shapes.

Figure 8
Precise extraction of various shaped spots. Examples of the four shapes of fluorescence spots from micrograph are shown in the top row. (a) Gaussian-like shape (low peak height), (b) Gaussian-like shape (high peak height), (c) irregular shape, and (d) ...

Discussion and conclusions

Our novel method to extract the spots in electron and fluorescence microscopic images uses the extended morphological filter through the top-hat transformation by RMP. We have successfully shown that the method is useful for extracting spots in biomedical images in which the conventional method is inadequate. The key concepts of our spot extraction method are the use of a straight line-segment SE and the rotation of the original image. By changing the length of SE, target spots of various sizes can be extracted. The method avoids the technical difficulties of traditional morphological processing and its performance is robust in the processing of biomedical images. The main advantages of our method are that it is computationally simple and easily modified for the extraction of target spots of different sizes and shapes, and that it can handle images in various conditions, e.g., aggregated target spots, poor SNR, and a background with large variations in intensity. The method yields directional information regarding the spatial distribution of spots within the cell as well as the frequency distribution of the size and intensity of the spots.

Our method is based on a line-segment SE with a 1-pixel width as the minimum separation distance and therefore allowed two or more target spots located close to each other to be clearly distinguished (Figure (Figure2).2). With conventional morphological top-hat transformation using the common SE shape (such as a disk or square), it is difficult to separate such spots. A similar difficulty arises when the "ball" SE is used. Since it has a radius that is larger than the inter-space distance between adjacent spots, it cannot fit within the space. The top surface obtained during opening with the rolling ball cannot reach the baseline allowing for separation of the spots.

To verify the optimality of the number of rotational directions (N) shown in Figure Figure3,3, we investigated how an artifacts-contaminated image (Figure (Figure3b)3b) could be restored by the RMP opening with increasing N. The experimental result (Figure (Figure3e)3e) showed that N = 36 was a better trade-off because the value of PSNR was low for N < 36 while for N > 36 the processing time became longer. In our method, a large computational cost, which is proportional to the size of the input images, is inevitable.

We compared our spot extraction method with the conventional top-hat and h-dome transformations. As seen in Figure Figure5,5, our method outperformed the others with respect to the three criteria (Figure (Figure5d).5d). For the proposed method, the performance in terms of F-measure rate was maintained at 100% among all background typed images. The precision rate of the conventional top-hat transformation was much lower due to its higher FP value (124 in type-A and 228 in type-B, respectively) in detection of the noise. Furthermore, conventional top-hat transformation could not separate adjacent spots, as it used the disk SE. Meanwhile, the recall rate of the h-dome transformation was much lower due to its higher FN value (113 in type-A and 28 in type-B, respectively). Thus, the number of undetected true spots was large.

We further investigated the change in the three measurements as a function of decreasing PSNR from 17.732 to 7.797 dB (Figure (Figure5e).5e). The F-measure rate was maintained at 100% until PSNR decreased to about 10 dB. Subsequently, when PSNR decreased further, the F-measure rate decreased as well due to a decreasing recall rate (thus, increasing the FN value); however, the precision rate was constantly 100%. These results showed that our method is accurate in spot detection.

In the measurement of gold particles in the electron micrograph (Figure (Figure6),6), the value of the averaged Feret's diameter of the extracted spots and the value of the nominal diameter of the gold particles were in close agreement. Thus, our method effectively extracted spots of the specified size with high accuracy.

Figure Figure77 shows the location of the small spots in the cell and the estimation of the spots intensities. Previously, Orlichenko reported that stimulation of epithelial cells with epithelial growth factor (EGF) resulted in a profound increase in the number of caveolar structures at the plasma membrane [43]. Our method was able to carry out precise quantitative measurements of the spatial and intensity distributions of the membrane domain with respect to caveolae.

Furthermore, our method allows effective extraction of various shaped spots. Since it is based on the top-hat transformation, the spots are extracted independently in terms of the shape of the surface relief, which is based on variations in the intensity value within a spot region (Figure (Figure8).8). In the conventional spot detection methods that rely on matched filtering, a 2-D Gaussian distribution is commonly used for the matched filter, assuming that a point-spread function of a signal spot has a 2-D Gaussian distribution. However, because most spots have an irregular topology, as in the example in Figure Figure8,8, accurate spot extraction is difficult using the matched filtering method.

Signal spots extracted by our method can be transformed into a 2-D Gaussian distribution as a normalization of spot shape. This allows the application of our method to the conventional automatic tracking system of individual fluorescent particles [44].

The RMP-based method enables a shape and intensity analysis for various types of biomedical images (2-D gel electrophoresis image, DNA microarray image, electron micrograph, X-ray mammographic image, etc.). It can be applied not only to spot extraction but also to a wide variety of important image processing techniques, such as segmentation, smoothing, and pattern extraction [45]. Overall, it provides a wide-ranging analytical approach to biological and biomedical informatics.

Authors' contributions

YK conceived the study, developed the algorithms, carried out the testing and fine-tuning of the algorithms using the programming language C/C + +, and wrote the first draft of the paper. NB provided useful comments on methodology and helped to revise the manuscript. NM carried out all of the LM, TEM, and molecular biology experiments, supervised the work, and edited and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Dr. S. Yuasa (National Center of Neurology and Psychiatry) for helpful discussions and encouragement, Prof. T. Kodama (Osaka University) for continuing guidance and encouragement. We also express our appreciation to Prof. E. Katayama (University of Tokyo) and Dr. K. Aoyama (FEI Japan). This work was supported in part by a Health Labor Science Research Grant (Nano-001) and Grants-in-Aid for Scientific Research from the Ministry of Education Culture, Sports, Science and Technology to N. Morone.

References

  • Shimomura O, Johnson FH, Saiga Y. Extraction, purification and properties of aequorin, a bioluminescent protein from the luminous hydromedusan, Aequorea. J Cell Comp Physiol. 1962;59:223–239. doi: 10.1002/jcp.1030590302. [PubMed] [Cross Ref]
  • Paddock S. Over the rainbow: 25 years of confocal imaging. Biotechniques. 2008;44(5):643–644. doi: 10.2144/000112798. [PubMed] [Cross Ref]
  • Giepmans BNG, Adams SR, Ellisman MH, Tsien RY. The fluorescent toolbox for assessing protein location and function. Science. 2006;312:217–224. doi: 10.1126/science.1124618. [PubMed] [Cross Ref]
  • Handley DA. In: Colloidal Gold: Principles, Methods, and Applications. Hayat MA, editor. Vol. 1. San Diego: Academic Press; 1989. The development and application of colloidal gold as a microscopic probe; pp. 1–32.
  • Hainfeld JF, Powell RD. New frontiers in gold labeling. J Histochem Cytochem. 2000;48:471–480. [PubMed]
  • Peng H. Bioimage informatics: a new area of engineering biology. Bioinformatics. 2008;24:1827–1836. doi: 10.1093/bioinformatics/btn346. [PubMed] [Cross Ref]
  • Bozinov D, Rahnenfuhrer J. Unsupervised technique for robust target separation and analysis of DNA microarray spots through adaptive pixel clustering. Bioinformatics. 2002;18:747–756. doi: 10.1093/bioinformatics/18.5.747. [PubMed] [Cross Ref]
  • Rahnenfuhrer J, Bozinov D. Hybrid clustering for microarray image analysis combining intensity and shape features. BMC Bioinformatics. 2004;5:47. doi: 10.1186/1471-2105-5-47. [PMC free article] [PubMed] [Cross Ref]
  • Mete M, Hennings L, Spencer HJ, Topaloglu U. Automatic identification of angiogenesis in double stained images of liver tissue. BMC Bioinformatics. 2009;10(Suppl 11):S13. doi: 10.1186/1471-2105-10-S11-S13. [PMC free article] [PubMed] [Cross Ref]
  • Berth M, Moser FM, Kolbe M, Bernhardt J. The state of the art in the analysis of two-dimensional gel electrophoresis images. Appl Microbiol Biotechnol. 2007;76:1223–1243. doi: 10.1007/s00253-007-1128-0. [PMC free article] [PubMed] [Cross Ref]
  • Gonzalez RC, Woods RE. Digital Image Processing. 3. New Jersey: Prentice Hall; 2008.
  • Suri JS, Wilson DL, Laxminarayan S, Eds. Handbook of Biomedical Image Analysis: Volume 2: Segmentation Models Part B. New York: Springer-Verlag; 2005.
  • Bankman IN, Ed. Handbook of Medical Imaging: Processing and Analysis. San Diego: Academic Press; 2000. [PubMed]
  • van der Heijden F, Apperloo W, Spreeuwers LJ. Numerical optimisation in spot detector design. Pattern recognition letters. 1997;18:1091–1097. doi: 10.1016/S0167-8655(97)00086-X. [Cross Ref]
  • Blanford RP, Tanimoto SL. Bright-spot detection in pyramids. Comput Vis Graph Image Process. 1988;43:133–149. doi: 10.1016/0734-189X(88)90058-8. [Cross Ref]
  • Olivo-Marin J-C. Extraction of spots in biological images using multiscale products. Pattern Recognition. 2002;35:1989–1996. doi: 10.1016/S0031-3203(01)00127-3. [Cross Ref]
  • Genovesio A, Liedl T, Emiliani V, Parak WJ, Coppey-Moisan M, Olivo-Marin J-C. Multiple particle tracking in 3-D + t microscopy: Method and application to the tracking of endocytosed quantum dots. IEEE Trans Image Process. 2006;15:1062–1070. doi: 10.1109/TIP.2006.872323. [PubMed] [Cross Ref]
  • Zhang B, Fadili J, Starck J-L, Olivo-Marin J-C. Multiscale variance-stabilizing transform for mixed-Poisson-Gaussian processes and its applications in bioimaging. Proceedings of the IEEE International Conference on Image Processsing. 2007;6:VI-233–VI-236.
  • Zeng L, Wu W. Motion Objects Detection Based on Wavelet Clustering. Proceedings of 2nd IEEE International Conference on Computer Science and Information Technology. 2009. pp. 562–566. full_text.
  • Bao W, Zhou R, Yang J, Yu D, Li N. Anti-aliasing lifting scheme for mechanical vibration fault feature extraction. Mech Syst Signal Process. 2009;23:1458–1473. doi: 10.1016/j.ymssp.2009.02.010. [Cross Ref]
  • Serra J. Image analysis and mathematical morphology. London: Academic Press; 1982.
  • Meyer F. Iterative image transformations for an automatic screening of cervical smears. J Histochem Cytochem. 1979;27:128–135. [PubMed]
  • Bright DS, Steel EB. Two-dimensional top hat filter for extracting spots and spheres from digital images. J Microsc. 1987;146:191–200.
  • Breen EJ, Joss GH, Williams KL. Locating objects of interest within biological images: The top hat box filter. J Comput-Assist Microsc. 1991;3:97–102.
  • Angulo J. Mathematical morphology operators for reading radioactivity DNA array images. Proceedings of the IASTED International Conference on Visualization, Imaging, and Image Processing. 2004. pp. 802–807.
  • Janson ME, Setty TG, Paoletti A, Tran PT. Efficient formation of bipolar microtubule bundles requires microtubule-bound γ-tubulin complexes. J Cell Biol. 2005;169:297–308. doi: 10.1083/jcb.200410119. [PMC free article] [PubMed] [Cross Ref]
  • Meyer F. In: Quantitative analysis of microstructures in material sciences, biology and medicine. Chermant JL, editor. Stuttgart: Rieder Verlag; 1977. Contrast features extraction; pp. 374–380.
  • Sternberg SR. Grayscale morphology. Comput Vis Graph Image Process. 1986;35:333–355. doi: 10.1016/0734-189X(86)90004-6. [Cross Ref]
  • Manders EM, Hoebe R, Strackee J, Vossepoel AM, Aten JA. Largest contour segmentation: a tool for the localization of spots in confocal images. Cytometry. 1996;23:15–21. doi: 10.1002/(SICI)1097-0320(19960101)23:1<15::AID-CYTO3>3.0.CO;2-L. [PubMed] [Cross Ref]
  • Giakoumis I, Nikolaidis N, Pitas I. Digital image processing techniques for the detection and removal of cracks in digitized paintings. IEEE Trans Image Process. 2006;15:178–188. doi: 10.1109/TIP.2005.860311. [PubMed] [Cross Ref]
  • Yang Y, Huang S, Rao N. An automatic hybrid method for retinal blood vessel extraction. Int J Appl Math Comput Sci. 2008;18:399–407. doi: 10.2478/v10006-008-0036-5. [Cross Ref]
  • Kapsalas P, Maravelaki-Kalaitzaki P, Zervakis M, Delegou ET, Moropoulou A. A morphological fusion algorithm for optical detection and quantification of decay patterns on stone surfaces. Constr Building Mater. 2008;22:228–238. doi: 10.1016/j.conbuildmat.2006.08.024. [Cross Ref]
  • Kim H, Morgan DE, Zeng H, Grizzle WE, Warram JM, Stockard CR, Wang D, Zinn KR. Breast tumor xenografts: diffusion-weighted MR imaging to assess early therapy with novel apoptosis-inducing anti-DR5 antibody. Radiology. 2008;248:844–851. doi: 10.1148/radiol.2483071740. [PubMed] [Cross Ref]
  • Ji Q, Engel J, Craine E. Texture analysis for classification of cervix lesions. IEEE Trans Med Imaging. 2000;19:1144–1149. doi: 10.1109/42.896790. [PubMed] [Cross Ref]
  • Stevenson RL, Arce GR. Morphological Filters: Statistical and Further Syntactic Properties. IEEE Trans Circuits and Systems. 1987;34:1292–1305. doi: 10.1109/TCS.1987.1086067. [Cross Ref]
  • Song J, Delp EJ. The analysis of morphological filters with multiple structuring elements. Comput Vis Graph Image Process. 1990;50:308–328. doi: 10.1016/0734-189X(90)90150-T. [Cross Ref]
  • Smal I, Loog M, Niessen W, Meijering E. Quantitative comparison of spot detection methods in fluorescence microscopy. IEEE Trans Med Imaging. 2010;29:282–301. doi: 10.1109/TMI.2009.2025127. [PubMed] [Cross Ref]
  • Vincent L. Morphological grayscale reconstruction in image analysis: Applications and efficient algorithms. IEEE Trans Image Process. 1993;2:176–201. doi: 10.1109/83.217222. [PubMed] [Cross Ref]
  • Jiang S, Zhou X, Kirchhausen T, Wong STC. Detection of molecular particles in live cells via machine learning. Cytometry A. 2007;71:563–575. [PubMed]
  • Soille P. Morphological Image Analysis: Principles and Applications. Berlin: Springer-Verlag; 1999.
  • Lacoste TD, Michalet X, Pinaud F, Chemla DS, Alivisatos AP, Weiss S. Ultrahigh-resolution multicolor colocalization of single fluorescent probes. Proc Natl Acad Sci USA. 2000;97:9461–9466. doi: 10.1073/pnas.170286097. [PubMed] [Cross Ref]
  • Smal I, Draegestein K, Galjart N, Niessen W, Meijering E. Particle filtering for multiple object tracking in dynamic fluorescence microscopy images: application to microtubule growth analysis. IEEE Trans Med Imaging. 2008;27:789–804. doi: 10.1109/TMI.2008.916964. [PubMed] [Cross Ref]
  • Orlichenko L, Huang B, Krueger E, McNiven MA. Epithelial growth factor-induced phosphorylation of caveolin 1 at tyrosine 14 stimulates caveolae formation in epithelial cells. J Biol Chem. 2006;281:4570–4579. doi: 10.1074/jbc.M512088200. [PubMed] [Cross Ref]
  • Cheezum MK, Walker WF, Guilford WH. Quantitative comparison of algorithms for tracking single fluorescent particles. Biophys J. 2001;81:2378–2388. doi: 10.1016/S0006-3495(01)75884-5. [PubMed] [Cross Ref]
  • Kimori Y, Oguchi Y, Ichise N, Baba N, Katayama E. A procedure to analyze surface profiles of the protein molecules visualized by quick-freeze deep-etch replica electron microscopy. Ultramicroscopy. 2007;107:25–39. doi: 10.1016/j.ultramic.2006.04.012. [PubMed] [Cross Ref]

Articles from BMC Bioinformatics are provided here courtesy of BioMed Central