vol_local_std — 3D feature op

数据种类:voxelvoxel

调用: import fullseye as fs; fs.ledger.vol_local_std(vol, size=5)(要直接调用实现,import volops; volops.vol_local_std(vol, size=5);从台账取用则 ops3d.get("vol_local_std"))

用法

立方窗口内的局部标准差(无偏)。

> 以下的详细说明为原文 —— 摘要与标题已翻译。

The 3-D counterpart of the registry op `local_std` (HALCON's

`deviation_image`). The variance is taken **after subtracting the global

mean** — `E[x^2] - E[x]^2` loses every significant digit on a volume whose

values sit far from zero, and the leftover shows up as a fake texture on a

uniform block. It is then unbiased twice: `n/(n-1)` on the variance and

`c4(n) on the square root, so the estimate of sigma` itself (not of

`sigma^2`) is unbiased.

*size* is the edge of the cubic window in voxels (odd, >= 3), so

`n = size**3` voxels enter each estimate. **The relative standard error of

a single voxel's estimate is `1/sqrt(2(n-1))`** — 3.1% for a 5x5x5 window,

1.1% for 9x9x9. Quote that figure whenever a noise or roughness number is

read off this volume.

Returns a `(D, H, W)` float64 volume in the input's own units (it is *not*

normalised, so values from different scans compare directly).

参考(示例数据・文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

可运行的示例(实际调用该算子并已验证的样例)

ct_porosity_and_fibre_morphometrypy -3.11 examples_3d/ct_porosity_and_fibre_morphometry.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_thickness · vol_orientation_coherence · vol_euler_number


*Provenance: volops.py — 3D 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

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