Bagaimana Anda menggunakan v.kernel GRASS?

12

Saya bingung bagaimana cara menggunakan v.kernel GRASS.

Saya memiliki lapisan vektor sekitar 2,5 juta poin. Saya ingin membuat peta panas menggunakan v.kernel untuk menunjukkan konsentrasi, karena saya memiliki instance variabel dengan titik-titik yang tumpang tindih, kadang-kadang tumpang tindih besar.

Saya sudah mendapatkan layer vektor ini dalam GRASS, dan ini akan menampilkannya dengan baik.

Saya sudah mencoba menggunakan perintah v.kernel GRASS berdasarkan apa yang saya lihat di sini dan di forum lain, dan saya tidak bisa membuatnya melakukan apa pun selain menghasilkan raster yang hanya berupa kotak merah muda.

Inilah perintah yang saya gunakan:

v.kernel --verbose input=master_grass7 output=master_grass7a_heatmap stddeviation=.0001

Saya telah memvariasikan stddeviasi untuk semua jenis nilai mulai dari 10.000 hingga 0,000001, dan itu tidak berpengaruh.

Saya sudah membaca dokumentasi v.kernel berulang kali dan tidak begitu mengerti apa maksudnya. Setidaknya, instruksinya ada pada konsep esoteris, tidak ada yang praktis. Saya juga sudah memeriksa kode sumbernya , dan saya juga tidak terlalu memahaminya. Ya, saya bisa membaca C. Masalahnya tergantung pada banyak hal yang didefinisikan di GRIS GIS.

Saya juga telah melakukan banyak pencarian Google, dan saya tidak dapat menemukan panduan yang komprehensif. Semua yang saya dapatkan adalah salinan halaman doc / man v.kernel yang tersebar atau orang-orang yang tampaknya dapat bekerja tanpa repot.

Saya juga telah memeriksa konsep estimasi kepadatan kernel (KDE), dan bahkan kemudian saya tidak melihat cara menggunakan perintah v.kernel. Perintah itu tampaknya merupakan interpretasi khusus dari KDE; sakelarnya tampaknya tidak sesuai dengan konsep KDE umum.

Jadi kembali ke pertanyaan utama di sini: bagaimana mungkin seseorang yang tidak akrab dengan pengembangan produk GRASS menggunakan perintah v.kernel? Apakah tersedia terjemahan bahasa biasa?

Aren Cambre
sumber

Jawaban:

9

The v.kernelalgoritma menghitung kepadatan vektor titik untuk setiap sel dari peta raster. Jika Anda sejauh ini hanya menggunakan objek vektor, kemungkinan besar Anda belum mengatur region(yang tidak hanya menentukan luasan, tetapi juga resolusi raster) secara memadai: Anda mungkin memiliki regionset hanya satu baris dan satu kolom, yang berarti v.kernelalgoritma hanya akan menghitung kepadatan kernel sebagai nilai tunggal di seluruh peta. Pengaturan wilayah seperti ini baik-baik saja jika Anda tidak menggunakan peta raster, karena peta vektor tidak peduli dengan pengaturan resolusi. Periksa pengaturan wilayah Anda menggunakan g.region -p, dan jika rowsdan colsdiatur ke 1, tingkatkan resolusi dengan menggunakan

g.region rows=<number of rows> cols=<number of columns>

atau

g.region res=<length>

di mana <length>panjang sel resolusi di unit peta.

Jake
sumber