替换光栅之间的值

Question

我把我的问题简化了一点，希望它有意义。

我有三个栅格，我正在使用。

library(raster)
LC <- raster("finalRRaster2018.tif")
nativeAET <- raster('kohala_aet_mm_ann_int.tif')
nonnativeAET <- raster('pred_AET.tif`)

栅格LC是一个栅格，其土地覆盖属性值1为本地土地覆盖，0为非本地土地覆盖。

本地蒸散量和非本地蒸散量分别是本地物种和非本地物种的蒸散量。两个光栅的属性都在[015000]之间

id要做的是将1的所有值替换为nativeet值，将0的所有值替换为nonnativeet值。

我的想法是将土地覆盖光栅（值为1或0）转换为AET光栅（值介于0和15000之间）。

我使用了以下代码

native <- values(LC) == 1
nonative <- values(LC) ==0
LC[native] <- values(nativeAET)[native]
LC[nonnative] <- values(nonnativeAET)[nonnative]

但是当我绘制图表时，我得到警告，要替换的项目数不是替换长度的倍数要替换的项目数不是替换长度的倍数

情节是这样的，

我的问题是，我希望生成的图看起来不那么厚。我觉得栅格没有正确对齐，但是当我检查栅格的范围和尺寸时，它们似乎都共享相同的投影。

非常感谢您的帮助！

Answer 1

假设更正了nonative（one"n"）中的错字，您的代码应该可以工作（并且速度很快）。

错误可能来自光栅的空间分辨率差异。因此，<代码>ncell（LC）！=ncell（nativeAET）和ncell（LC）！=ncell（非本地设置）。如果情况并非如此，则需要使用以下方法对光栅重新采样，例如：

nativeAET <- resample(LC, nativeAET, method='ngb')
nonnativeAET <- resample(LC, nonnativeAET, method='ngb')

如果LC光栅中只有两个值（0和1），则可以通过执行以下操作来加快代码速度：

AET <- nativeAET 
AET[values(LC) == 0] <- values(nonnativeAET)[values(LC) == 0] 

Answer 2

如果您提供一些代码生成的示例数据，这个问题将更容易回答。

示例数据：

library(raster)
nativeAET <- nonnativeAET <- LC <- raster(ncol=4, nrow=4)
values(nativeAET) <- 1:16
values(nonnativeAET) <- 21:36
values(LC) <- c(rep(0,8), rep(1,8))

方法1：

AETnative <- mask(nativeAET, LC, maskvalue=1, updatevalue=0)
AETnonnative <- mask(nonnativeAET, LC, maskvalue=0, updatevalue=0)
AET <- AETnative + AETnonnative

方法2：

AET <- mask(nativeAET, LC, maskvalue=1)
AET <- cover(AET, nonnativeAET)