vol_fft_bandpass — 3D frequency op

• 数据种类:voxel → voxel

• 调用: import fullseye as fs; fs.ledger.vol_fft_bandpass(vol, low, high, spacing=None)(要直接调用实现,import volfreq; volfreq.vol_fft_bandpass(vol, low, high, spacing=None);从台账取用则 ops3d.get("vol_fft_bandpass"))

用法

高斯带通 lowpass(high) - lowpass(low):保留两个尺度之间的结构。

要求 low < high,二者单位均为 cycles/voxel,或在给定 *spacing* 时为 cycles/mm。典型用途:分离出一个纹理尺度,或在减去之前取出周期性伪影频带。

传递函数为 exp(-|f|^2 / (2 high^2)) - exp(-|f|^2 / (2 low^2))。在 DC 处为 1 - 1 = 0,因此平均亮度被去除,只留下峰值位于 low 与 high 之间的高斯差分(DoG 型)频带。峰高达不到 1(因为是 2 个 Gaussian 之差, low 与 high 越接近,通过量越小)。返回值是与输入相同 (D, H, W) 的 float64,有符号。

参数:必须 low < high(二者均为正的有限值)。单位在 spacing=None 时为 cycles/voxel,指定 spacing 时为 cycles/mm。若按周期(voxel 或 mm)考虑, 1/high 是最小结构尺寸,1/low 是最大结构尺寸。

验证(ValueError):非 3-D / NaN、Inf / 超过 MAX_VOXELS / low、high 非正、非有限、2 次方下溢 / low >= high / spacing 无效。

用途:提取 1 个纹理尺度;取出周期性条纹、环形伪影频带并从输入中减去。与 vol_fft_lowpass(vol, high) - vol_fft_lowpass(vol, low) 等价。

参考(示例数据、文献)

• 示例数据目录(下载 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_highpass

*Provenance: volfreq.py — 3D 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

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