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:
summary(saf)
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
