kidwellj/uk_census_shapes_clipped
1
0
Fork 0
mirror of https://github.com/kidwellj/uk_census_shapes_clipped.git synced 2025-02-20 15:36:31 +00:00
Code Issues Packages Projects Releases Wiki Activity

added calculations for lsoa and oa

Browse source
This commit is contained in:
Jeremy Kidwell 2022-04-02 11:11:11 +01:00
parent 3180a86ee7
commit fa803798c4
1 changed files with 39 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

48
script.R
View file

@ -48,30 +48,46 @@ st_write(buildings_buffer, "data/OS_Open_Zoomstack_district_buildings_buffered.g
# with filled space for unoccupied layers
difference_ctry <- st_difference(countries_simplified_100m$geom, buildings_buffer)
## Write results to a file:
st_write(difference_ctry, "data/infuse_ctry_2011_simplified_100_buildings_overlay_.gpkg")
st_write(difference_ctry, "data/infuse_ctry_2011_simplified_100m_buildings_overlay.gpkg")
# Simplify layer
difference_ctry_simplified = st_simplify(difference_ctry, dTolerance = 100) # 100 m
st_write(difference_ctry_simplified, "data/infuse_ctry_2011_simplified_100_buildings_overlay_simplified.gpkg")
st_write(difference_ctry_simplified, "data/infuse_ctry_2011_simplified_100m_buildings_overlay_simplified.gpkg")
# Invert polygons in this new layer
inversion_ctry <- st_difference(difference_ctry, countries_simplified_100m$geom)
# Write results to a file:
st_write(inversion_ctry, "data/infuse_ctry_2011_simplified_100_buildings.gpkg")
st_write(inversion_ctry, "data/infuse_ctry_2011_simplified_100m_buildings.gpkg")
inversion_ctry_simplified = st_simplify(inversion_ctry, dTolerance = 100) # 100 m
st_write(inversion_ctry_simplified, "data/infuse_ctry_2011_simplified_100_buildings_simplified.gpkg")
st_write(inversion_ctry_simplified, "data/infuse_ctry_2011_simplified_100m_buildings_simplified.gpkg")
## Process local authorities simplified shapefiles -----
difference_la <- st_difference(local_authorities_simplified_100m$geom, buildings_buffer)
## Write results to a file:
st_write(difference_la, "data/infuse_dist_lyr_2011_simplified_100_buildings_overlay.gpkg")
st_write(difference_la, "data/infuse_dist_lyr_2011_simplified_100m_buildings_overlay.gpkg")
# Simplify layer
difference_simplified = st_simplify(difference_la, dTolerance = 100) # 100 m
st_write(difference_simplified, "data/infuse_dist_lyr_2011_simplified_100_buildings_overlay_simplified.gpkg")
st_write(difference_simplified, "data/infuse_dist_lyr_2011_simplified_100m_buildings_overlay_simplified.gpkg")
# Invert polygons in this new layer
inversion_la <- st_difference(difference_la, local_authorities_simplified_100m$geom)
# Write results to a file:
st_write(inversion_la, "data/infuse_dist_lyr_2011_simplified_100_buildings.gpkg")
st_write(inversion_la, "data/infuse_dist_lyr_2011_simplified_100m_buildings.gpkg")
inversion_simplified_la = st_simplify(inversion_la, dTolerance = 100) # 100 m
st_write(inversion_simplified_la, "data/infuse_dist_lyr_2011_simplified_100_buildings_simplified.gpkg")
st_write(inversion_simplified_la, "data/infuse_dist_lyr_2011_simplified_100m_buildings_simplified.gpkg")
## Process lsoa simplified shapefile -----
difference_lsoa <- st_difference(lsoa11_simplified_100m$geom, buildings_buffer)
## Write results to a file:
st_write(difference_lsoa, "data/infuse_lsoa_lyr_2011_simplified_100m_buildings_overlay.gpkg")
# Simplify layer
difference_lsoa_simplified = st_simplify(difference_lsoa, dTolerance = 100) # 100 m
st_write(difference_lsoa_simplified, "data/infuse_lsoa_lyr_2011_simplified_100m_buildings_overlay_simplified.gpkg")
## Process oa simplified shapefile -----
difference_oa <- st_difference(oa11_simplified_100m$geom, buildings_buffer)
## Write results to a file:
st_write(difference_oa, "data/infuse_oa_lyr_2011_simplified_100m_buildings_overlay.gpkg")
# Simplify layer
difference_oa_simplified = st_simplify(difference_oa, dTolerance = 100) # 100 m
st_write(difference_oa_simplified, "data/infuse_oa_lyr_2011_simplified_100m_buildings_overlay_simplified.gpkg")
## Process countries shapefile -----
difference_ctry <- st_difference(countries$geom, buildings_buffer)
@ -92,4 +108,18 @@ st_write(difference_la, "data/infuse_dist_lyr_2011_buildings_overlay.gpkg")
# Invert polygons in this new layer
inversion_la <- st_difference(difference_la, local_authorities$geom)
# Write results to a file:
st_write(inversion_la, "data/infuse_dist_lyr_2011_buildings.gpkg")
st_write(inversion_la, "data/infuse_dist_lyr_2011_buildings.gpkg")
## Process lsoa shapefile -----
difference_lsoa <- st_difference(lsoa11$geom, buildings_buffer)
# Write results to a file:
st_write(difference_lsoa, "data/infuse_lsoa_lyr_2011_buildings_overlay.gpkg")
difference_lsoa_simplified = st_simplify(difference_lsoa, dTolerance = 100) # 100 m
st_write(difference_lsoa_simplified, "data/infuse_lsoa_lyr_2011_buildings_overlay_simplified.gpkg")
## Process oa shapefile -----
difference_oa <- st_difference(oa11$geom, buildings_buffer)
# Write results to a file:
st_write(difference_oa, "data/infuse_oa_lyr_2011_buildings_overlay.gpkg")
difference_oa_simplified = st_simplify(difference_oa, dTolerance = 100) # 100 m
st_write(difference_oa_simplified, "data/infuse_oa_lyr_2011_buildings_overlay_simplified.gpkg")