我有一个包含 900K 随机采样坐标点的数据集,每个点都有一个值。我想在地图上制作一个网格,并为每个单元格分配单元格内所有点的平均值。最后,我想将网格绘制为热图。
网格的分辨率可以是 100 平方公里。小区可达 10000 平方公里。
首先,您的调用语法autoKrige
不正确,应该是:
kr <- autoKrige(mean_week1~1,project_df)
该formula
参数从不包含向量,它始终包含来自project_df
.
autoKrige
可以预测任何Spatial*
对象,例如SpatialGrid
或SpatialPoints
,只需将要预测的对象传递给 的new_data
参数autoKrige
。因此,如果您想预测不规则点,只需创建一个SpatialPoints
包含位置的对象并将其传递给autoKrige
.
关于计算时间,900.000 个点是一个大数据集。nmax
您可以使用(使用的最大点数)或maxdist
(使用点的最大距离)参数来限制用于预测的数据量。这些是 , 的一部分,并从包...
中传递给。krige
gstat