🐥 절기(월주 경계) 정밀 계산 도입 — 데이터 확보부터 3중 독립 검증까지

1. 배경 — 왜 이 조사가 시작됐나

신한라이프 평생운세 벤치마크 비교 리포트에서 무작위 생년월일시 10건을 대형 금융사의 무료 사주 서비스와 우리 엔진에 나란히 돌렸을 때, 년·일·시주는 10/10 전부 일치했지만 월주만 1건(1999-10-08 축시, 남) 갈렸다. 신한은 유월(酉月)로, 우리는 술월(戌月)로 계산해 정확히 한 칸씩 밀린 결과였다.

2. 근본 원인 — 고정 근사 날짜 방식의 한계

src/manseryeok.js의 월주 경계표(MONTH_BOUNDARIES)는 12절기 절입일을 매년 똑같은 고정 날짜로 근사한다(예: 한로 = 매년 10/8). 하지만 실제 절입 시각은 지구 공전 특성상 매년 몇 시간~하루 정도씩 미세하게 달라진다 — 1999년 실제 한로는 10/9였다. 경계 ±1일 근처 출생자만 영향을 받고(전체의 대략 2% 안팎으로 추정), 나머지 대다수 날짜는 원래도 정확했다.

3. 1단계 — KASI 정부 공식 데이터 확보와 실제 범위 발견

한국천문연구원(KASI)이 data.go.kr에 공개한 "특일 정보" API(데이터셋 15012690, get24DivisionsInfo 엔드포인트)는 24절기 절입 날짜·시각·태양황경을 정부 공식으로 제공한다. 문서만으로는 지원 연도 범위가 명시돼 있지 않아서, 연도를 하나씩 직접 호출해서 실측했다.

원래 기대했던 "1900~2050년 전체 커버"와는 다르게, 이 API는 "올해·최근 몇 년" 조회용으로 설계된 것으로 보이는 딱 29년 구간(2000~2028)만 지원했다. 이 구간은 실측 데이터를 PRECISE_TERMS_BY_YEAR 테이블에 그대로 내장해서, 매 요청마다 외부 API를 부르지 않고도 (외부 API 장애가 핵심 서비스에 영향 주지 않도록) 쓸 수 있게 했다.

4. KASI 원본 데이터에서 발견한 오류 2건

정부 공식 데이터라고 해서 무결하지는 않았다. 교차검증 과정에서 KASI 원본 데이터 자체의 오류 2건을 발견했다.

5. 2단계 — 자체 천문 계산식 구현과 1차 검증

KASI API가 지원하지 않는 1900~1999·2029~2050년을 채우기 위해, 태양의 겉보기 황경(apparent geocentric ecliptic longitude)을 직접 계산하는 천문 공식을 구현했다:

• 태양황경: Meeus 저정밀 공식(평균 황경 L0 + 평균 근점이각 M + 중심차 C, 章동·광행차 보정 포함)
• ΔT(TT-UT 보정): Espenak/Meeus 근대 구간 구분 다항식(1900~2050년대별로 다른 근사식 적용)
• 절입 시각 탐색: 목표 황경(15도 간격, 24절기)에 도달하는 순간을 ±7일 구간 안에서 이분탐색(60회 반복)으로 특정

이 계산식을 KASI 실측 데이터가 있는 2000~2028년 구간(24절기 × 29년 ≈ 695개 포인트)과 전수 대조해서 1차 검증했다.

6. 3단계 — NASA JPL 정밀 천체력으로 독립 재검증

5단계까지의 검증은 2000~2028년 구간에서만 이뤄졌다는 한계가 있었다 — 정작 자체 계산식이 실제로 채워야 하는 1900~1999·2029~2050년(외삽 구간)은 "공식이 그 안에서도 잘 맞을 것"이라는 가정에 의존했을 뿐, 독립적으로 확인된 적이 없었다. 이 공백을 메우기 위해 KASI와 전혀 무관한 제3의 권위있는 출처로 다시 검증했다.

NASA JPL이 발행하는 DE421 정밀 천체력(태양계 천체의 실제 수치적분 궤도력, 아크초 단위 정밀도)을 Python skyfield 라이브러리로 로컬에 내려받아 직접 계산했다(외부 API 호출이 아니라 다운로드한 궤도력 파일로 자체 계산). 1900~2050년 전체 151년 × 12절 = 1,812개 포인트를 전수 계산해서 대조했다.

2000~2028년 근방에서만 검증됐던 이전과 달리, 외삽 구간(1900~1999·2029~2050)을 포함한 전체 150년이 실제로 분 단위 정확도라는 게 KASI와 독립적인 제3의 출처로 다시 확인됐다.

2011년 대한 KASI 오류, 3중 확인

§4에서 "KASI 원본 오류로 추정"했던 2011년 대한 24시간 불일치 건을 JPL로 다시 확인했다.

다만 대한(大寒)은 애초에 월주 경계로 쓰는 12개 "절"에 포함되지 않는 절기라, 이 오류는 실제 서비스 계산에는 영향을 준 적이 없었다.

7. 최종 확인 — 원래 버그가 실제로 고쳐졌는가

이 조사를 촉발한 §1의 원래 케이스를 배포된 최종 코드로 다시 계산해서 확인했다.

추가로 1955-02-04(입춘 경계 케이스), getTodayPillar()(오늘의 운세 실시간 계산), 1900~2050년 604개 표본 사주 전수 계산에서 크래시 0건을 확인했다.

8. 남은 한계와 배포 현황

• 배포 완료: src/manseryeok.js의 PRECISE_TERMS_BY_YEAR에 1900~2050년 전체 반영, SAJU_CANON_VERSION = "2026-07-22-precise-terms-1900-2050", ./deploy.sh로 프로덕션 배포 + git 커밋 완료.
• 유일하게 남은 한계: 1961년 이전 한국 표준시는 UTC+8:30이었고(1961년부터 현재의 UTC+9), 1948~51·55~60·87~88년엔 서머타임도 있었다. 지금 계산은 전 구간을 KST=UTC+9로 고정 가정한다 — 이건 절기 계산 정확도 문제가 아니라 "그 시각을 몇 시(KST)로 환산하느냐"의 별개 이슈로, 아직 미해결 상태로 남아있다.

조사 타임라인

• 1단계
  신한라이프 벤치마크 비교에서 1999-10-08 월주 불일치 발견 → 원인이 고정 근사 절기표라는 것 확인
• 2단계
  KASI 정부 API 인증키 발급 → 연도별 실측으로 2000~2028년만 지원한다는 것 발견 → 그 구간 우선 반영·배포
• 3단계
  "2000년 이전도 마저 해결"이라는 요청에 따라 자체 천문 계산식(Meeus + ΔT) 구현 → KASI 실측 695포인트와 대조, 평균오차 6.97분 확인 → 1900~2050년 전체 반영·배포
• 4단계
  "다른 방법으로 데이터 확보 가능한지 조사" 요청에 따라 NASA JPL DE421 천체력으로 KASI와 무관한 독립 재검증 → 1,812포인트 중 4건만(전부 자정 경계) 불일치 → 외삽 구간까지 포함해 전체 검증 완료, 2011년 KASI 오류 3중 확인