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

用法