[R-sig-Geo] st_transform problem

[Deurloo, Marinus]

I have a shapefile from CBS/Kadaster in EPSG:28992 (Amersfoort / RD New) projection and I use st_transform to reproject to EPSG:4326.
Small R script example:


setwd("c:/testTransform")
library("sf")
polyPurmerend <- st_read(dsn = "c:/testTransform", layer = "Purmerend",stringsAsFactors = FALSE)
polyPurmerend_WGS84 <- st_transform(polyPurmerend,4326)
st_write(polyPurmerend_WGS84,dsn=paste0("","Purmerend_WGS84",".gpkg"),"polygon", layer_options="OVERWRITE=YES")

Unfortunately the resulting layer is at a different position.

[cid:7954d268-0209-4785-bc33-32c08ec236d0]

I cannot figure out why.

Rinus Deurloo
University of Amsterdam.




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20181114/03d4c055/attachment.html>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 35402 bytes
Desc: pastedImage.png
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20181114/03d4c055/attachment.png>