fullseye

correlation_score — REPRCONV score op

• データ種: voxel × voxel → score
• 呼び出し: import fullseye as fs; fs.ledger.correlation_score(voxel_a, voxel_b) (実装を直接呼ぶなら import reprconv; reprconv.correlation_score(voxel_a, voxel_b)、台帳から引くなら opsreprconv.get("correlation_score"))

使い方

2 つの voxel → 正規化相互相関の score volume。この型の唯一の入口。

score は refine_peak_newton が食う型だが、台帳のどの op も score を産まなかった(実測。生成器の種を置いてようやく到達していた)。 ここでは FFT による循環相互相関を返す:

score[s] = sum_x (a[x] - mean_a) * (b[x + s] - mean_b) / (N * std_a * std_b)

したがって b が a を s0 だけ np.roll したものなら、 ピークは厳密に s0 に立つ(閉形式の真値。テストがこれを使う)。 循環相関なので端は巻き込む —— 打ち切り相関ではないことを明記しておく。

Args: voxel_a: (D, H, W)。 voxel_b: (D, H, W)、voxel_a と同形。 Returns: (D, H, W) float64、値域は概ね [-1, 1] (完全一致で 1.0)。 Raises: ValueError: 3-D でない / 形が違う / 定数体積(標準偏差 0)/ 非有限。

参考(サンプルデータ・文献)

• サンプルデータ カタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
• 演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。
• アルゴリズムの正典(著者・年)と用途は上記ファミリ使い方ガイドに記載。

実行できる例(この op を実際に呼ぶ検証済みサンプル)

• poc_leak_localization — py -3.11 examples/poc_leak_localization.py
• representation_conversion — py -3.11 examples/representation_conversion.py
• representation_roundtrip — py -3.11 examples/representation_roundtrip.py

型が繋がる次の op(score を入力に取れる)

score_to_position · score_to_image2d

同カテゴリ(score)

score_to_position · score_to_image2d

---

Provenance: reprconv.py — REPRCONV operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。

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

This site is open source. Improve this page.