tps_set_ctrl_pts.m

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

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;