我在 R 中有两个栅格（Landsat slc-off 图像）。两者都缺少一些数据，但间隙位置完全偏移。例如，我创建了两个 r1及r2以下的栅格。
r1 <- raster(system.file("external/test.grd", package="raster"))
r1_mat <- as.matrix(r1)
r1[which(!is.na(as.matrix(r1)))[1:600]] <- NA
par(mfrow=c(3,1))
plot(r1)
r2 <- raster(system.file("external/test.grd", package="raster"))
r2[which(!is.na(as.matrix(r2)))[900:1400]] <- NA
plot(r2)
然而，第二张图像是在不同的大气条件下拍摄的，比如在太阳光照较好的情况下。要模拟这种效果（非常简单）：
r2 <- r2 + 200
现在我想用 r2 像素填充 r1 中的缺失数据，这些像素与间隙重叠，这很简单：
r3 <- cover(r1, r2)
plot(r3)
查看 的结果，本页图 2 中也显示的plot(r3)所谓“剥离效果”是显而易见的。推荐的解决方案之一是在进行间隙填充之前通过匹配它们的直方图来规范化这两个图像。该技术基于候选图像的累积分布函数，例如将在 Grass 7 中完成的。
我怎样才能在 R 中实现这一点？谢谢！