Convert output from CycleStreets.net into sf object
json2sf_cs( obj, cols = NULL, cols_extra = c("elevation_start", "elevation_end", "gradient_segment", "elevation_change", "provisionName", "quietness_segment"), smooth_gradient = FALSE, distance_cutoff = 50, gradient_cutoff = 0.1, n = 3 )
Arguments
| obj | Object from CycleStreets.net read-in with |
|---|---|
| cols | Columns to be included in the result, a character vector or |
| cols_extra | Additional columns to be added providing summaries of gradient and other variables |
| smooth_gradient | Identify and fix anomalous gradients? TRUE by default. See https://github.com/Robinlovelace/cyclestreets/issues/14 |
| distance_cutoff | Distance (m) used to identify anomalous gradients |
| gradient_cutoff | Gradient (%, e.g. 0.1 being 10%) used to identify anomalous gradients |
| n | The number of segments to use to smooth anomalous gradents. The default is 3, meaning segments directly before, after and including the offending segment. |
Examples
from = "Leeds Rail Station" to = "University of Leeds" # from_point = tmaptools::geocode_OSM(from) # to_point = tmaptools::geocode_OSM(to) from_point = c(-1.54408, 53.79360) to_point = c(-1.54802, 53.79618) # save result from the API call to journey.json # res_json = journey(from_point, to_point, silent = FALSE, save_raw = TRUE) # jsonlite::write_json(res_json, "inst/extdata/journey.json") f = system.file(package = "cyclestreets", "extdata/journey.json") obj = jsonlite::read_json(f, simplifyVector = TRUE) rsf = json2sf_cs(obj, cols = c("distances")) names(rsf)#> [1] "distances" "elevation_start" "elevation_end" #> [4] "gradient_segment" "elevation_change" "provisionName" #> [7] "quietness_segment" "geometry"rsf#> Simple feature collection with 5 features and 7 fields #> Geometry type: LINESTRING #> Dimension: XY #> Bounding box: xmin: -1.54802 ymin: 53.7936 xmax: -1.54324 ymax: 53.79618 #> Geodetic CRS: WGS 84 #> distances elevation_start elevation_end gradient_segment elevation_change #> 1 87 27 28 0.01149425 1 #> 2 234 28 33 0.02564103 6 #> 3 18 33 32 0.05555556 1 #> 4 131 32 34 0.01526718 2 #> 5 36 34 34 0.00000000 0 #> provisionName quietness_segment geometry #> 1 Unclassified road 0.2478632 LINESTRING (-1.54408 53.793... #> 2 Main road 0.1800000 LINESTRING (-1.54324 53.794... #> 3 Main road 0.5294118 LINESTRING (-1.54624 53.795... #> 4 Main road 0.2330961 LINESTRING (-1.54649 53.795... #> 5 Main road 0.2000000 LINESTRING (-1.54747 53.796...#> [1] "time" "busynance" "signalledJunctions" #> [4] "signalledCrossings" "name" "walk" #> [7] "elevations" "distances" "elevation_mean" #> [10] "elevation_start" "elevation_end" "elevation_max" #> [13] "elevation_min" "gradient_mean" "gradient_median" #> [16] "gradient_p75" "gradient_max" "gradient_segment" #> [19] "elevation_change" "provisionName" "quietness_segment" #> [22] "start" "finish" "startSpeed" #> [25] "start_longitude" "start_latitude" "finish_longitude" #> [28] "finish_latitude" "crow_fly_distance" "event" #> [31] "whence" "speed" "itinerary" #> [34] "clientRouteId" "plan" "note" #> [37] "length" "quietness" "west" #> [40] "south" "east" "north" #> [43] "leaving" "arriving" "grammesCO2saved" #> [46] "calories" "edition" "flow" #> [49] "geometry"#> Simple feature collection with 5 features and 9 fields #> Geometry type: LINESTRING #> Dimension: XY #> Bounding box: xmin: -1.54802 ymin: 53.7936 xmax: -1.54324 ymax: 53.79618 #> Geodetic CRS: WGS 84 #> time busynance elevations elevation_start elevation_end gradient_segment #> 1 24 351 NA 27 28 0.01149425 #> 2 105 1300 NA 28 33 0.02564103 #> 3 5 34 NA 33 32 0.05555556 #> 4 48 562 NA 32 34 0.01526718 #> 5 9 180 NA 34 34 0.00000000 #> elevation_change provisionName quietness_segment #> 1 1 Unclassified road 0.2478632 #> 2 6 Main road 0.1800000 #> 3 1 Main road 0.5294118 #> 4 2 Main road 0.2330961 #> 5 0 Main road 0.2000000 #> geometry #> 1 LINESTRING (-1.54408 53.793... #> 2 LINESTRING (-1.54324 53.794... #> 3 LINESTRING (-1.54624 53.795... #> 4 LINESTRING (-1.54649 53.795... #> 5 LINESTRING (-1.54747 53.796...json2sf_cs(obj, cols = c("distances"), smooth_gradient = TRUE, gradient_cutoff = 0.05, distance_cutoff = 50)#>#> Simple feature collection with 5 features and 8 fields #> Geometry type: LINESTRING #> Dimension: XY #> Bounding box: xmin: -1.54802 ymin: 53.7936 xmax: -1.54324 ymax: 53.79618 #> Geodetic CRS: WGS 84 #> distances elevation_start elevation_end gradient_segment elevation_change #> 1 87 27 28 0.01149425 1 #> 2 234 28 33 0.02564103 6 #> 3 18 33 32 0.05555556 1 #> 4 131 32 34 0.01526718 2 #> 5 36 34 34 0.00000000 0 #> provisionName quietness_segment geometry #> 1 Unclassified road 0.2478632 LINESTRING (-1.54408 53.793... #> 2 Main road 0.1800000 LINESTRING (-1.54324 53.794... #> 3 Main road 0.5294118 LINESTRING (-1.54624 53.795... #> 4 Main road 0.2330961 LINESTRING (-1.54649 53.795... #> 5 Main road 0.2000000 LINESTRING (-1.54747 53.796... #> gradient_smooth #> 1 0.01149425 #> 2 0.02564103 #> 3 0.02349869 #> 4 0.01526718 #> 5 0.00000000
