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.