gstat

0熱度

1回答

我想從R包gstat中使用函數variogramST來計算時空變差函數。有12年的數據在空間和時間不規則點有20'000個數據點（沒有全格或部分格）。我必須使用時空包中的STIDF來處理不規則的數據集。我想要一個時間半變異函數，參考點爲0,90,180,270天，直到幾年等。不幸的是，計算和內存問題都會發生。當命令 samplevariogram<-variogramST(formula=for

1熱度

1回答

gstat中使用交叉驗證的長度誤差R中的krigecv

今天我遇到了一個我找不到解決方案的問題。我想在R中使用gstat使用krigcv交叉驗證，並且我有一條錯誤消息： 0％model.frame.default（terms（公式），as（data，「data.frame」），na中的錯誤。行動= na.fail）：可變長度不同（找到'x'）這是我的代碼行： QCini.uk.cv < - krige.cv（QCini $ var1〜x + y +

0熱度

1回答

自動地圖包：負值和變差擬合

我前兩天開始使用自動映射封裝自動化變差函數建模，併發揮它周圍的（這樣做手工是繁瑣和不來太多的優勢......雖然我覺得奇怪的是，automap似乎真的很喜歡Matérn函數來模擬變差函數）。我想到的是數據集中存在負值，而這些數據集中沒有顯示這些數值。這裏我只是粘在一起的人物：NR。 1是使用autoFitVariogram和autoKrige後的輸出。它似乎沒有負面價值。由於我不喜歡圖形輸出，我

0熱度

1回答

R中的IDW參數

我想使用包中的idw命令使用R執行IDW插值。我有這樣的數據： #settings library(gstat) library(dplyr) library(sp) library(tidyr) id_rep <- rep(c(1,2), 20) f <- rep(c(930,930.2), each=20) perc <- rep(c(90, 80), each=10) x

1熱度

1回答

插值到柵格

我有一個包含LAT，LON和溫度數據的三列數據集我想生成一個基於24個數據記錄器數據點預測景觀溫度的柵格圖像。該數據集可以訪問這裏：DATA 這是我這一步嘗試： #Lets try to interpolate the data onto a raster library (raster) library (gstat) library (sp) #Temp and XY data t

1熱度

1回答

R'automap'如何創建用於AutoKrige的預測網格（例如meuse.grid）？

創建預測網格（用於new_data參數）與automap包中的autoKrige函數一起使用時遇到了很多困難。我已經嘗試下面在這個崗位（How to subset SpatialGrid using SpatialPolygon）的步驟，但出現以下錯誤：錯誤X @ COORDS [我，放下= FALSE]：（下標）邏輯標太長另外：警告消息： 1：在min（x）中：沒有非缺少參數min;返回

-1熱度

1回答

用R格子包畫出線條列表

我試圖在xy.plots上面顯示一些變種圖模型。 panel.plots除了我想添加到對應的子圖的行列表之外。 require(gstat) require(sp) data(meuse) names(meuse) #make directional variograms b<-variogram(log(zinc)~1, meuse, alpha = c(0, 45, 90, 135)

1熱度

1回答

局部塊克里格與gstat本地變差函數

我一直無法找到任何信息，具體到本地塊克里格與使用gstat包河還有就是免費從澳大利亞中心精準農業所謂VESPER，能夠局部變差要做到這一點，從我讀過的應該可以在R中，我可以使用一些幫助，把一個for循環，使gstat函數在本地工作。使用默茲數據集作爲一個例子，我已經能夠計算和適應全球變差函數的數據的集合： library(gstat) data(meuse) coordina

0熱度

3回答

當在R中使用gstat或automap包時重複數據

我想使用普通克里金來根據預測變量使用g中的gstat或automap包預測動物的空間預測數據。我有很多（超過100）重複座標點，我不能扔掉，因爲這些臺站多年來多次採樣。每次我運行普通克里格的代碼時，我都會得到一個LDL錯誤，這是由於重複點造成的。有誰知道如何解決這個問題，而不會拋出數據？我已經嘗試了automap包中的代碼，該代碼旨在糾正重複項，但是我無法使其工作。感謝您的幫助！ coordina

1熱度

1回答

R中的簡單變異函數，瞭解gstat :: variogram（）和對象gstat

R中有一個data.frame，它的變量代表位置，其觀測值是這些位置中某個變量的度量值。我想根據距離來測量某些位置的相關衰減，所以變異函數對我的研究特別有用。我想使用gstat庫，但我對某些參數有些困惑。據我理解（經驗）變差函數應該只需要基本的數據：變量的位置觀察這些變量，然後在其它參數，如maximun距離，方向，... 現在，gstat::variogram()功能需要作爲第一個輸入類