Skip to content
Prev 28855 / 29559 Next

Mixing points and sf objects

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)

Best,

Patrick