Home > Mobile >  Creating multiple isochrones with openrouteservice-r for a dataframe of locations
Creating multiple isochrones with openrouteservice-r for a dataframe of locations

Time:04-15

What I want to do:

I have an sf object with a number of point locations. For each location, I want to generate isochrones with different parameters, using the openrouteservice-r package.

The object looks like so:

library(sf)
library(dplyr)
library(openrouteservice)
dd <- data.frame(x = c(-89.33, -89.33, -89.33, -89.33, -89.55, -89.55, -89.55, -89.55),
                 y = c(43.12, 43.12, 43.12, 43.12, 43.12, 43.12, 43.12, 43.12),
                 bike_type = c("bike", "bike", "e-bike", "e-bike", "bike", "bike", "e-bike", "e-bike"),
                 time = c(300, 600, 300, 600, 300, 600, 300, 600))

df <- st_as_sf(dd, coords = c("x","y"))
df

I created a helper function that takes three parameters, runs them through ors_isochrones and returns the geometry (a polygon) of the isochrone:

# requires API key: https://openrouteservice.org/dev/#/login
# ors_api_key(key, service = "openrouteservice")

create_iso <- function(location, profile_1 = "e-bike", range_1){
  iso <- ors_isochrones(locations = st_coordinates(location), 
                                 profile = ors_profile(mode = profile_1),
                                 range = range_1,
                       output = "sf")
  iso$geometry
}

My thought was that now I could call the function from within mutate to the locations dataframe:

df  %>% 
  mutate(iso = create_iso(geometry, profile, time))

What works:

  1. Doing this for a single row of df:
> df[1,]  %>% 
    mutate(iso = create_iso(geometry, bike_type, time))
[1] 300
Simple feature collection with 1 feature and 2 fields
Active geometry column: geometry
Geometry type: POINT
Dimension:     XY
Bounding box:  xmin: -89.33044 ymin: 43.12018 xmax: -89.33044 ymax: 43.12018
CRS:           NA
  bike_type time                   geometry                            iso
1      bike  300 POINT (-89.33044 43.12018) POLYGON ((-89.34975 43.1200...
  1. Doing this for multiple rows of df with keeping the two other parameters fixed:
df[4:7,]  %>%
  mutate(iso = create_iso(geometry, "e-bike", 100))

What doesn't work

Doing it for multiple rows while providing all three parameters

df[4:7,]  %>% 
    mutate(iso = create_iso(geometry, bike_type, time))
[1] 600 300 600 300
Error in `stopifnot()`:
! Problem while computing `iso = create_iso(geometry, bike_type,
  time)`.
Caused by error in `match.arg()`:
! 'arg' must be of length 1
Run `rlang::last_error()` to see where the error occurred.
Called from: signal_abort(cnd, .file)

Full reprex

library(sf)
#> Warning: package 'sf' was built under R version 4.1.3
#> Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1; sf_use_s2() is TRUE
library(dplyr)
#> Warning: package 'dplyr' was built under R version 4.1.3
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union
library(openrouteservice)
#> Registered S3 method overwritten by 'geojsonlint':
#>   method         from 
#>   print.location dplyr
dd <- data.frame(x = c(-89.33044, -89.33044, -89.33044, -89.33044, -89.55, -89.55, -89.55, -89.55),
                 y = c(43.12018, 43.12018, 43.12018, 43.12018, 43.12019, 43.12019, 43.12019, 43.12019),
                 bike_type = c("bike", "bike", "e-bike", "e-bike", "bike", "bike", "e-bike", "e-bike"),
                 time = c(300, 600, 300, 600, 300, 600, 300, 600))

df <- st_as_sf(dd, coords = c("x","y"))
df
#> Simple feature collection with 8 features and 2 fields
#> Geometry type: POINT
#> Dimension:     XY
#> Bounding box:  xmin: -89.55 ymin: 43.12018 xmax: -89.33044 ymax: 43.12019
#> CRS:           NA
#>   bike_type time                   geometry
#> 1      bike  300 POINT (-89.33044 43.12018)
#> 2      bike  600 POINT (-89.33044 43.12018)
#> 3    e-bike  300 POINT (-89.33044 43.12018)
#> 4    e-bike  600 POINT (-89.33044 43.12018)
#> 5      bike  300    POINT (-89.55 43.12019)
#> 6      bike  600    POINT (-89.55 43.12019)
#> 7    e-bike  300    POINT (-89.55 43.12019)
#> 8    e-bike  600    POINT (-89.55 43.12019)

create_iso <- function(location, profile_1 = "e-bike", range_1){
  print(range_1)
  iso <- ors_isochrones(locations = st_coordinates(location), 
                                 profile = ors_profile(mode = profile_1),
                                 range = range_1,
                       output = "sf")
  iso$geometry
}

x <- df[4:7,]  %>% 
  mutate(iso = create_iso(geometry, bike_type, time))
#> [1] 600 300 600 300
#> Error in `stopifnot()`:
#> ! Problem while computing `iso = create_iso(geometry, bike_type, time)`.
#> Caused by error in `match.arg()`:
#> ! 'arg' must be of length 1

Created on 2022-04-14 by the reprex package (v2.0.1)

CodePudding user response:

To apply your custom function to each row of your data.frame, you could use rowwise():

library(dplyr)

df %>% 
  rowwise() %>% 
  mutate(iso = create_iso(geometry, bike_type, time))
  • Related