0

一段时间以来，我一直在努力解决 Kriging + GAM 问题，希望有人能够帮助我。

我正在尝试将污染物插值到一个区域中。为此，除了 GAM 之外，我还使用 gstat 中的普通克里金法，基于此方法（我找到它的链接如下）：

GAM <- gam(formula = pollutant ~ 1+ s(Long, Lat) + s(Individual, bs="re"), data = dataset)
GAMPredictions <- predict(mydataset, GAM, type="response")
residKrigMap <- krige(formula = residuals(GAM) ,
                  locations = mydataset, 
                  model = myVariogram,
                  newdata = MyGrid)
residKrigRstLayer <- as(residKrigMap, "RasterLayer")   
gamKrigMap <- GAMPredictions + residKrigRstLayer

所有克里金部分都做得很好，但是当我尝试这样做时： 

gamKrigMap <- GAMPredictions + residKrigRstLayer

这没用。而且我很确定这是因为我正在对网格进行插值。事实上，我正在使用一个 rds 文件，在该文件上我要插入的区域。但是您不能在这样的文件中进行预测，因为它不包含数据集中的信息。此外，我的数据集中有 760 个观察值，但网格中有超过 6000 个观察值，因此存在不匹配。到目前为止，对于我一直在研究的内容，我应该根据我的数据集创建一个网格。

所以我的问题是：

  • 如何使用 R 从我的数据集中创建网格？我看到的所有方法都包括我没有（或知道如何创建）的现有栅格或 shapefile。

  • 从那时起，我如何添加国家边界（插值是在北大西洋进行的，我想添加该地区周围的国家）

当然，如果您需要更多信息，请告诉我。

最好的，

席琳

链接：https://www.r-exercises.com/2018/03/31/advanced-techniques-with-raster-data-part-3-regression-kriging/? cf_chl_jschl_tk =6311eac67bd8c8522c68f83f62be077ef572ca2c-1588325277-0-AaBAYr6yFHMceWCKlsRhOxbdPIKI3oKq6c8PVrtfzQaBgUDqBu3-QG8E1z94b3ygHSpaaYMPtHYo4ThfcJWwsn2p1ZxcvutlBdn-nh0wKILVFc3xrLHUzNqxRqJlY9sjRcHDft9tn9PFPngo06AonPkEkCUIHtQ1mgjz3Kgt9zTp93OJyJEgAp0XRAvc1pzE8MTXv-fAwlkWki21vE9RYiZEgss5m5AYJ-fppvRUL2ckSTl0W765HVVeH3GySp6nuoTvmvaocZ_Q8axMB_za2Iiiwb2iar_h-9s8fOZXh5QbQc-9PojDlPBfU-nxy4E1gGHjFyfXNvnkghKYqZaz6fr4pQ7u7K3R-GglVx2vPlEmIfuIzHYdhAC4oJd84hKMvg )

4

0 回答 0