Skip to content

Instantly share code, notes, and snippets.

@cavedave
Created April 23, 2017 21:25
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save cavedave/105469579bd611f59567ee4445a37522 to your computer and use it in GitHub Desktop.
Save cavedave/105469579bd611f59567ee4445a37522 to your computer and use it in GitHub Desktop.
read_csv('GEOSTAT_grid_POP_1K_2011_V2_0_1.csv') %>%
rbind(read_csv('JRC-GHSL_AIT-grid-POP_1K_2011.csv') %>%
mutate(TOT_P_CON_DT='')) %>%
mutate(lat = as.numeric(gsub('.*N([0-9]+)[EW].*', '\\1', GRD_ID))/100,
lng = as.numeric(gsub('.*[EW]([0-9]+)', '\\1', GRD_ID)) * ifelse(gsub('.*([EW]).*', '\\1', GRD_ID) == 'W', -1, 1) / 100) %>%
filter(lng > 28, lng < 34) %>%
filter(lat > 28, lat < 34) %>%
group_by(lng=round(lng, 1), lat=round(lat, 1)) %>%
summarize(value = sum(TOT_P, na.rm=TRUE)) %>%
ungroup() %>%
group_by(lat) %>%
arrange(lng) %>%
mutate(next.value = lead(value)) %>%
ggplot(aes(lng, lat + 5*(value/max(value)), xend=lng+0.1, yend=lat + 5*(next.value/max(value)))) +
geom_segment(size=0.4, alpha=0.8, color='#5A3E37') +
ggthemes::theme_map() +
coord_equal(0.9)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment