Question

我需要根据以下逻辑构造一个新的栅格（称为ch，它基于两个栅格（分别称为tmin和tmax）。新栅格的每个元素都遵循此逻辑。tmin, tmax和ch都是尺寸相同的光栅砖，但我认为我首先会只尝试单层。

if (tmin > 7) ch <- 0
if (tmax < 7) ch <- 24
else
ch <- (7 - tmin)/(tmax - tmin)

我尝试通过以下代码遵循此answer的逻辑。 testout如预期那样由零或24组成。我无法弄清楚如何将逻辑的最后一部分纳入该过程。

library(raster)
r1 <- raster()
r2 <- raster()
set.seed(10)
values(r1) <- runif(ncell(r1), min = 0, max = 10)
values(r2) <- runif(ncell(r2), min = 0, max = 10)

myFun <- function(r1, r2) {
  0*(r1 > 7)+
  24*(r2 < 7)  
}

testout <- myFun(r1, r2)

Answer 1

示例数据

library(raster)
r1 <- r2 <- raster(ncol=10, nrow=10)
set.seed(10)
values(r1) <- runif(ncell(r1), min = 0, max = 10)
values(r2) <- runif(ncell(r2), min = 0, max = 10)

向量化功能

myFun <- function(tmin, tmax) {
    ch <- (7 - tmin)/(tmax - tmin)
    ch[tmin > 7] <- 0
    ch[tmax < 7 & tmin <= 7] <- 24
    ch
}

ch1 <- overlay(r1, r2, fun=myFun)

未记录的替代方法

ifel <- raster:::.ifel
ch2 <- ifel(r1 > 7, 0, ifel(r2 < 7, 24, (7 - r1)/(r2 - r1)))

R栅格从两个现有栅格构建新栅格

1 个答案: