| 양쪽 이전 판이전 판다음 판 | 이전 판 |
| gnss [2026/04/13 12:02] – GNSS sync flyingtext | gnss [2026/04/13 12:02] (현재) – GNSS sync flyingtext |
|---|
| === 위성 궤도 배치와 배치 전략 === | === 위성 궤도 배치와 배치 전략 === |
| |
| 전 지구적 피복 범위를 확보하기 위한 중궤도 위성의 배치 원리와 궤도 평면 설계를 다룬다. | [[위성 항법 시스템]](GNSS)의 전 지구적 피복 범위를 달성하기 위한 핵심 전략은 [[중궤도]](Medium Earth Orbit, MEO)를 활용한 최적의 [[위성군]](Constellation) 배치에 있다. 중궤도는 약 19,000km에서 24,000km 사이의 고도에 위치하며, 이는 [[저궤도]](Low Earth Orbit, LEO)에 비해 한 기의 위성이 커버할 수 있는 지상 면적이 넓으면서도 [[정지 궤도]](Geostationary Orbit, GEO)에 비해 신호 감쇠와 지연 시간이 적다는 공학적 이점을 지닌다. 전 지구 어디에서나 최소 4기 이상의 [[가시 위성]]을 확보하기 위해, GNSS 설계자들은 위성의 수, [[궤도 평면]]의 개수, 그리고 각 평면 사이의 상대적 위치 관계를 수학적으로 정의하는 [[Walker Delta]] 배치 방식을 주로 채택한다. |
| | |
| | 위성군의 기하학적 구조를 설계할 때 가장 우선적으로 고려되는 요소는 [[궤도 경사각]](Inclination)과 [[승교점 적경]](Right Ascension of the Ascending Node, RAAN)의 분배이다. 전 지구적 서비스를 목표로 하는 시스템은 위성이 적도 상공에만 머물지 않고 고위도 지역까지 충분히 통과할 수 있도록 약 55도에서 65도 사이의 궤도 경사각을 설정한다. 예를 들어, 미국의 [[지피에스]](GPS)는 55도의 경사각을 가진 6개의 궤도 평면에 위성을 배치하여 전 지구적 균일성을 확보하는 반면, 러시아의 [[글로나스]](GLONASS)는 고위도 지역의 가시성을 높이기 위해 64.8도의 상대적으로 높은 경사각을 가진 3개의 궤도 평면을 운용한다((Decadal evolution of GPS, GLONASS, and Galileo mean orbital elements, https://link.springer.com/article/10.1007/s10291-024-01708-5 |
| | )). 이러한 궤도 평면의 다각화는 특정 지역에서 위성이 일렬로 배치되어 발생하는 [[정밀도 저하율]](Dilution of Precision, DOP)의 급격한 상승을 방지하고, 기하학적 강건성을 제공한다. |
| | |
| | 배치 전략의 또 다른 핵심은 동일 궤도 평면 내 위성 간의 위상차와 서로 다른 평면 간 위성들의 상대적 위치를 결정하는 ’위상 고정(Phasing)’이다. 이는 $i: T/P/F$라는 기호로 표현되는 Walker 배치 이론에 근거하는데, 여기서 $T$는 총 위성 수, $P$는 궤도 평면 수, $F$는 인접 평면 간의 위상 인자를 의미한다. 유럽 연합의 [[갈릴레오]](Galileo)는 3개의 궤도 평면에 각 8기씩 총 24기의 위성을 배치하는 $56^\circ: 24/3/1$ 구조를 기본 골격으로 삼아, 위성 한 기가 궤도를 이탈하더라도 서비스 공백을 최소화할 수 있는 중첩 피복 전략을 구사한다((THE INTEROPERABLE GLOBAL NAVIGATION SATELLITE SYSTEMS SPACE SERVICE VOLUME, https://www.unoosa.org/res/oosadoc/data/documents/2021/stspace/stspace75rev_1_0_html/st_space_75rev01E.pdf |
| | )). 반면 중국의 [[베이두]](BeiDou)는 중궤도 위성뿐만 아니라 [[경사 동기 궤도]](Inclined Geosynchronous Orbit, IGSO)와 정지 궤도 위성을 혼합한 하이브리드 배치를 통해, 아시아-태평양 지역의 정밀도를 우선적으로 강화하면서 전 지구적 서비스를 병행하는 독자적인 전략을 취한다. |
| | |
| | 위성 궤도 배치는 단순히 초기 발사 단계에서 끝나는 것이 아니라, 지구 중력의 불균일성이나 태양 복사압 등에 의한 [[섭동]]으로 발생하는 궤도 이탈을 지속적으로 보정하는 유지 관리 전략을 포함한다. 각 위성은 사전에 설계된 궤도 슬롯(Slot) 내에 머물러야 하며, 이를 위해 지상 제어 부문은 주기적으로 위성의 위치를 측정하고 필요시 탑재된 추진기를 이용해 궤도를 수정한다. 이러한 정밀한 궤도 배치와 유지 전략은 사용자가 지상에서 수신하는 신호의 기하학적 배치 상태를 최적으로 유지하게 함으로써, 삼차원 위치 결정의 오차 범위를 최소화하고 시스템의 신뢰성을 보장하는 물리적 토대가 된다. |
| |
| === 위성 탑재 장비와 신호 송신기 === | === 위성 탑재 장비와 신호 송신기 === |
| === 지상 관제소의 역할 === | === 지상 관제소의 역할 === |
| |
| 위성의 궤도를 추적하고 원격 측정 데이터를 수집하여 시스템 전체를 운영하는 과정을 설명한다. | 제어 부문(Control Segment)은 우주 공간에 배치된 위성군이 설계된 기술적 사양을 준수하며 정확한 항법 정보를 송출할 수 있도록 관리하는 시스템의 중추적인 신경망이다. 지상 관제소의 일차적인 역할은 위성으로부터 수신한 [[원격 측정]](Telemetry) 데이터를 분석하여 각 위성의 기계적 상태와 궤도 환경을 실시간으로 감시하는 것이다. 위성 항법 시스템의 정확도는 위성이 자신의 위치와 시각 정보를 얼마나 정밀하게 알고 있느냐에 달려 있으므로, 제어 부문은 위성의 궤도와 시계 오차를 끊임없이 보정하는 일련의 루프(Loop)를 형성한다. 이를 위해 제어 부문은 전 세계에 전략적으로 분산 배치된 [[모니터 스테이션]](Monitor Station), 수집된 데이터를 통합 처리하는 [[마스터 관제소]](Master Control Station, MCS), 그리고 위성에 수정된 데이터를 전송하는 [[지상 안테나]](Ground Antenna)로 구성된다. |
| | |
| | 전 세계 각지에 위치한 모니터 스테이션은 가시권 내에 있는 모든 위성으로부터 항법 신호를 수신하여 [[의사거리]](Pseudorange)와 [[반송파 위상]](Carrier Phase) 관측값을 수집한다. 이 스테이션들은 위치 좌표가 정밀하게 기지된 지점에 설치되어 있어, 위성 신호에서 발생하는 미세한 오차를 역으로 산출할 수 있는 기준점 역할을 수행한다. 수집된 원시 데이터는 통신망을 통해 마스터 관제소로 전송되며, 이곳에서 시스템의 핵심적인 연산 과정이 이루어진다. 마스터 관제소는 수집된 관측 데이터를 바탕으로 [[칼만 필터]](Kalman Filter)와 같은 고도의 통계적 추정 알고리즘을 적용하여 위성의 실제 위치와 예정된 위치 사이의 편차를 계산하고, 향후 위성이 이동할 경로인 [[궤도력]](Ephemeris)을 예측한다. 이와 동시에 위성에 탑재된 [[원자시계]]의 드리프트(Drift) 현상을 분석하여 시스템 표준시와의 동기화를 위한 시계 보정 계수를 산출한다. |
| | |
| | 마스터 관제소에서 생성된 최신 궤도 정보와 시계 보정 데이터는 [[항법 메시지]](Navigation Message)의 형태로 가공된다. 이렇게 갱신된 데이터는 지상 안테나를 통해 해당 위성으로 무선 송신되는데, 이를 [[업로드]](Upload) 과정이라 한다. 위성은 지상으로부터 받은 최신 메시지를 자신의 메모리에 저장한 후, 이를 다시 지상의 사용자들에게 방송한다. 이러한 갱신 과정은 통상적으로 하루에 수차례 반복되며, 이를 통해 위성이 궤도 [[섭동]](Perturbation)이나 시계 오차로 인해 발생할 수 있는 누적 오차를 최소화하도록 돕는다. 만약 지상 관제소가 위성의 신호에 이상이 있음을 감지하면, 해당 위성의 상태 지수(Health bit)를 ’비정상’으로 변경하여 사용자들이 잘못된 정보를 이용하지 않도록 조치함으로써 시스템의 [[무결성]](Integrity)을 유지한다. |
| | |
| | 또한 지상 관제소는 위성의 물리적인 생존 가능성을 관리하는 임무도 병행한다. 위성의 태양 전지판이 태양을 올바르게 지향하고 있는지, 배터리의 충방전 상태가 양호한지, 그리고 자세 제어를 위한 연료가 충분한지를 상시 점검한다. 위성이 중력적 요인으로 인해 지정된 궤도 구역을 벗어날 조짐이 보이면, 관제소는 위성의 추진기를 가동하는 [[궤도 기동]] 명령을 하달하여 위성을 정해진 위치로 복귀시킨다. 이처럼 지상 관제소는 단순한 데이터 중계를 넘어, 우주 자산의 수명을 관리하고 전 지구적 서비스의 [[연속성]]을 보장하는 운영 주체로서 기능한다. 제어 부문의 이러한 정밀한 관리 체계가 뒷받침되지 않는다면, 우주 부문의 위성들은 시간이 지남에 따라 점차 정확도를 상실하여 항법 시스템으로서의 가치를 잃게 된다.((Global Positioning System Standard Positioning Service Performance Standard, https://www.gps.gov/sites/default/files/2020-05/2020-SPS-performance-standard.pdf |
| | )) |
| |
| === 궤도 결정과 시각 동기화 === | === 궤도 결정과 시각 동기화 === |
| === 수신기 구조와 안테나 === | === 수신기 구조와 안테나 === |
| |
| 미약한 위성 신호를 포착하여 디지털 신호로 변환하는 하드웨어적 구성 요소를 다룬다. | 사용자 부문에서 핵심적인 하드웨어인 수신기는 위성으로부터 송출된 극미약한 무선 주파수(Radio Frequency, RF) 신호를 포착하여 사용자가 해독 가능한 디지털 정보로 변환하는 역할을 수행한다. 위성 항법 신호는 약 20,000km 이상의 거리를 전파하며 확산되기 때문에, 지표면에 도달할 무렵의 신호 강도는 대략 $-160 \, \text{dBW}$ 정도로 매우 낮다. 이는 주변 환경의 [[열잡음]](Thermal noise)보다도 낮은 수준이므로, 수신기는 고감도의 안테나와 정밀한 RF 프런트엔드(RF Front-end) 설계를 통해 신호를 복원해야 한다. |
| | |
| | 안테나는 위성 신호를 수집하는 첫 번째 단계로, GNSS 신호의 특성인 [[우측 원편파]](Right-Hand Circularly Polarized, RHCP)를 효과적으로 수신할 수 있도록 설계된다. 대다수의 GNSS 안테나는 소형화와 제작 용이성을 위해 [[패치 안테나]](Patch antenna) 형식을 취하며, 지면에서 반사되어 들어오는 [[다중 경로]](Multipath) 신호를 차단하기 위해 특정 각도 이하의 신호를 억제하는 이득 패턴을 갖는다. 안테나의 물리적 중심과 전자기적 중심인 [[위상 중심]](Phase center)의 일치 여부는 cm 단위의 정밀 측위에서 매우 중요한 요소이다. 안테나 바로 뒷단에는 [[저잡음 증폭기]](Low Noise Amplifier, LNA)가 배치되어, 신호에 포함된 잡음의 지수를 최소화하면서도 미약한 신호를 후속 처리가 가능한 수준으로 증폭한다. |
| | |
| | RF 프런트엔드는 안테나를 통해 들어온 고주파 신호를 처리하기 쉬운 낮은 주파수 대역으로 변환하는 기능을 담당한다. 이 과정은 대개 [[슈퍼헤테로다인]](Superheterodyne) 방식이나 직접 변환(Direct conversion) 방식을 통해 이루어진다. [[국부 발진기]](Local Oscillator, LO)에서 생성된 신호를 [[혼합기]](Mixer)에서 위성 신호와 결합하여 [[중간 주파수]](Intermediate Frequency, IF)를 생성하며, 이 과정에서 [[대역 통과 필터]](Band Pass Filter)를 사용하여 인접 대역의 간섭 신호를 제거한다. 이후 [[아날로그-디지털 변환기]](Analog-to-Digital Converter, ADC)를 거쳐 아날로그 IF 신호는 디지털 샘플 데이터로 양자화된다. 양자화 비트 수가 높을수록 신호의 해상도가 정밀해지지만, 연산량과 전력 소모가 증가하므로 용도에 따라 적절한 비트 수가 선택된다((Advanced System Analysis and Survey on the GPS Receiver Front End, https://ieeexplore.ieee.org/ielx7/6287639/9668973/09726792.pdf |
| | )). |
| | |
| | 디지털로 변환된 신호는 기저대역(Baseband) 처리부로 전달된다. 기저대역부의 핵심 하드웨어 요소는 [[상관기]](Correlator)이며, 이는 수신기 내부에서 생성된 복제 코드와 수신된 디지털 신호를 비교하여 신호의 도달 시간을 측정한다. 현대의 GNSS 수신기는 수십 개에서 수백 개의 채널을 동시에 운용하며, 각 채널은 특정 위성의 신호를 추적하기 위해 독립적인 [[수치 제어 발진기]](Numerically Controlled Oscillator, NCO)와 루프 필터를 포함한다. 이러한 하드웨어 구조는 [[지연 고정 루프]](Delay Locked Loop, DLL)와 [[위상 고정 루프]](Phase Locked Loop, PLL)를 통해 위성 신호의 코드와 반송파 위상을 정밀하게 유지한다. |
| | |
| | 수신기 내부의 모든 처리 과정은 기준 시계의 통제를 받는다. 위성 항법은 시간 측정의 정밀도가 곧 위치 정확도로 직결되므로, 수신기 내부에는 일정한 주파수를 유지할 수 있는 [[온도 보상 수정 발진기]](Temperature Compensated Crystal Oscillator, TCXO) 또는 [[오븐 제어 수정 발진기]](Oven Controlled Crystal Oscillator, OCXO)가 탑재된다. 위성의 원자시계만큼 정밀하지는 않으나, 수신기는 위치 산출 과정에서 네 개 이상의 위성 신호를 이용하여 자신의 시계 오차를 미지수로 두고 계산함으로써 하드웨어적인 시계 오차를 수학적으로 보정한다((GLOBAL POSITIONING SYSTEM STANDARD POSITIONING SERVICE SIGNAL SPECIFICATION, https://www.gps.gov/technical/ps/1995-SPS-signal-specification.pdf |
| | )). |
| |
| === 신호 처리와 위치 산출 알고리즘 === | === 신호 처리와 위치 산출 알고리즘 === |
| |
| 코드 및 위상 관측값을 이용하여 수신자의 좌표를 계산하는 수학적 과정을 설명한다. | 수신기에 도달한 [[위성 항법 시스템]](GNSS) 신호로부터 위치를 산출하는 과정은 크게 신호 처리(Signal Processing) 단계와 항법 해 도출(Navigation Solution) 단계로 구분된다. 수신기는 안테나를 통해 수신된 초고주파(RF) 신호를 [[기저 대역]](Baseband) 신호로 변환한 후, 신호 획득(Acquisition)과 추적(Tracking) 과정을 거쳐 항법 메시지와 관측값을 추출한다. 신호 획득 단계에서는 수신된 신호와 수신기 내부에서 생성한 [[복제 코드]](Replica code) 사이의 상관관계를 분석하여 [[코드 지연]](Code delay)과 [[도플러 편이]](Doppler shift)를 대략적으로 찾아낸다. 이후 신호 추적 단계에서는 [[지연 잠금 루프]](Delay Lock Loop, DLL)와 [[위상 잠금 루프]](Phase Lock Loop, PLL)를 사용하여 코드와 위상의 변화를 정밀하게 추적하며, 이를 통해 [[의사거리]](Pseudorange)와 [[반송파 위상]](Carrier phase) 관측값을 생성한다. |
| | |
| | 수신기가 위치를 계산하기 위해 사용하는 가장 기본적인 수학적 모델은 의사거리 관측 방정식이다. 수신기 시간 $ t_u $와 위성 시간 $ t^s $ 사이의 동기화 오차를 고려할 때, 특정 위성 $ i $에 대한 의사거리 $ ^{(i)} $는 다음과 같은 비선형 방정식으로 표현된다. |
| | |
| | $$ \rho^{(i)} = \sqrt{(x^{(i)} - x_u)^2 + (y^{(i)} - y_u)^2 + (z^{(i)} - z_u)^2} + c(dt_u - dt^{(i)}) + I^{(i)} + T^{(i)} + \epsilon^{(i)} $$ |
| | |
| | 여기서 $ (x^{(i)}, y^{(i)}, z^{(i)}) $는 위성의 좌표, $ (x_u, y_u, z_u) $는 수신기의 미지 좌표이다. $ c $는 진공에서의 광속이며, $ dt_u $와 $ dt^{(i)} $는 각각 수신기와 위성의 시계 오차를 나타낸다. $ I^{(i)} $와 $ T^{(i)} $는 전리층 및 대류권 지연 오차이며, $ ^{(i)} $는 수신기 잡음 및 다중 경로 오차를 포함한 잔여 오차항이다. 관측 방정식은 수신기의 좌표와 시계 오차라는 4개의 미지수를 포함하고 있으므로, 최소 4기 이상의 위성으로부터 신호를 수신해야 유일한 해를 구할 수 있다. |
| | |
| | 실제 위치 산출 알고리즘에서는 비선형인 관측 방정식을 계산하기 용이하도록 [[테일러 급수]](Taylor series)를 이용하여 선형화한다. 수신기의 초기 추정 위치 $ (x_0, y_0, z_0) $를 기준으로 방정식을 전개하면, 관측값의 변화량과 미지수의 보정량 사이의 관계를 나타내는 행렬식 $ = + $를 얻게 된다. 여기서 $ $는 각 위성에 대한 방향 코사인 성분으로 구성된 [[야코비 행렬]](Jacobian matrix)이다. 정지 상태의 수신기는 주로 [[최소제곱법]](Least Squares Method)을 사용하여 잔차의 제곱합을 최소화하는 방향으로 위치를 결정하며, 이동 중인 사용자의 경우에는 이전 시점의 상태 정보를 바탕으로 현재 위치를 예측하고 보정하는 [[칼만 필터]](Kalman Filter) 알고리즘이 널리 사용된다((Sanz Subirana, J., Juan Zornoza, J. M., & Hernández-Pajares, M. (2013). GNSS Data Processing, Volume I: Fundamentals and Algorithms. European Space Agency. https://gssc.esa.int/navipedia/GNSS_Book/ESA_GNSS-Book_TM-23_Vol_I.pdf |
| | )). |
| | |
| | 고정밀 위치 결정이 필요한 경우에는 코드 관측값보다 오차가 훨씬 작은 반송파 위상 관측값을 활용한다. 그러나 반송파 위상은 신호의 파장 단위로 측정되기 때문에, 수신기가 신호를 처음 포착했을 때 몇 주기의 파장이 지나갔는지 알 수 없는 [[정수 모호성]](Integer Ambiguity) 문제가 발생한다. 이를 해결하기 위해 [[이중 차분]](Double Difference) 기법을 사용하여 수신기 및 위성의 시계 오차를 제거하고, [[LAMBDA 알고리즘]](Least-squares AMBiguity Decorrelation Adjustment, LAMBDA)과 같은 통계적 기법을 동원하여 정수 모호성을 결정한다. 정수 모호성이 정확히 해소(Fix)되면 수 밀리미터에서 수 센티미터 수준의 정밀한 위치 산출이 가능해진다. |
| |
| ===== 작동 원리와 신호 체계 ===== | ===== 작동 원리와 신호 체계 ===== |
| === 대류권 지연 효과 === | === 대류권 지연 효과 === |
| |
| 수증기와 기압 등 대기 하층부의 기상 상태로 인해 발생하는 신호 지연 현상을 기술한다. | [[대류권]](Troposphere)은 지표면으로부터 약 10~16km 고도까지 형성된 대기층으로, 전체 대기 질량의 약 80% 이상이 집중되어 있으며 기상 현상이 발생하는 주된 공간이다. [[위성 항법 시스템]](GNSS)의 신호가 이 영역을 통과할 때, 대기를 구성하는 중성 기체 분자와 수증기에 의해 신호의 속도가 지연되고 경로가 굴절되는 현상이 발생한다. 이를 대류권 지연(Tropospheric delay)이라 하며, 이는 전리층 지연과 달리 신호의 주파수에 의존하지 않는 [[비분산 매질]](Non-dispersive medium)의 특성을 가진다. 따라서 다중 주파수 수신기를 사용하더라도 주파수 조합을 통해 오차를 직접 제거할 수 없으며, 반드시 대기압, 온도, 습도 등의 기상 매개변수를 기반으로 한 수학적 모델을 통해 보정해야 한다. |
| | |
| | 대류권 지연의 크기는 전파가 통과하는 경로상의 [[굴절률]](Refractive index) $ n $에 의해 결정된다. 실제 계산에서는 굴절률과 1의 차이가 매우 미세하므로, 이를 수치적으로 다루기 쉬운 [[굴절도]](Refractivity) $ N $으로 변환하여 사용한다. 굴절도는 대기의 물리적 상태에 따라 다음과 같이 정의된다. |
| | |
| | $$ N = (n - 1) \times 10^6 = N_{\text{dry}} + N_{\text{wet}} $$ |
| | |
| | 여기서 $ N_{} $는 질소와 산소 등 건조 기체에 의한 건조 분력을, $ N_{} $은 수증기에 의한 습윤 분력을 의미한다. 이에 따라 전체 대류권 지연량은 성격이 다른 두 가지 구성 요소, 즉 [[정역학적 지연]](Hydrostatic delay)과 [[습윤 지연]](Wet delay)의 합으로 표현된다. |
| | |
| | 정역학적 지연은 대류권 전체 지연량의 약 90%를 차지하며, 지표면의 [[기압]]과 밀접한 관련이 있다. 건조 대기는 [[정역학적 평형]](Hydrostatic equilibrium) 상태에 있다고 가정할 수 있으므로, 지표 기압을 정밀하게 측정하면 모델을 통해 수 밀리미터 이내의 높은 정확도로 예측이 가능하다. 위성이 관측자의 머리 위에 위치할 때의 [[천정 정역학적 지연]](Zenith Hydrostatic Delay, ZHD)은 해수면 기준으로 약 2.3m에 달한다. |
| | |
| | 반면 습윤 지연은 대기 중의 [[수증기]] 분자에 의해 발생한다. 수증기는 건조 기체와 달리 시공간적 변동성이 매우 크고 수직 분포가 불균일하여 물리적 모델만으로 정확히 예측하기 어렵다. [[천정 습윤 지연]](Zenith Wet Delay, ZWD)은 통상 0~30cm 범위에 머물러 절대적인 크기는 정역학적 지연보다 작지만, 불확실성이 높아 GNSS 정밀 위치 결정에서 주요한 오차 원인이 된다. |
| | |
| | 위성의 고도각($ e $)이 낮아질수록 신호가 통과하는 대기 경로의 길이는 길어지며 지연량 또한 증가한다. 이를 보정하기 위해 천정 방향의 지연량에 [[사영 함수]](Mapping function)를 곱하여 실제 신호 경로상의 사선 지연(Slant delay)을 산출한다. |
| | |
| | $$ \Delta_{\text{trop}}(e) = ZHD \cdot m_h(e) + ZWD \cdot m_w(e) $$ |
| | |
| | 여기서 $ m_h(e) $와 $ m_w(e) $는 각각 정역학적 성분과 습윤 성분에 대한 사영 함수이다. 초기의 사영 함수는 단순한 삼각함수 형태였으나, 현대의 정밀 분석에서는 대기의 곡률과 온도 구배를 고려한 [[사스타모이넨 모델]](Saastamoinen model), [[홉필드 모델]](Hopfield model), 그리고 기상 관측 데이터를 기반으로 계수를 산출하는 [[닐 사영 함수]](Niell Mapping Function, NMF) 등이 널리 사용된다((Niell, A.E, “Global Mapping for the Atmospheric Delay at Radio Wavelenghts”, Journal of Geophysical Research 111 (B2), 3227-3246. 1996, https://www.sciepub.com/reference/122309 |
| | )). |
| | |
| | 대류권 지연 효과는 단순한 오차 요인에 그치지 않고, 역으로 대기 상태를 파악하는 중요한 정보원으로 활용된다. GNSS 관측 데이터에서 정역학적 지연을 제거하고 남은 습윤 지연량을 분석하면 대기 중의 전체 수증기량을 산출할 수 있는데, 이를 [[가청 수증기량]](Precipitable Water Vapor, PWV)이라 한다. 이러한 기법은 [[GNSS 기상학]]이라는 학문 분야로 발전하여 수치 예보 모델의 정확도를 높이고 국지성 호우를 예측하는 등 기상 관측 인프라의 핵심적인 역할을 수행하고 있다((Tropospheric Delay - Navipedia, https://gssc.esa.int/navipedia/index.php/Tropospheric_Delay |
| | )). |
| |
| ==== 기하학적 오차와 환경적 요인 ==== | ==== 기하학적 오차와 환경적 요인 ==== |
| === 위성 배치에 따른 정밀도 저하율 === | === 위성 배치에 따른 정밀도 저하율 === |
| |
| 수신자가 관측하는 위성들의 기하학적 분포가 위치 결정 정밀도에 미치는 상관관계를 설명한다. | [[위성 항법 시스템]](GNSS)을 이용한 위치 결정 과정에서 발생하는 오차는 신호의 전파 지연이나 수신기 잡음과 같은 물리적 요인뿐만 아니라, 관측 당시 위성들의 공간적 배치 상태에 의해서도 결정적인 영향을 받는다. 이를 정량화한 지표가 [[정밀도 저하율]](Dilution of Precision, DOP)이다. DOP는 위성과 수신기 사이의 기하학적 배치가 사용자 위치 해의 정밀도를 얼마나 저하시키는지를 나타내는 무차원 계수이다. 동일한 수준의 [[의사거리]](Pseudorange) 측정 오차가 존재하더라도, 가시 위성들이 하늘의 좁은 영역에 모여 있는 경우보다 넓게 고루 분산되어 있을 때 최종적인 위치 오차는 작아진다. 이는 [[삼변측량]]의 원리상 각 위성으로부터 측정된 거리 정보가 교차하는 기하학적 형태가 견고할수록 외부 잡음에 대한 저항력이 높아지기 때문이다. |
| | |
| | 수학적으로 DOP는 선형화된 관측 방정식의 [[기하학적 설계 행렬]](Geometric Design Matrix)로부터 유도된다. 수신기의 위치와 시계 편차를 포함한 상태 벡터의 변화량을 $ $, 관측된 의사거리 잔차 벡터를 $ $, 그리고 관측 잡음을 $ $이라 할 때, 선형화된 관계식은 $ = + $으로 표현된다. 여기서 행렬 $ $는 수신기로부터 각 위성을 향하는 단위 방향 벡터와 시계 편차 항으로 구성된다. [[최소자승법]](Least Squares Method)을 통해 추정된 상태 벡터의 [[공분산 행렬]](Covariance Matrix)은 측정 오차의 분산 $ ^2 $과 설계 행렬의 조합인 $ = (^T )^{-1} $의 곱으로 나타난다. 이 행렬 $ $의 대각 성분들이 바로 각 성분별 정밀도 저하율의 제곱에 해당한다((A closed-form formula to calculate geometric dilution of precision (GDOP) for multi-GNSS constellations, https://link.springer.com/article/10.1007/s10291-015-0440-x |
| | )). |
| | |
| | DOP는 분석하고자 하는 차원에 따라 여러 세부 지표로 분류된다. [[기하학적 정밀도 저하율]](Geometric DOP, GDOP)은 3차원 위치 좌표와 수신기 시계 편차를 모두 포함한 전체적인 정밀도를 의미하며, $ $ 행렬의 모든 대각 성분의 합에 루트를 씌워 계산한다. [[위치 정밀도 저하율]](Position DOP, PDOP)은 시계 성분을 제외한 3차원 공간상의 위치 정밀도만을 나타내며, 이를 다시 수평면 성분인 [[수평 정밀도 저하율]](Horizontal DOP, HDOP)과 수직 성분인 [[수직 정밀도 저하율]](Vertical DOP, VDOP)로 분리할 수 있다. 시계 오차의 정밀도만을 별도로 나타낼 때는 [[시간 정밀도 저하율]](Time DOP, TDOP)을 사용한다. |
| | |
| | 기하학적 관점에서 볼 때, 위성들이 형성하는 사면체의 부피가 클수록 DOP 값은 작아지며 이는 높은 정밀도를 의미한다. 예를 들어, 한 기의 위성이 사용자의 머리 위인 [[천정]](Zenith) 부근에 위치하고 나머지 위성들이 낮은 [[고도각]]에서 서로 넓은 [[방위각]] 간격을 두고 배치될 때 가장 이상적인 PDOP 값을 얻을 수 있다. 반대로 모든 위성이 특정 방향에 치우쳐 배치되면 기하학적 강인성이 떨어져 작은 측정 오차에도 위치 해가 크게 흔들리게 된다. |
| | |
| | 일반적으로 GNSS 관측에서 VDOP는 HDOP보다 항상 큰 값을 갖는 경향이 있다. 이는 수신기를 기준으로 하늘의 윗부분에만 위성이 존재하고, 지구 본체에 가려진 아래쪽 절반 공간에는 위성을 배치할 수 없다는 물리적 한계 때문이다. 이러한 기하학적 비대칭성으로 인해 수직 방향의 위치 결정은 수평 방향에 비해 본질적으로 오차에 취약한 구조를 갖는다. 따라서 정밀한 고도 정보가 요구되는 [[항공 항법]]이나 자율 주행 분야에서는 단순히 가시 위성의 수를 확보하는 것을 넘어, 낮은 DOP를 유지할 수 있는 위성 선택 알고리즘과 [[다중 GNSS]](Multi-GNSS) 통합 활용이 필수적으로 요구된다. |
| |
| === 다중 경로 오차와 수신기 잡음 === | === 다중 경로 오차와 수신기 잡음 === |
| |
| 건물이나 지면에 반사된 신호가 수신기에 도달하여 발생하는 간섭 현상과 장비 자체의 오차를 다룬다. | [[위성 항법 시스템]](GNSS)의 신호가 위성에서 송출되어 수신기에 도달하는 과정에서 발생하는 환경적 오차 중 가장 지배적인 요인 중 하나는 [[다중 경로]](Multipath) 현상이다. 다중 경로 오차는 위성으로부터 발신된 직접 신호(Direct signal) 외에 지면, 건물, 수면 등 수신기 주변의 지형지물에 반사된 하나 이상의 반사 신호(Reflected signal)가 수신 안테나에 함께 유입되면서 발생한다. 반사 신호는 직접 신호에 비해 항상 더 긴 경로를 이동하므로 수신기에 지연되어 도착하며, 직접 신호와 중첩되어 [[상관 함수]](Correlation function)의 기하학적 형상을 왜곡시킨다. 이러한 왜곡은 수신기가 신호의 도달 시간을 추적하는 [[지연 고정 루프]](Delay Lock Loop, DLL)에서 정합 지점의 편차를 유발하며, 결과적으로 [[의사거리]](Pseudorange) 측정치에 오차를 삽입하게 된다. |
| | |
| | 다중 경로 오차의 크기는 신호의 종류에 따라 상이하게 나타난다. 코드 관측값을 이용하는 경우에는 오차의 크기가 수 미터에서 환경에 따라 수십 미터에 달할 수 있으나, [[반송파 위상]](Carrier phase) 관측값의 경우 오차의 최대치가 반송파 파장의 4분의 1 수준인 수 센티미터 이내로 제한되는 특성을 갖는다. 또한 다중 경로 신호는 직접 신호와 비교하여 위상과 진폭이 다르기 때문에 두 신호가 보강 간섭 또는 상쇄 간섭을 일으키며 [[신호 대 잡음비]](Signal-to-Noise Ratio, SNR)의 주기적인 변동을 초래한다. 이러한 특성은 동적 환경에서 다중 경로를 검출하고 식별하는 주요 지표로 활용되기도 한다. |
| | |
| | 다중 경로를 억제하기 위한 기술적 대응은 하드웨어와 소프트웨어 양면에서 이루어진다. 하드웨어 측면에서는 [[우원편파]](Right-Hand Circular Polarization, RHCP) 특성을 갖는 GNSS 신호가 반사될 때 좌원편파(LHCP)로 변하는 성질을 이용하여 반사파를 차단하거나, 안테나 주변에 물리적인 차폐 구조물인 [[초크 링]](Choke ring)을 설치하여 지면으로부터의 반사파 유입을 최소화한다. 소프트웨어 측면에서는 수신기 내부의 [[상관기]](Correlator) 간격을 좁게 설정하는 [[협대역 상관기]](Narrow Correlator) 기술이나, 상관 함수의 미분 특성을 이용한 고도화된 신호 처리 알고리즘을 통해 반사 신호에 의한 상관 피크의 왜곡을 보정한다. |
| | |
| | 수신기 잡음(Receiver noise)은 다중 경로와 달리 수신기 내부의 전자 부품과 신호 처리 과정에서 발생하는 무작위 오차를 의미한다. 이는 주로 안테나와 수신기 전단부(Front-end)에서 발생하는 [[열잡음]](Thermal noise)에 기인하며, 신호의 [[양자화]](Quantization) 과정에서 발생하는 오차와 내부 [[발진기]](Oscillator)의 불안정성 등도 포함된다. 수신기 잡음은 백색 잡음(White noise)의 특성을 띠므로 장시간 관측을 통한 평균화로 어느 정도 상쇄가 가능하지만, 순간적인 정밀도를 결정하는 근본적인 한계치로 작용한다. |
| | |
| | 수신기 잡음에 의한 거리 측정 오차의 표준편차 $\sigma_{n}$은 대략적으로 다음과 같은 관계를 갖는다. |
| | |
| | $$ \sigma_{n} \propto \frac{1}{\sqrt{C/N_0}} \cdot \frac{B_L}{T_{int}} $$ |
| | |
| | 여기서 $C/N_0$는 반송파 대 잡음 밀도비이며, $B_L$은 추적 루프의 대역폭, $T_{int}$는 적분 시간을 의미한다. 즉, 수신 신호의 세기가 강할수록, 그리고 추적 루프의 대역폭이 좁을수록 잡음의 영향은 줄어든다. 그러나 대역폭을 지나치게 좁히면 수신기의 동적 특성이 저하되어 신호 추적을 놓칠 위험이 있으므로, 설계 시에는 정밀도와 추적 성능 사이의 [[트레이드오프]](Trade-off)를 고려해야 한다. 현대의 고성능 수신기는 수 밀리미터 수준의 잡음 수준을 유지하고 있으나, 다중 경로와 수신기 잡음의 복합적인 영향은 여전히 [[도심 협곡]](Urban canyon)이나 울창한 삼림 지대와 같은 가혹한 환경에서 GNSS 가용성을 저해하는 주요 기술적 과제로 남아 있다. ((Characterization of Signal Quality Monitoring Techniques for Multipath Detection in GNSS Applications, https://mdpi-res.com/d_attachment/sensors/sensors-17-01579/article_deploy/sensors-17-01579-v4.pdf?version=1500360568 |
| | )) |
| |
| ==== 오차 보정 기술 ==== | ==== 오차 보정 기술 ==== |
| === 차분 위성 항법 시스템 === | === 차분 위성 항법 시스템 === |
| |
| 위치를 알고 있는 기준국에서 생성한 보정 정보를 실시간으로 전송하여 오차를 제거하는 방식을 설명한다. | 차분 위성 항법 시스템(Differential Global Navigation Satellite System, DGNSS)은 위치 오차를 유발하는 공통적인 요인들을 상쇄하기 위해, 정확한 좌표를 이미 알고 있는 [[기준국]](Reference Station)에서 생성된 보정 정보를 활용하는 기술이다. GNSS 신호가 우주 공간에서 지상 수신기에 도달하는 과정에서는 [[전리층]] 및 [[대류권]] 지연, 위성 궤도 정보의 미세한 부정확성, 위성 시계 오차 등 다양한 오차 요인이 개입한다. 이러한 오차들은 특정 지역 내에 인접한 수신기들에게 거의 동일하게 작용하는 [[공통 오차]](Common-mode error)의 특성을 갖는다. DGNSS는 이러한 오차의 공간적 상관관계를 이용하여, 기준국에서 실시간으로 산출한 오차 보정치를 주변의 [[이동국]](Rover)에 전송함으로써 항법 정밀도를 획기적으로 개선한다. |
| | |
| | 기준국에서의 보정치 산출 과정은 기하학적 관계를 기반으로 한다. 기준국은 정밀 측량을 통해 자신의 지심 좌표 $(x_b, y_b, z_b)$를 미리 확보하고 있으며, 특정 시점에 궤도 정보로부터 계산된 위성의 위치 $(x_s, y_s, z_s)$를 바탕으로 위성과 자신 사이의 이론적인 [[기하학적 거리]](Geometric range)를 산출할 수 있다. 이때 기준국이 실제로 위성 신호를 수신하여 측정한 [[의사거리]](Pseudorange)와 이론적 거리 사이에는 차이가 발생하는데, 이를 [[의사거리 보정치]](Pseudorange Correction, PRC)라고 정의한다. 이를 수식으로 표현하면 다음과 같다. |
| | |
| | $$ PRC = \rho_{true} - \rho_{measured} $$ |
| | |
| | 여기서 $ %%//%%{true} $는 기준국과 위성 사이의 실제 기하학적 거리이며, $ %%//%%{measured} $는 수신기가 측정한 의사거리이다. 기준국은 가시 범위 내에 있는 모든 위성에 대해 이러한 보정치를 실시간으로 계산한다. 생성된 보정 정보는 [[RTCM]](Radio Technical Commission for Maritime Services)과 같은 국제 표준 데이터 포맷으로 부호화되어, 초단파(VHF), 초고주파(UHF), 또는 [[위성 데이터 링크]]나 인터넷망을 통해 주변의 이동국으로 즉시 송신된다. |
| | |
| | 보정 정보를 수신한 이동국은 자신이 개별 위성으로부터 측정한 의사거리에 기준국이 보내온 보정치를 가산하여 보정된 의사거리를 산출한다. 이동국이 기준국과 충분히 가까운 거리에 있다면, 두 장치가 겪는 대기 지연과 위성 관련 오차는 거의 동일하므로 보정 과정을 통해 이러한 계통적 오차가 대부분 제거된다. 이러한 방식은 일반적인 단독 측위 시 발생하는 수 미터에서 수십 미터의 오차를 수 미터 이내, 정밀한 경우 서브미터(sub-meter) 수준까지 줄여준다. |
| | |
| | 그러나 DGNSS의 보정 효과는 기준국과 이동국 사이의 거리인 [[기선]](Baseline)의 길이에 반비례하여 감소하는 특성을 보인다. 두 지점 사이의 거리가 멀어질수록 신호가 통과하는 대기층의 상태가 서로 달라지며, 위성을 바라보는 각도의 차이로 인해 오차의 공통성이 상실되는 [[공간적 상관성 저하]](Spatial decorrelation) 현상이 발생하기 때문이다. 일반적으로 수백 킬로미터 이내에서는 유효한 보정 성능을 유지하지만, 거리가 멀어질수록 정확도는 점차 저하된다. |
| | |
| | DGNSS는 활용하는 관측값의 종류에 따라 [[코드 차분]](Code-based differential) 방식과 [[반송파 위상 차분]](Carrier-phase differential) 방식으로 분류된다. 코드 차분 방식은 주로 선박의 항해, 차량 관제, 저정밀 지도 제작 등에 널리 쓰이며, 반송파 위상을 이용하는 방식은 [[실시간 이동 측위]](Real-Time Kinematic, RTK)로 발전하여 토목 측량이나 자율주행 자동차, 농업용 로봇 등 센티미터 단위의 극정밀도가 요구되는 분야의 핵심 기술이 되었다. 최근에는 단일 기준국의 한계를 극복하기 위해 광역에 분산된 여러 기준국의 데이터를 중앙 제어국에서 통합 처리하여 보정 모델을 생성하는 [[광역 차분 위성 항법 시스템]](Wide Area DGNSS, WADGNSS) 및 [[위성 기반 보정 시스템]](Satellite Based Augmentation System, SBAS)으로 기술적 확장이 이루어지고 있다. |
| |
| === 정밀 지점 위치 결정 기법 === | === 정밀 지점 위치 결정 기법 === |
| |
| 정밀한 위성 궤도와 시계 정보를 이용하여 단일 수신기로도 높은 정확도를 얻는 고급 해석 기법을 다룬다. | 정밀 지점 위치 결정(Precise Point Positioning, PPP)은 단일 [[GNSS]] 수신기를 이용하여 센티미터(cm)에서 데시미터(dm) 수준의 고정밀 좌표를 산출하는 고급 해석 기법이다. 기존의 [[차분 위성 항법 시스템]](Differential GNSS, DGNSS)이나 [[실시간 이동 측위]](Real-Time Kinematic, RTK) 기법이 인근의 기준국(Reference Station)으로부터 보정 정보를 전달받아 수신기와 기준국 간의 공통 오차를 상쇄하는 방식인 것과 달리, 정밀 지점 위치 결정은 위성 자체의 상태에 대한 정밀한 정보를 직접 활용한다. 이 기법은 전 세계 어디서든 기준국의 가용 여부에 구애받지 않고 독립적으로 고정밀 위치 결정을 수행할 수 있다는 점에서 독보적인 효용성을 가진다. |
| | |
| | 정밀 지점 위치 결정의 핵심은 [[국제 GNSS 서비스]](International GNSS Service, IGS)와 같은 전 지구적 관측망 운영 기구에서 제공하는 정밀 위성 궤도(Precise Orbit)와 정밀 위성 시계(Precise Clock) 정보를 활용하는 데 있다. 일반적인 [[항법 메시지]]에 포함된 궤도 정보는 수 미터 수준의 오차를 포함하고 있으나, 사후 처리를 통해 산출된 정밀 궤도 정보는 수 센티미터 수준의 정확도를 보장한다. 수신기는 이러한 정밀 데이터를 바탕으로 [[의사거리]](Pseudorange)와 [[반송파 위상]](Carrier Phase) 관측값을 결합하여 위치를 계산한다. 특히 반송파 위상 측정치는 잡음이 매우 적어 정밀도 향상에 필수적이지만, [[미지 정수]](Integer Ambiguity)라는 불확정 요소를 포함하고 있어 이를 정확히 추정하는 과정이 수반되어야 한다. |
| | |
| | 수학적으로 정밀 지점 위치 결정은 [[이중 주파수]](Dual-frequency) 관측값을 사용하여 [[전리층 지연]](Ionospheric delay)의 1차 항을 제거하는 전리층 제거 조합(Ionosphere-free combination) 모델을 기본으로 한다. 전리층 제거 조합된 의사거리 $ P_{IF} $와 반송파 위상 $ L_{IF} $의 관측 방정식은 다음과 같이 표현된다. |
| | |
| | $$ P_{IF} = \rho + c(dt - dT) + T_{trop} + \epsilon_P $$ $$ L_{IF} = \rho + c(dt - dT) + T_{trop} + \lambda_{IF}N_{IF} + \epsilon_L $$ |
| | |
| | 여기서 $ $는 위성과 수신기 사이의 기하학적 거리이며, $ c $는 광속, $ dt $와 $ dT $는 각각 수신기와 위성의 시계 오차이다. $ T_{trop} $는 [[대류권 지연]](Tropospheric delay)을 의미하며, $ %%//%%{IF} $와 $ N%%//%%{IF} $는 전리층 제거 조합에서의 파장과 미지 정수를 나타낸다. 정밀 지점 위치 결정에서는 위성 시계 오차 $ dT $를 외부에서 제공된 정밀 정보로 고정하므로, 수신기는 자신의 위치 좌표와 시계 오차, 그리고 대류권의 습윤 지연 성분 및 미지 정수를 미지수로 설정하여 최적해를 도출한다. |
| | |
| | 차분 기법을 사용하지 않기 때문에, 정밀 지점 위치 결정에서는 미세한 물리적 현상에 의한 오차까지 정밀하게 모델링하여 보정해야 한다. 여기에는 위성과 수신기의 상대적 운동에 따른 [[상대성 이론]]적 효과, 위성의 회전에 의해 발생하는 [[위상 권선]](Phase Wind-up) 현상, 위성과 수신기 안테나의 위상 중심 오차(Antenna Phase Center Offset, PCO) 및 변동(Variation, PCV)이 포함된다. 또한 지구의 물리적 변형인 [[고체 지구 조석]](Solid Earth Tide), [[해양 조석 하중]](Ocean Tide Loading), 그리고 [[극운동]](Pole Tide)에 의한 지각의 미세한 움직임까지 계산 모델에 반영되어야 센티미터 단위의 정확도에 도달할 수 있다. |
| | |
| | 정밀 지점 위치 결정 기법의 주요한 기술적 과제 중 하나는 수렴 시간(Convergence time)의 단축이다. 단일 수신기를 사용하므로 미지 정수가 안정적으로 결정될 때까지 통상 수십 분 이상의 연속적인 관측 데이터가 요구된다. 최근에는 이를 해결하기 위해 위성에서 발생하는 위상 바이어스(Phase Bias) 정보를 추가로 제공하여 미지 정수를 정수값으로 고정하는 PPP-AR(Ambiguity Resolution) 기술이 발전하고 있다. 이러한 기술적 진보는 정밀 지점 위치 결정을 실시간 정밀 항법과 [[자율 주행]], [[정밀 농업]] 등 고도의 신뢰성이 요구되는 민간 산업 분야로 확산시키는 동력이 되고 있다. ((IGS PPP Products, http://www.acc.igs.org/igsacc_ppp.html |
| | )) |
| |
| ===== 주요 위성 항법 시스템의 종류 ===== | ===== 주요 위성 항법 시스템의 종류 ===== |
