Robin Lovelace

Should now be fixed upstream in **spDataLarge**

Being installed on CI, fingers crossed: ``` Run Rscript -e 'remotes::install_github("geocompr/geocompkg", dependencies = TRUE, force = TRUE)' Using github PAT from envvar GITHUB_PAT Downloading GitHub repo geocompr/geocompkg@HEAD spDataLarge (071e59bef... ->...

Still failing, not sure why, and cannot reproduce locally.

One thing for sure: it's not due to spatial subsetting: ``` zone_cents_rail = zone_cents[desire_rail, ] Error in st_geos_binop("intersects", x, y, sparse = sparse, prepared = prepared, : st_crs(x) == st_crs(y)...

Working hypothesis: this is the culprit. ``` bristol_rail_points = rbind( st_sf(data.frame( Node = "Origin and destination locations", col = "black" ), geometry = zone_cents_rail$geometry), st_sf(data.frame( Node = "Public transport node",...

This seems to be the cause: ``` waldo::compare(sf::st_crs(bristol_stations), sf::st_crs(bristol_zones)) lines(old$wkt) vs lines(new$wkt) - "GEOGCRS[\"WGS 84\"," + "GEOGCS[\"WGS 84\"," - " ENSEMBLE[\"World Geodetic System 1984 ensemble\"," + " DATUM[\"WGS_1984\"," - "...

Reprex showing it working locally: ``` r # Aim: find out the cause of the build failing. library(tidyverse) remotes::install_github("nowosad/spDataLarge") #> Skipping install of 'spDataLarge' from a github remote, the SHA1...

Another reprex: ``` r library(spDataLarge) u = "https://github.com/Nowosad/spDataLarge/raw/master/data/bristol_stations.rda" f = basename(u) download.file(u, f) load(f) waldo::compare(sf::st_crs(bristol_zones), sf::st_crs(bristol_stations)) #> lines(old$wkt) vs lines(new$wkt) #> - "GEOGCS[\"WGS 84\"," #> + "GEOGCRS[\"WGS 84\"," #> -...

Latest error: ``` Quitting from lines 526-541 (13-transport.Rmd) Error: arguments have different crs In addition: Warning messages: 1: In st_centroid.sf(zones_od) : st_centroid assumes attributes are constant over geometries of x...

That was fixed, not its... ``` Quitting from lines 691-695 (13-transport.Rmd) Error in geos_op2_geom("difference", x, y, ...) : st_crs(x) == st_crs(y) is not TRUE Calls: local ... st_difference.sf -> geos_op2_df...