Skip to main content

Fractal frontiers in cardiovascular magnetic resonance: towards clinical implementation


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.

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
figure 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
figure 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
figure 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
figure 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
Fig. 5
figure 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
figure 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
figure 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)

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



D T :

Topological dimension

D E :

Euclidean dimension

FD :

Fractal dimension


Sarcomere gene mutation carriers without left ventricular hypertrophy

λ :



Region of interest


  1. Mandelbrot B. The fractal geometry of nature. Ed. Brickmann, J. Free. Co., San Fr. 1982.

  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.

    CAS  Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    CAS  Article  PubMed  Google 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 Central  Article  PubMed  Google 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 Central  Article  PubMed  Google 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.

    CAS  Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    Article  Google 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 Central  Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    CAS  Article  PubMed  Google Scholar 

  16. Avnir D. Applied Mathematics:is the Geometry of Nature Fractal? Science. 1998;279:39–40.

    CAS  Article  Google Scholar 

  17. Cross SS. The application of fractal geometric analysis to microscopic images. Micron. 1994;25:101–13.

    CAS  Article  PubMed  Google 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.

    CAS  Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    Article  PubMed  Google 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 Central  Article  PubMed  Google 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.

    CAS  Article  PubMed  Google 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 Central  Article  PubMed  Google Scholar 

  24. Cross SS. Fractals in pathology. J Pathol. 1997;182:1–8.

    CAS  Article  PubMed  Google 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 Central  CAS  Article  PubMed  Google 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 Central  CAS  Article  PubMed  Google 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 Central  Article  PubMed  Google 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 Central  Article  PubMed  Google 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.

    CAS  Article  Google 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.

    CAS  Article  PubMed  Google Scholar 

  31. Abramoff MD, Magalhães PJ, Ram SJ. Image Processing with ImageJ. Biophotonics International. 2004;11:36--42.

  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.

    PubMed  Google 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 Central  Article  PubMed  Google 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.

  36. Bossuyt PM. The STARD Statement for Reporting Studies of Diagnostic Accuracy: Explanation and Elaboration. Clin Chem. 2003;49:7–18.

    CAS  Article  PubMed  Google Scholar 

  37. Sado DM, Flett AS, Moon JC. Novel imaging techniques for diffuse myocardial fibrosis. Future Cardiol. 2011;7:643–50.

    Article  PubMed  Google Scholar 

  38. Eke A, Herman P, Kocsis L, Kozak LR. Fractal characterization of complexity in temporal physiological signals. Physiol Meas. 2002;23:R1–38.

    CAS  Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    Article  PubMed  Google 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.

    Article  PubMed  Google Scholar 

  42. Somasundaran P. Encyclopedia of surface and colloid science. Volume 4. Boca Raton, FL: CRC Press, Taylor & Francis Group; 2006:6675.

Download references




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.

Author information

Authors and Affiliations


Corresponding author

Correspondence to James C. Moon.

Additional information

Competing interests

Academic collaboration with Circle CVI.

The authors declare that they have no competing interests.

Authors’ contributions

All authors have contributed significantly to the submitted work: JCM and GC wrote the article. ALK, CTG, XG, PB, FZ, SEP, DAB and CL provided expert advice and critical review of the manuscript. All authors read and approved the final manuscript.

Rights and permissions

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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Captur, G., Karperien, A.L., Li, C. et al. Fractal frontiers in cardiovascular magnetic resonance: towards clinical implementation. J Cardiovasc Magn Reson 17, 80 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Cardiovascular magnetic resonance
  • Segmentation
  • Image processing