library(leaflet)
library(sp)
## Warning: package 'sp' was built under R version 3.3.2
library(rgdal)
## rgdal: version: 1.1-10, (SVN revision 622)
##  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/Vinod Basalalli/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/Vinod Basalalli/Documents/R/win-library/3.3/rgdal/proj
##  Linking to sp version: 1.2-3
library(maptools)
## Warning: package 'maptools' was built under R version 3.3.2
## Checking rgeos availability: TRUE
library(ggmap)
## Warning: package 'ggmap' was built under R version 3.3.3
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 3.3.2
library(HistData)
## Warning: package 'HistData' was built under R version 3.3.3
library(ggthemes)
## Warning: package 'ggthemes' was built under R version 3.3.2
data(Snow.deaths)
data(Snow.pumps)
data(Snow.streets)
data(Snow.polygons)
data(Snow.deaths)

london_map <- get_map(c(-.137,51.513), zoom=17)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=51.513,-0.137&zoom=17&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false
london <- ggmap(london_map)

Denoting the John Snow’s Pump Polygon

SnowMap()

Adding Pump polygons to the John Snow Data

SnowMap(polygons=TRUE, main="John Snow's Cholera Map with Pump Polygons")

Plot for showing density estimate with zoom

SnowMap(xlim=c(6.5,17.5), ylim=c(6,17), polygons=TRUE, density=TRUE,
        main="Annotated John Snow's Cholera Map")

Mapping of Deaths due to Cholera Outbreak & Pump Area in London

setwd("C:/Users/Vinod Basalalli/Documents/Vinod/R files/Spring 2017/AnlyII/Problem Set Mapping/SnowGIS")
deaths <- readShapePoints("Cholera_Deaths")
## Warning: use rgdal::readOGR or sf::st_read
df_deaths <- data.frame(deaths@coords)

coordinates(df_deaths)=~coords.x1+coords.x2

proj4string(df_deaths)=CRS("+init=epsg:27700") 
df_deaths = spTransform(df_deaths,CRS("+proj=longlat +datum=WGS84"))
df=data.frame(df_deaths@coords)


df1 <- data.frame(df$coords.x1, df$coords.x2)

Pump_area <- readShapePoints("Pumps")
## Warning: use rgdal::readOGR or sf::st_read
pump_area <- data.frame(Pump_area@coords)

coordinates(pump_area)=~coords.x1+coords.x2

proj4string(pump_area)=CRS("+init=epsg:27700") 
pump_area = spTransform(pump_area,CRS("+proj=longlat +datum=WGS84"))
pa=data.frame(pump_area@coords)


pa1 <- data.frame(pa$coords.x1, pa$coords.x2)

df2 <- cbind(Parameter = "Deaths", df1)
colnames(df2) <- c("Parameter", "coords.x1", "coords.x2")

pa2 <- cbind(Parameter = "Pump", pa1)
colnames(pa2) <- c("Parameter", "coords.x1", "coords.x2")
tot <- rbind(df2, pa2)

london + geom_point(mapping=aes(x=coords.x1, y=coords.x2, col=Parameter,size=Parameter , shape=Parameter), data=tot) + scale_size_discrete(range = c(2, 4)) + labs(title = " Mapping of Deaths due to Cholera Outbreak & Pump Area in London")
## Warning: Using size for a discrete variable is not advised.
## Warning: Removed 11 rows containing missing values (geom_point).