geodesy op• 데이터 종류: depth → signal
• 호출: import fullseye as fs; fs.ledger.dem_geoid_height(geoid, lat_deg, lon_deg, lat0_deg, lon0_deg, d_lat_deg, d_lon_deg)(구현을 직접 호출하려면 import demops; demops.dem_geoid_height(geoid, lat_deg, lon_deg, lat0_deg, lon0_deg, d_lat_deg, d_lon_deg), 원장에서 가져오려면 opsdem.get("dem_geoid_height"))
지오이드고 격자를 위도·경도로 쌍선형 보간하여 N [m]을 반환한다(signal).
GNSS가 반환하는 것은 타원체고 h, 지도와 설계도가 쓰는 것은 표고 H. 그 차이가지오이드고 N(H = h - N)이며, 일본 부근에서는 대략 +30~+40 m 있다. 이 op는 공개된 지오이드 모델의 격자(GSIGEO / EGM 계열 같은등간격 그리드)를 읽는 쪽의 도구로, 모델 자체는 갖지 않는다 —— 어느 판을썼는지는 호출자의 책임이며, 판이 바뀌면 표고는 수 cm 움직인다.
Args:
geoid: (H, W)의 지오이드고 [m]. geoid[0, 0]이
(lat0_deg, lon0_deg), 행이 늘면 위도가 d_lat_deg 늘고,
열이 늘면 경도가 d_lon_deg 는다(어느 부호든 좋다).
lat_deg / lon_deg: 조회할 점의 위도, 경도 [도]. 배열 가능.
lat0_deg / lon0_deg: 격자의 원점 [도].
d_lat_deg / d_lon_deg: 격자의 스텝 [도]. 0은 거부한다.
Returns:
(n,)의 지오이드고 [m].
닫힌 식으로 검사할 수 있는 것: N이 위도, 경도의 1차식인 격자에대해서는, 쌍선형 보간은 엄밀(기계 정밀도). 2차 항이 있을 때의 오차는 |∂²N| * d² / 8로 위에서 억제된다.
ValueError: 격자가 2-D가 아님 / 2x2 미만 / 비유한, 스텝이 0이거나 비유한, 위도가 ±90을 벗어남, 그리고 격자 밖의 점(해역이나 범위 밖을 조용히 끝 값으로메우지 않는다 —— 외삽한 지오이드고는 측량값이 아니다).
• 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
signal 를 입력으로 받는 것)dem_height_frame_convert · dem_height_frame_residual
geodesy)dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_geocentric_grid · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope · dem_height_frame_convert · dem_height_frame_residual
*Provenance: demops.py — DEM 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.