hessian3d — 3D feature op

• 数据种类:voxel → hessian

• 调用: import fullseye as fs; fs.ledger.hessian3d(vol, device='cpu')(要直接调用实现,import match3d; match3d.hessian3d(vol, device='cpu');从台账取用则 ops3d.get("hessian3d"))

• GPU:该算子有 GPU 路径(device="cuda")

用法

3D Hessian 的 6 个独立分量 (fzz,fyy,fxx,fzy,fzx,fyx)。可分离 conv3d(二阶/一阶×平滑)。

核: 2 阶 [1,−2,1] (增益 1)、1 阶 [−0.5,0,0.5] (增益 1)、其余轴为 [1,2,1]/4 的平滑(增益 1)。对角分量为「该轴的 2 阶 × 其他 2 轴的平滑」,交叉分量为「2 个轴的 1 阶 × 其余轴的平滑」。单位为 1/voxel² 的真实值(与 sobel3d 的 32 倍增益不同)。边缘 replicate。返回值为 6 个 torch tensor 组成的 list,各为 (D,H,W) float32,位于 device 上,顺序为(zz, yy, xx, zy, zx, yx)(轴 0=z, 1=y, 2=x)。需要 numpy 时用 .cpu().numpy()。输入为 numpy 等价物(转换为 float32)。

用途: curvature_maps 的主曲率(与 sobel3d 配合使用)、blob/管状结构的检测。

参考(示例数据、文献)

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

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

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

• diff_features — py -3.11 examples_3d/diff_features.py

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

fuse_to_voxel

同类别(feature)

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


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

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