Интерполяция обратным взвешиванием расстояний: от SYMAP до современных методов
Inverse distance weighting
Интерполяция методом обратных расстояний (IDW): детерминированный способ многомерной интерполяции, взвешенное среднее для оценки значений в точках. GIS.
Сравнивайте с английским: нажмите на абзац — оригинал откроется в окне. Кнопка EN под абзацем показывает его прямо в тексте.
Содержание
Введение
Тип детерминированного метода для многовариантной интерполяции
Type of deterministic method for multivariate interpolation
Инверсное взвешивание по расстоянию (IDW) — это тип детерминированного метода для многовариантной интерполяции по известному рассеянному набору точек. Значения, присваиваемые неизвестным точкам, рассчитываются как взвешенное среднее значений, доступных в известных точках. Этот метод также может использоваться для создания матриц пространственных весов при анализе пространственной автокорреляции (например, I Морана). Название этого метода обусловлено используемым взвешенным средним, поскольку для определения весов используется обратная величина расстояния до каждой известной точки ("мера близости").
Inverse distance weighting (IDW) is a type of deterministic method for multivariate interpolation with a known scattered set of points. The assigned values to unknown points are calculated with a weighted average of the values available at the known points. This method can also be used to create spatial weights matrices in spatial autocorrelation analyses (e. g. Moran's I). The name given to this type of method was motivated by the weighted average applied, since it resorts to the inverse of the distance to each known point ("amount of proximity") when assigning weights.
Историческая справка
В Гарвардской лаборатории компьютерной графики и пространственного анализа, начиная с 1965 года, собралась разнородная группа ученых, чтобы переосмыслить, в том числе, то, что сегодня известно как географические информационные системы. Говард Фишер, вдохновитель лаборатории, разработал усовершенствованную программу компьютерного картографирования под названием SYMAP, и с самого начала стремился улучшить методы интерполяции. Он демонстрировал свою работу над SYMAP первокурсникам Гарвардского колледжа, и многие из них участвовали в мероприятиях лаборатории. Один из первокурсников, Дональд Шепард, решил полностью переработать интерполяцию в SYMAP, что привело к публикации его знаменитой статьи в 1968 году. На алгоритм Шепарда также повлиял теоретический подход Уильяма Варнца и других сотрудников лаборатории, занимавшихся пространственным анализом. Он провел ряд экспериментов с показателем степени расстояния, остановившись на значении, близком к модели гравитации (степень 2). Шепард реализовал не только базовое обратное взвешивание расстояния, но и предусмотрел возможность использования барьеров (проницаемых и непроницаемых) при интерполяции. В то же время над интерполяцией работали и другие исследовательские центры, в частности, Университет Канзаса с их программой SURFACE II. Однако функциональность SYMAP оставалась передовой, несмотря на то, что она была разработана студентом.
At the Harvard Laboratory for Computer Graphics and Spatial Analysis, beginning in 1965, a varied collection of scientists converged to rethink, among other things, what are now called geographic information systems. The motive force behind the Laboratory, Howard Fisher, conceived an improved computer mapping program that he called SYMAP, which, from the start, Fisher wanted to improve on the interpolation. He showed Harvard College freshmen his work on SYMAP, and many of them participated in Laboratory events. One freshman, Donald Shepard, decided to overhaul the interpolation in SYMAP, resulting in his famous article from 1968. Shepard's algorithm was also influenced by the theoretical approach of William Warntz and others at the Lab who worked with spatial analysis. He conducted a number of experiments with the exponent of distance, deciding on something closer to the gravity model (exponent of 2). Shepard implemented not just basic inverse distance weighting, but also allowed barriers (permeable and absolute) to interpolation. Other research centers were working on interpolation at this time, particularly University of Kansas and their SURFACE II program. Still, the features of SYMAP were state of the art, even though programmed by an undergraduate.
Модифицированный метод Шепарда
Другая модификация метода Шепарда вычисляет интерполированное значение, используя только ближайших соседей в пределах сферы радиуса R (вместо использования всей выборки). Веса в этом случае незначительно изменяются. В сочетании с быстрой структурой пространственного поиска (например, kd-деревом), это становится эффективным методом интерполяции со сложностью N log N, подходящим для задач большого масштаба.
Another modification of Shepard's method calculates interpolated value using only nearest neighbors within R sphere (instead of full sample). Weights are slightly modified in this case:
When combined with fast spatial search structure (like kd tree), it becomes efficient N log N interpolation method suitable for large scale problems.