rgdal

To ensure that appropriate projection metadata are at every step associated with the coordinates, I’d suggest converting the points to a SpatialPointsDataFrame object as soon as possible. See ?”SpatialPointsDataFrame-class” for more on how to convert simple data.frames or matrices to SpatialPointsDataFrame objects. library(sp) library(rgdal) xy <- data.frame(ID = 1:2, X = c(118, 119), Y = … Read more