dlvr.it
Access door-to-door: An intercity efficiency and distributional analysis of the costs of travel by plane, train, and automobile
Recently published:
*
Yu L., Li, M., Dai, Z., Cui, M., and Levinson, D. (2026) Access door-to-door: An intercity efficiency and distributional analysis of the costs of travel by plane, train, and automobile. Transport Policy. Volume 184, August 2026, 104188 [doi]
Existing studies typically evaluate air travel accessibility by examining either air network performance or ground access to airports in isolation. This paper offers a complementary perspective by assessing the national air travel accessibility through a door-to-door framework and comparing it against multiple modes of intercity transport. We compare four scenarios: air-only, railway-only, highway-only, and an optimal-mode scenario. The first three rely exclusively on a single mode for intercity trips, whereas the optimal-mode scenario selects the lowest-cost option among air, rail, and direct driving for each origin–destination pair. The results show that air travel provides higher accessibility and more balanced spatial equity than rail or highway travel at higher cost thresholds. Air travel also delivers clear advantages in regions with significant geographical constraints, where land-based transport infrastructure is limited. Although the optimal-mode scenario generally enhances spatial equity, it reduces within-group equity in regions characterized either by highly developed urban cores (e.g., the Yangtze River Delta in East China) or by significant geographic constraints (e.g., the peninsula areas of Northeast China).
Figure 9: Difference in accessibility among three scenarios (air-only, railway-only, and highway-only) based on the cost-weighted accessibility metrics.
![use anyhow::bail;
use arrow_extendr::FromArrowRobj;
use extendr_api::prelude::*;
use geo_traits::{CoordTrait, PointTrait};
use geoarrow::array::{GeoArrowArrayAccessor, PointArray};
#[extendr]
fn geoarrow_centroid(x: Robj) -> anyhow::Result<Doubles> {
let res = PointArray::from_arrow_robj(&x)?;
let (sum_x, sum_y, n) = res.iter_values().filter_map(|xi| xi.ok()).fold(
(0.0_f64, 0.0_f64, 0usize),
|(sx, sy, n), v| {
let (x, y) = v.coord().unwrap().x_y();
(sx + x, sy + y, n + 1)
},
);
if n == 0 {
bail!("cannot compute centroid of empty array");
}
let res = [sum_x / n as f64, sum_y / n as f64]
.map(Rfloat::from)
.into_iter()
.collect::<Doubles>();
Ok(res)
}](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreifvhwk7acara7a64jwglprbjgszw5vd7dmfox55qpyhhe4c5i3sry)
![library(wk)
library(geoarrow)
devtools::load_all()
#> ℹ Loading extendrtest
pnts <- wk::xy(rnorm(100, -180, 180), runif(100, -90, 90))
x <- geoarrow::as_geoarrow_array(pnts)
geoarrow_centroid(x)
#> [1] -182.711392 -2.806439
sf::st_as_sfc(pnts) |>
sf::st_combine() |>
sf::st_centroid()
#> Geometry set for 1 feature
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: -182.7114 ymin: -2.806439 xmax: -182.7114 ymax: -2.806439
#> CRS: NA
#> POINT (-182.7114 -2.806439)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreifb7wrdfujkxs356k5dj6peukdmtj2hekllcpbheckr4uhihc4r5a)
