Hi You use several functions which are not base without mentioning their source, therefore you probably won't get much help. I wonder that you do not get error with stack.
AFAIK you will get a character vector in dir2. stack a character vector results in error > stack(letters) Error in rep.int(names(x), lapply(x, length)) : invalid 'times' value > Normally stack works on lists or data frames > test teorie kal1 kal2 kal3 1 2 2.04 1.80 1.851 2 4 4.07 3.50 3.753 3 6 5.86 5.23 5.663 4 8 7.87 7.09 7.598 5 10 10.01 8.87 9.428 6 12 NA 10.61 11.308 > stack(test) values ind 1 2.000 teorie 2 4.000 teorie 3 6.000 teorie 4 8.000 teorie 5 10.000 teorie 6 12.000 teorie 7 .... So you either have some other stack function or I am completly out of the topic. Regards Petr > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Jonsson > Sent: Monday, March 18, 2013 6:32 PM > To: r-help@r-project.org > Subject: [R] Why stacking rasters return NAs? > > 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. > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting- > guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
