Correspondence to:
M. Zamir, Department of Applied Mathematics, University of Western Ontario, London, Ontario, N6A 5B7, Canada. Fax
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Parametric Lindenmayer systems (L-systems) are formulated to generate branching tree structures that can incorporate the physiological laws of arterial branching. By construction, the generated trees are de facto fractal structures, and with appropriate choice of parameters, they can be made to exhibit some of the branching patterns of arterial trees, particularly those with a preponderant value of the asymmetry ratio. The question of whether arterial trees in general have these fractal characteristics is examined by comparison of pattern with vasculature from the cardiovascular system. The results suggest that parametric L-systems can be used to produce fractal tree structures but not with the variability in branching parameters observed in arterial trees. These parameters include the asymmetry ratio, the area ratio, branch diameters, and branching angles. The key issue is that the source of variability in these parameters is not known and, hence, it cannot be accurately reproduced in a model. L-systems with a random choice of parameters can be made to mimic some of the observed variability, but the legitimacy of that choice is not clear.
Key Words: arterial tree, L-system, branching, self-similarity, vasculature
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The branching structure of vascular systems has been the subject of much discussion and debate since it was first suggested that these systems may have fractal architecture (
There is no doubt that even here the entire spectrum of vascular systems cannot be treated in the same way. The vascular system of the kidney (
|
|
These considerations suggest that discussion of fractal properties of vascular systems, if it is to relate to function, should be limited to particular categories of systems rather than collectively to all. In this paper, the focus is on arterial systems in the systemic circulation whose main function is the distribution and delivery of blood, rather than filtration or other processing functions. These systems are found to have a predominantly open tree structure and include large and small vessels up to, but not including, the capillary bed. Capillary beds are not included because they are of a different category and their branching structure, characterized by interconnections and closed loops, requires separate treatment (
The term "arterial tree" shall be used generically to mean the branching structure of a main artery in the systemic circulation from its source to its final branches as they reach the capillary bed, but not including the latter. It may represent the aorta from its origin at the left ventricle and include its entire distribution or a smaller system such as that of a coronary artery. It has been found that these and other arterial trees in the systemic circulation have an open tree structure based on repeated bifurcations. This structure is a fractal structure in the sense that each level of the tree is obtained from that preceding it by dichotomous division. However, in the context of physiology, this fractal character is not useful unless the dichotomous divisions are allowed to have the range of properties that arterial bifurcations in the cardiovascular system are known to have (
![]() |
MATERIALS AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The fractal properties of an arterial tree may be discussed in terms of its "fractal dimension" (
Since its inception, L-system formalism has been used and validated as a tool. In plants, where it was first used, its power and utility was measured by its ability to generate the branching patterns of trees and bushes observed in nature. The test is a qualitative one, not quantitative, and the same approach is being used in this paper. We look at whether known patterns in arterial branching can be mimicked by patterns produced within the confines of L-system formalism.
Basic L-System Branching Model
L-system formalism is particularly suited to open tree structures because it has a "grammar" especially developed for repeated branching. Although the formalism was originally developed for botanical structures, its language carries directly to arterial structures. Following that language and notation ( and production rule p, is given by:
![]() |
(1) |
where F represents a line of unit length in the horizontal direction and X is an auxiliary symbol that has no graphical representation but plays an essential role in the branching process. The square brackets represent the departure from ([) and return to (]) a branch point, whereas the plus and minus signs represent turns through a given angle in the clockwise and anticlockwise directions, respectively. The first four stages of a tree produced by this system, denoted by n = 14, are then given by:
![]() |
(2) |
and the results are shown graphically in Fig 3. Note that the tree is perfectly symmetrical and has uniform properties. Asymmetry and variable properties are introduced later.
|
Fractal Character
The fractal character of the tree produced in the previous section is evident by the fact that its growth from one stage to the next is achieved by the same production rule, namely that defined by p in Equation 1. It is also evident by the self-similarity of the tree structure, namely the fact that any subtree taken from the tree as a whole has the same (dichotomous branching) structure as the whole. This self-similarity property can be expressed more rigorously within the L-system formalism by adopting a somewhat different L-system model for generating the tree, namely:
![]() |
(3) |
The first four stages of the tree are given by:
![]() |
(4) |
and the results are shown in Fig 4. An important difference between this L-system and that of the previous section is that the production rule here, although producing the same tree structure as that produced by Equation 1, involves graphical duplication of some branch segments at each production beyond n = 2. At n = 1, the first junction point is produced and thereafter the graphical "turtle" always returns to that point after tracing each of the two sides of the tree arising from that junction. If each of the duplications is treated as an overwriting of the existing segment, then the resulting tree is graphically the same as that produced by Equation 1.
|
An important advantage of this system, however, is that it provides a more explicit expression of the self-similarity property of the tree by making it possible to define any stage of the tree as an axiom, and then applying the same production rule to that axiom to produce the tree. In particular, taking the second stage of the tree in Equation 4 as an axiom, we obtain the following system:
![]() |
(5) |
The first three stages of the tree are given by:
![]() |
(6) |
and the results are shown in Fig 5. Similarly, taking the third stage of the tree in Equation 4 (or the second stage of the tree in Equation 6) as an axiom, we obtain the following system:
![]() |
(7) |
|
The first three stages of the tree are now given by:
![]() |
(8) |
and the results are shown in Fig 6. Comparison of the three L-systems shows that the production rules in Equation 5 and Equation 7 are identical in form to the production rule in Equation 3. Indeed, this bracketed form is the generator of dichotomy in the three cases.
|
The key difference between these two systems and that of the previous section is that production here proceeds in terms of B or C without any reference back to F once the symbols B and C have been defined. This feature is illustrated graphically in Fig 5 and Fig 6. The lowest building unit here is not F but B or C as the case may be and as illustrated in these figures. The fact that B or C can be any stage of the tree while the production rule remains unchanged provides an algorithmic expression of self-similarity and of fractal character. This rather elegant expression of self-similarity is not possible with Equation 1 since, in that system, F continues to appear explicitly in the production rule.
Although Equation 3 has the problem of duplication or overwriting of branch segments as discussed earlier, this is balanced by the advantage that production here can be "speeded up" considerably by taking a larger and larger subtree as axiom. This can have significant advantage both graphically, as illustrated in Fig 5 and Fig 6, and algorithmically. Indeed, in terms of F, the algorithmic expression for the third stage of the tree in Equation 8 would be the following:
![]() |
(9) |
compared with the much simpler expression in terms of C in Equation 8.
Elements of Arterial Branching
As stated earlier, arterial trees have been found to consist of primarily dichotomous branching, the basic structural unit of the tree being an arterial "bifurcation," whereby an arterial segment divides into two branches and then each of the branches goes on to do the same, and so on (1 and
2, which the two branches make with the direction of the parent segment (Fig 7). In what follows, the convention is used that all properties with subscript 1 refer to the branch with the larger diameter, i.e., the inequality d1 > d2 is always assumed to apply, except, of course, when the two branch diameters are equal. This allows the definition of a useful parameter, sometimes referred to as the bifurcation index or asymmetry ratio:
![]() |
(10) |
which has the convenient range of values 01.0. The index is a measure of the asymmetry of the two diameters at an arterial bifurcation, with = 1.0 when the two diameters are equal, and
0 when one of the two diameters is very much larger than the other.
|
One of the most fundamental laws underlying the structure of arterial trees is that relating the diameter d of an arterial segment to the flow rate q, which the segment is destined to carry. A classical result obtained by d3. The result is based on a compromise between the power required to drive steady flow through the vessel, which varies as d-4, and the rate of expenditure of metabolic energy required to maintain the volume of blood filling the vessel, which varies as d2. Data from the cardiovascular systems of humans and animals have given the law considerable support, though with a great deal of scatter (
dk have been considered, with values of k different from three (
Application of the cube law at an arterial bifurcation provides the basis for its role in arterial branching. Conservation of flow rate at an arterial bifurcation requires that flow rate in the parent vessel equals the sum of flow rates in the two branches (Fig 7), i.e., q0 = q1 + q2. Application of the cube law then converts this into a relation between the diameters of the three vessels, namely d03=d13+d23, and the following two diameter ratios can then be defined in terms of :
![]() |
(11) |
Another important parameter is the ratio of branch length to the length of the parent vessel at an arterial bifurcation. So far, there have been no theoretical grounds for determining this ratio, but biological data indicate that the length of a vessel segment is very weakly related to its diameter in the sense that vessel segments of smaller diameter tend to have smaller lengths, though exceptions occur in sufficiently large number to blur any consistent relation. At best, the data suggest a maximum length to diameter ratio of
35 and an approximate average of 10 (
for the two branches in an arterial bifurcation are then equal to the corresponding values of the diameter ratios, namely
![]() |
(12) |
Finally, to complete the structure of an arterial bifurcation, it is necessary to define the angles, 1 and
2, that the two branches make with the direction of the parent vessel (Fig 7). There are optimum values of these angles that minimize the total volume of the three vessel segments, given by
![]() |
(13) |
Combined with the cube law, the same angles also minimize the pumping power required to drive the flow through the bifurcation. It can be observed from these expressions that, in a nonsymmetrical bifurcation ( < 1.0), the branch with the larger diameter makes a smaller branching angle than does the other branch, i.e.,
1 <
2.
Parametric L-System Branching
The variables of arterial branching ,
,
1, and
2, can be incorporated into the L-system formalism by using a so- called "parametric L-system" (
1) and -(
2) in which the angles of clockwise and anticlockwise rotations are stated explicitly and individually through the values of
1 and
2 .
The basic parametric L-system for a dichotomous tree then takes the form:
![]() |
(14) |
where the tree production now requires that values of the parameters L, W, 1, and
2 be specified. In fact, the basic L-systems considered in previous sections also require these values, but when the parameters are not included explicitly in the algorithm, as in previous sections, it is assumed that their values are fixed throughout the tree structure. Indeed, in the trees shown in Fig 3 Fig 4 Fig 5 Fig 6, the values L = 1, W = 1, and
1 =
2 = 40° were implicitly used, though not stated. Also, it is convenient, particularly for computer-generated images, to allow the scales of L and W to be different from each other. In this way, L = 1, W = 1 above can be taken to mean that L is one unit of length (in centimeters, for example) while W is one unit of width (e.g., in millimeters); thus, the unit of length being different from the unit of width.
One of the most powerful features of parametric L-systems is that the values of these parameters can be changed along a tree structure, and the pattern of that change can be incorporated into the production algorithm. This makes it possible to make branches smaller than parents, both in length and diameter, and to give branches individual branching angles. These changes can be made consistent with the characteristics of arterial branching discussed in the previous section by taking,
![]() |
(15) |
It should be noted that although X is only an auxiliary symbol that has no graphical interpretation, it must nevertheless carry the parameters prescribed by the production rule because these become the parameters of the branches in the next production cycle. Values of L0,W0 in the axiom represent the length and width of the first line segment that must be specified. The lengths and widths of subsequent line segments are then adjusted by the scaling parameters 1,
2,
1, and
2, and their branching angles are determined by
1 and
2. The scaling occurs at each production cycle, and is therefore compounded as production progresses. This can be observed by writing down the first three stages of the tree, namely
![]() |
(16) |
It is to be recalled that, by convention, the parameters 1,
1, and
1 at each bifurcation are associated with the branch that has the larger diameter, whereas
2,
2, and
2 are associated with the branch that has the smaller diameter. Furthermore, in the above Equation 15 and Equation 16, the angle
1 has been associated with an anticlockwise turn (minus sign), whereas the angle
2 has been associated with a clockwise turn (plus sign). As a result of these conventions, the branch with the larger diameter at a bifurcation will always appear on the same side of the parent vessel, namely at an anticlockwise angle. These conventions are purely arbitrary, and can be reversed or changed in a number of ways. The conventions do affect the orientation and appearance of a tree structure, however, and must therefore be recalled when interpreting the resulting tree structure.
Finally, since the scaling parameters 1,
1,
2, and
2 and branching angles
1 and
2 are all functions of the bifurcation index
(Equation 11 HREF="#FD12">Equation 12Equation 13), the above Equation 15 will generate a different arterial tree for each value of
. For
= 1.0, all bifurcations are symmetrical and so is the resulting tree structure (Fig 8 a), although the lengths and diameters of vessel segments are diminishing in accordance with the rules of arterial branching. The effect of small asymmetry is shown in Fig 8 b, where
= 0.8, and the results for
= 0.6, 0.4, and 0.2 are shown in Fig 8ce.
|
![]() |
DISCUSSION AND CONCLUSIONS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Results in Fig 8 (ae) demonstrate clearly that parametric L-systems can generate a wide spectrum of tree structures that have the theoretical properties of arterial branching and that, by construction, have a fractal pattern.
At one end of this spectrum is a tree with symmetrical bifurcations ( = 1.0) in which the diameters of vessel segments diminish rapidly as the tree structure progresses (Fig 8 a). At symmetrical bifurcations, the diameter reduction ratios
1 and
2 are equal and have their maximum value of 0.794 (Equation 11), i.e., vessel diameters are reduced by
20% at each bifurcation. In this case, the tree progresses most "rapidly" toward its ultimate level at which the branches implement the delivery of blood to the capillary bed. Vessels with this branching pattern have been termed "delivering vessels" (
|
At the other end of the spectrum, is a tree with highly nonsymmetrical bifurcations (
0) in which the diameter of the leading branch at each bifurcation is only slightly reduced. In Fig 8 e, where
= 0.2, and hence from Equation 11
1
0.997, the reduction is only 0.3%. In this case, the leading artery retains its identity and much of the flow rate that it carries as the tree progresses. Vessels with this pattern of branching have been termed "distributing vessels" (
Thus, certain arterial branching patterns can be mimicked by fractal L-system trees, using an appropriate choice of branching parameters. However, the L-system trees have one important characteristic that is rarely shared by arterial trees in the cardiovascular system. In each panel of Fig 8, the value of is fixed throughout that particular tree structure (although values of the L-system parameters do change from one level of the tree to the next). This type of uniformity is rarely observed in the physiological system. The vasculature of the kidney shown in Fig 2, for example, is like the fractal tree shown in Fig 8 a in general appearance, but the first clearly lacks the uniformity inherent in the second. The primary source of this nonuniformity in arterial trees is the value of the bifurcation index
. Data from the vascular system have shown consistently that
takes on almost the full range of values (01.0) at each level of a tree and from one level to the next (
as shown in Fig 10.
|
Variability in the value of the asymmetry ratio, , can be incorporated into the L-system by assigning a random value to
between 0 and 1.0 at each branch point. Two examples of the resulting trees are shown in Fig 11, which raise two important questions. First, is the wide variability in the value of
observed in physiological data actually random, or is it determined by local conditions and other constraints that, therefore, cannot be modeled by purely random variation? Second, is it legitimate to refer to the trees in Fig 11 and to arterial trees that they mimic as fractal structures? Although the trees in Fig 11 are generated by an L-system in which the production rule has the same form at all levels of the tree, the values of the parameters involved are not predetermined. The second stage of one of the trees in Fig 11, for example, cannot be used to generate the third and subsequent stages of that tree as was done in Fig 5 and Fig 6. Indeed, the second stage itself cannot be regenerated because it would be randomly different at each production.
|
The question of whether arterial trees have fractal structures has important bearing on the more general question of scaling in biology (. Although a constant value of
near 0 or 1.0 produce branching patterns that have been observed in the physiological system, a random choice of
does not produce the pattern of variability found in the physiological system.
In the analogous problem of botanical tree structures, probabilities have sometimes been used to determine the value or nature of a branching variable in the production process ( a certain probability of occurring at each branch point. Data from the cardiovascular system have so far not shown any basis for assigning such probabilities. The choice of random values of
at each branch point is at least consistent with the range of variability shown by the data from the cardiovascular system, even if in the physiological setting the source of that variability is not known. Here, the values of
may be influenced by local anatomy, local flow requirements, and other constraints, and therefore the variability may not be purely random. On the other hand, in the botanical situation, there are indications that probabilities may well play a role in the growth of tree structures (
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
This work was supported by the Natural Sciences and Engineering Research Council of Canada.
Submitted: 13 February 2001
Revised: 6 July 2001
Accepted: 9 July 2001
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Hutchins, G.M., Miner, M.M., and Boitnott, J.K. 1976. Vessel caliber and branch-angle of human coronary artery branch-points. Circ. Res. 38:572-576[Abstract].
Horton, R.E. 1945. Erosional development of streams and their drainage basins: hydrophysical approach to quantitative morphology. Bull. Geol. Soc. Am. 56:275-370.
MacDonald, N. 1983. Trees and Networks in Biological Models. Chichester, John Wiley & Sons Inc, pp. 215.
Mandelbrot, B.B. 1977. Fractals: form, chance, and dimensions. San Francisco, Freeman, pp. 365.
Mayrovitz, H.N., and Roy, J. 1983. Microvascular blood flow: evidence indicating a cubic dependence on arteriolar diameter. Am. J. Physiol. 245:H1031-H1038
Moffat, D.B. 1979. The anatomy of the renal circulation. In Black D., Jones N.F, eds. Renal Disease. Oxford, Blackwell.
Murray, C.D. 1926a. The physiological principle of minimum work. I. The vascular system and the cost of blood volume. Proc. Natl. Acad. Sci. 12:207-214.
Murray, C.D. 1926b. The physiological principle of minimum work applied to the angle of branching of arteries. J. Gen. Physiol. 9:835-841
Nishida, T. 1980. KOL-simulating almost but not exactly the same development. Mem. Fac. Sci. Kyoto Univ. Ser. Biol. 8:97-122.
Peitgen, H.O., Jurgens, H., and Saupe, D. 1992. Chaos and Fractals: New Frontiers of Science. New York, Springer-Verlag, pp. 984.
Pollanen, M.S. 1992. Dimensional optimization at different levels of the arterial hierarchy. J. Theor. Biol. 159:267-270[Medline].
Prusinkiewicz, P., and J. Hanan. 1989. Lindenmayer Systems, Fractals, and Plants. Lecture Notes in Biomathematics. Vol. 79. Springer-Verlag, New York. 120 pp.
Prusinkiewicz, P., and Lindenmayer, A. 1990. The Algorithmic Beauty of Plants. New York, Springer-Verlag, pp. 228.
Rodbard, S. 1975. Vascular caliber. Cardiology. 60:4-49[Medline].
Roy, A.G., and Woldenberg, M.J. 1982. A generalization of the optimal models of arterial branching. Bull. Math. Biol. 44:349-360[Medline].
Sernetz, M., Gelleri, B., and Hofman, F. 1985. The organism as a bioreactor, interpretation of the reduction law of metabolism in terms of heterogeneous catalysis and fractal structure. J. Theor. Biol. 117:209-230[Medline].
Sherman, T.F. 1981. On connecting large vessels to small. The meaning of Murray's Law. J. Gen. Physiol. 78:431-453[Abstract].
Weibel, E.R. 1984. The Pathways for Oxygen. Cambridge, MA, Harvard University Press, pp. 425.
West, G.B., Brown, J.H., and Enquist, B.J. 1997. A general model for the origin of allometric scaling laws in biology. Science. 276:122-126
Woldenberg, M.J., and Horsfield, K. 1983. Finding the optimal length for three branches at a junction. J. Theor. Biol. 104:301-318[Medline].
Zamir, M. 1978. Nonsymmetrical bifurcations in arterial branching. J. Gen. Physiol. 72:837-845[Abstract].
Zamir, M. 1988. Distributing and delivering vessels of the human heart. J. Gen. Physiol. 91:725-735[Abstract].
Zamir, M. 1999. On fractal properties of arterial trees. J. Theor. Biol. 197:517-526[Medline].
Zamir, M., and Brown, N. 1982. Arterial branching in various parts of the cardiovascular system. Am. J. Anat. 163:295-307[Medline].
Zamir, M., and Phipps, S. 1987. Morphometric analysis of the distributing vessels of the kidney. Can. J. Physiol. Pharmacol. 65:2433-2440[Medline].
Zamir, M., Sinclair, P., and Wonnacott, T.H. 1992. Relation between diameter and flow in major branches of the arch of the aorta. J. Biomech. 25:1303-1310[Medline].
Zamir, M., Wrigley, S.M., and Langille, B.L. 1983. Arterial bifurcations in the cardiovascular system of a rat. J. Gen. Physiol. 81:325-335[Abstract].