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.