30.4 Интерполяция на рассеянных данных

Важное применение триангуляции Делоне заключается в том, что она может использоваться для интерполяции из рассеянных данных в произвольную совокупность точек. Для этого вычисляется N-симплекс известного набора точек с помощью delaunay или delaunayn. Затем идентифицируются симплексы, в которых находятся искомые точки. Наконец, вершины симплексов используются для интерполяции в искомые точки. Функции, выполняющие эту интерполяцию, это griddata, griddata3 и griddatan.

zi = griddata (x, y, z, xi, yi)

zi = griddata (x, y, z, xi, yi, method)

[xi, yi, zi] = griddata (…)

Генерирует регулярную сетку из нерегулярных данных с использованием интерполяции.

Функция определяется по z = f (x, y). Входные данные x, y, z являются векторами одинаковой длины или x, y являются векторами, а z является матрицей.

Точки интерполяции все (xi, yi). Если xi, yi являются векторами, то они преобразуются в 2-мерную сетку.

Метод интерполяции может быть "nearest", "cubic" или "linear". Если метод опущен, он по умолчанию равен "linear".

См. также: griddata3, griddatan, delaunay.

vi = griddata3 (x, y, z, v, xi, yi, zi)

vi = griddata3 (x, y, z, v, xi, yi, zi, method)

vi = griddata3 (x, y, z, v, xi, yi, zi, method, options)

Генерирует регулярную сетку из нерегулярных данных с использованием интерполяции.

Функция определяется по v = f (x, y, z). Точки интерполяции задаются xi, yi, zi.

Метод интерполяции может быть "nearest" или "linear". Если метод опущен, он по умолчанию равен "linear".

Необязательный аргумент options передаётся непосредственно в Qhull при вычислении триангуляции Делоне, используемой для интерполяции. См. delaunayn для получения информации о значениях по умолчанию и о том, как передавать различные значения.

См. также: griddata, griddatan, delaunayn.

yi = griddatan (x, y, xi)

yi = griddatan (x, y, xi, method)

yi = griddatan (x, y, xi, method, options)

Генерирует регулярную сетку из нерегулярных данных с использованием интерполяции.

Функция определяется по y = f (x). Точки интерполяции — все xi.

Метод интерполяции может быть "nearest" или "linear". Если метод опущен, он по умолчанию равен "linear".

Необязательный аргумент options передаётся непосредственно в Qhull при вычислении триангуляции Делоне, используемой для интерполяции. См. delaunayn для получения информации о значениях по умолчанию и о том, как передавать различные значения.

См. также: griddata, griddata3, delaunayn.

Пример использования функции griddata:

rand ("state", 1);
x = 2*rand (1000,1) - 1;
y = 2*rand (size (x)) - 1;
z = sin (2*(x.^2+y.^2));
[xx,yy] = meshgrid (linspace (-1,1,32));
zz = griddata (x, y, z, xx, yy);
mesh (xx, yy, zz);

интерполирует из случайного рассеяния точек в равномерную сетку. Результат выше можно увидеть на рисунке 30.6.

Рисунок 30.6: Интерполяция из рассеянных данных на регулярную сетку