tps_warp.m

来自「薄板样条」· M 代码 · 共 10 行

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);