Я должен рассчитать оценку плотности ядра (kde) из списка координат широты и долготы. Но один градус по широте - это не то же самое расстояние, что и один градус по долготе, это означает, что отдельные ядра будут овальными, особенно если точка дальше от экватора.
В моем случае все точки достаточно близки друг к другу, поэтому их преобразование в плоскую землю не должно вызывать многих проблем. Однако мне все еще интересно, как это должно быть правильно обработано на случай, если это не так.
spatial
kernel-smoothing
geostatistics
geography
Аарон де Виндт
источник
источник
Ответы:
Вы могли бы рассмотреть использование ядра, особенно подходящего для сферы, такого как плотность фон Мизеса-Фишера
где и - местоположения на единичной сфере, выраженные в трехмерных декартовых координатах.μ x
Аналогом пропускной способности является параметр . Вклад в местоположение из точки входа в местоположение на сфере, имеющий вес , поэтому равенκ x μ ω(μ)
Для каждого суммируйте эти вклады по всем входным точкам .x μi
Чтобы проиллюстрировать это, вотμi ω(μi) κ 6
R
код для вычисления плотности фон Мизеса-Фишера, генерации некоторых случайных местоположений и weights (12 из них в коде), и отображения карты результирующей плотности ядра для указанного значение (равно в коде).Точки показаны черными точками, размеры которых пропорциональны их весам . Вклад большой точки вблизи очевиден во всех северных широтах. Яркое желто-белое пятно вокруг него будет приблизительно круглым, если показано в подходящей проекции, такой как Орфографическая (земля из космоса).μi ω(μi) (100,60)
источник