Hi there,I am currently trying to calculate "weighted" spatial annual "global" values for precipitation using the object "Prec20". I need to do this for each of the 140 years (i.e. 140 layers) of global precipitation data that I have. The idea would be to somehow apply weights to each grid cell value for each year by using the cosine of its latitude (which means grid cells at the equator would have a weight of 1 (i.e. the cosine of 0 degrees is 1), and the poles would have a value of 1 (as the cosine of 90 is 1)). The idea would also then be to make a time series of these values, from Year 1 to Year 140, after all newly derived grid cell values are averaged for each year, creating 140 (weighted) averages)). The object "Prec20" looks like this:class : RasterBrick

dimensions : 64, 128, 8192, 140 (nrow, ncol, ncell, nlayers)

resolution : 2.8125, 2.789327 (x, y)

extent : -181.4062, 178.5938, -89.25846, 89.25846 (xmin, xmax, ymin, ymax)

coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

data source : C:/Users/Travis/Documents/Other documents/All netCDF files/netcdffiles/MaxPrecIPSLIPSL-CM5B-LR1pctCO2.nc

names : X1, X2, X3, X4, X5, X6, X7, X8, X9, X10, X11, X12, X13, X14, X15, ...

z-value : 1, 140 (min, max)

varname : onedaymax

Here is what I have done so far:ncfname20 <- "MaxPrecIPSLIPSL-CM5B-LR1pctCO2.nc"Prec20 <- brick(ncfname20,var="onedaymax")w20 <- cos(Prec20*(pi/180))

Would this approach apply the weights appropriately?

I would greatly appreciate any help with this!!!!Thanks,

[[alternative HTML version deleted]]

_______________________________________________

R-sig-Geo mailing list

[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo