diff --git a/code/censusdata.R b/code/censusdata.R
index 6180243..1057763 100644
--- a/code/censusdata.R
+++ b/code/censusdata.R
@@ -12,38 +12,41 @@ tyneandwear = northeast %>%
            lad_name =="South Tyneside")
 tm_shape(tyneandwear) + tm_polygons()
 
+northeast = northeast %>% 
+  select(geo_code)
+
 # Do either get_pct or get_pct_lines include intrazonal flows?
 # should we use get_od?
-od = get_pct(region = "north-east", purpose = "commute", geography = "msoa", layer = "l")
-lines = od %>% 
+lines = get_pct(region = "north-east", purpose = "commute", geography = "msoa", layer = "l")
+lines_tyneandwear = lines %>% 
   filter(lad_name1 == "Newcastle upon Tyne" | lad_name1 == "Sunderland" | 
            lad_name1 == "Gateshead" | lad_name1 == "North Tyneside" |  
            lad_name1 =="South Tyneside" | lad_name2 == "Newcastle upon Tyne" | 
            lad_name2 == "Sunderland" | lad_name2 == "Gateshead" | 
            lad_name2 == "North Tyneside" |  lad_name2 =="South Tyneside"
          )
-lines = lines[, c(2:length(lines), 1)]
-tm_shape(lines) + tm_lines()
+tm_shape(lines_tyneandwear) + tm_lines()
 
-lines_foot = lines %>% 
+lines_foot = lines_tyneandwear %>% 
   filter(foot > 0) %>% 
-  select(-c(car_driver:ebike_sico2), -bicycle)
-lines_foot = lines_foot %>% 
+  # select(-c(car_driver:ebike_sico2), -bicycle) %>% 
   select(geo_code1, geo_code2, foot)
 
-lines_bicycle = lines %>% 
+lines_bicycle = lines_tyneandwear %>% 
   filter(bicycle > 0) %>% 
-  select(-c(foot:ebike_sico2))
+  # select(-c(foot:ebike_sico2)) %>% 
+  select(geo_code1, geo_code2, bicycle)
 
-lines_car = lines %>% 
+lines_car = lines_tyneandwear %>% 
   filter(car_driver > 0) %>% 
-  select(-bicycle, -foot, -c(govtarget_slc:ebike_sico2))
+  # select(-bicycle, -foot, -c(govtarget_slc:ebike_sico2)) %>% 
+  select(geo_code1, geo_code2, car_driver)
 
 
 # Unjittered routes ------------------------------------------------------
 
 library(stplanr)
-# coords = od_coords(lines)
+# coords = od_coords(lines_tyneandwear)
 # from = coords[, 1:2]
 # to = coords[, 3:4]
 # routes_osrm = route_osrm(from = from, to = to, osrm.profile = "foot") # not working]
@@ -78,8 +81,7 @@ tm_shape(car_osrm) + tm_lines("car_driver")
 min_distance_meters = 500
 disag_threshold = 50
 
-northeast = northeast %>% 
-  select(geo_code)
+osm_foot = readRDS("data/osm_foot_2022-12-07.Rds")
 
 # Why do we get this error?
 # Error in UseMethod("st_write") : 
@@ -88,7 +90,7 @@ od_foot_jittered = odjitter::jitter(
   od = lines_foot,
   zones = northeast,
   # zone_name_key = "geo_code",
-  subpoints = osm_drive,
+  subpoints = osm_foot,
   disaggregation_threshold = disag_threshold,
   disaggregation_key = "foot",
   min_distance_meters = min_distance_meters
@@ -112,6 +114,7 @@ od_car_jittered = odjitter::jitter(
   od = lines_car,
   zones = northeast,
   zone_name_key = "geo_code",
+  subpoints = osm_drive,
   disaggregation_threshold = disag_threshold,
   disaggregation_key = "car_driver",
   min_distance_meters = min_distance_meters
diff --git a/code/get-osm.R b/code/get-osm.R
index 21f303c..97c3ab7 100644
--- a/code/get-osm.R
+++ b/code/get-osm.R
@@ -1,4 +1,4 @@
-
+library(tidyverse)
 # Cycling -----------------------------------------------------------------
 
 # Download osm data
@@ -123,3 +123,58 @@ dim(osm_highways)
 
 # tm_shape(osm_highways) + tm_lines()
 saveRDS(osm_highways, "data/osm_drive_2022-12-07.Rds") # 100 MB file
+
+# Walking -----------------------------------------------------------------
+
+# Download osm data
+et = c(
+  "maxspeed",
+  "lanes",
+  "sidewalk",
+  "ref",
+  "surface",
+  "segregated",
+  "state"
+)
+
+# Read-in road network data
+# box = sf::st_bbox(c(xmin = -1.2, xmax = 2.22, ymin = 54.73, ymax = 55.2))
+osm_lines1 = osmextract::oe_get_network(
+  place = "Tyne and Wear",
+  mode = "walking",
+  extra_tags = et
+  # , force_download = TRUE # keep it up-to date
+)
+osm_lines2 = osmextract::oe_get_network(
+  place = "Northumberland",
+  mode = "walking",
+  extra_tags = et
+  # , force_download = TRUE # keep it up-to date
+)
+osm_lines3 = osmextract::oe_get_network(
+  place = "Durham",
+  mode = "walking",
+  extra_tags = et
+  # , force_download = TRUE # keep it up-to date
+)
+osm_lines = bind_rows(osm_lines1, osm_lines2, osm_lines3)
+
+# Simple exploratory analysis
+osm_lines %>%
+  filter(highway == "track") %>%
+  sample_n(200) %>% 
+  qtm()
+osm_lines %>%
+  filter(highway == "trunk") %>%
+  sample_n(500) %>% 
+  qtm()
+
+to_exclude = "motorway|trunk|services|disused|emergency|escap|far|rest"
+
+osm_highways = osm_lines %>%
+  filter(!str_detect(string = highway, pattern = to_exclude))
+
+dim(osm_highways) 
+# [1] 210700     20
+saveRDS(osm_highways, "data/osm_foot_2022-12-07.Rds") # 100 MB file
+