Saya telah mengikuti tutorial online untuk sprigal kriging dengan keduanya geoR
dan gstat
(dan juga automap
). Saya dapat melakukan kriging spasial dan saya memahami konsep utama di baliknya. Saya tahu cara membuat semivariogram, cara menyesuaikan model dengan itu dan bagaimana melakukan kriging biasa.
Yang tidak saya mengerti adalah bagaimana bobot dari nilai yang diukur di sekitarnya ditentukan. Saya tahu mereka berasal dari semivariogram dan tergantung pada jarak dari lokasi prediksi dan pada pengaturan spasial dari titik-titik yang diukur. Tapi bagaimana caranya?
Adakah yang bisa membuat model kriging (non-bayesian) biasa dengan 3 titik pengukuran acak dan 1 lokasi prediksi? Itu akan mencerahkan.
prediction
spatial
weights
kriging
Pigna
sumber
sumber
Jawaban:
Pertama saya akan menjelaskan kriging biasa dengan tiga poin secara matematis. Asumsikan kita memiliki bidang acak intrinsik stasioner.
Kriging Biasa
Kami mencoba memprediksi nilainyaZ(x0) menggunakan nilai yang diketahui Z= ( Z(x1) , Z(x2) , Z(x3) ) Prediksi yang kita inginkan adalah dari bentuk
Perilaku bobot yang tepat sulit dilihat hanya dari persamaannya, tetapi orang dapat dengan sangat kasar mengatakan:
Namun saya akan fokus pada lokasi titik di pesawat. Saya menulis fungsi R kecil ini yang mengambil poin dari[ 0 , 1]2 dan plot bobot kriging (untuk fungsi kovarian eksponensial dengan nol nugget).
Anda dapat bermain dengannya menggunakan
clickppp
fungsi spatstat :Berikut ini beberapa contohnya
Poin berjarak sama darix0 dan dari satu sama lain
Poin yang dekat satu sama lain akan berbagi bobot
Titik terdekat "mencuri" bobot
Dimungkinkan untuk mendapatkan bobot negatif
Semoga ini memberi Anda merasakan bagaimana bobot bekerja.
sumber