I have a NetCDF (nc) file I'm using to create a Raster stack in R. The coord. ref. is listed as "NA" when I load it in, but plots as Polar Stereographic.
When I try set its CRS, it lists that it has changed in the Raster info, but still plots incorrectly. How do I change the CRS to WGS1984 and plot it as Equirectangular?
I need to crop it using lat/long coordinates etc, so can't really leave it as stereographic projection.
Data is available here https://drive.google.com/file/d/0Bz0W7Ut_SNfjeE9ObXpySzJ5UWs/view
Code used is:
library(ncdf4) library (raster) library(sp) nc<-nc_open('NOAA_Snow_JanJune2016.nc') print(nc) varname<-names(nc$var) r<-stack('NOAA_Snow_JanJune2016.nc',varname='snow_cover_extent') crs(r) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" plot(r[])
(I can provide images of plots etc if you don't want to run the above)