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