forked from mgimond/Spatial
-
Notifications
You must be signed in to change notification settings - Fork 0
/
01-intro.Rmd
108 lines (70 loc) · 9.91 KB
/
01-intro.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
# Introduction to GIS {#introGIS}
```{r echo = FALSE, message = FALSE}
#source('libs/Common.R')
```
## What is a GIS?
A Geographic Information System is a multi-component environment used to create, manage, visualize and analyze data and its spatial counterpart. It's important to note that most datasets you will encounter in your lifetime can all be assigned a spatial location whether on the earth's surface or within some arbitrary coordinate system (such as a soccer field or a gridded petri dish). So in essence, any dataset can be represented in a GIS: the question then becomes "does it need to be analyzed in a GIS environment?" The answer to this question depends on the purpose of the analysis. If, for example, we are interested in identifying the ten African countries with the highest conflict index scores for the 1966-78 period, a simple table listing those scores by country is all that is needed.
```{r 01-Africa-conflict, results='asis', message=FALSE, echo=FALSE}
library(xtable)
library(spdep)
data(afcon)
x <- data.frame(Country = afcon$name, Conflicts = afcon$totcon)
# Split dataframe into two parts
x <- x[ order(-x$Conflicts),]
x2 <- split(x, c(rep(1,floor(nrow(x)/2)), rep(2,floor(nrow(x)/2))))
x3 <- as.data.frame(x2)
colnames(x3) <- rep(c("Country", "Conflicts"),2)
rownames(x3) <-NULL
# Print output
knitr::kable(x3, booktabs=TRUE, caption="Index of total African conflict for the 1966-78 period [@Anselin1992a].")
```
*Data source: Anselin, L. and John O'Loughlin. 1992. Geography of international conflict and cooperation: spatial dependence and regional context in Africa. In The New Geopolitics, ed. M. Ward, pp. 39-75.*
A simple sort on the Conflict column reveals that `r as.vector(x[1:10,]$Country)` are the top ten countries.
What if we are interested in knowing whether countries with a high conflict index score are geographically clustered, does the above table provide us with enough information to help answer this question? The answer, of course, is no. We need additional data pertaining to the geographic location and shape of each country. A map of the countries would be helpful.
```{r map-africa, echo=FALSE, fig.cap = "Choropleth representation of African conflict index scores. Countries for which a score was not available are not mapped.", fig.width=7.5, fig.height=7.5}
knitr::include_graphics("img/Africa_conflicts.png")
```
Maps are ubiquitous: available online and in various print medium. But we seldom ask how the boundaries of the map features are encoded in a computing environment? After all, if we expect software to assist us in the analysis, the spatial elements of our data should be readily accessible in a digital form. Spending a few minutes thinking through this question will make you realize that simple tables or spreadsheets are not up to this task. A more complex data storage mechanism is required. This is the core of a GIS environment: a spatial database that facilitates the storage and retrieval of data that define the spatial boundaries, lines or points of the entities we are studying. This may seem trivial, but without a spatial database, most spatial data exploration and analysis would not be possible!
### GIS software
Many GIS software applications are available--both commercial and open source. Two popular applications are **ArcGIS** and **QGIS**.
#### ArcGIS
A popular commercial GIS software is [**ArcGIS**](http://desktop.arcgis.com/en/) developed by ESRI ([ESRI](https://en.wikipedia.org/wiki/Esri), pronounced *ez-ree*),was once a small land-use consulting firm which did not start developing GIS software until the mid 1970s. The ArcGIS desktop environment encompasses a suite of applications which include ArcMap, ArcCatalog, ArcScene and ArcGlobe. ArcGIS comes in three different license levels (basic, standard and advanced) and can be purchased with additional *add-on* packages. As such, a single license can range from a few thousand dollars to well over ten thousand dollars. In addition to software licensing costs, ArcGIS is only available for Windows operating systems; so if your workplace is a Mac only environment, the purchase of a Windows PC would add to the expense.
### QGIS
A very capable open source (free) GIS software is [**QGIS**](http://qgis.org). It encompasses most of the functionality included in ArcGIS. If you are looking for a GIS application for your Mac or Linux environment, QGIS is a wonderful choice given its multi-platform support. Built into the current versions of QGIS are functions from another open source software: **GRASS**. GRASS has been around since the 1980's and has many advanced GIS data manipulation functions however, its use is not as intuitive as that of QGIS or ArcGIS (hence the preferred QGIS alternative).
## What is Spatial Analysis?
A distinction is made in this course between GIS and spatial analysis.
In the context of mainstream GIS software, the term *analysis*
refers to data manipulation and data querying. In the context
of **spatial analysis**, the *analysis* focuses on the *statistical analysis* of
patterns and underlying processes or more generally, spatial analysis
addresses the question “what could have been the genesis of the
observed spatial pattern?” It’s an **exploratory process** whereby we
attempt to quantify the observed pattern then explore the processes
that may have generated the pattern.
For example, you record the location of each tree in a well defined
study area. You then map the location of each tree (a GIS task).
At this point, you might be inclined to make inferences about the
observed pattern. Are the trees clustered or dispersed? Is the tree
density constant across the study area? Could soil type or slope have
led to the observed pattern? Those are questions that are addressed
in spatial analysis using quantitative and statistical techniques.
```{r f01-ppp, echo=FALSE, fig.cap = "Distribution of Maple trees in a 1,000 x 1,000 ft study area.", fig.width=2, fig.height=2}
library(spatstat)
OP <- par( mar=c(0,0,0,0) )
plot(split(lansing)$maple, pch=16, cex=0.5, cols="#222222",main="")
par(OP)
```
What you will learn in this course is that popular GIS software like ArcGIS are great tools to create and manipulate spatial data, but if one wishes to go beyond the data manipulation and analyze patterns and processes that may have led to these patterns, other quantitative tools are needed. One such tool we will use in this class is **R**: an open source (freeware) data analysis environment.
R has one, if not the *richest* set of spatial data analysis and statistics tools available today. Learning the R programming environment will prove to be quite beneficial given that many of the operations learnt are transferable across many other (non-spatial) quantitative analysis projects.
[R](http://www.r-project.org/) can be installed on both Windows and Mac operating systems. Another related piece of software that you might find useful is [RStudio](https://www.rstudio.com/products/rstudio/download/) which offers a nice interface to R. To learn more about data analysis in R, visit the [ES218 course website](http://mgimond.github.io/ES218/).
## What's in an Acronym?
GIS is a ubiquitous technology. Many of you are taking this course in part because you have seen GIS listed as a "desirable"" or "required" skill in job postings. Many of you will think of GIS as a “map making” environment as do many ancillary users of GIS in the workforce. While "visualizing" data is an important feature of a GIS, one must not lose sight of *what* data is being visualized and for what purpose. O'Sullivan and Unwin [@Unwin1] use the term **accidental geographer** to refer to those *"whose understanding of geographic science is based on the operations made possible by GIS software"*. We can expand on this idea and define **accidental data analyst** as one whose understanding of data and its analysis is limited to the point-and-click environment of popular pieces of software such as spreadsheet environments, statistical packages and GIS software. The aggressive marketing of GIS technology has the undesirable effect of placing the *technology* before *purpose* and *theory*. This is not unique to GIS, however. Such concerns were shared decades ago when personal computers made it easier for researchers and employees to graph non-spatial data as well as perform many statistical procedures.
The different purposes of mapping spatial data have strong parallels to that of graphing (or plotting) non-spatial data. **John Tukey** [@Tukey1972] offers three broad classes of the latter:
* "*Graphs from which numbers are to be read off- substitutes for tables.*
* *Graphs intended to show the reader what has already been learned (by some other technique)--these we shall sometimes impolitely call propaganda graphs.*
* *Graphs intended to let us see what may be happening over and above what we have already described- these are the analytical graphs that are our main topic.*"
A *GIS world* analogy is proposed here:
* **Reference maps** (USGS maps, hiking maps, road maps). Such maps are used to navigate landscapes or identify locations of points-of-interest.
* **Presentation maps** presented in the press such as the NY Times and the Wall Street Journal, but also maps presented in journals. Such maps are designed to convey a very specific narrative of the author's choosing. (Here we'll avoid Tukey's harsh description of such visual displays, but the idea that maps can be used as *propaganda* is not farfetched).
* **Statistical maps** whose purpose it is to manipulate the raw data in such a way to tease out patterns otherwise not discernable in its original form. This usually requires multiple data manipulation operations and visualization and can sometimes benefit from being explored outside of a spatial context.
This course will focus on the last two spatial data visualization purposes with a strong emphasis on the latter (*Statistical maps*).