tps_compute_param.m

薄板样条

function [param] = tps_compute_param(PP,kernel,U,Pm,Q1,Q2,R,lambda,target)


TB = U*PP;
QQ = Q2*Q2';
A = inv(TB'*QQ*TB + lambda*kernel)*TB'*QQ;
tps = A*target;
affine = inv(R)*Q1'*(target-TB*tps);
param = [affine; tps];