frequency op• 데이터 종류: voxel → voxel
• 호출: import fullseye as fs; fs.ledger.vol_fft_lowpass(vol, cutoff, spacing=None)(구현을 직접 호출하려면 import volfreq; volfreq.vol_fft_lowpass(vol, cutoff, spacing=None), 원장에서 가져오려면 ops3d.get("vol_fft_lowpass"))
가우스 저역통과: `1/cutoff` 보다 거친 구조를 남깁니다(복셀 단위 또는 실치수).
> 아래 상세 설명은 원문입니다 —— 요약과 제목은 번역되어 있습니다.
mm with *spacing*), attenuates finer detail smoothly. Transfer
`exp(-f^2 / (2 cutoff^2))` — the DC level (mean intensity) passes
unchanged. Typical use: extract the illumination/thickness drift.
手順: `np.fft.rfftn で実 FFT → 各周波数 bin の大きさ |f| = sqrt(fz^2 + fy^2 + fx^2)`
に伝達関数 `exp(-|f|^2 / (2 cutoff^2)) を掛ける → irfftn` で戻す。入力と同じ
`(D, H, W)` の float64 を返す(spacing を渡しても shape は変わらない)。
引数:
• `cutoff: 正の有限値。spacing=None` なら cycles/voxel(Nyquist = 0.5)、
`spacing=(sz, sy, sx)(mm、または spacing_mm を持つ volio.VolumeMeta`)を
渡すと cycles/mm。異方 voxel では軸ごとに `fftfreq(n, d=spacing)` で物理
周波数に直すので、同じ物理構造が同じ扱いになる。
• 減衰は Gaussian で、`|f| = cutoff で exp(-1/2) ≈ 0.61`(brick-wall では
ない。リンギングは出ない)。
検証(`ValueError): 3-D でない / NaN・Inf / voxel 数が MAX_VOXELS(1 << 27`)
超 / `cutoff が非正・非有限・2 乗がアンダーフローするほど小さい / spacing` が
長さ 3 でない・非正。
注意: FFT は volume を周期的とみなす。向かい合う面の輝度差が大きいと wrap を
またいで漏れる(窓掛けは黙ってしない)。`vol_fft_highpass` は厳密にこの補集合
`1 - lowpass` で、両者の和は入力に一致する。
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL 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_highpass · vol_fft_bandpass
*Provenance: volfreq.py — 3D 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.