Metode open source untuk kriging?

42

Saya memiliki dataset titik yang ingin saya Krige , idealnya menggunakan paket perangkat lunak sumber terbuka. Jika memungkinkan, saya juga ingin memilih model semi-variogram selama proses untuk meningkatkan estimasi.

scw
sumber

Jawaban:

45

Tergantung pada jenis Kriging yang ingin Anda terapkan, ada paket berbeda untuk dipilih:

Kriging Biasa

Versi paling umum diterapkan misalnya di:

Kriging sederhana

Simple Kriging menggunakan rata-rata dari seluruh kumpulan data sementara Biasa Kriging menggunakan rata-rata lokal. Oleh karena itu, Simple Kriging bisa kurang akurat, tetapi umumnya menghasilkan hasil yang "lebih halus". Diimplementasikan di:

Kriging Universal

Universal Kriging memungkinkan pertimbangan drift dalam data. Implementasi termasuk dalam:

Jenis Kriging Lainnya

GRASS v.krige juga mendukung Block Kriging.

HPGL mengimplementasikan sejumlah besar metode Kriging yang kurang dikenal (lihat manual untuk informasi lebih lanjut tentang itu):

  • Indikator Kriging (IK)
  • Variatif Lokal Kriging Bervariasi (LVM Kriging)
  • CoKriging Sederhana (Markov Model 1 & 2)
  • Simulasi Indikator Berurutan (SIS)
  • Corellogram Local Varying Mean SIS (CLVM SIS)
  • SIS Bervariasi Lokal (SIS LVM)
  • Simulasi Gaussian Berurutan (SGS)
  • Simulasi Gaussian Terpotong (GTSIM) [dalam koleksi skrip Python]

SAGA menawarkan versi berbeda dari Biasa dan Universal Kriging.

Gstat krige juga mendukung Block dan Point Kriging.

underdark
sumber
1
Jawaban Kriging yang luar biasa!
Ragi Yaser Burhum
1
Tautan fauna Anda tidak berfungsi ...
Alex Leith
17

Sepertinya ada beberapa opsi dengan GRASS GIS. Lihatlah halaman GRASS Kriging Wiki: http://grass.osgeo.org/wiki/Kriging

Proyek Google Summer of Code pada tahun 2009 menghasilkan V.krige: http://grass.osgeo.org/wiki/V.krige_GSoC_2009

Paket GPL gstat harus bekerja sendiri atau dihubungkan dengan GRASS GIS. http://www.gstat.org/

Dylan Beaudette memiliki contoh yang bagus untuk melakukan kriging dengan GRASS. http://casoilresource.lawr.ucdavis.edu/drupal/node/438 (Blognya penuh dengan contoh-contoh hebat dan menarik menggunakan OpenSource GIS dan alat statistik!)

DavidF
sumber
14

The R-proyek memiliki sejumlah besar paket perangkat lunak statistik spasial , tetapi R memiliki kurva belajar agak curam.

dodobas
sumber
Orang selalu mengatakan itu, tetapi saya bertanya-tanya: curam relatif terhadap apa?
Matt Parker
Saya telah melihat komentar "kurva pembelajaran curam" yang dilemparkan ke R beberapa kali - itu tidak masuk akal bagi saya. Saya adalah satu tahun dalam hubungan saya dengan MATLAB ketika saya menemukan R. Saya menemukan R sangat mudah dipelajari sehingga saya memberi MATLAB salut satu jari dan segera berhenti menggunakannya dengan berat.
Sharpie
1
saya pikir itu karena orang jarang mencoba untuk memahami statistik, dan karena itu memiliki kurva belajar kurva curam bijaksana, jarang ada masalah mengambilnya
dodobas
2
Saya pikir sintaks bijak itu adalah salah satu bahasa yang lebih mudah dipelajari. Apa contoh bahasa statistik yang mudah dipelajari dari baris perintah. Saya pikir orang mengeluh karena itu bukan Excel.
TheSteve0
1
Ini adalah langkah dibandingkan dengan program berbasis GUI. Jika Anda menggunakan Windows sepanjang hidup Anda, dan program berbasis GUI, Anda akan berjalan ketika Anda melihat tampilan baris perintah. Perbandingan Excel memang membuat mereka berjalan. Tetapi R sangat mudah digunakan jika seseorang dapat menunjukkan trik dasar kepada Anda. Anda harus siap untuk mempelajari konsep-konsep baru seperti vektor, matriks, fungsi, loop, yang dalam dunia Excel / Windows tidak ada. Jika sebelumnya Anda menggunakan Linux, itu bukan kurva langkah.
Niculita Mihai
14

Jika Anda senang membaca raster Anda menjadi array numpy ( gdal bisa melakukan ini), maka Anda bisa menggunakan implementasi High Performance Geostatistics Library dari Python atau C / C ++.

HPGL mengimplementasikan algoritma berikut:

  1. Kriging Sederhana (SK)
  2. Kriging Biasa (OK)
  3. Indikator Kriging (IK)
  4. Variatif Lokal Kriging Bervariasi (LVM Kriging)
  5. CoKriging Sederhana (Markov Model 1 & 2)
  6. Simulasi Indikator Berurutan (SIS)
  7. Corellogram Local Varying Mean SIS (CLVM SIS)
  8. SIS Bervariasi Lokal (SIS LVM)
  9. Simulasi Gaussian Berurutan (SGS)
  10. Simulasi Gaussian Terpotong (GTSIM) [dalam koleksi skrip Python]

Saya belum menggunakannya sendiri tetapi telah mendengar hal-hal baik tentang itu, terutama sehubungan dengan kecepatan.

fmark
sumber
4

Saya ingat menggunakan SAGA untuk melakukan ini beberapa tahun yang lalu untuk beberapa hasil pemodelan banjir. Open Source dan layak untuk dilihat.

Simon
sumber
4

gvSIG (GIS gratis lain) memungkinkan kriging, menggunakan Sextante. Ini pada dasarnya sama dengan menggunakan SAGA, tetapi gvSIG memberikan pengalaman gis yang lebih 'khas' (seperti ESRI).

Rudi
sumber
3

Anda dapat mencoba model Kriging di Surfpack versi 1.1 (saya menulisnya ketika saya masih di tim DAKOTA), atau versi terbaru dan terhebat yang hadir dengan versi "stabil" DAKOTA (Surfpack adalah sub-paket DAKOTA) , ini memang universal Kriging dari perspektif fungsi korelasi daripada semi-variograms.

Baru-baru ini seorang pengguna, Joel Guerrero, membandingkannya secara langsung dengan sekelompok implementasi lain dan menyatakan bahwa "Selalu terkait dengan surfpack, kami membandingkannya dengan implementasi lain (termasuk yang komersial), dan sejauh ini mengungguli semuanya, ke titik yang kadang-kadang tampaknya melakukan sihir hitam "

Keith Dalbey
sumber
3

GSLIB (Perpustakaan Perangkat Lunak Geostatistik) adalah perangkat lunak yang digerakkan oleh file / perintah yang dikembangkan dari Universitas Stanford dan dirilis pada 1990-an, dengan beberapa pemeliharaan dekade terakhir. Kode sumber dapat diunduh dan dikompilasi secara bebas di Linux / Windows menggunakan kompiler Fortran. Ada sumber daya online dan buku tersedia.

Kriging adalah salah satu kekuatan perangkat lunak:

  • 1, 2 atau 3-D kriging kisi, validasi silang, jackknifing
  • SK, OK, Inggris, kriging dengan penyimpangan eksternal
  • kokriging
  • indikator kriging
Mike T
sumber