Research Article  Open Access
Amir Lashgari, Shahriar Ghamami, Saeedeh Shahbazkhany, Guillermo SalgadoMorán, Daniel GlossmanMitnik, "Fractal Dimension Calculation of a ManganeseChromium Bimetallic Nanocomposite Using Image Processing", Journal of Nanomaterials, vol. 2015, Article ID 384835, 9 pages, 2015. https://doi.org/10.1155/2015/384835
Fractal Dimension Calculation of a ManganeseChromium Bimetallic Nanocomposite Using Image Processing
Abstract
Bimetallic materials, which have the ability to convert heat change into mechanical movement, normally consist of two bonded strips of dissimilar metals that expand at different rates. We describe how we made a manganesechromium (MnCr) bimetallic nanocomposite using the centrifuge method and a lowtohigh approach. We conducted scanning electron microscope (SEM) imaging, energydispersive Xray spectroscopy (EDX) analysis, and Xray diffraction spectra of the nanocomposite to prove its identity. We examined how centrifuge speed, process time, and the use of an “intruder agent” affected the properties of the material. The fractal dimension is a significant factor that can be used to approximate the surface roughness, the texture segmentation, and an image of the studied compounds. We calculated the technique of fractal dimensions using imageprocessing values on a computer and histogram plot with the SEM image of the MnCr bimetallic nanocomposite using MATLAB software. We applied the Statistical Package for the Social Sciences software for statistics data extracted from the SEM image of the nanocomposite and obtained the following results: mean = 1.778, median = 1.770, max = 1.98, min = 1.60, skewness = 0.177, range = 0.38, and harmonic mean = 1.771 for fractal dimension of the SEM image.
1. Introduction
Metal matrix composites (MMCs) are an important engineering material due to their various applications. Lightweight MMCs, due to their significant hardness and strength, have found wide applications in the aerospace and automotive industries. However, these applications have remained limited due to the lack of stretch of the composites. Composite conversion nanocomposites increase the yield strength and remove these limitations [1, 2]. MMCs have steadily become more favored in recent years for their high strength and stress resistance properties [3–11].
Chromium (Cr) and chromium compounds are widely used in many industrial procedures such as textile dyeing, printing inks, electroplating, leather tanning, chromium plating, and pigments as serious manufacturing resources [12, 13]. Manganese (Mn) is a mineral created in large amounts in both plant and animal tissues. Manganese lowers the temperature at which austenite converts into ferrite, therefore avoiding cementite precipitation at ferrite grain boundaries, and purifies the resulting pearlitic structures. Alloying elements, such as manganese, also enable solutionhardening of ferrite, but this effect is limited compared with that of carbon, nitrogen, phosphorus, and even silicon [14, 15].
Image processing of microstructures involves several preparation steps, division, and postprocessing that finally lead to the extraction of quantitative data. Recently, computerbased imageprocessing methods have advanced quickly, making it possible to quantitatively describe complex colors, shapes, texture properties, and sizes. Imageprocessing methods are regularly used in tandem with mechanical and instrumental devices to change human manipulative work in the show of an assumed process [16, 17].
The term “fractal dimension” has been used to describe several aspects of particle morphology. The fractal dimension of an individual particle, which can be a main particle or a composite particle such as a floc or aggregate, denotes the indiscretion and selfsimilarity of the particle. The fractal dimension parameter is determined by computing the length of the predictable perimeter by taking steps of a known length along the border and including the number of ladders. The border length increases with decreasing step size because a smaller step size allows for the presence of more detail in the particle outline [18–20]. The word “fractal” is derived from the Latin “fractious” meaning fragmented; a fractal is defined as a mathematical object with a fractional (noninteger) dimension. It seems that wherever there is structural disorder and chaos, fractal geometry can be an efficient way of describing and analyzing it.
Fractal dimension is a very good mathematical method for the study of unequal geometric substances. Fractal dimension can reflect characteristics appearing in nature. The application of fractal system allows us to explain numerous states of fragmenting and branching occurring in ecological, biological, and other systems [21–24]. The output is a binarized image with a onepixelwide boundary appropriate to be used for boxcounting calculations to approximate the fractal dimension of the boundary. We developed a MATLAB (MathWorks, Natick, Massachusetts, USA) program based on theory to implement the boxcounting technique for fractal dimension calculations.
Several specific theoretical and numerical models have been proposed in order to assess the fractal dimension. Recent studies have examined dimensions, but the majority of them were only one or twodimensional. These dimensions are a clear limitation to the study of complex and more general strengthening geometries [25, 26].
Manganese and chromium have been used in this nanocomposite; potassium permanganate and potassium dichromate were used as the starting materials. Scanning electron microscope (SEM) images of the nanocomposites may be regular or irregular and their fractal dimension can be calculated. It is interesting to compare the fractal images of the nanocomposites. For example, in this paper, we investigated the size of the nanoparticles and calculated their fractal dimension. We also measured the average, minimum, and maximum of fractal dimensions and image’s histogram and normalized image’s histogram plot of the SEM image [27–29].
2. Experimental
2.1. Materials and Instruments
The starting materials were obtained from Merck (Berlin, Germany) and were used without further purification. The surface morphology of the product was characterized by using a LEO1430VP scanning electron microscope (Zeiss, Germany) with an accelerating voltage of 15 kV. Xray powder diffraction (XRD) measurements were performed using a Philips X’pert diffractometer (PANalytical, Almelo, Netherlands) with monochromatized CuKα radiation. The sizes of the selected samples were estimated using the Scherer method. For identification, the SEM samples were coated with gold.
2.2. Synthesis of Nanocomposite Manganese and Chromium
Potassium permanganate (KMnO_{4}) and potassium dichromate (K_{2}Cr_{2}O_{7}) were used as starting materials to prepare this nanocomposite. Initially, 0.2 g of potassium permanganate was transferred to a test tube and 0.01 g of clay was added as an intruder agent. Ethanoic acid, CH_{3}COOH (4 mL), was then added to the mixture and the tube was quickly transferred to a centrifuge and run at 3000 rpm for 45 minutes. The sample was then taken out of the centrifuge and put aside. Then, 0.5 g of potassium dichromate and 0.01 g of clay were added to another test tube. NH_{3} (5 mL) solution was added and the tube was placed in a centrifuge and run at 3000 rpm for 45 minutes.
2.3. Characterization of the Nanocomposite
The XRD technique was used to determine the composition of the milled powder. The morphology and weight and atomic percentage of the nanocomposite were observed using a SEM image and EDX analysis, respectively. The contents of the first tube were thoroughly stirred and transferred to a second tube. This second tube was returned to the centrifuge and run at 3000 rpm for 45 minutes. The test tube was then removed from the centrifuge and the precipitate was separated out using filter paper. After drying, the precipitate was imaged with SEM instrument, and the atomic percentages were achieved by EDX analysis. So, the XRD spectrum was obtained to show that a manganesechromium nanocomposite had been produced.
2.4. Calculation of Fractal Dimension and Image Analysis
Dimension is one of the most significant ideas in fractal geometry. Now, there are numerous descriptions of this quantity, including packing, dimension, boxcounting dimension, Hausdorff dimension, and modified boxcounting dimension. Since the mathematic design and approximation of the boxcounting dimension are straightforward, this quantity has become one of the most commonly used general dimensions [30, 31]. The calculations and analysis were done by programming with MATLAB and Statistical Package for the Social Sciences (SPSS; IBM, Armonk, New York, USA) software [32]. Moreover, by declaring that the calculation results were comparable, all of the images were processed using a similar crop size and threshold value. Using this method, we calculated the fractals dimension theoretically by applying the following formula:where is the fractal dimension and is the number of squares including a part of the considered fractal dimensions.
Considering the resulting numbers, we drew a logarithmic diagram in which the vertical and horizontal axes were and , respectively, and the coefficient of the box size was 1.3. We used a program written in MATLAB to calculate the fractal dimension. There are several methods for calculating the fractal definition including the selfsimilarity definition [33], the mass fractal dimension [34], the sparesphere fractal dimension [35], the Euclid dimension [36], the grid style [37], the analytical approach for fractal diminishing [38], and the boxcounting fractal dimension [39].
We used the selfsimilar method for determining the fractal dimension from the SEM image, which employed the boxcounting method for the fractal dimension.
For the boxcounting fractal dimension we suppose that is a nonempty and bounded subset of and we suppose that is the fewest number of collections with a maximum dimension of that can cover . The boxcounting dimensions under and above : if they are equal we say they are equal to the boxcounting dimension or the box dimension; then we defineThe circle of cubes in coordinates from implies cubes in the following form: If , then are integers. Finally, . If , then . If we limit , then . Finally, .
3. Results and Discussion
3.1. Analysis of the Effective Factors for Creating the Nanomaterial
It is important to note that we used centrifuges to make the nanomaterials with the lowtohigh method. After we mixed the solutions they were immediately put into the centrifuge, and an intruder agent was used to ensure that the particles were not connected. This intruder agent must be inert and must not participate in the main reaction. Centrifugal force ensures that the particles produced are not connected and remain on the nanoscale. The advantages of this method are its low cost, rapidity, and ease of application.
Here, we examine the preparation of a manganesechromium bimetallic nanocomposite. The XRD spectrum and SEM image and EDX analysis prepared for the sample both confirmed the creation of a nanocomposite. In the XRD spectrum (Figure 1), no distinguishing peaks are visible and these data did not identify the substance. This outcome is normal because composites are amorphous and their spectra are wide and uniform, unlike the spectra of crystalline materials that have clear and sharp peaks.
(a)
(b)
The SEM image (Figure 2) shows that the particles are smaller than 100 nm, which also confirms that this material is a nanomaterial. We were therefore able to confirm that we made a manganesechromium bimetallic nanocomposite. We then obtained 30 random pictures from the SEM image and calculated their fractal dimensions using MATLAB software [40]. For these 30 fractal dimensions we obtained image’s histogram and normalized image’s histogram, mean, median, max and min, range, skewness, and harmonic mean using SPSS software [41].
The EDX analysis (IDFix report) for the manganesechromium bimetallic nanocomposite is shown in Table 1 and Figure 3. As shown in the IDFix report, the gold present (11.89%) was used for SEM imaging and was not part of the nanocomposite. The primary elements of the nanocomposite were manganese (24.33%), chromium (3.84%), carbon (6.77%), nitrogen (2.61%), potassium (8.73%), and oxygen (36.34%); other elements (silicon (2.24%), sulfur (0.25%), copper (2.07%), and calcium (0.92%)) were related to the clay intruder compound (Figure 3). These results confirmed the correct synthesis of nanocomposite.
 
Element; ^{b}intensity; ^{c}weight percentage; ^{d}atomic percentage. 
3.2. Analysis Using MATLAB and SPSS Data
Thirty pictures were selected at random from the SEM image. We first calculated the fractal dimension for these 30 images using MATLAB software (Table 2). Mean, median, max and min, range, skewness, and harmonic mean were also obtained for these data (Table 2) using SPSS software (Table 3). In Table 3, the mean of the fractal dimension is 1.778, the median is 1.770, the minimum is 1.60, the maximum is 1.98, the range is 0.38, the skewness is 0.177, and the harmonic mean is 1.771.


In statistics and probability theory, the skewness reflects the degree of asymmetry of the probability distribution. Skewness is a measure of the presence or absence of distribution function symmetry. For a perfectly symmetrical distribution, the skewness is zero; for an asymmetric distribution with stretching toward higher quantities, the skewness is positive. Conversely, for an asymmetric distribution with stretching toward smaller quantities, the skewness is negative. Given that the amount of skewness for particles in this nanocomposite was positive (+0.177), this result reveals an asymmetrical distribution of fractal data with stretching toward higher quantities (Table 3). The harmonic mean is a measure of central tendency and is typically calculated when the average of the data is important. The harmonic mean was found to be 1.771 for this nanocomposite (Table 3).
A grayscale digital image consists of separate points of gray tones, or brightness, before incessantly variable tones. A normal image is separated into a number of specific points of brightness, and each of those points is defined via a digital data value. A pixel is the greatest fundamental element of some digital image. Each brightness point is a pixel of the digital image. Essentially, the image histogram shows the distribution of the pixel intensities in the image and is continuously used as a reference. Interactive thresholding can be very effective and afford fast, precise information [42, 43]. To establish the frequency distribution model for fractal data we plotted image’s histogram and normalized image’s histogram (Figure 4) by MATLAB software for the fractal data distribution.
The descriptive statistics of the fractal dimension was listed in Table 4. This table shows fractal dimensions of 30 random separated pieces of the SEM image with frequency and valid and cumulative percent of each piece. When dealing with statistical information, we come through numbers. Sometimes, the numbers are repeated. After a value repeated more than once, it is recognized as the frequency of a number. The total frequency of all modules less than the higher class boundary of a quantified class is named the cumulative frequency of that class. Cumulative frequency is found by addition of the frequency of a class interval and the frequencies of the former intervals up to that class interval. This is another way of expressing frequency distribution. Cumulative frequency is used to determine the number of explanations that lie below or above a specific value. The best advantage of using cumulative percentage is that it offers an easier method to compare different sets of data.

Figure 5 shows a diagram in which there are 30 particles on the horizontal axis and the vertical axis shows the fractal dimension of these particles, which ranges from 1.60 to 1.98. The differences in the fractal dimensions of the particles can be clearly seen, which reveals the nonhomogeneity of the produced particles.
Normal probability plots were generated to check for errors in the normal distribution. If the points of these plots were in a straight line, then normally distributed errors are assumed. According to Figure 6, therefore, which shows a nearly straight line, we conclude that the fractal dimension distribution of the nanocomposite particles is a normally distributed error. Plotting residuals is a useful way of elucidating the goodness of fit of the data. Plotting the residuals versus the corresponding fitted values is a useful model for understanding the common types of disproportion. If the fitted model is appropriate the curve should be symmetrical and the points should be uniformly distributed. Based on Figure 7, one can see that the points are not uniformly distributed, revealing an asymmetric curve that conforms to a skewness of 0.177.
4. Conclusion
Our calculations of the fractal dimensions of 30 randomly selected particles from the SEM image of the nanomaterial and the max, min, and range of the data can be used for analyzing the nanoparticles. We considered homogeneity and uniformity of the SEM image by selecting the 30 random particles from it. In fact, the fractal dimension can be obtained from the size and morphology of the particles. If the max and min of the data are close to each other or the range of data is small then we can conclude that the size and dispersion of the particles in the SEM image are similar and that these created nanoparticles show good homogeneity and uniformity. In the SEM image of this nanocomposite the range was large (about 0.38). For this reason, we concluded that the sizes and morphologies of the particles were different. This conclusion was consistent with the SEM image of the nanocomposite. Initially, it seemed that the picture was uniform. However, after doing the calculations it became clear that the produced particles were not homogenous; despite the uniform dispersion of the particles, their sizes and morphologies were different.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors gratefully acknowledge the financial support from the Research Council of Imam Khomeini International University and Mollasadra student home research. The authors also warmly acknowledge the support of Shahrbano Maleki, Islamic Azad University, Ardabil Branch, Faculty of Science, chemistry group, for his M.S. degree thesis on inorganic chemistry. Daniel GlossmanMitnik is a Researcher at CONACYT and CIMAV and acknowledges partial support from both institutions.
References
 M. Surender, B. Basu, and R. Balasubramaniam, “Wear characterization of electrodeposited NiWC composite coatings,” Tribology International, vol. 37, no. 9, pp. 743–749, 2004. View at: Publisher Site  Google Scholar
 N. Sombatsompop, K. Sukeemith, T. Markpin, and N. Tareelap, “A new experimental apparatus of electrocodeposited system for NiWC composite coatings,” Materials Science and Engineering A, vol. 381, no. 12, pp. 175–188, 2004. View at: Publisher Site  Google Scholar
 G. L. Povirk, A. Needleman, and S. R. Nutt, “An analysis of residual stress formation in whiskerreinforced AlSiC composites,” Materials Science and Engineering A, vol. 125, no. 2, pp. 129–140, 1990. View at: Publisher Site  Google Scholar
 G. L. Povirk, A. Needleman, and S. R. Nutt, “An analysis of the effect of residual stresses on deformation and damage mechanisms in AlSiC composites,” Materials Science and Engineering A, vol. 132, pp. 31–38, 1991. View at: Publisher Site  Google Scholar
 Y.L. Shen, A. Needleman, and S. Suresh, “Coefficients of thermal expansion of metalmatrix composites for electronic packaging,” Metallurgical and Materials Transactions A, vol. 25, no. 4, pp. 839–850, 1994. View at: Publisher Site  Google Scholar
 L. L. Shaw and D. B. Miracle, “Effects of an interfacial region on the transverse behavior of metalmatrix composites—a finite element analysis,” Acta Materialia, vol. 44, no. 5, pp. 2043–2055, 1996. View at: Publisher Site  Google Scholar
 P. J. Ward, H. V. Atkinson, P. R. G. Anderson et al., “Semisolid processing of novel MMCs based on hypereutectic aluminiumsilicon alloys,” Acta Materialia, vol. 44, no. 5, pp. 1717–1727, 1996. View at: Publisher Site  Google Scholar
 J. C. Michel and P. Suquet, “An analytical and numerical study of the overall behaviour of metalmatrix composites,” Modelling and Simulation in Materials Science and Engineering, vol. 2, pp. 637–658, 1994. View at: Publisher Site  Google Scholar
 N. Ramakrishnan, “An analytical study on strengthening of particulate reinforced metal matrix composites,” Acta Materialia, vol. 44, no. 1, pp. 69–77, 1996. View at: Publisher Site  Google Scholar
 D. B. Zahl, S. Schmauder, and R. M. McMeeking, “Mechanical behavior of residually stressed composites with ductile and brittle constituents,” Modelling and Simulation in Materials Science and Engineering, vol. 2, pp. 267–276, 1994. View at: Publisher Site  Google Scholar
 M. Suéry, C. Teodosiu, and L. F. Menezes, “Thermal residual stresses in particlereinforced viscoplastic metal matrix composites,” Materials Science and Engineering A, vol. 167, no. 12, pp. 97–105, 1993. View at: Publisher Site  Google Scholar
 L. Khezami and R. Capart, “Removal of chromium(VI) from aqueous solution by activated carbons: kinetic and equilibrium studies,” Journal of Hazardous Materials, vol. 123, no. 13, pp. 223–231, 2005. View at: Publisher Site  Google Scholar
 X. Guo, G. T. Fei, H. Su, and L. De Zhang, “Highperformance and reproducible polyaniline nanowire/tubes for removal of Cr(VI) in aqueous solution,” The Journal of Physical Chemistry C, vol. 115, no. 5, pp. 1608–1613, 2011. View at: Publisher Site  Google Scholar
 S. R. Yeduru, A. Backen, S. Fähler, L. Schultz, and M. Kohl, “Transformation behaviour of freestanding epitaxial NiMnGa films,” Journal of Alloys and Compounds, vol. 577, no. 1, pp. S353–S357, 2013. View at: Publisher Site  Google Scholar
 K. Hayashi, K. Miyata, F. Katsuki, T. Ishimoto, and T. Nakano, “Individual mechanical properties of ferrite and martensite in Fe0.16 mass% C1.0 mass% Si1.5 mass% Mn steel,” Journal of Alloys and Compounds, vol. 577, no. 1, pp. S593–S596, 2013. View at: Publisher Site  Google Scholar
 Q. Li, M. Wang, and W. Gu, “Computer vision based system for apple surface defect detection,” Computers and Electronics in Agriculture, vol. 36, no. 23, pp. 215–223, 2002. View at: Publisher Site  Google Scholar
 A. Hosseini, S. M. Jafari, H. Mirzaei, A. Asghari, and S. Akhavan, “Application of image processing to assess emulsion stability and emulsification properties of Arabic gum,” Carbohydrate Polymers, vol. 126, pp. 1–8, 2015. View at: Publisher Site  Google Scholar
 D. H. De Boer and M. Stone, “Fractal dimensions of suspended solids in streams: comparison of sampling and analysis techniques,” Hydrological Processes, vol. 13, no. 2, pp. 239–254, 1999. View at: Publisher Site  Google Scholar
 D. H. de Boer, “An evaluation of fractal dimensions to quantify changes in the morphology of fluvial suspended sediment particles during baseflow conditions,” Hydrological Processes, vol. 11, no. 4, pp. 415–426, 1997. View at: Publisher Site  Google Scholar
 B. H. Kaye, “Image analysis techniques for characterizing fractal structures,” in The Fractal Approach to Heterogeneous Chemistry: Surfaces, Colloids, Polymers, D. Avnir, Ed., Wiley, Chichester, UK, 1989. View at: Google Scholar
 P. A. Burrough, “Fractal dimensions of landscapes and other environmental data,” Nature, vol. 294, no. 5838, pp. 240–242, 1981. View at: Publisher Site  Google Scholar
 B. T. Milne, “Measuring the fractal geometry of landscapes,” Applied Mathematics and Computation, vol. 27, no. 1, pp. 67–79, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B. J. West and A. L. Goldberger, “Physiology in fractal dimensions,” The Journal of American Science, vol. 75, pp. 354–365, 1987. View at: Google Scholar
 C. J. Du and D. W. Sun, “Recent developments in the applications of image processing techniques for food quality evaluation,” Trends in Food Science & Technology, vol. 15, no. 5, pp. 230–249, 2004. View at: Publisher Site  Google Scholar
 R. J. Arsenault and M. Taya, “Thermal residual stress in metal matrix composite,” Acta Metallurgica, vol. 35, no. 3, pp. 651–659, 1987. View at: Publisher Site  Google Scholar
 P. J. Withers, W. M. Stobbs, and O. B. Pedersen, “The application of the eshelby method of internal stress determination to short fibre metal matrix composites,” Acta Metallurgica, vol. 37, no. 11, pp. 3061–3084, 1989. View at: Publisher Site  Google Scholar
 B. B. Mandelbrot, The Fractal Geometry of Nature, W. H. Freeman and Company, New York, NY, USA, 1982. View at: MathSciNet
 J.F. Gouyet and B. B. Mandelbrot, Physics and Fractal Structures, Springer, New York, NY, USA, 1996. View at: MathSciNet
 K. Falconer, Fractal Geometry: Mathematical Foundations and Applications, John Wiley & Sons, New York, NY, USA, 2003. View at: Publisher Site
 X. C. Jin, S. H. Ong, and Jayasooriah, “A practical method for estimating fractal dimension,” Pattern Recognition Letters, vol. 16, no. 5, pp. 457–464, 1995. View at: Publisher Site  Google Scholar
 A. K. Bisoi and J. Mishra, “On calculation of fractal dimension of images,” Pattern Recognition Letters, vol. 22, no. 67, pp. 631–637, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Zhang, R. Zhang, B. Hu, and S. Bai, “Fractal image processing and analysis by programming in MATLAB,” in Proceedings of the 8th WSEAS International Conference on Mathematics and Computers in Biology and Chemistry, Vancouver, Canada, June 2007. View at: Google Scholar
 B. Mandelbrot, “How long is the coast of Britain? Statistical selfsimilarity and fractional dimension,” Science Magazine, vol. 156, no. 3775, pp. 636–638, 1967. View at: Publisher Site  Google Scholar
 E. Perfect, “Estimating soil mass fractal dimensions from water retention curves,” Geoderma, vol. 88, no. 34, pp. 221–231, 1999. View at: Publisher Site  Google Scholar
 M. Borkovec, W. de Paris, and R. Peikert, “The fractal dimension of the Apollonian sphere packing,” Fractals, vol. 2, no. 4, pp. 521–526, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Khanbabaei, A. Karam, and G. Rostamizad, “Studying relationships between the fractal dimension of the drainage basins and some of their geomorphological characteristics,” International Journal of Geosciences, vol. 4, pp. 636–642, 2013. View at: Publisher Site  Google Scholar
 P. ParsonsWingerter, B. Lwai, M. C. Yang et al., “A novel assay of angiogenesis in the quail chorioallantoic membrane: stimulation by bFGF and inhibition by angiostatin according to fractal dimension and grid intersection,” Microvascular Research, vol. 55, no. 3, pp. 201–214, 1998. View at: Publisher Site  Google Scholar
 E. Guérin and E. Tosan, “Fractal inverse problem: approximation formulation and differential methods,” Fractals in Engineering: New Trends in Theory and Applications, pp. 271–285, 2005. View at: Publisher Site  Google Scholar
 M. K. Biswas, T. Ghose, S. Guha, and P. K. Biswas, “Fractal dimension estimation for texture images: a parallel approach,” Pattern Recognition Letters, vol. 19, no. 34, pp. 309–313, 1998. View at: Publisher Site  Google Scholar
 The MathWorks, MATLAB—Statistics Toolbox, Release 2012, The MathWorks, Natick, Mass, USA, 2012.
 SPSS, SPSS Statistics for Windows, Version 17.0, SPSS, 2008.
 G. F. Voort, ASM Handbook, Volume 9: Metallography and Microstructures, ASM International, Materials Park, Ohio, USA, 2004.
 L. Wojnar, Image Analysis, Applications in Materials Engineering, CRC Press, Boca Raton, Fla, USA, 1999.
Copyright
Copyright © 2015 Amir Lashgari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.