vol_gradient_magnitude — 3D feature op

• 資料種類:voxel → voxel

• 呼叫: 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)。

每個 g* 是沿一個軸的 scipy.ndimage.sobel 導數。響應集中在強度邊界上(階躍邊緣在介面處亮起,平坦區域約為 0)。回傳 (D, H, W) float64 的 volume。

計算: 依 axis 0 (z)、1 (y)、2 (x) 的順序施加 scipy.ndimage.sobel, 回傳 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_metrology — py -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 · edt_jfa_vector · vol_frangi · vol_local_std · vol_local_thickness


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

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