替代在使用栅格对象时提取函数

作者: 無口君
发布时间: 2024-12-31 06:21:41 (1月前)
转自：

2 条回复

0#
回复此人
满目山河 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 你可以使用 <code> zonal </code> 方法： </p> <pre> <code> library(raster) cat_r <- raster(ncol=3,nrow=3, xmn=-10, xmx=10, ymn=-10, ymx=10, vals=c(1,2,1,3,4,3,4,4,4 )) val_r <- setValues(cat_r, c(1,0,1,5,2,5,2,2,2)) zonal(val_r, cat_r, "sum") # zone sum #[1,] 1 2 #[2,] 2 0 #[3,] 3 10 #[4,] 4 8 </code> </pre> <P> 这相当于 </p> <pre> <code> s <- stack(cat_r, val_r) v <- values(s) tapply(v[,2], v[,1], sum) # 1 2 3 4 # 2 0 10 8 </code> </pre> </DIV>

编辑

登录后才能参与评论