Bug in identicalCRS()?

[Marc Choisy]

Let's consider the example below:

m <- matrix(c(0, 0, 1, 1), ncol = 2, dimnames = list(NULL, c("min", "max")))
crs1 <- CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
crs2 <- CRS("+init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
S1 <- Spatial(bbox = m, proj4string = crs1)
S2 <- Spatial(bbox = m, proj4string = crs2)

Shouldn't the call identicalCRS(S1, S2) return TRUE instead of FALSE?
Thanks for help.

[Edzer Pebesma]

Yes, this can be considered a false negative. sf uses a different mechanism:

> st_crs("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0") == st_crs("+init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
[1] TRUE

If you need it, you may consider backporting that.

[Marc Choisy]

No, no worries, thanks!