feature op• 数据种类:voxel → curvature
• 调用: import fullseye as fs; fs.ledger.curvature_maps(vol, device='cpu', mc=0.000625)(要直接调用实现,import match3d; match3d.curvature_maps(vol, device='cpu', mc=0.000625);从台账取用则 ops3d.get("curvature_maps"))
• GPU:该算子有 GPU 路径(device="cuda")
level-set 的主曲率 → 形状指数 S(Koenderink)与弯曲度。闭式(Kindlmann 2003)。
2D 轮廓曲率(1 个标量)的 线→面提升: 曲面有 2 个主曲率 κ1,κ2。mean = (κ1+κ2)/2 = (|g|²trH − gᵀHg)/|g|³,Gauss K = κ1κ2 = gᵀadj(H)g/|g|⁴(g=∇, H=Hessian)。S=(2/π)atan2(κ1+κ2, κ1−κ2) ∈[-1,1] 是 对强度、旋转不变的局部曲面类型(cup−1/rut/saddle0/ ridge+.5/cap+1)。按外向法线约定,亮的凸 blob=cap(+1)。返回值 (S, curvedness, mask, |g|),全为 torch。
单位制: 在此把 sobel3d 可分离 conv 的增益 32(deriv[-1,0,1]×平滑[1,2,1]²)除回去,因此 κ1,κ2/curvedness 为 真实的 1/voxel 单位(半径 R 的球壳上 curvedness=1/R),|g| 和 mask 阈值 mc 为 每 voxel 的真实梯度单位。旧版(至 2026-08-29)忘了除回, curvedness 为 1/32 倍、mc 为原始 sobel3d 单位(shape index S 是比值,不受影响)。若使用过旧 mc 值,请乘以 1/32 后传入(默认值 0.02→6.25e-4 也已等价换算)。
• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。
• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。
• diff_features — py -3.11 examples_3d/diff_features.py
curvature 作为输入)feature)sobel3d · hessian3d · 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.