Radially symmetric thin plate splines interpolating a circular contour map
This work provides theoretical foundations for radial thin plate splines on unbounded domains, which is incremental for researchers in approximation theory and scattered data interpolation.
The paper derives two types of radially symmetric thin plate spline profiles that minimize the Beppo Levi energy over the full semi-axis, achieving an L^p-approximation order of 3/2+1/p. Numerical examples and a novel basis function representation are provided.
Profiles of radially symmetric thin plate spline surfaces minimizing the Beppo Levi energy over a compact annulus $R_{1}\leq r\leq R_{2}$ have been studied by Rabut via reproducing kernel methods. Motivated by our recent construction of Beppo Levi polyspline surfaces, we focus here on minimizing the radial energy over the full semi-axis $0<r<\infty$. Using a $L$-spline approach, we find two types of minimizing profiles: one is the limit of Rabut's solution as $R_{1}\rightarrow0$ and $R_{2}\rightarrow\infty$ (identified as a `non-singular' $L$-spline), the other has a second-derivative singularity and matches an extra data value at $0$. For both profiles and $p\in\left[ 2,\infty\right] $, we establish the $L^{p}$-approximation order $3/2+1/p$ in the radial energy space. We also include numerical examples and obtain a novel representation of the minimizers in terms of dilates of a basis function.