Generate_U: Generate spatial data that matches population counts

In DanielBonnery/Strategy: X

Description

Generate spatial data that matches population counts

Usage

Generate_U(SpatialData, .id = NULL, .spatialobject, type = "random")

Arguments

SpatialData	: an object of class that includes
type	: argument to be passed to sp::spsample

Examples

data(parish110217popest,package="dataONS")
data("mtcty150217population",package="dataONS")
shapeData2<-dataONS::dataParishes_December_2011_Boundaries_EW_BFC()
yy<-unique(get(data(Output_Area_to_Parish_to_Local_Authority_District_December_2011_Lookup_in_England_and_Wales,package="dataONS"))[c("PAR11CD","LAD11NM")])
names(yy)<-tolower(names(yy))
shapeData<-sp::merge(shapeData2,yy,by="par11cd",duplicateGeoms = TRUE)
parish110217popest2<-parish110217popest[
 is.element(parish110217popest$PAR11CD,
              shapeData$par11cd)&
                  parish110217popest$year=="mid_2006",
                    c("PAR11CD","Population")]
                    names(parish110217popest2)<-tolower(names(parish110217popest2))
                    shapeData=sp::merge(shapeData,parish110217popest2,by="par11cd",duplicateGeoms = TRUE)
                    shapeData$population[is.na(shapeData$population)]<-mean(shapeData$population,na.rm=TRUE)
                    shapeData<-subset(shapeData,is.element(lad11nm ,c("Allerdale", "Barrow-in-Furness", "Carlisle", "Copeland", "Eden","South Lakeland")))

U<-Generate_U(shapeData,.id="par11cd",.spatialobject="st_areasha",type="random")
popbins<-quantile(shapeData$population,(seq_len(11)-1)/10) 
poppal <- colorBin(heat.colors(5), bins=popbins, na.color = "#aaff56",reverse = T)
library(leaflet)

leaflet(U) %>% 
  addPolygons(data=shapeData,
              stroke=TRUE,
              weight=1,
              color="black",
              fillOpacity=5,
              fillColor=~poppal(shapeData$population)) %>% 
  addTiles() %>% 
 addLegend(title = "Population count", pal=poppal, 
             values=shapeData$population,
              opacity=1, 
              na.label = "Not Available")

DanielBonnery/Strategy documentation built on Dec. 17, 2021, 4:03 p.m.