3

I have an Australia shapefile (can be found here http://www.abs.gov.au/AUSSTATS/abs@.nsf/DetailsPage/1259.0.30.001July%202011?OpenDocument) and some custom data which contains latitudes and longitudes for a number of people.

How can I join the shapefile to the custom data using latitude and longitude so that I can create a population distribution map and colour the segments of the map accordingly?

  • Can you give an example for what exactly your 'custom data' table looks like? – cengel Jul 9 '15 at 15:45
3

Provided I understood your question correctly, here is a suggestion. I made up some data for population, so if those are different for you the aggregate function call might need to be adapted.

library(rgdal)
library(sp)

# shapefile
tmpdir <- tempdir()
download.file("http://www.abs.gov.au/ausstats/subscriber.nsf/log?openagent&1259030001_ste11aaust_shape.zip&1259.0.30.001&Data%20Cubes&D39E28B23F39F498CA2578CC00120E25&0&July%202011&14.07.2011&Latest", "ste11aaust.zip")
unzip("ste11aaust.zip", exdir = tmpdir )
ste <- paste0(tmpdir,"/STE11aAust.shp")
AUS <- readOGR(ste, "STE11aAust")
unlink(tmpdir) 

# make up the population data 
set.seed(100)
lon <- runif(10, 1, 147) 
lat <- runif(10, -31, -20)
no_ppl <- sample(1:100, 10)

# create SpatialPointsDataframe with people 
people.spdf <- SpatialPointsDataFrame(cbind(lon, lat), data.frame(no_ppl), proj4string=CRS(proj4string(AUS)))

# find which polygons the points fall into and sum them up
ppl.sum <- aggregate(x = people.spdf["no_ppl"], by = AUS, FUN = sum)

# plot crude choropleth map (note that there are NAs for some states)
spplot(ppl.sum, "no_ppl")

enter image description here

  • +1 for aggregate. Alternately, you could use %over% and tapply(). I also have a point.in.poly() function in the spatialEco package that shortcuts this and adds the polygon ID's and optionally, attributes, to the points. – Jeffrey Evans Jul 9 '15 at 18:35

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.