search for: rasterfinal

Displaying 1 result from an estimated 1 matches for "rasterfinal".

2013 Jul 17
2
Using RasterBricks
Dear listers, I am trying to create a RasterLayer of the values of a rasterbrick object. The rasterbrick object has, for example, 100cells library(raster) r <- raster(ncol=10, nrow=10) r[]=1:ncell(r) s <- brick(r,r,r) s <- s * 1:3 Each cell of the rasterfinal will have the AREA UNDER CURVE formed by the values in each cell of the three original rasters. For example: s[4] has the values: 4,8,12 y<-s[4]; x<-0:2; You can see the AUC in this graph: plot(x,y, type='l') And calculate it using the function trapz (package:pracma) library(pracm...