Saya sangat baru dalam menggunakan data GIS dan hanya sedikit berpengalaman dengan R. Saya telah membaca tentang bagaimana menganalisis data spasial menggunakan buku PDF spatial-analyst.net, jadi saya tidak sepenuhnya hilang, tapi saya pikir saya bisa menggambarkan masalah saya dan orang mungkin menyarankan ide.
Saya memiliki dataset dengan sekitar 2000 pengukuran pada koordinat lat / long yang berbeda, meskipun saya mungkin akan membagi dataset ini karena data dikumpulkan selama 3 tahun dan kondisinya berubah dari waktu ke waktu. Mari kita sebut variabel yang diukur "IP."
Saya ingin membuat peta IP di area penuh yang dipertanyakan menggunakan Kriging atau metode interpolasi lainnya pada data sampel. Lalu saya ingin membuat histogram yang mengukur jumlah lahan di berbagai ember IP. Saya juga perlu membuat histogram yang menunjukkan jumlah sampel dalam setiap ember (perhatikan sampel bisa memiliki IP aktual yang lebih tinggi atau lebih rendah dari apa yang diprediksi kriging untuk lahannya).
Saya mengikuti cara memuat data ke dalam SpatialPointsDataFrame dan menjalankan analisis kriging, di mana saya mengalami masalah adalah bagaimana mengkonversi data itu ke dalam kerangka data yang disambungkan sehingga saya dapat melakukan analisis histogram.
Adakah saran untuk mengubah poin menjadi kisi?
pts
ada dalam contoh saya di atas. Cukup jalankan kode pada objek SpatialPointsDataFrame Anda!r
) untuk digunakan dalamhist
prosedur yang mirip, yang hanya masalah ekspresihist(getValues(r))
.Paket plotKML memiliki fungsi yang disebut
vect2rast
. Fungsi ini pada dasarnya memperluasrasterize
fungsi yang tersedia dalam paket raster. Keuntungan darivect2rast
; Namun, adalah bahwa ia tidak memerlukan input dari sisi pengguna yaitu secara otomatis menentukan ukuran sel grid dan kotak pembatas berdasarkan pada properti dari set data input. Ukuran sel kisi diperkirakan berdasarkan kepadatan / ukuran fitur di peta (nndist
fungsi dalam paket spatstat).sumber