Last Mile Problems

I have two-ish sensible routes on bike to/from work, depending on mood/weather/roadworks. Except for the last mile.

Next time I’m in the office I’ll grab the GPS log, process it outside of here, and show how the last mile ¹ takes far too long.

Update: this happened:

ride = read_csv("data.csv") %>% 
  mutate(secs = set_units(secs, "seconds")) %>% 
  mutate(dist = set_units(dist, "meters")) %>% 
  mutate(dist = set_units(dist, "km")) %>% 
  mutate(time = set_units(secs, "minutes")) 

ride %>% 
  ggplot(aes(x=time, y=dist)) + geom_point() +
    geom_smooth(method = "lm") +
    expand_limits(x=60)

dist is distance as the crow flies from the target location. Different paths to work are different lengths, but I want to compare door-to-door time.

The gradient of that line corresponds to a sort of average useful speed - time spent going away from target is penalised, only movement towards the target counts. I kind-of want to throw some components together - compare crow-flying distance to work against my average bus of about 70 mins each way. But I suspect with a very large variance…

This means I can actually blog about this, since I’ve already given away enough info for anyone to work out which office I’m at, but haven’t said where home is.

Between Cycle Superhighway 1² and National Route 66 I can get to Point A and I want to be at Work:

coordinates = tribble(~name, ~lat, ~lon,
                      "point A", 53.79314,-1.55170,
                      "Work", 53.79778,-1.53044
        ) %>% 
  select(name, lon, lat) #cheap fix for coords being lon-lat, not lat-lon as copied from OSM

coordinates %>% 
  st_as_sf(coords = c("lon", "lat"), crs = 4326) %>% 
  tm_shape() +
    tm_dots()

The title for today’s post derives from the crow-flying distance between those points:

coordinates %>% 
  st_as_sf(coords = c("lon", "lat"), crs = 4326) %>% 
  st_distance() %>% 
  knitr::kable()

0.000 [m]	1488.619 [m]
1488.619 [m]	0.000 [m]

It’s a matrix of pairwise distances, but 1489m is close enough to 1 mile to make comparisons to “last mile” in shipping.

I tinkered with the valhallr library last year and still have a local Docker image of West Yorkshire set up for routing. Docker is great - I can tell the image to go to sleep when I’m not playing with it.

After the initial setup, which was basically feeding a bbox of limit coordinates into a config file, waiting about 1 afternoon for it to work out the network, this is by far the easiest roll-your-own routing engine I’ve played with.

Note - 1 way systems mean that point A -> work is a different problem to work -> point A.

route_1 = route(coordinates[1,], coordinates[2,], costing="bicycle") 

route_2 = route(coordinates[2,], coordinates[1,], costing="bicycle") 

map_trip(route_1)

map_trip(route_2)

Ah. Point A is in a slightly awkward spot for the return trip, and valhalla doesn’t believe in dismounting and crossing the road as a pedestrian.

point_b = tribble(
  ~name, ~lat, ~lon,
  "point B",53.79698,-1.53084
) %>% 
  select(name, lon, lat)


route_3 = route(point_b, coordinates[1,], costing="bicycle") 

map_trip(route_3)

That looks much more sensible. However, I suspect that the OSM data that I fed into my valhalla server is way out of date, and the odd road might be shut for roadworks now.

Compare with cyclestreets there’s pretty good agreement. In the return trip it does believe in dismounting to cross awkward roads. I also really like the cycle-focused OSM tiles CS use - cycle paths, pubs and cafes are highlighted.

Both are better than what I tried the first time, which mostly involved getting stuck in the inner Leeds loop, roads/cycleways closed for roadworks, and random disappearing cycle paths.

I’ve played with the Google Maps routing service, but the programmers at Google seem to believe that any road that is bike-legal is fair game.

ORS has some interesting features. For the full route I could mark out a polygon of “avoid this area” around the death road if it tried to take me that way. I might try building their thing another time, or playing with their free API tokens. I doubt I’d go past the free tier.