match_phase_3d — 3D match_pose op

Data kinds: voxel × voxelshift

Call: import fullseye as fs; fs.ledger.match_phase_3d(a, b, device='cpu') (to call the implementation directly, import match3d; match3d.match_phase_3d(a, b, device='cpu'); from the registry, ops3d.get("match_phase_3d"))

GPU: this op has a GPU path (device="cuda")

Usage

3-D phase correlation (FFT). Returns the integer shift (dz, dy, dx) that aligns b to a.

> The detailed description below is the original text — the summary and the headings are translated.

Reddy & Chatterji の 3D 版。相互パワースペクトルの逆 FFT のピーク = 平行移動。テンプレート

不要・全 volume・O(N log N)。回転/スケールは別途(PCA / log-polar)。

引数: `a, b` は同形の 3-D 配列(違えば ValueError)。float32 に落として FFT する。

返り値: int の tuple `(dz, dy, dx)、各軸 (−N/2, N/2]` に折り返し済み。意味は

`np.roll(b, (dz,dy,dx), axis=(0,1,2)) ≈ a`(b をこれだけ動かすと a に重なる)。

• 循環相関なので、はみ出した部分は反対側から回り込む(窓掛けはしない)。シフトが volume の

半分を超えると符号が反転して見える。

• 位相のみ(`R/|R|`)なので振幅・コントラスト差に不変だが、ノイズが白色化されてピークが

埋もれることがある。全 0 の volume は 0 になり index 0 を返す。

• 整数精度。サブボクセルは `refine_translation_lk / refine_peak_newton` へ。

• 回転・スケールがあると効かない(`match_logpolar_z` → 回転補正 → 本 op の順)。

References (sample data, literature)

• Sample-data catalog (download URLs / licences) — 2-D uses skimage.data (BSD/public domain) plus synthetic images; 3-D lists download URLs for real data sources (Stanford, PDS, …).

• Operator provenance and references — the sources of the research/methods this op family came from.

Runnable examples (verified samples that actually call this op)

shape_desc_posepy -3.11 examples_3d/shape_desc_pose.py

Ops the type connects to (they accept shift as input)

fuse_to_voxel

Same category (match_pose)

match_pca · moment_axes · match_logpolar_z


*Provenance: match3d.py — 3D operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*

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