refine op• 資料種類:voxel × voxel × position → position
• 呼叫: import fullseye as fs; fs.ledger.refine_translation_lk(scene, template, init_pos, device='cpu', iters=30, tol=0.0001)(要直接呼叫實作,import match3d; match3d.refine_translation_lk(scene, template, init_pos, device='cpu', iters=30, tol=0.0001);從台帳取用則 ops3d.get("refine_translation_lk"))
• GPU:該運算子有 GPU 路徑(device="cuda")
用 Gauss-Newton 逆合成 Lucas-Kanade 做 3D 平移的次體素精化。
> 以下的詳細說明為原文 —— 摘要與標題已翻譯。
粗マッチ(整数 NCC / Fourier-Mellin / Hough)が与えた整数初期位置 `init_pos` を
出発点に、SSD `Σ|I(x+p) − T(x)|² を最小化してサブボクセル並進 p` へ収束させる。
逆合成(inverse-compositional, Baker–Matthews)方式のため steepest-descent 画像
`SD = ∇T と Hessian H = Σ SDᵀSD` を 反復前に一度だけ前計算し、各反復は
「scene の trilinear ワープ + 残差 + 3×3 線形解 `Δp = H⁻¹ Σ SDᵀ(I(x+p)−T)`」のみ。
並進の合成は `p ← p − Δp`。純並進ワープでは ∂W/∂p=I なので SD=∇T がそのまま使える。
座標系: `init_pos` と戻り値はいずれも テンプレート原点(corner, index 0,0,0) が
scene のどの (dz,dy,dx) に載るか。`sobel3d / grid_sample` の corner 規約に一致
(NCC(ncc_locate_3d)の中心規約とは T//2 だけ異なる点に注意)。
Parameters
----------
scene : (D,H,W) array_like
探索対象ボリューム。
template : (Td,Th,Tw) array_like
位置合わせするテンプレート(scene より小)。
init_pos : (3,) sequence
整数初期位置 (dz,dy,dx) = テンプレート原点の scene 座標。
device : str
"cpu" / "cuda" 等。device 非依存。
iters : int
最大反復数。
tol : float
‖Δp‖ がこの値を下回ったら収束打ち切り。
Returns
-------
pos : (3,) np.ndarray(float64)
精緻化されたサブボクセル位置 (dz,dy,dx)。
Notes
-----
• 滑らか(帯域制限)な密度場を仮定。整数初期値が真値の ±0.5〜1 voxel 内であれば
通常 5〜8 反復で ‖err‖ < 0.05 voxel(低ノイズ時)。実測(独立 cubic-spline GT):
ノイズ無し mean 0.008 / max 0.013 voxel(≈6 反復, ≈1.2ms/回)、NCC サブボクセル
baseline(mean 0.56 voxel)を約60×改善。
• `grad_scale=32` は分離 sobel3d(導関数[-1,0,1]×平滑[1,2,1]²)の固定スケール
(線形ランプで実測 32.0)。真の勾配へ正規化して Δp のスケールを正す。
• H には微小 Levenberg 正則化を加え、勾配の乏しい平坦テンプレートでの数値破綻を防ぐ。
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• refinement — py -3.11 examples_3d/refinement.py
position 作為輸入)refine_peak_newton · refine_lm · fuse_to_voxel
refine)refine_peak_newton · refine_lm · refine_rotation_z · icp_point2point_3d · icp_point2plane
*Provenance: match3d.py — 3D 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.