geodesy op• 데이터 종류: 없음 → points(인자만으로 정해지는 연산자 —— 이미지나 데이터 입력을 받지 않습니다)
• 호출: import fullseye as fs; fs.ledger.dem_datum_shift_3param(lat_deg, lon_deg, h_m, dx_m, dy_m, dz_m, a_from_m, f_from, a_to_m, f_to)(구현을 직접 호출하려면 import demops; demops.dem_datum_shift_3param(lat_deg, lon_deg, h_m, dx_m, dy_m, dz_m, a_from_m, f_from, a_to_m, f_to), 원장에서 가져오려면 opsdem.get("dem_datum_shift_3param"))
측지 성과(datum)의 변환 —— 지심 3 매개변수 평행 이동(points).
같은 “북위 35도 41분”이라도, 구 일본 측지계와 세계 측지계에서는 지상에서 수백 미터 떨어진 점을 가리킨다. 이 op는 그 전환을 지심 직교 좌표에서의 평행 이동(Molodensky-Badekas의 회전·축척 없는 경우)으로 수행한다::
ECEF(a_from, f_from) + (dx, dy, dz) -> 측지 좌표(a_to, f_to)
★파라미터에 기본값을 두지 않는다. 어느 성과에서 어느 성과로, 어느 평행 이동량으로옮기는지는 나라·지역·판마다 다르며, 말없이 가정하면 “예외는 나지 않지만 수백 미터어긋난 좌표”가 나온다 —— 그것이 이 op의 존재 이유이므로, 호출하는 쪽이 반드시 쓰게 한다.
Args:
lat_deg / lon_deg: 원래 성과에서의 위도·경도 [도]. 배열 가능.
h_m: 원래 성과에서의 타원체고 [m].
dx_m / dy_m / dz_m: 지심의 평행 이동량 [m] (원래 → 대상).
a_from_m / f_from: 원래 타원체의 장반경 [m]과 편평률.
a_to_m / f_to: 대상 타원체의 장반경 [m]과 편평률.
Returns:
(n, 3)의 (위도[도], 경도[도], 타원체고[m]).
닫힌 식으로 검사할 수 있는 것: 평행 이동 0이고 같은 타원체이면 항등 변환 (왕복의 하한은 dem_ecef_to_geodetic과 같다). 평행 이동만 주었을 때의지상 이동량은, 그 점의 ECEF 기저로 분해한 성분과 일치한다.
ValueError: 비유한, 위도가 ±90을 벗어남, 장반경이 양수가 아님, 편평률이 [0, 1)을 벗어남.
• dem_terrain_analysis 패밀리 가이드
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.
• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.
• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.
• poc_geodetic_benchmarks_real — py -3.11 examples/poc_geodetic_benchmarks_real.py
• poc_geodetic_height_frames — py -3.11 examples/poc_geodetic_height_frames.py
points 를 입력으로 받는 것)dem_ecef_to_geodetic · dem_geodetic_from_enu
geodesy)dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_geocentric_grid · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope · dem_geoid_height · dem_height_frame_convert
*Provenance: demops.py — DEM 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.