我正在尝试对与lat / long值相关联的z值执行高斯平滑处理。我拥有的每个纬度/长度都是1英里乘1英里网格的质心。这是一个数据示例:
latitude <- c(40.607674, 40.607048, 40.606419, 40.605785, 40.605152, 40.604515, 40.603874)
longitude <- c(-91.865349, -91.846184, -91.827026, -91.807861, -91.788704, -91.769547, -91.750381)
value <- 1:7
这些纬度/长度在WGS 1984坐标系中。需要注意的一件重要事情是我无法访问rgdal库,但是我可以访问ARCGIS(我没有空间分析师许可证)来将这些值投影到另一个坐标系。
到目前为止,我已经能够使用这些值在R中创建SpatialPointsDataFrame。我的目的是将SpatialPointsDataFrame转换为高斯平滑的矩阵表示,但我还没有成功地做到这一点。