Next Article in Journal
Ensemble Estimation of Information Divergence
Next Article in Special Issue
Intrinsic Computation of a Monod-Wyman-Changeux Molecule
Previous Article in Journal
Quantum Dynamics and Non-Local Effects Behind Ion Transition States during Permeation in Membrane Channel Proteins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Index Based on Binary Entropy to Confirm the Spatial Expansion Degree of Urban Sprawl

School of Geographic and Oceanographic Sciences, Nanjing University, Nanjing 210046, China
*
Authors to whom correspondence should be addressed.
Entropy 2018, 20(8), 559; https://doi.org/10.3390/e20080559
Submission received: 18 May 2018 / Revised: 25 June 2018 / Accepted: 25 July 2018 / Published: 27 July 2018
(This article belongs to the Special Issue Information Theory in Complex Systems)

Abstract

:
The phenomenon of urban sprawl has received much attention. Accurately confirming the spatial expansion degree of urban sprawl (SEDUS) is a prerequisite to controlling urban sprawl. However, there is no reliable metric to accurately measure SEDUS. In this paper, based on binary entropy, we propose a new index named the spatial expansion degree index (SEDI), to overcome this difficulty. The study shows that the new index can accurately determine SEDUS and, compared with other commonly used measures, the new index has an obvious advantage in measuring SEDUS. The new index belongs to the second-order metrics of point pattern analysis, and greatly extends the concept of entropy. The new index can also be applied to other spatial differentiation research from a broader perspective. Although the new index is influenced by the scaling problem, because of small differences between different scales, given that the partition scheme in the research process is the same, the new index is a quite robust method for measuring SEDUS.

1. Introduction

Urban sprawl is considered a threat to sustainable economic and social development in the world today [1]. It is reported that worldwide the built-up land area increased 58,000 km2 from 1970 to 2000 and, at this rate, by 2030, 1,527,000 km2 of built-up land will be added globally [2]. Meanwhile, the built-up land expansion often encroaches on suburban peripheries and results in cropland fragmentation [3,4], so its spatial form mostly appears ‘‘dispersed’’ [5]. The phenomenon of the high-speed and dispersed expansion of built-up land has become known as urban sprawl [5,6,7,8].
There are two different viewpoints about the impact of urban sprawl. One idea is that urban sprawl can relieve traffic congestion in urban central area, relax residential pressure, increase employment, and promote social fairness [5,9]. The opposite view is that development associated with urban sprawl not only decreases the amount of cultivated land and forest resources, endangers food and energy security, destroys the environment, and also affects biological diversity [1,10,11,12,13,14,15]. As time has passed, the phenomenon of urban sprawl has received much attention, and most scholars share the view that the negative aspects of the impact outweigh the positive ones [8].
Urban sprawl is considered a multidimensional phenomenon [5,8,16,17,18,19]. As the subject became popularized, some of the researchers began conceptualizing multifaceted analytical frameworks and establishing multiple indices by GIS analysis or descriptive statistical analysis to measure the sprawl [6,18,19,20]. Galster et al. conducted a pioneering study, which proposed eight dimensions—density, continuity, concentration, clustering, centrality, nuclearity, mixed uses, and proximity—to measure housing sprawl in 13 urbanized areas [18]. Tsai employed a set of four dimensions—size, density, degree of equal distribution and degree of clustering—to systematically identify a metropolitan form [20]; Hamidi and Ewing operationalized four dimensions—development density, land use mix, activity centering and street accessibility—to study the changes of urban sprawl between 2000 and 2010 in the United States [19]. However, if only considering urban sprawl as a typical geographical phenomenon and focusing solely on the external spatial form, two indispensable dimensions increasing amount of built-up land area and spatial expansion degree should be confirmed [7,21]. Measurement of increasing amount of built-up land area can be easily obtained by remote sensing imageries and techniques [8,13,22,23,24,25,26,27,28,29,30], however, there is no reliable metrics to measure the spatial expansion degree.
Accurately confirming the spatial expansion degree of urban sprawl (SEDUS) will not only provide monitoring information of land use/cover change and urban planning decision-making bases, but also make contributions to controlling urban sprawl [11,29,31,32] and the protection of the environment and ecosystems [7,29,30,31,32,33]. There are about 50 metrics that characterize the spatial degree of expansion of urban sprawl, and many of them have been implemented in the public domain statistical package FRAGSTATS3, which are generally known as spatial metrics or landscape metrics [6,20,34]. Spatial metrics have found important applications in quantifying SEDUS [6,16,20,23,31], however, many of them are correlated [6,34]. Often, different metrics employed may also result in contradictory conclusions [35].
Most of the spatial metrics used in the research of spatial heterogeneity require partitioning the study area, so partitioning methods can affect the values of the spatial metrics employed [14,29,34,36,37], which is the so-called modifiable area unit problem [29,36,37]. The problem includes two aspects: scale effect and zoning effect, respectively. The former refers to the phenomenon that the analytical results change with the aggregation of spatial data through different basic area units; the latter refers to the change caused by the effect of different aggregation methods on the same basic unit. Since scale effect to the spatial metrics varies from different metrics while the zoning effect has the commonness [36,37], the research in the influence of partition method on a new metric is often just focused on the scale effect [7,21,29,37].
Entropy methods, such as Shannon’s entropy and Renyi’s entropy [38], are perhaps the most widely used spatial metrics for measuring SEDUS with the integration of remote sensing and GIS [20,29,39,40,41,42,43]. In fact, the term of entropy is originally a thermodynamic item [44,45,46]. Claude Shannon introduced the concept of entropy into information theory and physics and provided a new scientific field [47]. Jaynes proposed the principle of maximum entropy on Shannon’s work to apply to subjective statistic inference; since then, entropy has occurred as a statistical inference method for supporting spatial locations [48,49]. Theil obtained a new understanding of Shannon’s entropy in the form of its relative value (i.e., relative entropy) and regarded it as a measure of a discrete probability distribution [50]. Yeh and Li used relative Shannon entropy (RS in the paper) to successfully measure the extent of urban sprawl. Henceforth, RS has become a main tool in measuring SEDUS [6,29,44]. However, for RS, there is still a considerable deficiency: built-up land with different configurations may have the same entropy value (Figure 1). That is, although RS is more robust than the other measures, developing new tools and metrics are highly demanded. Like many other spatial metrics, the modifiable area unit problem also exists in entropy methods [29].
In this paper, a new index named spatial expansion degree index (SEDI) is proposed to confirm the spatial expansion degree of urban sprawl. The structure of this paper is as follows: first, based on choosing the spatial analysis method and determining the dimensions of SEDUS, we elaborate on the building-up process of the new index and introduce the study area. Second, based on the case result, we discuss the obvious advantage of the new index compared with other commonly used measures, outline the characteristics of the new index, investigate the influence of partition method on the index, and propose its relativity in measuring SEDUS and future research directions. Finally, we provide our conclusions.

2. Materials and Methods

2.1. Building Spatial Expansion Degree Index (SEDI)

2.1.1. Choosing the Spatial Analysis Method and Determining Dimensions of SEDUS

In the building-up process of the new index, we adopted the second-order metrics of point pattern analysis as our spatial analysis method [51]. The so-called second-order metrics are also known as distance-based point pattern measures, which take an object of interest at a particular location in the study region as a point pattern and use the distances between points to characterize the geospatial object of interest [51]. That is, the spatial analysis method is on the basis of point-pair distance [51,52,53].
In this paper, by reference to partition methodology of the Geographical Analysis Machine (GAM) [52], we laid out a two-dimensional grid over the study region as a partitioning method, treated each grid centroid as the center of the subzone, and attached the message of the built-up land area in the corresponding subzone to each point (grid centroid). In this way, we set up a one-to-one correspondence relationship between the built-up land in the subzone and each grid centroid in the point pattern, that is, the term “two linked subzones” is equivalent to a “point-pair” in meaning. In addition, for the sake of convenient research, the built-up process of the new index goes from one point-pair to multiple point-pairs step-by-step.
To confirm the spatial expanding degree of urban sprawl (SEDUS), it is necessary to determine the dimensions of SEDUS. This paper is influenced by the ideal of Jaeger et al. work which focused only on the external spatial form of urban sprawl, regardless of the driving forces of urban sprawl [7,21]. As seen in Figure 2, the SEDUS changes with two dimensions: the distances between point-pairs and built-up land area configurations in subzones; intuitively, the greater the distances and the more balanced area configurations, the greater the SEDUS. That is, SEDUS is a function of the two dimensions.

2.1.2. Defining the Point-Pair Entropy Distance

The task consists of two sides: how to measure the two dimensions of SEDUS and how to integrate the two dimensions in only one point-pair configuration. First, we utilized Pythagoras’ theorem (the straight-line distance between two points) [51] to measure the magnitude of the distance. Second, we employed binary entropy function (recorded as B ( i , i ) ) to measure the area configuration. B ( i , i ) is the shorthand form of binary entropy function, when 0 log 2 0 is taken, the value of B ( i , i ) is equal to 0; when an even area configuration occurs, B ( i , i ) attains its maximum value [54]. These specific properties of B ( i , i ) match the process of the SEDUS well (Figure 3). Thirdly, we used the product (recorded as E D ( i , i ) ) of B ( i , i ) and the distances to characterize the comprehensive effect of SEDUS. The related mathematical expressions are written as:
E D ( i , i ) = k B ( i , i ) r ( i , i ) .
where:
  B ( i , i ) = P i log 2 P i P i log 2 P i .
where:
P i = p i / ( p i + p i ) ; P i = 1 P i = p i / ( p i + p i ) .
E D ( i , i ) is called point-pair entropy distance, which represents the SEDUS in only one point-pair configuration; k is a constant and is usually set to 1; r ( i , i ) represents the straight-line distance between two points; p i is the proportion of built-up land area in the ith subzone over the built-up land area in whole area ( p i = x i / 1 n x i ) [6,14,29]. From Equation (1), if the value of r ( i , i ) is fixed, the magnitude of E D ( i , i ) varies with the binary entropy function; when the expansion with a fully-even area configuration ( P i = P i ) occurs, the SEDUS reaches the highest value; given that the area configuration is fixed ( r ( i , i ) is variable), the greater the value of r ( i , i ) , the greater the magnitude of E D ( i , i ) . Therefore, E D ( i , i ) well reflects the law that the SEDUS changes with its dimensions.

2.1.3. Giving Weights to Point-Pair Entropy Distance and Optimizing r ( i , i )

Based on the point-pair entropy distance, we discuss the SEDUS with its dimensions in more than two point-pair configurations. First, as can be seen from Figure 4, if two point-pairs in a configuration occur, binary entropy function is incapable of describing this change. We also found in Figure 4 that, other things being equal, the greater the value of the sum of the proportion of the built-up land area, the higher the spatial expansion degree. To express this effect, we took Y ( i , i ) (recorded as Y ( i , i ) = p i + p i ) as the weight of B ( i , i ) , and named the adjusted E D ( i , i ) as the weighted point-pair entropy distance. In fact, Y ( i , i ) is always equal to constant 1 in the only one point-pair configuration, thus, E D ( i , i ) can be regarded as a special case. The mathematical expression of the weighted point-pair entropy distance is written as:
W E D ( i , i ) = k B ( i , i ) Y ( i , i ) r ( i , i ) .
Second, let us investigate SEDUS with more than two points distances. In the case, the so-called “problem of monotonous reaction to increased spreading of three urban patches” will emerge [21]. This means that, when an urban area is broken up into more parts and the distribution of these parts becomes more dispersed, the SEDUS will increase monotonously as Figure 1 demonstrates.
If computing SEDUS in more than a two-point configuration with no three being collinear, we can accomplish the task well by using the method of summing the resulting value of W E D ( i , i ) . However, when configurations with three points in a line occur, the method is powerless, as Figure 5 shows.
Jaeger et al. proposed a power function f ( x ) = c x γ ( 0 < γ < 1 ; c is a constant) to solve the problem (see [7]). We select the simplest square root function f ( x ) = 2 x to optimize r ( i , i ) . Thus, Equation (4) is adjusted as:
W E D ( i , i ) = k B ( i , i ) Y ( i , i ) 2 r ( i , i ) .
where W E D ( i , i ) denotes the optimized weighted point-pair entropy distance; and 2 r ( i , i ) denotes the optimized r ( i , i ) . Simple verification shows that Equation (5) is effective.

2.1.4. Constructing Spatial Expansion Degree Index (SEDI)

Suppose a zone (the whole study area) with n subzones, if two subzones ( i , i ) are randomly selected, n ( n 1 ) / 2 point-pairs in total can be obtained. We apply Equation (5) to all point-pairs and average all resulting values, the global average optimized weighted point-pair entropy distance can be developed as:
D n = 2 k ( B ( i , i ) Y ( i , i ) 2 r ( i , i ) ) / [ n ( n 1 ) ] .
Since Equation (6) contains an absolute amount, 2 r ( i , i ) , it is difficult to compare SEDUS in different regions. To solve the problem, by introducing the relative value, we further adjust Equation (6). The first work is to calculate the maximum of D n . Since the partition scheme in the process of research remains the same, the point-to-point distances must be fixed. Meanwhile, when equal probability occurs within the zone (if a fully-even area configuration occurs, then p i = 1 / n ) (see [44]), B ( i , i ) reaches the maximum 1 (See Equations (2) and (3)), at the same time, Y ( i , i ) obtains the maximum 2 / n , and D n obtains the maximum. The resulting maximum is written as:
M A X ( D n ) = k ( 2 / n ) D n = 4 k 2 r ( i , i ) / [ n 2 ( n 1 ) ] .
Then, SEDI is defined as the form of the relative value:
S E D I = D n / M A X ( D n ) = n W ( i , i ) 2 r ( i , i ) / [ 2 2 r ( i , i ) ] .
where:
W ( i , i ) = B ( i , i ) Y ( i , i ) = B ( i , i ) ( p i + p i ) .
W ( i , i ) is a weighted version of B ( i , i ) (see Equations (2) and (3)), which is only associated with p i . The values of the index range from 0 to 1, and the greater the value is, the greater the SEDUS.

2.2. Study Case

For providing a clear visual effect, instead of choosing the real data of urban sprawl, we selected a grid map containing six study subplots as a study case to verify the robustness of the new index. Meanwhile, to demonstrate the effect that SEDUS changes with the two dimensions mentioned above, we presented two contrastive patterns, where Pattern 1 denotes the increasing SEDUS with increasing dispersion of built-up land in constant area configurations (Figure 6); Pattern 2 shows an increasing SEDUS in more balanced built-up land area configurations (Figure 7).

3. Results

The ArcGIS 10.1 software from ESRI (Environmental Systems Research Institute, Inc. Redlands, California, USA) was employed to perform the data processing, whose flow is shown in Figure 8. In order to improve the computational efficiency and accuracy, we made a Python program (see Supplementary Materials) to perform the data processing. The resulting values of SEDI applied to Pattern 1 (Figure 6) were: (a) (0.1484), (b) (0.1759), (c) (0.1820), (d) (0.2099), (e) (0.2348), and (f) (0.2571), respectively; while the values in Pattern 2 (Figure 7) were: (a) (0.0000), (b) (0.0783), (c) (0.1296), (d) (0.1830), (e) (0.2481), and (f) (0.3076), respectively. The results matched well the tendency of urban sprawl in Pattern 1 and 2. Pattern 1 can correspond to many practical applications, for example, controlling urban sprawl requires optimizing the structure and layout of land use in the process of drawing up or amending an overall plan of land use, especially in the case of the total amount of built-up land area being constant. That is, based on resulting values of SEDI applied to different area configurations of the land use layout, a suitable plan for urban development can be chosen. Pattern 2 reproduces the phenomenon that a new town (i.e., satellite town) is emerging far away from the old central city. The new index SEDI can confirm SEDUS and provide accurate changing information of SEDUS for urban planning decision-making. Based on the information, targeted strategies of controlling urban sprawl can be proposed.

4. Discussion

4.1. Comparsion with Commonly Used Indices of SEDUS

We selected the patch shape index (PSI) from the landscape metrics [55], the Geary coefficient (GI) [20], the relative Shannon entropy index (RS) [29], and the leapfrog index (LPI) [31] to perform their capability to compare with the new index, because they are the most widely used in measuring SEDUS [6,20,29,31,55,56]. Their mathematical expressions and meanings are listed in Table 1. Applying these indices to Figure 6 and Figure 7, we recorded the resulting values in Table 2 and Table 3, respectively. As in Table 2 and Table 3, the values of PSI and GI show fluctuation trends, which do not match SEDUS in Figure 6 and Figure 7, so, we can conclude that the two indices are not suitable to characterize SEDUS. At the same time, we can see that RS and LPI can confirm the gradual increasing SEDUS in Figure 7, while the two indices cannot identify the trends in Figure 6. Since Pattern 2 (Figure 7) represents a commonplace urban sprawl pattern, which misleads us with the illusion that both indices, especially RS, can confirm SEDUS of any urban sprawl patterns. As in Table 2 and Table 3, the new index, SEDI, could confirm the two patterns, so, the new index has an apparent advantage in measuring SEDUS.

4.2. The Extreme Value and Certainty of SEDI

Dietzel et al. empirically tested a theory of spatiotemporal urban growth dynamics, which suggests that the patterns of SEDUS can be characterized into two types: diffusion and compaction, respectively [57]. This theory also suggests that the processes is continuously observable until the geographical area becomes completely urbanized. Here, we can regard diffusion pattern as the maximum of spatial expansion degree, and compaction pattern as the minimum. If the maximum and the minimum are mapped as integer 1 and integer 0, the changes in the value of SEDI can be understood as the dynamic change between the two states of diffusion and compaction. Each different value of SEDI can indicate a certain spatial expansion degree. If there is no extreme value, there is no standard to quantitatively describe spatial expansion degree, so the extreme value of SEDI is the qualitative prerequisite for the certainty of SEDI.

4.3. The New Index Deepening the Understanding of the Concept of Second-Order Metrics

The new measure SEDI belongs to the second-order metrics of point pattern analysis. The second-order metrics, such as F(d) function, G(d) function, and Ripley’s K function [51], characterize the general relationship of homogeneous point-pairs (homogeneous point-pairs can be understood as each point in a point-pair with a same size). The new index can examine the relationship of heterogeneous point-pairs, which greatly deepens the understanding of the concept of second-order metrics. In other words, conventional point-pairs can be taken as a special case of point-pair entropy in which each point has an equal area weight.

4.4. Effect of Partitioning Method on SEDI

Like most spatial metrics, the new index is also influenced by the scaling problem. Here, we demonstrate the scale effect by comparing the values of SEDI in different scenarios at three scales (Figure 9).
From Table 4, we can see that, given a fixed extent, as for the same scenario and different scales, the smaller the scale is, the greater the value of SEDI. However, as can be seen from Table 4, although the values of SEDI varied with different scales, the relative differences are small.

4.5. Relativity of SEDI in Measuring Urban Sprawl

As the literature states, if only considering urban sprawl as a typical geographical phenomenon, measuring urban sprawl requires two dimensions: the increasing amount of built-up land area and spatial expansion degree [7,21]. Meanwhile, from the study described above, measuring the spatial expansion degree also requires two dimensions the distances between point-pairs and area configurations, respectively. The new index can confirm well the spatial expansion degree, but it is only one side in measuring urban sprawl. The reason is that the index of SEDI based on “p” is a relative numerical indicator and “increasing amount of built-up land area” is an absolute one. The same thing applies to the index of RS.

4.6. Future Work

In the next step, the new index can be used individually or in combination with other metrics to analyze a region in different periods, or different regions in the same period. From a broader perspective, the new index cannot just limit to urban sprawl research, but can be extended to other spatial differentiation research, such as arable land distribution and the spread of infectious diseases.

5. Conclusions

Urban sprawl can result in excessive resource consumption, environmental pollution, and ecosystem destruction, which has received much attention and has become an important research topic. Accurately confirming the spatial expansion degree of urban sprawl (SEDUS) will not only provide monitoring information of land use/cover change for government officials and planners to propose targeted strategies of controlling urban sprawl, but also make contributions to the protection of the environment and ecosystems, thus, it has important theoretical and practical significance. Until now, about 50 metrics used to measure SEDUS have been proposed, however, there is no reliable metric among them. In this paper, and through previous work, we proposed developing a new index named the spatial expansion degree index (SEDI), to overcome the difficulty.
Based on the literature, this paper articulated that, regardless of driving forces of urban sprawl, SEDUS changes with two dimensions: the distances between linked subzones and built-up land area configurations, respectively, and confirmed that SEDUS is a function of the two dimensions. We utilized Pythagoras’ theorem to measure the magnitude of the different distances, employed binary entropy function to measure the other’s dimension, and used the product of the two dimensions to characterize the comprehensive effect of SEDUS. Meanwhile, the study adopted the second-order metrics of point pattern analysis as spatial analysis method, followed the research approach from a single point-pair to multiple point-pairs and, finally, the new spatial expansion degree index (SEDI) was used to accurately confirm that SEDUS was successfully built up.
This study shows that the new index is capable of identifying SEDUS in constant or changing built-up land area configurations. Meanwhile, we also found that the patch shape index (PSI) and the Geary coefficient (GI) are not suitable to characterize SEDUS, and the relative Shannon entropy index (RS) and leapfrog index (LPI) could identify the SEDUS in the changing area configurations, but they could not identify SEDUS in the constant area configurations. Therefore, the new index has an obvious advantage in measuring SEDUS.
The new index is a new addition to existing second-order metrics for point pattern analysis, and it has three excellent characteristics: extreme value, certainty, and point-pair heterogeneity. The extreme value is the key basis of quantitative SEDUS. The values of the index range from the minimum 0 to the maximum 1, and the greater the value is, the greater the SEDUS. Thus, if the extreme values are mapped as integer 1 and integer 0, the changes in the value of SEDI can be understood as the dynamic change between the two extreme states, that is, SEDUS has quantitative uncertainty. The point-pair heterogeneity of the new index, which is mainly reflected in binary entropy function, greatly deepening the understanding of the concept of the point-pair from homogeneity to heterogeneity.
In the next step, the new index can be applied to practical applications, such as measuring SEDUS or amending an overall plan of land use. From a broader perspective, the new index can be applied to other different spatial differentiation studies, such as arable land distribution and the spread of infectious diseases. Due to the index of SEDI being a relative numerical indicator, it has relativity in measuring urban sprawl. Meanwhile, this study shows that the new index is still influenced by partition schemes. However, given the same partition scheme in the research process, the new index is considered as a robust tool in measuring SEDUS for government officials and relevant theoretical research workers.

Supplementary Materials

The following are available online at https://www.mdpi.com/1099-4300/20/8/559/s1.

Author Contributions

Z.C. had the original idea for the study and was responsible for the writing, drawing and programming. Y.Z. was responsible for funding acquisition and dividing the work. X.J. was responsible for revising and polishing the manuscript. All authors have read and approved the final manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pielke, R.A. Land use and climate change. Science 2005, 310, 1625–1626. [Google Scholar] [CrossRef] [PubMed]
  2. Seto, K.C.; Michail, F.; Burak, G.; Reilly, M.K. A meta-analysis of global urban land expansion. PLoS ONE 2011, 6, e23777. [Google Scholar] [CrossRef] [PubMed]
  3. Visalatchi, A.; Padmanaban, R. Land Use and Land Cover Mapping and Shore Line Changes Studies in Tuticorin Coastal Area Using Remote Sensing. Int. J. Remote Sens. 2012, 1, 1–12. [Google Scholar]
  4. Monishiya, B.G.; Padmanaban, R. Mapping and change detection analysis of marine resources in Tuicorin and Vembar group of Islands using remote sensing. Int. J. Adv. For. Sci. Manag. 2012, 1, 1–16. [Google Scholar]
  5. Ewing, R. Is Los Angeles-style sprawl desirable? J. Am. Plan. Assoc. 1997, 1, 107–126. [Google Scholar] [CrossRef]
  6. Bhatta, B.; Saraswati, S.; Bandyopadhyay, D. Urban sprawl measurement from remote sensing data remote sensing data. Appl. Geogr. 2010, 30, 731–740. [Google Scholar] [CrossRef]
  7. Jaeger, J.A.; Bertiller, R.; Schwick, C.; Cavens, D.; Kienast, F. Urban permeation of landscapes and sprawl per capita: New measures of urban sprawl. Ecol. Indic. 2010, 10, 427–441. [Google Scholar] [CrossRef]
  8. Wilson, E.H.; Hurd, J.D.; Civco, D.L.; Prisloe, M.P.; Arnold, C. Development of a geospatial model to quantify, describe and map urban growth. Remote Sens. Environ. 2003, 86, 275–285. [Google Scholar] [CrossRef]
  9. Peter, G.; Harry, W.R. Are compact cities a desirable planning goal? J. Am. Plan. Assoc. 1997, 63, 95–106. [Google Scholar]
  10. Miller, M.D. The impacts of Atlanta’s urban sprawl on forest cover and fragmentation. Appl. Geogr. 2012, 34, 171–179. [Google Scholar] [CrossRef]
  11. Tan, M.; Li, X.; Xie, H.; Lu, C. Urban land expansion and arable land loss in China-a case study of Beijing-Tianjin-Hebei region. Land Use Policy 2005, 22, 187–196. [Google Scholar] [CrossRef]
  12. Cabral, P.; Feger, C.; Levrel, H.; Chambolle, M.; Basque, D. Assessing the impact of land-cover changes on ecosystem services: A first step toward integrative planning in Bordeaux, France. Ecosyst. Serv. 2016, 22, 318–327. [Google Scholar] [CrossRef]
  13. Padmanaban, R.; Bhowmik, A.K.; Cabral, P.; Zamyatin, A.; Almegdadi, O.; Wang, S. Modelling Urban Sprawl Using Remotely Sensed Data: A Case Study of Chennai City, Tamilnadu. Entropy 2017, 19, 164. [Google Scholar] [CrossRef]
  14. Sun, H.; Forsythe, W.; Waters, N. Modeling Urban Land Use Change and Urban Sprawl: Calgary, AB, Canada. Netw. Spat. Econ. 2007, 7, 353–376. [Google Scholar] [CrossRef]
  15. Stow, D.A.; Hope, A.; McGuire, D.; Verbyla, D.; Gamon, J.; Huemmrich, F.; Houston, S.; Racine, C.; Sturm, M.; Tape, K.; et al. Remote sensing of vegetation and land-cover change in Arctic Tundra Ecosystems. Remote Sens. Environ. 2004, 89, 281–308. [Google Scholar] [CrossRef] [Green Version]
  16. Arribas-Bel, D.; Nijkamp, P.; Scholten, H. Multidimensional urban sprawl in Europe: A self-organizing map approach. Comput. Env. Syst. 2011, 35, 263–275. [Google Scholar] [CrossRef] [Green Version]
  17. Sudhira, H.S.; Ramachandra, T.V.; Jagdish, K.S. Urban sprawl: Metrics, dynamics and modeling using GIS. Int. J. Appl. Earth. Obs. Geoinform. 2004, 5, 29–39. [Google Scholar] [CrossRef]
  18. Galster, G.; Hanson, R.; Ratcliffe, M.R.; Wolman, H.; Coleman, S.; Freihage, J. Wrestling sprawl to the ground: Defining and measuring an elusive concept. Hous. Policy Debate 2001, 12, 681–717. [Google Scholar] [CrossRef]
  19. Hamidi, S.; Ewing, R. A longitudinal study of changes in urban sprawl between 2000 and 2010 in the united states. Landsc. Urban Plan. 2014, 128, 72–82. [Google Scholar] [CrossRef]
  20. Tsai, Y. Quantifying urban form: Compactness versus ‘sprawl’. Urban Stud. 2005, 42, 141–161. [Google Scholar] [CrossRef]
  21. Jaeger, J.A.; Bertiller, R.; Schwick, C.; Kienast, F. Suitability criteria for measures of urban sprawl. Ecol. Indic. 2010, 10, 397–406. [Google Scholar] [CrossRef]
  22. Hasse, J.E.; Lathrop, R.G. Land resource impact indicators of urban sprawl. Appl. Geogr. 2003, 23, 159–175. [Google Scholar] [CrossRef]
  23. Herold, M.; Couclelis, H.; Clarke, K.C. The role of spatial metrics in the analysis and modeling of urban change. Comput. Env. Syst. 2005, 29, 339–369. [Google Scholar]
  24. Padmanaban, R. Integrating of urban growth modelling and utility management system using Spatiotemporal data mining. Int. J. Adv. Earth Sci. Eng. 2012, 1, 13–15. [Google Scholar]
  25. Jat, M.K.; Garg, P.K.; Khare, D. Monitoring and modelling of urban sprawl using remote sensing and GIS techniques. Int. J. Appl. Earth Obs. Geoinform. 2008, 10, 26–43. [Google Scholar] [CrossRef]
  26. Aithal, B.H.; Ramachandra, T.V. Visualization of Urban Growth Pattern in Chennai Using Geoinformatics and Spatial Metrics. J. Indian Soc. Remote Sens. 2016, 44, 617–633. [Google Scholar] [CrossRef]
  27. Jayaprakash, M.; Senthil Kumar, R.; Giridharan, L.; Sujitha, S.B.; Sarkar, S.K.; Jonathan, M.P. Bioaccumulation of metals in fish species from water and sediments in macrotidal Ennore creek, Chennai, SE coast of India: A metropolitan city effect. Ecotoxicol. Environ. Saf. 2015, 120, 243–255. [Google Scholar] [CrossRef] [PubMed]
  28. Gowri, V.S.; Ramachandran, S.; Ramesh, R.; Pramiladevi, I.R.; Krishnaveni, K. Application of GIS in the study of mass transport of pollutants by Adyar and Cooum Rivers in Chennai, Tamilnadu. Environ. Monit. Assess. 2008, 138, 41–49. [Google Scholar] [CrossRef] [PubMed]
  29. Yeh, A.G.O.; Li, X. Measurement and monitoring of urban sprawl in a rapidly growing region using entropy. Photogramm. Eng. Remote Sens. 2001, 67, 83–90. [Google Scholar]
  30. Ferdinent, J.; Padmanaban, R. Development of a Methodology to Estimate Biomass from Tree Height Using Airborne Digital Image. Int. J. Adv. Remote Sens. GIS 2013, 2, 49–58. [Google Scholar]
  31. Frankel, A.; Ashkenazi, M. Measuring urban sprawl: How can we deal with it? Environ. Plan. B 2008, 35, 56–79. [Google Scholar] [CrossRef]
  32. Frenkel, A.; Orenstein, D.E. Can urban growth management work in an era of political and economic change? J. Am. Plan. Assoc. 2012, 78, 16–33. [Google Scholar] [CrossRef]
  33. Dewan, A.M.; Humayun, K.; Nahar, M.K.; Rahman, Z. Urbanisation and environmental degradation in Dhaka Metropolitan Area of Bangladesh. Environ. Sustain. Dev. 2012, 11, 118–147. [Google Scholar] [CrossRef]
  34. Riitters, K.H.; O’Neill, R.V.; Hunsaker, C.T.; Wickham, J.D.; Yankee, D.H.; Timmins, S.P.; Jones, K.B.; Jackson, B.L. A factor analysis of landscape pattern and structure metrics. Landsc. Ecol. 1995, 10, 23–39. [Google Scholar] [CrossRef]
  35. Herold, M.; Goldstein, N.C.; Clarke, K.C. The spatiotemporal form of urban growth: Measurement, analysis and modeling. Remote Sens. Environ. 2003, 86, 286–302. [Google Scholar] [CrossRef]
  36. Openshaw, S. Optimal zoning systems for spatial interaction models. Environ. Plan. B 1977, 9, 169–184. [Google Scholar] [CrossRef]
  37. Jelinski, D.E.; Wu, J. The modifiable areal unit problem and implications for landscape ecology. Landsc. Ecol. 1996, 11, 129–140. [Google Scholar] [CrossRef]
  38. Fan, Y.; Yu, G.; He, Z.; Yu, H.; Bai, R.; Yang, L.; Wu, D. Entropies of the Chinese Land Use/Cover Change from 1990 to 2010 at a County Level. Entropy 2017, 19, 51. [Google Scholar] [CrossRef]
  39. Torrens, P.M. A toolkit for measuring sprawl. Appl. Spat. Anal. Policy 2008, 1, 5–36. [Google Scholar] [CrossRef]
  40. Tewolde, M.G.; Cabral, P. Urban sprawl analysis and modeling in Asmara, Eritrea. Remote Sens. 2011, 3, 2148–2165. [Google Scholar] [CrossRef]
  41. Yue, W.; Liu, Y.; Fan, P. Measuring urban sprawl and its drivers in large Chinese cities: The case of Hangzhou. Land Use Policy 2013, 31, 358–370. [Google Scholar] [CrossRef]
  42. Sabet Sarvestani, M.; Ibrahim, A.L.; Kanaroglou, P. Three decades of urban growth in the city of Shiraz, Iran: A remote sensing and geographic information systems application. Cities 2011, 28, 320–329. [Google Scholar] [CrossRef]
  43. Dewan, A.M.; Corner, R.J. Spatiotemporal Analysis of Urban Growth, Sprawl and Structure. In Dhaka Megacity: Geospatial Perspectives on Urbanization, Environment and Health; Dewan, A., Corner, R.J., Eds.; Springer Science and Business Media: Dordretch, The Netherlands, 2014; pp. 99–121. [Google Scholar]
  44. Cabral, P.; Augusto, G.; Tewolde, M.; Araya, Y. Entropy in urban systems. Entropy 2013, 15, 5223–5236. [Google Scholar] [CrossRef]
  45. Clausius, R. The Mechanical Theory of Heat: With its Applications to the Steam-Engine and to the Physical Properties of Bodies; J. Van Voorst: London, UK, 1867; p. 376. [Google Scholar]
  46. Uffink, J. Bluff your way in the second law of thermodynamics. Stud. Hist. Philos. Sci. Part B Stud. Hist. Philos. Mod. Phys. 2001, 32, 305–394. [Google Scholar] [CrossRef]
  47. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  48. Wilson, A. Entropy in urban and regional modelling: Retrospect and prospect. Geogr. Anal. 2010, 42, 364–394. [Google Scholar] [CrossRef]
  49. Leopold, L.B.; Langbein, W.B. The Concept of Entropy in Landscape Evolution; United States Geological Survey: Washington, DC, USA, 1962.
  50. Thomas, R.W. Information Statistics in Geography; Headey Brothers Ltd., The Invicta Press: Ashford, UK, 1981. [Google Scholar]
  51. O’Sullivan, D.; Unwin, D.J. Point Pattern Analysis. In Geographic Information Analysis, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NY, USA, 2010; pp. 121–156. ISBN 9780470549094. [Google Scholar]
  52. O’Sullivan, D.; Unwin, D.J. Practical Point Pattern Analysis. In Geographic Information Analysis, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NY, USA, 2010; pp. 157–186. ISBN 9780470288573. [Google Scholar]
  53. Perry, G.L.W.; Miller, B.P.; Enright, N.J. A comparison of methods for the statistical analysis of spatial point patterns in plant ecology. Plant. Ecol. 2006, 187, 59–82. [Google Scholar] [CrossRef]
  54. Wikipedia. Available online: https://en.wikipedia.org/wiki/Binary_entropy_function (accessed on 17 May 2018).
  55. Wu, J.G. Landscape Ecology-Pattern, Process, Scale and Hierarchy, 2nd ed.; Higher Education Press: Beijing, China, 2007; pp. 107–116. ISBN 978–7-04-020879-5. (In Chinese) [Google Scholar]
  56. Jiang, F.; Liu, S.; Yuan, H.; Zhang, Q. Measuring urban sprawl in Beijing with geo-spatial indices. J. Geoginform. Sci. 2007, 17, 469–478. [Google Scholar] [CrossRef]
  57. Dietzel, C.; Oguz, H.; Hemphill, J.J.; Clarke, K.C.; Gazulis, N. Diffusion and coalescence of the Houston metropolitan area: Evidence supporting a new urban theory. Environ. Plan. B 2005, 32, 231–236. [Google Scholar] [CrossRef]
Figure 1. The limitations of Relative Shannon Entropy (RS) in describing the spatial expanding degree of urban sprawl (SEDUS). Note: This graph contains three scenarios, and each scenario contains 16 grid subzones or cells, where the blue grid cell area (1 × 1) indicates that the built-up land area is equal to the total area of the reporting subzone (the size of a cell is set to 1); every yellow subzone indicates the land without buildings in it. Obviously, the SEDUS of the three subplots from left to right is gradually deepening. As for RS, its value ranges from 0 to 1, and the greater the value is, the greater the SEDUS [6,29]. From the resulting values (the values are 0.5000), each zone has the same expansion degree.
Figure 1. The limitations of Relative Shannon Entropy (RS) in describing the spatial expanding degree of urban sprawl (SEDUS). Note: This graph contains three scenarios, and each scenario contains 16 grid subzones or cells, where the blue grid cell area (1 × 1) indicates that the built-up land area is equal to the total area of the reporting subzone (the size of a cell is set to 1); every yellow subzone indicates the land without buildings in it. Obviously, the SEDUS of the three subplots from left to right is gradually deepening. As for RS, its value ranges from 0 to 1, and the greater the value is, the greater the SEDUS [6,29]. From the resulting values (the values are 0.5000), each zone has the same expansion degree.
Entropy 20 00559 g001
Figure 2. Spatial expanding degree of urban sprawl (SEDUS) is a function of the distances between point-pairs and built-up land area configurations in subzones. Note: Each blue grid cell can be understood as the basic unit of built-up land area, and each dotted green square box denotes the basic partition (i.e., subzone). (ac) denote two linked subzones with constant built-up land area configurations and different distances, which characterizes the process that the SEDUS is gradually deepening with increasing distances. (df) represent two linked subzones with different built-up land area configurations and the same distance, which reveals that the SEDUS grows greater in more balanced area configurations (the process could be viewed as a course that the area unit of built-up land was moved from one subzone to another).
Figure 2. Spatial expanding degree of urban sprawl (SEDUS) is a function of the distances between point-pairs and built-up land area configurations in subzones. Note: Each blue grid cell can be understood as the basic unit of built-up land area, and each dotted green square box denotes the basic partition (i.e., subzone). (ac) denote two linked subzones with constant built-up land area configurations and different distances, which characterizes the process that the SEDUS is gradually deepening with increasing distances. (df) represent two linked subzones with different built-up land area configurations and the same distance, which reveals that the SEDUS grows greater in more balanced area configurations (the process could be viewed as a course that the area unit of built-up land was moved from one subzone to another).
Entropy 20 00559 g002
Figure 3. The diagram in the description of the correspondences between the binary entropy function and SEDUS in one point-pair configuration. Note: The legends in the figure are same as in Figure 2. Scenarios (a), (b), (c), (d), and (e), respectively, represent different SEDUS. Scenarios (a) and (e) display the lowest SEDUS, while Scenarios (c) displays the highest. Scenarios (a), (c) and (e) respectively correspond to the three special points of binary entropy function: the starting point, the midpoint, and the terminal point, respectively. Scenario (b) represents a random point in the changing process of the expansion from the lowest degree to the highest, and the process just corresponds to the increasing zone of the binary entropy function, while the scenario (d) displays just the opposite.
Figure 3. The diagram in the description of the correspondences between the binary entropy function and SEDUS in one point-pair configuration. Note: The legends in the figure are same as in Figure 2. Scenarios (a), (b), (c), (d), and (e), respectively, represent different SEDUS. Scenarios (a) and (e) display the lowest SEDUS, while Scenarios (c) displays the highest. Scenarios (a), (c) and (e) respectively correspond to the three special points of binary entropy function: the starting point, the midpoint, and the terminal point, respectively. Scenario (b) represents a random point in the changing process of the expansion from the lowest degree to the highest, and the process just corresponds to the increasing zone of the binary entropy function, while the scenario (d) displays just the opposite.
Entropy 20 00559 g003
Figure 4. An example for illustrating the necessity of giving weights to point-pair entropy distances. Note: The legends in the figure are same as in Figure 2. Due to two scenarios with different built-up land area configurations ( p 1 = p 2 = 2 p 3 = 2 p 4 ) and the same distances ( r 12 = r 34 ), SEDUS in Scenario 1 is twice of that in Scenario 2. That is, giving no weight to point-pair entropy distance (other things being equal), when applying Equation (1) to the two scenarios, the resulting values of E D ( i , i ) are the same because P i = P i (see Equation (3)). However, when we operated the adjusted E D ( i , i ) (Equation (4)), we found that the resulting values aligned with the SEDUS in the two scenarios.
Figure 4. An example for illustrating the necessity of giving weights to point-pair entropy distances. Note: The legends in the figure are same as in Figure 2. Due to two scenarios with different built-up land area configurations ( p 1 = p 2 = 2 p 3 = 2 p 4 ) and the same distances ( r 12 = r 34 ), SEDUS in Scenario 1 is twice of that in Scenario 2. That is, giving no weight to point-pair entropy distance (other things being equal), when applying Equation (1) to the two scenarios, the resulting values of E D ( i , i ) are the same because P i = P i (see Equation (3)). However, when we operated the adjusted E D ( i , i ) (Equation (4)), we found that the resulting values aligned with the SEDUS in the two scenarios.
Entropy 20 00559 g004
Figure 5. An example for illustrating the necessity of optimizing r ( i , i ) in configurations with three collinear subzones. The legends in the figure are same as in Figure 2. As the figure shows, three subzones with equal built-up land area line up, where the two outer subzones are fixed and the third one can freely shift along the line (Scenario 1); when the third subzone is right in the middle (Scenario 2), the SEDUS is the highest; in other words, the closer the third subzone shifts to one of the two outer subzones, the lower the SEDUS will be. However, when applying Equation (4) to the two scenarios, due to p 1 = p 2 = p 3 ( p i is the proportion of built-up land area and does not refer to the patch’s code in the subzone) and r 13 + r 12 + r 23 = 2 d , their resulting values are the same. Therefore, it is necessary to optimize r ( i , i ) .
Figure 5. An example for illustrating the necessity of optimizing r ( i , i ) in configurations with three collinear subzones. The legends in the figure are same as in Figure 2. As the figure shows, three subzones with equal built-up land area line up, where the two outer subzones are fixed and the third one can freely shift along the line (Scenario 1); when the third subzone is right in the middle (Scenario 2), the SEDUS is the highest; in other words, the closer the third subzone shifts to one of the two outer subzones, the lower the SEDUS will be. However, when applying Equation (4) to the two scenarios, due to p 1 = p 2 = p 3 ( p i is the proportion of built-up land area and does not refer to the patch’s code in the subzone) and r 13 + r 12 + r 23 = 2 d , their resulting values are the same. Therefore, it is necessary to optimize r ( i , i ) .
Entropy 20 00559 g005
Figure 6. Increasing SEDUS with increasing dispersion of built-up land in constant area configurations (Pattern 1). Note: Each subplot contains 16 grid subzones or cells, where a blue grid cell (1 m × 1 m) indicates that the built-up land area in the cell is equal to the total area of the cell; each yellow cell indicates the land without buildings in it. The black spot in each grid cell indicates the geographic center of each subzone; the values inside parentheses are the relative coordinate values of its geographic center.
Figure 6. Increasing SEDUS with increasing dispersion of built-up land in constant area configurations (Pattern 1). Note: Each subplot contains 16 grid subzones or cells, where a blue grid cell (1 m × 1 m) indicates that the built-up land area in the cell is equal to the total area of the cell; each yellow cell indicates the land without buildings in it. The black spot in each grid cell indicates the geographic center of each subzone; the values inside parentheses are the relative coordinate values of its geographic center.
Entropy 20 00559 g006
Figure 7. Increasing SEDUS in more balanced built-up land area configurations (Pattern 2). Note: The red grid cell can be viewed as old urban area, the blue one can be regarded as the patch of newly-grown built-up land, and (af) can be understood as the process of urban sprawl over time; other legends or labels are the same as in Figure 6.
Figure 7. Increasing SEDUS in more balanced built-up land area configurations (Pattern 2). Note: The red grid cell can be viewed as old urban area, the blue one can be regarded as the patch of newly-grown built-up land, and (af) can be understood as the process of urban sprawl over time; other legends or labels are the same as in Figure 6.
Entropy 20 00559 g007
Figure 8. The plot of the data processing flow. Note: “P_INPUT” (or “P_NEAR”) corresponds to p or p* in Equations (2) and (3).
Figure 8. The plot of the data processing flow. Note: “P_INPUT” (or “P_NEAR”) corresponds to p or p* in Equations (2) and (3).
Entropy 20 00559 g008
Figure 9. SEDUS in different scenarios at different scales. Note: From the vertical configurations in the figure, the SEDUS from (a) to (f) is gradually deepening; the horizontal subplots indicate the same configurations at three scales. The legends or labels in the figure are same as in Figure 6.
Figure 9. SEDUS in different scenarios at different scales. Note: From the vertical configurations in the figure, the SEDUS from (a) to (f) is gradually deepening; the horizontal subplots indicate the same configurations at three scales. The legends or labels in the figure are same as in Figure 6.
Entropy 20 00559 g009
Table 1. The mathematical expressions and meanings of the commonly used indices of SEDUS.
Table 1. The mathematical expressions and meanings of the commonly used indices of SEDUS.
IndexMathematical ExpressionExplanation
PSI P S I = P / 4 A P denotes the perimeter of the polygon measured, and A is its area. The greater the value, the greater the SEDUS [26].
GI G I = ( N 1 ) 1 N 1 N W i j ( X i X j ) 2 2 ( 1 N 1 N W i j ) 1 N ( X i X ) 2 i, j refer to different subzones, respectively; X i , X j refer to the area proportions of built-up lands of subzones i, j over total area, respectively; X is the mean value of all subzones; W i j is the reciprocal of the distance between the centroids of subzones i, j. The lower the value of GI is, the higher the SEDUS [16].
RS R S = 1 n p i log 2 ( p i ) / log 2 n p i is the proportion of different land use in grid cells, and n is the total number of subzones. Values of the index range from 0 to 1. The greater the value, the greater the SEDUS [3].
LPI L P I = A o u t / A u A o u t denotes “leapfrog area”, A u is the total area of built-up land. The higher the value, the greater the SEDUS [20].
Table 2. Capability comparison of the new index and other commonly used indices in constant built-up land area configurations.
Table 2. Capability comparison of the new index and other commonly used indices in constant built-up land area configurations.
IndicesThe Resulting Values of the Different Scenarios in Figure 6
Subplot (a)Subplot (b)Subplot (c)Subplot (d)Subplot (e)Subplot (f)
PSI1.00001.25001.50002.00002.00002.0000
GI1.07700.82031.06971.04930.97680.9108
RS0.50000.50000.50000.50000.50000.5000
LPI1.00001.00001.00001.00001.00001.0000
SEDI0.14840.17590.18200.20990.23480.2571
Table 3. Capability comparison of the new index and other commonly used indices in changing built-up land area configurations.
Table 3. Capability comparison of the new index and other commonly used indices in changing built-up land area configurations.
IndicesThe Resulting Values of the Different Scenarios in Figure 7
Subplot (a)Subplot (b)Subplot (c)Subplot (d)Subplot (e)Subplot (f)
PSI1.00001.41421.44341.50001.56521.4289
GI0.83561.03641.06691.07681.05011.0036
RS0.00000.25000.39620.50000.58050.6462
LPI0.00000.50000.66670.75000.80000.8333
SEDI0.00000.07830.12960.18300.24810.3076
Table 4. The values of SEDI in different scenarios at different scales.
Table 4. The values of SEDI in different scenarios at different scales.
ScaleSubplot (a)Subplot (b)Subplot (c)Subplot (d)Subplot (e)Subplot (f)
1 m × 1 m0.14840.17590.18200.20990.23480.2571
0. 5 m × 0.5 m0.17140.19760.20320.22960.25360.2750
0.25 m × 0.25 m0.17570.20160.20700.23320.25700.2783

Share and Cite

MDPI and ACS Style

Chen, Z.; Zhou, Y.; Jin, X. A Novel Index Based on Binary Entropy to Confirm the Spatial Expansion Degree of Urban Sprawl. Entropy 2018, 20, 559. https://doi.org/10.3390/e20080559

AMA Style

Chen Z, Zhou Y, Jin X. A Novel Index Based on Binary Entropy to Confirm the Spatial Expansion Degree of Urban Sprawl. Entropy. 2018; 20(8):559. https://doi.org/10.3390/e20080559

Chicago/Turabian Style

Chen, Zhen, Yinkang Zhou, and Xiaobin Jin. 2018. "A Novel Index Based on Binary Entropy to Confirm the Spatial Expansion Degree of Urban Sprawl" Entropy 20, no. 8: 559. https://doi.org/10.3390/e20080559

APA Style

Chen, Z., Zhou, Y., & Jin, X. (2018). A Novel Index Based on Binary Entropy to Confirm the Spatial Expansion Degree of Urban Sprawl. Entropy, 20(8), 559. https://doi.org/10.3390/e20080559

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop