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"))

用法

兩個 voxel → 正規化互相關的 score 體。該型別的唯一入口。

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 以 np.roll 平移 s0 得到的, 峰值會嚴格地立在 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)/ 非有限。

參考(範例資料、文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者、年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• 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

型別可銜接的下一個運算子(可接受 score 作為輸入)

score_to_position · score_to_image2d

同類別(score)

score_to_position · score_to_image2d

*Provenance: reprconv.py — REPRCONV 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

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