[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