Blok spatio-temporal kriging dengan paket R gstat?

9

Dalam dokumentasi spsaya melihat bahwa coordsargumen SpatialPolygonsfungsi hanya menerima 2 kolom. Apakah tidak mungkin membuat poligon 3 dimensi? Saya melakukan spatio-temporal kriging menggunakan paket gstatdengan menambahkan waktu sebagai dimensi spasial ketiga dari a SpatialPointsDataFrame. Saya ingin memperluas itu untuk memblokir kriging dengan poligon sebagai blok. Apakah ada cara untuk melakukan itu?

Nima
sumber

Jawaban:

5

Dalam sp, SpatialPoints*, SpatialPixels*dan SpatialGrid*(dengan *dihilangkan atau diganti dengan DataFrame) melakukan dukungan lebih dari 2 dimensi spasial, seperti yang dilakukan OP, tetapi SpatialPolygons*dan SpatialLines*tidak. Dengan gstatAnda dapat melakukan kriging blok 3-D dengan blok 3-D (menggunakan block = c(10,10,10)), tetapi Anda tidak dapat melakukan ini untuk blok non-persegi panjang, seperti yang diinginkan OP. Sangatlah OK untuk mengganti waktu untuk dimensi ketiga, tetapi Anda dibatasi untuk variogram ST metrik.

library(gstat)
vignette("st")

memberi Anda lebih banyak opsi untuk model variogram, tetapi tidak untuk memprediksi nilai rata-rata blok (ini FYI, bukan jawaban untuk pertanyaan).

Satu-satunya jawaban untuk pertanyaan itu adalah melakukan simulasi kondisional 3D, dan nilai-nilai titik agregat di atas blok 3D Anda yang sewenang-wenang (2D poligon + batas waktu). Membosankan, tetapi mungkin; juga hanya di sepanjang jalur 3D, tidak di sepanjang jalur yang dijelaskan dalam sketsa ST ( krigeSTtidak melakukan simulasi - belum!).

Edzer Pebesma
sumber
3

Anda benar-benar perlu melakukan sedikit riset lebih lanjut tentang metodologi Anda dan membaca dokumentasi untuk memahami struktur objek kelas sp S4 dan interaksi objek sp dengan fungsi gstat yang relevan. Dalam sp Vignette ada penjelasan rinci tentang perbedaan antara SpatialPolygons (hanya topologi poligon) dan SpatialPolygonDataFrame (poligon dengan atribut) objek.

Apa yang Anda jelaskan tidak memblokir Kriging dan menggunakan waktu sebagai atribut tidak menghasilkan estimasi spasial-temporal. Metodologi konseptual yang Anda jelaskan sangat tidak valid. Menggunakan poligon atau poligon centroid melanggar asumsi Kriging dari bidang acak seragam, anisotropi, dan nonstasioneritas.

Berikut ini adalah sketsa gstat yang bagus pada model spasial-temporal menggunakan antarmuka ke paket ruangwaktu. Saya juga harus mencatat bahwa paket constrainedKriging dapat melakukan blok Kriging pada blok berbentuk sewenang-wenang menggunakan fungsi mean nonstasioner dan variogram stasioner isotropik lemah.

Yang mengatakan, untuk menjawab pertanyaan Anda, Anda bisa mengirimkan objek sp SpatialPointsDataFrame langsung ke model variogram / Kriging di gstat. Dalam jenis objek sp, atribut berada di slot "data" dan sudah melekat pada koordinat melalui struktur kelas S4 internal.

# COERCE meuse DATAFRAME TO sp SpatialPointsDataFrame OBJECT
require(gstat)
data(meuse)
coordinates(meuse) <- ~ x + y
head(meuse@data)

# CREATE SEMIVARIOGRAM USING THE zinc ATTRIBUTE
# NOTE: THERE IS NO ARGUMENT FOR A "4th DIM"
v <- variogram(log(zinc) ~ 1, meuse)
plot(v, type = "l")
Jeffrey Evans
sumber
0

Waktu akan menjadi dimensi keempat, z level atau tinggi akan menjadi yang ketiga, tins atau dems atau z level menangani tinggi dalam gis, waktu mungkin diwakili oleh database berversi. Anda bisa memetakan poligon ke kaleng atau berpotongan dengan dem, atau mengaitkan kolom az ke kolom xy yang mewakili titik-titik bentuk poligon. Mungkin Anda bisa menambahkan bidang waktu juga dengan rentang nilai dan perubahan nilai sepanjang waktu, ditambahkan sebagai kolom. Tentu saja ESRI mungkin memiliki alat untuk itu.

Lewis
sumber
Saya tidak berpikir Analis Geostatistik mendukung kriging ruang-waktu.
Paul Hiemstra