vol_gradient_magnitude — 3D feature op

資料種類:voxelvoxel

呼叫: import fullseye as fs; fs.ledger.vol_gradient_magnitude(vol)(要直接呼叫實作,import volops; volops.vol_gradient_magnitude(vol);從台帳取用則 ops3d.get("vol_gradient_magnitude"))

用法

3-D Sobel 梯度幅值 `sqrt(gz**2 + gy**2 + gx**2)`。

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

Each `g* is a scipy.ndimage.sobel` derivative along one axis. The

response localises at intensity boundaries (a step edge lights up on the

interface and is ~0 in flat regions). Returns a `(D, H, W)` float64 volume.

計算: `scipy.ndimage.sobel` を axis 0 (z), 1 (y), 2 (x) の順に掛け、

`sqrt(gz**2 + gy**2 + gx**2)` を返す。Sobel は微分 [-1,0,1] と平滑 [1,2,1] の

分離カーネル(直交 2 軸で各 4 倍)なので、軸に垂直な単位ステップ(0→1)に対する

応答は 1 ではなく 16(界面をはさむ 2 voxel で 16.0、平坦部は 0。scipy の既定

重み・境界 `mode='reflect')。値は正規化しない([0, 1]` には収まらない)。

spacing は受けず、voxel 単位の差分(異方 voxel でも軸ごとに補正しない)。

検証(`ValueError): 3-D でない / NaN・Inf を含む / voxel 数が MAX_VOXELS`

(`1 << 27`)を超える。

使いどころ: `vol_watershed の地形(landscape)入力、vol_local_maxima` で

界面の峰を拾う、`vol_stretch[0, 1]` に正規化して表示。

參考(範例資料・文獻)

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

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

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

vessel_metrologypy -3.11 examples_3d/vessel_metrology.py

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

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

同類別(feature)

sobel3d · hessian3d · curvature_maps · edt_jfa · vol_frangi · vol_local_std · vol_local_thickness · vol_orientation_coherence


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

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