Mapping the corona virus
The outbreak of the corona virus has caused a lot of concern and has raised alarms across the globe.
There is a near real-time reporting being mapped by the John Hopkins university.
The data is being published in a Google Spreadsheet at https://gisanddata.maps.arcgis.com/apps/opsdashboard/index.html#/bda7594740fd40299423467b48e9ecf6
The Google Spreadsheet is: https://docs.google.com/spreadsheets/d/1yZv9w9zRKwrGTaR-YzmAqMefw4wMlaXocejdxZaTs6w/htmlview?usp=sharing&sle=true
Here is a quick map of the confirmed cases.
Blog post to be updated. Made in R.
Intermediate update: Modified thanks to inputs at https://www.linkedin.com/posts/santosh-srinivas-ba2b2919bmapping-the-corona-virus-activity-6628120232328339456-fGkg
setwd("~/Dropbox/pandora/My-Projects/repos/diary/writing/")
# Mapping the corona virus
# read data from googlesheets
# devtools::install_github("tidyverse/googlesheets4")
library(googlesheets4)
# got from google
# https://docs.google.com/spreadsheets/d/1yZv9w9zRKwrGTaR-YzmAqMefw4wMlaXocejdxZaTs6w/htmlview?usp=sharing&sle=true
sheet_id <- "1yZv9w9zRKwrGTaR-YzmAqMefw4wMlaXocejdxZaTs6w"
sheets_url <- "https://docs.google.com/spreadsheets/d/1yZv9w9zRKwrGTaR-YzmAqMefw4wMlaXocejdxZaTs6w/htmlview?usp=sharing&sle=true"
sheets_get(sheet_id)
sheets <- sheets_get(sheets_url)
df <- read_sheet(sheets, sheet = 1)
# I need to map countries to latitude and longtude!!?? How?
# data pasta to rescue
# get the data from https://developers.google.com/public-data/docs/canonical/countries_csv
# read here: https://www.santoshsrinivas.com/data-pasta-is-a-must-have-to-working-with-r-and-the-web/
require('data.table')
# old habits! i love data.table
# pasted below using data pasta mentioned above!
dt_countries <- data.table::data.table(
country = c("AD","AE",
"AF","AG","AI","AL","AM","AN","AO","AQ",
"AR","AS","AT","AU","AW","AZ","BA",
"BB","BD","BE","BF","BG","BH","BI","BJ",
"BM","BN","BO","BR","BS","BT","BV","BW",
"BY","BZ","CA","CC","CD","CF","CG","CH",
"CI","CK","CL","CM","CN","CO","CR","CU",
"CV","CX","CY","CZ","DE","DJ","DK","DM",
"DO","DZ","EC","EE","EG","EH","ER",
"ES","ET","FI","FJ","FK","FM","FO","FR",
"GA","GB","GD","GE","GF","GG","GH","GI",
"GL","GM","GN","GP","GQ","GR","GS","GT",
"GU","GW","GY","GZ","HK","HM","HN","HR",
"HT","HU","ID","IE","IL","IM","IN","IO",
"IQ","IR","IS","IT","JE","JM","JO",
"JP","KE","KG","KH","KI","KM","KN","KP",
"KR","KW","KY","KZ","LA","LB","LC","LI",
"LK","LR","LS","LT","LU","LV","LY","MA",
"MC","MD","ME","MG","MH","MK","ML","MM",
"MN","MO","MP","MQ","MR","MS","MT","MU",
"MV","MW","MX","MY","MZ",NA,"NC","NE",
"NF","NG","NI","NL","NO","NP","NR","NU",
"NZ","OM","PA","PE","PF","PG","PH",
"PK","PL","PM","PN","PR","PS","PT","PW",
"PY","QA","RE","RO","RS","RU","RW","SA",
"SB","SC","SD","SE","SG","SH","SI","SJ",
"SK","SL","SM","SN","SO","SR","ST","SV",
"SY","SZ","TC","TD","TF","TG","TH","TJ",
"TK","TL","TM","TN","TO","TR","TT",
"TV","TW","TZ","UA","UG","UM","US","UY",
"UZ","VA","VC","VE","VG","VI","VN","VU",
"WF","WS","XK","YE","YT","ZA","ZM","ZW"),
latitude = c(42.546245,
23.424076,33.93911,17.060816,18.220554,
41.153332,40.069099,12.226079,-11.202692,
-75.250973,-38.416097,-14.270972,47.516231,
-25.274398,12.52111,40.143105,43.915886,13.193887,
23.684994,50.503887,12.238333,42.733883,
25.930414,-3.373056,9.30769,32.321384,
4.535277,-16.290154,-14.235004,25.03428,27.514162,
-54.423199,-22.328474,53.709807,17.189877,
56.130366,-12.164165,-4.038333,6.611111,
-0.228021,46.818188,7.539989,-21.236736,
-35.675147,7.369722,35.86166,4.570868,9.748917,
21.521757,16.002082,-10.447525,35.126413,
49.817492,51.165691,11.825138,56.26392,
15.414999,18.735693,28.033886,-1.831239,58.595272,
26.820553,24.215527,15.179384,40.463667,
9.145,61.92411,-16.578193,-51.796253,
7.425554,61.892635,46.227638,-0.803689,55.378051,
12.262776,42.315407,3.933889,49.465691,
7.946527,36.137741,71.706936,13.443182,9.945587,
16.995971,1.650801,39.074208,-54.429579,
15.783471,13.444304,11.803749,4.860416,
31.354676,22.396428,-53.08181,15.199999,45.1,
18.971187,47.162494,-0.789275,53.41291,
31.046051,54.236107,20.593684,-6.343194,33.223191,
32.427908,64.963051,41.87194,49.214439,
18.109581,30.585164,36.204824,-0.023559,
41.20438,12.565679,-3.370417,-11.875001,
17.357822,40.339852,35.907757,29.31166,19.513469,
48.019573,19.85627,33.854721,13.909444,
47.166,7.873054,6.428055,-29.609988,55.169438,
49.815273,56.879635,26.3351,31.791702,
43.750298,47.411631,42.708678,-18.766947,7.131474,
41.608635,17.570692,21.913965,46.862496,
22.198745,17.33083,14.641528,21.00789,
16.742498,35.937496,-20.348404,3.202778,
-13.254308,23.634501,4.210484,-18.665695,-22.95764,
-20.904305,17.607789,-29.040835,9.081999,
12.865416,52.132633,60.472024,28.394857,
-0.522778,-19.054445,-40.900557,21.512583,
8.537981,-9.189967,-17.679742,-6.314993,12.879721,
30.375321,51.919438,46.941936,-24.703615,
18.220833,31.952162,39.399872,7.51498,
-23.442503,25.354826,-21.115141,45.943161,
44.016521,61.52401,-1.940278,23.885942,-9.64571,
-4.679574,12.862807,60.128161,1.352083,
-24.143474,46.151241,77.553604,48.669026,
8.460555,43.94236,14.497401,5.152149,3.919305,
0.18636,13.794185,34.802075,-26.522503,
21.694025,15.454166,-49.280366,8.619543,15.870032,
38.861034,-8.967363,-8.874217,38.969719,
33.886917,-21.178986,38.963745,10.691803,
-7.109535,23.69781,-6.369028,48.379433,
1.373333,NA,37.09024,-32.522779,41.377491,
41.902916,12.984305,6.42375,18.420695,18.335765,
14.058324,-15.376706,-13.768752,-13.759029,
42.602636,15.552727,-12.8275,-30.559482,
-13.133897,-19.015438),
longitude = c(1.601554,
53.847818,67.709953,-61.796428,-63.068615,
20.168331,45.038189,-69.060087,17.873887,
-0.071389,-63.616672,-170.132217,14.550072,
133.775136,-69.968338,47.576927,17.679076,
-59.543198,90.356331,4.469936,-1.561593,25.48583,
50.637772,29.918886,2.315834,-64.75737,
114.727669,-63.588653,-51.92528,-77.39628,
90.433601,3.413194,24.684866,27.953389,
-88.49765,-106.346771,96.870956,21.758664,
20.939444,15.827659,8.227512,-5.54708,-159.777671,
-71.542969,12.354722,104.195397,-74.297333,
-83.753428,-77.781167,-24.013197,105.690449,
33.429859,15.472962,10.451526,42.590275,
9.501785,-61.370976,-70.162651,1.659626,
-78.183406,25.013607,30.802498,-12.885834,
39.782334,-3.74922,40.489673,25.748151,179.414413,
-59.523613,150.550812,-6.911806,2.213749,
11.609444,-3.435973,-61.604171,43.356892,
-53.125782,-2.585278,-1.023194,-5.345374,
-42.604303,-15.310139,-9.696645,-62.067641,
10.267895,21.824312,-36.587909,-90.230759,
144.793731,-15.180413,-58.93018,34.308825,
114.109497,73.504158,-86.241905,15.2,-72.285215,
19.503304,113.921327,-8.24389,34.851612,
-4.548056,78.96288,71.876519,43.679291,
53.688046,-19.020835,12.56738,-2.13125,-77.297508,
36.238414,138.252924,37.906193,74.766098,
104.990963,-168.734039,43.872219,-62.782998,
127.510093,127.766922,47.481766,-80.566956,
66.923684,102.495496,35.862285,-60.978893,
9.555373,80.771797,-9.429499,28.233608,
23.881275,6.129583,24.603189,17.228331,-7.09262,
7.412841,28.369885,19.37439,46.869107,
171.184478,21.745275,-3.996166,95.956223,
103.846656,113.543873,145.38469,-61.024174,
-10.940835,-62.187366,14.375416,57.552152,
73.22068,34.301525,-102.552784,101.975766,
35.529562,18.49041,165.618042,8.081666,167.954712,
8.675277,-85.207229,5.291266,8.468946,
84.124008,166.931503,-169.867233,174.885971,
55.923255,-80.782127,-75.015152,-149.406843,
143.95555,121.774017,69.345116,19.145136,
-56.27111,-127.439308,-66.590149,35.233154,
-8.224454,134.58252,-58.443832,51.183884,
55.536384,24.96676,21.005859,105.318756,29.873888,
45.079162,160.156194,55.491977,30.217636,
18.643501,103.819836,-10.030696,14.995463,
23.670272,19.699024,-11.779889,12.457777,
-14.452362,46.199616,-56.027783,6.613081,
-88.89653,38.996815,31.465866,-71.797928,
18.732207,69.348557,0.824782,100.992541,71.276093,
-171.855881,125.727539,59.556278,9.537499,
-175.198242,35.243322,-61.222503,177.64933,
120.960515,34.888822,31.16558,32.290275,NA,
-95.712891,-55.765835,64.585262,12.453389,
-61.287228,-66.58973,-64.639968,-64.896335,
108.277199,166.959158,-177.156097,-172.104629,
20.902977,48.516388,45.166244,22.937506,
27.849332,29.154857),
name = c("Andorra",
"United Arab Emirates","Afghanistan",
"Antigua and Barbuda","Anguilla","Albania","Armenia",
"Netherlands Antilles","Angola",
"Antarctica","Argentina","American Samoa","Austria",
"Australia","Aruba","Azerbaijan",
"Bosnia and Herzegovina","Barbados","Bangladesh",
"Belgium","Burkina Faso","Bulgaria","Bahrain",
"Burundi","Benin","Bermuda","Brunei","Bolivia",
"Brazil","Bahamas","Bhutan","Bouvet Island",
"Botswana","Belarus","Belize","Canada",
"Cocos [Keeling] Islands","Congo [DRC]",
"Central African Republic","Congo [Republic]",
"Switzerland","Côte d'Ivoire","Cook Islands",
"Chile","Cameroon","China","Colombia",
"Costa Rica","Cuba","Cape Verde","Christmas Island",
"Cyprus","Czech Republic","Germany","Djibouti",
"Denmark","Dominica","Dominican Republic",
"Algeria","Ecuador","Estonia","Egypt",
"Western Sahara","Eritrea","Spain","Ethiopia",
"Finland","Fiji",
"Falkland Islands [Islas Malvinas]","Micronesia","Faroe Islands","France",
"Gabon","United Kingdom","Grenada",
"Georgia","French Guiana","Guernsey","Ghana",
"Gibraltar","Greenland","Gambia","Guinea",
"Guadeloupe","Equatorial Guinea","Greece",
"South Georgia and the South Sandwich Islands",
"Guatemala","Guam","Guinea-Bissau","Guyana",
"Gaza Strip","Hong Kong",
"Heard Island and McDonald Islands","Honduras","Croatia","Haiti",
"Hungary","Indonesia","Ireland","Israel",
"Isle of Man","India",
"British Indian Ocean Territory","Iraq","Iran","Iceland","Italy",
"Jersey","Jamaica","Jordan","Japan","Kenya",
"Kyrgyzstan","Cambodia","Kiribati","Comoros",
"Saint Kitts and Nevis","North Korea",
"South Korea","Kuwait","Cayman Islands","Kazakhstan",
"Laos","Lebanon","Saint Lucia",
"Liechtenstein","Sri Lanka","Liberia","Lesotho",
"Lithuania","Luxembourg","Latvia","Libya","Morocco",
"Monaco","Moldova","Montenegro",
"Madagascar","Marshall Islands","Macedonia [FYROM]",
"Mali","Myanmar [Burma]","Mongolia","Macau",
"Northern Mariana Islands","Martinique",
"Mauritania","Montserrat","Malta","Mauritius",
"Maldives","Malawi","Mexico","Malaysia",
"Mozambique","Namibia","New Caledonia","Niger",
"Norfolk Island","Nigeria","Nicaragua",
"Netherlands","Norway","Nepal","Nauru","Niue",
"New Zealand","Oman","Panama","Peru",
"French Polynesia","Papua New Guinea","Philippines",
"Pakistan","Poland",
"Saint Pierre and Miquelon","Pitcairn Islands","Puerto Rico",
"Palestinian Territories","Portugal","Palau",
"Paraguay","Qatar","Réunion","Romania","Serbia",
"Russia","Rwanda","Saudi Arabia",
"Solomon Islands","Seychelles","Sudan","Sweden",
"Singapore","Saint Helena","Slovenia",
"Svalbard and Jan Mayen","Slovakia","Sierra Leone",
"San Marino","Senegal","Somalia","Suriname",
"São Tomé and Príncipe","El Salvador","Syria",
"Swaziland","Turks and Caicos Islands","Chad",
"French Southern Territories","Togo","Thailand",
"Tajikistan","Tokelau","Timor-Leste",
"Turkmenistan","Tunisia","Tonga","Turkey",
"Trinidad and Tobago","Tuvalu","Taiwan","Tanzania",
"Ukraine","Uganda",
"U.S. Minor Outlying Islands","United States","Uruguay","Uzbekistan",
"Vatican City",
"Saint Vincent and the Grenadines","Venezuela","British Virgin Islands",
"U.S. Virgin Islands","Vietnam","Vanuatu",
"Wallis and Futuna","Samoa","Kosovo","Yemen",
"Mayotte","South Africa","Zambia","Zimbabwe")
)
# this can be done using tidyverse
# i'm used to data.table after using for a decade! need to get more comfortable with tidyverse
# lets do the data.table magic for now!
dt <- as.data.table(df)
dt_countries
# need to make it consistent!
dt[,name:=`Country/Region`]
# damn missed this!
dt[name=="Mainland China", name:="China"]
# the merge works!
dt <- merge.data.table(dt, dt_countries, all.x=TRUE)
# none found! works!
dt[is.na(latitude)]
dt[is.na(longitude)]
# install.packages("config")
options(mapbox.accessToken = config::get("mapbox")$accessToken)
library(leaflet)
library(leaflet.mapboxgl)
# credit: https://rpubs.com/mattdray/basic-leaflet-maps
cols <- c("Confirmed", "Deaths", "Recovered")
# https://stackoverflow.com/a/15905564/644081
dt_agg <- dt[,lapply(.SD,sum,na.rm=TRUE),.SDcols=cols, by=list(name, latitude, longitude)]
# naming to make maps work
dt_agg[,lat:=latitude]
dt_agg[,lng:=longitude]
dt_agg <- dt_agg[,list(country=name, confirmed=Confirmed, lat, lng)]
# credit: https://bookdown.org/robinlovelace/geocompr/adv-map.html
# https://rstudio-pubs-static.s3.amazonaws.com/307862_b8c8460272dc4a2a9023d033d5f3ec34.html
pal <- colorNumeric("YlOrRd", domain = dt_agg$confirmed)
# credit: https://rstudio-pubs-static.s3.amazonaws.com/307862_b8c8460272dc4a2a9023d033d5f3ec34.html
dt_agg[,label:=sprintf(
"<strong>%s</strong> confirmed<br/>%s",
confirmed, country
) %>% lapply(htmltools::HTML)]
# Check options here!
# credit: https://bookdown.org/robinlovelace/geocompr/adv-map.html
# setting view to Hong Kong
leaflet(dt_agg) %>%
addMapboxGL(style = "mapbox://styles/mapbox/dark-v10") %>%
addCircleMarkers(radius = ~sqrt(confirmed), col = ~pal(confirmed), opacity = 0.9, , popup = ~label, fillOpacity = 1) %>%
setView(114.109497, 22.396428, zoom = 3)
