## Not run:
# 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)
# ## End(Not run)
Run the code above in your browser using DataLab