Question

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

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

dest: '<%= theme %>/app.js'

然后我想做一些点模式分析，因此使用了以下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)

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

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

接下来我要做的是从上面创建一个Geotiff图像。我想这样做，因为我想在某些矢量数据上叠加tiff图像。

我的问题是：

如何将上面的图像转换为Geotiff格式......？

Answer 1

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

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

使用R将图像转换为地理基准

1 个答案: