dem_skyline — GEOCAM `skyline` op

• Data kinds: depth → table

• Call: 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) (to call the implementation directly, 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); from the registry, opsgeocam.get("dem_skyline"))

Usage

> This operator's description has not been translated yet. The original text follows as it is.

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 未満、昇順。

Detailed usage guide

• geocam family guide

References (sample data, literature)

• Sample-data catalog (download URLs / licences) — 2-D uses skimage.data (BSD/public domain) plus synthetic images; 3-D lists download URLs for real data sources (Stanford, PDS, …).

• Operator provenance and references — the sources of the research/methods this op family came from.

• The canonical algorithm (author, year) and its uses are named in the family usage guide above.

Runnable examples (verified samples that actually call this op)

• poc_public_camera_heading — py -3.11 examples/poc_public_camera_heading.py

Ops the type connects to (they accept `table` as input)

render_skyline_view · camera_orientation_from_skyline

Same category (`skyline`)

skyline_extract · render_skyline_view

*Provenance: geocam.py — GEOCAM operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*

dem_skyline — GEOCAM skyline op