geodesy op• 資料種類:depth → image2d
• 呼叫: import fullseye as fs; fs.ledger.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees')(要直接呼叫實作,import demops; demops.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees');從台帳取用則 opsdem.get("dem_geodetic_slope"))
經緯度網格(等角度間隔)DEM 的坡度。格子單元尺寸隨緯度變化。
公開的 DEM 大多像「1 秒網格(mesh)」那樣按角度等間隔。把這種網格
直接以固定單元尺寸傳給 dem_slope,就會忽略經度方向的實際距離短了 cos(緯度) 倍這一點。若兩軸都用緯度方向的尺寸, 東西向坡度會偏小(北緯 60 度時減半);若兩軸都用經度方向的尺寸, 則反過來南北向偏大。偏向哪一邊取決於把什麼當作常數,所以正確做法不是記住「偏大/偏小」,而是按緯度計算尺寸。
★ 最初這個 docstring 寫的是「東西向會偏大」,但用測試量測後 變成了一半(偏小)。關於符號或方向的主張要先量測再寫。
這裡按緯度計算東西向的單元尺寸,把網格修正為等距離後再量測。
Args:
dem: (H, W)。第 0 列為北。
lat0_deg: 西北角的緯度 [度]。
d_lat_deg / d_lon_deg: 每個單元的緯度、經度步長 [度]。正值。
method / units: 與 dem_slope 相同。
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• 演算法的正典(作者、年份)與用途見上面的族使用指南。
• dem_geodesy_tour — py -3.11 examples/dem_geodesy_tour.py
image2d 作為輸入)—
geodesy)dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_geocentric_grid · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geoid_height · dem_height_frame_convert · dem_height_frame_residual
*Provenance: demops.py — DEM 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.