freeform op• 数据种类:image2d × image2d × image2d × bspline_surface → measurement
• 调用: import fullseye as fs; fs.ledger.surface_residual(x, y, z, tck)(要直接调用实现,import bspline_surf; bspline_surf.surface_residual(x, y, z, tck);从台账取用则 ops3d.get("surface_residual"))
• 经台账调用的返回值: fullseye.ledger.surface_residual(...) 只返回**声明的 out 类型 measurement 的值**(底层函数还会返回辅助信息)。需要被丢弃的部分时,使用 fullseye.ledger.surface_residual.raw(...),或直接调用 bspline_surf.surface_residual。
• 底层返回: {"rms","max","pv"} → pv float
返回散布数据与曲面 tck 的残差统计(形状误差=相对拟合的偏离)。
> 以下的详细说明为原文 —— 摘要与标题已翻译。
各サンプル位置で曲面高さを評価し、観測 z との差の RMS / 最大絶対値 / PV
(peak-to-valley = 最大 - 最小)を計測する。検査ではこれが自由曲面からの
ずれ量(打痕・うねり・欠肉)の定量指標になる。
Returns
-------
dict
`{"rms": float, "max": float, "pv": float}`。max は最大絶対残差、
pv は符号付き残差の最大 - 最小(片側だけの凸/凹も捉える)。
補足:
• 残差の符号は `resid = z - zhat(観測 − 曲面)。正 = 曲面より高い(盛り上がり)、負 = 低い(欠肉)。max は絶対値なので向きは分からず、向きが要るなら eval_bspline_surface で zhat` を取り自分で差を取る。
• x, y, z は任意 shape を受け、内部で `eval_bspline_surface(tck, x, y, grid=False) を通す(tck の種類検査もそこで行われ、曲線 tck や多項式 dict は ValueError`)。z と x の要素数が違うと numpy の形状エラーになる(専用の検査は無く、z が 1 要素だと黙って broadcast される)。
• 単位は z と同じ。`rms は全点の二乗平均平方根なので局所的な打痕は平均で薄まる。局所欠陥は max か pv` で見る。
• 典型: `fit_bspline_surface` → 本 op。フィットに使った点で評価すると smooth が小さいほど残差は 0 に近づく(過適合の指標にもなる)。
• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。
• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。
• bspline_freeform — py -3.11 examples_3d/bspline_freeform.py
measurement 作为输入)vol_gaussian_psf · fuse_to_voxel · fresnel_reflectance · snell_angle
freeform)fit_bspline_surface · eval_bspline_surface · fit_bspline_curve · eval_bspline_curve
*Provenance: bspline_surf.py — 3D 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.