[R-sig-Geo] Mixing points and sf objects

Roger Bivand Roger@B|v@nd @end|ng |rom nhh@no
Thu Mar 3 10:32:11 CET 2022

On Thu, 3 Mar 2022, Patrick Giraudoux wrote:

> Dear listers,
> Just for curiosity, using one sf object (depts) and one classical
> data.frame with two columns bearing coordinates, I wonder why this
> sequence works (points and department limits are correctly drawn):
> plot(chbrut[,5:6],asp=1)
> plot(st_geometry(depts),add=TRUE)
> .. but not this one (points of chbrut are not projected - or projected
> somewhere else, and are not visible on the map)
> plot(st_geometry(depts),reset=FALSE)
> points(chbrut[,5:6],pch=19,col=cols)

Sorry, cannot reproduce:

nc <- st_read(system.file("gpkg/nc.gpkg", package="sf"))
nc_xy <- st_transform(nc, "EPSG:32019")
pts <- as.data.frame(st_coordinates(st_sample(nc_xy, 100)))

plot(pts, asp=1)
plot(st_geometry(nc_xy), add=TRUE)



seem the same. reset= is irrelevant if the plot method for sfc objects is 
used, I think it applies to sf objects only, when multiple maps may be 
displayed with keys.

Hope this helps,


> Best,
> Patrick
> 	[[alternative HTML version deleted]]
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo using r-project.org
> https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-geo&data=04%7C01%7CRoger.Bivand%40nhh.no%7C5d59728741d046dc88ab08d9fcf204c5%7C33a15b2f849941998d56f20b5aa91af2%7C0%7C0%7C637818938783874996%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=KxQmoTdwqh0IhlHnHeXmtX0DTIf2kGceX2INux47%2FHM%3D&reserved=0

Roger Bivand
Emeritus Professor
Department of Economics, Norwegian School of Economics,
Postboks 3490 Ytre Sandviken, 5045 Bergen, Norway.
e-mail: Roger.Bivand using nhh.no

More information about the R-sig-Geo mailing list