Москва: проверенные б/у ноутбуки с гарантией и проверкой перед оплатой
+7 906 722-94-45

Главная → Статьи → Постобработка ГНСС и дрона в RTKLIB (PPK) на ноутбуке: файлы RINEX 3.04, метки .MRK и фазовый центр антенны

Постобработка ГНСС и аэрофотосъёмки БПЛА в RTKLIB-EX (PPK) на ноутбуке: от сырых UBX/RINEX до фиксированного решения Q=1

Как геодезисту выполнять кинематическую постобработку PPK в RTKLIB / RTKLIB-EX на полевом ноутбуке без сотовой связи в тайге и карьере: конвертация RINEX 3.04, связка меток затвора .MRK, калибровка фазового центра .ATX и отбраковка Float Q=2.

Постобработка ГНСС и аэрофотосъёмки БПЛА в RTKLIB-EX (PPK) на ноутбуке: от сырых UBX/RINEX до фиксированного решения Q=1

Режим реального времени (RTK по радиоканалу УКВ или модему NTRIP) прекрасен в городе, но в глубоком рудном карьере, в тайге на трассе трубопровода или при высотном облёте крупного массива леса радиолинк неизбежно рвётся за гребнем сопки. В такие моменты на выручку приходит кинематика с постобработкой (PPK — Post-Processed Kinematic): и базовая станция на штативе, и подвижный ровер (вешка или геодезический БПЛА) автономно пишут сырые фазовые измерения несущей на флеш-память с частотой 5–20 Гц, а вечером на полевом ноутбуке камеральщик сводит оба потока вперёд и назад по времени (Forward + Backward Filter), вытягивая сантиметровое фиксированное решение даже там, где обычный RTK-канал через NTRIP ушёл бы в срыв.

Золотым бесплатным ядром такой обработки во всём мире выступает пакет RTKLIB — а именно его модернизированная многочастотная ветка RTKLIB-EX (прямое продолжение знаменитой сборки demo5, переехавшее в основную ветку rtklibexplorer/RTKLIB, на базе которого также работает дружелюбная оболочка Emlid Studio).

Из каких файлов собирается задача PPK на диске ноутбука

Чтобы не запутаться в десятке однотипных расширений с карты памяти дрона и приёмника, придерживайтесь строгой структуры каталога на каждый день изысканий:

Файл входного набора Откуда берётся За что отвечает в расчёте RTKPOST / rnx2rtkp Критическое условие успеха
Ровер *.OBS (RINEX 3.04) Конвертация сырого лога дрона/ровера (*.UBX, *.BIN, *.RTCM3, *.CNB) в RTKCONV Содержит сырые псевдодальности и фазы несущих L1/L2/L5 (GPS, ГЛОНАСС, Galileo, BeiDou) Частота записи на БПЛА — не ниже 5 Гц (при полёте 12 м/с шаг 1 Гц равен 12 метрам интерполяции!)
База *.OBS + *.NAV (RINEX 3.04) Локальная база на пункте ГГС/ОМС либо скачанный суточный/часовой лог ближайшей референц-станции Вычитает ионосферные/тропосферные задержки и ошибки часов спутников; .NAV даёт точные эфемериды орбит Интервал записи базы 1 с (1 Гц) должен полностью накрывать время полёта с запасом +10 минут до взлёта и после посадки!
Журнал меток камеры Timestamp.MRK Папка полёта геодезического БПЛА (DJI Mavic 3E / Matrice) Хранит микросекундный момент срабатывания шторки затвора в шкале времени GPS Week + Tow и плечи смещения Lever Arm Без интерполяции по .MRK координаты EXIF фото имеют полусекундный сдвиг (5–7 м ошибки!)
Калибровка антенн igs20.atx Официальный каталог IGS (ANTEX) Компенсирует смещение электрического фазового центра антенны базы (PCO / PCV) от точки измерения высоты штатива (ARP) Перепутанный отсчёт высоты до низа крепления вместо бампера антенны даёт систематический промах по высоте H на 6–14 см!

Четыре настройки Options в RTKPOST, которые решают судьбу доли Fix (Q=1)

Открыв RTKPOST (или консольный пакетный решатель rnx2rtkp для обсчёта сразу десятка вылетов скриптом), проверьте четыре ключевые вкладки перед нажатием Execute:

  1. Setting 1 -> Filter Type = Combined (или Combined - no phase reset): решатель прогоняет фильтр Калмана сначала от начала полёта к посадке (Forward), затем в обратную сторону от конца к старту (Backward) и сшивает оба вектора! Если в середине залесённого маршрута ровер на 15 секунд зацепил кроны сосен и потерял часть созвездия, обратный проход «подтягивает» целые неоднозначности со стороны открытой поляны — доля Q=1 (Fix) вырастает с 78 % до 98,5–100 %!
  2. Setting 1 -> Frequencies = L1+L2 (или L1+L2+L5 для современных плат) + все четыре глобальные системы (GPS + GLO + Galileo + BeiDou): именно форк RTKLIB-EX / demo5 корректно разрешает целочисленные фазовые неоднозначности смешанных группировок и межчастотные задержки GLONASS FDMA на разных марках приёмников базы и ровера (опция GLONASS AR = Fix and Hold / Autocal).
  3. Setting 1 -> Elevation Mask = 15° и фильтр SNR: спутники, висящие низко над горизонтом (5–10°), проходят через максимальную толщу приземной тропосферы и дают жесточайшее переотражение от борта карьера, крыши автомобиля и металлического забора (Multipath). Поднятие маски возвышения до 15° (а в глубоком залесённом коридоре до 18–20°) парадоксальным образом не ухудшает, а резко стабилизирует сходимость Fix!
  4. Проверка опорных координат базы (Positions -> Base Station): обязательно убедитесь, в какой именно системе координат и эпохе заданы координаты вашей базовой точки (ITRF / WGS-84 против ГСК-2011), прежде чем переводить полученные центры фотографирования в местную систему координат МСК региона и считать итоговую среднюю квадратическую погрешность Mt для межевого или технического плана!

Как по графику RTKPLOT на ноутбуке отличить честный сантиметр от «ложного фикса»

Никогда не экспортируйте таблицу центров снимков в фотограмметрию вслепую — нажмите кнопку Plot и посмотрите цветовую дорожку траектории:

  • Зелёные точки (Q=1, Fix): целочисленные фазовые неоднозначности успешно разрешены, типичное СКО SDN / SDE составляет 0,008–0,018 м в плане и 0,015–0,030 м по высоте.
  • Жёлтые точки (Q=2, Float): плавающее решение с реальным дрейфом 20–80 см, а на коротких срывах — до метра! Если среди 600 кадров облёта 25 кадров оказались Q=2, ни в коем случае не назначайте им в Agisoft Metashape / RealityCapture вес точности 2 см — укажите в колонке точности их реальное СКО из файла events.pos (0,5–1,0 м), чтобы уравнивание связующих точек не выгнуло цифровую модель рельефа «седлом»!
  • Всегда закладывайте на земле минимум 3–5 независимых контрольных опознаков (Check Points), снятых вешкой и не участвующих в уравнивании блока — это единственная честная защита геодезиста от забытого на 10 см удлинителя вехи базовой станции!

Требования к экспедиционному ноутбуку камеральщика

Сам расчёт траектории rnx2rtkp на часовом полёте занимает на современном процессоре 5–25 секунд, но сразу следом на том же ноутбуке в гостинице или вахтовом балке проверяется экспресс-ортофотоплан и плотность перекрытия снимков до отъезда бригады за 400 км с объекта:

  • Лёгкий и крепкий корпус 14–16" с зарядкой USB-C Power Delivery 65–100 Вт от бортовой сети УАЗа / пикапа, встроенным слотом карт SD / microSD, 16–32 ГБ ОЗУ и NVMe SSD от 1 ТБ (сырые DNG/JPEG + RINEX за двухнедельный полевой выезд занимают 300–700 ГБ — подбор выносливых мобильных платформ см. в гиде по ноутбукам для командировок и полевой работы).

Чек-лист полевой бригады перед выключением базовой станции

  • [ ] Запись сырых логов на базе включена минимум за 10 минут до взлёта БПЛА/старта ровера и остановлена не ранее чем через 5–10 минут после посадки при чистом открытом небе.
  • [ ] В полевой журнал с фотографией рулетки записана точная наклонная/вертикальная высота антенны базы с точностью до 1 мм и способ измерения (до низа адаптера ARP или до метки на бампере).
  • [ ] Прямо на объекте до сбора штатива файлы базы и ровера слиты на ноутбук, прогнаны через RTKCONV + RTKPOST, получен Q=1 > 95 % и проверено отсутствие пропусков карточек памяти!

Частые вопросы

До какой максимальной длины базовой линии от базы до дрона работает PPK без потери точности? На двухчастотных приёмниках L1/L2 в спокойной ионосфере невязка линейно растёт примерно на 1 мм на каждый километр удаления (1 см + 1 ppm): для инженерно-топографических планов масштаба 1:500 старайтесь ставить локальную базу не далее 5–10 км от участка съёмки (на 25–35 км дифференциальная ионосфера начинает резко снижать процент Fix в полуденные часы солнечной активности!).

Почему количество меток вспышки в .MRK оказалось на 2 штуки меньше числа JPEG-фотографий на флешке дрона? Классическая ситуация при тестовом ручном щелчке затвором на земле до набора спутников или сбое записи крайнего кадра при выключении питания. Сопоставляйте кадры и строки .MRK не тупо по порядковому номеру с единицы, а с контролем дельты интервалов времени между соседними экспозициями (Δt), как это делают профильные скрипты сшивки PPK.