Barometrická altimetrie: principy, chyby a kompenzace v reálném čase

Proč barometrická altimetrie a kde selhává bez kompenzace

Barometrický výškoměr patří mezi nejenergeticky efektivní a rychlé senzory výšky pro drony. Tlak vzduchu klesá se stoupající nadmořskou výškou, proto z měření statického tlaku lze odvodit výšku nad referenční hladinou. Problémem je, že tlak se mění i časem (počasí) a lokálně (teplotní gradienty, ohřev elektroniky, proudění od vrtulí). Bez inteligentní kompenzace teploty a driftu se barometrická výška dokáže rozchýlit o desítky metrů během minut. Tento článek nabízí systémový návod, jak realizovat kompenzaci v reálném čase v praxi – od fyziky, přes umístění senzoru a filtraci, po fúzi s IMU/GNSS a online odhad zaujatosti (biasu).

Fyzikální základ: z tlaku na výšku

Hypsometrická rovnice (ISA, troposféra): při standardním teplotním gradientu L ≈ 0,0065 K/m platí:
h = (T₀/L) · [1 − (p/p₀)^R·L/g],
kde R je plynová konstanta pro suchý vzduch a g tíhové zrychlení.
Lokální linearizace: blízko pracovního bodu platí Δh ≈ −(RT/g) · (Δp/p). Prakticky: při hladině moře ~1 hPa ≈ 8–9 m, s rostoucí výškou klesá citlivost (Pa/m).
Teplota: výška závisí na teplotě vzduchu ve sloupci, nikoli jen na teplotě čipu. Nesprávné teplotní předpoklady vedou ke systematické chybě.

Zdroj chyb: z čeho vzniká drift

Počasí: barický trend (±2–3 hPa za několik hodin) → desítky metrů chyby bez opětovného referencování (QNH/QFE).
Teplota senzoru: samoohřívání a ohřev od elektroniky mění offset tlaku i jeho citlivost (typicky desítky Pa na 10 °C).
Dynamický tlak a proudění: „prop wash“, boční vítr, ventilace trupu → lokální podtlaky/nadtlaky.
Vibrace: mikromechanické šumy, akustické vlivy vrtulí → vysokofrekvenční šum.
IMU/GNSS chyby při fúzi: pokud se baro „potáhne“ za driftující vertikální rychlostí z IMU, vznikne pomalý posun.

Umístění a mechanika: statický tlak bez rušení

Statický port: pokud je možné, používejte externí statický otvor/komoru s labyrintem a filtrační membránou; vyhýbejte se přímým proudnicím od vrtulí.
Termální izolace: senzor oddělte od horkých komponent (regulátory, procesor) a stíněním snižujte teplotní gradienty.
Anti-vibrace: měkké uložení PCB a krátké hadičky/komory snižují akustické buzení.
Redundance: dva různé baro senzory v jiné pozici (interní + externí) pomáhají detekovat lokální chyby.

Teplotní kompenzace: od výrobních polynomů po online adaptaci

Výrobní koeficienty: moderní senzory mají interní kalibrační konstanty (offset, citlivost). Použijte je, ale počítejte s reziduem.
Vlastní model: v teplotní komoře změřte p_mer vs. T při referenčním tlaku a fitujte polynom 2.–3. řádu pro offset(T) a scale(T). Uložte koeficienty do FW.
Online korekce: během letu sledujte pomalý drift δp(T,t) pomocí inovační sekvence EKF (baro – predikce z IMU/GNSS) a adaptujte offset s dlouhou časovou konstantou (např. 100–300 s).
Warm-up fáze: prvních 60–120 s po zapnutí je kritických. Aplikujte zvýšené filtrování a do not trust flag pro autopilota až do teplotní stabilizace.

Filtrační řetězec: od surového tlaku po „flight-grade“ výšku

Oversampling/decimace: čtěte senzor rychle (100–200 Hz), zprůměrujte na 50–100 Hz, snižujete kvantizační šum.
IIR nízkoprůchodný filtr: α v rozsahu 0,05–0,2 podle dynamiky platformy (prudké manévry vyžadují vyšší pásmo).
Notch/median filtr: potlačení prop-harmonických špiček (např. 80–200 Hz) nebo median z 3–5 vzorků pro „spikes“.
Detektor proudění: pokud je přítomen diferenční tlak nebo zvýšené vibrace IMU, dočasně snižte důvěru baro měření (zvyšte R v EKF).

Fúze senzorů: EKF/komplementární přístup

Krátkodobá stabilita: IMU (svislé zrychlení integrované na rychlost/výšku) s vysokofrekvenčním pásmem.
Střednědobá reference: barometrická výška drží pomalé změny (sekundy až minuty), ale je náchylná na barické trendy.
Dlouhodobě: GNSS (geodetická výška) nebo radar/laser altimetr (nad zemí) slouží jako „anchor“ proti barickému driftu.
Odhad biasu: rozšířený KF obsahuje stav b_baro (zaujatost). Jeho derivaci svazujte na nulovou průměrnou vertikální rychlost v časových oknech („zero-mean constraint“).
Zero-velocity updates (ZUPT): při stabilním visení z Vision/IMU (vertikální rychlost ~0) korigujte baro bias tak, aby ∂h/∂t → 0.

QNH, QFE a reference v autonomii dronů

QNH: tlak redukovaný na hladinu moře → absolutní nadmořská výška. Vhodné pro konzistentní geo-referenci.
QFE: tlak na místě startu → výška nad bodem vzletu (h = 0 při startu). Praktické pro bezpečnostní výšky v lokální misi.
QNE (1013,25 hPa): standardní reference pro letové hladiny – méně běžná v nízkých UAS operacích.
Praxe: při autonomii používejte QFE pro relativní nadmořskou výšku a GNSS/DEM pro absolutní kontrolu nad terénem.

Reálné-časové řízení driftu: praktický algoritmus

Inicializace: 10–20 s průměrování tlaku v klidu, nastavit p_ref a T_ref; pokud je dostupné QNH/QFE ze země, použít.
Teplotní korekce: aplikovat polynom offset(T) a scale(T) na surové měření.
Filtrace: projít přes median filtr (pro odstranění spike) → IIR LPF → downsample na 50 Hz.
EKF fúze: stavy: h, v_z, b_baro. Měření: baro h, IMU a_z, GNSS h (nízká frekvence), laser/radar h_AGL (pokud je).
Adaptivní R: zvyšovat měřicí kovarianci baro při překročení prahu RMS vibrací IMU nebo Δp′.
Anchor aktualizace: každých 30–120 s porovnat trend baro vs. GNSS/AGL; přiřadit malou časovou konstantu k db/dt pro vyrovnání dlouhodobého driftu.
ZUPT/hover: pokud |v_z| < prah a držení polohy je aktivní, přibližovat b_baro k hodnotě, která nulová vertikální drift.

Šum a stabilita: jak nastavit časové konstanty

Allanova analýza: odhadněte bílý šum a bias nestabilitu senzoru. Podle ní nastavte LPF a Q/R v EKF.
Časová konstanta biasu: typicky 100–300 s (elektrické drony), rychlejší při výrazném teplotním gradientu, pomalejší v stabilní troposféře.
Zisk pro GNSS anchor: malý (např. 0,001–0,01 s⁻¹), aby krátkodobé výkyvy GNSS neovládly baro.

Kompenzace při manévrech a vliv dynamického tlaku

Detekce režimu: při prudkých změnách plynu nebo náklonu rostou chyby baro. V autopilotovi během „aggressive“ módu letu snižte váhu baro.
Model prop-wash: pokud znáte RPM/tah, lze aditivně korigovat baro o funkci f(RPM). Jednodušší je robustní filtrování a vhodné umístění portu.

Fúze s AGL (laser/radar) a DEM

AGL senzor: výborný do ~30–50 m nad zemí; nad touto hranicí jeho šum a spotřeba narůstají – baro přebírá roli.
Terrain following: kombinujte baro (absolutní trend) + AGL (relativní reference) + DEM (digitální model reliéfu) pro bezpečné clearance nad terénem.

Redundance a validace měření

Dual-baro voting: pokud je rozdíl > prah (např. 1,5–2,5 hPa), označte chybu a snižte váhu problematického kanálu.
Reasonableness checks: odvozená vertikální rychlost z baro by neměla dlouhodobě odporovat IMU/GNSS (testy konzistence v EKF).

Provozní praxe: kalibrace a misijní návyky

Před letem: reset QFE/QNH, krátké zahřátí (≥60 s), kontrola stability měření (σ < prah).
Během letu: pokud je dostupný pozemní barograf nebo meteostanice, přenášejte trend (telemetrie) a aktualizujte referenci dlouhým filtrem.
Po letu: porovnejte baro „h_landing“ s nulou; zaznamenejte drift biasu pro účely adaptace koeficientů.

Bezpečnostní limity a fallbacky

Geofence ve výšce: absolutní max/min výšky navázané na GNSS a AGL, ne pouze na baro.
RTH profil: návratová trajektorie používá výšku s „anchorem“ (GNSS/AGL) a konzervativní rezervu nad terénem.
Degradované módy: při baro chybě přepnout na GNSS+IMU s oznámenou degradací přesnosti výšky.

Metodika testování a validace

Bench test: stabilní tlak (uzavřená komora), teplotní cyklus 0–50 °C, měření offsetu(T) a hystereze.
Prop-wash test: spuštění motorů na stojanu, vyhodnocení spektra rušení a optimalizace filtru/notch.
Field test – statika: dlouhý hover (5–10 min) s AGL; sledování driftu baro po aplikaci ZUPT korekcí.
Field test – stoupání/klesání: kalibrace zpoždění a přefiltrování (fázové zpoždění) vůči referenci (AGL/RTK).
Metoo test: paralelní záznam baro ze dvou dronů, stejný profil letu, porovnání trendů a odezvy na frontu.

Implementační tipy pro firmware a autopilota

Časové značky: všechna měření časujte na společnou časovou bázi; kompenzujte latence filtrů (predikcí stavu).
Stavový stroj důvěry: flagy warming, valid, degraded, fault; řízení vah v EKF podle stavu.
Logování: ukládejte surový tlak, T, filtrovaný tlak, inovační termín EKF, odhad biasu, váhy měření; usnadňuje ladění.
Konfigurovatelnost: oddělené parametry pro LPF cut-off, zisk anchoru, prahy vibrací, časové konstanty biasu dle typu platformy (MR/FW/VTOL).

Co očekávat: realistická přesnost po kompenzaci

Krátkodobě (sekundy): šum ±0,1–0,3 m u kvalitního senzoru a vhodného LPF.
Střednědobě (minuty): drift ≤ ±1–2 m za stabilního počasí s aktivním bias trackingem.
Dlouhodobě (desítky minut až hodiny): bez anchoru (GNSS/AGL/Q