vol_fft_highpass — 3D frequency op

資料種類:voxelvoxel

呼叫: 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` 的精確補(兩者相加還原)。

> 以下的詳細說明為原文 —— 摘要與標題已翻譯。

to the input to float precision, proven in tests). Removes the DC level and

slow drift, keeps edges/texture. Output is signed (mean ~ 0).

伝達関数は `1 - exp(-|f|^2 / (2 cutoff^2))。DC(f = 0`)は係数 0 なので平均

輝度は完全に落ち、`|f| = cutoff1 - 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 の性質上、面どうしの輝度差は wrap で漏れる。

參考(範例資料・文獻)

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

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

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

deconv_fft_restorepy -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.