Saya harus menghitung estimasi kepadatan kernel 2d (KDE) dari daftar koordinat lintang dan bujur. Tetapi satu derajat dalam garis lintang tidak sama dengan satu derajat dalam garis bujur, ini berarti bahwa masing-masing kernel akan berbentuk oval, khususnya titik selanjutnya dari khatulistiwa.
Dalam kasus saya, semua titik cukup dekat satu sama lain sehingga mengubahnya menjadi tanah datar seharusnya tidak menimbulkan banyak masalah. Namun saya masih penasaran tentang bagaimana ini harus ditangani dengan benar jika ini tidak benar.
spatial
kernel-smoothing
geostatistics
geography
Aaron de Windt
sumber
sumber
Jawaban:
Anda mungkin mempertimbangkan untuk menggunakan kernel yang sangat cocok untuk bola, seperti kepadatan von Mises-Fisher
di mana dan x adalah lokasi pada unit sphere yang dinyatakan dalam koordinat Cartesian 3D.μ x
Analog dari bandwidth adalah parameter . Kontribusi ke lokasi x dari titik input di lokasi μ pada bola, memiliki bobot ω ( μ ) , oleh karena itu adalahκ x μ ω(μ)
Untuk setiap , jumlah kontribusi ini pada semua titik input μ i .x μi
Untuk menggambarkan, di sini adalahμi ω(μi) κ 6
R
kode untuk menghitung kepadatan von Mises-Fisher, menghasilkan beberapa lokasi acak dan bobot ω ( μ i ) (12 dari mereka dalam kode), dan menampilkan peta kepadatan kernel yang dihasilkan untuk nilai tertentu dari κ (sama dengan 6 dalam kode).sumber