od_to_sf: Convert OD data into geographic 'desire line' objects

In od: Manipulate and Map Origin-Destination Data

Convert OD data into geographic 'desire line' objects

Description

Convert OD data into geographic 'desire line' objects

Usage

od_to_sf(
  x,
  z,
  zd = NULL,
  odc = NULL,
  silent = FALSE,
  filter = TRUE,
  package = "sfheaders",
  crs = 4326
)

od_to_sfc(
  x,
  z,
  zd = NULL,
  silent = TRUE,
  package = "sfheaders",
  crs = 4326,
  filter = TRUE
)

Arguments

x	A data frame in which the first two columns are codes representing points/zones of origin and destination
z	Zones representing origins and destinations
zd	Zones representing destinations
odc	A matrix containing coordinates representing line start and end points
silent	Hide messages? FALSE by default.
filter	Remove rows with no matches in z? TRUE by default
package	Which package to use to create the sf object? sfheaders is the default.
crs	The coordinate reference system of the output, if not known in z. 4326 by default.

Examples

x = od_data_df
z = od_data_zones
desire_lines = od_to_sf(x, z)
desire_lines[1:3]
plot(desire_lines)
desire_lines_d = od_to_sf(od_data_df2, od_data_centroids2, od_data_destinations)
o1 = od_data_centroids2[od_data_centroids2[[1]] == od_data_df2[[1]][1], ]
d1 = od_data_destinations[od_data_destinations[[1]] == od_data_df2[[2]][1], ]
plot(desire_lines_d$geometry)
plot(od_data_centroids2$geometry, add = TRUE, col = "green")
plot(od_data_destinations$geometry, add = TRUE)
plot(o1, add = TRUE)
plot(d1, add = TRUE)
plot(desire_lines_d$geometry[1], lwd = 3, add = TRUE)
n = 7
on = od_data_centroids2[od_data_centroids2[[1]] == od_data_df2[[1]][n], ]
dn = od_data_destinations[od_data_destinations[[1]] == od_data_df2[[2]][n], ]
plot(desire_lines_d$geometry)
plot(on, add = TRUE)
plot(dn, add = TRUE)
plot(desire_lines_d$geometry[n], lwd = 3, add = TRUE)

od documentation built on Sept. 11, 2024, 9:04 p.m.