ozmaps: Australia Maps

#Data file links
#http://www.abs.gov.au/AUSSTATS/subscriber.nsf/log?openagent&1270055001_sa2_2016_aust_shape.zip&1270.0.55.001&Data%20Cubes&A09309ACB3FA50B8CA257FED0013D420&0&July%202016&12.07.2016&Latest
#http://www.abs.gov.au/AUSSTATS/subscriber.nsf/log?openagent&1270055001_sa4_2016_aust_shape.zip&1270.0.55.001&Data%20Cubes&C65BC89E549D1CA3CA257FED0013E074&0&July%202016&12.07.2016&Latest
#http://www.abs.gov.au/AUSSTATS/subscriber.nsf/log?openagent&1270055001_ste_2016_aust_shape.zip&1270.0.55.001&Data%20Cubes&65819049BE2EB089CA257FED0013E865&0&July%202016&12.07.2016&Latest
#http://www.aec.gov.au/Electorates/gis/gis_datadownload.htm

# Population data set
#http://www.abs.gov.au/ausstats/subscriber.nsf/log?openagent&32350ds0012_asgs2011_2016.xls&3235.0&Data%20Cubes&98E620F732FBA7A6CA2581870037A23F&0&2016&28.08.2017&Latest


library(tidyverse)

#remane state in population dataset
popData <- readr::read_csv('population.csv') %>%
  mutate(State = case_when(
    State == 'New South Wales' ~ 'NSW',
    State == 'Victoria' ~ 'VIC',
    State == 'Queensland' ~ 'QLD',
    State == 'South Australia' ~ 'SA',
    State == 'Western Australia' ~ 'WA',
    State == 'Tasmania' ~ 'TAS',
    State == 'Northern Territory' ~ 'NT',
    State == 'Australian Capital Territory' ~ 'ACT',
    TRUE ~ as.character(State)))

popData %>%
  group_by(State) %>%
  summarise(pop = sum(`Total Persons`)) -> statePopulation

popData %>%
  group_by(`SA4 code`, State) %>%
  summarise(pop = sum(`Total Persons`)) %>%
  rename(SA4_CODE16 = `SA4 code`) -> sa4Population

popData %>%
  group_by(`SA2 code`, State) %>%
  summarise(pop = sum(`Total Persons`)) %>%
  rename(SA2_CODE16 = `SA2 code`) -> sa2Population


## SA4 
#sa4Shp <- rgdal::readOGR('SA4')
#sa4Small <- rmapshaper::ms_simplify(sa4Shp, keep = 0.02)
# load("~/mapsR/sa4Small.Rda")
# 
# sa4_data <- sa4Small@data
# sa4_data$id <- row.names(sa4_data)
# sa4_data$SA4_CODE16 <- as.integer(as.character(sa4_data$SA4_CODE16))
# sa4_map <- ggplot2::fortify(sa4Small)
# sa4_map$group <- paste("g",sa4_map$group,sep=".")
# sa4_map$piece <- paste("p",sa4_map$piece,sep=".")
# 
# sa4_data %>%
#   select(id, SA4_CODE16, SA4_NAME16, AREASQKM16) %>%
#   rename(name = SA4_NAME16, Area_SqKm =AREASQKM16) %>%
#   right_join(sa4Population) %>%
#   right_join(sa4_map) -> sa4_map
# 
# ggplot(sa4_map) + geom_polygon(aes(long, lat, group = group), colour = 'grey')
# 

## SA2
#sa2Shp <- rgdal::readOGR('SA2')
# load("~/mapsR/sa2Small.Rda")
#sa2Small <- rmapshaper::ms_simplify(sa2Shp, keep = 0.05)
# 
# sa2_data <- sa2Small@data
# sa2_data$id <- row.names(sa2_data)
# sa2_data$SA2_CODE16 <- as.integer(as.character(sa2_data$SA2_5DIG16))
# sa2_map <- ggplot2::fortify(sa2Small)
# sa2_map$group <- paste("g",sa2_map$group,sep=".")
# sa2_map$piece <- paste("p",sa2_map$piece,sep=".")
# 
# sa2_data %>%
#   select(id, SA2_CODE16, SA2_NAME16) %>%
#   rename(name = SA2_NAME16) %>%
#   right_join(sa2Population) %>%
#   right_join(sa2_map)-> sa2_map
# 
# ggplot(sa2_map) + geom_polygon(aes(long, lat, group = group), colour = 'grey')
# 

## Electorates

# electShp <- rgdal::readOGR('elect')
#electSmall <- rmapshaper::ms_simplify(electShp, keep = 0.05)
#save(electSmall, file = "electSmall.Rda")
# load("~/mapsR/electSmall.Rda")
# 
# elect_data <- electSmall@data
# elect_data$id <- row.names(elect_data)
# elect_map <- ggplot2::fortify(electSmall)
# elect_map$group <- paste("g",elect_map$group,sep=".")
# elect_map$piece <- paste("p",elect_map$piece,sep=".")
# 
# # Incorporate population
# elect_data %>%
#   select(id, Elect_div, State, Area_SqKm, Total_Population) %>%
#   rename(name =  Elect_div,
#          pop = Total_Population) %>%
#   right_join(elect_map) -> elect_map
# 
# 
# ggplot(elect_map) + geom_polygon(aes(long, lat, group = group),
#                                  colour = 'grey')


##Australian map
#ausShp <- rgdal::readOGR('Aus')
#ausSmall <- rmapshaper::ms_simplify(ausShp, keep = 0.05)
#load("~/mapsR/ausSmall.Rda")
# 
# aus_data <- ausSmall@data
# aus_data$id <- row.names(aus_data)
# aus_map <- ggplot2::fortify(ausSmall)
# aus_map$group <- paste("g", aus_map$group, sep = ".")
# aus_map$piece <- paste("p", aus_map$piece, sep = ".")
# 
# # Switch states to the codes and attach population data to map
# aus_data %>%
#   mutate(State = case_when(
#     STE_NAME16 == 'New South Wales' ~ 'NSW',
#     STE_NAME16 == 'Victoria' ~ 'VIC',
#     STE_NAME16 == 'Queensland' ~ 'QLD',
#     STE_NAME16 == 'South Australia' ~ 'SA',
#     STE_NAME16 == 'Western Australia' ~ 'WA',
#     STE_NAME16 == 'Tasmania' ~ 'TAS',
#     STE_NAME16 == 'Northern Territory' ~ 'NT',
#     STE_NAME16 == 'Australian Capital Territory' ~ 'ACT',
#     TRUE ~ as.character(STE_NAME16))) %>%
#   select(id, State) %>%
#   right_join(statePopulation) %>%
#   right_join(aus_map) -> aus_map