dem_skyline — GEOCAM `skyline` op

• 資料種類:depth → table

• 呼叫: import fullseye as fs; fs.ledger.dem_skyline(dem, cell_size, observer_rc, eye_height=1.5, az_step_deg=1.0, max_distance=None, earth_curvature=True)(要直接呼叫實作,import geocam; geocam.dem_skyline(dem, cell_size, observer_rc, eye_height=1.5, az_step_deg=1.0, max_distance=None, earth_curvature=True);從台帳取用則 opsgeocam.get("dem_skyline"))

用法

> 該運算子的說明尚無譯文,以下照原文給出。

DEM の 1 点から見た全方位のスカイライン(方位ごとの地平線仰角)→ table。

`dem_horizon_angle` は「全セル × 1 方位」、これは「1 点 × 全方位」—— 固定カメラの向きを

決めるのに要るのは後者。各方位へ `cell_size / 2` 刻みで視線を進め(標高は双一次補間)、

`atan((z_j − z_0 − drop(d)) / d) の最大値をその方位の仰角にする。drop(d) = d² / (2 R_eff)`

は地球の丸みと大気屈折(有効半径 R / (1 − 0.13))で遠くの山が沈む分。視線が DEM の外に

出たらそこで止める(`max_distance` でも止まる)。DEM の外側は「DEM の最低標高の平面が

地平線まで続く」と仮定し、仰角の下限を地平線の沈み `−sqrt(2 h / R_eff)`(h = 観測点の

最低標高からの高さ)にする —— 有限の DEM の端でたまたま低い点を見て、実際には地球の丸みで

見えない −7° のような「穴」を空にしないため(`horizon_dip_deg` に返す)。

Args:

dem: (H, W) 標高 [m]。行 0 が北端。

cell_size: [m/セル]。

observer_rc: 観測点 (row, col)。float 可(セルの中に立てる)。

eye_height: 地面からのカメラ高さ [m]。

az_step_deg: 方位の刻み [度]。

max_distance: 視線の最大距離 [m] (None = DEM の端まで)。

earth_curvature: False なら平らな地球。

Returns:

table: `azimuth_deg (M,) / elevation_deg (M,) / distance_m` (M,)(仰角を決めた

地点までの距離)/ `observer_elevation_m / horizon_dip_deg / n`。方位は 0 以上

360 未満、昇順。

詳細使用指南

• geocam 族使用指南

參考(範例資料・文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者・年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• poc_public_camera_heading — py -3.11 examples/poc_public_camera_heading.py

型別可銜接的下一個運算子(可接受 `table` 作為輸入)

render_skyline_view · camera_orientation_from_skyline

同類別(`skyline`)

skyline_extract · render_skyline_view

*Provenance: geocam.py — GEOCAM 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

dem_skyline — GEOCAM skyline op

用法