Open Access

Fractal frontiers in cardiovascular magnetic resonance: towards clinical implementation

  • Gabriella Captur1, 2,
  • Audrey L. Karperien3,
  • Chunming Li4,
  • Filip Zemrak2, 5,
  • Catalina Tobon-Gomez6,
  • Xuexin Gao7,
  • David A. Bluemke8,
  • Perry M. Elliott1, 2,
  • Steffen E. Petersen2, 5 and
  • James C. Moon1, 2Email author
Journal of Cardiovascular Magnetic Resonance201517:80

Received: 8 April 2015

Accepted: 5 August 2015

Published: 7 September 2015


Many of the structures and parameters that are detected, measured and reported in cardiovascular magnetic resonance (CMR) have at least some properties that are fractal, meaning complex and self-similar at different scales. To date however, there has been little use of fractal geometry in CMR; by comparison, many more applications of fractal analysis have been published in MR imaging of the brain.

This review explains the fundamental principles of fractal geometry, places the fractal dimension into a meaningful context within the realms of Euclidean and topological space, and defines its role in digital image processing. It summarises the basic mathematics, highlights strengths and potential limitations of its application to biomedical imaging, shows key current examples and suggests a simple route for its successful clinical implementation by the CMR community.

By simplifying some of the more abstract concepts of deterministic fractals, this review invites CMR scientists (clinicians, technologists, physicists) to experiment with fractal analysis as a means of developing the next generation of intelligent quantitative cardiac imaging tools.


Cardiovascular magnetic resonance Segmentation Image processing

Fractals-irregularity and complexity in nature

The earliest formal references to fractal geometry were made by Leibniz [1] in the mid-1600s. Centuries later, the first fractal prototype was abstractly introduced (only in passing) by German mathematician Georg Cantor in 1883. But the word ‘fractal’ did not come into existence until at least one century later. Inspired by the Latin fractus, meaning “broken”, the term was first coined in 1975 by Benoit Mandelbrot [1] to describe complex patterns that were self-similar across infinite scales. A fractal object is defined as a rough, fragmented, or detailed geometric shape that can be subdivided into parts, each of which is a reduced copy or approximate copy of the whole, where their self-similarity may be exact, quasi, or statistical.

Theoretical mathematical fractals are indeed infinitely self-similar. We can generate limited practical graphical representations of them by repeating a pattern at different scales in a recursive or iterative loop or by recursion of algebraic equations. Algebraic fractals typically require thousands or millions of iterations before their fractal nature is realised, and thus are usually visualised using computer software. Not surprisingly, widespread appreciation of fractal complexity developed only after the advent of the computer in the 1980s and thanks to Mandelbrot’s work [1].

Natural quasi fractal objects, unlike theoretical fractals but much like graphical representations of fractals, are scale invariant across only a limited range of scales. We are surrounded by natural objects that iterate, branch or spiral, spanning a wide range of scales. Some large-scale examples in the physical world include recursing coastlines, branching tree and river networks, and spiralling galaxies (Fig. 1a) and hurricanes. Some small-scale examples in biology include the spirals of a nautilus and whorls of a seashell (Fig. 1b). Small-scale examples in the human body include the lattices of cancellous bone (Fig. 1c), neuronal dendrites, tumor growth patterns, and specifically for the cardiovascular system, branching vascular networks (Fig. 1d), endocardial trabeculae, and the quasi-fractal ordering of collagen and fibrosis in the diseased myocardium as seen by micro-histology [2].
Fig. 1

Exact self-similar elements cannot always be recognized in naturally occurring fractals. The spiral galaxy (a) is an example of a large-scale fractal in the physical world. Biology is full of fractal objects: from the whorls of molluscs (b), to the woven lattice of human cancellous bone (c); from the branching pulmonary arterial tree (d) to the trabeculated apex of the left ventricle (e). CMR = cardiovascular magnetic resonance

In cardiovascular magnetic resonance (CMR), much of what we see, report, measure and compute in everyday clinical practice also has some quasi-fractal property and is amenable to description and quantification by fractal mathematics, generating an index of their space-filling. To date however, much more emphasis on Fourier analysis and processing of CMR data has existed. Fractal analysis of magnitude images is a more recent application—although more than 100 [36] publications indexed in PubMed have described fractal analysis in magnetic resonance imaging of the brain, only 4 publications exist for CMR [710]. Summing up this biological complexity in medical images is clinically important, to guide treatment decisions and improve disease diagnosis, but attempting to do so using traditional mathematics (perimeter estimates or area under the curve) is unsatisfactory—it will tend to either oversimplify the motif’s detail and/or vary with the iteration being interrogated (Fig. 2). In general, the fractal approach is ideal for measuring complicated image details that are beyond simple calliper measurement, and permits results from different scanners to be meaningfully compared.
Fig. 2

a The first 3 iterations of the Koch coastline, an exact geometrical fractal. It can be quantified by its perimeter, its AUC or its FD. With each successive iteration of the Koch coastline the original pattern is repeated at a finer level, corresponding to how with increasingly greater magnification increasingly fine detail is revealed in fractals. By traditional methods, the AUC will converge on \( \frac{8}{5} \) and the perimeter of the curve after n iterations will be \( \left(\frac{4}{3}\right)n \) times the original perimeter (4 times more lines, \( \frac{1}{3} \) greater length per iteration), and since \( \left|\frac{4}{3}\right|>1, \) perimeter will tend to infinity. These exemplify the inherent problem with traditional mathematics: it is capable of providing only scale-dependent descriptors that give limited insight into the motif’s overarching complexity. The FD of the Koch curve, on the other hand, summarises its complexity independently of scale. At every iteration (from 1 to infinity) the FD is invariant at \( \frac{ \log 4}{ \log 3}\approx 1.26186 \). Biological quasi-fractals are measured by ‘sampling’ them with an imaging ‘camera’ relevant to a particular imaging modality. Different cameras have different resolutions, but in all cases increasing resolution is similar to accumulating iterations on a mathematical fractal. Natural quasi-fractals are self-similar across a finite number of scales only—a lower limit of representation is imposed by the limit of the screen (pixel resolution). For CMR cines, blurring (quite extreme in b) has the same effect as setting a lower resolution for the particular sequence, and this is equivalent to having fewer fractal iterations. With such manipulation, it can be seen that the area of the set changes little (here by 2 %), the perimeter a lot (by 43 %) and the FD less (by 8 %). This implies that high image resolution (and a fractal approach) may not add much value when attempting to measure the left ventricular volume; but image resolution (and a fractal approach) will make a considerable difference when intricate features like trabeculae are the features of interest: the perimeter length or other 1D approach will be less robust than the FD. AUC = area under the curve; d = length of segment; 1D = one-dimension/al; FD = fractal dimension; px = pixels. Other abbreviation as in Fig. 1

By summarising some of the fundamental principles underpinning the science of deterministic fractals, and by pointing to existing tools and approaches, this paper invites CMR scientists to experiment with fractal analysis as a means of developing an alternative breed of quantitative cardiac imaging tools.

How to measure

Geometrically a fractal would exist in between our more familiar topological dimensions (D T ): between the 1st and 2nd D T , or between the 2nd and 3rd, etc. An understanding of the concept of fractal dimensionality begins therefore with at least some understanding of D T and Euclidean dimensionality (D E ) (Fig. 3). Euclidean space refers to an object’s embedding space and encompasses dimensions that we define using Cartesian coordinates (real numbers e.g., x, y and z). Figure 3 explains why some objects will have D T  = D E , while others will have D T  < D E . Unlike the topological and Euclidean dimensions, the fractal dimension (FD) measures the detailed self-similarity of fractals—the space-filling capacity of a set of points embedded in space or its complexity. It is related to D E and D T by Eq 1:
Fig. 3

A line, square or cube all exist in Euclidean space with a certain number of dimensions described classically by D E  = 0 for a single point, 1 for a line (a), 2 for a plane (b) and 3 for a 3D object (d) [38]. The concept of topology is rooted in the idea of connectedness among points in a set. The null (empty) set in topology () has no points and its D T is by definition ‘-1’. A single point or a number of points makes up a ‘countable set’. In topology, a set’s D T is always 1 integer value greater than the particular D T of the simplest form that can be used to ‘cut’ the set into two parts [42]. A single point or a few points (provided they are not connected) are already separated, so it takes ‘nothing’ () to separate them. Thus the D T of a point is 0 (−1 + 1 = 0). A line (a) or an open curve can be severed by the removal of a point so it has D T  = 1. A topological subset such as b can have an interior, boundary and exterior. b has a closed boundary of points (like y). When its interior is empty, b is referred to as a boundary set. Its interior may instead be full of points (like x) that are not boundary points because separating them from the exterior is a neighbourhood of other points also contained in b. All points of the subset that are neither interior nor boundary will form the exterior of b. A line of D T  = 1 is required to split this topological set into 2 parts, therefore the D T of b = 2. Flat disks (c) have D T  = 2 because they can be cut by a line with a D T  = 1. A warped surface can be cut by a curved open line (of D T  = 1) so its D T  = 2 although its D E  = 3. Therefore, while lines and disks have D T  = D E , warped surfaces have D T one less than D E . D E = Euclidean dimension; D T = topological dimension

$$ {D}_T\ \le\ FD\ \le\ {D}_E $$

These definitions also apply to fractal analysis in CMR. The heart itself exists in three-dimensional (3D) space, but diagnostic images provide 2D data a large part of the time, from which we extract patterns. The pattern of a drawn endocardial contour, for example on a left ventricular short axis CMR cine slice, appears more complicated than a simple curved line so its FD will be > 1. Because it partly but not completely ‘fills’ 2D space however its FD will be < 2. Therefore the range of possible FD s for a quasi-fractal object like the endocardial contour extracted from a CMR sequence will be consistently a value between 1 and 2.

The mathematical details of a fractal analysis are generally taken care of by software, but this is typically preceded by some medical image preparation. It may be necessary to generate the needed image format (e.g., grayscale, binary or red-green-blue (RGB) data type) or to remove image complexity unrelated to the feature to be measured. For example, a short-axis cine slice may carry signal originating from the myocardium, blood-myocardial boundary, blood pool, and surrounding tissues, all of which are measurable, either separately or together. To be able to measure the quasi-fractal properties of an endocardial contour (the blood-myocardial boundary) some image transformation would be needed in order to extract its relevant pattern, in particular its binary outline. In a segmented image, derived according to a fixed thresholding rule, the meaning of each single pixel is reduced to the binary logic of existence (pixel present/foreground) and nonexistence (pixel absent/background). Typically, the FD of a binary filled object (e.g., the binary mask of the blood pool) is greater than that of its binary outlined counterpart (e.g., the edge image of the endocardial contour), and the FD of such binary images (whether filled or outlined) will be generally greater than the equivalent FD [11] of the original grayscale object (Fig. 4) [12].
Fig. 4

The 3D FD (between 2 and 3) of the grayscale cine is computed using the differential box-counting algorithm that takes 3D pixel intensity information into account. In the standard box-counting method applied to binary images as either outlines or filled silhouettes, intensity information is lost as foreground pixels are contrasted from the background pixels to derive the 2D FD (range 1 – 2). For the same original image and considering only the mantissa, it is usually the case that the binary FD is greater than the grayscale equivalent. Furthermore, the FD of the filled binary mask would usually be nearer to 2 when compared to the FD of the equivalent binary outline as the FD of the filled areas massively outweigh the FD of the edges. Abbreviations as in Fig. 2

Assuming the preprocessing approaches used (threshold, subtract background, dilate, trace, find edges, binarise or skeletonize, either automatically or manually) are appropriate for the type of image [13, 14], then it is reasonable to expect that the FD of the resultant region of interest (ROI) will closely approximate the real FD of the aspect of the physical object or process being investigated, at least over a range of image resolutions, and that it will encode potentially valuable biological information.

Once the ROI is extracted, the FD can be calculated using many analysis methods (Table 1). Each will compute a different type of FD but fundamentally they all measure the same property of the ROI—they are all meters of complexity. Even for a single method (e.g., box-counting) multiple algorithmic variants may exist (box-counting may use either a conventional, overlapping, folded or symmetric surface scanning approach [15]). The conventional procedure for box-counting (Fig. 5) rests on simple arbitrary scaling and can be applied to structures lacking strictly self-similar patterns. It works by systematically laying a series of grids of boxes of decreasing calibre onto the ROI and counting (at each level) the number of boxes that overlies pixel detail. The FD is derived from the slope of the logarithmic regression line graphing the relationship of box count and scale. The number of data points used to generate these log-log plots is related to the number of measuring steps. Theoretically, given a priori knowledge of the scaling rules, a mathematical fractal would generate data points that lie along a perfect straight line. The point of practical analysis, however, is to find the scaling rule in the first place. For anisotropic biological objects (like left ventricular endocardial contours) as well as for precisely generated fractal images analysed without knowledge of the scaling rule, the data points do not generally lie on a straight line, reflecting sampling limitations as well as limited self-similarity [16], thus the slope is estimated from the regression line for the log-log plot. The choice of image preparation routine and the details of the method used to gather the data for fractal analysis are important as they can either increase or decrease the correlation coefficient of the double logarithmic plot (more linear or more sigmoid fit respectively).
Table 1

List of fractal dimensions that are most commonly used

Dimension Synonym



Author, Year described



Generic term first introduced by Mandelbrot

Mandelbrot, 1983

Hausdorff Hausdorff-Beisicovitch


Uses image coverage by a number of countable spheres; widely used in pure mathematics but less suitable for use with natural fractals

Hausdorff, 1919

Beisicovitch, 1935

Mandelbrot, 1983

Falconer, 1990

Gulick, 1992

Minkowski-Bouligand Kolmogorov


Uses circle sweep like for D H ; easier to evaluate than D H ; outputs usually greater than or equal to D H

Mandelbrot, 1983

Smith, 1989

Schroeder, 1991

Calliper Perimeter-stepping, Divider, Richardson, Compass


Calculates the fractal complexity of a simple continuous perimeter

Richardson, 1961

Mandelbrot, 1967/83

Falconer, 1990

Smith, 1989

Peitgen, 1992

Box-counting Capacity


Uses a grid method to measure the fractal complexity of 2D and 3D noncontiguous outlines commonly encountered in biological structures

Mandelbrot, 1983

Falconer, 1990

Gulick, 1992

Peitgen, 1992



Typically used in the context of clusters and networks; can be applied to surfaces and biological objects

Caserta, 1990

Jelinek, 1998



Used for measuring the dimension of strange attractors in time series analysis.

Gulick, 1992



Uses dense packing by disjoint balls of differing small radii.

Falconer, 1990

Local connected set


Variant of box-counting applied to binary images where they are sampled pixel by pixel according to the local connectedness of each pixel

Landini, 1995



Uses dense packing by disjoint balls of differing small radii.

Falconer, 1990

Grayscale box-counting Differential box-counting Fourier Higuchi’s


Does not require image segmentation; suitable for being performed in an unsupervised manner and most amenable to automation.

Sarkar, 1994

Azemin, 2011

Higuchi, 1988

Fig. 5

Applying fractal analysis to a 2D cine CMR slice (a) at the mid-left ventricular level [9]. Trabecular detail is extracted by a region-based level-set segmentation [40], followed by binarisation (b) and edge-detection (c). Binarisation eliminates pixel detail originating from the blood pool. The edge image is covered by a series of grids (d). The total number of sized d boxes making up this exemplar grid is 16, and the number of boxes N(d) required to completely cover the contour, 14 (2 boxes overlie blank space). For this set, box-counting will involve the application of 86 grid sizes. The minimum size is set to 2 pixels. The maximum size of the grid series is dictated by the dimensions of the bounding box (discontinuous red line) where ‘bounding box’ refers to the smallest rectangle that encloses the foreground pixels. The box diameter for each successive grid is set to drop by d-1 pixels each time. Through the implementation of this 2D box-counting approach, a fractal output of between 1 and 2 is expected. The log-lot plot (e) produces a good fit using linear regression and yields a gradient equivalent to - FD (1.363). d = box dimension; Ln = natural logarithm; N(d) = number of boxes. Other abbreviations as in Figs. 1 and 2

The FD is not the only tool available in fractal geometry—others such as lacunarity also exist that provide a different layer of information relating more to the texture of objects [17]. Lacunarity (λ) mesures the size distribution of gaps (lacunae) in an image, providing a measure of heterogeneity [18]. It is the counterpart to the FD but the two are non-identical (Fig. 6). If an image has few, small, and regular gaps and is translationally and rotationally invariant, it will have low λ; if it has many large and irregular gaps with notable translational and rotational variance, it will have high λ. The translational invariance (spatial heterogeneity [19]) that is measured by lacunarity implies that: 1) λ is highly scale-dependent, meaning an image that appears highly heterogenous at low scale may appear much more homogenous at large scale producing two very different values of λ; and 2) λ (like the related box-counting fractal analysis) may be used to study non-fractal objects. λ and the FD are usually used complementarily, but for some biomedical applications lacunarity may be preferred (e.g., quantification of trabecular bone by MR [20] where the widely varying pattern of emptiness between spicules is the feature of interest, Fig. 1c), and in others the FD is preferred (e.g., endocardial contours with large central emptiness and edge detail, Fig. 5).
Fig. 6

It is possible to construct a family of fractals that share the same FD, but differ sharply in their overall texture so they have uncorrelated values for λ —likewise two objects may have the same λ but very different FD. In a, two 2D binary sets are presented that share the same λ but have different FD. For quantifying myocardial trabecular complexity in CMR cines, FD was chosen over λ for a number of reasons: 1) experiments on grayscale short-axis imaging sequences showed λ was confounded by signal from the central blood pool; 2) as λ measures translational invariance (imagine the binary edge-image rotated clockwise as per curved arrow in b), it is theoretically possible for a heavily but symmetrically trabeculated heart (b, left image) to have a lower value for λ than one with fewer, more irregularly spaced trabeculae (b, right image). On the contrary, if there are more trabeculae, whether regularly or irregularly spaced, FD will always be higher. As the sole objective of this tool was to quantify trabeculae, the extra information on spatial heterogeneity encoded in λ could only have distracted from the biological signal of interest; 3) λ is a very scale-dependent meter and potentially more susceptible to differences in image resolution across vendors and CMR centres compared to FD. λ = lacunarity. Other abbreviations as in Figs. 2 and 3

Previous use of fractal analysis in medicine

Fractal geometry has already found effective research application in the medical imaging field across several modalities (such as plain radiography, retinal photography, ultrasonography [21], computed tomography, MR and nuclear MR [22]). It has been used to study a wide variety of processes: the complex geometries of biological cell types [23]; tumor growth patterns [24]; gene expression [25]; retinopathy [26]; cellular differentiation in space and time [27]; bone and dental matrix composition [24, 28]; brain matter changes [29] etc. Fractal methods are popular and convenient because they lend themselves to automated computer-assisted image processing providing a precise and quantitative metric. Robust measurement of biological complexity in the medical imaging field is clinically important and worth pursuing because fractal indices have been shown to permit early diagnosis of disease (in osteoporosis [20]), predict likelihood of malignancy (in mediastinal nodes imaged by endobronchial ultrasound [21]), predict outcome (of lacunar strokes on the basis of retinal vessel complexity [26]) and measure treatment response (to radiochemotherapy in malignant head and neck tumors [30]).

Utility in CMR

Pertinent to CMR, and for certain applications (e.g., myocardial trabecular quantification), there are clear advantages in using the FD: because it is less susceptible to magnification, it works on different CMR sequences, with different voxel sizes acquired on different platforms; because it is independent of the size of the ROI, it works for small as well as large hearts. There are also potential limitations. For cine imaging, loss of image detail is a particular concern due to partial volume effects at the blood-myocardial boundary in the relatively extended through-plane voxel dimension. Left ventricular cine stacks may be prone to variable spatial resolution but we have previously shown how FD is robust to small changes in slice thickness (6 mm vs. 7 mm vs. 8 mm [10]). Future work should explore whether the higher spatial resolution of computerised tomography provides more suitable image data for fractal analysis than does CMR, especially with respect to vascular trees and probably also myocardial trabeculae provided blood-myocardial contrast is sufficient.

Experimenting with fractal analysis of images in the CMR domain, typically involves the in-house development of scripts written for a specific programming environment (e.g., MATLAB, ImageJ [31], Insight Toolkit [ITK] [32], etc.). It may be possible to repurpose already available tools in the form of commercial and open-source fractal plugins and codes [33]. Examples include: Fractalyse (ThèMA, F), Whinrhizo (Regent Instruments Inc.), Image Pro Plus (Media Cybernetics), FDSURFFT (in MATLAB) and Fraclac for ImageJ [34]. Our group started off with Fraclac and then moved to an in-house MATLAB implementation. We applied fractal analysis to CMR cine data for trabecular quantification. In left ventricular noncompaction (n = 30) compared to healthy volunteers (n = 75) fractal analysis (Fig. 7) revealed FD elevation in the apical half of the left ventricle [8] (1.392 ± 0.010 versus 1.235 ± 0.004). When we studied patients at our centre with hypertrophic cardiomyopathy (n = 107), fractal analysis showed abnormally increased apical FD not only in overt disease, but also in sarcomere gene mutation carriers without left ventricular hypertrophy (G + LVH-, 1.249 ± 0.07) compared to controls (1.199 ± 0.05) [9]. In a multi-centre setting high FD was further shown to predict hypertrophic cardiomyopathy sarcomere gene mutation carriage in G + LVH- (n = 73) [10]. Applied to 2547 participants in the population-representative MESA study, fractal analysis was able to provide ethnically-appropriate normal reference ranges for left ventricular endocardial complexity [35].
Fig. 7

Clinical application of a fractal analysis for trabecular quantification by CMR in LVNC. It is noteworthy how in healthy hearts, it is the mid-LV that holds the greatest fractal complexity (papillary muscles), a fact that is commonly overlooked as the more intricately trabeculated apex commonly distracts. LVNC = left ventricular noncompaction. Authorization for this adaptation has been obtained both from the owner of the copyright in the original work [8] and from the owner of copyright in the translation or adaptation (JCMR)

Whether to measure endocardial complexity or any other imaging feature of interest, all novel CMR fractal tests will invariably need to satisfy the usual STAndards for the Reporting of Diagnostic accuracy studies (STARD) [36]. To become useful clinical tools, they will need to pass the 15 developmental “check-points” [37]. Table 2 underscores how two efforts in this field are still some way off from clinical utility (e.g., the further developed of the two is at step 11—development of normal reference values).
Table 2

The 15 steps needed to turn a fractal tool in a clinically valid test (also considering STARD [39] criteria)

Developmental step

Fractal quantification of trabecular complexity [9]

Fractal quantification of the spatial distribution of pulmonary flow [7]

1. Technical development and theoretical basis of the test

Achieved – method first implemented in Java [8] and then in MATLAB [9] to improve computational efficiency; many segmentation algorithms tested before choosing a region-based level-set function [40]

Achieved – fractal dimension used as an index of pulmonary perfusion heterogeneity; image preparation included a coil inhomogeneity correction

2. Comparison with gold-standard or tissue biopsy (animal models and then human biopsy material)

Achieved – validated against episcopic mouse embryo datasets and using synthetically constructed phantoms with well-known FD: 1) regular geometrical objects (plane, cube surface, sphere surface) and 2) ideal monofractal signals (4th, 5th and 6th iteration of the Sierpinski carpet or 9th, 10th and 11th iteration of the Sierpinski gasket)

Part achieved – validated using 3 MR reference phantoms applied to each patient’s chest

3. Detection of changes in established disease compared with normals

Achieved – FD in left ventricular noncompaction compared to healthy volunteers

Not achieved

4. Correlation with other equivalent cardiac imaging markers

Achieved – correlated with perimeter and with noncompacted/compacted wall thickness ratio [41]

Achieved – comparison is made with relative dispersion and the geometric standard deviation

5. Correlation with other relevant biomarkers

Not achieved

Achieved – data correlated with pulmonary function test from spirometry and repeated for three different inspired oxygen concentrations (normoxia, hypoxia and hyperoxia)

6. Demonstration of the test in more than one condition

Achieved – noncompaction and also subclinical and overt hypertrophic cardiomyopathy (and hypertension, in press)

Not achieved

7. Demonstration of test sensitivity (early disease or change with age)

Achieved – in subclinical hypertrophic cardiomyopathy

Not achieved

8. Demonstration of ability to track changes over time

Not achieved

Not achieved

9. Demonstration of predictive or prognostic value of the test

Achieved – in combination with other CMR imaging markers, high FD shown to predict sarcomere gene mutation carriage in subclinical hypertrophic cardiomyopathy

Not achieved

10. Standardization of the test (reproducibility, different equipment, in non-research settings, quality control, limitations of test)

Achieved – intra- and inter-observer variability, inter-scanner reproducibility, field-strength and slice-thickness independence demonstrated; community roll-out started through open-source release of an OsiriX plugin and development of an equivalent commercial version (in cvi42, Circle Cardiovascular Imaging)

Not achieved

11. Development of robust age/ethnic normal reference ranges

Achieved – through analysis of the Multi-Ethnic Study of Atherosclerosis (in press); robust to multi-centre/multi-vendor implementation

Not achieved

12. Changes in biomarker remain tied to the disease after treatment

Not achieved

Not achieved

13. Demonstration of test as surrogate trial end point

Not achieved

Not achieved

14. Clinical use and regulatory approval of test

Not achieved

Not achieved

15. Prove that test use improves clinical outcome

Not achieved

Not achieved

Not achieved marks a developmental milestone that has not yet been reached/published to our knowledge

Nonetheless, on the broader frontier, there is reason for optimism with regard to developing useful CMR applications in the future. We think that potential, as yet untested applications could include such things as textural analysis to quantify scar in late gadolinium enhancement images; spatiotemporal analysis to track cardiac motion of cine objects; stochastic fractal models [38] to study nonperiodic fluctuations in physiological parameters in MR flow data; and fractal analysis in general to aid pattern recognition in pixel-wise parametric mapping.


Although the description of modern fractal analysis by Mandelbrot occurred more than 40 years ago and in spite of clinical practice bringing us face to face with multifarious fractal features daily, the CMR community is only beginning to evaluate potential applications of fractal analysis to cardiac imaging. This review reminds us of the accessibility of fractal mathematics and methods and aspires to attract more cardiac imagers to the library of efficient fractal analysis tools available, as well as invite them to innovate. A deeper fractal exploration of the human heart by CMR has the ability to teach us new facts relating to cardiac function, haemodynamics and tissue characterisation. With additional validation, software tools based on fractal analysis may ultimately prove to have clinical utility in the field of CMR.



Cardiovascular magnetic resonance




Topological dimension


Euclidean dimension


Fractal dimension


Sarcomere gene mutation carriers without left ventricular hypertrophy




Region of interest





JCM: Higher Education Funding Council for England and the UK National Institute for Health Research, University College London, Biomedical Research Centre; GC: NIHR BRC University College London. DAB: Intramural research program, National Institutes of Health.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

UCL Institute of Cardiovascular Science, University College London
Division of Cardiovascular Imaging, The Heart Hospital, part of University College London NHS Foundation Trust
Centre for Research in Complex Systems, School of Community Health, Charles Sturt University
Department of Radiology, University of Pennsylvania
Cardiovascular Biomedical Research Unit, Barts and the London School of Medicine and Dentistry, Queen Mary University of London
Division of Imaging Sciences & Biomedical Engineering, King’s College London
Circle Cardiovascular Imaging Inc., Panarctic Plaza
Radiology and Imaging Sciences, National Institutes of Health Clinical Center
Barts Heart Centre


  1. Mandelbrot B. The fractal geometry of nature. Ed. Brickmann, J. Free. Co., San Fr. 1982.Google Scholar
  2. Zouein FA, Kurdi M, Booz GW, Fuseler JW. Applying fractal dimension and image analysis to quantify fibrotic collagen deposition and organization in the normal and hypertensive heart. Microsc Microanal. 2014;20:1134–44.View ArticlePubMedGoogle Scholar
  3. Squarcina L, De Luca A, Bellani M, Brambilla P, Turkheimer FE, Bertoldo A. Fractal analysis of MRI data for the characterization of patients with schizophrenia and bipolar disorder. Phys Med Biol. 2015;60:1697–716.View ArticlePubMedGoogle Scholar
  4. Sandu A, Rasmussen Jr IA, Lundervold A, Kreuder F, Neckelmann G, Hugdahl K, et al. Fractal dimension analysis of MR images reveals grey matter structure irregularities in schizophrenia. Comput Med Imaging Graph. 2008;32:150–8.View ArticlePubMedGoogle Scholar
  5. Quodbach J, Moussavi A, Tammer R, Frahm J, Kleinebudde P. Assessment of disintegrant efficacy with fractal dimensions from real-time MRI. Int J Pharm. 2014;475:605–12.View ArticlePubMedGoogle Scholar
  6. Gupta L, Besseling RMH, Overvliet GM, Hofman PA, de Louw A, Vaessen MJ, et al. Spatial heterogeneity analysis of brain activation in fMRI. NeuroImage Clin. 2014;5:266–76.PubMed CentralView ArticlePubMedGoogle Scholar
  7. Arai TJ, Prisk GK, Holverda S, Sá RC, Theilmann RJ, Henderson AC, et al. Magnetic resonance imaging quantification of pulmonary perfusion using calibrated arterial spin labeling. J Vis Exp. 2011;30(51):2712.Google Scholar
  8. Captur G, Muthurangu V, Cook C, Flett AS, Wilson R, Barison A, et al. Quantification of left ventricular trabeculae using fractal analysis. J Cardiovasc Magn Reson. 2013;15:36.PubMed CentralView ArticlePubMedGoogle Scholar
  9. Captur G, Lopes LR, Patel V, Li C, Bassett P, Syrris P, et al. Abnormal cardiac formation in hypertrophic cardiomyopathy: fractal analysis of trabeculae and preclinical gene expression. Circ Cardiovasc Genet. 2014;7:241–8.View ArticlePubMedGoogle Scholar
  10. Captur G, Lopes LR, Mohun TJ, Patel V, Li C, Bassett P, et al. Prediction of sarcomere mutations in subclinical hypertrophic cardiomyopathy. Circ Cardiovasc Imaging. 2014;7:863–7.View ArticlePubMedGoogle Scholar
  11. Sarkar N, Chaudhuri BB. An efficient differential box-counting approach to compute fractal dimension of image. IEEE Trans Syst Man Cybern. 1994;24:115–20.View ArticleGoogle Scholar
  12. Aliahmad B, Kumar DK, Sarossy MG, Jain R. Relationship between diabetes and grayscale fractal dimensions of retinal vasculature in the Indian population. BMC Ophthalmol. 2014;14:152.PubMed CentralView ArticlePubMedGoogle Scholar
  13. Grosgeorge D, Petitjean C, Fares J. Automatic cardiac ventricle segmentation in MR images : a validation study To cite this version. 2010.Google Scholar
  14. Petitjean C, Dacher J-N. A review of segmentation methods in short axis cardiac MR images. Med Image Anal. 2011;15:169–84.View ArticlePubMedGoogle Scholar
  15. Uemura K, Toyama H, Baba S, Kimura Y, Senda M. Uchiyama a. Generation of fractal dimension images and its application to automatic edge detection in brain MRI. Comput Med Imaging Graph. 2000;24:73–85.View ArticlePubMedGoogle Scholar
  16. Avnir D. Applied Mathematics:is the Geometry of Nature Fractal? Science. 1998;279:39–40.View ArticleGoogle Scholar
  17. Cross SS. The application of fractal geometric analysis to microscopic images. Micron. 1994;25:101–13.View ArticlePubMedGoogle Scholar
  18. Manera M, Dezfuli BS, Borreca C, Giari L. The use of fractal dimension and lacunarity in the characterization of mast cell degranulation in rainbow trout (Onchorhynchus mykiss). J Microsc. 2014;256:82–9.View ArticlePubMedGoogle Scholar
  19. Soares F, Janela F, Pereira M, Seabra J, Freire MM. 3D lacunarity in multifractal analysis of breast tumor lesions in dynamic contrast-enhanced magnetic resonance imaging. IEEE Trans Image Process. 2013;22:4422–35.View ArticlePubMedGoogle Scholar
  20. Zaia A, Eleonori R, Maponi P, Rossi R, Murri R. MR imaging and osteoporosis: Fractal lacunarity analysis of trabecular bone. IEEE Trans Inf Technol Biomed. 2006;10:484–9.View ArticlePubMedGoogle Scholar
  21. Fiz JA, Monte-Moreno E, Andreo F, Auteri SJ, Sanz-Santos J, Serra P, et al. Fractal dimension analysis of malignant and benign endobronchial ultrasound nodes. BMC Med Imaging. 2014;14:22.PubMed CentralView ArticlePubMedGoogle Scholar
  22. Bauer WR, Hiller KH, Galuppo P, Neubauer S, Köpke J, Haase A, et al. Fast high-resolution magnetic resonance imaging demonstrates fractality of myocardial perfusion in microscopic dimensions. Circ Res. 2001;88:340–6.View ArticlePubMedGoogle Scholar
  23. Kam Y, Karperien A, Weidow B, Estrada L, Anderson AR, Quaranta V. Nest expansion assay: a cancer systems biology approach to in vitro invasion measurements. BMC Res Notes. 2009;2:130.PubMed CentralView ArticlePubMedGoogle Scholar
  24. Cross SS. Fractals in pathology. J Pathol. 1997;182:1–8.View ArticlePubMedGoogle Scholar
  25. Aldrich PR, Horsley RK, Ahmed YA, Williamson JJ, Turcic SM. Fractal topology of gene promoter networks at phase transitions. Gene Regul Syst Bio. 2010;4:75–82.PubMed CentralView ArticlePubMedGoogle Scholar
  26. Doubal FN, MacGillivray TJ, Patton N, Dhillon B, Dennis MS, Wardlaw JM. Fractal analysis of retinal vessels suggests that a distinct vasculopathy causes lacunar stroke. Neurology. 2010;74:1102–7.PubMed CentralView ArticlePubMedGoogle Scholar
  27. Tiraihi A, Tiraihi M, Tiraihi T. Self-organization of developing embryo using scale-invariant approach. Theor Biol Med Model. 2011;8:17.PubMed CentralView ArticlePubMedGoogle Scholar
  28. Amer ME, Heo M-S, Brooks SL, Benavides E. Anatomical variations of trabecular bone structure in intraoral radiographs using fractal and particles count analyses. Imaging Sci Dent. 2012;42:5–12.PubMed CentralView ArticlePubMedGoogle Scholar
  29. Doan NT, van Rooden S, Versluis MJ, Buijs M, Webb AG, van der Grond J, et al. An automated tool for cortical feature analysis: Application to differences on 7 Tesla T2 (*) -weighted images between young and older healthy subjects. Magn Reson Med. 2014;00:1–9.View ArticleGoogle Scholar
  30. Abramyuk A, Hietschold V, Appold S, von Kummer R, Abolmaali N. Radiochemotherapy-induced changes of tumour vascularity and blood supply estimated by dynamic contrast-enhanced CT and fractal analysis in malignant head and neck tumours. Br J Radiol. 2015;88:20140412.View ArticlePubMedGoogle Scholar
  31. Abramoff MD, Magalhães PJ, Ram SJ. Image Processing with ImageJ. Biophotonics International. 2004;11:36--42.Google Scholar
  32. Yoo TS, Ackerman MJ, Lorensen WE, Schroeder W, Chalana V, Aylward S, et al. Engineering and algorithm design for an image processing Api: a technical report on ITK--the Insight Toolkit. Stud Health Technol Inform. 2002;85:586–92.PubMedGoogle Scholar
  33. Mancardi D, Varetto G, Bucci E, Maniero F, Guiot C. Fractal parameters and vascular networks: facts & artifacts. Theor Biol Med Model. 2008;5:12.PubMed CentralView ArticlePubMedGoogle Scholar
  34. A. K. FracLac for ImageJ. 1999–2013.
  35. Captur G, Zemrak F, Muthurangu V, Petersen SE, Li C, Bassett P, et al. Fractal analysis of myocardial trabeculations in 2547 subjects: the Multi-Ethnic Study of Atherosclerosis. Radiology 2015:Epub ahead of print.Google Scholar
  36. Bossuyt PM. The STARD Statement for Reporting Studies of Diagnostic Accuracy: Explanation and Elaboration. Clin Chem. 2003;49:7–18.View ArticlePubMedGoogle Scholar
  37. Sado DM, Flett AS, Moon JC. Novel imaging techniques for diffuse myocardial fibrosis. Future Cardiol. 2011;7:643–50.View ArticlePubMedGoogle Scholar
  38. Eke A, Herman P, Kocsis L, Kozak LR. Fractal characterization of complexity in temporal physiological signals. Physiol Meas. 2002;23:R1–38.View ArticlePubMedGoogle Scholar
  39. Bossuyt PM, Reitsma JB, Bruns DE, Gatsonis CA, Glasziou PP, Irwig LM, et al. The STARD statement for reporting studies of diagnostic accuracy: explanation and elaboration. Ann Intern Med. 2003;138:W1–12.View ArticlePubMedGoogle Scholar
  40. Li C, Huang R, Ding Z, Gatenby JC, Metaxas DN, Gore JC. A level set method for image segmentation in the presence of intensity inhomogeneities with application to MRI. IEEE Trans Image Process. 2011;20:2007–16.View ArticlePubMedGoogle Scholar
  41. Zemrak F, Ahlman MA, Captur G, Mohiddin SA, Kawel-Boehm N, Prince MR, et al. The relationship of left ventricular trabeculation to ventricular function and structure over a 9.5-year follow-up. J Am Coll Cardiol. 2014;64:1971–80.View ArticlePubMedGoogle Scholar
  42. Somasundaran P. Encyclopedia of surface and colloid science. Volume 4. Boca Raton, FL: CRC Press, Taylor & Francis Group; 2006:6675.Google Scholar


© Captur et al. 2015