v.krigegrass - Dalam Talian di Awan

JADUAL:

NAMA

v.krige - Melakukan kriging biasa atau blok untuk peta vektor.

Kata-kata berkaitan

vektor, interpolasi, raster, kriging

SINOPSIS

v.krige
v.krige - membantu
v.krige input=nama ruang=nama [output=nama] [pakej=rentetan]
[model=rentetan[,rentetan,...]] [menyekat=integer] [pelbagai=integer] [nugget=integer]
[ambang=integer] [output_var=nama] [--menimpa] [--membantu] [--kata kerja] [--tenang]
[--ui]

Bendera:
--tulis ganti
Benarkan fail output untuk menulis ganti fail sedia ada

- membantu
Cetak ringkasan penggunaan

--verbose
Keluaran modul verbose

--senyap
Keluaran modul senyap

--ui
Paksa melancarkan dialog GUI

Parameter:
input=nama [diperlukan]
Nama peta vektor input
Nama peta vektor titik yang mengandungi data sampel

ruang=nama [diperlukan]
Nama lajur atribut dengan nilai berangka untuk diinterpolasi

output=nama
Nama untuk peta raster output
Jika ditinggalkan, akan _kriging

pakej=rentetan
Pakej R untuk digunakan
Pilihan: gstat
default: gstat

model=rentetan [, rentetan,...]
Model Variogram
Biarkan kosong untuk menguji semua model (memerlukan automap)
Pilihan: Nug, Exp, Sph, Gau, Exc, tikar, Ste, Cir, Lin, Bes, pen, Per, Hol, Log, Pow,
Spl, kaki, Err, Int

menyekat=integer
Saiz blok (blok persegi)
Saiz blok. Digunakan oleh blok kriging.

pelbagai=integer
Nilai julat
Dibetulkan secara automatik jika tidak ditetapkan

nugget=integer
Nilai nugget
Dibetulkan secara automatik jika tidak ditetapkan

ambang=integer
Nilai ambang
Dibetulkan secara automatik jika tidak ditetapkan

output_var=nama
Nama untuk peta raster varians keluaran
Jika ditinggalkan, akan _kriging.var

DESCRIPTION

v.krige membenarkan melaksanakan operasi Kriging dalam persekitaran GRASS GIS, menggunakan perisian R
fungsi di latar belakang.

NOTA

v.krige hanyalah bahagian hadapan kepada R. Pilihan dan parameter adalah sama yang ditawarkan oleh
pakej automap and gstat.

Kriging, seperti kaedah interpolasi lain, bergantung sepenuhnya pada ciri data input.
Analisis penerokaan data digalakkan untuk mengetahui outlier, trend, anisotropi,
pengagihan yang tidak sekata dan seterusnya memilih algoritma kriging yang akan memberikan yang paling banyak
hasil yang boleh diterima. Pengetahuan yang baik tentang set data adalah lebih berharga daripada ratusan
parameter atau perkakasan berkuasa. Lihat buku Isaaks dan Srivastava, lengkap dan jelas
walaupun agak ketinggalan zaman.

Kebergantungan
R perisian >= 2.x

rpy2
Python mengikat R. Nota! rpy versi 1 tidak disokong.

R pakej automap, gstat, rgrass7 and rgeos.
automap adalah pilihan (menyediakan muat variogram automatik). Pasang pakej melalui R
baris arahan (atau GUI pilihan anda):
install.packages("rgeos", dep=T)
install.packages("gstat", dep=T)
install.packages("rgrass7", dep=T)
install.packages("automap", dep=T)

Nota khususnya Debian GNU / Linux
Pasang tanggungan. Perhatian! python-rpy IS TIDAK SESUAI.:
aptitude install R python-rpy2
Untuk memasang pakej R, gunakan sama ada fungsi R yang disenaraikan di atas (sebagai akar atau sebagai pengguna), sama ada
pakej Debian [5], tambahkan ke senarai repositori untuk 32bit atau 64bit (ambil
talian yang sesuai):
deb http://debian.cran.r-project.org/cran2deb/debian-i386 ujian/
deb http://debian.cran.r-project.org/cran2deb/debian-amd64 ujian/
dan dapatkan pakej melalui aptitude:
pasang aptitude r-cran-gstat r-cran-rgrass7

Nota khususnya Windows
Susun GRASS GIS mengikut panduan ini. Anda juga boleh menggunakan Linux dalam mesin maya. Ataupun
pasang Linux dalam partition berasingan HD. Ini tidak menyakitkan seperti yang kelihatan,
terdapat banyak panduan melalui Internet untuk membantu anda.

Pengiraan masa isu-isu
Sila ambil perhatian bahawa walaupun bilangan titik data input tinggi dan/atau resolusi rantau tinggi
menyumbang kepada keluaran yang lebih baik, kedua-duanya juga akan memperlahankan pengiraan kriging.

CONTOH

Contoh Kriging berdasarkan peta ketinggian (set data Spearfish).

bahagian 1: rawak pensampelan daripada 2000 titik vektor daripada peta ketinggian yang diketahui. Setiap mata akan
terima nilai ketinggian daripada raster ketinggian, seolah-olah ia datang daripada tinjauan titik.
# kurangkan peleraian untuk contoh ini
g.region raster=elevation -p res=100
v.keluaran rawak=rand2k_elev npoints=2000
v.db.addtable map=rand2k_elev columns="elevation double precision"
v.what.rast map=rand2k_elev raster=elevation column=elevation
bahagian 2: mengeluarkan mata kekurangan ketinggian sifat-sifat. Titik sampel di sempadan
peta ketinggian tidak menerima sebarang nilai. v.krige tidak mempunyai tindakan pilihan untuk mengatasi no
nilai data, jadi pengguna mesti menyemaknya dan memutuskan perkara yang perlu dilakukan (buang mata, isi
dengan nilai titik terdekat, isikan dengan min global/tempatan...). Dalam yang berikut
baris kod, titik tanpa data dialih keluar daripada peta.
v.ekstrak input=rand2k_elev output=rand2k_elev_filt where="ketinggian bukan NULL"
Semak hasil baris sebelumnya ("bilangan atribut NULL" mestilah 0):
v.univar map=rand2k_elev_filt type=point column=elevation
bahagian 3: membina semula DEM melalui kriging. Cara paling mudah untuk berlari v.krige daripada CLI sedang menggunakan
muat variogram automatik (nota: memerlukan pakej automap R). Nama peta output adalah pilihan,
modul menciptanya secara automatik menambahkan "_kriging" pada nama peta input dan juga
menyemak untuk menulis ganti. Jika output_var ditentukan, peta varians juga dibuat.
Kesesuaian variogram automatik disediakan oleh automap pakej R. Model variogram yang diuji oleh
fungsi pemasangan ialah: eksponen, sfera, Gaussian, Matern, M.Stein's
penjajaran. Rangkaian model yang lebih luas tersedia daripada pakej gstat dan boleh diuji
pada GUI melalui pemplotan variogram. Jika model dinyatakan dalam CLI, juga ambang,
nilai nugget dan julat perlu disediakan, jika tidak, ralat akan dibangkitkan (lihat second
contoh v.krige perintah).
# muat variogram automatik
v.krige input=rand2k_elev_filt column=elevation \
output=rand2k_elev_kriging output_var=rand2k_elev_kriging_var
# tentukan model variogram, cipta peta varians juga
v.krige input=rand2k_elev_filt column=elevation \
output=rand2k_elev_filt_kriging output_var=rand2k_elev_filt_kriging_var \
model=Tikar ambang=2500 nugget=0 julat=1000
Atau jalankan wxGUI, untuk menyesuaikan variogram secara interaktif dan meneroka pilihan:
v.krige
Mengira ramalan kesilapan:
r.mapcalc "rand2k_elev_kriging_pe = sqrt(rand2k_elev_kriging_var)"
r.univar map=elevation
r.univar map=rand2k_elev_kriging
r.univar map=rand2k_elev_kriging_pe
Keputusan menunjukkan ralat yang tinggi, seperti teknik kriging (kriging biasa dan blok).
tidak dapat mengendalikan set data dengan arah aliran, seperti yang digunakan dalam contoh ini: ketinggian ialah
lebih tinggi di sudut barat daya dan lebih rendah di sudut timur laut. Universal kriging boleh memberi
keputusan yang jauh lebih baik dalam kes ini kerana ia boleh mengendalikan arah aliran. Ia boleh didapati dalam pakej R
gstat dan akan menjadi sebahagian daripada keluaran v.krige akan datang.

Gunakan v.krigegrass dalam talian menggunakan perkhidmatan onworks.net