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.