Kathmandu COVID-19 Live Map

library(tidyverse)

## -- Attaching packages ----------------------------------------------------------------- tidyverse 1.3.0 --

## v ggplot2 3.3.2     v purrr   0.3.4
## v tibble  3.0.1     v dplyr   1.0.0
## v tidyr   1.1.0     v stringr 1.4.0
## v readr   1.3.1     v forcats 0.5.0

## -- Conflicts -------------------------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

library(sf)

## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1

library(leaflet)
library(viridis)

## Loading required package: viridisLite

getwd()

## [1] "C:/Users/gocoo/OneDrive/Documents"

covid<-read.csv('cases.csv')
setwd("C:/Users/gocoo/OneDrive/Documents/")

neighborhoods <- read_sf("~/shpfile/local_unit.shp")
neighborhoods<-neighborhoods %>% filter(DISTRICT == "KATHMANDU"| DISTRICT == "LALITPUR"|DISTRICT == "BHAKTAPUR")

hud_grants <- read_csv("cases.csv") %>% 
  st_as_sf(coords = c("lon", "lat"), 
           crs = 4326, agr = "field", na.fail=F)

## Warning: Missing column names filled in: 'X1' [1]

## Parsed with column specification:
## cols(
##   X1 = col_double(),
##   lon = col_double(),
##   lat = col_double(),
##   covid_cases = col_double(),
##   label = col_character()
## )

pal <- colorFactor(
  palette = viridis_pal(begin = .95, end = .4, option = 'C')(136),
  domain = covid$covid_cases
)

hud_grants_clean <- hud_grants%>% mutate(popup_label = paste(paste0('<b>Cases: ', covid_cases, '</b>'),
                                                             paste0('Location: ', label), 
                                                             sep = '<br/>'))

leaflet() %>% 
  addTiles() %>% 
  addPolygons(data = neighborhoods,
              color = 'white',
              weight = 1.5,
              opacity = 1,
              fillColor = "black",
              fillOpacity = .7,
              highlightOptions = highlightOptions(color = "#FFF1BE", 
                                                  weight = 5),
              popup = ~DISTRICT)%>% 
addCircleMarkers(data = hud_grants_clean,
                 popup = ~popup_label,
                 stroke = F,
                 radius = 4,
                 fillColor = ~pal(covid_cases),
                 fillOpacity = 1)

## Warning: sf layer has inconsistent datum (+proj=longlat +ellps=evrst30 +no_defs).
## Need '+proj=longlat +datum=WGS84'

## Warning in validateCoords(lng, lat, funcName): Data contains 1 rows with either
## missing or invalid lat/lon values and will be ignored