vol_gradient_magnitude — 3D feature op

• Data kinds: voxel → voxel

• Call: import fullseye as fs; fs.ledger.vol_gradient_magnitude(vol) (implementation को सीधे call करना हो तो import volops; volops.vol_gradient_magnitude(vol), registry से लेना हो तो ops3d.get("vol_gradient_magnitude"))

उपयोग

3-D Sobel gradient magnitude sqrt(gz**2 + gy**2 + gx**2)।

हर g* एक axis के साथ scipy.ndimage.sobel derivative है। Response intensity की boundaries पर केंद्रित होता है (step edge interface पर चमकता है और flat क्षेत्रों में लगभग 0 रहता है)। (D, H, W) float64 volume लौटाता है।

गणना: scipy.ndimage.sobel को axis 0 (z), 1 (y), 2 (x) के क्रम में लगाया जाता है, और sqrt(gz**2 + gy**2 + gx**2) लौटाया जाता है। Sobel derivative [-1,0,1] और smoothing [1,2,1] का separable kernel है (2 orthogonal axes पर हर एक में 4 गुना), इसलिए axis के लंबवत unit step (0→1) का response 1 नहीं बल्कि 16 होता है (interface के दोनों ओर के 2 voxel पर 16.0, flat हिस्से में 0; scipy के default weights, boundary mode='reflect')। मान normalise नहीं होते ([0, 1] में नहीं आते)। spacing नहीं लिया जाता, अंतर voxel इकाई में है (anisotropic voxel में भी axis के हिसाब से correction नहीं)।

जाँच (ValueError): 3-D न हो / NaN या Inf हो / voxel की संख्या MAX_VOXELS (1 << 27) से ज़्यादा हो।

कहाँ इस्तेमाल करें: vol_watershed का landscape input, vol_local_maxima से interface की चोटियाँ चुनना, vol_stretch से [0, 1] में normalise करके दिखाना।

संदर्भ (sample data, साहित्य)

• Sample data catalog (DL URL / license) — 2-D के लिए skimage.data (BSD/public) + synthetic, 3-D के लिए real data sources (Stanford/PDS आदि) के DL URL।

• Operators का स्रोत और संदर्भ — इस op family के मूल research/methods के स्रोत।

चलाए जा सकने वाले उदाहरण (verified samples जो सच में यह op call करते हैं)

• vessel_metrology — py -3.11 examples_3d/vessel_metrology.py

वे अगले ops जिनसे type जुड़ता है (voxel को input के रूप में लेते हैं)

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

इसी category के (feature)

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


*Provenance: volops.py — 3D operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*

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