Convert flow data to SpatialLinesDataFrame

Origin-destination ('OD') flow data is often provided in the form of 1 line per flow with zone codes of origin and destination centroids. This can be tricky to plot and link-up with geographical data. This function makes the task easier.

od2line(flow, zones, destinations = NULL, zone_code = names(zones)[1],
  origin_code = names(flow)[1], dest_code = names(flow)[2],
  zone_code_d = NA, silent = FALSE)

od2line2(flow, zones)

Arguments

flow	A data frame representing the flow between two points or zones. The first two columns of this data frame should correspond to the first column of the data in the zones. Thus in cents(), the first column is geo_code. This corresponds to the first two columns of flow().
zones	A spatial object representing origins (and destinations if no separate destinations object is provided) of travel.
destinations	A SpatialPolygonsDataFrame or SpatialPointsDataFrame representing destinations of travel flows.
zone_code	Name of the variable in zones containing the ids of the zone. By default this is the first column names in the zones.
origin_code	Name of the variable in flow containing the ids of the zone of origin. By default this is the first column name in the flow input dataset.
dest_code	Name of the variable in flow containing the ids of the zone of destination. By default this is the second column name in the flow input dataset or the first column name in the destinations if that is set.
zone_code_d	Name of the variable in destinations containing the ids of the zone. By default this is the first column names in the destinations.
silent	TRUE by default, setting it to TRUE will show you the matching columns

Details

The function expects zone codes to be in the 1st column of the zones/destinations datasets and the 1st and 2nd columns of the flow data, respectively.

od2line2() is a faster implementation (around 6 times faster on large datasets) that returns a SpatialLines object, omitting the data and working only when there is no destinations dataset (i.e. when the geography of origins is the same as that of destinations).

See also

Other od: dist_google, od2odf, od_aggregate_from, od_aggregate_to, od_aggregate, od_coords2line, od_coords, od_dist, od_id, od_oneway, od_radiation, od_to_odmatrix, odmatrix_to_od, points2flow, points2odf, sp_aggregate

Examples

l <- od2line(flow = flow, zones = cents)
plot(cents)
plot(l, lwd = l$All / mean(l$All), add = TRUE)
# When destinations are different
head(flow_dests[1:5]) # check data
#>        Area.of.residence    WZ11CD All Work.mainly.at.or.from.home
#> 920573         E02002361 E33010350 109                           0
#> 920575         E02002361 E33010390  38                           0
#> 920578         E02002361 E33010468  10                           0
#> 920582         E02002361 E33012367  44                           0
#> 920587         E02002361 E33012225  34                           0
#> 920591         E02002361 E33012368   7                           0
#>        Underground..metro..light.rail..tram
#> 920573                                    0
#> 920575                                    0
#> 920578                                    0
#> 920582                                    0
#> 920587                                    0
#> 920591                                    0
head(destinations[1:5])
#> class       : SpatialPointsDataFrame 
#> features    : 6 
#> extent      : -1.540901, -1.492713, 53.80131, 53.81724  (xmin, xmax, ymin, ymax)
#> crs         : +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
#> variables   : 5
#> names       :    WZ11CD,   LAD11CD, COWZEW_SG,                     COWZEW_SGN, COWZEW_G 
#> min values  : E33010351, E08000035,         2, Manufacturing and distribution,       2c 
#> max values  : E33012352, E08000035,         7,                       Top jobs,       7c 
flowlines_dests <- od2line(flow_dests, cents, destinations = destinations)
plot(flowlines_dests)
l <- od2line(flow, zones_sf)
#> Creating centroids representing desire line start and end points.
#> Warning: st_centroid assumes attributes are constant over geometries of x
#> Warning: st_centroid does not give correct centroids for longitude/latitude data
plot(l["All"], lwd = l$All/mean(l$All))