[R] How to re-project ease( Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?

Jonsson amen.alyaari at Bordeaux.inra.fr
Fri Jan 18 18:11:11 CET 2013

I have nc files  for global soil moisture,here is one file 

from the metadata ,the projection is cylindrical and the resolution is 25
km(it is based on authalic sphere based on International 1924 ellipsoid).As
I want to compare with other data, I have to make them identical.
- my other data are in WGS84 with 0.25*0.25 degree resolution and extent
-180, 180, -90, 90.
So I want to re-project  the data I have here
EASE grid, cylindrical,25 km ,1383 pixel  586 lines
WGS84,0.25*0.25 degree,1440 pixel 720 lines:           
                    A =
        I tried this:
varname = "Soil_Moisture") 
        class       : RasterLayer 
        dimensions  : 586, 1383, 810438  (nrow, ncol, ncell)
        resolution  : 0.2603037, 0.2916659  (x, y)
        extent      : -180, 180, -85.4581, 85.4581  (xmin, xmax, ymin, ymax)
      coord. ref. : +proj=longlat +datum=WGS84 
        names       : Retrieved.soil.moisture.value 
           zvar        : Soil_Moisture 
       a <- spTransform(d, CRS ("+proj=longlat +ellps=WGS84 +datum=WGS84
but I got this error:
     Error in function (classes, fdef, mtable)  : 
              unable to find an inherited method for function ‘spTransform’
for signature ‘"RasterLayer", "CRS"’        

View this message in context: http://r.789695.n4.nabble.com/How-to-re-project-ease-Equal-Area-Scalable-Earth-grid-with-a-25-km-cylindrical-projection-to-WGS84-0-tp4655977.html
Sent from the R help mailing list archive at Nabble.com.

More information about the R-help mailing list