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

**|**Nonlinear Biomed Phys**|**v.4; 2010**|**PMC3017044

Formats

Article sections

Authors

Related links

Nonlinear Biomed Phys. 2010; 4: 7.

Published online 2010 December 6. doi: 10.1186/1753-4631-4-7

PMCID: PMC3017044

Wlodzimierz Klonowski: lp.waw.bibi@nolkw; Robert Stepien: lp.waw.bibi@pets; Pawel Stepien: lp.waw.bibi@neipetsp

Received 2010 October 27; Accepted 2010 December 6.

Copyright ©2010 Klonowski et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), 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.

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.

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.

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 *D _{f }*from the scaling that

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, *D _{f }*, 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

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 *g _{min }= 0 *to

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, *G*_{m }, for *m = 1,...,M*. Normalizing the numbers by using the largest of those values, **G**_{mmax}, we produce the series of real numbers *NGS *_{ℇ}*[0 1] *that we call *'horizontal landscape'*

$$\begin{array}{cc}\mathit{N}\mathit{G}{\mathit{S}}_{\mathit{m}}={\mathit{G}}_{\mathit{m}}/{\mathit{G}}_{\mathit{m}\mathit{m}\mathit{a}\mathit{x}}& (\mathit{m}=\mathit{1},\mathrm{...},\mathit{M})\end{array}$$

(1)

and we call Higuchi's fractal dimension of this *NGS *series *D*_{h}.

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, *G*_{n }, for *n = 1,...,N*. Normalizing the numbers by using the largest of those values, *G*_{nmax }, we produce the series of real numbers *NGS*_{ℇ}*[0,1] *that we call '*vertical landscape*'

$$\begin{array}{cc}\mathit{N}\mathit{G}{\mathit{S}}_{\mathit{n}}={\mathit{G}}_{\mathit{n}}/{\mathit{G}}_{\mathit{n}\mathit{m}\mathit{a}\mathit{x}}& (\mathit{n}=\mathit{1},\mathrm{...},\mathit{N})\end{array}$$

(2)

and we call Higuchi's fractal dimension of this *NGS *series *D*_{v}.

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.

Having the contour of an image specified in any rectangular coordinate system, i.e. by the set of pairs (*x _{i }, y_{i }*) such that pairs

We calculate arithmetic averages, *x _{0 }*and

$${\mathit{r}}_{\mathit{i}}^{\mathit{2}}={({\mathit{x}}_{\mathit{i}}-{\mathit{x}}_{\mathit{0}})}^{\mathit{2}}+{({\mathit{y}}_{\mathit{i}}-{\mathit{y}}_{\mathit{0}})}^{\mathit{2}}$$

(3)

The series *r _{i }*is a 1-D '

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].

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).

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.).

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.

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 *D*_{f }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 *D _{f }*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

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].

The authors declare that they have no competing interests.

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

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.

- 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**

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. |