R パッケージ gstat の関数 variogramST を使用して、時空間バリオグラムを計算しようとしています。
空間と時間の不規則なポイントに 20,000 のデータ ポイントを持つ 12 年間のデータがあります (完全なグリッドまたは部分的なグリッドはありません)。不規則なデータセットには、時空パッケージの STIDF を使用する必要があります。0、90、180、270 日、最大数年などの基準点を持つ時間セミバリオグラムが欲しいのですが、残念ながら、計算とメモリの問題の両方が発生します。コマンドが
samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1)
追加の引数なしで実行すると、セミオバリオグラムは、セミバリオグラムの基準点に関して非常に短い期間のみを考慮しているため、固有のデータ構造を適切にキャプチャしていないようです。
ユーザーが自由に使用できるこの関数の引数は他にもありますが、それらを正しくパラメーター化する方法がわかりません: tlag、tunit、twindow。具体的には、それらがどのように相互作用し、上記の目標をどのように達成するのか疑問に思っています. だから私は次のコードを試しました
samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag= ...., tunit=... , twindow= ...)
次のコードの結果は、32G バイトの RAM コンピューターのメモリの問題により機能しません。
samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag=90*(0:20), tunit="days")
それ以外の場合は、おそらく欠陥がある可能性があります。さらに、コードの最後の行も、計算時間の点で実行不可能に見えます。
希望の時間間隔を目指して、gstat パッケージから variogramST 関数を正しく指定する方法を誰かが知っていますか?
ありがとう