Abstract

The area-perimeter allometric scaling is an important approach for researching fractal cities, and the basic ideas and models have been researched for a long time. However, the fractal parameters based on this scaling relation have not been efficiently utilized in urban studies. This paper is devoted to developing a description method of urban evolution using the fractal parameter sets based on the area-perimeter measure relation. The novelty of this methodology is as follows: first, the form dimension and boundary dimension are integrated to characterize the urban structure and texture; second, the global and local parameters are combined to characterize an urban system and individual cities; third, an entire analytical process based on the area-perimeter scaling is illustrated. Two discoveries are made in this work: first, a dynamic proportionality factor can be employed to estimate the local boundary dimension; second, the average values of the local fractal parameters are approximately equal to the corresponding global fractal parameters of cities. By illustrating how to carry out the area-perimeter scaling analysis of Chinese cities in Yangtze River Delta in the case of remote sensing images with low resolution, we propose a possible new approach to exploring fractal systems of cities.

1. Introduction

A scientific study is always involved with two processes: the first process is to describe how a system works using mathematics or measurement, and the second process is to gain understanding of the causality behind the system using observation, experience, or artificially constructed experiments [1, 2]. Fractal geometry is a powerful tool of describing urban systems because fractal dimension is a kind of characteristic parameter to describe scale-free phenomena [39]. By fractal dimension description, we can get insight into the spatial dynamics of urban evolution [1013]. There are two simple and convenient ways of understanding city fractals. One is measurement process (see, e.g., [1419]), and the other is the geometric measure relation (see, e.g., [12, 20, 21]). If a measurement result such as boundary length, urban area, and land-use density depends on the scale to measure (the length of yardstick), we may face a fractal object. On the other hand, if the proportional relation between two measures such as urban area and urban perimeter is involved with a scaling exponent indicating a pair of fractional parameters, we will meet a fractal phenomenon. The former way can be used to realize any types of fractals, while the latter way can be employed to identify the fractal boundary of a region. Sometimes an urbanized area is treated as a Euclidean space, but the urban boundary can be regarded as a fractal line. In this case, the boundary dimension of a city provides a good way of looking at city size and shape [4]. So far, we have more than five empirical approaches to estimating the fractal dimension of geographical boundaries [2227].

An urban area bears an analogy to a random Koch island; thus, the fractal dimension of its boundary can be estimated with the area-perimeter scaling [4, 22, 28, 29]. The geometrical measure relation between urban area and perimeter is in essence an allometric scaling relation, which is similar to the relation between urban area and population [30]. If the population size of a city is compared to the weight of an animal, then the urban area can be compared to the volume of the animal and the urban perimeter to the surface area (the area of the whole skin) [20]. In many cases, it is difficult to investigate city population, but it is easy to measure urban area and the corresponding perimeter by using the remote sensing images and the geographical information system (GIS) techniques. The area-perimeter allometric scaling is a simple approach to revealing the spatiotemporal evolution dynamics of urban systems. However, two problems remain to be solved. First, in many cases, it is impossible to compute the fractal dimension of each city in an urban system because of inadequate remote sensing data or the lower resolution of remote sensing images. Second, the traditional formula of the boundary dimension based on the area-perimeter scaling is not exact enough to guarantee the effective results.

The processes and patterns of urban evolution follow scaling laws [3134]. Fractal geometry is an effective tool to explore scaling in cities. In recent years, the new formulae for revising the boundary dimension calculations through the area-perimeter scaling have been proposed [23]. By means of these formulae, we can correct the errors of the boundary dimension and estimate the form dimension of cities. Using the boundary dimension and form dimension, we can characterize urban shapes and the spatiotemporal evolution of urban systems. This paper is devoted to developing the description method of urban patterns and processes using the fractal dimension sets based on the area-perimeter allometric scaling. It tries to solve several problems as follows: how to combine the global fractal parameters with the local fractal parameters for spatial analysis; how to make use of remote sensing images of low resolution for urban fractal studies; how to understand the influence of urban sprawl on fractal dimension change. The rest of the work is organized as follows. In Section 2, several new formulae of fractal dimension estimation are introduced and clarified for spatiotemporal analysis of urban form and growth. In Section 3, the new models and formulae are applied to the system of cities and towns in Yangtze River Delta, China, to make an empirical analysis. In Section 4, based on the mathematical models and the empirical case, the new analytical process of urban evolution based on the urban area-perimeter scaling is presented and illustrated, and several related questions are discussed. Finally, the paper reaches its conclusions by outlining its major viewpoints.

2. Fractal Parameters

2.1. Basic Formulae

In theory, a city figure can be divided into two parts: one is the urban boundary, and the other is the urban area within the boundary. The former is termed urban envelope () and can be described with the boundary dimension [21], while the latter is named urban area () and can be characterized by the form dimension [4, 35]. The form dimension is a structural dimension, while the boundary dimension is a textural dimension [20, 36, 37]. The form dimension and boundary dimension compose the shape dimension of cities. In fact, the family of shape indexes includes the well-known form ratio, which is defined by area and perimeter of a city [38]. The form dimension can be measured by the area-scale relation, while the boundary dimension can be determined by the perimeter-scale relation [23]. In this sense, urban shape can be characterized by both the form and the boundary dimensions.

Based on digital maps or remote sensing images, fractal cities are always defined in a 2-dimensional space. In practice, the urban area can be regarded as a Euclidean plane with a dimension , and accordingly, the urban boundary is treated as a fractal line [4, 23]. Thus, the boundary dimension represents the basic fractal dimension of urban shape and can be estimated with the regression analysis based on the method of the ordinary least squares (OLS). By the fractal measure relation [28, 29, 39], the urban area and perimeter follow a power law as below:where refers to the perimeter of the urban envelope and refers to the corresponding urban area. As for the parameters, denotes the boundary dimension of traditional meaning and is related with the proportionality coefficient. In this context, the boundary dimension should be termed initial boundary dimension because it represents the conventional concept of the fractal dimension of urban boundary. Equation (1) is in fact an allometric scaling relation, and the scaling exponent iswhere denotes the Euclidean dimension of the embedding space of urban form. The boundary dimension is often estimated by the formula . From (1), it follows thatwhich is an approximate formula of the boundary dimension estimation. By analogy with squares and empirical analysis, the proportionality parameter is always taken as fixed value; that is, [40, 41]. In fact, the fixed coefficient can be replaced by a dynamic parameter, which can be estimated with regression analysis. Thus, (2) provides a simple approach to estimating the boundary dimension especially when spatial data are short for computing the fractal dimension.

2.2. Adjusting Formulae

It is easy to evaluate the boundary dimension of a fractal region such as the Koch island and urban envelope [21, 36]. However, the method based on the geometric measure relation always overestimates the fractal dimension value [22, 30]. In order to lessen the errors resulting from (1) and (2), a formula is derived as follows [23]: where represents the revised boundary dimension of a city (a textural dimension). Equation (4) suggests a linear relation between and . If we calculate the value of using (3) or log-linear regression analysis, we can estimate the value by means of (4).

In reality, a city is a complex spatial system, and urban area does not correspond to a 2-dimensional region. In this case, (1) is not enough to describe the geometric measure relation between urban area and perimeter. According to the studies of Cheng [42], Imre [43], and Imre and Bogaert [44], (1) can be generalized as follows:where denotes the fractal dimension of urban form within the urban envelope (a structural dimension). The form dimension can be estimated with the following formula [23]:which suggests a hyperbolic relation between and . From (4) and (6), it follows thatwhich suggests a hyperbolic relation between and . Accordingly, . If , then we have , and vice versa. If so, we will have a Euclidean object. Based on (5), the allometric scaling exponent expressed by (2) can be rewritten as where refers to the revised scaling exponent of the area-perimeter allometry. An allometric exponent is usually a ratio of one fractal dimension to another fractal dimension. In many cases, it is allometric scaling exponents rather than fractal dimensions that play an important role in spatial analysis of urban systems [30, 45, 46].

Using the mathematical models, allometric scaling relations, and the fractal parameter formulae, we can describe and analyze the spatial development and evolution of urban systems in the real world. For an urban system, the fractal dimensions and the related allometric scaling exponents can be classified as global fractal parameters and local fractal parameters. The global fractal parameters can be reckoned with the cross-sectional data and used to describe a system of cities as a whole, while the local fractal parameters can be figured out by using the data of individual cities and used to describe each city as an element in the urban system.

3. Empirical Analysis

3.1. Study Area and Method

The area-perimeter scaling and the adjusting formulae of fractal dimension can be applied to the actual cities by means of remote sensing data. Yangtze River Delta in China is taken as a study area to make an empirical analysis. The region includes 68 cities and towns, which can be regarded as an urban system (Figure 1). The remote sensing images used in this research came from National Aeronautics and Space Administration (NASA), including Landsat MSS, TM, and ETM images from 1985, 1996, and 2005, which were downloaded from the Earth Science Data Interface (ESDI) at the Global Land Cover Facility (GLCF) center of University of Maryland (http://glcf.umd.edu/). Firstly, these images were transformed to the GCS_WGS_1984 geographic coordinate system and Asia_Lambert_Conformal_Conic projected coordinate system. Then, the supervised classification method (SCM) was employed to extract the built-up area and the boundary of each city in a given year. The results can be manually corrected in ESRI ArcGIS to ensure better accuracy. The SCM is one of the popular remote sensing image classification techniques [47]. The basic idea of this technique is to select “training areas” as representative samples to identify the land cover classes in an image. The classifier is then used to attach labels to all the image pixels according to the trained parameters. The maximum likelihood classification (MLC) is commonly used to generate the trained parameters, such as mean vectors and variance-covariance matrix. Then, the land cover is classified based on the spectral signature that is defined by these trained parameters. In this work, we use the software of ERDAS IMAGINE which has a well-defined classifier for the SCM. We classify the Landsat images into five classes of land cover: city, forest, agriculture, water, and other. The land cover of city can be obtained and then the rough boundary of cities can be extracted on the base of the classification.

The next step is the manual error correction of the land cover of cities obtained in the first step. We check the boundary by comparing with the existing maps. The Google Map and Chinese Administrative Division Maps for various years are the major sources for manual error correction. The classified land cover of city is geocoded with these maps. So we can drop the pixels which are not located within the city administrative boundary. Afterwards, we can use a range of critical landmark features (e.g., river, lake, and mountain) as benchmarks to correct the outlier pixels and the pixels that are wrongly classified. Based on the error correction, the much more accurate city boundary and urban land use can be extracted in ArcGIS. As such, we can extract the accurate built-up area which falls into each urban envelope, and then calculate areal and perimetric values through the Calculate Geometry function in the Attribute Table of ArcGIS. Finally, we have 3 datasets of urban areas and perimeters for the 68 cities and towns in three years (Table 1).

Using the relations and formulae, we can evaluate the scaling exponents and fractal dimensions. If the observational data have been processed, the procedure of parameter estimation is as follows.

Step 1 (preliminary estimation of fractal parameters). Using the geometric measure relation and formulae including (1) and (2), we can calculate the initial boundary dimension and the proportionality coefficient by means of scatterplots and the least squares regression.

Step 2 (revision of the global fractal parameters). Using (4) and (6), we can transform the initial global boundary dimension into the revised global boundary dimension and global form dimension , which are used to describe the urban system in the study area.

Step 3 (estimation of the local fractal parameters). Using the approximate formula, (3), and the values from Step 1, we can estimate the initial local boundary dimension for each city.

Step 4 (revision of the local fractal parameters). Using (4) and (6), we can transform the initial local boundary dimension into the revised local boundary dimension and local form dimension , which are used to describe the individual cities in the study area.

Step 5 (related calculations). Fractal parameters can be associated with traditional spatial measurements. For example, the local boundary dimension is a function of the compactness of urban form. The compactness ratio of a city can be figured out with urban area and perimeter.

3.2. Results

The regression analysis based on OLS can be employed to estimate the allometric parameters and the global fractal parameters. The log-log scatterplots show that the relations between urban area and perimeter follow the allometric scaling law (Figure 2). The slopes of the trend lines give the values of the scaling exponent in different years. Taking natural logarithms on both sides of (1) yields a linear relation as below:where and . Using the least squares calculation, we can fit (9) to the datasets displayed in Table 1. The results including the values of the proportionality coefficient , scaling exponent , and the goodness of fit are tabulated as below (Table 2). By (2), (4), and (6), we can further estimate the initial boundary dimension , the revised boundary dimension , and the form dimension . The values of boundary dimension range from 1.4 to 1.5. Accordingly, the revised boundary dimension is about 1.2, and the corresponding form dimension is around 1.7.

If we had the remote sensing images with enough high resolution of all the cities, we would calculate the form and boundary dimension for each city [16, 27]. If so, the local fractal parameters could be evaluated with the regression analysis or other algorithms such as the major axis method (MAM). Unfortunately, we have only the images of the region instead of each city. In this case, the local parameters can be estimated by (3). The key is to determine the value of . There are two possible approaches to estimating the value. One is empirical method. As indicated above, a fixed coefficient is proposed by empirical analysis [40, 41]. However, if we take , the boundary dimension of the cities in our study area will be significantly overestimated. For example, for 1985, the average value of the 68 cities and towns is about 2.258, which is greater than the Euclidean dimension of the embedding space and thus unreasonable. Generally speaking, the boundary dimension comes between 1 and 1.5, and the average value is often near 1.25 [20, 30, 48]. The other is the regression method. By the least squares calculation of cross-sectional datasets, we can estimate the values using (9). This is an unfixed parameter: the value changes over time (Table 2). Using the variable value to replace the fixed value in (3), we can estimate the boundary dimension of each city. Then, using (4) and (6), we can further estimate the revised boundary dimension and the corresponding form dimension for each city. All these results represent the local fractal parameters of the urban system in the study region (Table 3).

A discovery is that the average values of the local fractal parameters are approximately equal to the corresponding global fractal parameters of the whole cities and towns. Comparing the means of the fractal dimensions listed in Table 3 with the related parameter values displayed in Table 2, we can find the consistency of the global estimation with the local average (Table 4). The global parameter values can be used to analyze the development of the system of cities and towns at the macrolevel (the whole), while the local parameter values can be employed to analyze the evolution process of the individual cities at the microlevel (the parts or elements).

3.3. Global and Local Analysis

Generally speaking, a mathematical model reflects the macrostructural properties of a system, while the model parameters reflect the microinteraction of elements. It is necessary to examine the mathematical expressions of area-perimeter relations and the corresponding fractal parameters. Analyzing the model forms and parameter values will result in a number of findings.

First, the allometric scaling degenerated from power law to linear relations from 1985 to 2005. A real allometric scaling relation takes on a double logarithmic relation, but it sometimes degenerates and changes to single logarithmic relations including exponential relation and logarithmic relation (semidegeneration), or even to a linear relation (full degeneration) [49]. In 1985, the area-perimeter relation followed the power law, and the goodness of fit () of the power function was significantly higher than those of the linear function and the single logarithmic functions. However, in 1996, the case was different, and the goodness of fit of the linear function is a little higher than that of the power function. In 2005, the goodness of fit of the linear function is significantly higher than that of the power function (Table 5). A complex system such as a city and a system of cities can be modeled by more than one mathematical form. Different mathematical models reflect different states of system evolution. Fractals suggest the optimized structure of nature. In theory, fractal models represent the optimum state of city development as a fractal object can fill its space in the best way. The degeneration of the fractal measure relation indicates some disorder problems in the processes of urbanization and urban evolution. If the power law degenerates into a linear relation, we can treat it as quasiallometric scaling in light of fractal theory in order to bring the parameter values into comparison (Table 4).

Second, both the boundary dimension and the form dimension approach certain constants. This seems to suggest some conversation law of fractal parameter evolution. The average revised boundary dimension approaches a constant: for 1985, for 1996, and for 2005. In short, the mean of the textual dimension is close to 1.2, namely, . Accordingly, the average form dimension approaches another constant: for 1985, for 1996, and for 2005. In short, the mean of the structural dimension is close to 1.7, namely, (Table 4). The average value of the form dimension lends further support to the suggestion that the structural dimension of urban space changes around 1.7 [4, 23, 48]. The is an interesting value for the form dimension of cities. By computer simulation, the fractal dimension of spatial aggregates such as cities proved to approach 1.7 averagely [50]; by empirical analysis, the form dimension mean of cities is near 1.7 [4]. It has been demonstrated that the reasonable range of the form dimension comes between 1.5 and 2 [48]. If , the urban field is not efficiently utilized due to underfilling of space; on the contrary, if , the urban field is excessively utilized owing to overfilling of space. In theory, the best value of the form dimension is close to 1.7 [23, 35]. In terms of (8), the corresponding boundary dimension is about .

Third, the fractal models and the fractal parameters reflect the social and economic state of China. For a long time, China is socialist country based on centrally planned economy rather than market economy. Chinese city development is always associated with its political and economic conditions. Since the introduction of the economic reforms and opening-up policies at the end of 1978 and with the gradual establishment of a socialist market economic system from 1992, namely, after Deng’s South Tour Speeches, the command economy based on top-down mechanism and the market economy based on bottom-up mechanism combined with one another to form a mixed economy (Table 6). China has witnessed a rapid urban development [51]. If we relate the urban evolution with the political and economic background of China, we can understand the changes of fractals and fractal dimension of Chinese cities. At the beginning of reform and opening up of China, that is, in 1985, the fractal dimension values are in disorder to some extent, which can be reflected by the relation between the boundary dimension and compactness ratio (see Figure 3 and Section 4). Moreover, several values are abnormal: the value is less than 1 or even less than 0, and the value is greater than 2. The values range from −0.330 (an abnormal value) to 1.525, and values vary from 0.398 to 2.053. However, in 2005, both the and values fluctuate between 1 and 2. The values range from 1.046 to 1.422, and the values vary from 1.542 to 1.916 (Table 3). The boundary dimension of cities is proved to be a function of the compactness of urban form. The fractal dimension evolution from chaos to order can be mirrored by the relations between the boundary dimension and the compactness ratio (Figure 3).

4. Questions and Discussion

An urban analytical process based on the area-perimeter scaling is illustrated by using the 68 cities and towns in the study area. Through the empirical analysis, the following questions have been clarified. First, the fractal dimension of urban shape includes form dimension and boundary dimension, which are associated with one another. However, the boundary dimension of cities and patches of a city are indeed overestimated by using the traditional area-perimeter scaling formula. In Table 3, many values are close to 2, and this is unreasonable because the dimension of a fractal boundary is often less than 3/2 [48]. By the ideas from multifractals and spatial correlation, we can demonstrate that is a threshold value of urban shape dimension, and there is a reasonable numerical range for the form dimension; that is, [48]. However, in previous literature, the boundary dimension is often overestimated using (1). Thus, the boundary dimension used to be expressed as . Here represents the Euclidean dimension of urban region. This suggests that the urban area is regarded as a 2-dimensional measure. However, the average dimension of urban area in a system is not actually equal to 2 [12, 23, 42, 44]. According to (5), the revised boundary dimension can be expressed as , where the form dimension is near 1.7 where average is concerned [4, 30]. Second, using the adjusting formula of fractal dimension, we can not only revise the boundary dimension indicating urban texture but also estimate the form dimension indicative of urban structure. The two formulae were previously derived [23], but the systematic empirical analysis is made for the first time. Third, in terms of the area-perimeter scaling, the fractal dimensions of urban shape fall into two types: the global parameters indicative of a system of cities and the local parameter indicating the individual cities in the urban system. If we estimate the global parameters by using the least squares calculations, the proportionality coefficient values in the local fractal parameter formula can be obtained by the constants of the global models. In other words, if we estimate the values for an urban system by the regression analysis based on (1), we can compute the boundary dimension of each city in the urban system by substituting the value in (3). This implies that we can use a dynamics proportionality factor instead of a fixed coefficient to estimate the local fractal parameters of cities. The analytical process of fractal cities can be demonstrated by a block diagram (Figure 4).

The main limitation of this study lies in the resolution of the remote sensing images. The resolution of the remote sensing images of the cities and towns is not high enough to guarantee the spatial data quality of individual cities. In this case, we cannot abstract large datasets of urban area and perimeter based on different scales for each human settlement. For a city, we have only two data points: urban area and perimeter. As a result, we cannot calculate the boundary dimension and form dimension for individual cities directly. In other words, we cannot evaluate the local fractal parameters for individual cities by the least squares calculation due to absence of datasets for each city. What we can do is to evaluate the global fractal parameters for the urban systems by means of the area-perimeter scaling based on cross-sectional data and then estimate the local fractal parameters for each city by the approximate formula. The shortcomings of the empirical analysis are as below. First, there exist many outliers in the results. For example, the boundary dimension of Deqing city is negative, which is absurd as a fractal dimension must be greater than 0 and less than the Euclidean dimension of the embedding space. A fractal city based on the remote sensing images is actually defined in a 2-dimensional space [4]. Thus, the Euclidean dimension of the embedding space for the fractal city is . The boundary dimension of Shangyu is greater than 2, going beyond the Euclidean dimension of its embedding space. Second, the comparability of the fractal dimension values in different year is doubtful. Generally speaking, owing to space filling, the form dimension increases year by year, until the limit is reached (). However, what with the data quality and the urban sprawl, the form dimension values of the cities and towns randomly fluctuated. Third, the numerical relation between the boundary and the compactness ratio degenerated. The boundary dimension of a city can be associated with the compactness of its urban form. The mathematical relation between the compactness ratio and the boundary dimension can be derived as below [20]:where Co denotes the compactness ratio of urban form, represents the boundary dimension ( or ), is the length of urban perimeter, and is a coefficient. However, based on the datasets in Table 3, the exponential form of (10) is replaced by a hyperbolic relation such aswhere and represent the intercept and slope if the reciprocal of the fractal dimension, and is treated as an independent variable. Despite these shortcomings, we can make use of the advantages and bypass the disadvantages of the remote sensing data in virtue of the formulae of the fractal dimension estimation.

A popular misunderstanding of the boundary dimension should be pointed out here. In previous literature, the initial boundary dimension is always associated with the stability index of urban evolution [52]. If the parameter , the urban boundary is regarded as resulting from Brownian motion and thus the urban development is considered to be unstable. This viewpoint is wrong. The reasons are as follows. First, it is the self-affine record dimension rather than the self-similar trail dimension that can be directly related with Brownian motion [28]. However, the boundary dimension of urban shape is a self-similar trail dimension instead of a self-affine record dimension [35]. Second, as indicated above, the boundary dimension is always numerically overestimated by the old formula. In fact, the properly estimated boundary dimension is seldom greater than or equal to 1.5. An urban boundary can be treated as a fractal line, and the dimension of a fractal line is usually less than 1.5. The boundary dimension should be revised using the adjusting formula (6). The revised boundary dimension values often come between 1 and 1.5 [23]. Finally, a random process does not indicate an unstable process. If the self-affine record dimension is close to 1.5, it suggests randomicity rather than instability of urban growth.

5. Conclusions

In this paper, a case study is employed to illustrate a systematic analytical process of fractal cities using the area-perimeter allometric scaling. In particular, we try to develop an approach to make the best use of the advantages and bypass the disadvantages of spatial data by means of the new fractal dimension formulae. The innovation of this paper is with two aspects: First, the global fractal parameters and the local fractal parameters are integrated to make urban spatial analyses, and the numerical link between the global and local parameters is revealed by statistical average. Second, the form dimension and the boundary dimension are associated with one another to describe urban shapes. Third, the area-perimeter scaling is employed to estimate the proportionality coefficient for the approximate formula of boundary dimension. An interesting finding is the numerical link between the global and local parameters.

Based on the empirical analysis using the methodology developed in this work, the main conclusions can be reached as follows. First, the relationships between urban area and perimeter follow the allometric scaling laws; the scaling exponent is the ratio of the urban boundary dimension to the urban form dimension. However, the boundary dimension value used to be overestimated because the form dimension () is mistaken for Euclidean dimension of urban area (). By using the revision formulae of fractal dimension, we can correct the results of fractal dimension estimation. What is more, the form dimension can be estimated through the boundary dimension. Thus, the form dimension indicative of urban structure and the boundary dimension indicating urban texture can be combined with each other to characterize urban evolution. Second, a fractal system of cities should be characterized by both the global and the local scaling parameters. The fractal parameters of an urban system, including form dimension and boundary dimension, fall into two groups: global parameters and local parameters. The global parameters reflect the spatial properties of a system of cities, while the local parameters reflect the geographical feature of individual cities. The average values of the local parameters are approximately equal to the corresponding values of the global parameters. This suggests that the global parameters can be decomposed into local parameters. By means of the average values, we can associate the global parameters with local parameters and further connect the macrolevel of an urban system with the microlevel of elements in the urban system. Third, the fractal measure relations between urban area and urban envelope may degenerate because of disorder competition of cities. Fractals follow scaling laws, which can be expressed as a set of power functions. However, a power-law relation between urban area and perimeter sometimes degenerates into a semilogarithmic relation or even a linear relation. Fractal structure based on power laws is a kind of optimized structure of natural and human systems. A fractal object can fill its space in the most efficient way. Fractal relation degeneration suggests some latent problems of urban evolution. If a city or urban system does not follow the allometric scaling law, the fractal structure may be broken and the structure of the city or system of cities should be improved by rational city planning.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research was supported financially by the National Natural Science Foundation of China (Grant no. 41171129). The support is gratefully acknowledged.