Re: Using CRS() method on SpatialPointDataFrame to project coordinates
On Thu, 20 Sep 2018, MacQueen, Don wrote:
> First, it looks like prcp_sites_ll does not have CRS information, because
> @projargs is NA.
I used coordinates(prcp_sites_ll) <- c('easting','northing') to transform
the original dataframe to a SpatialPointsDataFrame and thought that because lat-lon
was not projected the @projards=NA was expected. Now I have learned that ...
> You'll need something like
> proj4string(prcp_sites_ll) <- CRS('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0')
I missed this step in my readings.
> Otherwise, your spTransform() command looks ok to me.