[논문 리뷰] Geodesics on an ellipsoid of revolution

Charles F. F. Karney|arXiv (Cornell University)|2009. 03. 01.

Historical Geography and CartographySocial Sciences참고 문헌 25인용 수 21

한 줄 요약

이 논문은 타원체의 회전축을 기준으로 하는 직접 및 역지오데식 문제를 해결하기 위한 매우 정확하고 반복이 없는 알고리즘을 제시하며, 나노미터 이내의 정밀도를 달성한다. 보조 구면과 야코비 타원 함수를 사용한 새로운 수식을 도입함으로써 삼각측량, 지ap로프로젝션, 해양 경계, 다각형 면적 계산에 대한 강력한 해법을 가능하게 하며, 수렴 보장과 높은 수치 안정성을 확보한다.

ABSTRACT

Algorithms for the computation of the forward and inverse geodesic problems for an ellipsoid of revolution are derived. These are accurate to better than 15 nm when applied to the terrestrial ellipsoids. The solutions of other problems involving geodesics (triangulation, projections, maritime boundaries, and polygonal areas) are investigated.

연구 동기 및 목표

• 타원체의 회전축에 대해 직접 및 역지오데식 문제를 수치적으로 안정적이고 매우 정확하게 해결하기 위한 것이다.
• 기존 알고리즘(특히 비ン센티의 방법)의 한계를 극복하기 위해 고정된 순도 의존성 제거 및 거의 반대점에서의 수렴 실패 문제 해결을 위한 것이다.
• 지오데식 해법을 긴장 타원체로 일반화하고 삼각측량, 지도 투영, 해양 경계, 그리고 구면 다각형 면적 계산 응용 분야로 확장하기 위한 것이다.
• 지오데시브 라이브러리에서 사용 가능한 감소 길이 및 기타 고급 지오데틱 개념을 통합한 종합적인 이론적 및 계산 프레임워크를 제공하기 위한 것이다.

제안 방법

• 등각 사영을 통해 타원체 문제를 해결 가능한 구면 삼각형 문제로 변환하는 보조 구면 방법을 사용하여 지오데식을 수식화한다.
• 야코비 타원 함수(sn, dn, am)와 그 적분(F, E, G)을 사용하여 지오데식 호의 길이, 방위각 및 위치를 고정밀도로 표현한다.
• 편평함 $f$, 제3편평함 $n$, 또는 이심률 제곱 $e^2$에 대한 급수 전개를 사용하여 직접 문제를 반복 없이 해결함으로써 15 nm 이내의 정확도를 확보한다.
• 뉴턴 방법을 사용하여 역문제를 해결하며, 초기 추정치를 신중히 선택함으로써 거의 반대점인 경우에도 몇 번의 반복 내에 수렴을 보장한다.
• 지오데식 삼각형에서 핵심 불변량인 감소 길이 $m_{12}$ 개념을 도입하여 정확한 삼각측량 및 면적 계산을 가능하게 한다.
• 지오데식 거리와 방위각을 사용하여 기하학적 투영 및 구면 다각형 면적 공식을 타원체 경우로 일반화함으로써 정밀도를 유지한다.

실험 결과

연구 질문

• RQ1 타원체의 회전축에 대해 직접 및 역지오데식 문제를 나노미터 이내 정밀도와 보장된 수렴을 확보하여 어떻게 해결할 수 있는가?
• RQ2 지오데식 알고리즘을 긴장 타원체로 확장하기 위해 어떤 수정이 필요한가? 그리고 타원체의 경우와 긴장 타원체의 경우 해법은 어떻게 다를까?
• RQ3 지오데식 기반 방법을 사용하여 구면 다각형 면적을 고정밀도로 계산하고 해양 경계를 정의하기 위해 어떻게 활용할 수 있는가?
• RQ4 타원 적분과 야코비 함수는 직접 지오데식 문제의 반복 없는 해법을 가능하게 하기 위해 어떤 역할을 하는가?
• RQ5 감소 길이 개념은 삼각측량 및 지오데식 네트워크 계산의 정확도와 안정성을 어떻게 향상시킬 수 있는가?

주요 결과

• 직접 지오데식 문제는 $f$, $n$, 또는 $e^2$에 대한 급수 전개를 사용하여 반복 없이 해결되며, 타원체의 경우와 긴장 타원체의 경우 모두 15 nm 이내의 정확도를 확보한다.
• 뉴턴 방법을 사용하여 역지오데식 문제는 거의 반대점인 경우에도 몇 번의 반복 내에 안정적으로 수렴하며, 비ン센티 알고리즘의 핵심 결함을 극복한다.
• 감소 길이 $m_{12}$는 지오데식 네트워크에서 정확한 삼각측량 및 면적 계산을 가능하게 하는 기본 불변량이다.
• 긴장 타원체의 경우, 초기 방위각 $\alpha_1$에 대해 경도 차이 $\lambda_{12}$는 더 이상 단조적이지 않으며, 뉴턴 방법의 시작점 선택을 수정할 필요가 있다.
• 지오데식 다각형 면적은 $2\pi n - \sum E_{i-1,i}$로 계산되며, $n$은 전체 동쪽 방향 둘레 수를 나타내며, 고정밀도 면적 통합을 보장한다.
• 긴장 경우를 위한 지오데식 면적 적분은 $G(\phi, \alpha^2, k)$와 복소 타원 함수를 사용하여 재구성되며, $c^2$는 $\tan^{-1}(\sqrt{-e^2})/\sqrt{-e^2}$를 통해 정의되어 실수 표현을 유지한다.