log_zero_crossings — 3D `edges` op

• 数据种类:voxel → voxel

• 调用: import fullseye as fs; fs.ledger.log_zero_crossings(vol, sigma: 'float' = 1.5, rel_thresh: 'float' = 0.001) -> 'np.ndarray'(要直接调用实现,import edges3d; edges3d.log_zero_crossings(vol, sigma: 'float' = 1.5, rel_thresh: 'float' = 0.001) -> 'np.ndarray';从台账取用则 ops3d.get("log_zero_crossings"))

用法

Laplacian-of-Gaussian 的零交叉边缘。

Parameters

vol : (D,H,W) array

灰度 voxel。

sigma : float

LoG 的标准差(> 0)。

rel_thresh : float

抑制平坦区数值噪声造成的伪交叉的相对阈值。仅当交叉对的 LoG 差 |a-b|

超过 rel_thresh * max|L| 时才采用。

Returns

edge_mask : (D,H,W) bool

在符号变化的相邻对中,将 |LoG| 较小的一侧(更接近零的一侧)置位的 mask。

Notes

在各轴方向上检测两种过零点:

(1) 格间交叉 —— 相邻对 (a,b) 的符号积 a*b < 0(两侧都非零且符号反转)。为了抑制厚度,分配给对中 |LoG| 较小一侧的 1 个 voxel。

(2) 格点对齐交叉 —— 中央 voxel 严格为零(L==0)且两邻异号(符号在该 voxel 处穿过零)的情况。在 SDF 风格的格点对齐面上,LoG 恰好在格点上为 0 的真实边缘,a*b 总是 0,(1) 会漏检,因此把该 voxel 作为交叉置位。常数零区域两邻也为 0,Lm*Lp==0,因而被排除(不会误检)。

两种情况下,平坦区的数值噪声都由相对阈值(rel_thresh * max|L|)抑制。

参考(示例数据、文献)

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

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

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

• edges_3d — py -3.11 examples_3d/edges_3d.py

类型可衔接的下一个算子(可接受 `voxel` 作为输入)

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

同类别(`edges`)

gradient3d · canny3d · link_edges · edge_points

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

log_zero_crossings — 3D edges op

用法