사용자 도구

사이트 도구


정적_측위

문서의 이전 판입니다!


정적 측위

정적 측위의 개념과 정의

정적 측위(Static Positioning)는 범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 활용한 측위 방식 중 하나로, 관측 기간 동안 수신용 안테나를 지면에 고정된 특정 지점에 위치시켜 위성 신호를 수집하는 기법을 의미한다. 이는 수신기가 이동하며 위치를 산출하는 동적 측위(Kinematic Positioning)와 대비되는 개념이다. 현대 측량학에서 정적 측위는 가장 높은 수준의 정밀도를 확보할 수 있는 방법론으로 간주되며, 주로 국가 기준점 설치, 지각 변동 연구, 대규모 토목 구조물의 변위 측정 등 밀리미터(mm) 단위의 정확도가 요구되는 분야에서 핵심적인 역할을 수행한다.

정적 측위의 본질은 충분한 관측 시간을 확보함으로써 위성 신호에 포함된 각종 오차 성분을 통계적으로 처리하고 상쇄하는 데 있다. 위성으로부터 수신되는 신호는 전리층대류권 지연, 위성 궤도 오차, 시계 오차 등 다양한 물리적 요인에 의해 왜곡된다. 정적 측위에서는 이러한 오차를 효과적으로 제거하기 위해 주로 코드 데이터보다 정밀한 반송파 위상(Carrier Phase) 관측값을 사용하며, 두 대 이상의 수신기를 동시에 운용하는 상대 측위(Relative Positioning) 방식을 취한다. 이때 형성되는 두 수신기 간의 좌표 차이 벡터를 기선(Baseline)이라 하며, 후처리 과정을 통한 기선 해석을 통해 두 지점 사이의 정밀한 상대적 위치 관계가 규명된다.

정적 측위의 정밀도는 관측 시간 $ T $와 위성의 기하학적 배치 상태인 DOP(Dilution of Precision)에 밀접하게 의존한다. 관측 데이터의 수가 증가함에 따라 무작위 오차(Random Error)는 통계적 수렴 특성에 의해 감소하게 된다. 일반적으로 기선 벡터의 정밀도 $ $는 수신기 성능과 기선 길이에 따라 결정되며, 다음과 같은 형태의 실험식으로 표현되는 경우가 많다.

$$ \sigma = \sqrt{a^2 + (b \times L)^2} $$

위 식에서 $ a $는 수신기 및 안테나의 기계적 특성에 따른 고정 오차를 의미하며, $ b $는 기선 길이에 비례하여 증가하는 오차 계수, $ L $은 두 관측점 사이의 기선 길이를 나타낸다. 정적 측위는 짧게는 수십 분에서 길게는 수일에 걸친 장기 관측을 수행함으로써, 이러한 오차 항들을 최소화하고 최소제곱법(Least Squares Method) 등의 수학적 최적화 기법을 통해 신뢰도 높은 좌표 해를 도출한다.

결과적으로 정적 측위는 단순히 고정된 지점의 좌표를 얻는 행위를 넘어, 지구 물리적 현상을 정량화하는 고정밀 측정 도구로서의 가치를 지닌다. 이는 측지계의 구축과 유지에 필수적인 데이터 소스를 제공하며, 위성 신호의 미세한 변화를 추적하여 판 구조론적 이동이나 지반 침하를 감시하는 등 현대 우주 측지학의 근간을 이루는 기술적 토대가 된다. 수신기를 고정함으로써 얻어지는 데이터의 중첩 효과는 정적 측위만이 가지는 고유한 강점이며, 이는 실시간성이 강조되는 다른 측위 기법들이 도달하기 어려운 극도의 정밀도를 보장하는 원천이 된다.

정적 측위의 학술적 정의

정적 측위(Static Positioning)는 범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 활용한 측량 방법론 중 최고 수준의 정밀도를 확보할 수 있는 기법으로 정의된다. 학술적으로 정적 측위는 미지의 좌표를 결정하고자 하는 지점에 GNSS 수신 안테나를 고정된 상태로 설치하고, 충분한 시간 동안 위성 신호를 연속적으로 수신하여 누적된 시계열 데이터를 후처리(Post-processing) 또는 실시간 연산을 통해 해석하는 과정을 의미한다. 이 방식의 본질은 안테나를 공간적으로 고정함으로써 수신기의 이동에 따른 동역학적 변수를 제거하고, 오직 시간의 흐름에 따른 위성 신호의 변화만을 관측 데이터로 확보하는 데 있다.

정적 측위의 학술적 토대는 반송파 위상(Carrier Phase) 관측값의 정밀한 해석에 있다. 일반적인 코드 기반 측위가 수 미터 단위의 오차를 허용하는 것과 달리, 정적 측위는 위성 신호의 반송파 파장을 직접 측정하여 밀리미터(mm) 단위의 분해능(resolution)을 추구한다. 관측 과정에서 발생하는 전리층대류권 지연 오차, 위성의 궤도 및 시계 오차 등은 단일 시점의 관측으로는 완전히 제거하기 어려우나, 정적 측위에서는 장시간의 관측을 통해 이러한 오차 성분을 통계적으로 모델링하거나 소거한다. 특히 두 대 이상의 수신기를 동시에 운용하는 상대 측위(Relative Positioning) 방식을 취할 경우, 공통 오차를 제거하는 차분(Differencing) 기술을 적용하여 기선 벡터의 정밀도를 극대화할 수 있다.

수학적 관점에서 정적 측위는 관측 방정식의 중복성을 높여 최적의 해를 도출하는 과정이다. 관측 시간 $t$ 동안 수집된 방대한 데이터는 최소제곱법(Least Squares Method)이나 칼만 필터(Kalman Filter)와 같은 추정 알고리즘을 통해 처리된다. 이때 미지수로 설정되는 수신기의 위치 좌표 $(x, y, z)$는 시간이 경과함에 따라 관측 방정식의 수가 미지수의 수보다 많아지는 과결정 상태가 되며, 이는 결과적으로 해의 신뢰도와 정밀도를 높이는 결과로 이어진다. 특히 반송파 위상 측정에서 발생하는 정수 모호성(Integer Ambiguity)을 정확히 결정하기 위해서는 위성의 기하학적 배치 상태인 정밀도 저하율(Dilution of Precision, DOP)이 변화할 수 있는 충분한 관측 시간이 필수적으로 요구된다.

이러한 특성으로 인해 정적 측위는 측지학 및 지각 변동 연구에서 핵심적인 역할을 수행한다. 국가 기준점의 설치나 좌표계의 정의와 같이 절대적인 위치 정확도가 요구되는 분야에서는 정적 측위가 표준적인 방법으로 채택된다. 또한, 수 밀리미터 단위의 미세한 움직임을 추적해야 하는 판 구조론 기반의 지각 이동 관측이나 대형 구조물의 변위 모니터링에서도 정적 측위는 시공간적 제약을 극복하고 고정밀 데이터를 제공하는 유일한 수단으로 간주된다. 즉, 정적 측위는 단순히 정지 상태에서의 측정을 넘어, 시간적 적분을 통해 공간적 불확실성을 최소화하는 고도의 통계적 측위 기법이다.

정적 측위와 동적 측위의 비교

범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 이용한 위치 결정 방식은 수신기의 운동 상태와 데이터 처리 목적에 따라 크게 정적 측위와 동적 측위(Kinematic Positioning)로 구분된다. 두 방식의 근본적인 차이는 관측 시간 동안 안테나의 좌표 고정 여부에 있으며, 이는 오차 제거 기법, 데이터 처리 알고리즘, 그리고 최종적으로 도출되는 좌표의 정밀도에 결정적인 영향을 미친다. 정적 측위가 고정된 지점에서 장시간 관측을 통해 지구 중심 좌표계상의 절대적 위치를 정밀하게 규명하는 데 목적이 있다면, 동적 측위는 이동하는 수신기의 궤적을 실시간 또는 사후에 추적하는 데 중점을 둔다.

정적 측위는 미지의 지점에 안테나를 고정하고 위성의 기하학적 배치가 충분히 변화할 때까지 반송파 위상(Carrier Phase) 관측값을 수집한다. 이 방식은 확률론적 관점에서 충분한 표본 데이터를 확보함으로써 다중 경로 오차(Multipath Error)나 대기 지연과 같은 불규칙 오차 성분을 통계적으로 상쇄하는 효과를 가진다. 특히 두 대 이상의 수신기를 동시에 운용하는 상대 측위 방식을 적용할 경우, 위성과 수신기 사이의 공통 오차를 제거하여 수 mm에서 수 cm 수준의 고정밀 기선 해를 도출할 수 있다. 이는 국가 기준점 구축이나 지각 변동 감시와 같이 미세한 변위를 탐지해야 하는 분야에서 표준적인 방법론으로 자리 잡고 있다.

반면 동적 측위는 수신기가 이동하는 상태에서 연속적으로 위치를 결정하는 방식이다. 동적 측위의 핵심은 이동 시작 전 또는 이동 중에 정수 모호정(Integer Ambiguity)을 신속하게 결정하는 초기화(Initialization) 과정에 있다. 대표적인 기법인 실시간 동적 측위(Real-Time Kinematic, RTK)는 기준국으로부터 보정 정보를 실시간으로 수신하여 수 cm 급의 정밀도를 확보한다. 그러나 동적 측위는 정적 측위에 비해 관측 시간이 극히 짧기 때문에 위성 신호의 단절이나 급격한 주변 환경 변화에 취약하다. 특히 수신기가 교량 아래나 고층 건물 밀집 지역을 통과할 때 발생하는 사이클 슬립(Cycle Slip) 현상은 정수 모호정의 재결정을 요구하며, 이는 측위의 연속성과 신뢰성을 저하시키는 주요 원인이 된다.

두 방식의 정밀도 차이는 수학적으로 관측 방정식의 중복도와 관련이 있다. 정적 측위는 관측 시간 $T$ 동안 수집된 수많은 관측 방정식을 최소제곱법(Least Squares Method)으로 처리하여 미지수를 산출하므로, 시간적 중복성이 정밀도를 향상시키는 요인이 된다. 반면 동적 측위는 각 시점(Epoch)마다 위치를 추정해야 하므로, 주로 칼만 필터(Kalman Filter)와 같은 재귀적 추정 알고리즘을 사용하여 이전 시점의 정보를 바탕으로 현재의 위치를 갱신한다. 측위 정밀도 $\sigma$는 일반적으로 다음과 같은 관계를 갖는다.

$$ \sigma = \text{DOP} \times \sigma_{\text{obs}} / \sqrt{n} $$

위 식에서 정밀도 저하율(Dilution of Precision, DOP)은 위성의 배치 상태를 나타내며, $\sigma_{\text{obs}}$는 관측값의 표준편차, $n$은 독립적인 관측 횟수를 의미한다. 정적 측위는 관측 횟수 $n$을 극대화하여 $\sigma$를 최소화하는 전략을 취하는 반면, 동적 측위는 $n$이 제한적인 상황에서 최적의 DOP를 확보하고 알고리즘을 통해 $\sigma_{\text{obs}}$를 제어하는 데 집중한다.

응용 분야 측면에서 정적 측위는 높은 정밀도와 신뢰성이 요구되는 지적 측량, 공공 측량의 기준점 설치, 대형 구조물의 안전 진단 등에 활용된다. 데이터 처리 또한 관측 후 정밀 궤도력(Precise Ephemeris)을 활용하는 후처리(Post-processing) 방식이 일반적이다. 이와 대조적으로 동적 측위는 실시간성이 강조되는 자율주행 자동차의 항법, 무인 항공기(Unmanned Aerial Vehicle, UAV)의 경로 제어, 그리고 신속한 지형 정보 수집이 필요한 일반 측량 분야에서 광범위하게 사용된다. 결국 정적 측위와 동적 측위는 상호 배타적인 기술이 아니라, 요구되는 정밀도와 작업의 효율성 사이의 균형점에 따라 선택되는 상호 보완적인 체계라고 할 수 있다.

정밀도와 관측 시간의 상관관계

정적 측위에서 관측 시간은 산출되는 좌표의 신뢰도와 정밀도를 결정짓는 가장 핵심적인 변수 중 하나이다. 위성으로부터 수신되는 신호는 진공이 아닌 대기권을 통과하며 다양한 굴절과 지연을 겪고, 수신기 주변의 지형지물에 의한 간섭을 받게 된다. 이러한 오차 요인들은 시간에 따라 변하는 확률적 특성을 지니므로, 관측 시간을 충분히 확보함으로써 통계적 기법을 통해 이를 효과적으로 상쇄할 수 있다. 이는 오차론의 관점에서 데이터의 표본 크기를 늘려 모수에 대한 추정치의 정확도를 높이는 과정과 맥을 같이 한다.

통계적 원리에 따르면, 관측 데이터에 포함된 무작위 오차(Random error)는 대수의 법칙에 의해 관측 횟수가 증가할수록 그 평균값이 0으로 수렴하는 경향을 보인다. 정적 측위의 정밀도 $ $와 관측 시간 $ T $ 사이의 관계는 일반적으로 다음과 같은 반비례 모델로 설명된다.

$$ \sigma \approx \frac{k}{\sqrt{T}} $$

여기서 $ k $는 수신기의 성능, 위성의 배치 상태, 주변 환경 등에 의해 결정되는 상수이다. 위 식은 관측 시간이 길어질수록 정밀도가 향상됨을 수학적으로 보여주지만, 동시에 정밀도 향상의 효율이 시간의 제곱근에 비례하여 감소한다는 점도 시사한다. 즉, 초기 단계에서는 관측 시간의 짧은 연장만으로도 정밀도가 급격히 개선되나, 일정 수준 이상의 정밀도에 도달한 후에는 추가적인 관측이 주는 이득이 점차 줄어들게 된다.

특히 다중 경로(Multipath) 오차의 경우, 위성과 수신기, 그리고 반사체 사이의 기하학적 관계에 의존하는 특성을 가진다. 범지구 위성 항법 시스템(GNSS) 위성은 지구를 주회하며 끊임없이 위치가 변하므로, 장시간 관측을 수행하면 다중 경로 신호의 입사각과 위상이 변화하게 된다. 이러한 변화는 다중 경로 오차를 백색 잡음(White noise)과 유사한 형태로 분산시켜, 최소제곱법을 이용한 데이터 처리 과정에서 결과값에 미치는 영향을 최소화한다. 또한 전리층 지연이나 대류권 지연과 같은 대기 오차 역시 짧은 주기적 변동성을 내포하고 있어, 장기 관측을 통해 이러한 변동 성분을 평활화(Smoothing)하는 효과를 얻을 수 있다.

기선(Baseline)의 길이에 따른 정밀도 확보 전략에서도 관측 시간은 결정적인 역할을 한다. 기선의 길이가 길어질수록 두 수신기가 공유하는 대기 상태의 공통성이 낮아지기 때문에, 상대 측위를 통해 제거되지 않는 잔차 오차가 증가한다. 이를 극복하기 위해서는 더 많은 양의 관측 데이터를 확보하여 정수 모호정(Ambiguity)을 안정적으로 결정해야 한다. 일반적으로 수 킬로미터 이내의 단기선 측량에서는 수십 분 내외의 관측으로도 충분한 정밀도를 얻을 수 있으나, 수십에서 수백 킬로미터에 이르는 장기선 측량이나 국가 기준점 구축을 위한 정밀 측량에서는 수 시간에서 수일 이상의 연속 관측이 요구되기도 한다1).

결론적으로 정적 측위에서 관측 시간의 증가는 단순히 데이터의 양을 늘리는 것이 아니라, 시간에 따라 변화하는 다양한 물리적 오차 성분을 통계적으로 여과하는 과정이다. 이를 통해 위성 배치 기하학(Satellite Geometry)의 변화를 충분히 반영하고 신호의 신호 대 잡음비(SNR)를 실질적으로 개선함으로써, 밀리미터(mm) 단위의 고정밀 좌표 산출이 가능해진다.

정적 측위의 역사와 발전 과정

위성 측량의 역사는 1960년대 초반 미국 해군이 주도하여 개발한 트랜싯(Transit) 시스템, 일명 해군 위성 항법 시스템(Navy Navigation Satellite System, NNSS)의 구축과 함께 본격적으로 시작되었다. 초기 위성 측위의 핵심 원리는 위성이 이동함에 따라 발생하는 신호의 주파수 변화인 도플러 효과(Doppler effect)를 관측하는 것이었다. 지상에 고정된 수신기는 단일 위성이 궤도를 따라 이동하며 방출하는 전파의 주파수 이동량을 측정함으로써 위성과 수신기 사이의 상대적 거리를 계산하였다. 당시의 정적 측위는 위성이 관측 지점 상공을 통과하는 약 10분에서 15분 내외의 시간 동안 데이터를 수집하였으며, 수 차례의 위성 통과를 반복 관측하여 위치를 결정하였다. 이러한 초기 방식은 관측 시간이 길고 실시간 처리가 불가능하다는 한계가 있었으나, 장시간의 누적 데이터를 통해 수 미터 단위의 정밀도를 확보함으로써 지오이드(Geoid) 모델 구축 및 국가 간 기준망(Reference Network) 연결에 획기적인 전기를 마련하였다.

1970년대 후반부터 전개된 범지구 위성 항법 시스템(Global Positioning System, GPS)의 등장은 정적 측위의 패러다임을 근본적으로 변화시켰다. 미국 국방부가 개발한 GPS는 기존의 도플러 방식에서 탈피하여 위성으로부터의 신호 도달 시간을 측정하는 의사거리(Pseudorange) 관측과 반송파 위상(Carrier Phase) 관측을 도입하였다. 특히 정적 측위의 정밀도를 비약적으로 향상시킨 것은 반송파 위상 관측법의 발전이었다. 이는 파장이 수십 센티미터에 불과한 반송파 자체를 정밀한 자(scale)로 활용하는 방식으로, 정적 관측 시 밀리미터(mm) 단위의 정밀한 위치 결정이 가능해졌다. 이 시기에는 두 대 이상의 수신기를 동시에 운용하여 공통적인 오차 요인을 상쇄하는 상대 측위(Relative Positioning) 기법이 정립되었으며, 이는 현대 측량학에서 고정밀 좌표를 얻기 위한 표준적인 정적 측위 방법론으로 자리 잡았다.

21세기에 접어들어 정적 측위는 미국의 GPS에 의존하던 단일 시스템 체계를 넘어, 다양한 국가의 위성군을 통합 활용하는 범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS) 시대로 진입하였다. 러시아의 글로나스(GLONASS), 유럽 연합의 갈릴레오(Galileo), 중국의 베이두(BeiDou) 등이 차례로 가동됨에 따라 관측 가능한 위성의 수가 비약적으로 증가하였다. 다중 위성군(Multi-Constellation)의 활용은 위성 배치 상태를 나타내는 지표인 정밀도 저하율(Dilution of Precision, DOP)을 획기적으로 개선하였으며, 신호 차단이 빈번한 도심지나 산악 지형에서도 안정적인 정적 측위를 수행할 수 있는 환경을 조성하였다. 또한, 국제 GNSS 서비스(International GNSS Service, IGS)와 같은 국제 협력 기구가 제공하는 정밀 궤도 정보와 시계 보정 데이터는 초장기선(Very Long Baseline) 해석을 가능하게 하여, 지각 변동이나 대륙 이동과 같은 미세한 지구물리학적 변화를 밀리미터 단위로 추적하는 연구 분야의 발전을 이끌었다.

현대의 정적 측위 기술은 전통적인 기선 해석 방식을 넘어 정밀 단일 측위(Precise Point Positioning, PPP) 기술의 고도화로 이어지고 있다. PPP는 기준국과의 상대적 위치 관계를 규명하지 않고도 단일 수신기에서 수집된 관측값에 정밀 오차 모델을 직접 적용하여 고정밀 좌표를 산출하는 방식이다. 이는 광역적인 측위 환경에서 기준국 설치의 제약을 극복하게 해주며, 최근에는 인공지능(Artificial Intelligence)을 활용한 신호 필터링 기술과 다중 주파수(Multi-frequency) 신호 처리 기법이 결합되어 관측 시간 단축과 정밀도 향상을 동시에 꾀하고 있다. 이와 같은 발전 과정은 정적 측위가 단순한 지형 측량을 넘어 국가 위치 기준 체계의 유지, 거대 구조물의 안전 진단, 지구 환경 모니터링 등 정밀 좌표가 요구되는 현대 과학기술의 필수적인 기반 기술로 확립되는 과정이라 할 수 있다.

초기 위성 측량 시스템

위성 측량의 태동은 1950년대 후반 인공위성의 궤도 추적 과정에서 발견된 물리적 현상을 역으로 이용하려는 시도에서 비롯되었다. 1957년 소련의 스푸트니크 1호(Sputnik 1) 발사 당시, 미국의 존스 홉킨스 대학교 응용물리학연구소(Applied Physics Laboratory, APL) 연구진은 지상 수신기가 포착한 위성 신호의 주파수가 위성의 이동에 따라 변화하는 도플러 효과(Doppler Effect)를 관측하였다. 연구진은 위성의 정확한 궤도 역학 정보를 알고 있다면, 지상 수신 지점에서 관측된 도플러 천이(Doppler Shift)의 변화 양상을 분석하여 수신기의 역좌표를 산출할 수 있다는 점에 주목하였다. 이러한 발상은 세계 최초의 위성 항법 시스템인 트랜싯(Transit) 시스템, 즉 해군 위성 항법 시스템(Navy Navigation Satellite System, NNSS)의 개발로 이어졌다.

트랜싯 시스템의 핵심적인 운용 원리는 위성이 관측자의 상공을 통과할 때 발생하는 상대적 속도 변화를 이용하는 것이다. 위성이 수신기에 접근할 때는 수신 주파수가 송신 주파수보다 높게 측정되고, 멀어질 때는 낮게 측정된다. 이때 주파수의 변화율이 가장 가파른 지점은 위성이 수신기와 가장 가까운 거리(Closest Approach)에 위치한 시점이다. 이 이론적 배경 하에서 수신기는 위성으로부터 전송되는 궤도 정보인 천문력(Ephemeris)과 시간에 따른 주파수 측정값을 결합하여 자신의 위치를 결정한다. 초기 시스템은 전리층에 의한 신호 굴절 오차를 보정하기 위해 150MHz와 400MHz의 두 가지 주파수를 동시에 사용하였으며, 이는 현대 범지구 위성 항법 시스템(GNSS)의 다중 주파수 활용 모태가 되었다2).

정적 측위의 관점에서 초기 트랜싯 시스템은 현대적 시스템과 비교하여 뚜렷한 한계를 지니고 있었다. 우선, 트랜싯 위성은 고도 약 1,100km의 저궤도(Low Earth Orbit, LEO)를 선회하는 5~6기의 위성으로 구성되었기 때문에, 관측자가 다음 위성을 포착하기까지 짧게는 수십 분에서 길게는 몇 시간까지 대기해야 하는 불연속성의 문제가 발생하였다3). 또한, 한 번의 위성 통과(Pass) 시 수집되는 데이터만으로는 수 미터에서 수십 미터에 달하는 오차가 발생하였으므로, 정밀한 좌표를 얻기 위해서는 동일 지점에서 수일간 반복적으로 위성 신호를 수집하여 통계적으로 처리하는 정적 측위 방식이 필수적이었다.

특히 초기 시스템은 수신기의 이동 속도가 주파수 측정에 직접적인 오차를 유발하는 구조적 취약성을 지니고 있었다. 수신기가 이동할 경우 발생하는 도플러 효과가 위성의 이동에 의한 효과와 중첩되어 위치 산출의 정확도를 현저히 저하시켰기 때문이다. 이러한 이유로 트랜싯 시스템은 항행 중인 함정보다는 지면에 고정된 지점의 좌표를 결정하는 측지학적 목적이나 잠수함의 위치 교정 등에 주로 활용되었다. 비록 실시간 3차원 위치 결정은 불가능하였으나, 전 지구적 좌표계를 기반으로 대륙 간 기선 거리를 측정할 수 있는 토대를 마련하였다는 점에서 초기 위성 측량 시스템은 정적 측위 기술 발전의 중요한 분기점으로 평가된다.

범지구 위성 항법 시스템의 도입

초기 위성 측량의 주류였던 트랜싯(Transit) 시스템은 도플러 효과(Doppler Effect)를 이용한 위치 결정 방식을 채택하였으나, 위성 수가 부족하여 연속적인 관측이 불가능하고 고도 측정의 정확도가 현저히 떨어진다는 한계가 있었다. 이를 극복하기 위해 1970년대 초반 미국 국방부(Department of Defense, DoD)는 실시간으로 3차원 위치, 속도, 시간을 결정할 수 있는 전지구 위치 결정 시스템(Global Positioning System, GPS)의 개발에 착수하였다. GPS는 본래 군사적 목적으로 설계되었으나, 그 과정에서 도입된 고정밀 반송파 위상(Carrier Phase) 관측 기술은 정적 측위 분야에 혁신적인 변화를 가져왔다.

GPS의 도입으로 정적 측위는 기존의 도플러 주파수 편이 측정 방식에서 탈피하여, 위성에서 송신하는 반송파의 위상 차이를 직접 측정하는 방식으로 진화하였다. GPS 위성은 L1($1575.42 \text{ MHz}$)과 L2($1227.60 \text{ MHz}$) 대역의 주파수를 사용하여 신호를 송신하며, 수신기는 이 신호의 위상을 추적하여 위성과 수신기 사이의 거리를 측정한다. 정적 측위에서는 관측점의 위치가 고정되어 있다는 특성을 활용하여, 장시간 동안 축적된 반송파 위상 데이터를 통계적으로 처리함으로써 정수 모호정(Integer Ambiguity)을 해결하고 밀리미터 단위의 정밀도를 달성할 수 있게 되었다.

GPS를 필두로 한 현대적 위성 항법 시스템은 정적 측위에 상대 측위(Relative Positioning) 기법을 보편화시켰다. 이는 미지의 지점과 이미 좌표를 알고 있는 기준점(Reference Station)에서 동시에 위성 신호를 수신하여 그 차이를 분석하는 방식이다. 이 과정에서 발생하는 이중 차분(Double Difference) 관측 방정식은 다음과 같은 형태로 표현된다.

$$\Delta\nabla\Phi_{AB}^{jk} = \rho_{AB}^{jk} + \lambda\Delta\nabla N_{AB}^{jk} + \epsilon$$

여기서 $\Delta\nabla\Phi_{AB}^{jk}$는 두 수신기($A, B$)와 두 위성($j, k$) 사이의 이중 차분 위상 관측값이며, $\rho_{AB}^{jk}$는 기하학적 거리의 차이, $\lambda$는 반송파의 파장, $\Delta\nabla N_{AB}^{jk}$는 이중 차분된 정수 모호정, $\epsilon$은 잔여 오차를 의미한다. 이러한 수치적 접근을 통해 위성의 시계 오차(Clock Error)와 수신기의 시계 오차를 소거할 수 있게 되었으며, 이는 정적 측위의 신뢰성을 비약적으로 향상시키는 계기가 되었다.

또한 GPS의 도입은 전 지구적으로 통일된 측지계(Geodetic Datum)인 WGS84(World Geodetic System 1984)의 확립과 밀접하게 연계되었다. 과거의 정적 측위는 각 국가나 지역별로 서로 다른 타원체와 원점을 사용하여 성과 간의 호환성이 부족하였으나, GPS는 지구 중심을 원점으로 하는 좌표계를 제공함으로써 범지구적인 정적 측위 성과 통합을 가능하게 하였다. 이는 국제 측지학 및 지구물리학 연맹(International Union of Geodesy and Geophysics, IUGG)과 같은 기구를 통해 국제 지구 기준 좌표계(International Terrestrial Reference System, ITRS)가 고도화되는 물리적 토대가 되었다.

현대적 GNSS의 도입은 단순히 기술적 수단의 변화를 넘어 정적 측위의 패러다임을 효율성 중심으로 재편하였다. 과거 수일이 소요되던 고정밀 기선 해석 작업은 GPS의 다중 위성 배치와 고도화된 후처리 소프트웨어의 등장으로 수 시간 이내로 단축되었다. 이러한 변화는 국가 기준점 체계의 정비, 지각 변동의 정밀 감시, 그리고 대규모 토목 공사의 정밀 제어 측량 등 현대 측량학의 전 분야에 걸쳐 핵심적인 동력으로 작용하고 있다. 4) 5)

다중 위성군 활용 기술의 진화

초기 범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 활용한 정적 측위는 미국의 GPS(Global Positioning System)에 전적으로 의존하는 단일 위성군 체계였다. 그러나 2010년대 이후 러시아의 글로나스(GLONASS)의 현대화, 유럽연합의 갈릴레오(Galileo) 및 중국의 베이두(BeiDou) 시스템의 전 지구적 서비스 개시로 인해 현대 정적 측위는 다중 위성군(Multi-GNSS)을 통합 활용하는 방향으로 진화하였다. 이러한 변화는 가용 위성 수의 비약적인 증가를 불러왔으며, 이는 곧 정적 측위의 정밀도와 신뢰성, 그리고 관측 효율성의 향상으로 직결되었다.

다중 위성군 환경에서 정적 측위의 가장 큰 기술적 이점은 정밀도 저하율(Dilution of Precision, DOP)의 개선이다. 단일 위성군만을 사용할 때보다 수신기가 포착할 수 있는 위성의 개수가 2배에서 3배 이상 증가함에 따라, 위성 배치에 따른 기하학적 오차가 최소화된다. 특히 고층 건물이 밀집한 도심지나 깊은 계곡과 같은 신호 차폐 지역에서도 충분한 수의 위성을 확보할 수 있게 되어, 과거에는 정적 측위가 불가능했던 지역에서도 정밀한 좌표 산출이 가능해졌다. 또한, 서로 다른 궤도 경사각을 가진 위성군들이 결합되면서 관측 데이터의 중복성(Redundancy)이 높아져, 신호 잡음이나 다중 경로 오차에 대한 저항력도 강화되었다.

기술적으로 다중 위성군을 통합하기 위해서는 각 시스템이 사용하는 상이한 좌표계(Coordinate System)와 시간 체계(Time System)를 정밀하게 결합하는 과정이 필수적이다. GPS는 WGS84를 기준 좌표계로 사용하는 반면, 글로나스는 PZ-90, 갈릴레오는 GTRF(Galileo Terrestrial Reference Frame), 베이두는 CGCS2000을 사용한다. 이들 간의 미세한 차이를 보정하기 위해 국제 지구 자전 및 좌표 시스템 서비스(International Earth Rotation and Reference Systems Service, IERS)에서 정의한 변환 매개변수가 적용된다. 시간 체계 역시 원자시를 기반으로 하지만 각 시스템의 기원(Epoch)과 윤초 처리 방식이 다르므로, 이를 단일한 기준 시계로 환산하는 알고리즘이 현대 수신기와 후처리 소프트웨어의 핵심 요소가 되었다.

또한, 다중 위성군 활용은 반송파 위상 관측의 핵심 과제인 정수 모호정(Integer Ambiguity) 결정 시간을 단축시킨다. 다양한 주파수 대역의 신호를 복합적으로 활용함으로써, 모호정을 빠르게 확정할 수 있는 와이드레인(Wide-lane) 조합의 생성이 용이해지기 때문이다. 특히 갈릴레오와 베이두가 제공하는 3중 주파수(Triple-frequency) 신호는 기존의 2중 주파수 체계보다 오차 제거 능력이 탁월하여, 장기선 정적 측위에서도 밀리미터 단위의 정밀도를 신속하게 확보할 수 있게 한다. 이러한 기술적 진화는 정밀 지점 측위(Precise Point Positioning, PPP) 기법과 결합되어, 기준국 없이 단일 수신기만으로도 고정밀 좌표를 얻는 차세대 정적 측위 모델의 기반이 되었다.

현대적 정적 측위 기술은 단순히 위성 수를 늘리는 수준을 넘어, 각 위성군의 신호 특성을 고려한 가중치 부여 모델과 시스템 간 편차(Inter-System Bias, ISB) 보정 기술을 포함한다. 서로 다른 하드웨어 특성으로 인해 발생하는 신호 지연의 차이를 수학적으로 모델링하여 제거함으로써, 이질적인 위성군들이 마치 하나의 거대한 가상 위성군처럼 작동하도록 하는 것이다. 이러한 통합 활용 기술의 진화는 지각 변동 모니터링, 국가 기준점 체계 유지 등 초정밀 위치 정보가 요구되는 현대 사회의 핵심 인프라 기술로 자리 잡고 있다.

정적 측위의 물리적 및 수학적 원리

정적 측위(Static Positioning)는 위성 항법 시스템의 신호 중 반송파(Carrier Wave)의 위상을 주된 관측량으로 활용하여, 지표면에 고정된 수신기의 정밀 좌표를 결정하는 기술이다. 이 방식의 물리적 기초는 위성에서 송신된 전자기파의 위상과 수신기 내부에서 생성된 기준 신호 사이의 위상차를 측정하는 데 있다. 코드 측위(Code-based Positioning)가 신호의 시간 지연을 이용해 미터 단위의 오차를 갖는 것과 대조적으로, 반송파 위상 관측은 파장의 수 밀리미터 수준까지 분해능을 확장할 수 있어 고정밀 측량의 핵심이 된다.

반송파 위상 관측의 수학적 모델은 수신기 $ i $가 위성 $ j $로부터 수신한 위상 측정치 $ _i^j $를 기하학적 거리와 각종 오차 요인의 합으로 기술한다. 관측 방정식은 일반적으로 다음과 같이 정의된다.

$$ \phi_i^j = \frac{1}{\lambda} \rho_i^j + N_i^j + f \cdot (\delta t_i - \delta T^j) - I_i^j + T_i^j + \epsilon_i^j $$

여기서 $ $는 반송파의 파장, $ _i^j $는 위성과 수신기 사이의 실제 기하학적 거리, $ N_i^j $는 정수 모호정(Integer Ambiguity)이다. $ t_i $와 $ T^j $는 각각 수신기와 위성의 시계 오차를 나타내며, $ I_i^j $와 $ T_i^j $는 전리층(Ionosphere) 및 대류권(Troposphere) 지연 효과를 의미한다. 이러한 방정식은 위성 항법의 기하학적 원리와 신호 전파 특성을 포괄하는 기초가 된다.6)

정적 측위에서 가장 중요한 수학적 난제는 정수 모호정의 결정이다. 수신기가 위성 신호를 처음 잠금(lock-on)하는 순간, 수신기는 위상차의 소수 부분은 측정할 수 있으나 위성과 수신기 사이에 존재하는 전체 파장의 개수 $ N $은 알 수 없다. 이 미지수는 정수 값을 가지며, 관측이 중단되지 않는 한 일정하게 유지된다. 정적 측위에서는 장시간의 관측 데이터를 최소제곱법(Least Squares Method)이나 칼만 필터(Kalman Filter)와 같은 추정 기법으로 처리하여 이 정수 값을 추정하며, 모호정이 정확히 결정되었을 때 비로소 밀리미터 단위의 정밀도를 얻을 수 있다.

고정밀 좌표 산출을 위해 정적 측위는 단일 수신기를 이용하는 방식보다 두 대 이상의 수신기를 동시에 운용하는 상대 측위(Relative Positioning) 기법을 주로 사용한다. 이는 수신기 간의 상대적인 위치 관계를 나타내는 기선(Baseline) 벡터를 해석하는 과정이다. 이때 공통적으로 발생하는 오차를 제거하기 위해 차분(Differencing) 기술이 적용된다.

단일 차분(Single Difference)은 두 수신기가 하나의 위성을 동시에 관측할 때 발생하며, 이를 통해 위성 시계 오차와 궤도 오차를 상당 부분 상쇄한다. 이중 차분(Double Difference)은 두 수신기가 두 개의 위성을 관측한 단일 차분 값들의 차이를 구하는 방식이다. 이중 차분 과정을 거치면 수신기 시계 오차가 제거되며, 관측 방정식에서 정수 모호정 항만이 남게 되어 수학적 모델이 단순화된다.

$$ \nabla \Delta \phi_{AB}^{pq} = (\phi_B^q - \phi_A^q) - (\phi_B^p - \phi_A^p) $$

위 식에서 $ $는 이중 차분 연산자를 의미하며, 수신기 A, B와 위성 p, q 사이의 위상차를 결합하여 위치 매개변수를 산출한다. 이러한 수치적 결합은 정밀한 기선 해석의 핵심적인 도구가 된다.7)

마지막으로 삼중 차분(Triple Difference)은 서로 다른 두 시점 사이의 이중 차분 값의 변화를 분석한다. 삼중 차분은 정수 모호정 항을 소거하므로, 관측 도중 위성 신호가 일시적으로 차단되어 발생하는 사이클 슬립(Cycle Slip) 현상을 탐지하고 수정하는 데 매우 효과적이다. 또한 정수 모호정의 초기 근사치를 결정하는 단계에서 중요한 역할을 수행한다. 이러한 물리적·수학적 원리들이 결합하여 정적 측위는 지각 변동 감시나 국가 기준점 측량과 같은 초정밀 분야에서 신뢰할 수 있는 위치 정보를 제공한다.

반송파 위상 관측 원리

반송파 위상 관측은 범지구 위성 항법 시스템(GNSS)에서 제공하는 반송파(Carrier Wave)의 위상을 측정하여 위성과 수신기 사이의 거리를 결정하는 기법이다. 이는 코드 기반의 의사거리(Pseudorange) 측정 방식보다 훨씬 높은 정밀도를 제공하는데, 그 이유는 신호의 해상도 차이에 있다. 일반적으로 코드 기반 측위에서 사용하는 C/A 코드(Coarse/Acquisition Code)의 칩 길이는 약 300m에 달하는 반면, L1 반송파의 파장은 약 19cm, L2 반송파는 약 24cm에 불과하다. 수신기는 이 파장의 수백 분의 일에 해당하는 위상 변화를 감지할 수 있으므로, 이론적으로 밀리미터 단위의 거리 측정이 가능해진다.

반송파 위상 관측의 기본 원리는 위성에서 송신된 반송파의 위상과 수신기 내부에서 생성된 복제 반송파의 위상 차이를 측정하는 것이다. 위성에서 방출된 신호가 수신기에 도달하면, 수신기는 도달한 신호의 위상이 한 주기 내에서 어느 지점에 있는지를 판별한다. 그러나 위성에서 수신기까지의 전체 거리는 단순히 이 위상 차이만으로 결정되지 않는다. 위성과 수신기 사이에는 수많은 파장의 정수 배가 존재하는데, 수신기가 신호를 처음 추적하기 시작한 시점에서는 이 정수 배의 값을 알 수 없다. 이를 정수 모호정(Integer Ambiguity)이라 하며, 반송파 위상 관측을 통한 정밀 측위의 핵심 과제는 이 모호정 값을 정확히 결정하는 데 있다.

수학적으로 시각 $t$에서의 반송파 위상 관측 방정식은 다음과 같이 표현할 수 있다.

$$\phi(t) = \frac{1}{\lambda} \rho + N + f \cdot (\delta t_r - \delta t_s) - I + T + \epsilon$$

여기서 $\phi(t)$는 수신기가 측정한 위상값(주기 단위), $\lambda$는 반송파의 파장, $\rho$는 위성과 수신기 사이의 실제 기하학적 거리이다. $N$은 미지의 정수 모호정이며, $f$는 신호의 주파수, $\delta t_r$과 $\delta t_s$는 각각 수신기와 위성의 시계 오차를 나타낸다. 또한 $I$와 $T$는 각각 전리층 지연대류권 지연에 의한 오차 항이며, $\epsilon$은 다중 경로 오차 및 수신기 잡음 등을 포함한 잔차 항이다. 이 방정식에서 알 수 있듯이, 관측된 위상값은 기하학적 거리뿐만 아니라 정수 모호정과 다양한 오차 요인이 결합된 형태이다.

정적 측위에서는 수신기가 고정된 위치에서 장시간 데이터를 수집하므로, 위성의 궤도 운동에 따른 기하학적 배치의 변화를 이용하여 정수 모호정을 통계적으로 추정할 수 있다. 관측이 진행되는 동안 수신기가 신호를 놓치지 않고 연속적으로 추적한다면 $N$값은 상수로 유지된다. 그러나 물리적 장애물이나 신호 간섭으로 인해 신호 추적이 일시적으로 중단될 경우, 정수 배의 연속성이 깨지는 사이클 슬립(Cycle Slip) 현상이 발생한다. 이 경우 수신기는 다시 모호정을 결정해야 하며, 이는 측위 정밀도를 저하시키는 주요 원인이 된다.

결과적으로 반송파 위상 관측 원리는 극히 짧은 파장을 가진 전자기파를 거리 측정의 척도로 삼음으로써 미시경제학적 관점에서의 정밀한 계측과 유사하게 오차의 범위를 극소화한다. 정적 측위는 이러한 반송파의 특성을 극대화하기 위해 장시간의 관측을 수행하며, 이를 통해 정수 모호정을 확정하고 대기 지연 등 확률적으로 변동하는 오차들을 상쇄시킨다. 이러한 정밀 관측 원리는 측지학적 기준망 설정이나 지각 변동의 미세한 흐름을 관측하는 데 있어 필수적인 이론적 토대가 된다.

반송파의 파장과 분해능

정적 측위에서 고정밀 좌표를 산출하기 위해 반송파(Carrier Wave)를 주요 관측량으로 활용하는 근본적인 이유는 그 물리적 특성에 기인한 높은 해상도(Resolution)에 있다. 범지구 위성 항법 시스템(GNSS)에서 송신하는 신호는 크게 L 밴드(L-band) 전자기파에 실린 코드(Code) 데이터와 그 정보를 실어 나르는 반송파로 구분된다. 코드 기반 측위가 수백 미터에서 수십 미터에 달하는 칩(Chip) 길이에 의존하여 수 미터 수준의 오차를 허용하는 반면, 반송파 위상 관측은 전자기파의 위상 자체를 측정함으로써 밀리미터 단위의 정밀도를 확보할 수 있는 물리적 토대를 제공한다.

GPS(Global Positioning System)를 기준으로 할 때, 주요 반송파 주파수는 $L_1$($1575.42 \text{ MHz}$), $L_2$($1227.60 \text{ MHz}$), 그리고 현대화된 신호인 $L_5$($1176.45 \text{ MHz}$) 등으로 구성된다. 진공에서의 빛의 속도를 $c \approx 299,792,458 \text{ m/s}$라고 할 때, 각 주파수 $f$에 대응하는 파장 $\lambda$은 다음의 관계식에 의해 결정된다. $$ \lambda = \frac{c}{f} $$ 이 식에 따라 계산하면 $L_1$ 반송파의 파장은 약 $19.03 \text{ cm}$, $L_2$는 약 $24.42 \text{ cm}$, $L_5$는 약 $25.48 \text{ cm}$이다. 이러한 파장의 길이는 측위 시스템의 분해능을 결정하는 기초적인 척도가 된다. 일반적으로 고성능 GNSS 수신기는 반송파 위상의 약 1% 내외를 측정할 수 있는 능력을 갖추고 있으므로, 물리적으로 약 $2 \text{ mm}$ 수준의 관측 분해능을 달성할 수 있게 된다.

반송파의 주파수가 높을수록 파장은 짧아지며, 이는 이론적으로 더 정밀한 해상도를 의미한다. 그러나 실제 정적 측위의 정밀도는 단순히 파장의 길이에만 의존하지 않고, 전리층 지연과 같은 대기 환경 변수와의 상호작용에 의해 복합적으로 결정된다. 서로 다른 주파수를 가진 두 개 이상의 반송파(예: $L_1$과 $L_2$)를 동시에 사용하는 이중 주파수 관측은 전리층의 분산(Dispersion) 특성을 이용하여 대기 지연 오차를 효과적으로 제거할 수 있게 한다. 이는 정수 모호정(Integer Ambiguity)을 신속하고 정확하게 결정하는 데 기여하며, 결과적으로 기선 해석의 최종적인 정확도를 높이는 핵심 요인이 된다.

또한, 반송파의 신호 대 잡음비(Signal-to-Noise Ratio, SNR)는 위상 측정의 불확실성을 결정하는 중요한 변수이다. 신호의 강도가 충분하지 않거나 주변 지형지물에 의한 다중 경로(Multipath) 간섭이 발생할 경우, 반송파 위상의 미세한 변화를 감지하는 정밀도가 저하되어 분해능의 손실을 초래한다. 따라서 정적 측위에서는 파장이 짧은 고주파 대역의 신호를 안정적으로 수신할 수 있는 환경을 조성하고, 이를 수학적으로 처리하여 물리적 파장의 한계를 극복하는 과정이 필수적으로 요구된다. 이러한 반송파의 물리적 특성과 수신기의 측정 분해능 사이의 관계는 정적 측위가 지각 변동 감시나 정밀 기준점 측량과 같이 극도의 정확도를 요구하는 분야에서 표준 기술로 자리 잡게 된 핵심적인 근거가 된다. 8)

정수 모호정 결정 이론

반송파 위상(Carrier Phase) 관측을 이용한 정적 측위에서 가장 핵심적인 과제는 관측 시작 시점에서 수신기와 위성 사이의 전체 파장 수를 나타내는 미지수인 정수 모호정(Integer Ambiguity)을 정확히 결정하는 것이다. 범지구 위성 항법 시스템(GNSS) 수신기는 위성 신호의 소수점 이하 위상은 밀리미터 단위의 정밀도로 측정할 수 있으나, 신호가 방출되어 수신기에 도달하기까지 거쳐온 전체 파장의 개수 $N$은 직접적으로 측정할 수 없다. 이 $N$값이 정수라는 물리적 제약 조건을 활용하여 이를 확정하는 과정을 정수 모호정 결정이라 하며, 이 과정이 성공적으로 수행되어야만 비로소 센티미터 또는 밀리미터 수준의 고정해(Fixed Solution)를 얻을 수 있다.

정수 모호정 결정을 위한 수학적 모델은 일반적으로 이중 차분(Double Difference) 관측 방정식을 기초로 한다. 이중 차분을 통해 수신기 시계 오차와 위성 시계 오차를 제거하면, 위상 관측 방정식은 다음과 같이 정수 모호정 항을 포함하는 형태로 정리된다.

$ = + N + $

여기서 $ $는 이중 차분된 위상 관측량, $ $는 위성과 수신기 사이의 기하학적 거리 차이, $ $는 반송파의 파장, $ N $은 이중 차분된 정수 모호정, $ $은 잔여 오차를 의미한다. 정수 모호정 결정 이론은 크게 실수해 추정, 정수해 검색, 그리고 검정의 세 단계로 구성된다.

첫 번째 단계인 실수해(Float Solution) 추정에서는 $ N $이 정수라는 제약을 무시하고 이를 실수 변수로 취급하여 최소제곱법(Least Squares Method)이나 칼만 필터(Kalman Filter)를 통해 최적값을 산출한다. 이 과정에서 정수 모호정의 실수 추정치 $\hat{a}$와 그에 따른 공분산 행렬(Covariance Matrix) $Q_{\hat{a}}$가 함께 도출된다. 그러나 대기 지연이나 위성 궤도 오차 등의 영향으로 인해 이 실수해는 실제 정수값에서 벗어나게 된다.

두 번째 단계는 실수해 근처에서 실제 정수 조합을 찾는 정수해 검색(Integer Search) 단계이다. 초기에는 FARA(Fast Ambiguity Resolution Approach)와 같은 기법이 사용되었으나, 현재는 Teunissen에 의해 제안된 LAMBDA(Least-squares AMBiguity Decorrelation Adjustment) 기법이 표준적으로 사용된다9). GNSS 관측값은 위성 배치 전략상 변수 간의 상관관계가 매우 높기 때문에, 공분산 행렬의 타원체가 매우 길게 늘어진 형태를 띠어 검색 효율이 저하된다. LAMBDA 기법은 정수 특성을 보존하는 Z-변환(Z-transformation)을 통해 변수 간의 상관관계를 제거(Decorrelation)하고, 검색 공간을 구형에 가깝게 변형하여 최적의 정수 벡터를 신속하게 탐색한다. 이때 목적 함수는 다음과 같은 정수 최소제곱(Integer Least Squares) 문제를 해결하는 것이다.

$$ \min_{a \in \mathbb{Z}^n} (\hat{a} - a)^T Q_{\hat{a}}^{-1} (\hat{a} - a) $$

여기서 $a$는 후보 정수 벡터이며, 이 식을 최소화하는 $a$가 최적 정수해로 선택된다.

마지막 단계는 선택된 정수해의 통계적 신뢰성을 평가하는 검정(Validation) 단계이다. 가장 널리 사용되는 방법은 Ratio Test로, 검색 과정에서 발견된 가장 작은 오차 제곱합과 두 번째로 작은 오차 제곱합의 비율을 비교하는 것이다10).

$ Ratio = > k $

여기서 $ _1 $은 최적해의 잔차 제곱합, $ _2 $는 차선해의 잔차 제곱합이며, 임계치 $k$는 통계적 유의 수준에 따라 보통 2 또는 3으로 설정된다. 이 비율이 임계치를 초과할 경우에만 해당 정수해를 확정하고, 그렇지 못할 경우 실수해를 그대로 사용하거나 관측 시간을 연장하여 데이터를 추가로 확보해야 한다. 정수 모호정이 올바르게 결정되면 위상 관측값의 정밀도가 거리 단위로 환산되어 측량의 정확도가 비약적으로 향상된다.

기선 해석과 상대 측위

상대 측위(Relative Positioning)는 두 대 이상의 수신기를 동시에 운용하여 수신기 간의 상대적인 위치 관계를 결정하는 기법이다. 이 방식은 단일 수신기를 이용하는 절대 측위에서 발생하는 위성 궤도 오차, 시계 오차, 그리고 대기 지연 오차의 상당 부분이 인접한 두 수신기에서 공통적으로 나타난다는 점에 착안한다. 관측된 데이터에서 이러한 공통 오차 성분을 상쇄함으로써 정밀도를 획기적으로 높일 수 있으며, 이때 결정되는 두 지점 사이의 3차원 상대 좌표 차이 벡터를 기선(Baseline)이라 한다. 기선 해석(Baseline Processing)은 수집된 반송파 위상 관측 데이터를 수학적으로 처리하여 이 기선 벡터를 산출하는 일련의 과정을 의미한다.

기선 해석의 핵심은 이중 차분(Double Difference) 관측 방정식을 구성하여 미지수를 해결하는 데 있다. 두 대의 수신기 $A, B$가 두 개의 위성 $i, j$를 동시에 관측했을 때, 이중 차분 방정식은 다음과 같이 표현된다.

$$ \nabla\Delta\phi_{AB}^{ij} = \frac{1}{\lambda}\nabla\Delta\rho_{AB}^{ij} + \nabla\Delta N_{AB}^{ij} + \epsilon $$

여기서 $\nabla\Delta\phi_{AB}^{ij}$는 이중 차분된 위상 관측값이며, $\lambda$는 반송파의 파장, $\nabla\Delta\rho_{AB}^{ij}$는 위성과 수신기 사이의 기하학적 거리 차이의 조합이다. $\nabla\Delta N_{AB}^{ij}$는 이 과정에서 반드시 해결해야 할 정수 모호정(Integer Ambiguity)을 나타내며, $\epsilon$은 잔여 오차를 의미한다. 이 식을 통해 수신기 간의 상대적 좌표 증분량인 $(\Delta X, \Delta Y, \Delta Z)$를 미지수로 설정하고 최소제곱법(Least Squares Method)을 적용하여 최적해를 구한다.

기선 해석 과정은 크게 세 단계로 구분된다. 첫째, 정수해(Fixed Solution)를 구하기 전 단계로, 정수 모호정을 실수(float) 값으로 추정하는 실수해(Float Solution) 산출 단계이다. 둘째, 추정된 실수값으로부터 통계적 기법을 동원하여 가장 가능성 있는 정수값을 결정하는 모호정 결정 단계이다. 이때 LAMBDA 기법(Least-squares Ambiguity Decorrelation Adjustment)과 같은 알고리즘이 주로 사용된다. 셋째, 결정된 정수 모호정을 고정시킨 상태에서 최종적인 기선 벡터를 산출하는 고정해 산출 단계이다. 만약 관측 환경이 불량하여 정수 모호정이 정확히 결정되지 않을 경우 실수해를 최종 결과로 사용하기도 하지만, 이 경우 정밀도는 크게 저하된다.

산출된 기선 벡터의 품질은 분산-공분산 행렬(Variance-Covariance Matrix)과 RATIO 값 등을 통해 평가된다. RATIO 값은 가장 가능성 있는 정수 모호정 조합과 두 번째로 가능성 있는 조합 사이의 확률적 대비를 나타내는 지표로, 통상적으로 2.0 또는 3.0 이상의 값을 가질 때 정수해의 신뢰도가 높은 것으로 판단한다. 이렇게 결정된 개별 기선들은 이후 기망 조정(Network Adjustment) 과정을 거쳐 전체적인 기준점 체계 내에서 일관성을 확보하게 된다. 기선 해석은 단순히 거리와 방향을 구하는 것을 넘어, 지구 중심 좌표계ITRF(International Terrestrial Reference Frame) 또는 지역 좌표계상의 정밀한 위치 정보를 제공하는 정적 측위의 이론적 토대가 된다.

단일 차분과 이중 차분

수신기 간 또는 위성 간의 관측값 차이를 통해 공통 오차를 제거하는 수식을 고찰한다.

삼중 차분을 통한 초기치 결정

시간적 차이까지 고려하여 오차를 최소화하고 초기 근사치를 산출하는 과정을 다룬다.

오차의 원인과 보정 기술

정적 측위의 정밀도를 결정짓는 핵심 기제는 위성 신호의 생성, 전파, 수신 과정에서 개입하는 다양한 오차 요인을 물리적·수학적으로 격리하고 제거하는 데 있다. 오차 요인은 발생원에 따라 크게 위성 관련 오차, 대기 전파 경로상의 오차, 수신기 및 관측점 주변 환경에 의한 오차로 분류할 수 있다. 위성 관련 오차의 대표적인 사례는 위성 궤도 오차(Satellite Orbit Error)와 위성 시계 오차(Satellite Clock Error)이다. 위성의 실제 위치와 방송 궤도 정보 사이의 불일치는 위치 결정의 기하학적 정확도를 저해하며, 위성에 탑재된 원자시계의 미세한 편차는 거리 측정값에 직접적인 편향을 유도한다. 이를 보정하기 위해 국제 GNSS 서비스(International GNSS Service, IGS)에서 제공하는 정밀 궤도(Precise Orbit) 및 시계 정보를 활용하며, 이는 방송 위성 항법 메시지보다 수 센티미터 수준의 높은 정밀도를 보장한다.11)

신호가 대기권을 통과하며 발생하는 지연 현상은 정적 측위에서 가장 복잡한 보정 대상 중 하나이다. 전리층 지연(Ionospheric Delay)은 고도 약 50km에서 1,000km 사이의 자유 전자에 의해 발생하며, 신호의 주파수에 따라 지연 정도가 달라지는 분산 매질의 특성을 갖는다. 따라서 L1과 L2 등 서로 다른 주파수를 사용하는 이중 주파수(Dual-frequency) 관측값을 선형 결합한 전리층 제거 조합(Ionosphere-free combination)을 통해 오차의 대부분을 상쇄할 수 있다. 반면, 대류권 지연(Tropospheric Delay)은 비분산 매질에서 발생하므로 주파수에 따른 분리가 불가능하다. 대신 표준 대기 모델이나 사스타모이넨 모델(Saastamoinen model)과 같은 수학적 모델을 적용하거나, 정적 측위의 장시간 관측 데이터를 기반으로 천정 방향 지연(Zenith Total Delay)을 미지수로 설정하여 최소제곱법으로 추정하는 방식을 취한다.

수신기 주변의 지형지물에 의해 신호가 반사되어 들어오는 다중 경로 오차(Multipath Error)는 관측 환경에 종속적인 오차이다. 이는 수학적 모델링이 매우 까다롭기 때문에, 실무에서는 초크 링 안테나(Choke Ring Antenna)와 같은 특수 하드웨어를 사용하거나 신호 차폐가 없는 개활지를 선정함으로써 물리적으로 억제한다. 또한 수신기 시계 오차(Receiver Clock Error)와 안테나 내부의 위상 중심 변동(Phase Center Variation, PCV) 역시 정밀한 좌표 산출을 위해 반드시 고려되어야 할 항목이다. 안테나의 물리적 중심과 전기적 신호가 수렴하는 위상 중심은 일치하지 않으며, 위성의 고도각과 방위각에 따라 변하므로 이에 대한 보정 모델 적용이 필수적이다.

정적 측위에서 이러한 공통 오차를 제거하는 가장 강력한 수학적 도구는 상대 측위(Relative Positioning)에 기반한 차분(Differencing) 기술이다. 두 대의 수신기가 동일한 위성을 동시에 관측할 때 발생하는 관측값의 차이를 구하는 단일 차분(Single Difference)은 위성 시계 오차를 상쇄하며, 두 위성과 두 수신기 사이의 관측값 차이를 다시 결합하는 이중 차분(Double Difference)은 수신기 시계 오차까지 완전히 제거한다. 이 과정에서 반송파 위상(Carrier Phase) 관측값에 포함된 정수 모호정(Integer Ambiguity)을 정확히 결정하는 것이 정적 측위의 최종적인 정밀도를 좌우한다. 이중 차분된 관측식은 다음과 같은 일반적인 형태로 표현된다.

$$ \nabla\Delta \Phi = \nabla\Delta \rho + \lambda \nabla\Delta N + \nabla\Delta \epsilon $$

여기서 $ $는 이중 차분된 위상 관측값, $ $는 위성과 수신기 사이의 기하학적 거리 차이, $ $는 반송파의 파장, $ N $은 이중 차분된 정수 모호정, $ $은 잔여 오차 항을 의미한다. 정적 측위는 충분한 시간 동안 데이터를 누적함으로써 이러한 잔여 오차를 통계적으로 최소화하며, 최종적으로 밀리미터 단위의 기선 성분을 도출한다.

대기 지연 오차

범지구 위성 항법 시스템(GNSS)의 신호가 전리층(Ionosphere)과 대류권(Troposphere)을 통과하며 발생하는 대기 지연은 정적 측위의 정밀도를 결정짓는 가장 핵심적인 오차 요인이다. 위성에서 송신된 전자기파 신호는 진공 상태인 외기권을 지나 지구 대기권에 진입하면서 매질의 밀도와 물리적 특성에 따라 굴절(Refraction)과 속도 지연을 겪는다. 이러한 현상은 신호의 도달 시간을 실제보다 길게 측정하게 함으로써 위성과 수신기 사이의 거리에 오차를 유발하며, 이를 적절히 보정하지 않을 경우 수 미터에서 수십 미터에 이르는 좌표 오차가 발생할 수 있다.

전리층은 지상 약 50km에서 1,000km 사이의 영역으로, 태양의 자외선 및 X선에 의해 대기 분자가 이온화되어 자유 전자가 밀집된 층이다. 전리층에서의 신호 지연은 매질 내의 전전자수(Total Electron Content, TEC)에 비례하며, 신호의 주파수에 따라 굴절률이 변하는 분산(Dispersive) 매질의 특성을 가진다. 전리층 지연량 $d_{ion}$은 1차 근사적으로 다음과 같은 관계식을 따른다.

$$d_{ion} \approx \frac{40.3 \cdot TEC}{f^2}$$

여기서 $f$는 신호의 주파수를 의미한다. 이러한 주파수 의존적 특성 덕분에 L1L2 등 두 개 이상의 서로 다른 주파수를 수신할 수 있는 장비를 사용하면, 전리층 제거 선형 조합(Ionosphere-Free combination)을 통해 전리층 오차의 99% 이상을 제거할 수 있다. 단일 주파수 수신기의 경우에는 GPS 위성이 송출하는 클로부차 모델(Klobuchar model)이나 유럽의 네이퀵(NeQuick) 모델과 같은 경험적 모델을 사용하여 오차를 보정한다12).

대류권은 지표면에서 약 50km 고도까지의 중성 대기층을 의미하며, 전리층과 달리 주파수에 따라 굴절률이 변하지 않는 비분산(Non-dispersive) 매질이다. 따라서 대류권 지연은 다중 주파수 관측을 통해서도 제거되지 않으며, 물리적 모델링에 의존해야 한다. 대류권 지연은 크게 건조 대기의 질량에 의한 정역학적 지연(Hydrostatic delay)과 대기 중 수증기에 의한 습윤 지연(Wet delay)으로 구분된다. 정역학적 지연은 전체 대류권 지연의 약 90%를 차지하며 지표의 기압과 기온 정보를 통해 1mm 이내의 정밀도로 예측이 가능하다. 반면 습윤 지연은 수증기의 불규칙한 시공간적 분포로 인해 예측 모델의 오차가 상대적으로 크다.

대류권 오차 보정에는 사스타모이넨 모델(Saastamoinen model)이나 홉필드 모델(Hopfield model)이 표준적으로 사용된다. 천정 방향의 총 지연량(Zenith Total Delay, ZTD)은 위성의 고도각에 따른 투영 함수(Mapping Function)를 통해 실제 신호 경로상의 지연량으로 변환된다13). 정적 측위에서는 동일한 지점에서 장시간 데이터를 수집하므로, 대류권 지연량을 고정된 상수가 아닌 미지의 매개변수로 설정하여 최소제곱법이나 칼만 필터를 통해 좌표값과 동시에 추정함으로써 밀리미터 단위의 정밀도를 확보한다.

위성 및 수신기 관련 오차

정적 측위의 정밀도를 확보하기 위해서는 위성 시스템과 수신기 하드웨어 자체에서 기인하는 다양한 물리적 오차 요인을 정밀하게 모델링하고 보정해야 한다. 이러한 오차는 크게 위성 측의 요인인 궤도 및 시계 오차와 수신기 측의 요인인 시계 및 안테나 관련 오차로 구분된다.

위성 궤도 오차(Ephemeris Error)는 위성에서 송신하는 항법 메시지에 포함된 위성의 예측 위치와 실제 우주 공간에서의 위성 위치 사이의 차이를 의미한다. 지상 관제국은 위성의 궤도를 추적하여 예측된 궤도 정보를 위성에 업로드하며, 위성은 이를 방송 궤도(Broadcast Ephemeris)의 형태로 사용자에게 전달한다. 그러나 태양 복사압, 지구 중력장의 불균일성, 달과 태양의 인력 등 복합적인 섭동 요인으로 인해 방송 궤도는 수 미터 수준의 오차를 포함하게 된다. 고정밀 좌표 산출이 목적인 정적 측위에서는 국제 GNSS 서비스(International GNSS Service, IGS) 등 전문 기관에서 사후 분석을 통해 제공하는 정밀 궤도 데이터를 활용하여 이 오차를 수 센티미터 이내로 줄인다.

위성 시계 오차(Satellite Clock Error)는 위성에 탑재된 원자시계가 가리키는 시간과 GNSS 표준 시간 사이의 불일치로 인해 발생한다. 위성 시계는 극도로 안정적이지만 미세한 주파수 표류(Drift)가 존재하며, 1나노초($10^{-9}$초)의 시계 오차는 약 30cm의 거리 오차를 유발하므로 정밀한 보정이 필수적이다. 또한, 알베르트 아인슈타인상대성 이론에 따른 보정이 수반되어야 한다. 위성의 빠른 이동 속도로 인한 특수 상대성 이론적 효과와 지구 중력 차이에 의한 일반 상대성 이론적 효과를 종합하면, 지상의 시계보다 위성의 시계가 하루에 약 38마이크로초($\mu s$) 정도 빠르게 흐르게 된다. 시스템 설계 단계에서 위성 주파수를 미세하게 조정하여 이를 상쇄하지만, 타원 궤도에 따른 미세 변화량은 수신기에서 별도의 수식을 통해 실시간으로 보정한다.

수신기 시계 오차(Receiver Clock Error)는 수신기에 내장된 수정 발진기(Quartz Oscillator)의 불안정성에서 기인한다. 수신기 시계는 위성의 원자시계에 비해 정밀도가 현저히 낮아 매우 큰 오차를 발생시키지만, 정적 측위에서는 이를 수학적 기법으로 해결한다. 네 개 이상의 위성을 동시에 관측하여 수신기 시계 오차를 미지수로 포함해 산출하거나, 두 대의 수신기가 동일한 위성을 관측할 때 발생하는 관측값의 차이를 이용하는 이중 차분(Double Differencing) 기법을 통해 수신기 시계 오차 항을 완전히 소거한다.

안테나 위상 중심(Antenna Phase Center, APC) 오차는 신호가 수신되는 전기적 중심점이 안테나의 물리적·기하학적 중심과 일치하지 않아 발생하는 하드웨어적 변위이다. 안테나의 전기적 중심은 입사되는 위성 신호의 고도각방위각에 따라 미세하게 변화하는데, 이를 위상 중심 변동(Phase Center Variation, PCV)이라 하며 기하학적 중심과의 고정된 거리 차이를 위상 중심 오프셋(Phase Center Offset, PCO)이라 한다. 정밀 정적 측위 시에는 사용된 안테나 모델의 고유 교정값을 적용하여 수 밀리미터에서 수 센티미터에 이르는 위상 중심 오차를 보정함으로써 최종적인 좌표의 신뢰도를 높인다.

이외에도 위성과 수신기의 하드웨어 내부 회로에서 신호가 처리될 때 발생하는 기기적 편향(Hardware Bias) 등이 존재한다. 이러한 요인들은 단독 측위보다 두 지점 이상의 관측값을 비교 분석하는 상대 측위 방식을 취할 때 대부분 상쇄되거나 무시할 수 있는 수준으로 제어된다. 따라서 정적 측위는 이러한 하드웨어적 한계를 극복하기 위해 충분한 관측 시간을 확보하고 정밀 보정 모델을 적용하는 과정을 거친다.

주변 환경에 의한 오차

정적 측위의 정밀도를 저해하는 요인 중 관측점 주변의 지형지물에 의해 발생하는 오차는 위성이나 대기 상태와 무관하게 관측 환경의 물리적 특성에 따라 결정된다. 이러한 오차는 크게 다중 경로 오차(Multipath Error)와 신호 가림 및 회절(Diffraction) 현상으로 구분된다. 위성에서 송신된 전자기파 신호는 수신 안테나에 직접 도달하는 직접파(Direct Signal) 외에도, 주변의 건물, 지표면, 수면 등에 반사되어 도달하는 반사파(Reflected Signal)를 포함하게 된다. 수신기는 이 두 신호의 합성파를 처리하게 되는데, 이때 반사파의 이동 거리가 직접파보다 길기 때문에 위상 지연이 발생하며 이는 최종적인 거리 측정치에 오차를 유발한다.

다중 경로 오차는 반송파 위상(Carrier Phase) 관측량과 코드 관측량 모두에 영향을 미치나, 그 영향의 크기와 특성은 상이하다. 코드 관측에서는 오차의 크기가 수십 미터에 달할 수 있는 반면, 반송파 위상 관측에서는 파장의 4분의 1 수준인 수 센티미터 이내로 제한되는 경향이 있다. 그러나 밀리미터 단위의 정밀도를 추구하는 정적 측위에서 이러한 미세한 오차는 정수 모호정(Integer Ambiguity) 결정을 어렵게 하거나 좌표 산출의 신뢰도를 떨어뜨리는 주요 원인이 된다. 수신된 신호 $ S(t) $는 직접파와 $ n $개의 반사파의 합으로 다음과 같이 표현될 수 있다.

$$ S(t) = A_0 \cos(\phi_0) + \sum_{i=1}^{n} \alpha_i A_0 \cos(\phi_0 + \Delta\phi_i) $$

여기서 $ A_0 $와 $ _0 $는 직접파의 진폭과 위상이며, $ _i $와 $ _i $는 각 반사파의 상대적 감쇄 계수와 위상 변화량을 의미한다. 다중 경로 오차는 시간에 따른 위성 배치의 변화에 따라 주기적인 특성을 보이므로, 정적 측위에서는 관측 시간을 충분히 길게 확보하여 이를 통계적으로 평균화함으로써 오차를 완화한다.14)

관측 지점 주변의 고층 건물, 수목, 지형 등에 의한 신호 가림 현상은 가시 위성의 수를 감소시켜 기하학적 정밀도 저하율(Dilution of Precision, DOP)을 악화시킨다. 특히 신호가 완전히 차단되지 않더라도 구조물의 모서리 등을 통과하며 발생하는 회절 현상은 신호 대 잡음비(Signal-to-Noise Ratio, SNR)를 저하시키고 위상 측정의 불확실성을 증대시킨다.15) 회절된 신호는 직접파에 비해 강도가 약하고 위상 편이가 발생하여 수신기 내부의 상관기(Correlator)가 신호를 추적하는 과정에서 오차를 유발한다. 이는 도심지나 산악 지형에서의 정적 측위 시 좌표의 수직 성분 정밀도를 저하시키는 결정적인 요인이 된다.

이러한 주변 환경 오차를 최소화하기 위해 하드웨어 측면에서는 초크 링 안테나(Choke Ring Antenna)를 사용한다. 초크 링 안테나는 동심원 형태의 금속 고랑을 통해 지표면이나 주변 구조물에서 반사되어 하부에서 올라오는 신호를 물리적으로 차단하고 감쇄시키는 역할을 한다. 또한 수신기 내부적으로는 다중 경로에 내성이 있는 신호 처리 알고리즘을 적용하거나, 관측 데이터 후처리 과정에서 프레넬 존(Fresnel Zone) 분석을 통해 신호 품질이 낮은 위성을 배제하기도 한다. 결과적으로 정적 측위의 성패는 이러한 국지적 환경 요인을 사전에 분석하여 최적의 관측점을 선정하는 선점 단계의 효율성에 크게 의존한다.

정적 측위의 실행 단계와 방법론

정적 측위의 성공적인 수행은 체계적인 관측 계획 수립에서 시작된다. 가장 먼저 고려해야 할 요소는 관측점의 선점(Reconnaissance)과 위성 배치 상황의 분석이다. 관측점은 위성 신호의 수신을 방해하는 상부 차폐물이 없고, 수신된 신호가 주변 지형지물에 반사되어 발생하는 다중 경로 오차(Multipath Error)를 최소화할 수 있는 개활지를 선정해야 한다. 특히 정밀한 좌표 산출을 위해서는 정밀도 저하율(Dilution of Precision, DOP) 수치가 낮은 시간대를 파악하는 것이 필수적이다. 이를 위해 범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS) 운용 기관에서 제공하는 위성 궤도력 정보를 활용하여, 관측 예정 시간에 가시 위성 수가 충분한지 그리고 위성들이 천구(Celestial Sphere)에 고르게 분포되어 있는지를 사전에 시뮬레이션한다.

현장 관측 단계에서는 장비의 정밀한 거치와 기록의 정확성이 요구된다. 수신기를 삼각대 위에 설치할 때에는 정준구심을 엄격히 시행하여 안테나 위상 중심(Antenna Phase Center)이 지상 기준점의 수직 상단에 위치하도록 한다. 이때 안테나 높이의 측정 오차는 최종 수직 좌표의 오차로 직결되므로, 줄자 등을 이용하여 다각도에서 높이를 측정하고 이를 관측 기록지에 상세히 기입한다. 관측이 시작되면 정해진 데이터 수신 간격(Epoch)에 따라 반송파 위상(Carrier Phase) 관측값을 수집한다. 두 대 이상의 수신기를 동시에 운용하는 상대 측위 방식에서는 수신기 간의 동기화가 필수적이며, 기선의 길이에 따라 수십 분에서 수 시간 이상의 충분한 관측 시간을 확보하여 전리층대류권 지연 오차를 통계적으로 상쇄할 수 있도록 한다.

수집된 원시 데이터는 실내에서 전문 소프트웨어를 이용한 데이터 후처리 과정을 거친다. 먼저 각 제조사의 고유 형식으로 저장된 데이터를 공용 형식인 수신기 독립 교환 형식(Receiver Independent Exchange Format, RINEX)으로 변환한다. 이후 수행되는 기선 해석(Baseline Analysis)은 두 관측점 간의 상대적인 위치 벡터를 산출하는 핵심 과정이다. 이 과정에서 가장 중요한 수학적 과제는 정수 모호정(Integer Ambiguity)을 결정하는 것이다. 반송파의 파장 수를 정확히 결정함으로써 센티미터 수준의 정밀도를 확보하며, 해석 결과의 신뢰도를 평가하기 위해 분산 분석과 잔차 검토를 병행한다. 해석이 완료된 기선 벡터들은 하나의 유기적인 망을 형성하게 된다.

마지막 단계인 망 조정(Network Adjustment)은 개별적으로 해석된 기선들을 결합하여 전체 측량망의 최적해를 구하는 과정이다. 각 기선 벡터에 포함된 미세한 관측 오차를 수학적으로 배분하기 위해 최소제곱법(Least Squares Method)이 주로 사용된다. 고정된 기지점의 좌표를 기준으로 자유망 조정구속망 조정을 순차적으로 실시하여 망의 내적 정밀도와 외적 정확도를 검증한다. 조정 계산 결과 출력되는 좌표값과 오차 타원(Error Ellipse) 분석을 통해 허용 오차 범위 내에 들어오는지를 최종 확인하며, 이 모든 과정을 거쳐 도출된 측량 성과는 국가 기준점 체계에 부합하는 정밀 좌표로서의 효력을 갖게 된다.16)

관측 계획 및 선점

정적 측위의 성패와 최종적인 좌표의 정밀도는 단순히 수신기를 가동하는 시간에 비례하는 것이 아니라, 관측 전 단계에서 수행되는 체계적인 관측 계획 수립과 철저한 현장 선점에 의해 결정된다. 이는 범지구 위성 항법 시스템(GNSS)의 신호 특성과 수신 환경의 물리적 제약 조건을 사전에 분석하여 오차 발생 가능성을 최소화하기 위한 필수적인 공정이다.

관측 계획의 일차적인 목표는 위성 궤도 정보인 알마낙(Almanac) 데이터를 활용하여 최적의 관측 시간대를 선정하는 것이다. GNSS 위성은 지구 주위를 공전하므로 특정 관측점 상공에 나타나는 위성의 수와 기하학적 배치는 시간에 따라 끊임없이 변화한다. 이러한 위성의 배치 상태가 측위 정밀도에 미치는 영향을 수치화한 지표가 기하학적 정밀도 저하율(Dilution of Precision, DOP)이다. 일반적으로 위치 정밀도 저하율(Position Dilution of Precision, PDOP) 값이 낮을수록 위성들이 상공에 고르게 분산되어 있음을 의미하며, 이는 기하학적 강도를 높여 관측 오차의 전파를 억제하는 효과를 가진다. 따라서 정밀도가 요구되는 정적 측위에서는 가용 위성 수가 최대화되고 PDOP 값이 최소화되는 시간대를 사전에 예측하여 관측 일정을 수립해야 한다.

현장 선점(Reconnaissance)은 이론적으로 검토된 관측점 후보지를 직접 방문하여 실제 수신 환경을 조사하고 최적의 안테나 설치 위치를 확정하는 과정이다. 선점 시 가장 우선적으로 고려해야 할 요소는 상공 시계의 확보이다. GNSS 신호는 직진성이 강하여 수목, 건물, 산악 지형 등의 장애물에 의해 쉽게 차단되거나 회절된다. 특히 정수 모호정(Integer Ambiguity)의 안정적인 결정을 위해서는 안테나 주위의 고도각 15도 이상 범위에 차폐물이 없는 지점을 선정하는 것이 원칙이다. 만약 상공 일부분이 가려질 경우, 해당 방향의 위성 신호가 단절되면서 싸이클 슬립(Cycle Slip)이 발생하여 데이터 처리에 난항을 겪을 수 있다.

또한, 선점 과정에서는 다중 경로 오차(Multipath Error)를 유발할 수 있는 요인을 엄격히 배제해야 한다. 이는 위성으로부터 직접 도달하는 신호 외에 주변의 대형 건축물, 금속성 구조물, 혹은 수면 등에 반사된 신호가 수신기에 혼입되어 거리 측정에 오차를 일으키는 현상이다. 이를 방지하기 위해 반사체로부터 충분한 거리를 이격하고, 강력한 전자기파를 방출하여 신호 간섭을 일으킬 수 있는 통신탑이나 고압선 인근은 피해야 한다. 지반의 안정성 역시 중요한 고려 사항이다. 정적 측위는 장시간 동일한 지점에서 데이터를 수집하므로, 관측 기간 동안 안테나의 미세한 침하나 진동이 발생하지 않는 견고한 지반을 선택해야 한다.

복수의 수신기를 동시에 운용하는 상대 측위(Relative Positioning) 방식에서는 각 관측점 간의 동시 관측 시간 확보가 핵심이다. 기선(Baseline)의 길이에 따라 요구되는 최소 관측 시간이 달라지는데, 수십 킬로미터 이상의 장기선일수록 대기 지연 등 공통 오차의 상관성이 낮아지므로 이를 통계적으로 상쇄하기 위해 더 긴 관측 세션이 요구된다. 선점이 완료되면 관측점의 위치와 주변 상황, 장애물 분포 등을 상세히 기록한 선점도를 작성한다. 이는 실제 관측 시 안테나 설치의 일관성을 유지하고, 사후 데이터 처리 단계에서 특정 위성 신호의 품질 저하 원인을 분석하는 중요한 근거 자료로 활용된다. 이러한 일련의 사전 준비 과정은 측량의 효율성을 극대화하고 최종 산출물의 신뢰도를 학술적 수준으로 끌어올리는 토대가 된다.

현장 관측 및 데이터 수집

현장 관측 및 데이터 수집은 정적 측위의 정밀도를 결정짓는 실질적인 데이터 획득 단계로, 사전에 수립된 관측 계획에 따라 현장에서 범지구 위성 항법 시스템(GNSS) 수신기를 운용하는 과정을 포함한다. 이 단계에서의 미세한 실수는 후처리 과정에서 보정하기 어려운 치명적인 오차를 유발할 수 있으므로, 엄격한 표준 지침의 준수가 요구된다. 현장 작업은 크게 수신기 설치, 안테나 높이 측정, 데이터 로그 기록 및 관측 야장 작성으로 구분된다.

수신기 설치의 핵심은 정확한 구심(Centering)과 정준(Leveling)에 있다. 관측점의 중심인 표석 또는 대표정(Mark) 위에 삼각대(Tripod)를 거치하고, 정밀 구심 장치를 사용하여 수신 안테나의 연직축이 점의 중심과 일치하도록 조정한다. 이때 수평 기포관을 활용하여 안테나의 수평 상태를 유지하는 정준 과정을 수행하는데, 이는 안테나의 지향성을 확보하고 안테나 위상 중심(Antenna Phase Center, APC) 오차를 최소화하기 위함이다. 또한, 주변의 대형 건축물, 금속 울타리, 수면 등 위성 신호의 반사를 유도할 수 있는 지형지물로부터 충분한 거리를 확보하여 다중 경로 오차(Multipath Error)를 방지해야 한다.

안테나 높이 측정은 수신된 좌표를 지면의 기준점 좌표로 변환하기 위한 필수적인 절차이다. 일반적으로 수신 안테나의 하단 기준면인 안테나 기준점(Antenna Reference Point, ARP)까지의 수직 높이를 측정한다. 현장 상황에 따라 안테나의 측면까지의 사거리(Slant Height)를 측정한 후, 안테나의 반경 정보를 이용한 기하학적 계산을 통해 수직 높이로 변환하기도 한다. 안테나 높이의 측정 오차는 그대로 수직 위치 오차로 전이되므로, 관측 시작 전과 종료 후에 각각 측정하여 그 값의 일치 여부를 확인하는 것이 원칙이다.

데이터 로그 기록은 수신기가 포착한 위성 신호를 내부 메모리나 외부 저장 장치에 저장하는 과정이다. 정적 측위에서는 보통 15초 또는 30초 단위의 에포크(Epoch) 간격을 설정하여 데이터를 수집하며, 정밀도 요구 수준과 기선의 길이에 따라 수 시간 이상의 관측 세션을 유지한다. 수집되는 데이터는 위성별 의사거리(Pseudorange)와 반송파 위상 측정값, 도플러 효과(Doppler Effect) 측정값 등을 포함하며, 후처리의 호환성을 위해 수신기 독립 교환 형식(Receiver Independent Exchange Format, RINEX)으로 변환될 수 있는 원시 데이터 형태로 보존된다.

마지막으로 현장 작업자는 관측 야장(Field Note)을 충실히 작성하여 데이터의 신뢰성을 보증해야 한다. 야장에는 관측점 번호, 수신기 및 안테나의 일련번호, 안테나 높이, 관측 시작 및 종료 시각뿐만 아니라, 관측 당시의 기상 상태와 주변 환경의 특이 사항을 상세히 기록한다. 이러한 메타데이터는 이후 데이터 후처리 과정에서 특정 시간대의 이상치를 식별하거나 오차 모델을 적용할 때 중요한 판단 근거가 된다. 정적 측위의 모든 현장 절차는 표준화된 작업 규정에 따라 수행되어야 하며, 이는 측량 결과의 법적·기술적 공신력을 확보하는 토대가 된다.

데이터 후처리 및 망 조정

현장에서 수집된 범지구 위성 항법 시스템(GNSS)의 원시 데이터는 그 자체로 최종적인 좌표 성과가 될 수 없으며, 반드시 실내에서의 데이터 후처리(Data Post-processing) 과정을 거쳐야 한다. 정적 측위의 후처리는 크게 개별 기선의 상대적 위치 관계를 결정하는 기선 해석(Baseline Processing)과, 해석된 기선들을 결합하여 전체 측량망의 최적해를 산출하는 망 조정(Network Adjustment) 단계로 구분된다. 이 과정은 관측 시 발생한 확률적 오차를 수학적으로 정제하고, 기하학적 모순을 해결하여 측량 성과의 신뢰성을 확보하는 데 목적이 있다.

기선 해석은 두 대 이상의 수신기에서 동시에 관측된 반송파 위상(Carrier Phase) 데이터를 바탕으로 수신기 간의 좌표 차이인 기선 벡터를 산출하는 작업이다. 분석 시에는 현장에서 기록된 수신기 독립 교환 형식(Receiver Independent Exchange Format, RINEX) 파일을 활용하며, 위성의 위치 정보인 궤도력(Ephemeris)을 결합한다. 이때 실시간으로 송출되는 방송 궤도(Broadcast Ephemeris) 대신 국제 GNSS 서비스(International GNSS Service, IGS) 등에서 제공하는 정밀 궤도(Precise Ephemeris)를 사용하면 위성 위치 오차를 획기적으로 줄일 수 있다17). 기선 해석의 핵심은 수신기와 위성 사이의 정수 배 파장 수를 의미하는 정수 모호정(Integer Ambiguity)을 정확히 결정하는 데 있으며, 이를 위해 다양한 통계적 기법과 필터링 알고리즘이 적용된다.

개별 기선들에 대한 해석이 완료되면, 이들을 하나의 유기적인 체계로 묶는 망 조정 단계에 진입한다. 현장 관측값에는 불가피하게 우연 오차(Random Error)가 포함되어 있어, 여러 기선이 폐합 회로를 형성할 때 수치적 불일치가 발생하게 된다. 이러한 모순을 해결하기 위해 오차론의 핵심인 최소제곱법(Least Squares Method)이 동원된다. 최소제곱법은 관측값의 잔차(Residual) 제곱합에 가중치를 곱한 값이 최소가 되도록 미지수를 결정하는 원리이다. 기본적인 선형 모델은 다음과 같이 표현된다.

$$ V = AX - L $$

여기서 $ V $는 잔차 벡터, $ A $는 관측값과 미지수의 기하학적 관계를 나타내는 설계 행렬, $ X $는 보정해야 할 미지수 벡터, $ L $은 관측값과 근사값의 차이를 나타내는 벡터이다. 이에 대한 최적해 $ $는 가중 행렬 $ P $를 고려하여 다음과 같은 정규 방정식을 통해 산출된다.

$$ \hat{X} = (A^T PA)^{-1} A^T PL $$

망 조정은 통상적으로 두 단계로 나누어 수행한다. 먼저 기지점의 좌표를 고정하지 않고 관측값 사이의 내부적 일치성만을 검토하는 자유망 조정(Free Network Adjustment)을 실시하여, 특정 기선에 포함된 착오(Blunder)나 이상치를 식별한다18). 이후 신뢰성이 검증된 기선들을 바탕으로 이미 좌표가 알려진 국가 기준점에 결합하는 구속망 조정(Constrained Network Adjustment)을 수행하여 최종 성과를 도출한다.

최종 산출된 좌표의 품질을 평가하기 위해 카이제곱 검정(Chi-square Test)과 같은 통계적 검정이 수반된다. 이는 조정된 결과가 통계적으로 유의미한 범위 내에 있는지, 그리고 설정된 가중치가 적절했는지를 판별하는 지표가 된다. 만약 검정 결과가 허용 범위를 초과할 경우, 데이터 품질이 불량한 기선을 제거하거나 재관측을 수행하여 측량망의 정밀도를 보정해야 한다.

정적 측위의 주요 응용 분야

정적 측위(Static Positioning)는 범지구 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 활용하여 얻을 수 있는 가장 높은 수준의 정밀도를 제공하기 때문에, 밀리미터(mm) 단위의 오차 허용 범위가 요구되는 국가적 기반 시설 구축 및 학술적 연구 분야에서 중추적인 역할을 수행한다. 이 방식은 수신기를 장시간 고정하여 관측함으로써 대기 지연이나 위성 궤도 오차 등의 가변적 요인을 통계적으로 상쇄할 수 있다는 이점을 가진다.

가장 대표적인 응용 분야는 국가 위치 결정의 근간이 되는 국가기준점(National Control Point) 및 측량 기준망의 구축이다. 국토지리정보원과 같은 국가 기관은 전 국토의 수평 위치와 높이 기준을 설정하기 위해 통합기준점, 삼각점, 수준점 등을 설치하고 관리한다. 이러한 기준점들은 지적재조사나 각종 토목 설계의 기초가 되며, 정적 측위를 통해 산출된 고정밀 좌표는 국가 좌표계의 신뢰성을 담보하는 핵심 요소가 된다19).

지구물리학적 관점에서의 지각 변동(Crustal Deformation) 및 지반 침하(Land Subsidence) 감시 역시 정적 측위의 핵심적인 활용 영역이다. 판 구조론에 따른 지각판의 이동 속도는 연간 수 센티미터 수준에 불과하므로, 이를 정밀하게 추적하기 위해서는 초정밀 정적 측위가 필수적이다. 특히 지진 발생 전후의 미세한 지각 변위나 화산 활동에 따른 지표면의 팽창 및 수축을 관측함으로써 재난 예측 및 연구에 기여한다20). 또한, 대규모 지하수 개발이나 도심지 굴착 공사로 인해 발생하는 지반의 수직적 변위를 장기적으로 추적하여 시설물의 붕괴 위험을 사전에 인지하는 데 사용된다.

토목 및 건축 분야에서는 대형 구조물의 정밀 안전 진단(Structural Health Monitoring, SHM)에 정적 측위가 도입된다. 장대교량(Long-span Bridge), (Dam), 초고층 빌딩과 같은 거대 구조물은 자중, 풍하중, 온도 변화 등에 의해 미세한 변형이 발생한다. 이러한 구조물의 특정 지점에 GNSS 안테나를 고정하고 주기적으로 정적 측위를 수행함으로써, 구조물의 기하학적 형태 변화를 정밀하게 모니터링할 수 있다. 이는 육안 점검으로 확인하기 어려운 구조적 결함을 수치적으로 파악하여 대형 사고를 방지하는 근거를 제공한다21).

정적 측위의 정밀도($ $)는 일반적으로 기선 거리($ L $)에 비례하며, 다음과 같은 수식으로 표현되는 경우가 많다. $$ \sigma = \pm(a + b \times L) $$ 여기서 $ a $는 수신기 자체의 고정 오차를, $ b $는 기선 거리에 따른 오차 계수를 의미한다. 이러한 높은 정밀도 덕분에 정적 측위는 해양 측량, 경계 측량, 그리고 정밀한 지도 제작을 위한 항공 사진 측량의 지상 기준점(Ground Control Point, GCP) 설치 등 광범위한 산업 분야에서 표준적인 측위 방법론으로 자리 잡고 있다.

국가 기준점 및 측량 기준망 구축

국토의 위치 기준이 되는 삼각점 및 통합 기준점 설치에서의 역할을 다룬다.

지각 변동 및 지반 침하 감시

정적 측위(Static Positioning)는 반송파 위상(Carrier Phase) 관측치를 활용하여 밀리미터(mm) 단위의 고정밀 좌표 산출이 가능하므로, 지구 물리적 현상에 따른 미세한 지각의 움직임이나 국지적인 지반의 수직 변위를 추적하는 데 핵심적인 역할을 수행한다. 특히 판 구조론(Plate Tectonics)에 근거한 지각의 수평 이동과 지질학적 요인 혹은 인위적 요인에 의한 지반 침하(Land Subsidence)를 장기적으로 감시하기 위해서는 연속적이고 정밀한 정적 측위 데이터의 확보가 필수적이다. 이러한 정밀 감시는 지진이나 화산 활동과 같은 지질 재해의 전조 현상을 파악하는 기초 연구로 활용된다.

지각 변동 감시는 전 지구적 차원의 기준 좌표계(Reference Frame) 유지와 밀접하게 관련된다. 지각은 고정된 상태가 아니라 인접한 (Plate)과의 상호작용에 의해 끊임없이 이동하며, 그 속도는 연간 수 센티미터(cm) 수준에 달한다. 이러한 미세한 변화를 정밀하게 측정하기 위해 전 세계적으로 범지구 위성 항법 시스템(GNSS) 상시관측소가 운영되고 있다. 수집된 정적 측위 데이터는 국제 지구 회전 및 기준 좌표계 서비스(International Earth Rotation and Reference Systems Service, IERS)에 의해 처리되어 국제 지구 기준 좌표계(International Terrestrial Reference Frame, ITRF)와 같은 지구 중심 좌표계의 갱신에 활용된다. 이를 통해 특정 지역의 지각 속도(Crustal Velocity) 벡터를 산출할 수 있으며, 이는 지진 발생 가능성 평가나 지구 내부의 역학적 구조 연구의 기초 자료가 된다. 22)

지반 침하 감시 분야에서 정적 측위는 특정 시설물이나 지역의 안정성을 평가하는 정밀 도구로 사용된다. 도심지의 지하 공간 개발, 과도한 지하수 추출, 혹은 연약 지반압밀 현상 등으로 발생하는 지표면의 하강은 구조물의 붕괴나 범람 위험을 초래한다. 정적 측위를 활용한 감시 체계는 대상 지역에 고정된 관측점을 설치하고 주기적인 관측을 수행함으로써 수직 변위의 경향성을 파악한다. 관측된 위치 정보는 시간의 흐름에 따른 시계열 분석(Time Series Analysis)을 통해 처리되며, 다음과 같은 선형 모델을 기본으로 한다.

$$X(t) = X_0 + V(t - t_0) + \sum_{i=1}^{n} [a_i \sin(\omega_i t) + b_i \cos(\omega_i t)] + \epsilon$$

위 식에서 $X(t)$는 시간 $t$에서의 좌표, $X_0$는 기준 시점 $t_0$에서의 초기 좌표이며, $V$는 일정한 속도로 진행되는 지각 변동 혹은 침하 속도를 의미한다. 우항의 삼각함수 항은 기온 변화나 수문학적 하중(Hydrological Loading) 등에 의해 발생하는 계절적 변동 성분을 모델링한 것이며, $\epsilon$은 측정 오차를 나타낸다. 정적 측위는 관측 시간을 길게 확보할수록 이러한 확률적 오차와 계절적 잡음을 효과적으로 분리해낼 수 있어, 순수한 지반 변위량만을 정밀하게 추출하는 데 유리하다. 특히 최소제곱법(Least Squares Method)을 이용한 매개변수 추정 과정에서 관측 시간의 증가는 해의 신뢰도를 높이는 결정적 요인이 된다.

최근에는 정밀 단독 측위(Precise Point Positioning, PPP) 기술의 발전으로 인해 기준국과의 거리에 제약을 받지 않고도 단일 수신기만으로 고정밀 지각 변동 감시가 가능해졌다. 또한, 영상 레이더 간섭계(Interferometric Synthetic Aperture Radar, InSAR) 기술과 정적 측위 데이터를 결합하여 광역적인 지반 변위 지도와 특정 지점의 정밀 변위 시계열을 동시에 확보하는 융합 연구가 활발히 진행되고 있다. 이는 원격 탐사의 광역성과 정적 측위의 고정밀성을 결합한 것으로, 도시 전체의 안정성을 진단하는 데 효과적이다. 이러한 체계적인 감시는 국가 기간 시설의 안전 진단 및 재난 예방을 위한 공간 정보 인프라의 핵심 요소로 자리 잡고 있다. 23)

대형 구조물의 정밀 안전 진단

대형 구조물의 정밀 안전 진단에서 정적 측위는 구조물의 기하학적 형상 변화를 밀리미터(mm) 단위로 추적하여 잠재적인 붕괴 위험을 사전에 탐지하는 핵심 기술로 활용된다. 교량, 댐, 초고층 빌딩과 같은 대규모 사회기반시설은 외부 하중, 온도 변화, 지반 거동 등에 의해 미세한 변형이 상시 발생하며, 이러한 변형이 허용 임계치를 초과할 경우 구조적 결함으로 이어진다. 정적 측위는 장시간의 관측 데이터를 축적함으로써 일시적인 진동을 배제하고 구조물의 영구 변위(Permanent Displacement)나 장기적인 추세 변화를 정밀하게 산출하는 데 최적화되어 있다.

현수교사장교와 같은 장대교량의 경우, 주탑(Pylon)의 기울기나 상판(Deck)의 수직 처짐을 감시하기 위해 정적 측위 기반의 구조물 건강 모니터링(Structural Health Monitoring, SHM) 시스템이 구축된다. 교량의 특정 지점에 GNSS 수신기를 고정하고, 구조물 외부의 안정된 지점에 설치된 기준점(Reference Station)과의 상대 측위를 수행함으로써 기선(Baseline)의 변화를 측정한다. 이를 통해 계절적 온도 변화에 따른 신축 거동이나 풍하중에 의한 주탑의 미세 변위를 분석할 수 있으며, 이는 유한요소법(Finite Element Method, FEM)을 이용한 구조 해석 모델의 타당성을 검증하는 실측 자료가 된다.

댐 구조물에서의 정적 측위는 수압(Hydrostatic Pressure)에 의한 제체의 변형과 기초 지반의 침하를 감시하는 데 중추적인 역할을 한다. 댐은 저수량 변화에 따라 막대한 하중을 받게 되며, 이로 인해 발생하는 미세한 수평 및 수직 변위는 구조물의 안전성과 직결된다. 정적 측위는 다중 위성군(Multi-GNSS)을 활용하여 관측 가용성을 높이고, 이중 차분(Double Difference) 기법을 통해 대기 지연 및 위성 궤도 오차를 극소화함으로써 고정밀 좌표를 산출한다. 이때 산출되는 변위 벡터 $\Delta \mathbf{x}$는 다음과 같이 정의된다.

$$\Delta \mathbf{x} = \mathbf{x}_{current} - \mathbf{x}_{baseline}$$

여기서 $\mathbf{x}_{current}$는 현재 시점의 정밀 좌표이며, $\mathbf{x}_{baseline}$은 구조물 준공 직후나 안정화 단계에서 설정된 기준 좌표이다. 이 벡터의 크기와 방향을 시계열적으로 분석함으로써 댐의 비정상적인 거동 여부를 판별한다.

초고층 빌딩의 경우에는 건물의 수직도 유지와 기초 부등 침하 방지가 안전 진단의 핵심이다. 고층 구조물은 상층부로 갈수록 위성 신호 수신 환경은 양호해지나, 바람에 의한 미세 진동이 상시 존재하므로 순수한 정적 변위를 분리해내는 필터링 기술이 병행되어야 한다. 정적 측위 데이터는 가속도계(Accelerometer)나 경사계(Inclinometer)에서 얻은 동적 데이터와 결합하여 구조물의 통합적인 거동 특성을 규명하는 데 기여한다. 특히 도심지 내의 지반 침하나 인근 공사로 인한 구조물 영향 평가 시, 정적 측위는 장기적인 좌표 안정성을 제공하므로 신뢰도 높은 안전 진단 지표로 기능한다24)25).

최근에는 실시간 이동 측위(Real-Time Kinematic, RTK) 기술의 발전에도 불구하고, 신뢰성이 최우선인 안전 진단 분야에서는 여전히 후처리(Post-processing) 방식의 정적 측위가 선호된다. 이는 충분한 관측 시간을 확보함으로써 정수 모호정(Integer Ambiguity) 결정의 정확도를 높이고, 정밀 궤도(Precise Ephemeris) 정보를 적용하여 오차 요인을 완벽하게 제거할 수 있기 때문이다. 이러한 정밀 안전 진단 체계는 구조물의 수명을 연장하고 예방적 유지관리를 가능케 하여 대형 참사를 방지하는 사회적 안전망의 역할을 수행한다.

1)
NRTK, PPP or Static, That Is the Question. Testing Different Positioning Solutions for GNSS Survey, https://www.mdpi.com/2072-4292/13/7/1406
2) , 3)
Hospodka, J., “Doppler shift satellite navigation - NAVSAT-TRANSIT and adherents”, https://doi.org/10.14311/MAD.2013.02.03
5)
NOAA National Geodetic Survey, “The National Geodetic Survey and the Evolution of the National Spatial Reference System”, https://geodesy.noaa.gov/PUBS_LIB/Evolution_of_NSRS.pdf
6) , 7)
Global Positioning System Carrier Phase: Description and Use, https://repository.library.noaa.gov/view/noaa/51023/noaa_51023_DS1.pdf
8)
Global Positioning Systems Directorate, “NAVSTAR GPS Space Segment/User Segment L1C Interface (IS-GPS-800E)”, https://www.gps.gov/technical/icwg/IS-GPS-800E.pdf
9)
Teunissen, P. J. G., “The least-squares ambiguity decorrelation adjustment: a method for fast GPS integer ambiguity estimation”, https://link.springer.com/article/10.1007/BF00863419
10)
Verhagen, S., “The GNSS integer ambiguities: estimation and validation”, https://repository.tudelft.nl/islandora/object/uuid:998b965b-43d9-4820-8041-38507c30999a
12)
Regiomontan: A Regional High Precision Ionosphere Delay Model and Its Application in Precise Point Positioning, https://pmc.ncbi.nlm.nih.gov/articles/PMC7284771/
13)
Modeling wide‑area tropospheric delay corrections for fast PPP ambiguity resolution, https://gfzpublic.gfz-potsdam.de/pubman/item/item_5010427_1/component/file_5010766/5010427.pdf
15)
Performance assessment of GNSS diffraction models in urban areas, https://navi.ion.org/content/68/2/369
16)
Methods of practice and guidelines for using survey-grade global navigation satellite systems (GNSS) to establish vertical datum in the United States Geological Survey, https://www.usgs.gov/publications/methods-practice-and-guidelines-using-survey-grade-global-navigation-satellite-systems
17)
차상헌 외, 2등 측지기준점 GPS 관측데이터의 기선벡터 추정, https://www.koreascience.or.kr/article/JAKO200821036731209.page?lang=ko
18)
남기범 외, 국가기준점 망조정을 위한 GPS 3등기준점 기선해석, https://www.koreascience.kr/article/CFKO200716419439669.page?lang=ko
19)
국가기준점 성과 재고시가 지적재조사 측량성과에 미치는 영향분석 연구, https://www.dbpia.co.kr/journal/articleDetail?nodeId=NODE11190007
20)
정밀절대측위를 이용한 양산지역의 지각변위 해석, https://www.dbpia.co.kr/journal/articleDetail?nodeId=NODE02210404
21)
GPS 측량에 의한 해안구조물의 정적변위측정에 관한 연구, https://www.koreascience.or.kr/article/JAKO200011921803103.page?lang=ko
22) , 23)
준실시간 지각변동 모니터링 체계구축을 위한 초고정밀 GPS 관측데이터 연속처리 기술 연구, https://data.doi.or.kr/10.23000/TRKO201800005180?lang=ko
24)
GNSS integrated displacement and attitude determination for structural health monitoring of long-span bridges, https://satellite-navigation.springeropen.com/articles/10.1186/s43020-025-00174-9
25)
Analysis of structural monitoring with multi-GNSS positioning: comparison between PPP and static relative strategies, https://link.springer.com/article/10.1007/s12518-025-00649-z
정적_측위.1776237223.txt.gz · 마지막으로 수정됨: 저자 flyingtext