Polyharmonic Spline - Definition

Definition

y(\mathbf{x}) \, = \, \sum_{i=1}^N w_i \, \phi(||\mathbf{x} - \mathbf{c}_i||) + \mathbf{v}^T \, \begin{bmatrix} 1 \\ \mathbf{x} \end{bmatrix}

where

  • is a real-valued vector of nx independent variables,
  • are N vectors of the same size as (often called centers) that the interpolated curve shall pass
  • are the N weights of the basis functions.
  • are the nx+1 weights of the polynomial.
  • The linear polynomial with the weighting factors improves the interpolation close to the "boundary" and especially the extrapolation "outside" of the centers . If this is not desired, this term can also be removed (see also figure below).

The basis functions of polyharmonic splines are radial basis functions of the form:

\begin{matrix} \phi(r) = \begin{cases} r^k & \mbox{with } k=1,3,5,\dots, \\ r^k \ln(r) & \mbox{with } k=2,4,6,\dots \end{cases} \\ r = ||\mathbf{x} - \mathbf{c}_i||_2 = \sqrt{ (\mathbf{x} - \mathbf{c}_i)^T \, (\mathbf{x} - \mathbf{c}_i) } \end{matrix}

Other values of exponent k are not useful (such as ), because a solution of the interpolation problem might no longer exist. To avoid problems at r=0 (since ln(0) = -∞), the polyharmonic splines with the natural logarithm might be implemented as:

\phi(r) = \begin{cases} r^{k-1} \ln(r^r) & \mbox{for } r < 1 \\ r^k \ln(r) & \mbox{for } r \ge 1 \end{cases}

The weights and are determined such that the function passes through given points (i=1,2,...,N) and fulfill the orthogonality conditions:

0 = \sum_{i=1}^N w_i, \;\; 0 = \sum_{i=1}^N w_i \, c_{j,i} \;\;\; (j=1,2,...,nx)

To compute the weights, a symmetric, linear system of equations has to be solved:

\begin{bmatrix} \mathbf{A} & \mathbf{V}^T \\ \mathbf{V} & \mathbf{0} \end{bmatrix} \; \begin{bmatrix} \mathbf{w} \\ \mathbf{v} \end{bmatrix} \; = \; \begin{bmatrix} \mathbf{y} \\ \mathbf{0} \end{bmatrix}\;\;\;\;

where

A_{i,j} = \phi(||\mathbf{c}_i - \mathbf{c}_j||), \;\;\; \mathbf{V} = \begin{bmatrix} 1 & 1 & \cdots & 1 \\ \mathbf{c}_1 & \mathbf{c}_2 & \cdots & \mathbf{c}_{N} \end{bmatrix}, \;\;\; \mathbf{y} = ^T

Under very mild conditions (essentially, that at least nx+1 points are not in a subspace; e.g. for nx=2 that at least 3 points are not on a straight line), the system matrix of the linear system of equations is nonsingular and therefore a unique solution of the equation system exists.

Once the weights are determined, interpolation requires to just evaluate the top most formula for the provided .

Many practical details to implement and use polyharmonic splines are given in the book of Fasshauer. In Iske polyharmonic splines are treated as special cases of other multiresolution methods in scattered data modelling.

Read more about this topic:  Polyharmonic Spline

Famous quotes containing the word definition:

    “Scientific method is the way to truth, but it affords, even in
    principle, no unique definition of truth. Any so-called pragmatic
    definition of truth is doomed to failure equally.”
    —Willard Van Orman Quine (b. 1908)

    “... if, as women, we accept a philosophy of history that asserts that women are by definition assimilated into the male universal, that we can understand our past through a male lens—if we are unaware that women even have a history—we live our lives similarly unanchored, drifting in response to a veering wind of myth and bias.”
    —Adrienne Rich (b. 1929)

    “The physicians say, they are not materialists; but they are:MSpirit is matter reduced to an extreme thinness: O so thin!—But the definition of spiritual should be, that which is its own evidence. What notions do they attach to love! what to religion! One would not willingly pronounce these words in their hearing, and give them the occasion to profane them.”
    —Ralph Waldo Emerson (1803–1882)