kidwellj/hacking_religion_textbook
1
0
Fork 0
mirror of https://github.com/kidwellj/hacking_religion_textbook.git synced 2025-01-09 08:52:32 +00:00
Code Issues Packages Projects Releases Wiki Activity

adding osm tooling

Browse source
This commit is contained in:
Jeremy Kidwell 2024-02-29 18:54:23 +00:00
parent d5ca1e3749
commit c4eb678555
1 changed files with 19 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

27
hacking_religion/chapter_3.qmd
View file

@ -27,7 +27,6 @@ A good starting point is to aquire some "adminstrative" data. This is a way of r
library(sf) |> suppressPackageStartupMessages() library(sf) |> suppressPackageStartupMessages()
library(here) |> suppressPackageStartupMessages() library(here) |> suppressPackageStartupMessages()
library(tidyverse) |> suppressPackageStartupMessages() library(tidyverse) |> suppressPackageStartupMessages()
# better video device, more accurate and faster rendering, esp. on macos. Also should enable system fonts for display
library(ragg) |> suppressPackageStartupMessages() library(ragg) |> suppressPackageStartupMessages()
library(tmap) |> suppressPackageStartupMessages() library(tmap) |> suppressPackageStartupMessages()
@ -168,6 +167,9 @@ The sf() library has a host of tools for geospatial data analysis, including the
```{r} ```{r}
local_authorities$churches_count <- lengths(st_covers(local_authorities, os_openmap_pow)) local_authorities$churches_count <- lengths(st_covers(local_authorities, os_openmap_pow))
local_authorities$churches_percent <- prop.table(local_authorities$churches_count) local_authorities$churches_percent <- prop.table(local_authorities$churches_count)
bbox_brum$churches_count <- lengths(st_covers(bbox_brum, os_openmap_pow))
bbox_brum$churches_percent <- prop.table(bbox_brum$churches_count)
``` ```
Now let's visualise this data using tmap, which (now that we have that new column) we can achieve using `tm_fill` specifying the name of our new column: Now let's visualise this data using tmap, which (now that we have that new column) we can achieve using `tm_fill` specifying the name of our new column:
@ -232,7 +234,20 @@ os_openmap_important_buildings <- st_read(here("example_data", "os_openmap_impor
# OSM # OSM
library(osmdata) library(osmdata)
Working off great tutorial here: https://jcoliver.github.io/learn-r/017-open-street-map.html uk_pow <- getbb(place_name = "Birmingham") %>%
opq() %>%
add_osm_feature(key = "amenity", value = "place_of_worship") %>%
osmdata_sf()
length(uk_pow$osm_points)
#Test for counts in OS dataset
bb_birmingham <- getbb(place_name = "Birmingham united kingdom", format_out = "sf_polygon")
bb_birmingham <- st_transform(bb_birmingham, crs = 27700) # CRS for WGS84
lengths(st_covers(bb_birmingham, os_openmap_pow))
# Working off great tutorial here: https://jcoliver.github.io/learn-r/017-open-street-map.html
uk_major <- getbb(place_name = "Birmingham") %>% uk_major <- getbb(place_name = "Birmingham") %>%
opq(timeout = 50) %>% opq(timeout = 50) %>%
@ -259,11 +274,6 @@ street_plot <- street_plot +
street_plot street_plot
uk_pow <- getbb(place_name = "Birmingham") %>%
opq() %>%
add_osm_feature(key = "amenity", value = "place_of_worship") %>%
osmdata_sf()
rest_plot <- street_plot + rest_plot <- street_plot +
geom_sf(data = uk_pow$osm_points, geom_sf(data = uk_pow$osm_points,
inherit.aes = FALSE, inherit.aes = FALSE,
@ -273,8 +283,9 @@ rest_plot <- street_plot +
# Print map # Print map
rest_plot rest_plot
bbox_brum <- getbb(place_name = "Birmingham", format_out = "sf_polygon")
tm_shape(bb) + tm_borders()
``` ```