tps_warp.m

薄板样条

function [target, bending] = tps_warp(PP,kernel,U,Pm,param)


% param = [affine; tps]
[m,d] = size(Pm);
basis = [Pm U*PP];
target = basis*param;
tps = param(d+1:end,:);
bending = trace(tps'*kernel*tps);