r - R-将SpatialLines转换为栅格

Question

在R中，我们可以使用功能raster将SpatialLinesDataFrame转换为rasterToCountour：

library(raster)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
x <- rasterToContour(r)
class(x)

[1] "SpatialLinesDataFrame"
attr(,"package")
[1] "sp"

spplot(x)

在R中，是否有相反的方法？像 contourToRaster一样？

我们可以简单地获取与直线上每个点关联的字段值，但是我正在寻找更通用的东西，可以在直线之间进行插值并在定义的域上生成完整的栅格。

Answer 1

library(raster)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
x <- rasterToContour(r)

您可以栅格化值。在这种情况下，请先从因子标签中提取它们。

x$value <- as.numeric(as.character(x$level))
rr <- rasterize(x, r, "value")

然后提取单元格值并进行插值

xyz <- rasterToPoints(rr)

（如果您想跳过rasterize和rasterToPoints（如mikoontz建议），则可以改而做

#g <- geom(x) 
#xyz = cbind(g[, c("x", "y")], x$value[g[,1]])

以更复杂的模型为代价）

现在进行插值，例如使用Tps

library(fields) 
tps <- Tps(xyz[,1:2], xyz[,3])
p <- raster(r)
p <- interpolate(p, tps)

m <- mask(p, r)
plot(m)