Abstract

Irregular boundary lines can be characterized by fractal dimension, which provides important information for spatial analysis of complex geographical phenomena such as cities. However, it is difficult to calculate fractal dimension of boundaries systematically when image data are limited. An approximation estimation formula of boundary dimension based on square is widely applied in urban and ecological studies. But the boundary dimension is sometimes overestimated. This paper is devoted to developing a series of practicable formulae for boundary dimension estimation using ideas from fractals. A number of regular figures are employed as reference shapes, from which the corresponding geometric measure relations are constructed; from these measure relations, two sets of fractal dimension estimation formulae are derived for describing fractal-like boundaries. Correspondingly, a group of shape indexes can be defined. A finding is that different formulae have different merits and spheres of application, and the second set of boundary dimensions is a function of the shape indexes. Under condition of data shortage, these formulae can be utilized to estimate boundary dimension values rapidly. Moreover, the relationships between boundary dimension and shape indexes are instructive to understand the association and differences between characteristic scales and scaling. The formulae may be useful for the prefractal studies in geography, geomorphology, ecology, landscape science, and especially, urban science.

1. Introduction

Fractal systems can be characterized by fractal dimension, and the basic and important approach to understanding fractal dimension is the geometric measure relations. Euclidean geometric measure relations come from the principle of dimension consistency. A measure (e.g. length) is not proportional to another measure (e.g. area) unless they share the same spatial dimension [13]. From the principle of dimensional homogeneity, we can derive Euclidean geometric measure relations, which can be generalized to fractal geometric measure relation [35]. From fractal measure relations, we can derive fractal dimension and allometric scaling exponents [3, 68]. An allometric scaling relation can be regarded as a generalized fractal measure relation. Among various geometric measure relations, the common one is the area-perimeter scaling relation, which was used to obtain the boundary dimension of self-similar shapes embedded into a 2-dimensional space [3, 4]. In urban studies, the form dimension and boundary dimension can be derived from the fractal measure relations [9, 10]. Form dimension of cities includes box dimension and radial dimension, which are defined on the basis of the relations of urban area and the linear sizes of box or radius of concentric circles [9, 11]. This work focuses on boundary dimension, which can be associated with form dimension in theory.

If we have enough data, we can calculate various fractal parameters. Taking urban research as an example, we can research spatial distribution and structure using box dimension and research urban growth using radial dimension. If there is not enough information for urban morphology, the boundary dimension of a city can be calculated [9, 12, 13]. However, sometimes, we only know the urban area and urban envelope. Urban envelope represents closed urban boundary lines, and urban area represents the region within the boundary curve [9, 14]. In this case, we can only estimate the boundary dimension by means of the datasets of urban area and perimeter length. Referring to a Euclidean shape, we can construct a series of estimation formulae of fractal dimension. The basic reference shapes are standard circle and square. The formula of boundary dimension based on square has been constructed by Olsen et al. [15], and the formula was widely applied in urban and ecological studies [16, 17]. The formula is simple and easy to understand and has strong practicability. However, it has two drawbacks. First, this formula is mainly applicable to objects extended in the form of squares. Second, the formula sometimes overestimates the boundary dimension [8, 18]. Therefore, we need not only the fractal dimension estimation formulae based on other reference shapes but also a new fractal dimension estimation formula with reference square. This paper is devoted to deriving two sets of approximation formulae of fractal dimension estimation for the fractal-like boundary dimension of irregular shapes. The reference figures include regular triangle, square, regular hexagon, and standard circle. As a contrast model, the generator of Koch snowflake curve is also employed as one of the reference shapes. From the geometric measure relations based on these reference shapes, a series of formulae are derived to approximately estimate the boundary dimension of various irregular shapes such urban envelopes.

2. Fractal Measure Relations

2.1. The First Set of Formulae

The so-called fractal dimension values based on the approximation formulae are actually fractal indicators, which can be used to replace fractal dimension under the condition of absence adequate data. A basic postulate is that the boundary line of an irregular region is a closed prefractal curve. A prefractal is a fractal-like object, which is not a real fractal [19, 20]. The “length” of true fractal line is infinite. If the irregular boundary such as urban envelope is a real fractal curve, we cannot derive any simple formula for fractal indicators and shape indexes. Using the ideas from prefractals, we can find a number of approximation formulae of boundary dimension from a given reference shape. The reference shapes are some types of regular geometric figures defined in a 2-dimensional Euclidean space, including standard circle, regular triangle, square, regular hexagon, and regular six-pointed star (Figure 1). Triangle can be regarded as the basic shape in Euclidean geometry. All geometric figures, including squares, rectangles, trapezoids, circles, ellipses, and irregular shapes, can be reduced to triangles. So, the first formula of fractal dimension estimation for irregular boundaries should be derived from a regular triangle. For an equilateral triangle with a side length r, the area A and perimeter P can be expressed as follows:

Thus, the geometric measure relation between area A and perimeter P can be obtained by combining equation (1) and equation (2). Eliminating the side length r yields

Suppose that the three Euclidean sides are replaced by three fractal lines and the fractal dimension of these lines is D. In this case, a regular shape changes to an irregular shape (Figure 2). Thus, the Euclidean geometric measure relation, equation (3), should be substituted by a fractal geometric measure relation:where D denotes the fractal dimension of the boundary line. From equation (4), a formula of estimating the fractal dimension of boundary lines can be derived as below:

If the shape of a natural system such as a city is similar to a triangle or a system has three growing directions, the fractal dimension of the system’s boundary line can be estimated by equation (5).

The second formula can be constructed on the basis of square. A square is simple and regular, and it is easy to calculate its area A and perimeter P if the side length r is known. The area and perimeter formulae are as follows:

Combining equations (6) and (7) yields the geometric measure relation between the area A and perimeter P as below:

If the sides of the square are replaced by the fractal lines with fractal dimension D, the geometric measure relation, equation (8), will be replaced by

From equation (9), a fractal dimension estimation formula can be derived as follows:

This formula is familiar to many scholars who like geographical and ecological fractals because it was once derived by Olsen et al. [15] in another way. If the shape of a natural system is similar to a square, or a system has four growing directions, the boundary dimension of the system shape can be estimated by equation (10).

The regular hexagons can be best closed to each other in a geographical region. Therefore, the hexagonal networks were applied to the well-known central place theory [21]. The area A and perimeter P of a regular hexagon with a side length r can be calculated by the following formulae:

From the equations (11) and (12), we can derive a geometric measure relation such as

Substituting the Euclidean sides of the hexagon with fractal boundary lines, we can turn equation (13) into a fractal measure relation as follows:

From equation (14), a fractal dimension estimation formula can be obtained as below:

If the shape of a natural system is similar to a hexagon or a system has six growing directions, the boundary dimension of the system can be estimated by equation (15).

The standard circle is treated a simple and perfect shape in Euclidean geometry. Many shape indexes of geography are based on this kind of circle [2225]. The area A and perimeter P of a circle with a radius r can be given by

Integrating equation (16) into equation (17) yields the geometric measure relation between the circular area A and circumference P as follows:

Replacing the Euclidean perimeter with a fractal curve results in a fractal measure relation as below:

Thus, the fractal dimension of the boundary can be calculated by the following formula:

If the shape of a natural system is similar to a circle or if a system is of isotropic growth, the boundary dimension of the system can be estimated by equation (20).

All the above-given formulae are based on Euclidean figures. For comparison, it is advisable to construct a calculation formula based on fractal generators. Koch snowflake curve is one of classical models for fractal lines. We can design the formula using the generator of Koch snowflake curve, a regular six-pointed star. For fractal generator of Koch snowflake curve with a side length r, the area A and perimeter P are as follows:

Thus, the geometric measure relation between the area A and perimeter P can be derived as

The generator of Koch curve is not a fractal line, but the second step is a prefractal figure. Substituting the straight line segments with fractal lines yields a fractal measure relation as below:in which D refers to the fractal index of irregular curve. The formula of fractal dimension estimation based on the Koch snowflake generator can be derived from equation (23) as follows:

If the shape of a natural system is similar to a Koch snowflake or a system has six protruding growth directions, the boundary dimension of the system can be estimated by equation (24).

2.2. The Second Set of Formulae

The aforementioned formulae are suitable for the very irregular boundary lines; otherwise, the fractal dimension may be overestimated. In order to estimate the boundary dimension of general fractal-like line, we should improve the formulae. Based on regular triangles, equation (4) can be replaced by the following fractal measure relation:which can be re-expressed as

From equations (25) and (26), a triangle-based fractal dimension formula can be derived as below:

Compared with equation (5), equation (27) can give more realistic dimension values of fractal lines. Based on squares, equation (9) can be substituted bywhich can be rewritten as

From equation (29), a square-based fractal dimension formula can be derived as follows:

Compared with equation (10), equation (30) can yield more realistic values of boundary dimension. Based on regular hexagons, equation (14) can be replaced withwhich is equivalent to

From equations (31) and (32), a hexagon-based fractal dimension formula can be derived as follows:

Compared with equation (15), equation (33) can produce more realistic values of boundary dimension. Based on the standard circle, equation (19) can be replaced bywhich can be converted into

From equation (35), a circle-based fractal dimension formula can be derived as below:

Compared with equation (20), equation (36) can give more realistic fractal dimension values of boundary lines. Based on the regular six-pointed star, equation (23) can be substituted withwhich can be converted into

From equation (38), a Koch-snowflake-based fractal dimension formula can be derived as below:

Compared with equation (24), equation (39) can give more realistic values of fractal indexes of boundary lines.

The two sets of fractal dimension estimation formulae represent two sets of fractal boundary indexes. For the convenience of readers, the two sets of formulae are tabulated as follows (Table 1). Applying these formulae to a simple irregular shape (Figure 2), we can obtain two sets of fractal indexes values for boundary lines (Table 2). Based on the first set of formulae, the fractal dimension estimation results are marked as boundary dimension D(1), and based on the second set of formulae, the fractal dimension estimation results are marked as boundary dimension D(2). The values of D(2) are less than those of D(1).

2.3. Relations of Boundary Dimension to Shape Indexes

The area-perimeter measure relation is in essence a problem of scaling in complex systems. The formulae of boundary dimension represent various definitions of the scaling exponents of different shapes. The traditional thinking of mathematical modeling and quantitative analysis is based on characteristic scales [5, 2629]. Therefore, a number of shape indexes have been derived from the area-perimeter relations to describe natural morphology such as urban patterns [2225]. The typical shape indexes are the circularity ratios, which are defined on the basis of area and perimeter [23]. However, complex systems such as cities have no significant characteristic scales in many aspects [26]. In this case, the scaling concept is employed to substitute the notion of typical scales such as characteristic lengths. Correspondingly, we utilize scaling exponents instead of shape indexes to characterize the form features of irregular patterns. The scaling exponents are based on the ideas from fractal geometry, while the shape indexes are based on the notion of Euclidean geometry. Despite the difference between the scaling exponents and shape indexes, there is inherent association of boundary dimension with shape indexes. The relations between the circularity ratios and the reciprocal of the boundary dimension have been proved to be the exponential function [22]. Based on area and perimeter, a series of shape indexes similar to the circularity ratios can be derived from the aforementioned geometric measure relations, which are listed in Table 1. Applying these formulae to Figure 2 yields a series of values for shape indexes. The relations between the new shape indexes and the second set of boundary dimension can be demonstrated to satisfy the following exponential function:where the parameters are a = 1/P2 and b = 2 ln (P). For example, for Figure 2, the relationship between the second type of boundary dimension and the shape index is as follows:which can be verified by the observational data (Table 2). It is easy to testify that a = 1/295.1157^2 = 0.00001148 and b = 2 ∗ ln (295.1157) = 11.3747.

If the boundary line is a true fractal line, then the perimeter will be infinite. In this case, equation (40) will be invalid. This suggests that the fractal dimensions and shape indexes in Table 1 will be invalid for true fractal boundary lines. In the real world, the boundary lines are fractal-like lines rather than true fractal lines. A real fractal bears no scaling limitation, and the geometric measure relation can be reflected by a straight line of infinite length on a log-log plot. In contrast, a prefractal possesses fractal nature only within certain scaling range. Therefore, the perimeters of prefractal closed curves bear certain values, and thus the formulae derived above are valid. Moreover, equation (40) lends further support to the inference that the relationships between the reciprocal of boundary dimension and the circularity ratios meet an exponential function. However, the relations between the shapes indexes and the first set of boundary dimension cannot be described with the exponential function exactly. Despite this, the inherent correlation between characteristic scales and scaling can be reflected by the relationships between boundary dimension and shape indexes.

3. Empirical Analysis

3.1. Fractal Dimension Estimation Results

To show the effects of the two sets of fractal dimension estimation formulae, we can apply them to the cities in Beijing, Tianjin, and Hebei region, China (for short, Jing-Jin-Ji region). There are 13 main cities in the study area. As a preparation, it is necessary to extract urban boundary lines using a proper method. In urban geography, the boundary curve of a city is termed urban envelope, and the region within the urban envelope is termed urban area [9, 14]. There are at least four scientific approaches to identifying and delineating urban envelopes for these cities [30], including the city clustering algorithm (CCA) [31], the automatic identification method of urban settlement boundaries [32], the fractal-based method [33], and the approach to derive “natural cities” by clustering street nodes/blocks [34]. In this paper, CCA is employed to delineate urban boundary lines on interpreted remote sensing images in different years. The urban boundary determined by this method corresponds to an urban agglomeration, which approximately corresponds to an urbanized area. The urban envelopes give urban perimeters, and the corresponding urban areas can be counted (Table 3). These datasets can be used to verify the geometric measure relation between urban area and perimeter and evaluate the fractal boundary indexes [8, 30].

It is easy to calculate the boundary fractal dimension of the 13 cities in Beijing, Tianjin, and Hebei region using each formula. As indicated above, the fractal dimension estimation values based on the first set of formulae are marked as boundary dimension D(1), and the fractal dimension estimation values based on the second set of formulae are marked as boundary dimension D(2). Based on the first set of formulae, the estimated values of the boundary fractal dimension are sometimes very high or very low or even greater than 2 or less than 0. In contrast, if we utilize the second set of formulae to estimate the boundary dimension, the results are relatively reasonable. All the values range from 1 to 2. In other words, no value is greater than 2 or less than 0. The results of 2000, 2005, and 2010 are listed in Tables 46. The shape indexes are listed separately for reference (Table 7). The boundary fractal dimension based on the first set of formulae has weak correlation with the shape index. However, there is negative correlation between the fractal dimension values based on the second set of formulae and the shape indexes (Supplementary File 1).

Several points of explanations should be provided for the exceptional values in the fractal dimension estimation results. In theory, the boundary dimension defined in a 2-dimensional embedding space is supposed to come between 0 and 2 [3, 9, 12, 13, 35, 36]. The reasonable values vary from 1 to 1.5. However, the following causes often result in overestimation of boundary dimension. First, the boundary dimension calculation is based on the geometric measure relation deriving from regular real fractals in the mathematical world. A real fractal has no scaling range, or, in other words, the scaling range of a real fractal is infinite. Applying the fractal measure relations proceeding from regular real fractals to the random prefractals gives rise to significant bias in many cases [18]. Second, if the image resolution is high enough, the length of the boundary line may be very long, but the area within the boundary curve is certain. This phenomenon can be illustrated by a regular fractal termed Koch lake (see Appendix). The Koch lake can be treated as models of lakes, islands, urban region, and so on. Third, compared with the second set of approximate formulae, the first set of approximate formulae enlarges the ratio of the perimeter logarithm to the area logarithm of a shape relatively. For example, for the formulae based on square, in equation (10), the circumference is reduced to a quarter of the original length, while in equation (30), the area is enlarged to 16 times the original one. As a result, the boundary dimension value of a shape based on equation (10) is significantly greater than the value based on equation (30). Generally speaking, the estimated value of a boundary dimension is not less than 1. However, if a figure is near a Euclidean shape, the fractal dimension estimation result may be an outlier and less than 1 because the formulae are designed for prefractals rather than for Euclidean shapes.

3.2. Fractal Dimension Adjustment and Transformation

Fractal lines fall into two types: one is boundary lines and the other is what is called space-filling curves. In theory, the well-known boundary line is Koch snowflake curve, and the well-known space-filling curve is Peano curve [3]. In reality, fractal boundary lines include coast lines, urban boundaries, lake boundaries, and national boundaries, and space-filling curves include rivers, traffic networks, and hierarchical boundaries in central place systems. Generally speaking, the average values of fractal dimension of boundary lines come between 1 and 1.5, while the fractal dimension values of space-filling curves vary from 1.5 to 2. If a fractal dimension of urban boundary is overestimated, it can be adjusted by a simple formula. Furthermore, the adjusted boundary dimension can be converted into the fractal dimension of urban form. The adjusted boundary dimension is as follows [18]:where Db denotes the adjusted boundary dimension and Dl refers to the originally estimated fractal dimension (note: the subscript for Dl is the Latin letter l, not the Arabic numeral 1). Using this formula, we can adjust the boundary dimension values (Table 8). The boundary dimension and the form dimension satisfy a hyperbolic relation as below [18]:where Df refers to the form dimension of a city. Using this formula, we can transform the adjusted boundary dimension values into the fractal dimension of urban form (Table 9).

The approximate fractal dimensions are in fact fractal indexes. An index can condense many data into a number, describing the characteristics of a system and simplifying the analytical process. Given the area and perimeter of an irregular shape such as an urban envelope, we can calculate its fractal indexes using these sets of formulae. All these approximate formulae can be applied to estimation of boundary dimension of the cities in the Yangtze River Delta, China. The datasets in 1985, 1996, and 2005 have been published [8]. The results show that the fractal index values of city boundaries in the Yangtze River Delta are significantly lower than those in Beijing, Tianjin, and Hebei region. The reason may be that the resolution of remote sensing images of Beijing, Tianjin, and Hebei cities is higher than that of Yangtze River Delta cities (Supplementary File 2).

4. Discussion

The aforementioned results are based on power-law relations, and a power law represents a geometric measure relation and reflects a proportional relationship. A power function has two parameters: one is the proportionality coefficient and the other is the power exponent. In the framework of Euclidean geometry, the power exponent is always a known constant and bears little useful information. Thus, we can construct various shape indexes based on proportionality coefficients. A proportional constant is always a dimensionless parameter reflecting a ratio of one measure to another measure. On the contrary, in the framework of fractal geometry, the proportionality coefficient bears little information, but the power exponent is unknown parameter and possesses spatial information. A simple system has characteristic scale and can be described with the mathematical method based on Euclidean geometry, while a complex system has no characteristic scale and cannot be effectively described by conventional mathematical methods. In this case, it is necessary to replace the characteristic scale concept with scaling idea. The power exponent is known as scaling exponent. Fractal geometry is a powerful tool for scaling analysis of complex systems, and the fractal dimension is an important scaling exponent. Based on the notion of fractals, various fractal indexes can be defined to characterize fractal-like phenomena.

A set of formulae for estimating the boundary dimension of irregular shapes has been derived from the geometric measure relations between the area and perimeter of certain reference figures. The reference shapes include regular triangle, square, regular hexagon, circle, and the generator of Koch snowflake curve. From difference reference shape, we can obtain different formulae; from the same reference shape, we can obtain at least two different formulae based on different conditions (Table 1). Different formulae have different spheres of application. In practice, we should select the proper formula according to the shape and irregularity of studied objects. As a matter of fact, a number of methods can be used to calculate the boundary dimension. The common methods include divider method [37, 38], box counting method [39, 40], and so on. Longley and Batty [12, 13] developed four methods to measure the fractal dimension of fractal lines. Given enough image data, it is not problematic to calculate the boundary dimension [6, 8, 4143]. The formulae proposed in this paper are suitable for fractal dimension estimation of irregular boundaries under the condition of data shortage. Concretely speaking, we need to use these formulae in three cases. (1) Limitation of data: the amount of data is small, and the existing data do not support the calculation of fractal dimension by the least squares calculation. (2) Approximation of results: an approximate estimation of fractal dimension can meet the needs of special research. (3) Comparability of datasets: image data of different years or places have the same quality. For example, we do not have any data except the numbers of boundary lengths of a city’s and the area within the boundary lines in different years. In this instance, we can estimate the boundary dimension of the city and analyze its growing process and pattern. The results of fractal dimension estimation are not real fractal dimension, but a kind of characteristic indexes to describe the fractal-like boundaries.

A fractal measure relation can be treated as an allometric scaling relation. These scaling relations are widely applied to urban research. In fact, the geometric measure relation between urban area and urban perimeter bears analogy with the allometric scaling relation between urban population and urban area, which can be expressed aswhere S denotes the population size of a city. Thus, the allometric scaling relation between urban population size S and urban boundary length P can be derived as follows:which can be reexpressed aswhere μ = SPα represents the proportionality coefficient and α = Dp/D is the scaling exponent. This suggests that urban population size is in a proportion to α power of the urban perimeter.

Further, the geometric measure relation can be generalized to the traffic network of a city. Suppose that the urban area is A and the total length of traffic lines is L. According to the principle of dimension consistency, we havewhere Dw denotes the fractal dimension of traffic networks. For comparability, the proportionality coefficient is assumed to be 1, and then equation (47) can be transformed into

From equation (48), a fractal index can be derived as below:

This is the formula of fractal dimension of traffic network in an urbanized area. This also implies that the effect of measures’ dimension on the result can be eliminated by taking logarithms of measures when constructing a fractal index.

Fractal measures are significant in the research on complex landscape, which bear no characteristic scale and cannot be characterized by the common indexes in theory. Cities represent complex human landscape. The shape indexes and boundary dimension are basic measures in urban studies, and both characteristic scales and scaling are important concepts in urban geography. Comparably speaking, scaling concept is more important. Cities are complex spatial systems, and many aspects of urban systems have no characteristic scales. Scaling in cities has attracted more and more attention of scholars [26, 4451]. Fractal geometry is one of the powerful tools in scaling analysis and has been applied to urban studies, which resulted in a number of interesting achievements. The series of approximation formulae of boundary dimension provide simple approaches to scaling analysis of cities. The shortcomings of this study lie in two aspects. First, the formulae of boundary dimensions and the corresponding shape indexes have not yet been derived from regular pentagon. Compared with the regular triangle, square, regular hexagon, and standard circle, the regular pentagon is more complex because it can be associated with fractals. Second, the systematic positive study has not been made. The empirical evidences shown in this paper are only for methodology rather than for urban studies. Furthermore, the formulae are only suitable for random prefractals rather than real fractals (Appendix). Due to the limitation of space of the paper, the pending questions will be answered in future studies.

5. Conclusions

Boundary line represents a perspective of spatial patterns and landscape analysis of complex systems. In this paper, various possible formulae for estimating boundary dimension are systematically studied and compared. The aim of this work is to provide two sets of practical formulae for approximate estimation of boundary dimension of fractal-like phenomena. The main conclusions can be reached as follows. First, the approximate formulae of boundary dimension can be derived from geometric measure relation in light of scaling thinking. Traditional mathematical modeling and quantitative analysis are based on characteristic scales, and a number of shape indexes are derived from these relations to describe various shapes; complex systems bear no characteristic scales, so a number of scaling exponents are derived to characterize various patterns. Second, the approximate estimation formulae of boundary dimension are not unique, but diverse. On the one hand, different fractal boundary dimension can be defined based on different reference shapes; on the other hand, under different conditions, different formulae can be derived from the same geometric measure relation. Therefore, in practice, proper approximate formula should be selected according to the shape and irregularity of natural morphology. Third, the approximate boundary dimension are essentially scaling exponents for describing complex curves. These formulae are used to estimate fractal parameters of boundary lines only in case of data shortage. The estimation results using the approximation formulae are actually fractal indexes instead of the calculated values of real fractal dimension. These formulae are suitable for prefractal curves rather than real fractal lines. If we have enough image data, we should use normal methods to calculate the fractal dimension with higher confidence level. Sometimes, even if the data are sufficient, these formulae can be used to estimate the fractal parameters quickly when the accuracy requirement is not so high.

Appendix

Fractal Indexes of Koch Lake

The approximate formulae developed in this paper is suitable for fractal-like curves rather than real fractal lines. Let us apply the fractal index formulae to a regular fractal, i.e., Koch lake, which comprises of three Koch curves (Figure 3). Suppose that the initiator of each Koch curve is a line segment of unit length, namely, L1 = 1. Thus, the initiator of Koch lake is an equilateral triangle with a perimeter P1 = 3L = 3. Correspondingly, the area of the triangle is . According to the knowledge of geometry and trigonometry, the perimeter of Koch lake can be expressed aswhere m = 1, 2, 3, … denotes steps of fractal generation. Clearly, if m ⟶ ∞, then Pm ⟶ ∞. In contrast, the area of Koch lake is limited. In light of the knowledge of geometry and trigonometry, the area of Koch lake iswhere m = 1, 2, … refers to steps of fractal development. Under the limit condition, the area of Koch lake approximates a constant, that is,

This is one of the characteristics of fractals: infinite filling in a finite space. In this case, regardless of the formula in Table 1, the fractal dimension is a variable dependent on m instead of a constant. For example, based on square, the fractal indexes and shape index are as follows:

This suggests that both the approximate fractal dimension formulae and shape index are meaningless for real fractal curves.

Data Availability

All data, models, and code generated or used during the study are provided within the article and supplementary information files.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This research was sponsored by the National Natural Science Foundation of China (grant nos. 41590843 and 41671167). The support is gratefully acknowledged. The author would like to thank Ms Yuqing Long of Peking University for providing the essential data on China’s cities.

Supplementary Materials

File 1: datasets for the fractal indexes based on urban area-perimeter relations of Beijing-Tianjin-Hebei region. File 2: datasets for the fractal indexes based on urban area-perimeter relations of Yangtze River Delta. (Supplementary Materials)