tps_set_ctrl_pts.m

薄板样条

function [PP,kernel] = tps_set_ctrl_pts(ctrl_pts)


[n,d] = size(ctrl_pts);
K = tps_compute_kernel(ctrl_pts, ctrl_pts);
Pn = [ones(n,1) ctrl_pts];
PP = null(Pn');  % or use qr(Pn)
kernel = PP'*K*PP;