kidwellj/uk_census_shapes_clipped
1
0
Fork 0
mirror of https://github.com/kidwellj/uk_census_shapes_clipped.git synced 2024-10-31 23:42:19 +00:00
Code Issues Packages Projects Releases Wiki Activity

added postcode_areas

Browse source
This commit is contained in:
Jeremy Kidwell 2022-04-02 11:45:43 +01:00
parent fa803798c4
commit fbccb533a0
2 changed files with 74 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

43
README.md
View file

@ -5,3 +5,46 @@ Choropleth visualisations of data can often end up visually misleading in that t
It is time consuming to produce this data, so I have done so programatically and uploaded the resulting geopackage files to a public repository for easy re-use for reproducible research.
The files in this repository are only the code used to produce those files. Separate licenses apply to that data etc.
Data derived from these computations can be found [on zenodo](https://zenodo.org/deposit/6395804).
Please note, final clipping doesn't seem to work in R on the current version, so I am finishing these using `difference` in QGIS. Output with versions used is as follows:
```
QGIS version: 3.16.15-Hannover
QGIS code revision: e7fdad6431
Qt version: 5.14.2
GDAL version: 3.2.1
GEOS version: 3.9.1-CAPI-1.14.2
PROJ version: Rel. 6.3.2, May 1st, 2020
Processing algorithm…
Algorithm 'Difference' starting…
Input parameters:
{ 'INPUT' : 'infuse_dist_lyr_2011_simplified_100m.gpkg|layername=infuse_dist_lyr_2011_simplified_100m', 'OUTPUT' : 'infuse_dist_lyr_2011_simplified_100m_buildings_simplified.gpkg', 'OVERLAY' : 'infuse_dist_lyr_2011_simplified_100m_buildings_overlay_simplified.gpkg|layername=infuse_dist_lyr_2011_simplified_100m_buildings_overlay_simplified' }
Execution completed in 262.19 seconds
Results:
{'OUTPUT': 'infuse_dist_lyr_2011_simplified_100m_buildings.gpkg'}
Loading resulting layers
Algorithm 'Difference' finished
QGIS version: 3.16.15-Hannover
QGIS code revision: e7fdad6431
Qt version: 5.14.2
GDAL version: 3.2.1
GEOS version: 3.9.1-CAPI-1.14.2
PROJ version: Rel. 6.3.2, May 1st, 2020
Processing algorithm…
Algorithm 'Difference' starting…
Input parameters:
{ 'INPUT' : 'infuse_ctry_2011_simplified_100m.gpkg|layername=infuse_ctry_2011_simplified_100m', 'OUTPUT' : 'infuse_ctry_2011_simplified_100m_buildings_simplified.gpkg', 'OVERLAY' : 'infuse_ctry_2011_simplified_100m_buildings_overlay_simplified.gpkg|layername=infuse_ctry_2011_simplified_100_buildings_overlay_simplified' }
Execution completed in 1548.94 seconds
Results:
{'OUTPUT': 'infuse_ctry_2011_simplified_100m_buildings_simplified.gpkg'}
Loading resulting layers
Algorithm 'Difference' finished
```

32
script.R
View file

@ -7,32 +7,43 @@ oa11 <- st_read("./data/infuse_oa_lyr_2011.shp")
lsoa11 <- st_read("./data/infuse_lsoa_lyr_2011.shp")
# download local authorities data for whole UK - level 1 (least detailed)
local_authorities <- st_read("data/infuse_dist_lyr_2011.shp")
# Load buildings shapefile:
buildings <- st_read("./data/OS_Open_Zoomstack_district_buildings.gpkg")
# open country shapefiles
countries <- st_read("data/infuse_ctry_2011.shp")
# load postcode areas shapefile
postcode_areas <- st_read("data/GB_Postcodes/PostalArea.shp")
# load postcode districts shapefile
postcode_districts <- st_read("data/GB_Postcodes/PostalDistrict.shp")
# load postcode areas shapefile
postcode_sectors <- st_read("data/GB_Postcodes/PostalSector.shp")
# Load buildings shapefile:
buildings <- st_read("./data/OS_Open_Zoomstack_district_buildings.gpkg")
# Generate simplified versions at 10m:
oa11_simplified_10m <- st_simplify(oa11, dTolerance = 10) # 10 m
lsoa11_simplified_10m <- st_simplify(lsoa11, dTolerance = 10) # 10 m
local_authorities_simplified_10m <- st_simplify(local_authorities, dTolerance = 10) # 10 m
countries_simplified_10m <- st_simplify(countries, dTolerance = 10) # 10 m
postcode_areas_simplified_10m <- st_simplify(postcode_areas, dTolerance = 10) # 10 m
st_write(oa11_simplified_10m, "data/infuse_oa_lyr_2011_simplified_10m.gpkg")
st_write(lsoa11_simplified_10m, "data/infuse_lsoa_lyr_2011_simplified_10m.gpkg")
st_write(local_authorities_simplified_10m, "data/infuse_dist_lyr_2011_simplified_10m.gpkg")
st_write(countries_simplified_10m, "data/infuse_ctry_2011_simplified_10m.gpkg")
st_write(postcode_areas_simplified_10m, "data/PostalArea_simplified_10m.gpkg")
# Generate simplified versions at 100m:
oa11_simplified_100m <- st_simplify(oa11, dTolerance = 100) # 100 m
lsoa11_simplified_100m <- st_simplify(lsoa11, dTolerance = 100) # 100 m
local_authorities_simplified_100m <- st_simplify(local_authorities, dTolerance = 100) # 100 m
countries_simplified_100m <- st_simplify(countries, dTolerance = 100) # 100 m
postcode_areas_simplified_100m <- st_simplify(postcode_areas, dTolerance = 100) # 10 m
st_write(oa11_simplified_100m, "data/infuse_oa_lyr_2011_simplified_100m.gpkg")
st_write(lsoa11_simplified_100m, "data/infuse_lsoa_lyr_2011_simplified_100m.gpkg")
st_write(local_authorities_simplified_100m, "data/infuse_dist_lyr_2011_simplified_100m.gpkg")
st_write(countries_simplified_100m, "data/infuse_ctry_2011_simplified_100m.gpkg")
st_write(postcode_areas_simplified_100m, "data/PostalArea_simplified_100m.gpkg")
## Initial work on data with buildings shapefile for large scale visualation -----
# Create modest buffer around quite small building polygons for the sake of visualisation at
@ -89,6 +100,15 @@ st_write(difference_oa, "data/infuse_oa_lyr_2011_simplified_100m_buildings_overl
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 postcode areas simplified shapefile -----
difference_postcode_areas <- st_difference(postcode_areas_simplified_100m$geom, buildings_buffer)
## Write results to a file:
st_write(difference_postcode_areas, "data/PostalArea_simplified_100m_buildings_overlay.gpkg")
# Simplify layer
difference_postcode_areas_simplified = st_simplify(difference_postcode_areas, dTolerance = 100) # 100 m
st_write(difference_postcode_areas_simplified, "data/PostalArea_simplified_100m_buildings_overlay_simplified.gpkg")
## Process countries shapefile -----
difference_ctry <- st_difference(countries$geom, buildings_buffer)
# Write results to a file:
@ -123,3 +143,11 @@ difference_oa <- st_difference(oa11$geom, buildings_buffer)
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")
## Process postcode areas shapefile -----
difference_postcode_areas <- st_difference(postcode_areas$geom, buildings_buffer)
## Write results to a file:
st_write(difference_postcode_areas, "data/PostalArea_buildings_overlay.gpkg")
# Simplify layer
difference_postcode_areas_simplified = st_simplify(difference_postcode_areas, dTolerance = 100) # 100 m
st_write(difference_postcode_areas_simplified, "data/PostalArea_buildings_overlay_simplified.gpkg")