Full Quantitative Methods and Socio Economic Applications in GIS 2ed. Edition Wang Ebook All Chapters
Full Quantitative Methods and Socio Economic Applications in GIS 2ed. Edition Wang Ebook All Chapters
Full Quantitative Methods and Socio Economic Applications in GIS 2ed. Edition Wang Ebook All Chapters
com
https://ebookgate.com/product/quantitative-
methods-and-socio-economic-applications-in-
gis-2ed-edition-wang/
https://ebookgate.com/product/socio-economic-applications-of-
geographic-information-science-innovations-in-gis-9-1st-edition-
david-kidner/
https://ebookgate.com/product/dirichlet-and-related-
distributions-theory-methods-and-applications-1st-edition-kai-
wang-ng/
https://ebookgate.com/product/quantitative-eeg-analysis-methods-
and-clinical-applications-1-har-cdr-edition-shanbao-tong/
https://ebookgate.com/product/gis-and-crime-mapping-mastering-
gis-technol-applications-mgmnt-1st-edition-spencer-chainey/
Quantitative Methods In Linguistics 1st Edition Keith
Johnson
https://ebookgate.com/product/quantitative-methods-in-
linguistics-1st-edition-keith-johnson/
https://ebookgate.com/product/methods-of-nonlinear-analysis-
applications-to-differential-equations-2ed-edition-pavel-drabek/
https://ebookgate.com/product/globalization-and-transformations-
of-local-socio-economic-practices-routledge-advances-in-
sociology-1st-edition-ulrike-schuerkens/
https://ebookgate.com/product/advances-in-web-based-gis-mapping-
services-and-applications-1st-edition-songnian-li/
https://ebookgate.com/product/statistical-methods-in-the-
atmospheric-sciences-2ed-edition-wilks-d-s/
Second Edition
Quantitative Methods
and Socio-Economic
Applications in
GIS
Fahui Wang
Second Edition
Quantitative Methods
and Socio-Economic
Applications in
GIS
Second Edition
Quantitative Methods
and Socio-Economic
Applications in
GIS
Fahui Wang
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742
© 2015 by Taylor & Francis Group, LLC
CRC Press is an imprint of Taylor & Francis Group, an Informa business
This book contains information obtained from authentic and highly regarded sources. Reasonable
efforts have been made to publish reliable data and information, but the author and publisher cannot
assume responsibility for the validity of all materials or the consequences of their use. The authors and
publishers have attempted to trace the copyright holders of all material reproduced in this publication
and apologize to copyright holders if permission to publish in this form has not been obtained. If any
copyright material has not been acknowledged please write and let us know so we may rectify in any
future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced,
transmitted, or utilized in any form by any electronic, mechanical, or other means, now known or
hereafter invented, including photocopying, microfilming, and recording, or in any information stor-
age or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copy-
right.com (http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222
Rosewood Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that pro-
vides licenses and registration for a variety of users. For organizations that have been granted a photo-
copy license by the CCC, a separate system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are
used only for identification and explanation without intent to infringe.
Visit the Taylor & Francis Web site at
http://www.taylorandfrancis.com
and the CRC Press Web site at
http://www.crcpress.com
In loving memory of Katherine Z. Wang
To Lei and our three J’s (Jenny, Joshua, and Jacqueline)
Contents
List of Figures .......................................................................................................... xv
List of Tables ...........................................................................................................xix
Foreword .................................................................................................................xxi
Preface.................................................................................................................. xxiii
Author ..................................................................................................................xxvii
List of Major GIS Datasets and Program Files .....................................................xxix
List of Quick References for Spatial Analysis Tasks ............................................xxxi
vii
viii Contents
xv
xvi List of Figures
FIGURE 3.3 Zhuang and non-Zhuang place names in Guangxi, China. ............. 54
FIGURE 3.4 Dialog window for summarization.................................................. 55
FIGURE 3.5 Zhuang place name ratios in Guangxi by the FCA method. ........... 57
FIGURE 3.6 Kernel density of Zhuang place names in Guangxi. ....................... 57
FIGURE 3.7 Spatial interpolation of Zhuang place names in Guangxi by
the IDW method. ............................................................................. 58
FIGURE 3.8 Population change rate in Baton Rouge 2000–2010........................ 63
FIGURE 3.9 Flow chart for implementing the TDW method. ............................. 63
FIGURE 4.1 Constructing Thiessen polygons for five points............................... 71
FIGURE 4.2 Reilly’s law of retail gravitation. ................................................. 72
FIGURE 4.3 Proximal areas for the Cubs and White Sox. .................................. 77
FIGURE 4.4 Probability of choosing the Cubs by the Huff model. ..................... 81
FIGURE 4.5 Proximal areas for public hospitals in Louisiana. ...........................84
FIGURE 4.6 Service areas for public hospitals in Louisiana by Huff model. ..... 86
FIGURE 4.7 Probability of visiting LSUHSC-Shreveport Hospital by Huff
model. .............................................................................................. 87
FIGURE A4.1 Interface for implementing the Huff model. ................................. 91
FIGURE 5.1 Basic floating catchment area method in Euclidean distance..........96
FIGURE 5.2 Two-step floating catchment area method in travel time. ............... 98
FIGURE 5.3 Conceptualizing distance decay in G2SFCA. ............................... 101
FIGURE 5.4 Flow chart for implementing the 2SFCA in ArcGIS..................... 104
FIGURE 5.5 Accessibility to primary care physician in Chicago region by
2SFCA. .......................................................................................... 105
FIGURE 5.6 Accessibility to primary care physician in Chicago region by
2SFCA. .......................................................................................... 106
FIGURE 5.7 Comparison of accessibility scores by the 2SFCA and gravity-
based methods. .............................................................................. 109
FIGURE A5.1 Interface for implementing G2SFCA method. ............................ 112
FIGURE 6.1 Regional growth patterns by the density function
approach. ....................................................................................... 119
FIGURE 6.2 Excel dialog window for regression. ............................................. 122
FIGURE 6.3 Excel dialog window for Format Trendline. ................................. 123
FIGURE 6.4 Illustration of polycentric assumptions. ........................................ 127
List of Figures xvii
FIGURE 6.5 Population density surface and job centers in Chicago. ................ 130
FIGURE 6.6 Density versus distance exponential trend line (census tracts)...... 133
FIGURE 6.7 Density versus distance exponential trend line (survey
townships). ..................................................................................... 136
FIGURE 7.1 Scree plot and variance explained in principal components
analysis. .......................................................................................... 147
FIGURE 7.2 Major steps in principal components analysis and factor
analysis........................................................................................... 148
FIGURE 7.3 Dendrogram for a cluster analysis example. .................................. 150
FIGURE 7.4 Conceptual model for urban mosaic. ............................................. 153
FIGURE 7.5 Districts and subdistricts in Beijing. .............................................. 154
FIGURE 7.6 Spatial patterns of factor scores in Beijing. ................................... 158
FIGURE 7.7 Social areas in Beijing. .................................................................. 159
FIGURE 8.1 Mean centers and ellipses for racial–ethnic groups in the
Chicago area. ................................................................................. 167
FIGURE 8.2 SaTScan dialog windows for point-based spatial cluster
analysis. ..................................................................................... 171
FIGURE 8.3 A spatial cluster of Zhuang place names in Guangxi, China. ....... 172
FIGURE 8.4 ArcGIS dialog window for computing Getis–Ord General G. ..... 173
FIGURE 8.5 Clusters of homicide rates based on local Moran’s Ii. ................... 183
FIGURE 8.6 Clusters of homicide rates based on Gi∗. ........................................ 184
FIGURE 8.7 GeoDa dialog window for defining spatial weights. ..................... 186
FIGURE 8.8 GeoDa dialog window for regression. ........................................... 187
FIGURE 8.9 Standard residuals in the GWR model. ......................................... 188
FIGURE 8.10 Spatial variations of coefficients from the GWR model. ............ 189
FIGURE 9.1 Female breast cancer death rates in Illinois for 2003–2007. ......... 194
FIGURE 9.2 Example of assigning spatial order values to areas. ...................... 197
FIGURE 9.3 Example illustrating REDCAP....................................................... 200
FIGURE 9.4 Interface windows in REDCAP. ....................................................204
FIGURE 9.5 Late-stage breast cancer rates in zip code areas in the Chicago
region in 2000................................................................................206
FIGURE 9.6 Distribution of late-stage breast cancer rates in the Chicago
region in 2000. ...............................................................................207
xviii List of Figures
xix
xx List of Tables
xxi
xxii Foreword
book concludes with a new chapter on spatial statistics that also builds on clustering
with respect to local and more global structures.
Section III of the book treats more advanced topics: regionalization which again
relates to clustering, various linear methods of land use modeling and optimiza-
tion, and finally, a new chapter on modeling of traffic based on Monte Carlo tech-
niques. A very nice feature of the book is the wealth of examples that are included.
Different places such as Louisiana State, Baton Rouge, Louisiana, Columbus, Ohio
and Chicago, Illinois in the United States, and Beijing and Guangxi Province of
China, among other locations are used, while the applications of these various tech-
niques are to population density, retail trade, health care, hospital provision, social
area analysis, the spatial incidence of cancer, and traffic.
In this book, Fahui Wang shows just how far GIS has progressed. Virtually
everything that was developed prior to the GIS age can now be applied, further
developed, and interpreted through the GIS lens, and in some respects, his examples
illustrate the great range and diversity of potential applications which are the marks
of a mature technology. In fact, GIS is becoming part of the routine tool kit that any
analyst would use in studying data that varies across space. What is intriguing about
the treatment here is that the edifice of theory of urban and regional systems that
draws on locational analysis and social physics, although still advancing, is becom-
ing increasingly integrated with GIS, and it is treatment of the subject area such as
that developed here that shows how relevant these tools are to contemporary urban
policy. In fact, it is this focus on policy that marks the book. Those reading it will
find that the author weaves together explanations of the theories involved with their
translation into tools and models and their estimation using straightforward statistics
with notions about different ways of applying these models to real problems that
have strong policy implications. At the end of the day, it is not only understanding
cities better that is the quest for the tools introduced here but also understanding
them in deep enough ways so that effective policies can be advanced and tested that
will provide more sustainable and resilient cities—one of the challenges of near-
and medium-term futures in the socioeconomic domain. This book sets a standard,
shows how this can be achieved, and charts the way forward.
Michael Batty
Centre for Advanced Spatial Analysis
University College, London, United Kingdom
Preface
One of the most important advancements in recent social science research (includ-
ing applied social sciences and public policy) has been the application of quantita-
tive or computational methods in studying the complex human or social systems.
Research centers in computational social sciences have flourished in major uni-
versity campuses including Harvard University (http://www.iq.harvard.edu/),
Stanford University (https://css-center.stanford.edu/), UCLA (http://ccss.ucla.edu/),
University of Washington (http://julius.csscr.washington.edu/), and George Mason
University (http://www.css.gmu.edu/). Many conferences have also been organized
around this theme (http://computationalsocialscience.org/). Geographic Information
System (GIS) has played an important role in this movement because of its capability
of integrating and analyzing various data sets, in particular spatial data. The Center
for Spatially Integrated Social Science (CSISS) at the University of California Santa
Barbara, funded by the National Science Foundation (1999–2007), has been an
important force in promoting the usage of GIS technologies in social sciences. The
Centre for Advanced Spatial Analysis (CASA) at the University College, London,
UK, is also known for its leading efforts in applied GIS and geo-simulation research
with a focus on cities. The growth of GIS has made it increasingly known as geo-
graphic information science (GISc), which covers broader issues such as spatial data
quality and uncertainty, design and development of spatial data structure, social, and
legal issues related to GIS, and many others.
Many of today’s students in geography and other social science-related fields (e.g.,
sociology, anthropology, business management, city and regional planning, public
administration) share the same excitement in GIS. But their interest in GIS may fade
away quickly if the GIS usage is limited to managing spatial data and mapping. In
the meantime, a significant number of students complain that courses on statistics,
quantitative methods, and spatial analysis are too dry and feel irrelevant to their
interests. Over years of teaching GIS, spatial analysis, and quantitative methods, I
have learned the benefits of blending them together and practicing them in case stud-
ies using real-world data. Students can sharpen their GIS skills by applying some
GIS techniques to detecting hot spots of crime, or gain better understanding of a
classic urban land use theory by examining their spatial patterns in a GIS environ-
ment. When students realize that they can use some of the computational methods
and GIS techniques to solve real-world problems in their own field, they become bet-
ter motivated in class. In other words, technical skills in GIS or quantitative methods
are learned in the context of addressing subject issues. Both are important for today’s
competitive job market.
This book is the result of my efforts of integrating GIS and quantitative (com-
putational) methods, demonstrated in various policy-relevant socioeconomic appli-
cations. The applications are chosen with three objectives in mind. The first is to
demonstrate the diversity of issues where GIS can be used to enhance the studies
related to socioeconomic issues and public policy. Applications spread from typical
xxiii
xxiv Preface
themes in urban and regional analysis (e.g., trade area analysis, regional growth
patterns, urban land use and transportation) to issues related to crime and health
analyses. The second is to illustrate various computational methods. Some methods
become easy to use in automated GIS tools, and others rely on GIS to enhance the
visualization of results. The third objective is to cover common tasks (e.g., distance
and travel time estimation, spatial smoothing and interpolation, accessibility mea-
sures) and major issues (e.g., modifiable areal unit problem, rate estimate of rare
events in small populations, spatial autocorrelation) that are encountered in spatial
analysis.
One important feature of this book is that each chapter is task-driven. Methods
can be better learned in the context of solving real-world problems. While each
method is illustrated in a special case of application, it can be used to analyze dif-
ferent issues. Each chapter has one subject theme and introduces the method (or a
group of related methods) most relevant to the theme. For example, spatial regres-
sion is used to examine the relationship between job access and homicide patterns;
systems of linear equations are analyzed to predict urban land use patterns; linear
programming is introduced to solve the problem of wasteful commuting and allo-
cate healthcare facilities; and Monte Carlo technique is illustrated in simulating
urban traffic.
Another important feature of this book is the emphasis on implementation of
methods. All GIS-related tasks are illustrated in the ArcGIS platform, and most sta-
tistical analyses are conducted by SAS. In other words, one may only need access to
ArcGIS and SAS in order to replicate the work discussed in the book and conduct
similar research. ArcGIS and SAS are chosen because they are the leading software
for GIS and statistical analysis, respectively. Some specific tasks such as spatial clus-
tering and spatial regression use free software that can be downloaded from the
Internet. Most data used in the case studies are publicly accessible. Instructors and
advanced readers may use the data sources and techniques discussed in the book to
design their class projects or craft their own research projects. A link to the website
is provided for downloading all data and computer programs used in the case studies
(see “List of Major GIS Datasets and Program Files”). I plan to run a blog under my
homepage (http://ga.lsu.edu/faculty/fahui-wang/) in support of the book.
The book has 12 chapters. Section I includes the first three chapters, covering
some generic issues such as an overview of data management in GIS and basic spa-
tial analysis tools (Chapter 1), distance and travel time measurement (Chapter 2),
and spatial smoothing and interpolation (Chapter 3). Section II includes Chapters
4 through 8, covering some basic quantitative methods that require little or no pro-
gramming skills: trade area analysis (Chapter 4), accessibility measures (Chapter
5), function fittings (Chapter 6), factor analysis (Chapter 7), and spatial statistics
(Chapter 8). Section III includes Chapters 9 through 12, covering more advanced
topics: regionalization (Chapter 9), a system of linear equations (Chapter 10), linear
programming (Chapter 11), and Monte Carlo simulation (Chapter 12). Sections I
and II may serve an upper-level undergraduate course. Section III may be used for
a graduate course. It is assumed that readers have some basic GIS and statistical
knowledge equivalent to one introductory GIS course and one elementary statistical
course.
Preface xxv
Each chapter focuses on one computational method except for the first chapter.
In general, a chapter (1) begins with an introduction to the method, (2) discusses a
theme to which the method is applied, and (3) uses a case study to implement the
method using GIS. Some important issues, if not directly relevant to the main theme
of a chapter, are illustrated in appendixes. Many important tasks are repeated in dif-
ferent projects to reinforce the learning experience (see “Quick Reference for Spatial
Analysis Tasks”).
My interest in quantitative methods has been very much influenced by my doc-
toral advisor, Jean-Michel Guldmann, in the Department of City and Regional
Planning of The Ohio State University. I learned linear programming and solving a
system of linear equations in his courses on static and dynamic programming. I also
benefited a great deal from the mentorship of Donald Haurin in the Department of
Economics of The Ohio State University. The topics on urban and regional density
patterns and wasteful commuting can be traced back to his inspiring urban econom-
ics course. Philip Viton, also in the Department of City and Regional Planning of
The Ohio State University, taught me much of the econometrics. I only wish I could
have been a better student then.
It has been 8 years since the publication of the previous edition of this book.
Several friends and many users, especially my students, have found some errors and
given me valuable feedback. Several tasks in the previous edition had to be imple-
mented in the deprecated ArcInfo workstation environment (and its associated AML
program). The software, particularly ArcGIS, has advanced with increasing capaci-
ties and more user-friendly interfaces. Numerous people have urged me to automate
some popular tasks in Python. I have also collected a few new tricks from my own
research that I hope to share with a broader audience. The reasons for preparing a
new version just kept piling on. It came to a point that I could not find any more
excuses for not doing it when my sabbatical request was approved by Louisiana State
University in early 2013.
Inevitably I have lost a step or two over the years. The added administrative duty
did not help either. I truly felt that the revision took more work than the previous edi-
tion. I used most of my summer and the sabbatical in the fall of 2013 on the book,
and often found myself working on the book in my office very early in the mornings
of weekdays and weekends for much of the spring semester of 2014. In this version,
several popular tasks such as trade area delineation, spatial accessibility measures
(2SFCA), mixed-level regionalization (MLR), and Garin–Lowry model are auto-
mated as a convenient toolkit in ArcGIS; and some common tasks such as REDCAP
regionalization, wasteful commuting measure, and Monte Carlo simulation are now
implemented in more user-friendly programs. All case studies have been tested mul-
tiple times, and instructions were based on ArcGIS 10.2. I would like to believe that
it was the amount of revision work that merited the investment of so much of time.
It could be an understatement to say that more than 60% of the materials (including
case studies) are new. I hope that the readers will be convinced likewise that it is a
worthy cause.
I am so grateful for the generous help from several individuals in preparing this
second edition. Yujie Hu coded the program TSME used in the case study in Chapter
12 and coauthored Chapter 12. He also tested all case studies and corrected numerous
xxvi Preface
errors. Haojie Zhu implemented the Python programs for three toolkits (trade area
delineation, spatial accessibility measures, and Garin–Lowry model). Both Yujie
and Haojie are PhD students in the Department of Geography and Anthropology,
Louisiana State University. Lan Mu in the Department of Geography, University of
Georgia, implemented the mixed-level regionalization method in Python and coau-
thored Appendix 9B in Chapter 9. Xinyue Ye in the Department of Geography, Kent
State University, coauthored Appendix 10D on cellular automata for urban land use
modeling in Chapter 10. Carmi J. Neiger of Elmhurst College tested the case stud-
ies in Chapters 1 through 4 and edited the four chapters. I thank Michael Batty for
graciously writing the new Foreword on short notice.
Finally, I would like to thank the editorial team at Taylor & Francis: acquisi-
tion editor, Irma Britton; production coordinator, Laurie Schlags; and many others,
including typesetters, proofreaders, cartographers, and computer specialists. Thank
you all for guiding me through the whole process.
This book intends to mainly serve students in geography, urban and regional
planning, public policy, and related fields. It can be used in courses such as (1) spa-
tial analysis, (2) location analysis, (3) quantitative methods in geography, and (4)
applications of GIS in business and social science. The book can also be useful for
social scientists in general with research interests related to spatial issues. Some in
urban economics may find the studies on urban structures and wasteful commut-
ing relevant, and others in business may think the chapters on trade area analysis
and accessibility measures useful. The case study on crime patterns may interest
criminologists, and the one on regionalization and cancer analysis may find the
audience among epidemiologists. Some basic GIS knowledge (e.g., one introductory
GIS course) will help readers navigate through the case studies smoothly, but is not
required.
Additional material including the datasets and program files is available from the
CRC website: http://www.crcpress.com/product/isbn/9781466584723.
Author
Fahui Wang is James J. Parsons Professor and
Chair of the Department of Geography and
Anthropology, Louisiana State University. He
earned a BS degree in geography from Peking
University, China, an MA degree in economics and
a PhD degree in city and regional planning from
The Ohio State University. His research interests
include GIS applications in human geography
(urban, economic, and transportation), city and
regional planning, and public policy (crime and
health). His work has been supported by the National
Science Foundation, U.S. Department of Energy,
U.S. Department of Health and Human Services
(Agency for Healthcare Research & Quality and
the National Cancer Institute), U.S. Department of
Housing and Urban Development, U.S. Department of Justice (National Institute of
Justice and Office of Juvenile Justice and Delinquency Prevention), and the National
Natural Science Foundation of China. He is on the editorial boards of several inter-
national journals including Annals of the Association of American Geographers,
Annals of GIS, Applied Geography, and the Chinese Geographical Science. He has
published over 100 refereed articles.
xxvii
List of Major GIS Datasets
and Program Files*
Case Study Data Folder GIS Dataset Program File Study Area
1 BatonRouge BR.gdb, Census. East Baton Rouge
gdb Parish
2 Louisiana LA_State.gdb, Louisiana
LAState.gdb
3A China_GX GX.gdb Guangxi, China
3B BatonRouge BR.gdb East Baton Rouge
Parish
4A Chicago ChiRegion.gdb Huff Model.tbx Chicago Region
4B Louisiana LA_State.gdb, Huff Model.tbx Louisiana
LAState.gdb
5 Chicago ChiRegion.gdb Accessibility.tbx Chicago Region
6 Chicago ChiUrArea.gdb monocent.sas, polycent.sas, Chicago Urban Area
mono_twnshp.sas
7 Beijing BJSA.gdb PCA_FA_CA.sas, BJreg.sas Beijing, China
8A Chicago ChiUrArea.gdb Chicago Urban Area
8B China_GX GX.gdb Guangxi, China
8C Chicago ChiCity.gdb PCA_FA.sas Chicago City
9 Chicago/ Chizip shapefiles poisson_zip.sas, Mixed- Chicago Region
Chizip Level Regionalizaiton
Tools.tbx
10 SimuCity SimuCity.gdb Garin-Lowry Model.tbx Hypothetical city
11A Columbus Columbus.gdb WasteComm.R, LP.sas Columbus, Ohio
11B BatonRouge BR.gdb East Baton Rouge
Parish
12 BatonRouge BRMSA.gdb TSME.exe, Features To Baton Rouge MSA
Text File.tbx
* http://www.crcpress.com/product/isbn/9781466584723
xxix
List of Quick References for
Spatial Analysis Tasks
Section First Figure Sections (Steps)
Task Introduced (Step) Illustration Repeated
2SFCA method (automated toolkit) Appendix 5B Figure A5.1
Areal weighting Interpolation Section 3.5.1 (1–3)
Attribute join Section 2.4.1 (5) Figure 2.4 Section 3.3.1 (3),
Section 4.3.2 (6)
Attribute query Section 1.3.1 (1) Figure 1.1 Section 3.3.1 (4)
Dissolve Section 1.3.2 (10) Figure 1.8
Euclidean distance matrix Section 2.4.1 (3) Section 3.3.1 (2),
computation Section 4.3.2 (5)
Geocoding by address locator Section 2.4.1 (1) Figure 2.3
Geocoding by geographic coordinates Section 2.4.1 (1) Figure 2.2
Geographically weighted regression Section 8.7.2 (8)
(GWR)
Getis–Ord general G Section 8.7.1 (3)
Huff model (automated toolkit) Appendix 4B Figure A4.1
IDW spatial interpolation Section 3.3.2 (10) Section 6.5.1 (2)
Kernel density estimation (KDE) Section 3.3.2 (9) Section 6.5.1 (2)
Local Gi∗ Section 8.7.1 (4)
Local Moran’s Ii (LISA) Section 8.7.1 (4)
Mean center Section 8.2 (1)
Moran’s I Section 8.7.1 (3)
Network dataset building Section 2.4.2 (8) Section 10.3 (1, 3)
OLS regression in ArcGIS Section 8.7.2 (5) Section 6.5.1 (4–5)
Projection of spatial data Section 1.3.1 (3) Figure 1.3
Proximal area by Euclidean distance Section 4.3.1 (2–3)
Proximal area by travel time Section 4.4.1 (1–3)
Spatial join Section 1.3.2 (13) Figure 1.10
Spatial query Section 1.3.1 (2) Figure 1.2 Section 6.5.1 (1)
Spatial regression in GeoDa Section 8.7.2 (7)
Standard deviational ellipse Section 8.2 (3)
Standard distance Section 8.2 (2)
Summarize attribute Section 3.3.1 (5) Figure 3.4 Section 4.3.1 (4),
Section 4.3.2 (7)
Travel time O-D matrix computation Section 2.4.2 Section 10.3 (1, 3),
(9–12) Section 11.2.2 (3)
Weighted centroids Section 4.3.1 (1)
xxxi
Section I
GIS and Basic Spatial
Analysis Tasks
1 Getting Started with
ArcGIS
Data Management and Basic
Spatial Analysis Tools
3
4 Quantitative Methods and Socioeconomic Applications in GIS
on command line interface. Prior to the release of ArcGIS, ESRI had an entry-level
GIS package with a GUI called ArcView. ArcView is now referred to as ArcGIS for
Desktop Basic.
ArcGIS contains three major modules: ArcCatalog, ArcMap, and ArcToolbox.
ArcCatalog views and manages spatial data files. ArcMap displays, edits, and ana-
lyzes the spatial data as well as attribute data. ArcToolbox contains various task-
specific tools for data analysis such as data conversion, data management, spatial
analysis, and spatial statistics. ArcCatalog and ArcMap can be launched sepa-
rately from the Start button on the Windows taskbar, and either module can also be
accessed by clicking its distinctive symbol inside the other modules. ArcToolbox is
accessed in either ArcMap or ArcCatalog. ArcMap is used most often, and thus is
the module being referred to in the text unless otherwise noted.
Traditionally, a GIS uses either the vector or the raster data model to manage
spatial data. A vector GIS uses geographically referenced points to construct spatial
features of points, lines and areas; and a raster GIS uses grid cells in rows and col-
umns to represent spatial features. The vector data structure is used in most socio-
economic applications, and also in most case studies in this book. The raster data
structure is simple and widely used in environmental studies, physical sciences,
and natural resource management. Most commercial GIS software can convert
from vector into raster data, or vice versa. In ArcGIS, the tools are available under
ArcToolbox > Conversion Tools.
Earlier versions of ESRI’S GIS software used the coverage data model. Later,
shapefiles were developed for the ArcView package. Since the release of ArcGIS
8, the geodatabase model has become available, and represents the current trend of
object-oriented data models. The object-oriented data model stores the geometries
of objects (spatial data) also with the attribute data, whereas the traditional cover-
age or shapefile model stores spatial and attribute data separately. Spatial data used
for case studies in this book are provided almost entirely in geodatabases (with rare
exceptions such as for network datasets). Spatial and attribute data in socioeconomic
analysis often come from different sources, and a typical task is to join them together
in GIS for mapping and analysis. This involves attribute data management as dis-
cussed below.
in a table is done in ArcMap: on the open table > right-click the field, choose Field
Calculator to access the tool. In addition, basic statistics for a field can be obtained
in ArcMap by right-clicking the field and choosing Statistics. Step 5 in Section 1.3
demonstrates how to create a new field and update its values.
In GIS, an attribute join is often used to link information in two tables based
on a common field (key). The table can be an attribute table associated with a par-
ticular geodataset or a standalone table. In an attribute join, the names of the fields
used to link the two tables do not need to be identical in the two tables to be com-
bined, but the data types must match. There are various relationships between tables
when linking them: one-to-one, many-to-one, one-to-many, and many-to-many. For
either one-to-one or many-to-one relationship, a join is used to combine two tables
in ArcGIS. However, when the relationship is one-to-many or many-to-many, a join
cannot be performed. ArcGIS uses a relate to link two tables while keeping the two
tables separated. In a relate, one or more records are selected in one table, and the
associated records are selected in another table. Table 1.1 summarizes the relation-
ships and corresponding tools in ArcGIS.
A join or relate is performed in ArcMap. Under Table of Contents (TOC), right-
click the spatial dataset or the table that is to become the target table, and choose
Joins and Relates > Join (or Relate) > In the Join Data dialog window, choose “Join
attributes from a table,” and select (1) the common field in the target table, (2) the
source table, and (3) the common field in the source table. It is important to under-
stand the different roles of the target and source tables in the process. As a result, the
target table is expanded by adding the fields from the source table.
A join is temporary in the sense that no new data are created and the join rela-
tionship is lost once the project is exited without being saved. To preserve the result-
ing layer (or combined table), right-click it under Table of Contents > Data > Export
Data (or Export) and to save it in a new spatial dataset (table). If it is desirable to
only save the joined attribute table of a spatial dataset, in the open attribute table,
click the Table Options dropdown icon > Export. The latter can be termed “per-
manent join.” Alternatively, one may use the tools in ArcToolbox to implement
attribute joins: ArcToolbox > Data Management Tool > Joins > choose Add Join
for a temporary join or Join Field for a permanent join. These ArcToolbox tools
TABLE 1.1
Types of Relationships in Combining Tables
Relationship Match Join or Relate in ArcGIS
One-to-one One record in the target table → one record in Join
the source table
Many-to-one Multiple records in the target table → one Join
record in the source table
One-to-many One record in the target table → multiple Relate
records in the source table
Many-to-many Multiple records in the target table → multiple Relate
records in the source table
Getting Started with ArcGIS 7
are particularly useful for automating the process in programming (see step 8 in
Section 5.4.1).
Once the attribute information is joined to a spatial layer, mapping the data is
also done in ArcMap. Right-click the layer and choose Properties > Symbology to
invoke the dialog window. In the dialog, one can select a field to map, choose colors
and symbols, and plan the layout. From the main menu bar, choose View > Layout
View to see what the draft map looks like. Map elements (title, scale, north arrow,
legends, and others) can be added to the map by clicking Insert from the main
menu bar.
of each county (e.g., FIPS code, name, administrator) are assigned to schools that fall
within the county boundary. If multiple objects in the source layer correspond to one
object in the target layer, two operations may be performed: “summarized join” and
“distance join.” A summarized join summarizes the numeric attributes of features
in the source layer (e.g., average, sum, minimum, maximum, standard derivation, or
variance) and adds the information to the target layer. A distance join identifies the
closest object (out of many) in the source layer from the matching object in the target
layer, and transfers all attributes of this nearest object (plus a distance field showing
how close they are) to the target layer. For example, one may spatially join a point
layer of geocoded crime incidents (source layer) with a polygon layer of census tracts
(target layer), and generate aggregated crime counts in those census tracts.
There are a variety of spatial joins between different spatial features (Price,
2004, pp. 287–288). Table 1.2 summarizes all types of spatial joins in ArcGIS.
Spatial joins are accessed in ArcMap similar to attribute joins: right-click the source
layer > choose Joins and Relates > Join. In the Join Data dialog window, choose
“Join data from another layer based on spatial location” instead of “Join attributes
from a table.” In Section 1.3.2, step 13 illustrates the use of a spatial join, in particu-
lar, a case of a summarized join.
Map overlays may be broadly defined as any spatial analysis involving modifying
features from different layers. The following reviews some of the most commonly
used map overlay tools: Clip, Erase, Intersect, Union, Buffer, and Multiple Ring
Buffer. A Clip truncates the features from one layer using the outline of another. An
Erase trims away the features from the input layer that fall inside the erase layer.
An Intersect overlays two layers and keeps only the areas that are common to both.
A Union also overlays two layers but keeps all the areas from both layers. A Buffer
creates areas by extending outward from point, line, or polygon features over a speci-
fied distance. A Multiple Ring Buffer generates buffer features based on a set of
distances. In ArcGIS 10.2, the above map overlay tools are grouped under different
toolsets in ArcToolbox > Analysis Tools: Clip is under the Extract toolset; Erase,
Intersect and Union are under the Overlay toolset; and Buffer and Multiple Ring
Buffer are grouped under the Proximity toolset. Other common map overlay tools
used in case studies elsewhere in the book include Point Distance (see step 3 in
Section 2.4.1) and Near (see step 2 in Section 4.3.1). In Section 1.3.2, step 7 uses the
multiple ring buffer tool, and step 8 uses the intersect tool.
One may notice the similarity among spatial queries, spatial joins, and map over-
lays. Indeed, many spatial analysis tasks may be accomplished by any one of the
three. Table 1.3 summarizes the differences between them. A spatial query only finds
the information on the screen and does not create new datasets (unless one chooses
to export the selected records or features). A spatial join always saves the result in a
new layer. There is an important difference between spatial joins and map overlays.
A spatial join merely identifies the location relationship between spatial features of
input layers, and does not change existing spatial features or creates new features. In
the process of a map overlay, some of the input features are split, merged, or dropped
in generating a new layer. In general, a map overlay operation takes more computation
time than a spatial join, and a spatial join takes more time than a spatial query.
TABLE 1.2
Types of Spatial Joins in ArcGIS
Source Target
Layer (S) Layer (T) Simple Join Distance Join Summarized Join
Point Point For each point in T, find its closest point in S For each point in T, find all the points in S closer to this
and transfer attributes of that closest point to T point than to any other point in T and transfer the points’
summarized attributes to T
Getting Started with ArcGIS
Line For each point in T, find the closest line in S For each point in T, find all lines in S that intersects it and
and transfer attributes of that line to T transfer the lines’ summarized attributes to T
Polygon For each point in T, find the polygon in S containing For each point in T, find its closest polygon in S
the point, and transfer the polygon’s attributes to T and transfer attributes of that polygon to T
Point Line For each line in T, find its closest point in S and For each line in T, find all the points in S that either
transfer the point’s attributes to T intersect or lie closest to it and transfer the points’
summarized attributes to T
Line For each line in T, find the line in S that it is part of For each line in T, find all the lines in S that intersect it and
and transfer the attributes of source line to T transfer the summarized attributes of intersected lines to T
Polygon For each line in T, find the polygon in S that it falls For each line in T, find the polygon in S that it For each line in T, find all the polygons in S crossed by it
completely inside, and transfer the polygon’s is closest to and transfer the polygon’s and transfer the polygons’ summarized attributes to T
attributes to T attributes to T
Point Polygon For each polygon in T, find its closest point in S For each polygon in T, find all the points in S that fall inside
and transfer the point’s attributes to T it and transfer the points’ summarized attributes to T
Line For each polygon in T, find its closest line in S For each polygon in T, find all the lines in S that intersect it
and transfer the line’s attributes to T and transfer the lines’ summarized attributes to T
Polygon For each polygon in T, find the polygon in S that it For each polygon in T, find all the polygons in S that
falls completely inside and transfer the attributes intersect it and transfer the summarized attributes of
of source polygon to T intersected polygons to T
9
10 Quantitative Methods and Socioeconomic Applications in GIS
TABLE 1.3
Comparison of Spatial Query, Spatial Join, and Map Overlay
Basic Spatial Whether a New Whether New Computation
Analysis Tools Function Layer Is Created Features Are Created Time
Spatial query Finds information based No (unless the No Least
on location relationship selected features
between features from are exported to a
different layers and new dataset)
displays on screen
Spatial join Identifies location Yes No Between
relationship between
features from different
layers and transfers the
attributes to target layer
Map overlay Overlays layers to create Yes Yes (splitting, merging Most
new features and saves or deleting features
the result in a new layer to create new)
Throughout this book, all computer file names, field names, and also some com-
puter commands are in the Courier New font. The following is a step-by-step
guide for the project.
In many cases, GIS analysts need to extract a study region that is contained in a
larger area. This project begins with deriving the census tracts in East Baton Rouge
Parish (“parish” is a county-equivalent unit in Louisiana). The task is accomplished
by steps 1 and 2.
Step 1. Using an attribute query to extract the parish boundary: In ArcMap, click
the Add Data tab to add all three features in geodatabase Census.gdb and the
feature BRCenter in geodatabase BR.gdb to the project.
From the main menu, choose Selection > Select By Attributes >. In the dialog
window (Figure 1.1), make sure that the Layer is County; in the bottom box, input
NAMELSAD10 = 'East Baton Rouge Parish'*; and click Apply (or OK). The
parish is selected and highlighted on the map. The selection is based on the attributes
of a GIS dataset, and thus it is an attribute query.
In the layers window under Table Of Contents (TOC), right-click the layer
County > Data > Export Data >. In the dialog window, make sure that Selected
* For convenience and accuracy, input a field name such as “NAMELSAD10” by double clicking the
field from the top window. The parish name “East Baton Rouge Parish” is case sensitive.
12 Quantitative Methods and Socioeconomic Applications in GIS
features are in the top dropdown bar; under Output feature class, navigate to geoda-
tabase BR.gdb and name it as EBRP; and click OK.
Click Yes to add the exported data to the map as a new layer.
Step 2. Using a spatial query to extract census tracts in the parish: From the main
menu of ArcMap, choose Selection > Select By Location >. In the dialog window
(Figure 1.2), check Tract under Target layer(s), choose EBRP as Source layer, choose
“are within the source layer feature” under Spatial selection method for target layer
feature(s), and click OK. All census tracts within the parish are selected. The selection
is based on the spatial relationship between two layers, and thus it is a spatial query.
Similar to step 1, return to the layers window, and right-click the layer
Tract > Data > Export Data > export the selected features to a new feature BRTrt
in geodatabase BR.gdb. Use the Add Data tab to add BRTrt to the map.
Step 3. Transforming the spatial dataset to a projected coordinate system: In the
layers window, right-click the layer BRTrt > Properties >. In the Layer Properties
window, click Source from the top bar. It shows the Extent of the feature in decimal
degrees (dd) and its Geographic Coordinate System (GCS_North_American_1983).
Getting Started with ArcGIS 13
Repeat on other layers (Tract, County and State) to confirm that they are all
in a geographic coordinate system, as are all TIGER data. Now check the layer
BRCenter. It has the Extent measured in meters (m) and a Projected Coordinate
System (NAD_1983_UTM_Zone_15N).*
For reasons discussed previously, we need to transform the census tracts dataset
BRTrt in an unprojected geographic coordinate system to one in a projected coor-
dinate system. For convenience, we use the projection system that is predefined in
the layer BRCenter.
Open ArcToolbox by clicking the icon on the menu bar above the ArcMap
palette. In ArcToolbox, go to Data Management Tools > Projections and
Transformations > double-click the tool Project. In the Project dialog window (as
shown on the left of Figure 1.3), choose BRTrt for Input Dataset or Feature Class;
* In ArcMap, as we move the cursor in a map, the bottom right corner displays the corresponding coor-
dinates in the current map units, in this case, Decimal Degrees. In other words, maps in the data (or
layout) view are in geographic coordinates, which is set by the coordinate system of the first layer
added to the project.
14 Quantitative Methods and Socioeconomic Applications in GIS
name the Output Dataset or Feature Class as BRTrtUtm (make sure that it is in
geodatabase BR.gdb); under Output Coordinate System, click the graphic icon
to activate the Spatial Reference Properties dialog window >. In this window (as
shown on the right of Figure 1.3), click the Add Coordinate System tab and choose
Import; browse to select BRCenter and click Add*; and click OK to close the dia-
log >. Click OK again to execute the Project tool. Use the Add Data tab to add
the feature BRTrtUTM to the map.
Check the coordinate system of the new dataset BRTrtUtm.
Step 4. Calculating area in census tracts: Under File on the top-level menu, click
New to exit the project, click OK when the New Document dialog appears, and click
No to “Save changes to the Untitled.” Now first add the feature BRTrtUtm obtained
from step 3 and then the other layers. Maps now change to the projected coordinate
system in meters.
Under Table Of Contents, right-click the layer BRTrtUtm >. Open Attribute
Table >. In the table, click the Table Options tab (the first at the upper left cor-
ner) > Add Field >. In the Add Field dialog window, enter Area for Name, and
choose Double for Type, and click OK. A new field Area is now added to the table
as the last column.
In the open table, scroll to the right side of the table, right-click the field
Area > choose Calculate Geometry, and click Yes on the warning message >. In the
Calculate Geometry dialog window, choose Area for Property and “Square Meters
[sqm]” for Units, and click OK to update areas (click Yes again on the warning
message).
Step 5. Calculating population density in census tracts: Still on the open attribute
table of BRTrtUtm, follow the above procedure to add a new field PopuDen (choose
* Alternatively, expand Layers in the middle box > choose “NAD_1983_UTM_Zone_15N,” which is the
projection used in the layer BRCenter.
Getting Started with ArcGIS 15
Float for Type) >. Right-click the field PopuDen > choose Field Calculator >. In the
dialog window (Figure 1.4), under “PopuDen =,” input the formula [DP0010001]/
([Area]/1000000), and click OK to update the field.* Note that the map projection
unit is meter, and thus the area unit is square meter. The dominator in the formula
converts area unit in square meter into square kilometer, the numerator DP0010001
is total population, and the unit of population density is persons per square kilometer.
(In project instructions in other chapters of this book, the bracket signs for each vari-
able may be dropped in a formula for simplicity, for example, the above formula may
be written as PopuDen = DP0010001/(Area/1000000).)
It is often desirable to obtain some basic statistics for a field. Still on the open
table, right-click the field PopuDen > choose Statistics. The result includes count,
minimum, maximum, sum, mean, standard deviation, and a frequency distribution.
Close the table window.
Step 6. Mapping population density across census tracts: In the layers window,
right-click the layer BRTrtUtm > choose Properties. In the Layer Properties win-
dow, click Symbology > Quantities > Graduated Colors to map the field PopuDen
* Similar to step 1, field names such as [DP0010001] and [Area] in the formula are entered by double-
clicking the field names from the top box to save time and minimize typos.
16 Quantitative Methods and Socioeconomic Applications in GIS
(as shown in Figure 1.5). Click Apply to see the changes. Experiment with different
classification methods, number of classes, and color schemes.
From the main menu bar, choose View > Layout View to preview the map. Also
from the main menu bar, choose Insert > Legend (Scale Bar, North Arrow, and
others) to add map elements. When you are satisfied with the map, go to the main
menu bar > File > Export Map to save the map product in a chosen format. The
population density map shown in Figure 1.6 uses a manual classification for density
classes.
Save the project with the final map by going to the main menu bar > File > Save as.
They rested there, and in the Moe ihola lakou a ao aela, hele
morning Kekuaokalani came to maila o Kekuaokalani a ma ke
the presence of Hoapili and alo o Hoapili a me Naihe, olelo
Naihe. They then asked, “Are we aku laua: “O ka holo keia o
going?” Kekuaokalani said “Yes,” kakou?” ae aku o Kekuaokalani,
but his assent was not real. He “Ae.” Aka, aole ona ae io, he
was bent on war. Wherefore manao kaua kona. Noia mea,
Naihe and Hoapili said: “The olelo aku o Naihe me Hoapili:
navel is cut then, companion.” “Wehe i ka piko la e na
hoahanau.”
OF HEMA. NO HEMA.
Hema was one of the most Oia kekahi koa kaulana loa i loko
famous warriors in the days of o ko Kamehameha mau la, a he
Kamehameha, and a man who kanaka makau ole i ko hai koa a
was not afraid of the bravery me ka ikaika. O Hema, aole oia i
and strength of others. Hema ao ia i ke koa a me ke kaua, aole
was not trained to be a warrior, i maa ma na hoouka kaua, aole
or in the art of war. He was not no hoi oia he koa. Aka, he
accustomed to the waging of aipuupuu o Hema na
wars, and was never a warrior. Kamehameha, o kana hana o ka
He was a steward of lawe a me ka malama i wahi ai
Kamehameha, and his duties na ke ’lii i na la a pau loa. I ka
were to take and keep in charge wa hele mao a mao, i ka wa
rations for the king every day. kaua, oia ka mea lawe ai, i ka
While going abroad in war times wa e kaua ai o Kamehameha a
he was the bearer of food, and pololi, alaila lawe aku o Hema i
when a battle was in progress ka ai a me kahi ia.
and Kamehameha became
hungry, then Hema would bring
him provisions.
In all these works pertaining to a Ma keia mau hana a pau loa i pili
steward Hema was an expert, i ka aipuupuu, ua makaukau loa
and was satisfactory to o Hema, a ua kupono i ko
Kamehameha, but he had not Kamehameha makemake, aka,
received a chiefly term, nor aole i loaa kona inoa alii, a me
fame, nor was he a favorite, but ke kaulana a me ka punahele. A
on the day that he chose to be a i kona la i lalau ai i ke koa a luku
warrior and destroyed the aku i na enemi o Kamehameha,
enemies of Kamehameha, that ia la oia i lilo ai i alii, a punahele
day he became a chief and a na Kamehameha, a pau kona
favorite of Kamehameha, and lawe ana i ka oihana aipuupuu.
abandoned his stewardship.
OF NALU. NO NALU.
[486]
Hala.
Of Kekuawahine. No Kekuawahine.
OF MAKAIOULU. NO MAKAIOULU.
Makaioulu was one of Oia kekahi koa kaulana o
Kamehameha’s celebrated Kamehameha, i ke kaua ana
warriors in the war between Oahu nei me Kamehameha ma
Oahu and Kamehameha, at Nuuanu. O Kaalamakaoikuwa ko
Nuuanu. Kaalamakaoikuwa was Oahu nei koa, e noho ana ia ma
the warrior of Oahu, and resided Luahenewai, ma Waikiki kai.
at Luahenewai, Waikiki-kai. Hele maila o Makaioulu me kona
Makaioulu with a companion hoahele, he koa no, o
warrior, Naaimokuokama by Naaimokuokama, a halawai laua
name, came and met me Kaalamakaoikuwa, ia wa,
Kaalamakaoikuwa. Instantly hopu ia maila o Makaioulu a paa,
Makaioulu was seized and held holo o Naaimokuokama. Ia
prisoner, while Naaimokuokama Naaimokuokama e holo ana,
ran away. And while he was kahea mai o Makaioulu: “O ka
running, Makaioulu called after holo ka kau, haalele oe ia’u?” ia
him: “Do you run away and kahea ana, ku o
leave me?” When Naaimokuokama. Olelo aku o
Naaimokuokama heard this he Makaioulu, “O ko ihe e paa ala i
stopped. Makaioulu then said: ko lima, pahu oe a pololei i ku’u
“Hurl that spear that you hold in piko.” Ia wa, pahu o
your hand straight to my navel.” Naaimokuokama i kana ihe. Ma
Naaimokuokama did so, and at keia pahu ana, e alo ae ana o
the same time that he did, Makaioulu loaa o
Makaioulu dodged to one side, Kaalamakaoikuwa i ka ihe, a ku
the spear striking ihola, a make loa, pakele aku la
Kaalamakaoikuwa and killing o Makaioulu, a hele aku la laua
him. Makaioulu escaped, and mai laila aku, a hiki ma
they went thence to Puowaina. Puowaina.
Here were ten soldiers, who Ilaila e noho ana he umi koa, ike
when they saw the two, pointed maila lakou ia laua, hou like
their ten spears at them at the maila lakou i na ihe he umi. Ia
same time. And while they were lakou e hou ana, huli ae la o
doing this, Makaioulu turned Makaioulu i hope, a olelo aku i
rearward and said to his kona kokoolua: “E! hou lakou la i
companion: “Say, if they hurl ku’u alo, kunihi oe; i hou lakou i
their spears to my front, you ku’u aoao, ma ka akau oe o’u e
stand sideways; and if they hurl ku ai, mai imo oe o ku oe.”
at my side you stand at my right. Alaila, huli aku la o Makaioulu
Do not wink or you will be hit.” hele imua, a hiki ia io lakou la he
Makaioulu then faced about and umi, ia lakou e hou ana, alo ae
went forward. And when he la o Makaioulu hala, a hala ka
came near to the ten men they lakou ihe, e huelepo iho ana o
hurled spears at him. He dodged Makaioulu i kana laau palau, loaa
and the spears missed him. As eono kanaka, pakele aku eha.
their spears missed him,
Makaioulu made a sweep with
his war club which caught six
men, four escaping.