Я борюсь с проблемой Kriging + GAM и надеюсь, что кто-нибудь сможет мне помочь.
Я пытаюсь сделать интерполяцию загрязнителя в область. Для этого я использую обычный кригинг из gstat, в дополнение к GAM, на основе этого метода (ссылка, где я нашел его ниже):
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? Все методы, которые я вижу, включают в себя существующий растр или шейп-файл, которого у меня нет (или я не знаю, как его создать).
с тех пор, как я могу добавить границы страны (интерполяция выполняется в северной атланти c, и я хочу добавить страну, которая окружает область)
Конечно, дайте мне знать, если вам нужна дополнительная информация.
1024 * Best 1026 * Селин 1028 * ссылка:
https://www.r-exercises.com/2018/03/31/advanced-techniques-with-raster-data-part-3-regression-kriging/? cf_chl_jschl_tk = 6311eac67bd8c8522c68f83f62be077ef572ca2 c -1588325277-0-AaBAYr6yFHMceWCKlsRhOxbdPIKI3oKq6c8PVrtfzQaBgUDqBu3-QG8E1z94b3ygHSpaaYMPtHYo4ThfcJWwsn2p1ZxcvutlBdn- nh0wKILVFc3xrLHUzNqxRqJlY9sjRcHDft9tn9PFPngo06AonPkEkCUIHtQ1mgjz3Kgt9zTp93OJyJEgAp0XRAvc1pzE8MTXv-fAwlkWki21vE9RYiZEgss5m5AYJ-fppvRUL2ckSTl0W765HVVeH3GySp6nuoTvmvaocZ_Q8axMB_za2Iiiwb2iar_h-9s8fOZXh5QbQ c -9PojDlPBfU-nxy4E1gGHjFyfXNvnkghKYqZaz6fr4pQ7u7K3R-GglVx2vPlEmIfuIzHYdhAC4oJd84hKMvg )