1. Введение

В данном техническом отчёте приводятся следующие результаты:

Первые 3 результата приведены для одного набора данных, для последнего используется другой набор данных. Данные обрабатывались после миссии, чтобы реализовать необходимый сценарий проверки (в том числе эмуляцию работы в реальном времени).

2. Описание траекторий движения для первых трех результатов

Движение происходило в условиях города на легковом автомобиле и длилось от 10 до 30 минут. На большинстве проездов также присутствовала стоянка в начале продолжительностью около 10 минут. Ниже на графиках представлены траектории и профили скорости для всех 8 проездов. Проезды именуются по номеру дня и номеру проезда ("лога") в этот день.

trajectory day 1 log 1
trajectory day 2 log 1
trajectory day 3 log 1
trajectory day 3 log 2
trajectory day 4 log 1
trajectory day 4 log 2
trajectory day 4 log 3
trajectory day 5 log 1

3. Проверка точности определения углов ориентации

Для оценки точности определения углов ориентации системой Acrux было проведено сравнение с решением, полученным с помощью высокоточного инерциального измерительного блока Оптолинк БЧЭ 501.

3.1. Методика сравнения

При сравнении решения двух инерциальных навигационных систем нужно иметь ввиду, что оси этих систем никогда не будут идеально совпадать. Это означает, что разность углов ориентации будет содержать постоянное и неизвестное смещение. Чтобы его оценить и сравнить углы ориентации в абсолютном смысле, нужно провести серию экспериментов при фиксированном закреплении инерциальных измерительных блоков. Именно такой подход использовался в нашем сравнении.

Сравнение было проведено следующим образом:

  1. Решение Acrux было получено в режиме обычной эксплуатации в реальном времени: использовалось 2 антенны и спидометр, автомобильная модель движения, были доступны поправки типа RTCM3 с близлежащей базовой станции

  2. Решение по Оптолинку было получено с помощью стандартного алгоритма объединения измерений ИИБ и позиционных решений (стандартной точности) навигационного приемника (метод “слабой интеграции”). Более никаких данных или ограничений на движение не использовалось. Это позволило получить решение независимое от решения Acrux

  3. Временная синхронизация решения обеспечивалась привязкой измерений Оптолинк к системному времени Acrux

  4. Использовались 8 наборов данных записанных в течении 5 дней

  5. По этим данным были оценены смещения углов ориентации за счет несоосности блоков и далее были проанализированы скомпенсированные разности углов

3.2. Характеристика ИИБ Оптолинк БЧЭ 501

Инерциальный измерительный блок Оптолинк БЧЭ 501 имеет следующие характеристики датчиков:

  • Гироскопы:

    • среднеквадратичное значение смещения нуля при включении - 0.03 °/час

    • шум - 0.003 °/√час

    • нестабильность смещения нуля - 0.003 °/час на временах порядка часа

  • Акселерометры:

    • среднеквадратичное значение смещения нуля при включении - 0.002 м/с2

    • шум - 0.003 м/с/√час

    • нестабильность смещения нуля - 0.001 м/с2 на временах порядка часа

Параметры шума и стабильности были оценены по методике дисперсии Аллана и примерно соответствует заявляемым производителем. Эти параметры использовались в модели ИИБ для получения решения.

Решение, полученное с помощью ИИБ Оптолинк, имеет ожидаемую среднеквадратичную точность не хуже, чем:

  • 0.01° по крену и тангажу

  • 0.05° по курсу

Заявленная точность Acrux же следующая:

  • 0.05° по крену и тангажу

  • 0.15° по курсу

Т. к. точность решения с помощью ИИБ Оптолинк превосходит точность решения Acrux в 3 или более раз, то по принципам метрологии оно может выступать в качестве эталонного. Из этого мы и исходим в дальнейшем.

3.3. Графики разности углов ориентации

Для всех 8 наборов проездов были посчитаны разности углов ориентации между решением Acrux и решением с помощью ИИБ Оптолинк. Для приведения решений к одинаковым моментам времени использовалась интерполяция ориентации типа SLERP.

Ниже представлены графики разности углов ориентации. Черными пунктирными линиями отмечены моменты времени, после которых считается, что решения достигли оптимальной точности (“сошлись”). Эти значения были выбраны в данном случае "на глаз" по виду графиков разности.

attitude difference day 1 log 1
attitude difference day 2 log 1
attitude difference day 3 log 1
attitude difference day 3 log 2
attitude difference day 4 log 1
attitude difference day 4 log 2
attitude difference day 4 log 3
attitude difference day 5 log 1

По этим графикам можно сделать следующие выводы:

  1. Разность стабилизируется по прошествии некоторого времени. Это связано главным образом с тем, что оптимальная точность Acrux достигается в процессе движения

  2. Средняя разность между углами приблизительно совпадает на всех логах

  3. Разности углов тангажа и крена более стабильны, чем разности угла курса. Это объясняется принципом оценки этих углов, оценки тангажа и крена точнее и стабильней в принципе во всех системах инерциальной навигации

3.4. Количественная оценка точности

Для оценки точности Acrux было проделано следующее:

  1. Были посчитаны средние разности углов на каждом логе после того, как решения "сошлись". Смещения для углов крена, тангажа и курса за счет несоосности ИИБ были определены, как средние значений за каждый проезд

  2. Временные серии разностей в каждом проезде были скомпенсированы на вычисленные смещения (путем вычитания)

  3. Были вычислены среднеквадратичные значение разностей после компенсации, как характеристика точности

На графике ниже представлены средние значения разностей углов:

mean difference

Оценки смещения углов за счет несоосности блоков, как средние представленных значений, получились следующими:

  • 0.63° по крену

  • 1.06° по тангажу

  • -0.83° по курсу

Среднеквадратичные значения разностей углов ориентации после компенсации указанных смещений изображены на графике:

rms

3.5. Вывод

На всех проездах среднеквадратичные разности после компенсации несоосности не превосходили:

  • 0.04° по крену и тангажу

  • 0.12° по курсу

А общее среднеквадратичное значение (как среднеквадратичное от значений в каждый проезд) оказалось следующим:

  • 0.032° по крену (заявленная 0.05°)

  • 0.036° по тангажу (заявленная 0.05°)

  • 0.094° по курсу (заявленная 0.15°)

Данные числа свидетельствуют о том, что система Acrux обеспечивала заявленную точность по определению углов ориентации транспортного средства.

4. Проверка точности одометрического решения

Под одометрическим решением понимается решение, полученное с помощью инерциального измерительного блока и одометрического датчика (спидометра, одометра). Прибор выдает одометрическое решение при пропаже сигналов ГНСС (туннели, закрытые парковки и т. п.)

4.1. Методика оценки точности

Ошибка одометрического решения нарастает по мере движения транспортного средства и обычно характеризуется, как процент от пройденного расстояния. Скорость роста ошибки зависит от начальной ошибки ориентации, точности оценки параметров ИИБ и вообще корректности моделирования его ошибок, способа обработки информации.

Проверка была организована следующим образом:

  1. Использовались данные записанные при движении автомобиля в условиях города, которые включали ИИБ, ГНСС, курс, получаемый по двум антеннам, измерения скорости из CAN-шины автомобиля и поправки RTCM3 с близлежащей базовой станции

  2. Эталонная траектория была получена, как решение Acrux с использованием всей доступной информации

  3. Для оценки точности одометрии сигналы ГНСС отключались (включая курс с помощью двух антенн) после прохождения половины времени лога. То есть первую половину времени система работает в обычном режиме и оценивает свои параметры, а далее работает в одометрическом режиме, где используется только ИИБ и измерения скорости

  4. Ошибки одометрии вычисляются, как разности эталонной и одометрической траекторий после отключения сигналов ГНСС

4.2. Графики одометрической ошибки

Для каждого проезда были посчитаны одометрические ошибки, как разности координат одометрической и эталонной траекторий после отключения сигналов ГНСС. Графики начинаются с середины времени лога в момент пропажи ГНСС сигналов и перехода системы в одометрический режим. В этот момент одометрическая ошибка равняется нулю. Также на графиках изображены \$\pm 1 \sigma\$ границы точности решения, выдаваемые навигационным фильтром Acrux.

dropout day 1 log 1
dropout day 2 log 1
dropout day 3 log 1
dropout day 3 log 2
dropout day 4 log 1
dropout day 4 log 2
dropout day 4 log 3
dropout day 5 log 1

По этим графикам можно сделать следующие выводы:

  1. Ошибки укладываются в \$\pm 1 \sigma\$ границы с запасом или не выходят за них существенно. Это говорит о консистентном поведении фильтра и достаточно консервативном подходе к оценке точности решения

  2. Скорость роста вертикальной ошибки существенно меньше скорости роста горизонтальных ошибок, что объясняется более высокой точностью оценок углов крена и тангажа по сравнению с точностью оценки угла курса

4.3. Количественная оценка точности

В качестве характеристики точности были вычислены средние относительные горизонтальные и вертикальные ошибки. Под относительной ошибкой понимается ошибка, делённая на пройденное расстояние после момента пропажи ГНСС сигналов. Горизонтальная ошибка вычисляется, как корень из суммы квадратов ошибок в северном и восточных направлениях. Усреднение производится для расстояний больше 100 метров, дабы исключить случайные колебания относительной ошибки при малых значениях знаменателя (пройденного расстояния).

На графике ниже представлены указанные средние ошибки в %:

odometry performance

4.4. Вывод

На всех проездах одометрические ошибки не превосходили:

  • 0.9% для горизонтальной

  • 0.1% для вертикальной

Медианные значения оказались следующими:

  • 0.17% для горизонтальной (заявленное значение 0.5%)

  • 0.03% для вертикальной

Результаты (количественные и графики) свидетельствуют о достаточно точной и консистентной работе Acrux в режиме одометрии. Это является косвенным подтверждением того, что система получает оценки углов ориентации консистентные с их предсказанной точностью и спецификацией системы.

5. Постобработка: улучшение одометрического решения

Сигналы ГНСС обычно пропадают на некотором ограниченном участке траектории: например, когда автомобиль проезжает туннель. В этом случае, до и сразу после имеется точное позиционное решение.

В режиме постобработки алгоритм оптимизации “сшивает” точки до и после пропажи ГНСС сигналов, используя одометрические измерения, при этом точность решения по сравнению с режимом работы в реальном времени может существенно улучшиться.

5.1. Графики одометрической ошибки

Используется следующий сценарий обработки данных: ГНСС сигналы отключаются на 50% времени лога и включаются вновь на 95% времени лога. Тем самым моделируется ситуация пропажи ГНСС на участке траектории с последующим восстановлением.

Ниже приведены графики позиционных ошибок и стандартные отклонения, выдаваемые алгоритмом, с начала пропажи сигналов ГНСС:

dropout optimization day 1 log 1
dropout optimization day 2 log 1
dropout optimization day 3 log 1
dropout optimization day 3 log 2
dropout optimization day 4 log 1
dropout optimization day 4 log 2
dropout optimization day 4 log 3
dropout optimization day 5 log 1

По этим графикам можно заметить следующее:

  • Согласно графикам СКО, одометрическая ошибка достигает максимума где-то в середине интервала пропажи сигналов ГНСС. Это можно понять интуитивно, рассмотрев процесс интегрирования в обратном по времени направлении, начиная с момента восстановления ГНСС.

  • Использование постобработки (оптимизации) позволяет получить одометрического решение сопоставимое со стандартной точностью ГНСС приемников (2-3 метра) на интервалах времени превышающих 15 минут

  • На логе day_2_log_1 решение после восстановления ГНСС получилось не слишком точным и консистентным. Это связано с непродолжительным отрезком приема ГНСС в конце проезда в неблагоприятными условиях (отсутствие движения и отражения ГНСС сигналов). Поэтому на конечном отрезке решение типа RTK_FIXED не определяется, а решение типа RTK_FLOAT оказалось смещенным.

5.2. Сравнение с решением в режиме фильтра

Для режима фильтра результаты в почти аналогичном сценарии приведены в предыдущем разделе. Стоит отметить, что на данном наборе логов фильтр в одометрии обычно давал результаты существенно лучше средних ожидаемых. Наиболее существенное уменьшение ошибки в режиме постобработки можно продемонстрировать для лога day_5_log_1. Для него графики ошибок с границами СКО приведены ниже:

optimization vs filter odometry

На графиках видно, что ожидаемые пиковые горизонтальные ошибки (СКО) уменьшаются приблизительно в 10 раз, а реальные ошибки (на данном конкретном логе) уменьшилась еще более существенно.

5.3. Вывод

Режим оптимизации в постобработке позволяет кардинально изменить свойства одометрического решения при пропаже ГНСС сигналов и уменьшить его ошибку на порядок или более. Это является одним из наиболее впечатляющих и полезных свойств алгоритма оптимизации для постобработки логов.

6. Постобработка: работа в условиях помех сигналам ГНСС

В этом разделе приводится отчет по использовании постобработки (режим processing_mode: optimization) для получения решения в условиях помех сигналам ГНСС. Они характеризуется в среднем пониженным соотношением сигнал/шум, а также нерегулярными сильными его падениями с существенным или катастрофическим увеличением шума всех измеряемых величин. В некоторые промежутки времени сигналы ГНСС полностью пропадают.

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

6.1. Характеристики принимаемых сигналов ГНСС

Изменение характеристики сигналов ГНСС во времени изображены на графиках ниже:

raw stats 1
raw stats 2
raw stats 3

Для антенны типа ФАР можно отметить, что число сигналов примерно в два раза меньше (т. к. она принимает сигналы только в диапазоне L1), однако при её использовании “просадок” и полных потерь сигнала меньше. Между двумя другими антеннами принципиальной разницы не наблюдается.

В таких условиях использование RTK не имеет смысла, тем более невозможно проверить точность и корректность такого решения. Поэтому постобработка запускалась с настройкой rtcm_mode: off.

6.2. Построение решения на карте

Единственным способом проверить точность позиционного решения в таких условиях является наложение на карту, что было проделано в программе QGIS с использованием слоя Open Street Maps. Ниже приведены три таких изображения со следующими пояснениями:

  • изображен примерно один и тот же пространственный район, в котором условия приема ГНСС были наихудшими

  • точки решения разделены по временем интервалом 2 секунды

  • области полной пропажи ГНСС можно идентифицировать по оранжевым точкам, помеченными, как “Acrux DEAD_RECKONING”

  • изображения можно увеличить по щелчку мыши

map 1
Рис. 1. Проезд 1
map 2
Рис. 2. Проезд 2
map 3
Рис. 3. Проезд 3

Сложность условий приема ГНСС можно понять по хаотическому поведению и пропажи решения внутреннего приемника (красные точки). Также можно увидеть, что решения Acrux в постобработке по крайней мере качественно отражает движения и вписывается в контуры улиц. Также из-за принципа получения решения, оно оказывается гладким и лишенным любых “нефизичных” скачков или перестроений. При детальном анализе поведения решения относительно слоя карты, можно оценить максимальную ошибку решения, как 10 метров.

6.3. Сравнение с решением реального времени

Результаты работы Acrux в реальном времени можно оценить, как условно в 2-3 раза хуже в смысле максимального отклонения позиции. Также решение реального времени неминуемо подвержено скачкообразным коррекциям в моменты восстановления сигналов ГНСС. Чтобы оставить изображения карт и траекторий читаемыми, результаты такого решения не приводятся.

6.4. Вывод

Постобработка навигационных данных Acrux позволяет получить решение приемлемой точности (ошибка менее 10 метров на всей траектории) в условиях помех сигналам ГНСС. Для реализации точной навигации в таких условиях требуется либо использовать высокоточные и дорогостоящие ИИБ и одометр, либо привлечение измерений, основанных на других принципах и подходах (компьютерное зрение, лазерное сканирование). Навигационное решение типа ГНСС/ИНС в этом случае необходимо для получения хорошей начальной оценки геолокации.