tests/sp.R

In r-spatial/stars: Spatiotemporal Arrays, Raster and Vector Data Cubes

options(rgdal_show_exportToProj4_warnings = "none")
if (suppressPackageStartupMessages(require(sp, quietly = TRUE))) {
demo(meuse)
# remove +init=xxx
meuse@proj4string = CRS(
"+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs"
)
meuse.grid@proj4string = CRS(
"+proj=sterea +lat_0=52.1561605555556 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs"
)

suppressPackageStartupMessages(library(stars))
(m = st_as_stars(meuse))
(gr = st_as_stars(meuse.grid))

m2 = as(m, "Spatial")
summary(m2)
gr2 = as(gr, "Spatial")
summary(gr2)

all.equal(meuse, m2, check.attributes = FALSE)
fullgrid(meuse.grid) = TRUE
all.equal(meuse.grid, gr2, check.attributes = FALSE)

#all.equal(meuse, m2)
#all.equal(meuse.grid, gr2)
}