Сглаживание кривой

Привет!

Совсем замучился с проблемой )) Помогите, пожалуйста!

Есть двумерный очень большой массив точек - координаты (x,y) мест. В целом они идут по "красивой" кривой, но иногда возникают "клубки" (как на картинке - черным цветом).

Я не понимаю, какой алгоритм применять для того, чтобы "клубки" убрать и просто продолжить мою кривую (на картинке - красная линия).

Очень надеюсь на подсказку.

Миниатюры

 

0

Сообщение от Excalibur921 Придумайте данные чтобы ломанная была с клубками как в вашей картинке в первом посте. Атангенс вроде нельзя брать отбалды, там 4 четверти поэтому на графиках бред).

Мне не нужно придумывать данные, они мне даны, и мне необходимо их сгладить. У меня там дальше есть модифицированный график углов - первая картинка, действительно, была очень неправильная ))

Сообщение от Igor3D Не понял какой график чему соответствует. Но в любом случае что-то не так, не может угол так резко прыгать если первый вес 2. Исправьте и покрутите весы, попробуйте 3 и даже 4 (вместо 2)

Первый график - исходные данные, второй - полученные. Я, действительно, что-то перемудрил, по факту второй график совпадает полностью с первым, и изменение w0 в любую сторону во сколько угодно раз вообще ничего не дает (картинка прилагается).

Вот мой код:
У меня есть подозрение, что считать следующую точку, прибавляя к предыдущей просто скорость, не совсем то, что надо (координата + скорость? разве это соразмерные величины?). Хочется на что-нибудь умножить.

Миниатюры

   

0

Сообщение от spirart У меня есть подозрение, что считать следующую точку, прибавляя к предыдущей просто скорость, не совсем то, что надо (координата + скорость? разве это соразмерные величины?). Хочется на что-нибудь умножить.

На шаг по времени = 1. Когда скорость вычисляли на время не делили, а здесь не множим - время сократилось

Сообщение от spirart w[i] <- cos.angle(c(v$x[i-1], v$y[i-1]), c(v$x[i], v$y[i]))

А что это за косинус такой? Может он у Вас 100 и больше, тогда и выходит что ничего не меняется. Попечатайте его. Он должен вычисляться как скалярное произведение нормированных скоростей, т.е это косинус угла между отрезками
И здесь

Сообщение от spirart #смешиваю скорости v$x[i-1] <- (v$x[i-1]*w[1] + v$x[i]*w[i]) / (w[1] + w[i]) v$y[i-1] <- (v$y[i-1]*w[1] + v$y[i]*w[i]) / (w[1] + w[i])

Присваивать это [i - 1] ячейке незачем. И зачем хранить эту таблицу, нужен только накапливаемый вектор и текущий, т.е. 2 переменные

В общем, отладку никто не отменял - на любом языке программирования

0

Сообщение от Igor3D На шаг по времени = 1. Когда скорость вычисляли на время не делили, а здесь не множим - время сократилось

Да, я послал вопрос и почти сразу понял, почему так ))

Сообщение от Igor3D А что это за косинус такой? Может он у Вас 100 и больше, тогда и выходит что ничего не меняется.

Косинус угла между отрезками - ровно то, что Вы и написали, просто вынесено в отдельную функцию.
И возвращает ровно косинус угла между векторами, пример вывода:

Сообщение от Igor3D И зачем хранить эту таблицу, нужен только накапливаемый вектор и текущий, т.е. 2 переменные

Да, может, таблица и не нужна ))

Я пытаюсь сделать следующее: загружаю исходные точки из файла (по сути - это координаты от gps-приемника), убираю клубки, записываю новый, сглаженный массив: gps-приемник дает координаты с большой ошибкой, когда скорость его = 0 (получаются такие клубки), при движении - все в пределах допустимого. Хочу рисовать красивый маршрут на карте.

0

Не знаю что у вас тут за формулы и зачем они.
Если сделать как я говорил выше то будет примерно следующее:

0

Сообщение от spirart Я пытаюсь сделать следующее: загружаю исходные точки из файла (по сути - это координаты от gps-приемника), убираю клубки, записываю новый, сглаженный массив:

Я понял, но с отладкой помочь не могу, т.к. Вашего "козырного" языка никогда не видел На С/С++ все это делается за минуты (если графики под рукой), а тут как-то тяжеловато идет.

Сообщение от spirart gps-приемник дает координаты с большой ошибкой, когда скорость его = 0 (получаются такие клубки), при движении - все в пределах допустимого. Хочу рисовать красивый маршрут на карте

Тогда есть резон осредняться "по длине" как я писал выше. Только в качестве критерия лучше брать угол. Схема:

- для каждой точки просматриваем массив вперед и назад на заданную длину. Попавшие в диапазон точки принимаются для осреднения с весом 1, крайние - с частичным, пропорционально захваченной длине. Имеем для точки среднее значение.

- вычисляем сумму квадратов отклонений углов (а не позиций) и используем ее как вес для взвешивания между исходной и осредненной. Для крайних точек (не вошедших в диапазон) используем их частичные веса

Добавлено через 2 минуты

Сообщение от Excalibur921 Если сделать как я говорил выше то будет примерно следующее:

Да, и я ожидал примерно того же. Но "выкинуть" пики не так уж просто - напр они могут иметь какие-то относительно гладкие участки. Нужна "непрерывная" схема

0

К примеру так Даны 4 точки ABCD. Находим абсолютный угол AB,AC,AD.Если разность углов AB и AC<порог и разность углов AB и AD<порог то точка A принадлежит графику.

Сообщение от spirart это координаты от gps-приемника

С этого нужно было начинать.
ТС может вы вообще не так подходите к вопросу?
Сейчас вы показываете координаты GPS типа X Y координаты на местности, а если построить два графика типа
1) X ось это время Y ось это X координаты
2) X ось это время Y ось это Y координаты
Тогда клубков не будет вообще никогда, появиться просто шум на графиках остановки GPS .В любом случае сильно упростит анализ, и тема эта должна быть очень известная. А значит велик не нужен.
И может такой вопрос задать в теме GPS по типу Метод сглаживания GPS данных?
К примеру: http://habrahabr.ru/post/140274/
И вообще поискать типа Убрать шум на графике, сглаживание GPS данных.

0

Сообщение от Excalibur921 А значит велик не нужен. И может такой вопрос задать в теме GPS по типу Метод сглаживания GPS данных? К примеру: http://habrahabr.ru/post/140274/ И вообще поискать типа Убрать шум на графике, сглаживание GPS данных.

Я тоже рекомендовал "открыть книгу" и все такое Ведь задача, вроде бы, не то чтобы "известная", а даже "известнейшая". Но вот, порыпавшись по ссылкам - что-то не выходит, по крайней мере "сходу". Взять хотя бы ту ссылку что Вы привели - ну не так уж просто. Надо осознать параметры (которых не так уж мало), поэкспериментировать (читай - неск дней), а потом.. да вот хз. Может да, действительно хорошее, грамотное решение. Но может и нет - придется убедиться что Кальман - нормальный мужик, но здесь не катит. Я не "придираюсь", просто везде так же - да, есть, но (боже мой) как "обильно" (громоздко)

В общем, такие ссылки хорошо сливать другим, но не юзать самому Что впрочем совершенно нормально - автор не обязан следовать собственным рекомендациям - они ведь пишутся для других В конце-концов велик - двигатель прогресса, а если попробовать легко - так чего ж нет.

0

Сообщение от Igor3D что-то не выходит, по крайней мере "сходу"

Эта ссылка просто пример =). А вообще направление поиска должен выбрать ТС.
А эта проблема была всегда… когда хочешь что-то повторить, вроде ну вот же 10 листов описание 30 промежуточных формул ( зачем? в итоге 3 для практики)а как понять что куда и как =). Такое чувство, что когда математический уровень программиста проходит какой-то порог, он перестает адекватно писать формулы простыми обозначениями и простым языком пояснять расчет.. И в других форумах бегло глянул вроде как упоминают Калмана для этой задачи наверно ж не зря, упоминают же там GPS.
http://delirium-00.livejournal... tml?page=1

Сообщение от Igor3D Надо осознать параметры (которых не так уж мало),

А можно не курить суть а сделать параметры слайдерами и играть в реальном времени смотреть какие оставить значения настроек константами а алгоритм Калмана пусть будет черным ящиком.
Или

Сообщение от Excalibur921 1) X ось это время Y ось это X координаты 2) X ось это время Y ось это Y координаты

И к таким графикам уже применять:

Сообщение от Excalibur921 К примеру так Даны 4 точки ABCD. Находим абсолютный угол AB,AC,AD.Если разность углов AB и AC<порог и разность углов AB и AD<порог то точка A принадлежит графику.

Хотя это суть как и здесь почти, только в обе стороны:

Сообщение от Igor3D для каждой точки просматриваем массив вперед и назад на заданную длину.

А может и так хватит и углов нет. https://ru.wikipedia.org/wiki/Скользящая_средняя
“Нашёл что это как правило МНК, Калмен, Сплайн, ну и нечто логическое, аля статик навигейшн с СИРФ3... Но это всё не подтвержденные/общие данные.”
Похоже каждый делает как может(вариантов очень много) а коммерческие наверно секретны. На одном форуме говорят сплайн Акимы, создавал давно вопрос этот и глухо.

А если не записывать координаты когда автомобиль остановлен? И все…И ничего не нужно сглаживать =).

0