3D Point set warping by thin-plate/rbf function

3D Point set warping by thin-plate/rbf function

Radial basis function/Thin-plate spline 3D point set warping.
p3do = rbfwarp3d( p3d, ps, pd, varargin )
input:
p3d: 3d point set
ps: 3d source landmark [n*3]
pd: 3d destin landmark [n*3]
method:
'gau',r - for Gaussian function ko = exp(-|pi-pj|/r.^2);
'thin' - for Thin-plate function ko = (|pi-pj|^2) * log(|pi-pj|^2)
output:
p3do: output point set

Bookstein, F. L.
"Principal Warps: Thin Plate Splines and the Decomposition of Deformations."
IEEE Trans. Pattern Anal. Mach. Intell. 11, 567-585, 1989.