zones <- rgdal::readOGR(dsn = "/home/robin/npct/pct-bigdata/msoas.geojson", layer = "OGRGeoJSON")
proj4string(zones) <- proj4string(cents)
zones <- zones[cents,]
plot(zones)
points(cents)
# use_data(zones, overwrite = TRUE)
Run the code above in your browser using DataLab