John Snow Data
library(ggplot2)
library(rgdal)
## Warning: package 'rgdal' was built under R version 3.3.3
## Loading required package: sp
## Warning: package 'sp' was built under R version 3.3.3
## rgdal: version: 1.2-5, (SVN revision 648)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 2.0.1, released 2015/09/15
## Path to GDAL shared files: C:/Users/Kanchana/Documents/R/win-library/3.3/rgdal/gdal
## Loaded PROJ.4 runtime: Rel. 4.9.2, 08 September 2015, [PJ_VERSION: 492]
## Path to PROJ.4 shared files: C:/Users/Kanchana/Documents/R/win-library/3.3/rgdal/proj
## Linking to sp version: 1.2-4
library(maptools)
## Warning: package 'maptools' was built under R version 3.3.3
## Checking rgeos availability: FALSE
## Note: when rgeos is not available, polygon geometry computations in maptools depend on gpclib,
## which has a restricted licence. It is disabled by default;
## to enable gpclib, type gpclibPermit()
library(ggmap)
## Warning: package 'ggmap' was built under R version 3.3.3
library(sp)
library(HistData)
## Warning: package 'HistData' was built under R version 3.3.3
Get_London <- get_map(c(-.137,51.513), zoom=16)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=51.513,-0.137&zoom=16&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false
London <- ggmap(Get_London)
setwd("C:/Users/Kanchana/Documents/SnowGIS_SHP")
deaths <- readShapePoints("Cholera_Deaths")
## Warning: use rgdal::readOGR or sf::st_read
pumps <- readShapePoints("Pumps")
## Warning: use rgdal::readOGR or sf::st_read
tmp.deaths <- data.frame(deaths@coords)
tmp.pumps <- data.frame(pumps@coords)
tmp <- rbind(tmp.deaths, tmp.pumps)
tmp$type <- c(rep('death', times=dim(tmp.deaths)[1]),
rep('pump', times=dim(tmp.pumps)[1]))
coordinates(tmp)=~coords.x1+coords.x2
proj4string(tmp)=CRS("+init=epsg:27700")
tmp = spTransform(tmp, CRS("+proj=longlat +datum=WGS84"))
tmp <- data.frame(tmp@coords, type=tmp@data$type)
##Map of London:
##---------------------------------------
snow.plot <- London
snow.plot
## Plotting John Snow Dataset on London map:
##----------------------------------
London + geom_point(mapping = aes(x=coords.x1, y=coords.x2, col=type),data=tmp)
## Deaths due to Cholera using Density Lines:
##---------------------------------------------
snow.plot2 <- snow.plot + geom_density2d(data = tmp[tmp$type ==
"death", ], aes(x = coords.x1, y = coords.x2),
size = 0.3)
snow.plot2
## Plotting the heat-map for the deaths due to cholera:
##-------------------------------------------------------------------
snow.plot3 <- snow.plot2 + stat_density2d(data = tmp[tmp$type == "death", ], aes(x = coords.x1, y = coords.x2,
fill = ..level.., alpha = ..level..), size = 0.01, bins = 10, geom = "polygon") + scale_fill_gradient(low = "blue",
high = "red", guide = FALSE) + scale_alpha(range = c(0, 0.5), guide = FALSE)
snow.plot3
