Dear all,
I tried to project a Raster file (orignal projection: UTM) to a file
in the Goode projection.
The re-projection with proj4 was obviosly not correct, neither was a
re-projection from a geographical lon/lat file.
Does anyone know if this is a general problem with R/proj4?
Otherwise, I might have produced a faulty raster file. I used the USGS
Global Ecosystem Data (http://edc2.usgs.gov/glcc/globdoc2_0.php),
provided as raw binary, and had to create the Raster myself:
ny <- 17347
nx <- 40031
goge <-
raster(xmn=-20015500,xmx=20015500,ymn=-8673500,ymx=8673500,nrows=ny,ncols=nx,crs="+proj=goode
+a=6370997")
values(goge) <- binfromFTP
Although the created Raster File looked good to me naturally the error
could still occur in this step.
Afterwards I just used:
projectRaster(infile,goge, method="ngb")
Any hints are very much appreciated!
Kind regards,
Birgit