| 양쪽 이전 판이전 판다음 판 | 이전 판 |
| 동경측지계 [2026/04/14 22:21] – 동경측지계 sync flyingtext | 동경측지계 [2026/04/14 22:26] (현재) – 동경측지계 sync flyingtext |
|---|
| ==== 베셀 타원체의 채택과 기하학적 정의 ==== | ==== 베셀 타원체의 채택과 기하학적 정의 ==== |
| |
| 동경측지계의 기술적 근간은 [[독일]]의 수학자이자 천문학자인 [[프리드리히 빌헬름 베셀]](Friedrich Wilhelm Bessel)이 1841년에 발표한 [[베셀 타원체]](Bessel 1841 Ellipsoid)를 [[참조 타원체]](Reference Ellipsoid)로 채택한 데서 시작된다. 19세기 당시 [[측지학]](Geodesy) 분야에서 지구의 형상을 정의하기 위한 시도는 유럽 각국에서 활발히 진행되었으며, 베셀은 유럽과 러시아 등지에서 수행된 10개의 [[호측량]](Arc measurement) 결과에 [[최소제곱법]](Least Squares Method)을 적용하여 지구의 기하학적 형태를 가장 잘 나타내는 회전 타원체의 제원을 도출하였다. 일본이 메이지 시대에 근대적 측량 체계를 수립하며 베셀 타원체를 선정한 이유는 당시 이 모델이 동아시아 지역의 [[지오이드]](Geoid) 기복과 지형적 특성을 가장 정밀하게 반영하는 것으로 평가받았기 때문이다. | [[동경측지계]]의 기술적 근간은 [[독일]]의 수학자이자 천문학자인 [[프리드리히 빌헬름 베셀]](Friedrich Wilhelm Bessel)이 1841년에 발표한 [[베셀 타원체]](Bessel 1841 ellipsoid)를 [[참조 타원체]](reference ellipsoid)로 채택한 데서 시작된다. 19세기 당시 [[측지학]](geodesy) 분야에서 지구의 형상을 정의하기 위한 시도는 유럽 각국에서 활발히 진행되었으며, 베셀은 유럽, 러시아, 인도 등지에서 수행된 10개의 [[호측량]](arc measurement) 성과를 종합하였다. 그는 관측 데이터에 포함된 오차를 최소화하기 위해 [[최소제곱법]](least squares method)을 적용하였고, 이를 통해 지구의 기하학적 형태를 가장 잘 나타내는 회전 타원체의 제원을 도출하였다. [[일본 제국]]이 [[메이지 유신]] 이후 근대적 측량 체계를 수립하며 베셀 타원체를 선정한 이유는 당시 이 모델이 동아시아 지역의 [[지오이드]](geoid) 기복과 지형적 특성을 가장 정밀하게 반영하는 표준 모델로 간주되었기 때문이다. |
| |
| 기하학적 관점에서 베셀 타원체는 지구의 자전축을 회전축으로 하는 [[회전 타원체]](Oblate Spheroid)로 정의된다. 타원체의 형상을 결정하는 핵심 요소는 적도 반지름에 해당하는 [[장반경]](Semi-major axis, $ a $)과 극 반지름에 해당하는 [[단반경]](Semi-minor axis, $ b $)이다. 베셀이 제시한 수치적 정의에 따르면 장반경과 [[편평률]](Flattening, $ f $)의 값은 다음과 같다. | 기하학적 관점에서 베셀 타원체는 지구의 자전축을 회전축으로 하는 [[회전 타원체]](oblate spheroid)로 정의된다. 타원체의 형상을 결정하는 핵심 요소는 적도 반지름에 해당하는 [[장반경]](semi-major axis, $ a $)과 극 반지름에 해당하는 [[단반경]](semi-minor axis, $ b $)이다. 베셀이 제시한 수치적 정의에 따르면 장반경과 [[편평률]](flattening, $ f $)의 값은 다음과 같이 설정된다. |
| |
| $$ a = 6,377,397.155 \, \text{m} $$ $$ f = \frac{a - b}{a} = \frac{1}{299.1528128} $$ | $$ a = 6,377,397.155 \, \text{m} $$ $$ f = \frac{a - b}{a} = \frac{1}{299.1528128} $$ |
| |
| 이러한 수치적 정의로부터 단반경 $ b $는 $ b = a(1 - f) $의 관계식을 통해 약 6,356,078.963m로 계산된다. 또한, 측지 계산에서 빈번하게 사용되는 [[제1이심률]](First Eccentricity, $ e $)의 제곱은 다음과 같은 수식을 통해 도출된다. | 이러한 수치적 정의로부터 단반경 $ b $는 $ b = a(1 - f) $의 관계식을 통해 약 6,356,078.963m로 산출된다. 또한, 측지 계산 및 좌표 변환 과정에서 빈번하게 사용되는 [[제1이심률]](first eccentricity, $ e $)의 제곱은 다음과 같은 수식을 통해 도출된다. |
| |
| $$ e^2 = \frac{a^2 - b^2}{a^2} = 2f - f^2 $$ | $$ e^2 = \frac{a^2 - b^2}{a^2} = 2f - f^2 $$ |
| |
| 베셀 타원체의 채택은 단순한 수치적 선택을 넘어, 지표면상의 위치를 수학적 좌표로 변환하기 위한 기하학적 틀을 제공한다는 점에서 중요한 의미를 지닌다. [[동경측지계]]는 이 타원체를 [[경위도 원점]](Geodetic Datum Origin)인 도쿄 아자부다이에서 지오이드 면에 수평이 되도록 설정함으로써, 해당 지역 내에서의 측량 오차를 최소화하도록 설계되었다. 이는 지구 전체의 질량 중심을 기준으로 하는 현대의 [[지구중심 측지계]](Geocentric Datum)와 달리, 특정 지역의 수평면과 타원체면을 일치시키는 [[지역측지계]](Local Geodetic Datum)의 전형적인 접근 방식이다. | 베셀 타원체의 채택은 단순한 수치적 선택을 넘어, 지표면상의 위치를 수학적 좌표로 변환하기 위한 기하학적 틀을 제공한다는 점에서 중요한 의미를 지닌다. 동경측지계는 이 타원체를 [[일본 경위도 원점]](Japanese Geodetic Datum Origin)인 도쿄 아자부다이(麻布台)에서 지오이드 면에 수평이 되도록 설정함으로써, 해당 지역 내에서의 측량 오차를 최소화하도록 설계되었다. 이는 지구 전체의 질량 중심을 기준으로 하는 현대의 [[지구중심 측지계]](geocentric datum)와 달리, 특정 지역의 수직선과 타원체 법선을 일치시키는 [[지역측지계]](local geodetic datum)의 전형적인 접근 방식이다. 따라서 원점에서 멀어질수록 지오이드와 타원체면 사이의 괴리가 커질 수 있으나, 한반도를 포함한 인근 지역에서는 높은 국지적 정밀도를 유지할 수 있었다. |
| |
| 수학적으로 베셀 타원체 위의 임의의 점에 대한 위치는 [[지리 좌표계]](Geographic Coordinate System)인 위도($ $)와 경도($ $)로 표현된다. 이때 타원체 면상의 곡률 반경은 위도에 따라 달라지며, 자오선 곡률 반경($ M $)과 횡곡률 반경($ N $)은 각각 다음과 같은 수식으로 정의되어 정밀한 거리 및 면적 계산의 기초가 된다. | 수학적으로 베셀 타원체 위의 임의의 점에 대한 위치는 [[지리 좌표계]](geographic coordinate system)인 위도($ $)와 경도($ $)로 표현된다. 이때 타원체 면상의 곡률 반경은 위도에 따라 가변적이며, 자오선 곡률 반경($ M $)과 횡곡률 반경($ N $)은 각각 다음과 같은 수식으로 정의된다. |
| |
| $$ M = \frac{a(1 - e^2)}{(1 - e^2 \sin^2 \phi)^{3/2}} $$ $$ N = \frac{a}{\sqrt{1 - e^2 \sin^2 \phi}} $$ | $$ M = \frac{a(1 - e^2)}{(1 - e^2 \sin^2 \phi)^{3/2}} $$ $$ N = \frac{a}{\sqrt{1 - e^2 \sin^2 \phi}} $$ |
| ==== 세계측지계 도입의 필요성 ==== | ==== 세계측지계 도입의 필요성 ==== |
| |
| 20세기 후반 [[위성항법시스템]](Global Navigation Satellite System, GNSS)의 비약적인 발전과 보급은 전통적인 [[측지학]]의 패러다임을 근본적으로 변화시켰다. 과거의 측량 체계가 국가나 특정 지역 단위의 독립적인 기준을 중시하는 [[지역측지계]](Local Geodetic Datum) 중심이었다면, 현대는 지구 전체를 하나의 통합된 좌표계로 관리하는 [[세계측지계]](World Geodetic System)의 시대로 이행하였다. 이러한 변화의 중심에는 [[GPS]](Global Positioning System)를 필두로 한 위성 측위 기술의 보편화와 국제적 데이터 호환성 확보라는 기술적 요구가 자리 잡고 있다. | 20세기 후반 [[위성항법시스템]](Global Navigation Satellite System, GNSS)의 비약적인 발전과 보급은 전통적인 [[측지학]]의 패러다임을 근본적으로 변화시켰다. 기존의 측량 체계가 국가나 특정 지역 단위의 독립적인 기준을 중시하는 [[지역측지계]](Local Geodetic Datum) 중심이었다면, 현대는 지구 전체를 하나의 통합된 좌표계로 관리하는 [[세계측지계]](World Geodetic System)의 시대로 이행하였다. 이러한 변화의 중심에는 [[GPS]](Global Positioning System)를 필두로 한 위성 측위 기술의 보편화와 국제적 데이터 호환성 확보라는 기술적 요구가 자리 잡고 있다. |
| |
| 동경측지계와 같은 지역측지계는 특정 지역의 [[지오이드]](Geoid)에 가장 잘 부합하도록 설정된 [[베셀 타원체]](Bessel 1841 Ellipsoid)를 [[참조 타원체]]로 사용한다. 그러나 이러한 체계는 지구 전체의 형상을 근사하기보다는 해당 지역의 지형적 특성을 반영하는 데 치중하였기에, 타원체의 중심이 [[지구 질량 중심]](Center of Mass)과 일치하지 않는 구조적 한계를 지닌다. 반면, 현대의 위성 항법 기술은 지구의 질량 중심을 원점으로 하는 [[WGS84]](World Geodetic System 1984) 또는 [[ITRF]](International Terrestrial Reference Frame)와 같은 세계측지계를 기반으로 운용된다. 따라서 지역측지계 기반의 지도를 위성 항법 데이터와 직접 결합할 경우, 한반도 지역에서는 약 400미터에서 500미터에 달하는 상당한 수평 위치 편차가 발생하게 된다. | [[동경측지계]]와 같은 지역측지계는 특정 지역의 [[지오이드]](Geoid)에 가장 잘 부합하도록 설정된 [[베셀 타원체]](Bessel 1841 Ellipsoid)를 [[준거 타원체]]로 사용한다. 그러나 이러한 체계는 지구 전체의 형상을 근사하기보다는 해당 지역의 지형적 특성을 반영하는 데 치중하였기에, 타원체의 중심이 [[지구 질량 중심]](Center of Mass)과 일치하지 않는 구조적 한계를 지닌다. 반면, 현대의 위성 항법 기술은 지구의 질량 중심을 원점으로 하는 [[WGS84]](World Geodetic System 1984) 또는 [[ITRF]](International Terrestrial Reference Frame)와 같은 세계측지계를 기반으로 운용된다. 따라서 지역측지계 기반의 지도를 위성 항법 데이터와 직접 결합할 경우, 한반도 일대에서는 남동 방향으로 약 400m에서 500m에 이르는 상당한 수평 위치 편차가 발생하게 된다. |
| |
| 이러한 좌표의 불일치는 단순히 수치상의 차이를 넘어 항행 안전과 직결되는 중대한 문제를 야기한다. 특히 [[국제 민간 항공 기구]](International Civil Aviation Organization, ICAO)와 [[국제 해로 기구]](International Hydrographic Organization, IHO)는 항공 및 해상 안전을 위해 전 세계적으로 통일된 측지 기준의 사용을 강력히 권고하고 있다. 예를 들어, [[항공기]]가 착륙을 위해 GPS 정보를 활용할 때 지상의 활주로 위치가 지역측지계로 작성되어 있다면, 위성이 지시하는 위치와 실제 지형 사이에 치명적인 오차가 발생하여 사고로 이어질 위험이 있다. 이에 따라 항공 분야에서는 이미 오래전부터 WGS84를 표준 좌표계로 채택하여 운용하고 있으며, 이를 국가 전체의 측량 기준으로 확장하는 것은 국제 표준과의 정합성을 확보하기 위한 필수 과제가 되었다. | 이러한 좌표의 불일치는 단순히 수치상의 차이를 넘어 항행 안전과 직결되는 중대한 문제를 야기한다. 특히 [[국제 민간 항공 기구]](International Civil Aviation Organization, ICAO)와 [[국제 해로 기구]](International Hydrographic Organization, IHO)는 항공 및 해상 안전을 위해 전 세계적으로 통일된 측지 기준의 사용을 강력히 권고하고 있다. 예를 들어, [[항공기]]가 착륙을 위해 GPS 정보를 활용할 때 지상의 활주로 위치가 지역측지계로 작성되어 있다면, 위성이 지시하는 위치와 실제 지형 사이에 치명적인 오차가 발생하여 사고로 이어질 위험이 있다. 이에 따라 항공 분야에서는 이미 오래전부터 WGS84를 표준 좌표계로 채택하여 운용하고 있으며, 이를 국가 전체의 측량 기준으로 확장하는 것은 국제 표준과의 정합성을 확보하기 위한 필수 과제가 되었다. |
| |
| 기술적으로 세계측지계 도입은 서로 다른 타원체 간의 기하학적 관계를 규명하는 것에서 시작된다. 지역측지계 좌표를 세계측지계로 변환하기 위해서는 타원체의 중심 이동량인 $ X, Y, Z $와 회전 및 척도 계수를 포함하는 수학적 모델이 동원된다. 이러한 변환 과정은 [[지리 정보 시스템]](Geographic Information System, GIS) 내에서 다양한 출처의 데이터를 통합하는 데 핵심적인 역할을 한다. 과거에는 국가별로 상이한 좌표계를 사용하여 국경 인접 지역의 데이터 통합이 어려웠으나, 세계측지계로의 전환을 통해 범지구적인 지리 공간 정보의 상호 운용성이 확보되었다. | 기술적으로 세계측지계 도입은 서로 다른 타원체 간의 기하학적 관계를 규명하는 것에서 시작된다. 지역측지계 좌표를 세계측지계로 변환하기 위해서는 타원체의 중심 이동량인 $ X, Y, Z $와 회전 및 척도 계수를 포함하는 [[좌표 변환]] 모델이 동원된다. 이러한 변환 과정은 [[지리 정보 시스템]](Geographic Information System, GIS) 내에서 다양한 출처의 데이터를 통합하는 데 핵심적인 역할을 한다. 과거에는 국가별로 상이한 좌표계를 사용하여 국경 인접 지역의 데이터 통합이 어려웠으나, 세계측지계로의 전환을 통해 범지구적인 지리 공간 정보의 [[상호 운용성]]이 확보되었다. |
| |
| 또한, 4차 산업혁명의 핵심 기술인 [[자율 주행]], [[드론]] 항행, [[지능형 교통 체계]](Intelligent Transport Systems, ITS) 등은 센티미터(cm) 단위의 초정밀 위치 정보를 요구한다. 동경측지계와 같이 지역적으로 왜곡이 존재하는 체계로는 이러한 정밀도를 유지하기 어렵다. 따라서 지각 변동을 반영하고 지구 중심과의 정밀한 일치를 보장하는 세계측지계의 도입은 국가 공간 정보 인프라를 현대화하고, 고부가가치 위치 기반 서비스를 창출하기 위한 기술적 토대를 마련하는 과정이라 할 수 있다.((세계측지계의 체계적 적용방안에 관한 연구, https://www.koreascience.or.kr/article/JAKO200821036731213.page?lang=ko | 또한, [[4차 산업혁명]]의 핵심 기술인 [[자율 주행]], [[드론]] 항행, [[지능형 교통 체계]](Intelligent Transport Systems, ITS) 등은 센티미터(cm) 단위의 초정밀 위치 정보를 요구한다. 동경측지계와 같이 지역적으로 왜곡이 존재하는 체계로는 이러한 정밀도를 유지하기 어렵다. 따라서 지각 변동을 반영하고 지구 중심과의 정밀한 일치를 보장하는 세계측지계의 도입은 국가 공간 정보 인프라를 현대화하고, 고부가가치 [[위치 기반 서비스]]를 창출하기 위한 기술적 토대를 마련하는 과정이라 할 수 있다.((세계측지계의 체계적 적용방안에 관한 연구, https://www.koreascience.or.kr/article/JAKO200821036731213.page?lang=ko |
| )) | )) |
| |
| ==== 좌표 변환 계수와 수치적 변환 방법 ==== | ==== 좌표 변환 계수와 수치적 변환 방법 ==== |
| |
| 동경측지계 데이터를 세계측지계로 정밀하게 이전하기 위해 사용되는 수학적 모델을 소개한다. | 동경측지계(Tokyo Datum)와 [[세계측지계]](World Geodetic System) 사이의 좌표 변환은 단순히 서로 다른 [[참조 타원체]](Reference Ellipsoid) 간의 기하학적 차이를 보정하는 것을 넘어, 타원체의 중심 위치와 회전 상태를 일치시키는 복합적인 수리적 과정을 포함한다. 동경측지계는 [[베셀 타원체]](Bessel 1841 Ellipsoid)를 사용하며 그 중심이 [[지구 질량 중심]]에서 약 수백 미터 이격되어 있는 반면, 세계측지계는 지구 중심을 원점으로 하는 [[GRS80]] 또는 [[WGS84]] 타원체를 사용한다. 따라서 두 체계 간의 정밀한 데이터 이전을 위해서는 3차원 공간상의 위치 관계를 규명하는 변환 모델의 수립이 필수적이다. |
| | |
| | 가장 보편적으로 사용되는 수학적 모델은 [[7매개변수 변환]](7-parameter transformation)으로도 불리는 [[부르사-울프 모델]](Bursa-Wolf Model)이다. 이 모델은 두 좌표계 사이의 관계를 세 개의 평행 이동 요소($\Delta X, \Delta Y, \Delta Z$), 세 개의 회전 요소($\epsilon_x, \epsilon_y, \epsilon_z$), 그리고 하나의 축척 계수($s$)로 정의한다. 3차원 직각좌표계 $(X, Y, Z)$에서의 변환식은 다음과 같은 행렬 형태로 표현된다. |
| | |
| | $$ \begin{bmatrix} X_{target} \\ Y_{target} \\ Z_{target} \end{bmatrix} = \begin{bmatrix} \Delta X \\ \Delta Y \\ \Delta Z \end{bmatrix} + (1+s) \begin{bmatrix} 1 & \epsilon_z & -\epsilon_y \\ -\epsilon_z & 1 & \epsilon_x \\ \epsilon_y & -\epsilon_x & 1 \end{bmatrix} \begin{bmatrix} X_{source} \\ Y_{source} \\ Z_{source} \end{bmatrix} $$ |
| | |
| | 여기서 회전각이 매우 작다는 가정하에 회정 행렬의 삼각함수 항을 선형화하여 근사한 형식을 취한다. 대한민국 [[국토지리정보원]]은 과거 동경측지계 기반의 성과를 세계측지계로 전환하기 위해 공통점(Common Points)을 이용한 [[최소제곱법]](Least Squares Method)을 적용하여 국가 표준 변환 계수를 산출한 바 있다((국토지리정보원, 1/1,000 수치지형도 좌표변환 표준작업지침(ver1.0), https://www.ngii.go.kr/kor/contents/view.do?board_code=contents_data&sq=251 |
| | )). |
| | |
| | 지리좌표인 [[경위도]]($\phi, \lambda, h$) 상태에서 직접 변환을 수행할 때는 [[몰로덴스키 공식]](Molodensky Equations)이 활용된다. 이 방식은 3차원 직각좌표로의 번거로운 변환 과정을 생략하고 타원체 제원의 차이($\Delta a, \Delta f$)와 중심점 편차만을 이용하여 좌표를 계산한다. 그러나 몰로덴스키 공식은 회전 요소와 축척 변화를 엄밀하게 반영하지 못하는 한계가 있어, 광역적인 정밀 측량보다는 실시간 항법 등 계산 효율성이 강조되는 분야에서 주로 사용된다((Deakin, R.E. (2006). A Note on the Bursa-Wolf and Molodensky-Badekas Transformations, https://www.sciepub.com/reference/141611 |
| | )). |
| | |
| | 하지만 이러한 수치적 모델만으로는 지역적으로 발생하는 비선형적 왜곡(Distortion)을 완벽히 제거하기 어렵다. 과거의 [[삼각 측량]] 과정에서 누적된 오차나 지각 변동에 의한 미세한 위치 변화는 단순한 기하학적 변환식으로 설명되지 않기 때문이다. 이를 해결하기 위해 [[격자형 수치변환 모델]](Grid-based Transformation)이 도입되었다. 이 방법은 국토를 일정 간격의 격자로 분할하고, 각 격자점에서의 변환 잔차(Residuals)를 보간법으로 계산하여 적용함으로써 변환 정밀도를 극대화한다. 한국에서는 [[NTv2]](National Transformation version 2)와 유사한 격자 변환 방식을 채택하여 지적도 및 국가 기본도의 좌표 전환에 활용하고 있다((GIS 기본도 및 DB의 세계측지계 좌표변환 정확도 분석에 관한 연구, 한국지형공간정보학회지 제16권 제3호, https://scienceon.kisti.re.kr/aiq/mlsh3/pdf.do?cn=JAKO200824359115481&title=GIS+%EA%B8%B0%EB%B3%B8%EB%8F%84+%EB%B0%8F+DB%EC%9D%98+%EC%84%B8%EA%B3%84%EC%B8%A1%EC%A7%80%EA%B3%84+%A2%8C%ED%91%9C%EB%B3%80%ED%99%98+%EC%A0%95%ED%99%95%EB%8F%84+%EB%B6%84%EC%84%9D%EC%97%90+%EA%B4%80%ED%95%9C+%EC%97%B0%EA%B5%AC |
| | )). |
| | |
| | 실제 수치적 변환 과정에서는 변환 전후의 좌표계가 정의하는 [[투영법]]의 차이도 고려해야 한다. 동경측지계 기반의 [[가우스-크뤼거 투영법]] 좌표를 세계측지계 기반의 평면직각좌표로 변환할 때는, 먼저 평면좌표를 타원체상의 경위도로 역투영한 뒤, 위에서 언급한 변환 모델을 거쳐 새로운 타원체상의 경위도를 얻고, 이를 다시 해당 투영법에 맞게 투영하는 다단계 절차를 밟는다. 이 과정에서 발생하는 수치적 오차를 최소화하기 위해 배정밀도 부동소수점 연산이 권장된다. |
| |
| === 상사 변환과 몰로덴스키 변환 기법 === | === 상사 변환과 몰로덴스키 변환 기법 === |
| |
| 7매개변수 변환과 몰로덴스키 공식 등 좌표계 간 이동 및 회전을 계산하는 알고리즘을 다룬다. | [[동경측지계]]에서 [[세계측지계]]로 좌표를 전환하는 과정은 서로 다른 두 [[참조 타원체]] 사이의 기하학적 관계를 수학적으로 모델링하는 작업이다. 두 체계는 타원체의 제원뿐만 아니라 중심점의 위치와 좌표축의 방향이 서로 다르므로, 이를 일치시키기 위한 정밀한 변환 알고리즘이 요구된다. 가장 대표적인 방법으로는 [[3차원 직교좌표계]]를 이용한 [[상사 변환]](Similarity Transformation)과 [[경위도]] 좌표를 직접 이용하는 [[몰로덴스키 변환]](Molodensky Transformation) 기법이 있다. |
| | |
| | 상사 변환의 대표적인 형태인 7매개변수 변환은 [[헬머트 변환]](Helmert Transformation)으로도 알려져 있으며, 두 좌표계 간의 평행 이동량, 회전각, 그리고 축척 계수를 보정하는 방식이다. 이 방법은 먼저 기존의 경위도 좌표와 [[타원체고]]를 [[지심 직교좌표계]](Geocentric Cartesian Coordinate System)의 $X, Y, Z$ 성분으로 변환한 뒤 수행된다. 변환 모델에 포함되는 7개의 매개변수는 세 방향의 평행 이동 성분인 $\Delta X, \Delta Y, \Delta Z$, 각 축을 중심으로 한 회전 성분인 $\epsilon_x, \epsilon_y, \epsilon_z$, 그리고 두 좌표계 사이의 거리 비율을 조정하는 축척 계수 $s$로 구성된다. 이를 행렬식으로 표현하면 다음과 같다. |
| | |
| | $$ \begin{bmatrix} X_{WGS} \\ Y_{WGS} \\ Z_{WGS} \end{bmatrix} = \begin{bmatrix} \Delta X \\ \Delta Y \\ \Delta Z \end{bmatrix} + (1+s) \begin{bmatrix} 1 & \epsilon_z & -\epsilon_y \\ -\epsilon_z & 1 & \epsilon_x \\ \epsilon_y & -\epsilon_x & 1 \end{bmatrix} \begin{bmatrix} X_{Tokyo} \\ Y_{Tokyo} \\ Z_{Tokyo} \end{bmatrix} $$ |
| | |
| | 위 식에서 회전각이 매우 미세할 경우 [[삼각함수]]의 선형화 근사를 통해 계산의 효율성을 높인다. 7매개변수 변환은 공간적 왜곡을 최소화하면서도 물리적인 기하 구조를 잘 반영한다는 장점이 있으나, 변환 계수를 산출하기 위해 두 좌표계의 성과를 모두 가진 다수의 [[공통점]](Common Points)이 확보되어야 하며 [[최소제곱법]](Least Squares Method)을 통한 엄밀한 조정 과정이 수반되어야 한다. |
| | |
| | 반면 몰로덴스키 변환은 지심 직교좌표로의 중간 변환 과정 없이 위도, 경도, 타원체고의 변화량을 직접 계산하는 방식이다. 이 기법은 두 타원체 사이의 장반경 차이($\Delta a$)와 편평률 차이($\Delta f$)를 변환 공식에 직접 포함하여 지리적 좌표의 변화량을 산출한다. 표준 몰로덴스키 공식은 다음과 같은 구조를 가진다. |
| | |
| | $$ \Delta \phi'' = \frac{[-\Delta X \sin \phi \cos \lambda - \Delta Y \sin \phi \sin \lambda + \Delta Z \cos \phi + (a \Delta f + f \Delta a) \sin 2\phi]}{(M + h) \sin 1''} $$ $$ \Delta \lambda'' = \frac{[-\Delta X \sin \lambda + \Delta Y \cos \lambda]}{(N + h) \cos \phi \sin 1''} $$ |
| | |
| | 여기서 $M$은 [[자오선 곡률 반경]]을, $N$은 [[동서 곡률 반경]]을 의미하며, $h$는 타원체고이다. 몰로덴스키 변환은 계산 과정이 상대적으로 단순하여 과거 계산 자원이 한정적이었던 시기에 널리 사용되었으며, 특히 고도 정보가 불분명한 지역에서 [[수평 위치]]를 변환할 때 유용하다. 하지만 상사 변환에 비해 좌표축의 회전이나 국지적인 왜곡을 정밀하게 반영하는 데 한계가 있어, 현대의 고정밀 [[국가측지계]] 변환 업무에서는 7매개변수 모델이나 이를 보완한 [[격자 기반 변환]](Grid-based Transformation) 기법이 주로 채택된다. |
| | |
| | 대한민국에서는 동경측지계에서 세계측지계로 전환하기 위해 [[국토지리정보원]]에서 산출한 국가 좌표 변환 계수를 사용한다. 이 계수는 전국에 분포한 [[국가기준점]]의 성과를 분석하여 도출되었으며, 지역적 편차를 극복하기 위해 전국을 단일 계수로 처리하는 방식과 더불어 왜곡량을 보정하는 수치 모델이 병행 적용된다. 이러한 변환 기법의 적용을 통해 과거 [[베셀 타원체]] 기반의 데이터를 현대적인 [[GRS80]] 타원체 체계로 오차 없이 통합할 수 있다. ((국토지리정보원, 세계측지계 변환 기술지침, https://www.ngii.go.kr/kor/content.do?sq=545 |
| | )) |
| |
