vol_fft_highpass — 3D frequency op

• 資料種類:voxel → voxel

• 呼叫: import fullseye as fs; fs.ledger.vol_fft_highpass(vol, cutoff, spacing=None)(要直接呼叫實作,import volfreq; volfreq.vol_fft_highpass(vol, cutoff, spacing=None);從台帳取用則 ops3d.get("vol_fft_highpass"))

用法

高斯高通 —— 1 - lowpass 的精確補(兩者相加還原)。

該和在浮點精度內成立,已在測試中證明。去除 DC 分量和緩慢漂移,保留邊緣/紋理。輸出帶符號(平均 ~ 0)。

傳遞函數為 1 - exp(-|f|^2 / (2 cutoff^2))。DC(f = 0)係數為 0,因此平均輝度被完全去除,在 |f| = cutoff 處為 1 - 0.61 ≈ 0.39,頻率越高越接近 1。回傳值是與輸入相同 (D, H, W) 的 float64,且 帶符號(含負值)。若要傳給以 [0, 1] 為前提的後級(顯示、vol_window_level 等),請用 vol_stretch 正規化,或以 vol + highpass 的形式使用。

參數: cutoff 為正有限值。spacing=None 時為 cycles/voxel,傳入 spacing ((sz, sy, sx) 或 volio.VolumeMeta)時為 cycles/mm。

驗證(ValueError): 與 vol_fft_lowpass 相同(不是 3-D / NaN、Inf / 超過 MAX_VOXELS / cutoff、spacing 不合法)。

使用時機: 去除照明不均、厚度漂移、背景梯度(等同於減去 vol_fft_lowpass 的結果)。由於 FFT 將 volume 視為週期性的,各面之間的輝度差會經 wrap 洩漏。

參考(範例資料、文獻)

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

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

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

• deconv_fft_restore — py -3.11 examples_3d/deconv_fft_restore.py

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

voxel_to_mips · voxel_to_mesh · signed_distance_field · to_points · sobel3d · hessian3d · curvature_maps · edt_jfa

同類別(frequency)

vol_fft_lowpass · vol_fft_bandpass

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

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