PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nonbiophysBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleNonlinear Biomedical PhysicsJournal Front Page
 
Nonlinear Biomed Phys. 2010; 4: 7.
Published online 2010 December 6. doi:  10.1186/1753-4631-4-7
PMCID: PMC3017044

Simple fractal method of assessment of histological images for application in medical diagnostics

Abstract

We propose new method of assessment of histological images for medical diagnostics. 2-D image is preprocessed to form 1-D landscapes or 1-D signature of the image contour and then their complexity is analyzed using Higuchi's fractal dimension method. The method may have broad medical application, from choosing implant materials to differentiation between benign masses and malignant breast tumors.

1. Aims

Fractal and symbolic methods can be very useful for quantitative assessment and classification of images, based on analysis of experimental data such as microscopic images. Our philosophy is that to be applicable a method should preferably be really simple and easily understandable to non-specialists in the field. Presented methods are very simple and they both draw from multiple disciplines and have multidisciplinary applications.

2. Methods

T. Mattfeldt applied nonlinear deterministic methods from chaos theory to pattern analysis of tumor cells. He compared histological texture in 20 cases of mastopathy with 20 cases of mammary cancer. Epithelial texture plays a central role in histopathological diagnosis and grading of malignancy. T. Mattfeldt pre-processed microscopic 2-dimensional images of tumor cells' epithelium into 1-dimensional 'signals' (so called 'landscapes') and then by embedding these signals in a phase space using 'time-delay' method; he found that correlation dimension differs considerably between benign and malignant mammary gland tumors [1]. We have proposed to use a similar simple method for pre-processing of the surface's 2-D image to construct from any 2-D image two 1-D landscapes, but in the second step we use much simpler and more appropriate in this case Higuchi's fractal dimension method for analysis of the obtained landscapes. It evaluates the total 'length' L(k) of the curve defined by every k-th point and then determines the fractal dimension Df from the scaling that L(k) is proportional to kDf (cf. [2] - [3]).

R.M.Rangayyan and T.M.Nguyen used fractal analysis of contours of breast masses in mammograms to differentiate between malignant and benign tumors. They computed fractal dimension of contours of breast masses obtained from mammographic images calculated either directly from the 2-D contour or from a 1-D 'signature' derived from the contour applying either the ruler method or the box counting method [4]. Again, we propose to use Higuchi's method for analysis of 'signatures' - the method is simpler and leads to comparable results.

Higuchi's fractal dimension, Df , is calculated directly from the data series, without embedding the data in a phase space. It is, in fact, fractal dimension of the curve representing the series, and so it is always between 1 and 2, since a simple curve has, of course, dimension equal 1 and a plane has dimension equal 2. The fractional part of Df is a measure of the series complexity. It should not be misled with fractal dimension of an attractor in the system's phase space.

2.1. Analysis of 'landscapes' obtained from an image

A digitized image is a pattern stored as a rectangular data matrix. Grayscale images are matrices where the matrix elements can take on values from gmin = 0 to gmax = (2b-1), where b denotes the number of bits (gmax = 255 for b = 8). The rendering on a video screen is a presentation of the values from black (0) to white (2b - 1). Most color images are overlays of three monochrome images.

Stepping through a gray value image length of N pixels and height of M pixels row by row we calculate the sum of the gray values in each row, Gm , for m = 1,...,M. Normalizing the numbers by using the largest of those values, Gmmax, we produce the series of real numbers NGS [0 1] that we call 'horizontal landscape'

NGSm=Gm/Gmmax(m=1,...,M)
(1)

and we call Higuchi's fractal dimension of this NGS series Dh.

Similarly, stepping through the same image length of N pixels and height of M pixels column by column we calculate the sum of the gray values in each column, Gn , for n = 1,...,N. Normalizing the numbers by using the largest of those values, Gnmax , we produce the series of real numbers NGS[0,1] that we call 'vertical landscape'

NGSn=Gn/Gnmax(n=1,...,N)
(2)

and we call Higuchi's fractal dimension of this NGS series Dv.

Both landscapes are then analyzed using Higuchi's fractal dimension method. ([5,6]). Normalization in (1) and (2) is convenient but not really necessary since fractal dimension is invariant with respect to scaling of the data.

2.2. Analysis of 'signatures' obtained from an image contour

Having the contour of an image specified in any rectangular coordinate system, i.e. by the set of pairs (xi , yi ) such that pairs i-1, i, i+1 correspond to consecutive points on the contour for any i = 1,...,I; the first point i = 1 may be chosen arbitrary and the point i = I+1 coincides with the point i = 1 i.e. the contour is a closed planar curve.

We calculate arithmetic averages, x0 and y0 of coordinates of all contour points and transform rectangular coordinates into polar ones; it is enough to calculate r-coordinate

ri2=(xix0)2+(yiy0)2
(3)

The series ri is a 1-D 'signature' of the 2-D contour and we analyze signatures of breast masses contours using Higuchi's method.

3. Results

In [5] we have demonstrated that differences in fractal dimension of horizontal and vertical landscapes characterize surface texture (Figure (Figure1)1) and that one may quantitatively characterize surface roughness analyzing landscapes obtained from images (Figure (Figure2);2); such quantitative comparison may be used if the images are obtained using the same magnification. As examples we took surface images provided on Internet by T. Randen [7].

Figure 1
Example of assessment of surface texture using fractal analysis of landscapes.
Figure 2
Example of assessment of surface roughness using fractal analysis of landscapes.

We have used the same method for roughness assessment of implant materials based on analysis of SEM images (cf. [6], we analyzed images by C. Giordano et al. [8]). For example, we can compare quality of surfaces of implant materials for orthopedic prostheses - titanium-coated untreated and treated with different chemicals (Figure (Figure3,3, ,44).

Figure 3
Images of titanium surface after BSP biomimetic treatment and corresponding landscapes - horizontal (row by row) and vertical (column by column); magnifications 700× (a-c), 3500× (d-f); fractal dimension calculated in a window of 128 points, ...
Figure 4
Assessment of surface roughness of implant materials for orthopedic prostheses - TI - titanium-coated untreated, ETC - chemically etched titanium, commercially available, BSP - surface after new biomimetic treatment; magnifications 700×.

Surface roughness plays an important role in cell adhesion to the surface, so quality of materials used for implants depends on their surface properties. The greater is fractal dimension the better it is as implant material. Experiment with culturing cells showed that cell adhesion is really the best for BSP-treated. surface significantly increasing cell proliferation [8]. BSP surface shows also evident multifractal properties - fractal dimension strongly decreases with magnification (Figure (Figure3.),3.), while for untreated titanium coated surface fractal dimension does not change with magnification; when a cell culture grows on such a surface fractal dimension decreases in comparison with that of 'naked' surface [6].

Our method may also help to distinguish between different types of cancer (Figure (Figure55.).

Figure 5
Microscopic images of two types of breast cancer, Carcinoma Ductale (left) and Carcinoma Lobulare (right) and the landscapes of these images below (horizontal landscapes in the middle and vertical landscapes in the bottom) with their Higuchi's fractal ...

We have also applied analysis of Higuchi's fractal dimension to the contours of breast masses (cf. [4]). Signatures of contours (cf. Eq. (2)) of benign masses show significantly higher values of Higuchi's fractal dimension than those of malignant breast tumors (Figure (Figure6.6. and Table Table1.).1.). We lack numerical mammographic images of very high quality to test our method of landscapes calculated for cases of breast tumors, benign masses, cysts etc. in comparison with normal breast tissues.

Figure 6
Contours of a benign mass (upper left) and of malignant breast tumor (upper right) their signatures (cf. Eq. (3), middle row) and the signatures' Higuchi's fractal dimension (lower row). Signature of a malignant tumor shows lower fractal dimension than ...
Table 1
Higuchi's fractal dimension of the signatures of benign masses and of malignant breast tumors calculated either from the whole signature at once ('Global'), or calculated in moving window shifted in each step one element of a signature (Figure 6. c. and ...

While the contour of a benign breast mass (Figure (Figure6a.)6a.) seems to be more regular than the contour of a malignant breast tumor (Figure (Figure6b.)6b.) fractal dimension of the malignant breast tumors is lower than fractal dimension of benign breast masses (Table (Table1.).1.). If the contours are blown up one may observe that these of benign masses show many small irregular than those of malignant tumors. That is why signature of a benign mass shows many small 'fluctuations' while that of a malignant tumor does not (cf. Figure Figure6c.6c. and Figure Figure6d.6d. respectively) so leading to differences in their fractal dimension. Problem of calculation of the length of coast-line considered by Mandelbrot is quite analogous; in fact, fractal dimension turned out to be the best characteristics that actually gives possibility to compare properties of different coast-lines [9].

We analyzed 37 cases of benign masses and 20 cases of malignant tumors. Table Table1.1. shows that there exist significant differences in mean values of Df between signatures of benign mass contours and those of malignant tumors contours - these of benign masses are significantly higher. Also there are only small differences in Df if calculated from the whole signature at once ('Global' columns) or if calculated in a sliding window and only then averaged ('Window' columns), so our method enables really quick data analysis of the whole signature at once for check-up examinations. Standard deviations of all mean Df values are small and the ranges (MeanDf-bengn ±_SDbenign) and (MeanDf-malignant ±_SDmalignant) do not overlap. There are some outliers in each group - they would need further more detailed examination.

Conclusions

Fractal dimension of landscapes obtained from surface images does change with the surface properties. The smoother is a surface, i.e. the smaller are its unevenness at any particular scale, the greater is fractal dimension of any landscape obtained from an image of this surface at given magnification. If a surface shows anisotropic roughness properties (texture) then fractal dimensions of its horizontal and vertical landscapes differ from one another. We demonstrated that our method may be applied for choosing better implant materials. The same method may be used in histology to help distinguish between different types of cancer.

Fractal analysis of signatures of contours of breast masses may help in differentiating between mammographic images of benign masses and malignant tumors in screening medical examinations. The method is quick - one may analyze the whole signature at once to calculate Higuchi's fractal dimension of the signature.

Our method draws from multiple disciplines and may find multidisciplinary applications. The same fractal data-processing method may be used for extraction, fusion, and visualization of multi-modal information from (nano)sensors. as well as in hybrid modeling of living organism - the method is computationally effective and may be applied in real-time.

For more detailed fractal analysis of Rangayyan's data cf. [10].

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

The authors contributed equally to this article. All authors read and approved the final manuscript.

Acknowledgements

We thank Prof. R.M.Rangayyan (University of Calgary, Alberta, Canada) for data of contours of breast masses, used for calculations in Figure Figure6.6. and Table Table1.1. We also thank Prof. D.Maslinska (MRC PAS, Warsaw) for images of breast cancer (Figure (Figure5.).5.). This material was partially presented during the Fifth International Symposium 'Fractals in Biology and Medicine', Locarno, Switzerland, March 2009. This work was partially supported by IBBE PAS statutory activity 4.4/st/10.

References

  • Mattfeldt T. In: Fractals in Biology and Medicine. Losa GA, Merlini D, Nonnenmacher TF, Weibel ER, editor. II. Birkhäuser, Basel, Boston, Berlin; 1997. Spatial Pattern Analysis using Chaos Theory: A Nonlinear Deterministic Approach to the Histological Texture of Tumours; pp. 50–72.
  • Klonowski W. Nonlinear Biomedical Physics. 5. Vol. 1. BioMed Central, London; 2007. From conformons to human brains: an informal overview of nonlinear dynamics and its applications in biomedicine.http://www.nonlinearbiomedphys.com/content/1/1/5 [PMC free article] [PubMed]
  • Klonowski W. Signal and Image Analysis Using Chaos Theory and Fractal Geometry. Machine Graphics & Vision. 2004;9:403–431.
  • Rangayyan RM, Nguyen TM. Fractal Analysis of Contours of Breast Masses in Mammograms. Journal of Digital Imaging. 2007;20(3):223–237. doi: 10.1007/s10278-006-0860-9. [PMC free article] [PubMed] [Cross Ref]
  • Klonowski W, Olejarczyk E, Stepien R. A new simple fractal method for nanomaterials science and nanosensors. Materials Science-Poland. 2005;23(3):607–612.
  • Klonowski W, Olejarczyk E, Stepien R. In: Computer Recognition Systems. Kurzynski M, Wozniak M, Puchala E, Zolnierek A, editor. Springer Verlag, Berlin, Heidelberg; 2005. SEM Image Analysis for Roughness Assessment of Implant Materials; pp. 553–560. full_text.
  • Brodatz’s Database cf. http://www.ux.uis.no/~tranden/brodatz.html
  • Giordano C. et al. Titanium for osteointegration: Comparison between a novel biomimetic treatment and commercially exploited surface. J Appl Biomaterials & Biomechanics. 2004;2:35–44. [PubMed]
  • Mandelbrot BB. How long is the coast of Britain? Statistical self-similarity and fractal dimension. Science. 1967;155:636–638. doi: 10.1126/science.156.3775.636. [PubMed] [Cross Ref]
  • Stepien R, Stepien P. Analysis of Contours of Tumor Masses in Mammograms by Higuchi's Fractal Dimension. Biocyb Biomed Eng. 2010. in press .

Articles from Nonlinear Biomedical Physics are provided here courtesy of BioMed Central