e01sg:: Interpolation (NAG Toolbox)

nag_interp_2d_scat_shep (e01sg) constructs a smooth function

Q (x, y)

which interpolates a set of

m

scattered data points

(x_{r}, y_{r}, f_{r})

, for

r = 1, 2, \dots, m

, using a modification of Shepard's method. The surface is continuous and has continuous first partial derivatives.

The basic Shepard (1968) method interpolates the input data with the weighted mean

Q (x, y) = \frac{\sum_{r = 1}^{m} w_{r} (x, y) q_{r}}{\sum_{r = 1}^{m} w_{r} (x, y)},

where

q_{r} = f_{r}

w_{r} (x, y) = \frac{1}{d_{r}^{2}}

and

d_{r}^{2} = {(x - x_{r})}^{2} + {(y - y_{r})}^{2}

The basic method is global in that the interpolated value at any point depends on all the data, but this function uses a modification (see Franke and Nielson (1980) and Renka (1988a)), whereby the method becomes local by adjusting each

w_{r} (x, y)

to be zero outside a circle with centre

(x_{r}, y_{r})

and some radius

R_{w}

. Also, to improve the performance of the basic method, each

q_{r}

above is replaced by a function

q_{r} (x, y)

, which is a quadratic fitted by weighted least squares to data local to

(x_{r}, y_{r})

and forced to interpolate

(x_{r}, y_{r}, f_{r})

. In this context, a point

(x, y)

is defined to be local to another point if it lies within some distance

R_{q}

of it. Computation of these quadratics constitutes the main work done by this function.

The radii

R_{w}

and

R_{q}

are chosen to be just large enough to include

N_{w}

and

N_{q}

data points, respectively, for user-supplied constants

N_{w}

and

N_{q}

. Default values of these arguments are provided by the function, and advice on alternatives is given in Choice of and .

References

Parameters

Compulsory Input Parameters

Optional Input Parameters

Output Parameters

Error Indicators and Warnings

Accuracy

Further Comments

Timing

Choice of Nw and Nq

If non-default values are required they must be supplied to nag_interp_2d_scat_shep (e01sg) through positive values of nw and nq. Increasing these arguments makes the method less local. This may increase the accuracy of the resulting interpolant at the expense of increased computational cost. The default values

nw = \min (19, m - 1)

and

nq = \min (13, m - 1)

have been chosen on the basis of experimental results reported in Renka (1988a). In these experiments the error norm was found to vary smoothly with

N_{w}

and

N_{q}

, generally increasing monotonically and slowly with distance from the optimal pair. The method is not therefore thought to be particularly sensitive to the argument values. For further advice on the choice of these arguments see Renka (1988a).

Example

function e01sg_example


fprintf('e01sg example results\n\n');

% Scattered Grid Data
x = [11.16; 12.85; 19.85; 19.72; 15.91;  0.00; 20.87;  3.45; 14.26; ...
     17.43; 22.80;  7.58; 25.00;  0.00;  9.66;  5.22; 17.25; 25.00; ...
     12.13; 22.23; 11.52; 15.20;  7.54; 17.32;  2.14;  0.51; 22.69; ...
      5.47; 21.67;  3.31];
y = [ 1.24;  3.06; 10.72;  1.39;  7.74; 20.00; 20.00; 12.78; 17.87; ...
      3.46; 12.39;  1.98; 11.87;  0.00; 20.00; 14.66; 19.57;  3.87; ...
     10.79;  6.21;  8.53;  0.00; 10.69; 13.78; 15.03;  8.37; 19.63; ...
     17.13; 14.36; 0.33];
f = [22.15; 22.11;  7.97; 16.83; 15.30; 34.60;  5.74; 41.24; 10.74; ...
     18.60;  5.47; 29.87;  4.40; 58.20;  4.73; 40.36;  6.43;  8.74; ...
     13.71; 10.25; 15.74; 21.60; 19.31; 12.11; 53.10; 49.43;  3.25; ...
     28.63;  5.52; 44.08];

% Generate interpolant
nw = int64(0);
nq = int64(0);
[iq, rq, ifail] = e01sg(x, y, f, nw, nq);

% Interpolation points
u = [20.00;  6.41;  7.54;  9.91; 12.30];
v = [ 3.14; 15.44; 10.69; 18.27;  9.22];

% Interpolate at interpolation points 
[q, qx, qy, ifail] = e01sh(x, y, f, iq, rq, u, v);

fprintf('Interpolated values Q and its derivatives at (u,v)\n'); 
fprintf('     u      v      q      qx     qy\n');
for i = 1:size(u,1)
  fprintf('%7.2f%7.2f%7.2f%7.2f%7.2f\n', u(i), v(i), q(i), qx(i), qy(i));
end

e01sg example results

Interpolated values Q and its derivatives at (u,v)
     u      v      q      qx     qy
  20.00   3.14  15.89  -1.28  -0.63
   6.41  15.44  34.05  -3.62  -3.56
   7.54  10.69  19.31  -2.84   0.81
   9.91  18.27  13.68  -1.59  -4.71
  12.30   9.22  14.56  -0.68  -0.78

On entry,	$m < 6$ ,
or	$0 < nq < 5$ ,
or	$nq > \min (40, m - 1)$ ,
or	$nw > \min (40, m - 1)$ ,
or	$liq < 2 \times m + 1$ ,
or	$lrq < 6 \times m + 5$ .

NAG Toolbox: nag_interp_2d_scat_shep (e01sg)

▸▿ Contents

Purpose

Syntax

Description