Многомерная интерполяция

Интерполяция функций более чем одной переменной

В численном анализе многомерная интерполяция — это интерполяция функций более чем одной переменной ( многомерные функции ); когда переменные являются пространственными координатами , она также известна как пространственная интерполяция .

Интерполируемая функция известна в заданных точках , а задача интерполяции состоит в получении значений в произвольных точках . ${\displaystyle (x_{i},y_{i},z_{i},\dots )}$ ${\displaystyle (x,y,z,\dots )}$

Многомерная интерполяция особенно важна в геостатистике , где она используется для создания цифровой модели рельефа из набора точек на поверхности Земли (например, высот точек при топографической съемке или глубин при гидрографической съемке ).

Регулярная сетка

Сравнение некоторых 1- и 2-мерных интерполяций.
Черные и красные / желтые / зеленые / синие точки соответствуют интерполированной точке и соседним образцам соответственно.
Их высоты над землей соответствуют их значениям.

Для значений функции, известных на регулярной сетке (имеющей заранее определенный, не обязательно равномерный интервал), доступны следующие методы.

Любое измерение

• Интерполяция по методу ближайшего соседа
• n-линейная интерполяция (см. би- и трилинейная интерполяция и полилинейный полином )
• n-кубическая интерполяция (см. би- и трикубическая интерполяция )
• Кригинг
• Обратное взвешивание расстояний
• Интерполяция естественного соседа
• Сплайн-интерполяция
• Интерполяция радиальной базисной функции

2 измерения

• интерполяция Барнса
• Билинейная интерполяция
• Бикубическая интерполяция
• Поверхность Безье
• Повторная выборка Ланцоша
• триангуляция Делоне

Повторная выборка растровых изображений — это применение двумерной многомерной интерполяции при обработке изображений .

Три метода применены к одному и тому же набору данных, из 25 значений, расположенных в черных точках. Цвета представляют интерполированные значения.

• 
  Ближайший сосед
• 
  Билинейный
• 
  Бикубический

См. также Падуанские точки для полиномиальной интерполяции по двум переменным.

3 измерения

• Трилинейная интерполяция
• Трикубическая интерполяция

См. также передискретизация растрового изображения .

Сплайны тензорного произведения дляНразмеры

Сплайны Катмулла-Рома можно легко обобщить на любое количество измерений. Статья о кубических сплайнах Эрмита напомнит вам, что для некоторого 4-вектора , который является функцией только x , где — значение в интерполируемой функции. Перепишите это приближение как ${\displaystyle \mathrm {CINT} _{x}(f_{-1},f_{0},f_{1},f_{2})=\mathbf {b} (x)\cdot \left(f_{-1}f_{0}f_{1}f_{2}\right)}$ ${\displaystyle \mathbf {b} (x)}$ ${\displaystyle f_{j}}$ ${\displaystyle j}$

${\displaystyle \mathrm {CR} (x)=\sum _{i=-1}^{2}f_{i}b_{i}(x)}$

Эту формулу можно напрямую обобщить на N измерений: ^[1]

${\displaystyle \mathrm {CR} (x_{1},\dots ,x_{N})=\sum _{i_{1},\dots ,i_{N}=-1}^{2}f_{i_{1}\dots i_{N}}\prod _{j=1}^{N}b_{i_{j}}(x_{j})}$

Обратите внимание, что подобные обобщения можно сделать и для других типов сплайн-интерполяций, включая эрмитовы сплайны. Что касается эффективности, то общая формула фактически может быть вычислена как композиция последовательных операций -типа для любого типа тензорных сплайнов, как объяснено в статье о трикубической интерполяции . Однако факт остается фактом: если в одномерном -подобном суммировании есть члены, то в -мерном суммировании будут члены . ${\displaystyle \mathrm {CINT} }$ ${\displaystyle n}$ ${\displaystyle \mathrm {CR} }$ ${\displaystyle n^{N}}$ ${\displaystyle N}$

Нерегулярная сетка (разбросанные данные)

Схемы, определенные для разбросанных данных на нерегулярной сетке, являются более общими. Они все должны работать на регулярной сетке, обычно сводясь к другому известному методу.

• Интерполяция по методу ближайшего соседа
• Триангулированная нерегулярная сеть на основе естественного соседа
• Линейная интерполяция на основе триангулированной нерегулярной сети (тип кусочно-линейной функции )
  • n- симплексная (например, тетраэдрическая) интерполяция (см. барицентрическую систему координат )
• Обратное взвешивание расстояний
• ABOS - аппроксимация на основе сглаживания
• Кригинг
• Градиентный кригинг (GEK)
• Тонкая пластина-шлиц
• Полигармонический сплайн (тонкопленочный сплайн является частным случаем полигармонического сплайна)
• Радиальная базисная функция ( полигармонические сплайны являются частным случаем радиальных базисных функций с полиномиальными членами низкой степени)
• Сплайн наименьших квадратов
• Интерполяция естественного соседа

Сетка — это процесс преобразования нерегулярно распределенных данных в регулярную сетку ( сетчатые данные ).

Смотрите также

• Сглаживание
• Поверхностная подгонка

Примечания

1. Две иерархии сплайн-интерполяций. Практические алгоритмы для многомерных сплайнов высшего порядка

Внешние ссылки

• Пример кода C++ для нескольких 1D, 2D и 3D сплайн-интерполяций (включая сплайны Кэтмелла-Рома).
• Многомерная интерполяция и аппроксимация Эрмита, проф. Чандраджит Баджаджа, Университет Пердью
• Библиотека Python, содержащая методы 3D- и 4D-сплайн-интерполяции.