1

我从点数据的形状文件开始。然后我将此格式转换为 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 格式...？

4

1 回答 1

0

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

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

于 2017-01-17T10:43:24.113 回答