[R] Why stacking rasters return NAs?

Jonsson amen.alyaari at Bordeaux.inra.fr
Mon Mar 18 18:32:26 CET 2013


I have several rasters that I want to do some calculations ,basically
calculating the moving average. 

            dir2 <- list.files("D:\\2010+2011", "*.bin", full.names = TRUE)
           saf=stack(dir2)
          movi <-  overlay(stack(saf),fun=function(x) movingFun(x, fun=mean,
n=3,  na.rm=TRUE))
             Error in .overlayList(x, fun = fun, filename = filename, ...) : 
         cannot use this formula, probably because it is not vectorized
 I then checked the data but found that all values were returnd as NA and
this may explain why i am getting the error. 
        
        saf
       class       : RasterStack 
       dimensions  : 720, 1440, 1036800, 601  (nrow, ncol, ncell, nlayers)
       resolution  : 0.25, 0.25  (x, y)
       extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
      coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0
+no_defs 
     names       : Vol_025_H//00_1_wgs84, Vol_025_H//00_1_wgs84,
Vol_025_H//00_1_wgs84,    Vol_025_H//00_1_wgs84, Vol_025_H//00_1_wgs84,
Vol_025_H//00_1_wgs84, , ... 
    min values  :                    NA,                    NA,                   
NA,                       NA,                    NA,                    NA,                   
NA,                     NA,                    NA,                    NA,                   
NA,                    NA,                    NA,                    NA,                   
NA, ... 
    max values  :                    NA,                    NA,                   
NA,                    NA,                    NA,                    NA,                   
NA,                    NA,                    NA,                    NA,                   
NA,                       NA,                    NA,                    NA,                   
NA, ... 


I wonder why this is happening, I checked the files separably(summary) and
everything was right!as you can see bellow:

             ol_025_H14_2011092000_1_wgs84 Vol_025_H14_2011092100_1_wgs84  
Vol_025_H14_2011092200_1_wgs84 Vol_025_H14_2011092300_1_wgs84
Vol_025_H14_2011092400_1_wgs84
      Min.                           0.00000                      0.0000000                           
0.0000000                      0.0000000                      0.0000000
     1st Qu.                        0.31883                      0.3163167                     
0.3146436                      0.3113111                      0.3064551
     Median                      9999.00000                   9999.0000000                  
9999.0000000                   9999.0000000                   9999.0000000
     3rd Qu.                     9999.00000                   9999.0000000                  
9999.0000000                   9999.0000000                   9999.0000000
      Max.                        9999.00000                   9999.0000000                  
9999.0000000                   9999.0000000                   9999.0000000
     NA's                           0.00000                      0.0000000                       
0.0000000                      0.00000

I am gratful to anyhelp



--
View this message in context: http://r.789695.n4.nabble.com/Why-stacking-rasters-return-NAs-tp4661706.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list