Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Genomics**|**v.11(Suppl 2); 2010**|**PMC2975414

Formats

Article sections

- Abstract
- Background
- Methods
- Results
- Discussion
- Conclusions
- Competing interests
- Authors' contributions
- References

Authors

Related links

BMC Genomics. 2010; 11(Suppl 2): S9.

Published online 2010 November 2. doi: 10.1186/1471-2164-11-S2-S9

PMCID: PMC2975414

Jinshan Tang: ude.nrocla@gnatj; Shengwen Guo: moc.liamg@oug.newgnehs; Qingling Sun: moc.oohay@nusgnilgniq; Youping Deng: ude.msu@gned.gnipuoy

The 2009 International Conference on Bioinformatics & Computational Biology (BioComp 2009)

Hamid R Arabnia, Mary Qu Yang, Ruzena Bajcsy, Andrzej Niemierko, Brian D Athey, Weizhong Li and A Keith Dunker

http://www.biomedcentral.com/content/pdf/1471-2164-11-S2-info.pdfThe 2009 International Conference on Bioinformatics & Computational Biology (BioComp 2009)

13–16 July 2009

Las Vegas, NV, USA

Copyright ©2010 Tang et al; licensee 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.

This article has been cited by other articles in PMC.

Ultrasound imaging technology has wide applications in cattle reproduction and has been used to monitor individual follicles and determine the patterns of follicular development. However, the speckles in ultrasound images affect the post-processing, such as follicle segmentation and finally affect the measurement of the follicles. In order to reduce the effect of speckles, a bilateral filter is developed in this paper.

We develop a new bilateral filter for speckle reduction in ultrasound images for follicle segmentation and measurement. Different from the previous bilateral filters, the proposed bilateral filter uses normalized difference in the computation of the Gaussian intensity difference. We also present the results of follicle segmentation after speckle reduction. Experimental results on both synthetic images and real ultrasound images demonstrate the effectiveness of the proposed filter.

Compared with the previous bilateral filters, the proposed bilateral filter can reduce speckles in both high-intensity regions and low intensity regions in ultrasound images. The segmentation of the follicles in the speckle reduced images by the proposed method has higher performance than the segmentation in the original ultrasound image, and the images filtered by Gaussian filter, the conventional bilateral filter respectively.

Ultrasound imaging technology has wide applications in cattle reproduction and has been used to monitor individual follicles and determine the patterns of follicular development [1-6]. The adoption of ultrasound imaging technology in cattle reproduction can provide an effective way to understand a number of issues on bovine reproductive cycle and its concurrent disorders [4]. For example, with the help of ultrasound imaging technology, it is now known that follicular growth occurs in wave-like patterns during each estrous cycle [7]. Ultrasound imaging technology also provides a tool for understanding the influence of dominant follicles on medium and small follicles [7].

In the applications of ultrasound imaging to monitoring individual follicles and determining the patterns of follicular development, the acquisition of the measurements of the individual follicles such as diameters, areas and perimeters is very important. In order to acquire the measurements of an individual follicle, image segmentation techniques are often used to extract the individual follicles. However, speckles in ultrasound images affect the segmentation and finally affect the measurement of the follicles. Speckle noise, seen as a granular structure, is caused by the interaction between the ultrasound waves and the scatters within the tissue [8]. The inherent nature of speckles makes its removal difficult. Speckle noise is not an additive noise, but is considered as a kind of multiplicative noise [9][10]. Many speckle reduction technologies have been proposed. In [11], a Laplacian pyramid-based nonlinear diffusion (LPND) is presented for medical ultrasound imaging. In the proposed method, the image is first decomposed into multi-layer Laplacian pyramid and speckles are removed by nonlinear diffusion filtering of bandpass ultrasound images in Laplacian pyramid domain. In [12], a nonlinear multiscale wavelet diffusion for speckle reduction is proposed. Speckles are suppressed by employing the iterative multiscale diffusion on the wavelet coefficients. In [9], a speckle reduction algorithm---speckle reducing anisotropic diffusion (SRAD) is proposed. The proposed algorithm has good performance in the preservation of edges and speckle reduction.

In this paper, we will investigate using bilateral filter to reduce the speckles in ultrasound images for cattle follicle segmentation. It is well known that bilateral filter has good performance in noise reduction and edge preservation. However, current existing bilateral filters are mainly used for additive noise reduction. It is not effective when it is applied to speckles, which are generally modelled as multiplicative noise. In order to solve this issue, we propose an adaptive bilateral filter, which can reduce the speckles effectively.

Bilateral filter was developed by Tomasi and Manduchi [13]. The basic idea of bilateral filter is to replace a pixel value in an image by a weighted mean of its neighbors, which the weights depend on both the spatial distance and the intensity distance [14][15]. There are many types of bilateral filters depending on the choice of weighting functions. What we develop in this paper is based on the Gaussian bilateral filter [13][16]. For Gaussian bilateral filter, it can be expressed mathematically as [13][17]

(1)

where is the output pixel value, *J* (*Y*) is the input pixel values, *X* and *Y* are the coordinates vectors,* σ
_{d}*

(2)

In the above equation, when *X* and *Y* are 2-D vectors, the bilateral filter is called 2-D bilateral filter, which can be used to reduce the noise in 2-D images.

Bilateral filter is a good choice for image de-noising because it is stable and simple. The effectiveness of bilateral filter lies in the combined use of the domain filter, which is used to enforce spatial closeness by weighting pixel values with coefficients that fall off with distance [18], and the range filter, which assigns greater coefficients to those neighbouring pixels with light intensity that is more similar to the centre pixel value [18]. In bilateral filter, the choice of the parameters *σ
_{d}*

Generally speaking, noise can be modelled by an additive model or a multiplicative model. Additive noise model is the simpler case of the two and can be described by a linear model

*J*(*X*) = *I*(*X*) + *n*(*X*) (3)

where *J*(·) is the noised image, *I*(·) is the original image and *n*(·) is the noise. Multiplicative noise is generally expressed by a multiplicative model

*J* (*X*) = *I* (*X*) * *n*(*X*) (4)

It is well known that multiplicative noise appears much worse in bright image regions than dark regions since it multiplies the gray intensities.

Speckle noise is generally treated as multiplicative noise and can be modelled using equation (4). Thus, compared with other types of noise, speckle noise is generally difficult to be removed. Our research below shows that the conventional bilateral filter described in equation (1) and (2) generally gets bad results when it is applied to speckle reduction directly. Thus, the bilateral filter described in (1) and (2) needs improvement or enhancement so that it can be applied to reduce the speckles in images effectively. In order to do this, we will first analyze the behavior of of the bilateral filter in equation (1) in a homogenous region for both additive noise and multiplicative noise, then we will propose an adaptive bilateral filter for speckle reduction.

Let *J*(*Y*) and *J*(*X*) be two different pixels from image *J.* If *J* is corrupted by additive noise, then we can use equation (3) to compute the difference between these two pixels

||*J*(*Y*) - *J*(*X)*|| *= ||I*(*Y*) + *n(Y*) -*I*(*X*) - *n*(*X*)|| (5)

If both *J*(*Y*) and *J*(*X*) are from the same homogenous region, then we have *I*(*Y*) = *I*(*X*), thus

||*J*(*Y*) - *J*(*X)|| = ||n*(*Y*) - *n*(*X*)|| (6)

Equation (6) means that the difference between any two pixels from the same homogenous region is only related to the difference of the noise. If *J* is corrupted by multiplicative noise, then we can use equation (4) to compute the difference between these two pixels. From equation (4), we have

||*J*(*Y*) - *J*(*X)|| =* ||*I*(*Y*) * *n*(*Y*) - *I(X*) * *n*(*X*)|| (7)

Similarly, if both *J*(*Y*) and *J*(*X*) are from the same homogenous region, then we have *I*(*Y*) = *I*(*X*), thus

||*J*(*Y*) - *J* (*X*)|| = ||*I*(*X*)||||*n*(*Y*) - *n*(*X*)|| (8)

From equation (8), we can understand that the difference between two pixels in the same homogenous region(in the image corrupted by multiplicative noise) is not only related to the difference of the noise. It also depends on the intensity of the region. As is seen in equation (8), the difference is big when the intensity of the region is big while the difference is small when the intensity of the image is small.

The above analysis shows the bilateral filter described in (1) and (2) is not suitable for removing speckle noise, which is multiplicative noise. The reason lies in the difference of the corrupted image has different distributions in different homogenous regions. For example, if *σ
_{r}*

(9)

Equation (9) shows that the normalized difference is only related to the noise and doesn’t depend on the intensities of the region. Thus, the proposed adaptive bilateral filter can be expressed as follows

(10)

where

(11)

Bilateral filter is famous because it is non-iterative, however, the non-iterative bilateral filter doesn’t yield good results. In order to improve its effectiveness, we use iterative bilateral filter. The basic idea of iterative bilateral filter is to use the filtered image obtained by equation (10) as the input of equation (1) and implement it many times, the mathematical expression is as follows:

(12)

where

(13)

Where . Experiments show that iterative bilateral filter gives much better results than the non-iterative bilateral filter.

In order to analyze and monitor the reproduction of cattle, the acquisition of some quantitative parameters is very important. These parameters include diameters, areas and perimeters of the follicles. These parameters can be used to monitor the development and maturity of follicles. In order to get these parameters, we need to segment the follicles.

Many image segmentation methods have been proposed, which includes histogram based methods, edge detection based methods, region based methods, active contour model based methods, etc. Active contour model based methods have drawn a lot of attention in the past decade because of their significant advantages. In this paper, we adopt active contour model based method for the segmentation of the follicles. An active contour or a snake [19] is defined as a controlled continuity contour that is attracted to salient image features. However, there are some disadvantages related to the original model. Thus, many improved active models have been proposed based on the original model. The gradient vector flow (GVF) model is one of them [20][21]. GVF model is designed to overcome one of the disadvantages of original model, i.e. the original model is sensitive to the initialization of the snake. In GVF model, GVF fields are computed by another diffusion process, which can be implemented by minimizing the following energy function [20][21]:

(14)

where *g* is a decreasing function of the edge-force magnitude and is defined as follows:

(15)

Here *k* is a non-negative smoothing parameter for the field (*u*, *v*)*.* The functional described by equation (15) smoothes the force field (*u*, *v*) only when the edge strength is low. Solving the energy functional optimization problem in (14), we can obtain the generalized gradient vector flow, which can be used as external forces that attract the snake to the follicle boundary [20][21].

GVF provides external forces for a snake model, we also need internal forces to smooth the contour. In this paper, we use B-spline to represent the contour instead of the real internal forces. B-spline has been used in snake model in several applications and get pretty good results [22][23][24]. Let the control points be denoted by *P*_{0} through *P _{m}*. The knot-value sequence is a non-decreasing sequence of knot values

*Q _{i}* (

where 3 *≤ i ≤ m* and *t _{i} ≤ t ≤ t_{i+}*

(17)

(18)

When *p*=4, we obtained the blending function of cubic splines. The GVF snake with B-spline is called B-spline GVF snake [23][24][25].

For the segmentation of the follicles, we initialize the B-spline GVF snake using a circle inside each follicle. The circle is represented by B-spline and the number of control points is set to 48 in this paper. Then, starting from the initial contour, the GVF is used to drive the contour to the boundary of the follicle. The evolution of the contours is the same as that in the B-spline GVF snake in single scale proposed by [24].

To test the effectiveness of the proposed bilateral filter, we used both conventional bilateral filter and the proposed bilateral filter to process the synthetic image with speckles and compare the results. Fig.1(a) is the original image and Fig. 1 (b) is the corrupted image by speckles with mean 0 and variance 0.075. In order to demonstrate the effectiveness of the proposed filter and evaluate its performance in speckle reduction and edge preservation, we employed three measures in the experiments for comparison. These three measures are: normalized mean square error (NMSE), noise suppression measure *α* and edge preservation parameter *β*[26]. The NMSE is defined as [26]

(19)

where I_{0} and I are the original image and the corrupted image, respectively, N is the pixel number of the image I_{0} (or) I, and are the means of I and I_{0}, respectively. The NMSE generally represents the difference between the original image and the processed image. The noise reduction measure is defined as [26]

(20)

where

(21)

The edge preservation parameter is given by [26]

(22)

where Δ is the Laplacian operator. Higher *α* and *β* represent better performance in noise reduction and edge preservation.

The conventional bilateral filter and the proposed bilateral filter were applied to process the speckled images. In both filters, *σ
_{d}* was fixed to be 3 and

Fig.1(c) and 1(f) are the best results obtained by the conventional bilateral filter with *σ
_{r}* = 0.3 and the proposed bilateral filter with

All of the above experiments show that the proposed bilateral filter can achieve better performance in noise removal and edge preservation than the conventional bilateral filter.

In this subsection, we will compare the proposed bilateral filter with Gaussian filter and the conventional bilateral filter in speckle reduction using real ultrasound images. Fig.Fig.55 shows the original image and the results obtained by the three filters. Although Gaussian filter may reduce the speckles in the images as seen in Fig.5 (b), the edges and details are very blurred. The useful details in the processed image (see Fig. 5(c)) obtained by the conventional bilateral filter are retained, but there are still many speckles. In Fig.5 (d), we know that the proposed filter can reduce speckles effectively while preserve useful edges and details.

To compare and evaluate the three filters quantitatively, we used them to reduce the speckles in real ultrasound image and then calculated the contrast of the homogenous region and edges in the image. A good filter should preserve the edges and reduce speckles in the image, which means the contrast in homogenous region should be low while the contrast in edges should be high. The contrast measure used in this paper is the measure adopted in [27], which is defined as

(16)

where *c*(*x, y*)*,* the local contrast at pixel (*x, y*), is the Laplacian operation

*c*(*x*, *y*) = 4*I*(*x*, *y*) - {*I*(*x* - 1, *y*) + *I*(*x*, *y -* 1) + *I*(*x* + 1, *y*) + *I*(*x*, *y +* 1)} (17)

where *I*(*x, y*) is the pixel intensity value at pixel (*x*, *y*) of an image, *w* is a region or a set of edge points, and *m* is the number of the pixels in the region or edge points.

Fig.Fig.66 illustrates the contrast values in the preselected homogenous regions and the preselected sets of edge points of 12 follicle images. For the homogenous regions, Fig. 6 (a) shows that the contrast values of the regions obtained by Gaussian filter are smaller than those obtained by the conventional bilateral filter or the proposed bilateral filter. Besides, the proposed bilateral filter obtained the smallest contrast values (all are less than 0.04). These results show that the proposed bilateral filter can achieve the best performance in speckle reduction in homogenous regions. For the set of edge points, Fig.6 (b) shows that Gaussian filter has the smallest contrast values, which indicates that most of the edges have been smeared out. Although the conventional bilateral filter has higher contrast values in the set of edge points, the proposed filter has the biggest contrast values, which means it has higher performance in edge preservation.

After the images were processed, we applied B-spline snake [28] to extract the boundaries of the follicles. Fig.Fig.77 shows the experimental results. Fig.7 (a),(b),(c),and (d) show the boundaries of the follicles extracted by B-spline snake from the original images, the images processed by Gaussian filter, the contional bilateral filter and the proposed filter, respectively. Fig.7(a) shows that the final contour is away from the boundary of the follicle due to the speckles. Although there are less speckles in Fig.7(b), the final contour is also away from the real boundary because the edges are blurred by Gaussian filter. The result of Fig.7(c) is very close to the real boundary of the follicle than the contour in Fig.7(a) and Fig.7(b), but it is still affected by speckles. Fig.7 (d) shows that B-spline snake can accurately locate the real boundary of the follicle filtered by the proposed algorithm.

In order to evaluate the segmentation results, we adopted the segmentation metric, Pratt's quality measurement metric (FOM), which is defined as [29]

(18)

where I_{A} is the number of boundary pixels delineated by an automatic segmentation method, I_{I} is the number of boundary pixels delineated by the technicians. *d(i)* is the Euclidean distance between a boundary pixel of ground truth or delineated by the technicians and the nearest boundary pixel extracted by automatic segmentation, and γ is a scaling constant(0.05 in our experiments).

Fig.Fig.88 shows the FOM values of the 12 images processed by different filters. We can see that the Gaussian filter could improve segmentation, the conventional bilateral filter and the proposed filter achieved better FOM values than Gaussian filter. However, the proposed bilateral filter outperformed the other two filters.

Bilateral filter is a powerful technique in image de-noising due to its stability, and simplicity. The basic idea of bilateral filter is to replace a pixel value by a weighted average of its neighbours in both space and range (pixel values). However, the conventional bilateral filter performs poorly on ultrasound images due to the speckles. From the multiplicative noise model, we investigated a normalized scheme based on the conventional bilateral filter so as to remove the speckles effectively while preserving useful details. For bilateral filter, the parameters including *σ
_{d}*

We presented a normalized bilateral filter for speckle reduction in ultrasound images for follicles segmentation. We compared the conventional bilateral filter with the proposed filter using synthetic speckled images and demonstrated its good performance in speckle reduction and edge preservation. Besides, we also tested the proposed filter, the conventional bilateral filter and Gaussian filter using real ultrasound images of cattle follicles. The contrast values of homogenous regions and edge points set demonstrated the proposed algorithm could achieve the best performance. The segmentation experiments also proved that B-spline snake can accurately find the boundary of the follicles from the filtered images by the proposed method. Experimental results validated the effectiveness and the accuracy of the proposed filter in noise reduction and edge preservation for follicle segmentation.

The authors declare that they have no competing interests.

JT developed the algorithm and wrote non-results part of the paper. SG implemented the algorithm and wrote the result part. QS attended to develop the algorithm. YD and DZ helped data analysis. All authors read and approved the final manuscript.

This work has been supported by the United States Department of Agriculture (award No. 2007-38814-18488). The authors appreciate Dr. Evelin Cuadra and Ms. Melissa C. Mason's work in the acquisition of the cattle follicle images. Publication of this supplement was made possible with support from the International Society of Intelligent Biological Medicine (ISIBM).

This article has been published as part of BMC Genomics Volume 11 Supplement 2, 2010: Proceedings of the 2009 International Conference on Bioinformatics & Computational Biology (BioComp 2009). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/11?issue=S2.

- Knopf L, Kastelic JP, Schallenberger E, Ginther OJ. Ovarian follicular dynamics in heifers: test of twowave hypothesis by ultrasonically monitoring individual follicles. Domes Anim Endoc. 1989;6:111–119. doi: 10.1016/0739-7240(89)90040-4. [PubMed] [Cross Ref]
- Lavoir M, Fortune JE. Follicular dynamics in heifers after injection of PGF2 during the first wave of follicular development. Theriogenology. 1990;33:270. doi: 10.1016/0093-691X(90)90694-O. [Cross Ref]
- Lucy MC, Savio JD, Badinga L, De La Sota RL, Thatcher W. Factors that affect ovarian follicular dynamics in cattle. J Anim Sci. 1992;70:3615–3626. [PubMed]
- Ribadu AY, Nakao T. Bovine Reproductive Ultrasonography: A Review. J. Reprod. Dev. 1999. pp. 13–28.
- Pieterse MC, Taverne MA, Kruip TA, Willemse AH. Detection of corpora lutea and follicles in cows: a comparison of transvaginal ultrasonography and rectal palpation. Vet. 1990;56:552–554. [PubMed]
- Driancourt MA, Thatcher WW, Terequi M, Andriew D. Dynamics of ovarian follicular development in cattle during the estrous cycle, early pregnancy and in response to PMSG. Dom Anim Endocr. 1991;8:21–37. [PubMed]
- Rajamahendran R, Ambrose DJ, Burton B. Clinical and research applications of real-time ultrasonography in bovine reproduction: A review. Can Vet J. 1994;35:563–572. [PMC free article] [PubMed]
- Wachowiak S, Smolikova R, Zurada J, Elmaghraby A. Intelligent Engineering Systems Through Artificial Neural Networks: Smart Engineering System Design: Neural Networks, Fuzzy Logic, Evolutionary Programming, Data Mining and Complex Systems. Vol. 10. ASME Press, New York; 2000. A Neural Approach to Speckle Noise Modeling; pp. 837–842.
- Yu Y, Acton ST. Speckle Reducing Anisotropic Diffusion. IEEE Trans. Image Process. 2002;11:1260–1270. doi: 10.1109/TIP.2002.804276. [PubMed] [Cross Ref]
- Sun Q, Hossack J, Tang J, Acton ST. Speckel reduing Anisotropic Diffusion for 3D Ultrasound images. Computerized Medical Imaging and Graphics. 2004;28:461–470. doi: 10.1016/j.compmedimag.2004.08.001. [PubMed] [Cross Ref]
- Zhang F, Yoo YM, Koh LM, Kim Y. Nonlinear Diffusion in Laplacian Pyramid Domain for Ultrasonic Speckle Reduction. IEEE Transaction on Medical Imaging. 2007;26:200–211. doi: 10.1109/TMI.2006.889735. [PubMed] [Cross Ref]
- Yue Y, Croitoru MM, Bidani A, Zwischenberger JB, Clark JW. Nonlinear Multiscale Wavelet Diffusion for Speckle Suppression and Edge Enhancement in Ultrasound Images. IEEE Transactions on Medical Imaging. 2006;25:297–311. doi: 10.1109/TMI.2005.862737. [PubMed] [Cross Ref]
- Tomasi C, Manduchi R. Bilateral filtering for gray and color images. Proc. Int. Conf. Computer Vision. 1998. pp. 839–846.
- Phelippeau H, Talbot H, Akil M, Bara S. Shot noise adaptive bilateral filter. Proceedings of the 9th International Conference on Signal Processing. 2008. pp. 864–867. full_text.
- Barash D. Fundamental relationship between bilateral filtering, adaptivesmoothing, and the nonlinear diffusion equation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 2002;24:844–847. doi: 10.1109/TPAMI.2002.1008390. [Cross Ref]
- Butt I, Rajpoot N. Multilateral Filtering: A Novel Framework for Generic Similarity-based Image Denoising. Proceedings International Conference on Image Processing (ICIP'2009)
- Zhang M. Bilateral Filter in Image Processing, Master Theis. Louisiana State University; 2009.
- Hu Q, He Q, Zhou J. Multi-Scale Edge Detection with Bilateral Filtering in Spiral Architecture. Proceedings of the Pan-Sydney Area Workshop on Visual Information Processing. 2004;6:29–32.
- Kass M, Witkin A, Terzopolous D. Snakes: Active contour models. International Journal of Computer Vision. 1987;1(4):321–331. doi: 10.1007/BF00133570. [Cross Ref]
- Xu C, Prince JL. Snakes, shapes, and gradient vector flow. IEEE Trans. Image Processing. 1998;7(3):359–369. doi: 10.1109/83.661186. [PubMed] [Cross Ref]
- Xu C, Prince JL. Generalized gradient vector flow external force for active contours. Signal Processing. 1998;71(2):131–139. doi: 10.1016/S0165-1684(98)00140-6. [Cross Ref]
- Foley J, Dam A, Feiner S, Hughes J. Computer graphics: principles and practice. Addison Wesley Publishing; 1996.
- Tang J, Millington S, Acton S, Crandall J, Hurwitz S. Surface extraction and thickness measurement of the articular cartilage from MR images using directional gradient vector flow snake. IEEE Tr. On Biomedical Engineering. 2006;52(5):896–907. doi: 10.1109/TBME.2006.872816. [PubMed] [Cross Ref]
- Tang J, Acton S. vessel boundary tracking for intravital microscopy via multi-scale gradient vector flow snakes. IEEE Tr. On Biomedical Engineering. 2004;51(2):316–324. doi: 10.1109/TBME.2003.820374. [PubMed] [Cross Ref]
- Brigger P, Hoeg J, Unser M. B-spline snakes: a flexible tool for parametric contour detection. IEEE Tr. on Image Processing. 2000;9(9):1484–1496. doi: 10.1109/83.862624. [PubMed] [Cross Ref]
- Satter F, Floreby L. Image enhancement based on nolinear multiscale method. IEEE Tr. on Medical Imaging. 1997;6(6):888–895. doi: 10.1109/83.585239. [PubMed] [Cross Ref]
- Tang J, Sun Q. A 3-D anisotropic diffusion filter for speckle reduction in 3-D ultrasound images. Proc. SPIE. 2009. 7239:72390T-72390T-9.
- Klein A, Lee F, Amini A. Quantitative coronary angiography with deformable spline models. IEEE Tr. on Med. Imaging. 1997;16(5):468–482. doi: 10.1109/42.640737. [PubMed] [Cross Ref]
- Pratt W. Digital image processing. Willy, New York; 1978.

Articles from BMC Genomics are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |