If your original data is not in the desired projection, changing the CRS attributes will not do the transformation for you. In the case of vector data, you will need to use spTransform(). For raster data, you might be able to use projectRaster() but I would recommend using gdalwarp().
I wrote a brief blog post about this question, here: https://geoanalytic.github.io/coordinate-reference-systems-in-r/ Regards, David Currie, P.Eng [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
