使用 R 将图像转换为 geotiff

Question

我从点数据的形状文件开始。然后我将此格式转换为 ppp 格式以使用一些 spatstat 函数。

下面列出了执行此操作的代码：

sp_points_df_al <- rgdal::readOGR(dsn = shape_path, layer = "aland_points")
# convert the spatial points data.frame to just plain ol spatial points
sp_points_al <- as(sp_points_df_al, "SpatialPoints")
# now convert the spatial points to ppp 
ppp_al <- as(sp_points_al, "ppp")
# estimate the window from the points data
wind_rr <- spatstat::ripras(ppp_al)
ppp_al_constr <- spatstat::ppp(ppp_al$x, ppp_al$y, window = wind_rr)

然后我想做一些点模式分析，所以使用了以下 spatstat 函数：

image <- spatstat::density.ppp(ppp_al_constr, sigma = 0.004, dimyx=c(512, 512))

然后我使用以下内容显示图像：

bias_palette <- colorRampPalette(c("blue", "magenta", "red", "yellow", "white"), bias=2, space="Lab")
spatstat::plot.im(k, col=bias_palette(256), ribbon = FALSE)

接下来我要做的是从上面创建一个 Geotiff 图像。我想这样做是因为我想将 tiff 图像覆盖在一些矢量数据之上。

我的问题是：

如何将上面的图片转换成 Geotiff 格式...？

Answer 1

spatstat::density.ppp 的结果是类im 的对象。如果你对这个对象执行as.matrix，你会得到一个包含像素条目的矩阵。然后，您可以将其转换为任何其他类型的图像。图像对象的as.data.frame(Frame(..)) 也将检索边界框。

或者，转到处理spatstat 对象的maptools 包，并查找将spatstat 类im 对象转换为所需类的函数。