hessian_vesselness — REPRCONV derivative op

• 数据种类:hessian → voxel

• 调用: import fullseye as fs; fs.ledger.hessian_vesselness(hessian, alpha=0.5, beta=0.5, c=None, black_ridges=False)(要直接调用实现,import reprconv; reprconv.hessian_vesselness(hessian, alpha=0.5, beta=0.5, c=None, black_ridges=False);从台账取用则 opsreprconv.get("hessian_vesselness"))

• GPU:无(仅 CPU 参考实现)

用法

3-D Hessian → Frangi 管状度(vesselness)voxel,值域 [0, 1]。

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

固有値を |λ1| ≤ |λ2| ≤ |λ3| に並べ、Ra = |λ2|/|λ3|(板と管を分ける)、Rb = |λ1|/sqrt(|λ2 λ3|)(塊を除く)、S = sqrt(λ1²+λ2²+λ3²)(背景の雑音を除く)から

V = (1 − exp(−Ra²/2α²)) · exp(−Rb²/2β²) · (1 − exp(−S²/2c²))

を返す(Frangi et al. 1998, MICCAI)。明るい管(black_ridges=False)では λ2 > 0 か λ3 > 0 の voxel を 0 にする(暗い管なら符号を逆に)。c の既定は Frangi の勧めどおり S の最大値の半分。

vol_frangi は体積から Hessian を自分で作る(尺度ごとに)のに対し、こちらは hessian 型を受ける —— Hessian の出どころ(hessian3d / 外部の微分)を選べ、同じ Hessian から hessian_determinant など複数の量を引き出せる。1 尺度だけ(多尺度の最大は呼ぶ側で取る)。

Args:

hessian: hessian3d の 6 成分 (zz, yy, xx, zy, zx, yx)。

alpha, beta: Frangi の感度(既定 0.5 / 0.5)。> 0。

c: 構造の強さの感度。None で max(S)/2。> 0。全面 0 の Hessian では 0 を返す。

black_ridges: True で暗い管(背景より暗い)を拾う。

Returns:

(D, H, W) float64、[0, 1]。

Raises:

ValueError: 6 成分でない / 非有限 / alpha・beta・c が正でない。

参考(示例数据、文献)

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

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

• 算法的正典(作者、年份)与用途见上面的族使用指南。

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

• type_connectivity_tour — py -3.11 examples/type_connectivity_tour.py

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

correlation_score

同类别(derivative)

gradient_magnitude · gradient_coherence · hessian_determinant


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

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