# 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)
zones_sf = sf::st_as_sf(zones)
# }
Run the code above in your browser using DataLab