13.6.2.2.2. 원시 데이터 기반 GDOP(Geometric Dilution of Precision) 재계산 및 위성 기하학 배열 품질 저하(Satellite Geometry Degradation) 구역 식별

13.6.2.2.2. 원시 데이터 기반 GDOP(Geometric Dilution of Precision) 재계산 및 위성 기하학 배열 품질 저하(Satellite Geometry Degradation) 구역 식별

GNSS(Global Navigation Satellite System)의 신뢰성은 획득한 위성의 절대적인 숫자결정론(Number of Satellites)에 의해서만 좌우되지 않는다. PX4-Autopilot의 확장 칼만 필터(EKF2)는 위성들이 천구 상에 어떻게 퍼져 있는지를 수치화한 **기하학적 정밀도 저하율(GDOP, Geometric Dilution of Precision)**에 극도로 민감하게 반응한다.

도시 협곡(Urban Canyon)이나 깊은 산악 지대에서 비행할 때, ULog에 기록된 위성 수는 12개 이상으로 충분해 보임에도 불구하고 기체가 심각한 위치 이탈(Position Drift)을 겪는 현상의 90% 이상은 바로 이 부분적인 위성 기하학 배열(Satellite Geometry)의 붕괴에서 기인한다. ULog의 vehicle_gps_position_s 토픽을 활용하여 GDOP를 재계산하고 취약 구역을 식별하는 고도화된 프로파일링 기법은 무인기 비행 안전 확보의 필수 불가결한 과정이다.

0.1 GDOP 파생 인자(HDOP, VDOP)의 수학적 본질과 오토파일럿의 수용

GNSS 수신기 칩셋(예: u-blox F9P)은 수신된 위성들의 방위각(Azimuth)과 앙각(Elevation Angle)을 기반으로 공분산 행렬의 대각 성분을 취하여 오차 증폭 계수를 계산한다. 이렇게 산출된 수평적(HDOP, Horizontal DOP) 및 수직적(VDOP, Vertical DOP) 정밀도 저하율은 vehicle_gps_position_s.hdop, vehicle_gps_position_s.vdop 필드에 미터 단위나 단위 없는 스칼라 계수로 매칭되어 PX4 내부로 유입된다.

GDOP = \sqrt{HDOP^2 + VDOP^2 + TDOP^2}
(TDOP: Time DOP, 시계 오차)

PX4의 EKF 시스템은 내부적으로 이 계수들을 위치 측정치 잡음 공분산 매트릭스(Position Observation Noise Covariance Matrix, R)에 곱하여 증폭시킨다. 즉, HDOP가 2.0에서 4.0으로 두 배 뛰면, EKF는 일시적으로 GPS의 위치 데이터를 두 배 더 불신(Distrust)하게 되고, 상대적으로 IMU(관성 센서)의 적분 가중치를 끌어올린다.

2. PlotJuggler를 이용한 기하학 배열 품질 저하 구역 식별법

특정 비행 경로에서 기체가 궤적 추종(Trajectory Following)에 지속해서 실패하는 ‘블랙홀(Blackhole)’ 구간이 존재한다면, ULog 데이터를 통해 GDOP 기반 공간 맵핑을 수행해야 한다.

  1. 동기화 플로팅(Synchronized Plotting):
    • PlotJuggler 화면의 상단에 vehicle_local_position_s.x, y (EKF 융합 위치)와 vehicle_gps_position_s.lat, lon (GPS 원시 위치)의 차이(Error)를 그리고,
    • 하단 화면에 vehicle_gps_position_s.hdop 및 satellites_used를 겹쳐 그린다.
  2. 은폐된 결함(Hidden Flaw) 발견:
    분석 도중 위성 수는 14개로 넉넉하게 유지됨에도 불구하고, 특정 타임스탬프에서 hdop 값만 1.5에서 3.5로 수직 스파이크(Spike)를 치는 구간이 발견될 수 있다.
  3. 원인 규명:
    이는 드론이 높은 건물 측면에 바짝 붙어 비행함에 따라, 한쪽 반구(예: 동쪽 방향)의 위성 신호가 모조리 차단되고, 반대쪽 반구(서쪽)의 위성들만 집중적으로 수신되는 ‘반쪽 하늘(Half-sky)’ 현상 때문이다. 기하학적으로 위성들이 일렬로 줄을 선 셈이 되므로, GPS 기반 삼각측량의 교차점이 길쭉한 타원형(Error Ellipse)으로 길게 찢어져 HDOP가 폭발한 것이다.

3. 통계적 한계를 넘어서: 커스텀 알람 레이어 설계

산업용 드론 운용 시 이러한 현상을 사후에 ULog로만 분석하는 것은 소 잃고 외양간 고치기 격이다.
지상 관제 시스템(QGroundControl 등) 개발 시, 단순 위성 수(N > 10)를 기준으로 이륙 조건(Pre-flight Check)을 승인하는 관행은 폐기되어야 한다. 대신 MAVLink 메시지를 서브스크라이빙하여 아래와 같은 로직을 GCS 또는 컴패니언 컴퓨터(Companion Computer)에 구현해야 한다.

graph TD;
    A[Telemetry 수신: GPS_RAW_INT] --> B{eph, epv 및 hdop 추출}
    B --> C{hdop > 2.5 5초 지속?}
    C -- 예(Yes) --> D[지상 관제국 팝업 경고: "지형 지물 차폐 강함"]
    C -- 아니오(No) --> E{위성 수 < 10?}
    E -- 예(Yes) --> F[시스템 결함: GPS 위상 소실 주의]
    D --> G[EKF VIO 융합 모듈 강제 우선순위 격상 유도]
    G --> H[안전 비행 유지]

원시 데이터의 기하학적 맹점(Blind Spot)을 이해하면, 오토파일럿이 센서 융합의 딜레마에 빠져 허우적대기 전에 미연에 사고를 방지할 수 있는 진정한 의미의 예측적 유지보수(Predictive Maintenance) 시스템을 구축할 수 있다.