--- title: Local RBF Interpolation Overview url: https://www.emergentmind.com/topics/local-radial-basis-function-interpolation type: topic --- # Local RBF Interpolation Overview Searching arXiv for recent and foundational work on local radial basis function interpolation. Local radial basis function interpolation denotes a family of meshless interpolation schemes in which the approximation at a target location is constructed from a restricted neighborhood, subdomain, or stencil rather than from all data sites globally. The common motivation is computational and numerical: classical RBF interpolation uses the ansatz \[ s(x)=\sum_{j=1}^n y_j\,\phi(\|x-x_j\|), \] or, in polynomially augmented form, \[ f(x) = \sum_{j=1}^{N} \lambda_j \, \phi(\|x - x_j\|) + P_k(x), \] and enforcing interpolation yields a dense linear system that is expensive to solve directly [1006.2453], [1707.09881]. Local formulations alter this picture in different ways: by using compactly supported radial basis functions so that \(\phi(r_{ij})=0\) beyond a finite radius [1707.09881], by decomposing the domain into overlapping patches and solving only small local systems combined through partition of unity weights [1811.05193], [1408.0183], by forming pointwise neighborhoods for each query location [1908.03690], or by using local stencils as in RBF-generated finite differences [1812.04742]. Although these variants differ in algebraic structure, they share the aim of replacing the dense, globally coupled interpolation problem by sparse, localized, or multilevel computations while retaining the meshless character of RBF approximation. ## 1. Classical formulation and the case for localization The standard interpolation problem begins with scattered or gridded samples \(\{(x_i,h_i)\}_{i=1}^N\) or \(\{x_j,f_j\}\), and seeks an interpolant expressed as a linear combination of radial kernels evaluated at pairwise distances [1006.2453], [1707.09881]. In one common block formulation with polynomial augmentation, \[ \begin{bmatrix} \Phi & P \\ P^T & 0 \end{bmatrix} \begin{bmatrix} \lambda \\ a \end{bmatrix} = \begin{bmatrix} h \\ 0 \end{bmatrix}, \] where \(\Phi_{ij}=\phi(\|x_i-x_j\|)\) and \(P\) contains polynomial basis terms [1707.09881]. This framework is meshless and depends on distances rather than on a mesh or triangulation, which is one reason RBF methods are used for scattered \(d\)-dimensional interpolation [1707.09881]. The principal difficulty is that global kernels such as Gaussian, multiquadric, inverse multiquadric, and thin-plate spline induce full dense matrices [1606.07652], [1707.09881]. For large point sets, the storage and solution cost of these matrices is prohibitive. One paper states that for \(N \approx 10^6\), a dense \(N \times N\) matrix is clearly infeasible, while the dense/global RBF costs are roughly \(O(N^2)\) in memory and \(O(N^3)\) for a direct solve [1707.09881]. Another paper emphasizes that even matrix-vector multiplication for dense RBF interpolation without fast summation methods costs \(\mathcal{O}(N^2)\) [1606.07652]. Local RBF interpolation is therefore best understood as a collection of strategies for reducing global coupling. In some approaches, locality is imposed directly through the kernel support; in others, it is imposed through the computational domain decomposition, neighborhood selection, or transformed basis. This suggests that “local” in the RBF literature is not a single construction but a structural principle: the influence of each basis function, patch, or stencil is restricted in space or in algebraic support. ## 2. Compact support, sparsity, and neighborhood-limited influence A direct route to locality uses compactly supported radial basis functions. In that setting, the basis is defined on a normalized support interval by \[ \phi(r)= \begin{cases} (1-r)^q P(r), & 0 \le r \le 1,\\ 0, & r>1, \end{cases} \] with a scaling factor \(a>0\) used in practice so that the support radius is not fixed to \(1\) [1707.09881]. The decisive local property is explicit: if two points are farther apart than the support radius, then \[ \phi(r_{ij})=0. \] Each basis function therefore influences only nearby points, and the interpolation matrix becomes sparse rather than dense [1707.09881]. This sparse-versus-dense distinction is central to the large-data literature. Global RBFs such as thin-plate splines produce a full dense matrix \(\Phi\), whereas CSRBFs produce a sparse matrix because many entries are exactly zero [1707.09881]. The practical consequence is not merely reduced memory; sparse direct or iterative solvers become viable, and matrix-vector operations are cheaper. The same source emphasizes that CSRBFs are especially advantageous when the data set is large, the points are scattered or unstructured, a mesh or triangulation would be expensive or impractical, the problem is higher-dimensional, memory limits make dense RBF matrices infeasible, or local interpolation behavior is desired [1707.09881]. Compact support does not eliminate all numerical issues. The same analysis warns that large coordinate span harms numerical stability because the polynomial block \(P\) depends strongly on the absolute positions of points. If point coordinates are far from the origin or span a large interval, polynomial terms such as \(x\), \(y\), \(x^2\), \(xy\), and \(y^2\) can vary over many orders of magnitude, causing large ranges in entries of \(P^T P\), widely spread eigenvalues, ill-conditioning, and reduced numerical robustness [1707.09881]. A practical implication, stated there explicitly, is that translation and scaling of the point set matter. This localized-support view of local RBF interpolation should be distinguished from computational locality induced by fast summation. A band-limited fast multipole method for RBF interpolation introduces near-field direct interactions and far-field low-rank approximations, but the underlying interpolation remains global [1606.07652]. The locality in that case is “local in the FMM sense,” not local in the sense of compact support or patchwise independence. ## 3. Partition of unity and patchwise local approximants A second major line of local RBF interpolation is the partition of unity (PU) method. Here the domain is covered by overlapping subdomains \(\{\Omega_j\}_{j=1}^d\), each local point set \({\cal X}_{N_j}={\cal X}_N\cap\Omega_j\) defines a small interpolation problem, and the global interpolant is assembled as \[ {\cal I}(\boldsymbol{x})=\sum_{j=1}^d R_j(\boldsymbol{x})\,W_j(\boldsymbol{x}), \] where \(W_j\) are compactly supported PU weights satisfying \[ \sum_{j=1}^d W_j(\boldsymbol{x})=1,\qquad \operatorname{supp}(W_j)\subseteq \Omega_j \] [1811.05193]. This construction replaces one large dense solve by many smaller local systems \(A_j\boldsymbol{\alpha}_j=\boldsymbol{f}_j\) [1811.05193]. The PU framework is explicitly described as effective when large data sets are considered, precisely because only linear systems of relatively small size need to be solved [1811.05193]. In two-dimensional interpolation on \([0,1]^2\), a related cell-based PU algorithm uses local RBF interpolants \[ {\cal R}_j(x)=\sum_{k=1}^{m_j} c_k\,\phi(\|x-x_k\|_2) \] and blends them by \[ {\cal I}(x)=\sum_{j=1}^d {\cal R}_j(x)\,W_j(x) \] [1408.0183]. The paper emphasizes that if each local interpolant matches the data on the nodes in its subdomain, then the PU interpolant preserves interpolation because the weights sum to one [1408.0183]. Theoretical analysis in the PU setting depends on a regular covering: bounded overlap, an interior cone condition for each subdomain, and uniform control of local fill distances by the global fill distance [1811.05193], [1408.0183]. The PU weights are required to be \(k\)-stable in Wendland’s sense [1811.05193]. These assumptions allow the local errors to control the global approximation error. An important refinement replaces fixed spherical patches by variable-shaped ellipsoidal subdomains paired with anisotropic kernels [1811.05193]. In that approach, each subdomain \(\Omega_j\) is an ellipsoid specified by semi-axes \[ \boldsymbol{\delta}_j=(\delta_1^j,\dots,\delta_M^j), \] and each local shape parameter becomes a diagonal positive definite matrix \[ E_j=\operatorname{diag}(\varepsilon_1^j,\dots,\varepsilon_M^j). \] The rationale is that for track data and other anisotropically distributed point clouds, spherical patches can waste points in irrelevant directions and fail to reflect the local geometry [1811.05193]. Local parameter selection is then driven by a LOOCV error estimate \[ e_j=e_j(\boldsymbol{\varepsilon}_j,\boldsymbol{\delta}_j), \] optimized patch by patch via the Nelder–Mead simplex method implemented in MATLAB’s `fminsearch.m` [1811.05193]. The reported numerical evidence is strong: for artificial track data with \(N=1000\), classical PU yields RMSE \(=3.27\times 10^{-3}\), whereas PU-LOOCV yields RMSE \(=8.02\times 10^{-5}\) [1811.05193]. A distinct computational contribution within PU methods is the use of geometric search structures. In the cell-based 2D algorithm, the domain is partitioned into square cells with side length equal to the subdomain radius, so that each search inspects only a fixed \(3\times 3\) neighborhood of cells [1408.0183]. The paper reports CPU-time reductions relative to an earlier strip-based method; for \(n=66049\) and \(d=16384\), \(t_{\text{cell}}=2.6\) versus \(t_{\text{strip}}=6.5\) [1408.0183]. This illustrates a recurring theme in local RBF interpolation: localization is often as much about search and data structures as about the kernel itself. ## 4. Pointwise local interpolation and adaptive neighborhood selection Another form of locality is pointwise neighborhood interpolation: for each target location, a small local data set is selected and a separate RBF interpolant is built only on that neighborhood. A concrete example is the adaptive local MQ-RBF method for missing-value estimation in geographical data [1908.03690]. For each interpolated point, the algorithm computes its distance to all known data points, selects the 20 points with the smallest distances, and uses these \(N_{\text{loc}}=20\) samples as the local interpolation data set [1908.03690]. The local interpolant uses the multiquadric basis \[ \phi(r)=\sqrt{r^2+c^2}, \] with Euclidean distance \(r=\|x-x_j\|_2\) [1908.03690]. The distinctive feature is adaptive shape-factor selection based on the distribution of the local neighbors. The method defines the expected density over the full data set, \[ D_{\exp}=\frac{N_{dp}}{(X_{\max}-X_{\min})(Y_{\max}-Y_{\min})}, \] the local density in the neighborhood, \[ D_{loc}=\frac{N_{loc}}{(x_{\max}-x_{\min})(y_{\max}-y_{\min})}, \] and the density statistic \[ D(s_0)=\frac{D_{loc}}{D_{\exp}}. \] This quantity is normalized by a fuzzy membership function to \(\mu_D\), after which the shape factor \(c\) is chosen through a piecewise triangular membership rule with five preset levels \(c_1,\dots,c_5\) [1908.03690]. The workflow is explicitly local and repeated independently for each missing point. The method was evaluated on three DEM data sets from Beijing, Chongqing, and Longyan, each with randomly masked \(10\%\) of observed samples as missing values. The reported data sizes are: Beijing, 1,111,369 known values and 123,592 missing values; Chongqing, 1,074,379 known values and 97,525 missing values; Longyan, 1,040,670 known values and 119,050 missing values [1908.03690]. Relative to \(k\)NN interpolation and AIDW, the adaptive RBF method has the best accuracy in terms of RMSE, but is the least efficient computationally [1908.03690]. The source further notes that because the test DEM points are relatively evenly distributed, the advantages of adaptive shape tuning may not be fully realized, and may show more advantage on unevenly distributed data [1908.03690]. This pointwise model differs from PU in a consequential way. PU constructs a globally defined interpolant by blending local approximants; the adaptive local MQ-RBF method instead performs target-by-target imputation. A plausible implication is that pointwise local RBF interpolation is especially natural for missing-value recovery and localized prediction tasks, whereas PU is more naturally suited to constructing a continuous global field over the entire domain. ## 5. Local stencils, finite differences, and transformed bases Local RBF interpolation also appears in discretization schemes where the objective is not a global interpolant but the generation of local differentiation weights. In RBF-generated finite differences for the 2D Helmholtz equation, each node \(x_i\) is assigned a stencil \[ S_i=\{x_j^i\}_{j=1}^{n_i},\qquad x_1^i=x_i, \] formed from the \(n_i\) nearest neighbors [1812.04742]. On this stencil, the local interpolant is \[ \widetilde{u}(x)=\sum_{j=1}^{n_i}\alpha_j^i\,J_0\!\left(k_i\|x-x_j^i\|\right), \] where \(J_0\) is the Bessel function of the first kind used as an oscillatory radial basis [1812.04742]. The local coefficients are found from the stencil interpolation matrix \[ \mathbf{J}_{k_i}=\left(J_0\!\left(k_i\|x_l^i-x_j^i\|\right)\right)_{1\le l,j\le n_i}, \] and the differential-operator weights are then derived from this local interpolant [1812.04742]. The paper stresses that local interpolation can itself be ill-conditioned. For the Bessel RBF stencils, condition numbers can exceed \(10^{15}\), especially for larger wavenumbers, symmetric stencil geometries, and regular square or hexagonal node patterns [1812.04742]. The stabilization device is a diagonal perturbation, \[ \widetilde{\mathbf{J}_{k_i}}=\mathbf{J}_{k_i}+\beta I, \] with \(\beta\) chosen so that the perturbed matrix has a prescribed target condition number \(\kappa_0\) [1812.04742]. The reported convergence rates are about third order and sixth order depending on stencil choice and grid structure, with approximately \(\|u-\widetilde{u}\|_\infty \sim C h^{6.2}\) on a square grid and rates around \(h^{3.1}\) to \(h^{3.9}\) on a hexagonal grid in some cases [1812.04742]. A broader multilevel interpretation of locality appears in the hierarchical-basis solver for RBF interpolation [1104.2504]. There, the kernel may remain globally supported, but an orthonormal hierarchical basis \(P=[L\;T]\) is constructed so that the detail basis \(T\) annihilates all polynomials up to a chosen degree and is adapted to the node geometry via an octree-like multiresolution partition [1104.2504]. The authors state that the transformed system behaves like a localized, multilevel approximation rather than a single dense global solve. This produces a transformed matrix \(K_W=T^HKT\) whose condition number is essentially invariant to domain scaling; for the biharmonic example with \(m=3\), the transformed matrix had condition number about \(762\) across multiple scales, while the original system condition number reached values like \(10^{24}\) [1104.2504]. Combined with fast summation, the total solve is reported as roughly \(O(N\log^2 N)\) under the paper’s assumptions [1104.2504]. This suggests a useful conceptual distinction. Some local RBF methods are local in the strong sense that only nearby nodes participate in each solve. Others remain globally supported at the kernel level but reveal effective locality after basis transformation or multilevel decomposition. Both belong to the broader efficient-localized RBF tradition, but they address dense coupling in different ways. ## 6. Geometry-aware locality, structured grids, and parameter tuning Locality based on Euclidean distance can fail on domains with nontrivial geometry. A recent example is rescaled localized radial basis function interpolation based on geodesic distance, RL-RBF-G [2403.03665]. Standard localized RBF interpolation on non-matching meshes uses compactly supported Wendland \(C^2\) basis functions, \[ \phi(t,r)= \max\left\{1-\frac{t}{r},0\right\}^4\left(1+4\frac{t}{r}\right), \] with support radius chosen locally as \[ r_j = \alpha\, d\bigl(\mathbf x_j,\mathbf z_j\bigr), \] where \(\mathbf z_j\) is the \(M\)-th nearest neighbor of \(\mathbf x_j\) [2403.03665]. The interpolant is rescaled by interpolating the constant function \(g(\mathbf{x})=1\) and dividing: \[ f_{\mathrm{res}}(\mathbf x)=\frac{\hat f(\mathbf x)}{\hat g(\mathbf x)}. \] This RL-RBF formulation is already local because the basis functions are compactly supported and each support radius is chosen from local point spacing [2403.03665]. RL-RBF-G modifies the distance notion itself. Instead of using Euclidean distance alone, it employs a thresholded hybrid distance based on a discrete geodesic \(g_h\), \[ d(\mathbf x,\mathbf y;R)= \begin{cases} +\infty, & g_h(\mathbf x,\mathbf y)>R,\\ g_h(\mathbf x,\mathbf y), & \beta h_{\max}+\|\mathbf x-\mathbf y\|