# load data for urban proxies and regions
proxies <- load_proxies_belgium()
regions <- flexurba::units_belgium
# convert Belgian provinces ('GID_2') to a zonal grid
gridded_regions <- convert_regions_to_grid(regions, proxies$pop, "GID_2")
terra::plot(gridded_regions)
Run the code above in your browser using DataLab