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