tps_fit — 3D deform op

• 資料種類:points × points → deformation

• 呼叫: import fullseye as fs; fs.ledger.tps_fit(src_ctrl, dst_ctrl, lam=0.0)(要直接呼叫實作,import deform3d; deform3d.tps_fit(src_ctrl, dst_ctrl, lam=0.0);從台帳取用則 ops3d.get("tps_fit"))

用法

由控制點對應擬合 3D Thin-Plate-Spline。

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

移動側の制御点 `src_ctrl を固定側 dst_ctrl` へ写す TPS 係数を、鞍点系

を最小二乗(`numpy.linalg.lstsq`)で解いて求める。λ=0 なら制御点上で厳密に

内挿(`tps_warp(model, src_ctrl) == dst_ctrl`)、λ>0 で平滑化する。

引数:

src_ctrl: (K,3) 制御点(変形の始点、TPS のカーネル中心 p_i)。

dst_ctrl: (K,3) 対応する目標点(変形の終点 v_i)。

lam: 正則化係数 λ≥0。カーネル行列 K の対角へ λ を加える。大きいほど

変形は滑らか(制御点への当てはめは緩む)。

返り値:

model: dict。キーは

"ctrl" (K,3) カーネル中心 p_i、

"w" (K,3) 非線形(曲げ)係数、

"a" (4,3) アフィン係数([平行移動; 線形部]、a[0]=c, a[1:4]=Aᵀ)、

"lam" 使用した λ。

例外:

ValueError: 形状不一致、制御点数不足、非有限値、または制御点数が

`TPS_MAX_CTRL`(10,000)超(密な (K+4)² 系は O(K³) のため。

対応点を間引いてから渡す — TPS は疎な制御点で滑らかな変形を

表現するのが本来の使い方)。

參考(範例資料・文獻)

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

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

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

• nonrigid_deform — py -3.11 examples_3d/nonrigid_deform.py

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

fuse_to_voxel · tps_warp

同類別(deform)

tps_warp · register_nonrigid · register_cpd_rigid

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

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