Analysis of Urban Growth Pattern Using Gis, Kolkata
Analysis of Urban Growth Pattern Using Gis, Kolkata
To cite this article: B. Bhatta (2009): Analysis of urban growth pattern using remote sensing and
GIS: a case study of Kolkata, India, International Journal of Remote Sensing, 30:18, 4733-4746
This article may be used for research, teaching, and private study purposes. Any
substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing,
systematic supply, or distribution in any form to anyone is expressly forbidden.
The publisher does not give any warranty express or implied or make any representation
that the contents will be complete or accurate or up to date. The accuracy of any
instructions, formulae, and drug doses should be independently verified with primary
sources. The publisher shall not be liable for any loss, actions, claims, proceedings,
demand, or costs or damages whatsoever or howsoever caused arising directly or
indirectly in connection with or arising out of the use of this material.
International Journal of Remote Sensing
Vol. 30, No. 18, 20 September 2009, 4733–4746
Analysis of urban growth pattern using remote sensing and GIS: a case
study of Kolkata, India
B. BHATTA*
Senior Systems Engineer, Computer Aided Design Centre, Department of Computer
Science & Engineering, Jadavpur University, Kolkata — 700 032, India
Urban growth pattern has a direct influence on urban development processes both
at the city and neighbourhood levels, but on the other hand growth pattern is also
affected by land and housing policy, current built-up area, available developable
area, change in household count, population, economy, and others. Remote sen-
sing and GIS techniques can be used to analyse and model the urban growth
pattern. This study was carried out to identify and analyse the urban growth
pattern of Kolkata, India. Several statistical methods were applied to recognize
and analyse this pattern, which shows increasing dispersed development of the city
with the declining population growth rate. Apart from the derived results, this
study also shows the potentiality of remote sensing and GIS techniques for the
analysis of urban growth patterns.
1. Introduction
Rapid urbanization in the world is quite alarming, especially in developing countries
such as India (Kumar et al. 2007). Urbanization is a process whereby productive
agricultural lands, forests, surface water bodies and groundwater prospects are being
irretrievably lost (Pathan et al. 1989, 1991). This is mainly due to uncontrolled
population growth resulting in serious problems with regards to scarcity of food,
informal settlements, environmental pollutions, destruction of ecological structure,
unemployment, and so on (Maktav and Erber 2005).
Recent generation models of urban dynamics have dealt with diverse themes.
According to Batty (2000), there are currently some 20 or more applications of models
to cities, such as for the diffusion or migration of resident populations (Portugali et al.
1997), competitive location of economic activities (Benati 1997), joint expansion of
urban surfaces and traffic networks (Batty and Xie 1997), generic urban growth
(Clarke et al. 1997) and urban land-use dynamics (Deadman et al. 1993, Batty and
Xie 1994, Phipps and Langlois 1997, White and Engelen 1997, White et al. 1998),
amongst others.
Urban sprawl has become one of the most important issues facing most cities at
the onset of the new millennium. Housing developments and increases in imper-
vious surfaces encroach on lands that were formerly agricultural lands, grasslands,
wetlands, water bodies, and forests. Changes to the landscape occur every day
with significant implications for taxation, quality of life, water quality, agricul-
tural viability, wildlife habitat and social equity. Impervious surfaces are
*Email: basu_bhatta@rediffmail.com
International Journal of Remote Sensing
ISSN 0143-1161 print/ISSN 1366-5901 online # 2009 Taylor & Francis
http://www.tandf.co.uk/journals
DOI: 10.1080/01431160802651967
4734 B. Bhatta
2. Study area
Kolkata (formerly Calcutta), the capital of the Indian state of West Bengal, was
considered as the study area. It is the main business, commercial and financial hub
of eastern India and the north-eastern states. The city is more than 300 years old and it
served as the capital of India during the British Governance until 1911. One of the
largest metropolitan regions of the world, Kolkata serves as the global example of the
potentials and problems facing mega-cities in developing nations (Banerjee 2005).
Kolkata is located at 22 33¢ N 88 20¢ E on the east bank of River Hooghly (Ganges
Delta) at an elevation of approximate 9 m (Shuttle Radar Topography Mission,
National Aeronautics and Space Administration, February 2000). The city is located
154 km (96 miles) upstream from the Bay of Bengal.
Kolkata city, under the jurisdiction of the Kolkata Municipal Corporation (KMC),
now has an area of 185 km2 subdivided into 141 wards. The Kolkata urban agglom-
eration, however, is continuously expanding: the urban agglomeration (Kolkata
Metropolitan Area) is spread over 1750 km2. However, this study is limited to the
KMC area only.
The area is spread linearly along the banks of the River Hooghly in a north–south
direction. The east–west dimension of the proper city stretches from the Hooghly
River in the west to roughly the Eastern Metropolitan Bypass in the east.
l Landsat Thematic Mapper (TM) image (path 138, row 44) from 14 November
1990
l Landsat Enhanced Thematic Mapper Plus (ETM+) image (path 138, row 44)
from 17 November 2000
l Indian Remote Sensing (IRS) Resourcesat 1 Linear Imaging Self-scanning
Sensor (LISS) III image (path 108, row 56) from 5 February 2005.
Spectral details of the aforementioned imageries are given in table 1. It is worth
mentioning that the thermal bands of Landsat TM and ETM+ were not considered
for the analysis.
The vector map of Kolkata (KMC area) was used for clipping and subsetting the
satellite images. Other ancillary data such as topographical maps, land-use/cover
maps, municipal ward maps, and census data were also used either as references or
for analytical purposes.
Downloaded by [Duke University Libraries] at 10:53 11 July 2012
The images were obtained as standard products, i.e. geometrically and radiome-
trically corrected. However, due to the several standards and references used by the
image supplying agencies, the overlay of the images did not match with considerable
accuracy. To solve this problem, images were co-registered so that the overlay
matched with the sub-pixel accuracy. Nearest-neighbour resampling was used to
transform the images retaining original pixel values.
Images from different sensors have differences in their spatial resolution. One
approach to encounter this problem is to resample the higher resolution images so
that they match the resolution of the lowest spatial resolution image. But resampling
of an image to change the pixel size either averages the neighbouring pixel values (in
the case of bi-linear or bi-cubic) or pixel dropout/duplication occurs (in the case of
nearest neighbour). For this reason, images were kept without changing the pixel size,
despite the varying accuracy level of classification with the different spatial, spectral
and radiometric resolutions. This variation was accepted not to compromise with the
spatial detail and pixel value.
In the next step, co-registered images were clipped by the vector map of the KMC
area. The clipped images were then classified using a nonparametric parallelepiped
classifier to extract the built-up area (along with other impervious surfaces). It is
worth mentioning that if the study is concerned only with urban growth, that is only
the classes of urban and non-urban use are under consideration, detailed land-use
4 0.5–0.6 1 0.45–0.52
5 0.6–0.7 2 0.52–0.60 2 0.52–0.59
6 0.7–0.8 3 0.63–0.69 3 0.62–0.68
7 0.8–1.1 4 0.76–0.90 4 0.77–0.86
5 1.55–1.75 5 1.55–1.70
7 2.08–2.35
4736 B. Bhatta
maps are no longer necessary and a simple binary classification of remotely sensed
data is enough.
The classified images were then compared with the topographical maps, available
land-use/cover maps for the respective time period and with current ground checks to
assess the accuracy. An overall accuracy between 71% and 84% was achieved and
agreed to proceed with.
Kolkata is highly heterogeneous in nature. The north–south expansion is
roughly divided into North, Central and South Kolkata. The North Kolkata
locality is the oldest part of the city. South Kolkata grew mostly after indepen-
dence and consists of elite localities. Central Kolkata houses the central business
district (CBD) area. Much of the city (mainly in the eastern part) was originally
a vast wetland, reclaimed over the decades to accommodate the city’s burgeon-
ing population. West side contains a very large dock area, industries, residence,
and wetlands as well. Due to the heterogeneous nature of the city, in this study,
Downloaded by [Duke University Libraries] at 10:53 11 July 2012
it has been subdivided into five different subparts considering the distribution of
different land-use/cover classes as shown in figure 1. Table 2 summarizes the
characteristics of these zones and their associated ward numbers. It is important
to note that the conventional zoning concept used by local people or authorities
has not been considered; rather, land-use/cover was emphasized for demarcating
the zones.
Built-up areas for each zone and for each temporal image were calculated by
clipping the classified images with the respective vector zone, and then multiplying
the number of pixels in each zone by the pixel size.
The dates of remote sensing data were different from the dates of the decennial
census dates (2001, 1991, 1981, . . .). However, the image of November 1990 can
be associated with the census data of 1991, and the image of November 2000 can
be associated with the census data of year 2001. But, the 1975 and 2005 images
very much deviated from the census dates. Therefore, it was necessary to inter-
polate the census data, so that the population figures corresponded to the dates of
the images.
City population growth was assumed to be exponential, and the annual growth rate
was calculated using the following equation (Angel et al. 2005):
where Pt is the population of census year t, Pt+10 is the population of census year t+10,
and g is the annual population growth rate between census years t and t+10.
The population of 1975 (P1975) was interpolated using the following formula:
where P1971 is the population of 1971 = 3 714 700, and g is the annual population
growth rate between 1971 and 1981 0.0102445990.
The population of 2005 (P2005) was interpolated using the following formula:
where P2001 is the population of 2001 = 4 572 876, and g is the annual population
growth rate between 1991 and 2001 0.0042000523.
Analysis of urban growth pattern 4737
Downloaded by [Duke University Libraries] at 10:53 11 July 2012
Total
Part name Wards Characteristics area (km2)
Northern zone(N) 1–6, 13, 14, 29–32, 34, 35 Developed, mainly residential, 16
concrete dominated
mix in nature
Central zone(C) 7–12, 15–28, 36–56, 59–78, Fully developed, business 56
81–93, 105 cum residential, mainly
concrete in nature
Western zone(W) 79, 80, 128–141 Industrial and residential 39
mix, large dock area,
some wards are developing,
mix in nature
Southern zone(S) 94–104, 111–127 Newly developed 46
(or developing), mainly
residential, mix in nature
Eastern zone(E) 33, 57, 58, 106–110 Mainly wetland (mixed with 28
agriculture and trees),
developing
Figure 2. Observed census data (1971, 1981, 1991 and 2001) and interpolated data (1975 and
2005) of population.
The resulted population data and built-up/non built-up area data were then used
for analysis.
Zone
Year N C W S E KMC
N, Northern zone; C, Central zone; W, Western zone; S, Southern zone; E, Eastern zone; KMC,
Kolkata Municipal Corporation.
Zone
Year N C W S E
N, Northern zone; C, Central zone; W, Western zone; S, Southern zone; E, Eastern zone.
the non-built-up areas in these zones as well. Therefore, it is better to relate the
proportion of population and proportion of built-up area for each zone instead of
regarding the population as being a function of the built-up area.
The proportion of population and proportion of built-up area were calculated
Downloaded by [Duke University Libraries] at 10:53 11 July 2012
by dividing the population and built-up area of the respective zone by the total
population and total area of KMC, respectively. If we subtract the proportion of
population from the proportion of built-up area, the result will vary within the
range of –1 to 1 as shown in table 4, where 0 (or near to 0) indicates average
conditions. Higher negative values indicate population crowding, which may cause
serious environmental problems, traffic congestion, and a minimum of social (and
other) facilities. Higher positive values indicate higher per capita consumption of
built-up area.
The per capita land consumption refers to utilization of all the land develop-
ment initiatives such as the commercial, industrial, educational, and recreational
establishments along with the residential establishments per person. Since most of
the initiatives pave the way for the creation of jobs and subsequently help with
earning livelihoods, the development of land is seen as a direct consequence of
this and hence one can conclude that the per capita land consumption is inclusive
of all the associated land development (Sudhira et al. 2004). Higher per capita
consumption of built-up area indicates a better environment and extended
facilities.
increase in the number of working persons will not influence the built-up area.
Therefore, we need to account the percentage of area on which new construction
can be made (which is a portion of total non-built-up area). Due to lack of
information, in this study, the percentage of non-built-up area has been consid-
ered to model the relationship between economy, developable area and built-up
growth. The following simple model shows the trend of built-up change in relation
to the change in worker percentage and percentage of available area to be built-up
for a specific zone:
Wi · Bi
Ei ¼ (4)
Pn
Wi
i¼1
where Ei is the change in the percentage of built-up area between 1991 and 2001 for the
ith zone, Wi is the change in the percentage of working persons between 1991 and 2001
for the ith zone, and Bi is the percentage of available area for the ith zone in 1991 on
which new construction can be made. This can be calculated by the percentage of
non built-up area minus the percentage of area on which construction is not
possible or permissible. n is the total number of zones.
The above model shows the trend of change in the percentage of built-up areas
(Ei) between 1991 and 2001 for each zone as shown in table 6. The result supports
the reality as shown in table 5. Some variation occurs due to the consideration of
entire non built-up area as developable area, which is not true. However, this
model requires tests on other cities also. It has a major limitation that it assumes
Wi can not be 0 or less than 0. The zones should not be considered where Wi is 0
or less than 0.
Zone
N C W S E
N, Northern zone; C, Central zone; W, Western zone; S, Southern zone; E, Eastern zone.
Analysis of urban growth pattern 4743
Zone
N C W S E
N, Northern zone; C, Central zone; W, Western zone; S, Southern zone; E, Eastern zone.
Shannon’s entropy was computed to detect and quantify the urban sprawl phenom-
enon. Shannon’s entropy (Hn) has been calculated by:
X
n
Hn ¼ Pi loge ðPi Þ (5)
i¼1
where Pi is the proportion of the variable in the ith zone, i.e. the proportion of built-
up area (as a percentage) in each zone, calculated by: built-up area percentage in ith
zone/total of built-up area percentage of all zones. n = total number of zones = 5.
The value of entropy ranges from 0 to loge(n) (i.e. 1.609). A value of 0 indicates that
the distribution of built-up areas is very compact, while values closer to loge(n) reveal
that the distribution of built-up areas is dispersed. Higher values of entropy indicate
the occurrence of sprawl.
Table 7 shows that the entropy values are very much higher than the halfway mark
of loge(n) (i.e. 0.805) and it is increasing with time. It clearly shows that the city is
becoming more dispersed despite the declining population growth rate.
5. Conclusions
The discussion in the preceding sections has directed attention to the magnitude and
pattern of change in urban growth of Kolkata for the last few decades, which should
be very helpful in terms of guiding future planning and policy-making for the city.
The analysis shows a clear indication of urban sprawl or dispersion in growth for
the city of Kolkata; it can be concluded that cities in developing countries do not
always become more compact with the declining rate of population growth.
The key issue facing decision makers at local, national and international levels is
not whether or not urban expansion will take place, but rather what is likely to be the
scale of urban expansion and what needs to be done now to adequately prepare for it.
The models introduced and discussed in the preceding sections can be utilized to
provide some guidance as to the scale of change likely to result if trends of the past
However, the merits of restricting urban expansion and encouraging infill and
intensification of existing urban areas—even in the cities in industrialized coun-
tries—are by no means clear, nor is it self-evident that these are desired by the
majority of urban residents. Even in the best of circumstances, compact city policies
may have a marginal effect on the overall level of urban land consumption (Angel
et al. 2005).
The case for densification and intensification in the cities of developing countries
(like Kolkata)—where densities are, on average, three times higher than densities in
the cities of industrialized and developed countries—is even less clear (Angel et al.
2005). Therefore the merits of urban densification postulated for Kolkata or the like
seem far less convincing, where most of the city area is overcrowded and lacks
sufficient urban facilities, where property rights are not strictly enforced, where the
city has grown almost without prior city planning, where the city lacks future town-
planning schemes and even lacks current land-use/cover maps, where enforcement is
intermittent and often corrupt, and where a large part of the citizenry cannot afford
minimum standard shelter.
Needless to say, it is more expensive to provide trunk urban infrastructure in built-up
areas—especially in areas developed by the informal sector (as is the case for Kolkata)—
than to provide such services, or at least to protect the right-of-way needed for such
services before building takes place. While there are many reasons for neglecting to study
the past and present of urban growth pattern to prepare for the inevitable future growth
of cities, the absence of even minimal preparation for urban expansion—on both the
activist and the regulatory fronts—is, no doubt, an inefficient, inequitable and unsus-
tainable practice, imposing great economic and environmental costs on societies that can
ill afford them. This paper aims to present several models using remote sensing/GIS
techniques to study the past and understand the present. Profiting from this knowledge,
the future planning can be done accordingly. Our understanding of the complexity of
urban growth pattern essentially provides us with some of the tools necessary to meet it
in an efficient, equitable and sustainable manner in the years to come.
In summary, this study opened with an observation about the important role of
analytic models of urban growth pattern, proceeded to use some standard classical
and some new models to derive testable hypotheses about factors that influence urban
extent and growth. It focuses on the scope of research/application for the city planners
and other stakeholders in developing nations. The theory and models of urban spatial
growth pattern that are supported by preceding findings should prove useful for
Analysis of urban growth pattern 4745
devising policy responses to the problems associated with preparing for urban expan-
sion. However, there still remain some important policy issues for which analytical
support has to be made available. Future analytical research, using the global dataset,
should shed some light on these issues in the near future.
Acknowledgment
The author wishes to acknowledge the Director, Computer Aided Design Centre,
Jadavpur University, Kolkata for providing necessary facilities to execute this study.
The author is also grateful to Earth Science Data Interface, Global Land Cover
Facility, for providing free Landsat data.
References
ACIOLY, C.C. and DAVIDSON, F., 1996, Density in urban development. Building Issues, 8,
Downloaded by [Duke University Libraries] at 10:53 11 July 2012
pp. 3–25.
ALMEIDA, C.M., MONTEIRO, A.M.V., MARA, G., SOARES-FILHO, B.S., CERQUEIRA, G.C.,
PENNACHIN, C.S.L. and BATTY, M., 2005, GIS and remote sensing as tools for the
simulation of urban land-use change. International Journal of Remote Sensing, 26,
pp. 759–774.
ANGEL, S., SHEPPARD, S.C. and CIVCO, D.L., 2005, The Dynamics of Global Urban Expansion.
Transport and Urban Development Department. (Washington DC: The World Bank),
pp. 1–200. Available online at: http://www.citiesalliance.org/doc/resources/upgrading/
urban-expansion/worldbankreportsept2005.pdf (last accessed 13 February 2008).
ARNOLD, C.L. JR and GIBBONS, J.C., 1996, Impervious surface coverage—the emergence of a
key environmental indicator. Journal of the American Planning Association, 62,
pp. 243–258.
BANERJEE, A., 2005, Population growth, environment and development: some issues in
sustainability of the mega city of Kolkata (Calcutta), West Bengal. In Proceedings
of the National Seminar on Population Environment and Nexus, Population
Environment Centre, IIPS, Deonar, Mumbai, 21 October 2005. Available online at:
http://www.iipsenvis.nic.in/paper/fp_anuradhab.pdf (last accessed 13 February 2008).
BARNES, K.B., MORGAN III, J.M. and ROBERGE, M.C., 2000, Impervious surfaces and the quality
of natural and built environments. White paper prepared for the project Chesapeake Bay
from Space. Available online at: http://chesapeake.towson.edu/landscape/impervious/
download/Impervious.pdf (last accessed 13 February 2008).
BARNES, K.B., MORGAN III, J.M., ROBERGE M.C. and LOWE, S., 2001, Sprawl development: its patterns,
consequences, and measurement. White paper, Towson University. Available online at:
http://chesapeake.towson.edu/landscape/urbansprawl/download/Sprawl_white_paper.pdf
(last accessed 13 February 2008).
BATTY, M., 2000, GeoComputation using cellular automata. In Geocomputation, S. Openshaw
and R.J. Abrahart (Eds), pp. 95–126 (New York: Taylor & Francis).
BATTY, M. and XIE, Y., 1994, From cells to cities. Environment and Planning B, 21, pp. 31–48.
BATTY, M. and XIE, Y., 1997, Possible urban automata. Environment and Planning B, 24,
pp. 175–192.
BENATI, S., 1997, A cellular automaton for the simulation of competitive location. Environment
and Planning B, 24, pp. 205–218.
Census of India, 2001, Table-4, Table-13 and Table-B21, Directorate of Census Operations,
West Bengal.
CLARKE, K.C., HOPPEN, S. and GAYDOS, L., 1997, A self-modifying cellular automaton model of
historical urbanization in the San Francisco Bay area. Environment and Planning B, 24,
pp. 247–261.
DEADMAN, P., BROWN, R.D. and GIMBLETT, P., 1993, Modelling rural residential settlement
patterns with cellular automata. Journal of Environment Management, 37, pp. 147–160.
4746 B. Bhatta
KUMAR, J.A.V., PATHAN, S.K. and BHANDERI, R.J., 2007, Spatio-temporal analysis for mon-
itoring urban growth—a case study of Indore city. Journal of Indian Society of Remote
Sensing, 35, pp. 11–20.
LATA, K.M., RAO, C.H.S., PRASAD, V.K., BADARIANTH, K.V.S., RAHGAVASAMY, V. and RAO,
C.H.S., 2001, Measuring urban sprawl: a case study of Hyderabad. GIS Development, 5,
pp. 26–29.
LI, X. and YEH, A.G.O., 2004, Analyzing spatial restructuring of landuse patterns in a fast
growing region using remote sensing and GIS. Landscape and Urban Planning, 69,
pp. 335–354.
MAKTAV, D. and ERBER, F.S., 2005, Analysis of urban growth using multi-temporal satellite
data in Istanbul, Turkey. International Journal of Remote Sensing, 26, pp. 797–810.
PATHAN, S.K., JOTHIMANI, P., PENDHARKAR, S.P. and KUMAR, D.S., 1989, Urban land use
mapping and zoning of Mombay Metropolitan region using remote sensing data.
Journal of Indian Society of Remote Sensing, 17, pp. 11–22.
PATHAN, S.K., SHUKLA, V.K., PATEL, R.G. and MEHTA, K.S., 1991, Urban land use mapping—a
Downloaded by [Duke University Libraries] at 10:53 11 July 2012
case study of Ahmedabad city and its environs. Journal of Indian Society of Remote
Sensing, 19, pp. 95–112.
PHIPPS, M. and LANGLOIS, A., 1997, Spatial dynamics, cellular automata, and parallel proces-
sing computers. Environment and Planning B, 24, pp. 193–204.
PORTUGALI, J., BENENSON, I. and OMER, I., 1997, Spatial cognitive dissonance and sociospatial
emergence in a self-organizing city. Environment and Planning B, 24, pp. 263–285.
RICHARDSON, H.W., BAE, C.C. and BAXAMUSA, M.H., 2000, Compact cities in developing
countries: assessment and implications. In Compact Cities: Sustainable Urban Forms
for Developing Countries, J. Mike and R. Burgess (Eds), pp. 1–25 (London and New
York: Spon Press).
SCHUELER, T.R., 1994, The importance of imperviousness. Watershed Protection Techniques, 3,
pp. 100–111.
SUDHIRA, H.S., RAMACHANDRA, T.V. and JAGDISH, K.S., 2004, Urban sprawl: metrics, dynamics
and modelling using GIS. International Journal of Applied Earth Observation and
Geoinformation, 5, pp. 29–39.
WHITE, R.W. and ENGELEN, G., 1997, Cellular automaton as the basis of integrated dynamic
regional modeling. Environment and Planning B, 24, pp. 235–246.
WHITE, R.W., ENGELEN, G. and ULJEE, I., 1998, Vulnerability Assessment of Low-Lying Coastal
Areas and Small Islands to Climate Change and Sea Level Rise—Phase 2: Case Study St
Lucia, Kingston, Jamaica. RIKS Publication, Report to the United Nations
Environment Programme, Caribbean Regional Co-ordinating Unit.
YEH, A.G.O. and LI, X., 2001, Measurement and monitoring of urban sprawl in a rapidly
growing region using entropy. Photogrammetric Engineering and Remote Sensing, 67,
pp. 83–90.