camera_orientation_from_skyline — GEOCAM orientation op

• 資料種類:signal × table → table

• 呼叫: import fullseye as fs; fs.ledger.camera_orientation_from_skyline(skyline_rows, K, sky_table, yaw_step_deg=1.0, pitch_range=(-15.0, 15.0), roll_range=(-10.0, 10.0), coarse_step_deg=2.5)(要直接呼叫實作,import geocam; geocam.camera_orientation_from_skyline(skyline_rows, K, sky_table, yaw_step_deg=1.0, pitch_range=(-15.0, 15.0), roll_range=(-10.0, 10.0), coarse_step_deg=2.5);從台帳取用則 opsgeocam.get("camera_orientation_from_skyline"))

用法

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

写真のスカイライン(列ごとの行)と DEM のスカイラインから (yaw, pitch, roll) を決める → table。

写真の境界画素 (u, v) をカメラ光線にし、仮の姿勢で世界へ回すと、各列は (方位, 仰角) の点になる。

その仰角と、DEM のスカイラインをその方位で引いた仰角の差の RMS が姿勢のコスト。yaw を

`yaw_step_deg 刻みで一周、pitch・roll を coarse_step_deg` 刻みで格子探索し、最小の点から

Nelder–Mead で連続値に詰める(scipy)。Baatz ら 2012 のスカイライン照合を、位置既知・

1 台のカメラに絞った形。

曖昧さを返す: yaw ごとの最良コストの曲線(`yaw_profile_deg / yaw_profile_residual_deg`)

と、最良と 2 番目の谷の差 `margin_deg`(2 番目の谷 = 最良から 20° 以上離れた最小)。平地・

対称な尾根・半円しか写らない写真では谷が複数あり、margin が小さい。**margin が

`coarse の残差程度なら向きは決まっていない**と読む(黙って 1 つを返さない: ambiguous`)。

Args:

skyline_rows: (W,) `skyline_extract` の出力(列ごとの行)。

K: (fx, fy, cx, cy)。

sky_table: `dem_skyline` の出力。

yaw_step_deg / coarse_step_deg: 粗い格子の刻み。pitch_range / roll_range: 探索範囲 [度]。

Returns:

table: `yaw_deg / pitch_deg / roll_deg / residual_deg`(仰角残差 RMS)/

`margin_deg / runner_up_yaw_deg / ambiguous(bool)/ yaw_profile_deg` (M,) /

`yaw_profile_residual_deg (M,) / n_columns`。

詳細使用指南

• 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

同類別(orientation)

—


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

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.