지평선 합성으로 고고천문학 연구 지원

초록

본 논문은 디지털 고도 모델(DEM)을 이용해 관측지점에서 보이는 자연 지평선을 수치적으로 합성하는 간단한 절차를 제시한다. 지구 곡률과 대기 굴절을 보정한 후, 지정된 방위각마다 LOS(시선 경로) 고도 프로파일을 추출하고 최대 고도를 지평선 고도로 정의한다. SRTM 90 m DEM을 활용한 사례와 현장 측량(시돌림) 결과를 비교해 평균 0.1° 수준의 RMS 오차를 보이며, 10 km 이상 떨어진 지평선에 대해 신뢰할 수 있음을 확인한다.

상세 분석

이 연구는 고고천문학에서 건축물·무덤·제단 등의 방위 해석에 필수적인 자연 지평선 정보를, 전통적인 현장 측량 대신 디지털 고도 모델을 활용해 자동화·표준화하는 방법을 제시한다. 핵심 아이디어는 관측지점 P₀에서 특정 방위각 ϕ를 따라 일정 간격(Δl)으로 지표면 위의 점들을 샘플링하고, 각 점 B의 고도 h_B와 거리 d를 이용해 지구 곡률을 보정한 시선 고도 α를 계산한다(식 1). 여기서 사용된 곡률 보정식은 타원체 반경 R을 기반으로 한 1/2 d/R 항을 포함해, 300 km 이내에서는 0.01° 이하의 오차를 보장한다.

대기 굴절 보정은 Bomsford(1980)의 모델을 채택했으며, 굴절 상수 κ를 기압·온도·온도 구배에 따라 정의한다(식 2·3). 일반적인 대기 조건에서는 κ≈2.3″ km⁻¹(≈0.00064° km⁻¹)이며, 100 km 거리에서는 약 0.06°의 상승 효과를 제공한다. 이는 곡률 보정보다 약 7배 작지만, 최종 정확도 한계에 크게 기여한다.

DEM 활용 측면에서는 SRTM 90 m 데이터를 기본으로 삼았으며, 수직 RMS 오차 σ_z≈1.8 m가 1.5 km 거리에서 약 0.1°의 고도 오차로 전환된다는 계산을 제시한다(σ_α≈180√2 π σ_z/d). 수평 정확도는 20 m 이하이며, 이는 방위각 오차 σ_ϕ≈0.1°(d=10 km) 수준을 만든다. 따라서 10 km 이상 떨어진 지평선에 대해 방위와 고도 모두 0.1° 이내의 정밀도를 기대할 수 있다.

알고리즘은 다음 단계로 구성된다: (1) 관측지점 고도 z₀ 추출, (2) 방위각 ϕ와 거리 d_i에 따라 지오데식 전진 문제를 수치적으로 해결해 점 P_i 좌표 계산, (3) DEM에서 해당 점의 고도 z_i 추출, (4) 식 A4(곡률 보정 포함)와 굴절 보정을 적용해 시선 고도 α_i 산출, (5) α_i 중 최대값을 해당 방위의 지평선 고도로 정의. 이 과정을 0°–360° 전 방위에 걸쳐 Δϕ 간격으로 반복하면 전체 지평선 프로파일 α(ϕ)를 얻는다.

검증 사례로는 이탈리아 Aquileia와 St. Martin(Artegna)에서 직접 시돌림 측정과 합성 결과를 비교하였다. 멀리(>10 km) 떨어진 서쪽 지평선은 두 방법 간 차이가 0.1° 이하로 일치했으며, 가까운(≈1–2 km) 동쪽 지평선은 DEM의 평활화 효과로 고도가 낮게 추정돼 0.5°~1.5° 차이가 발생했다. 이는 DEM 해상도가 지형 급경사를 충분히 포착하지 못함을 보여준다. 따라서 10 km 미만의 근거리 지평선은 고해상도(30 m 이하) DEM이나 현장 측량이 필요하다.

마지막으로 스코틀랜드 Maes Howe 사례를 제시해, 합성된 지평선이 고대 구조물의 일출·일몰 정렬 해석에 어떻게 활용될 수 있는지를 설명한다. 전체적으로 이 방법은 전 세계 어디서든 공개 DEM을 이용해 빠르고 일관된 지평선 데이터를 제공함으로써, 고고천문학 연구의 효율성과 재현성을 크게 향상시킨다.