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 正則化を加え、勾配の乏しい平坦テンプレートでの数値破綻を防ぐ。
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL 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 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.