- •Введение
- •Основные возможности пакета Surfer
- •1.2. Основные возможности модуля геостатического анализа (Geostatical Analyst) Гис-пакета ArcGis 9.3.
- •В чем важность ArcGis Geostatistical Analyst?
- •Кому пригодится ArcGis Geostatistical Analyst?
- •4. Каковы функции ArcGis Geostatistical Analyst?
- •Исследование данных и Построение поверхностей
- •12. Что такое Исследовательский анализ пространственных данных?
- •13. Что такое методы интерполяции?
- •14. Как выбрать подходящий метод интерполяции?
- •15. Включены ли в Geostatistical Analyst методы интерполяции, основанные на сплайнах?
- •16. Что такое автокорреляция?
- •17. Что такое кригинг?
- •18. Что такое вариаограмма (полувариограмма)?
- •19. Что такое моделирование ошибок?
- •20. Может ли Geostatistical Analyst выявить ошибки в моих данных?
- •21. Нужно ли дополнительное программное обеспечение для оценки оптимальных параметров прогнозных моделей, включенных в Geostatistical Analyst?
- •22. Можно ли проводить в Geostatistical Analyst поблоковый кригинг?
- •1.3. Подходы к построению карт
- •Подходы к анализу пространственно распределенных данных
- •1.4. Детерминистские методы интерполяции
- •Визуализация интерполяции по методу глобального полинома
- •Визуализация интерполяции по методу локальных полиномов
- •Визуализация радиальных базисных функций
- •Метод взвешенных (обратных) расстояний
- •Метод радиальных базисных функций
- •Когда использовать радиальные базисные функции
- •1.5. Геостатические методы интерполяции
- •2. Методы построения сеточных моделей
- •2.1. Метод Крайгинга
- •Метод радиальных базисных функций
- •2.3. Метод обратных расстояний
- •2.4. Метод Шепарда
- •2.5. Метод минимальной кривизны
- •2.6. Метод полиномиальной регрессии
- •2.7. Триангуляция с линейной интерполяцией
- •2.8. Метод ближайшего соседа
- •2.9. Метод естественного соседа
- •2. 10. Сравнительная характеристика основных методов построения сеточной функции
Метод радиальных базисных функций
В этом методе искомая функция находится как линейная комбинация набора радиальных базисных функций:
, (2.3)
где а -константа, i - индекс точки измерений, i - неизвестные коэффициенты, Ri(x,y) - базисные функции, зависящие от расстояния точки (х, у) до i-ой точки наблюдения.
Существуют несколько типов базисных функций:
Inverse Multiquadric
Multilog
Мультиквадратичная (Multiquadric), наиболее часто используется;
Natural Cubic Spline
Thin Plate Spline ,
где R2 – фактор сглаживания, чем больше будет параметр, тем более сглаженные будут контура. Разумные значения показателя находятся в интервале от среднего межточечного расстояния выборки до половины этого среднего значения.
2.3. Метод обратных расстояний
Фактически большинство интерполяционных методов используют формулу (2.1) для расчета значений функций в произвольных точках. Различаются принципы, используемые для задания весовых коэффициентов wi и, соответственно, выражения для них расчета. В методе обратных расстояний используется обратно пропорциональная зависимость весовых коэффициентов от некоторой степени расстояния между расчетной точкой и точкой наблюдения (di):
, (2.4)
здесь к - число, обычно принимаемое 1, 2 или 3.
Такой подход имеет вполне очевидный смысл - более удаленные точки в меньшей степени определяют значение в расчетной точке и наоборот.
2.4. Метод Шепарда
Метод Шепарда аналогичен методу обратных расстояний. Он также использует обратные расстояния при вычислении весовых коэффициентов. Интерполяция же осуществляется по формуле (2.3) при к=2. Отличие состоит в том, что при построении интерполяционной функции в локальных областях используется метод наименьших квадратов. Это уменьшает вероятность появления на сгенерированной поверхности ложных структур вокруг точек наблюдений.
2.5. Метод минимальной кривизны
Поверхность, построенная с помощью этого метода, аналогична тонкой упругой пластине, проходящей через все экспериментальные точки данных с минимальным числом изгибов. Он генерирует наиболее гладкую поверхность, которая проходит настолько близко к экспериментальным точкам, насколько это возможно, но экспериментальные точки не обязательно принадлежать приближающей поверхности. Окно Parameters позволяет задать критерий сходимости.
Параметр Max Residuals (Максимальная невязка) измеряется в тех же единицах, что и экспериментальные данные. Как показывает опыт, наиболее подходящее значение этого параметра равно, примерно, 10% от точности исходных данных. Если исходные экспериментальные данные измерены с точностью 1.0 единиц измерения, то значение Max Residuals рекомендуется положить равным 0.1. Итерации продолжаются до тех пор, пока максимальная невязка для всей итерации не станет меньше значения параметра Max Residuals.
Параметр Max Iterations (Максимальное число итераций) разумно выбрать в интервале от N до 2N, где N - число узлов генерируемой сети.