Алгоритм пересчета значений карты в регулярную сетку
В программе Depth Imaging для расчета значений поля в узлах регулярной сети
применяется функция surface
свободного программного обеспечения GMT (Generic Mapping Tools).
http://www.soest.hawaii.edu/soest/gmt.html
разработанного Полем Весселем и Валтером Смитом
(Paul Wessel and
Walter H.F.Smith).
Алгоритм основан на итерационном решении уравнения
(1 - T) * L (L (z)) + T * L (z) = 0
где T - фактор регуляризации (изменяется от 0 до 1)
L - оператор-лапласиан.
Значение T = 0 дает решение "минимальной кривизны", которое является эквивалентным решениям в пакетах SuperMISP и ISM.
Минимальная кривизна может вызвать нежелательные осцилляции и ложные локальные максимумы или минимумы (Смотри Смит и Wessel,1990).
Вы можете увеличить T > 0, чтобы подавить эти эффекты.
Опыт свидетельствует о том, что T ~ 0.25 обычно хорош для потенциальных эксплуатационных данных (грави-магнитометрия).
Можно увеличить T (T ~ 0.35) для сильно изрезанных полей.
T = 1 дает гармоническую поверхность (никакие максимумы или минимумы невозможны кроме как в точках заданных изначальных (измеренных) значений поля z(x,y).
Ссылки:
Wessel, P., and W. H. F. Smith, 1995, The
Generic Mapping Tools (GMT) version 3.0
Technical Reference & Cookbook, SOEST/NOAA.
Wessel, P., and W. H. F. Smith, 1995, New
Version of the Generic Mapping Tools Released, EOS Trans. AGU, 76,
p. 329.
Wessel, P., and W. H. F. Smith, 1995, New
Version of the Generic Mapping Tools Released,
http://www.agu.org/eos_elec/95154e.html, Copyright
1995 by the American Geophysical Union.
Wessel, P., and W. H. F. Smith, 1991, Free
Software Helps Map and Display Data, EOS Trans. AGU, 72, p. 441.
Smith, W. H. F, and P. Wessel, 1990, Gridding with
continuous curvature splines in tension,
Geophysics, 55, 293-305.