snake_evolve — SEGMENTATION contour op

• 資料種類:image2d × pairs → table

• 呼叫: import fullseye as fs; fs.ledger.snake_evolve(image, init_points, *, alpha: 'float' = 0.1, beta: 'float' = 0.1, gamma: 'float' = 1.0, external: 'str' = 'edge', w_line: 'float' = 0.0, w_edge: 'float' = 1.0, sigma: 'float' = 1.0, kappa: 'float' = 1.0, gvf=None, gvf_mu: 'float' = 0.2, max_px_move: 'Optional[float]' = None, n_iter: 'int' = 300, spline_order: 'int' = 3, tol: 'float' = 0.0, resample_every: 'int' = 0, record_every: 'int' = 0) -> 'Dict[str, object]'(要直接呼叫實作,import segcontour; segcontour.snake_evolve(image, init_points, *, alpha: 'float' = 0.1, beta: 'float' = 0.1, gamma: 'float' = 1.0, external: 'str' = 'edge', w_line: 'float' = 0.0, w_edge: 'float' = 1.0, sigma: 'float' = 1.0, kappa: 'float' = 1.0, gvf=None, gvf_mu: 'float' = 0.2, max_px_move: 'Optional[float]' = None, n_iter: 'int' = 300, spline_order: 'int' = 3, tol: 'float' = 0.0, resample_every: 'int' = 0, record_every: 'int' = 0) -> 'Dict[str, object]';從台帳取用則 opssegmentation.get("snake_evolve"))

用法

離散 snake(Kass–Witkin–Terzopoulos 1988)的半隱式更新 x_t = (A + γI)^(-1)(γ x_(t-1) − f_x)(式 (19)(20))(snake_evolve)。

> 以下的詳細說明為原文 —— 摘要與標題已翻譯。

閉じた snake(周期境界)の点 `init_points`(N, 2)= [row, col]、N ≥ 5。内部エネルギー

E_int = ½ Σ (α |v_(i+1) − v_i|² + β |v_(i+1) − 2v_i + v_(i−1)|²) = ½ xᵀ A x(陰的)。外力(陽的)は `external` で選ぶ:

`"edge"` = ポテンシャル P = w_line · G_σ*I + w_edge · |∇(G_σ*I)|²、E_ext = −κ Σ P(v_i)(Kass の E_line と E_edge)、

`"potential" = image` をそのまま P として使う(第 2 実装と同じ P を渡すため)、

`"gvf" = 力そのものが −κ (v, u)(GVF、gvf に gvf_field の返りを渡すか、gvf_mu` で内部で作る)、

`"none" = 外力 0。P は spline_order` 次の補間スプライン(値と微分が整合)で点の位置へ。

`max_px_move` を与えると 1 反復の移動を max_px_move · tanh(Δ) に制限(skimage と同じ。近接勾配の性質は失う)。

`resample_every` > 0 なら k 反復ごとに弧長で等間隔に打ち直す(凹部へ入るとき点が疎になるのを防ぐ。エネルギーの

単調性はその時点で途切れる)。点は画像の範囲 [0, H−1] × [0, W−1] に切り詰める。

返り値: `points(N, 2)、mask(折れ線の内側)、energy / energy_internal / energy_external`(反復ごと、

先頭 = 初期。GVF では外力のエネルギーが定義されないので total と external は NaN)、`n_increase`(エネルギーが

増えた反復の数、GVF なら None)、`max_increase、lipschitz`(κP のヘッセ行列のスペクトルノルムの格子上の最大 =

外力の勾配の Lipschitz 定数の見積り)、`gamma_ge_lipschitz`(γ ≥ L なら降下補題で単調非増加が保証される)、

`radius_factor(外力 0 の円の 1 反復の縮み γ/(γ + λ_1) の閉形式)、history`(record_every ごとの点)、

`n_iter、converged`。

詳細使用指南

• halcon_segmentation 族使用指南

背景知識指南(該運算子背後的物理與規範)

• active_contours_and_level_sets — 変分・動的輪郭とレベルセット — どの輪郭がどこで止まり、何を保証するか

參考(範例資料・文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者・年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• poc_active_contours — py -3.11 examples/poc_active_contours.py

型別可銜接的下一個運算子(可接受 table 作為輸入)

class_ndim_norm

同類別(contour)

gvf_field · chan_vese_energy · chan_vese_evolve · morph_chan_vese · morph_geodesic_ac · edge_stop_g · level_set_reinit · drle_evolve


*Provenance: segcontour.py — SEGMENTATION 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

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