forked from r-spatial/stars
-
Notifications
You must be signed in to change notification settings - Fork 0
/
sp.R
26 lines (22 loc) · 774 Bytes
/
sp.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
options(rgdal_show_exportToProj4_warnings = "none")
library(sp)
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)