robust_spline_profile — ROUGHNESS `prepare` op

• 数据种类:signal → signal

• 调用: import fullseye as fs; fs.ledger.robust_spline_profile(z, dx, lambda_c, c=4.4, robust=True, max_iter=None, tol=1e-06)(要直接调用实现,import roughness; roughness.robust_spline_profile(z, dx, lambda_c, c=4.4, robust=True, max_iter=None, tol=1e-06);从台账取用则 opsroughness.get("robust_spline_profile"))

• 经台账调用的返回值: fullseye.ledger.robust_spline_profile(...) 只返回声明的 out 类型 signal 的值(底层函数还会返回辅助信息)。需要被丢弃的部分时,使用 fullseye.ledger.robust_spline_profile.raw(...),或直接调用 roughness.robust_spline_profile。

• 底层返回: (roughness, mean_line, info)

用法

稳健样条滤波器(截面):不被离群点拖动的平均线与粗糙度。

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

平均線 w を Σ δ_i (z_i - w_i)² + β Σ (Δ² w)_i² の最小化で求める(2 階差分の罰則つき最小二乗 = 離散の三次平滑化スプライン、Whittaker 1923 / Eilers 2003 / Krystek, Measurement 18 (1996) 9-15)。重み δ は biweight で反復して更新する。ガウスと違って 端に自然境界条件が付く ので、端を削らなくても平均線が端で外へ振れにくい。

β = 1 / (16 sin⁴(π dx/λc)) と置くと、線形(robust=False)の内部の伝達関数が H(λ) = 1 / (1 + (sin(π dx/λ) / sin(π dx/λc))⁴) になり、λ = λc でちょうど 50 % (閉形式)。dx/λ が小さい極限では 1/(1+(λc/λ)⁴)(張力 0 の三次スプラインフィルタの形)。

戻り値は (roughness, mean_line, info)。引数と fail-closed は robust_gaussian_profile と同じ(order / truncate は無い)。重みを持つデータが 2 点未満なら拒否する。

详细使用指南

• surface_roughness 族使用指南

参考(示例数据、文献)

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

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

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

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

• poc_robust_profile_filter — py -3.11 examples/poc_robust_profile_filter.py

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

robust_gaussian_profile · profile_params

同类别(`prepare`)

surface_form_remove · surface_filter · robust_gaussian_profile · robust_gaussian_surface

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

robust_spline_profile — ROUGHNESS prepare op

用法