Я извлекаю данные о высоте маршрута из цифровой модели высоты, используя
my.elev <- over(new.points, mygrid)
- new.points — это объект SpatialPoints с координатами (долгота/широта) около 7000 точек, преобразованными в CRS mygrid.
- mygrid — это SpatialGridDataFrame с более чем 8 миллионами элементов
(подробнее в моем предыдущем вопросе)
Имея несколько значений NA в my.elev, я отладил свой код и обнаружил, что почти все точки в new.points повторяются более одного раза (в моем маршруте несколько сегментов пересекаются два раза):
- первое вхождение имеет соответствующее значение my.elev правильное
- второе имеет значение NA (или иногда совсем другое значение)
Я могу легко решить проблему, устранив повторяющиеся значения в new.points, но мне интересно, почему функция over не возвращает одно и то же значение для одной и той же точки.