Математическое моделирование уровенной поверхности океана по спутниковым данным на основе двумерной B-сплайн аппроксимации тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат технических наук Гомонов, Александр Дмитриевич

• Гомонов, Александр Дмитриевич
• кандидат технических наук
• 2011

Оглавление диссертации кандидат технических наук Гомонов, Александр Дмитриевич

ВВЕДЕНИЕ.

ГЛАВА I УРОВЕННАЯ ПОВЕРХНОСТЬ МИРОВОГО ОКЕАНА. ВОССТАНОВЛЕНИЕ УРОВЕННОЙ ПОВЕРХНОСТИ ОКЕАНА ПО

ДАННЫМ СПУТНИКОВЫХ НАБЛЮДЕНИЙ.

§1.1. Уровенная поверхность Мирового океана. Основные понятия предметной области.

§ 1.2. Влияние уровенной поверхности Мирового океана на формирование зон повышенной биопродуктивности.

§ 1.3. Временная устойчивость уровенной поверхности в заданных акваториях Мирового океана.

§ 1.4. Обзор методов восстановления уровенной поверхности Мирового океана по спутниковым данным.

ГЛАВА II ОБЩАЯ ПОСТАНОВКА ЗАДАЧИ ВОССТАНОВЛЕНИЯ АЛЬТИМЕТРИЧЕСКОЙ ПОВЕРХНОСТИ МИРОВОГО ОКЕАНА В

ЗАДАННОЙ АКВАТОРИИ.

§ 2.1. Особенности организации спутниковых альтиметрических наблюдений.

§ 2.2. Методы представления поверхностей.

§ 2.3. Постановка задачи восстановления уровенной поверхности Мирового океана в заданной акватории по данным спутниковых альтиметрических наблюдений.

§ 2.4. Обоснование выбора аппроксимационных функций для моделирования полей уровенной поверхности Мирового океана.

§ 2.5. Математическая модель уровенной поверхности Мирового океана в терминах Б-сплайнов. Обоснование критерия качества для поиска ее параметров.

ГЛАВА III ЧИСЛЕННЫЙ МЕТОД ВОССТАНОВЛЕНИЯ ПОЛЕЙ АЛЬТИМЕТРИИ В ЗАДАННОЙ АКВАТОРИИ МИРОВОГО ОКЕАНА ПО

ДАННЫМ СПУТНИКОВЫХ ИЗМЕРЕНИЙ.

§ 3.1. Минимизация целевого функционала как задача минимизации функционала квадратичной функции.

§ 3.2. Численный метод восстановления альтимстрических данных, нерегулярно представленных в плоской области. Особенности его реализации.

§ 3.3. Расчет течений по данным спутниковой альтиметрии.

ГЛАВА IV. ПРОГРАММНАЯ РЕАЛИЗАЦИЯ ЧИСЛЕННОГО МЕТОДА ВОССТАНОВЛЕНИЯ АЛЬТИМЕТРИЧЕСКИХ ПОЛЕЙ В ЗАДАННОЙ АКВАТОРИИ МИРОВОГО ОКЕАНА. РЕЗУЛЬТАТЫ ПРАКТИЧЕСКИХ

РАСЧЕТОВ.

§ 4.1. Описание программно-исследовательского комплекса.

§ 4.1.1. Структура, основные модули и функциональные возможности программно-исследовательского комплекса.

§ 4.1.2. Интерфейс программно-исследовательского комплекса.

§ 4.2. Результаты практических расчетов и их характеристика.

§ 4.2.1. Построение поверхностей по моделируемым данным некоторых математических функций и их характеристика.

§ 4.2.2. Построение уровенных поверхностей Мирового океана по спутниковым данным и их характеристика.

Введение диссертации (часть автореферата) на тему «Математическое моделирование уровенной поверхности океана по спутниковым данным на основе двумерной B-сплайн аппроксимации»

Актуальность

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

Подобного рода задачи могут возникать, например, в океанологии, метеорологии и т.п. В данной работе решается задача восстановления уровенной поверхности океана, отличающейся, как известно, временной изменчивостью, на основе накапливающихся во времени в заданной акватории данных спутниковых наблюдений. Несмотря на имеющиеся исследования, проводившиеся в этой связи в океанологии, эта задача сохраняет свою актуальность. Кроме этого она имеет достаточно важное практическое значение.

В ряде исследований [8, 20, 21, 46, 52, 53, 58, 63, 64, 78] было доказано, что существует однозначная взаимосвязь между значениями уровенной поверхность океана, атмосферным давлением, температурой, соленостью воды, направлениями и интенсивностями течений в океане. Следовательно, построив уровенные поверхности океана и отслеживая тенденции их изменений во времени, можно вырабатывать обоснованные прогнозы рыбопромысловой обстановки.

Поэтому в конечном итоге проблема определения перспективных районов промысла может быть решена посредством построения уровенных поверхностей океана и их аномалий (альтиметрии).

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

Развитие спутниковой альтиметрии в последнее десятилетие открыло более эффективные возможности в описании океанологических особенностей промысловых условий в различных регионах Мирового океана. Так, применяя методы дистанционного зондирования, можно оперативно получать достаточно большие массивы данных за сравнительно короткие промежутки времени практически по всему Мировому океану. Кроме того, по сравнению с традиционными контактными измерениями дистанционные методы зондирования поверхности океана являются более подробными в пространстве, но главное их преимущество — всепогодность, так как для зондирования океана используется микроволновый диапазон, в котором радиоволны проходят через облака с малыми затуханиями.

Основная характеристика динамики морской поверхности, которую измеряют при помощи спутниковых альтиметров, — аномалии высоты морской поверхности океана, то есть вычисляется аномалия силы тяжести по превышениям геоида над эллипсоидом. Эти работы начались в 1985 г. с запуска американского спутника GEO SAT. Позднее были запущены спутники ERS-1, ERS-2, TOPEX/POSEIDON, Jason-1, GFO и NVS. Самым последним был запущен совместно NASA и Францией 20 июня 2008 г. спутник OSTM/Jason-2 [http://www.nasa.gov/missionpages/ostm/news/ostmb20090805b. html].

Можно отметить и первые удачные опыты оперативного использования альтиметрических данных в реальном масштабе времени для решения промысловых задач в районах Центрально-Восточной Атлантики и Юго-Восточной Атлантики [46], а также на рыбных промыслах в Севоро-Восточной Атлантике и Баренцевом море [52, 53, 63, 78].

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

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

Необходимо также учитывать, что в этих данных присутствуют ошибки самих измерений и внесенными специализированными поправками к аль-тиметрическому сигналу. Поправки обусловлены тем, что сигнал от спутникового альтиметра проходит через атмосферу (атмосферные поправки) и отражается от мгновенного уровня морской поверхности (поправки на состояние и уровень подстилающей поверхности). Самая большая поправка связана с приливами. Существует несколько моделей приливных поправок: CSR 3.0, FES 95.2, GOT 99. Эти модели убирают приливную составляющую из альти-метрических данных, однако они не всегда хорошо работают в прибрежных зонах и в динамически активных районах Мирового океана [45].

Следует отметить, что на настоящий момент времени реализовано несколько математических моделей и программ, позволяющих строить уровен-ные поверхности океана [3, 48, 58]. Однако во всех них заложен слишком большой период осреднения данных (минимум 10 суток), что является абсолютно неприемлемым для определения быстро изменяющихся гидродинамических структур в океане. Кроме того, существующие модели, как правило, работают в диапазоне широт ± 72°, что тоже в свою очередь является существенным их недостатком, а главное - они не учитывают дополнительные требования, предъявляемые к моделируемым поверхностям. Например, скорость изменения поверхности в некоторых ее точках в заданном направлении.

Учитывая вышеизложенное, разработка математической модели для построения уровенной поверхности Мирового океана и их аномалий по спутниковым данным за более короткие периоды (менее 10 суток) в диапазоне широт ±81,5° с учетом дополнительных требований, предъявляемых к ней, а также разработка алгоритма и программной реализации являются весьма актуальными и многообещающими с практической точки зрения.

Предметом исследования данной работы является уровенная поверхность Мирового океана.

Методы исследования: математическое моделирование, численные методы решения задач оптимизации с применением ЭВМ, теория аппроксимации функций, вычислительный эксперимент.

Источники исследования: работы отечественных и зарубежных исследователей в области теории аппроксимации функций и решения задач оптимизации, а также спутниковые вдольтрековые альтиметрические данные.

Цель работы

Целями данной работы являются разработка математической модели уровенной поверхности океана методом аппроксимации спутниковых вдольтрековых альтиметрических данных двумерными ^-сплайнами и реализация на ее основе программных средств, устойчиво работающих в различных широтах.

Задачи исследования

1. Выполнение обзора существующих методов построения уровенной поверхности Мирового океана и их сравнительный анализ.

2. Разработка методики определения временного периода устойчивости уровенных поверхностей Мирового океана.

3. Построение математической модели уровенной поверхности Мирового океана.

4. Разработка метода подбора параметров математической модели уровенной поверхности в заданной акватории океана по спутниковым вдольтрековым альтиметрическим данным.

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

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

7. Построение карт уровенных поверхностей Мирового океана в заданных акваториях двумерными 5-сплайнами в диапазоне широт ±81,5°.

8. Проведение анализа адекватности моделируемых карт аномалий высоты поверхности океана реальным данным спутниковых измерений.

Научная новизна содержится в следующих результатах диссертационной работы:

1. Предложена математическая модель уровенной поверхности Мирового океана на базе двухмерных ^-сплайнов.

2. Предложен целевой функционал для определения параметров модели, учитывающий рассогласование моделируемой поверхности со спутниковыми вдольтрековыми альтиметрическими данными и априорную информацию о характере поведения восстанавливаемого поля.

3. Предложена методика определения временного периода устойчивости уровенных поверхностей Мирового океана.

4. Предложена методика определения минимально достаточного времени накопления вдольтрековых альтиметрических данных для обеспечения необходимой степени достоверности моделируемой поверхности.

Достоверность результатов исследования подтверждается результатами вычислительных экспериментов, осуществлявшихся для различных акваторий Мирового океана по спутниковым альтиметрическим данным.

Практическая значимость работы заключается в реализации предложенных методов и алгоритмов в виде программного комплекса с отображением результатов моделирования уровенной поверхности океана в электронной картографической системе. Разработанный программный комплекс позволяет осуществлять моделирования морских поверхностей в диапазоне широт ±81,5°. Кроме того, предложенная методика показала хорошие результаты для восстановления полей распределений температур, соленостей, хлорофилла и т.д. Результаты моделирования уровенной поверхности океана позволяют рассчитать в реальном масштабе времени направления морских течений, что реализовано как один из режимов работы программного комплекса. Это позволяет определить зоны повышенной биопродуктивности в океане. В целом результаты работы могут способствовать более полному пониманию физических процессов, происходящих в океане, что дает возможность формировать обоснованные рекомендации рыбопромысловому флоту о наиболее предпочтительных местах ведения промысла и может способствовать повышению эффективности его работы.

Положения, выносимые на защиту:

Методика компьютерного моделирования уровенной поверхности океана по спутниковым данным, включающая в себя:

- математическую модель уровенной поверхности Мирового океана на базе двухмерных 5-сплайнов;

- целевой функционал для определения параметров модели, учитывающий рассогласование моделируемой поверхности со спутниковыми вдольтрековыми альтиметрическими данными и априорную информацию о характере поведения восстанавливаемого поля;

- методику определения временного периода устойчивости уровенных поверхностей Мирового океана;

- методику определения минимально достаточного времени накопления вдольтрековых альтиметрических данных для обеспечения необходимой степени достоверности моделируемой поверхности;

- программный комплекс, реализующий компьютерное моделирование уровенной поверхности Мирового океана по данным спутниковой альтиметрии.

Апробация работы

По теме диссертации опубликовано 14 печатных работ, среди которых две публикации в ведущих рецензируемых изданиях, рекомендованных в действующем перечне ВАК, и две — в нерецензируемых изданиях. Доклады доложены и получили одобрение на восьми международных научно-практических и всероссийских научно-технических конференциях. Основные положения работы защищены двумя авторскими свидетельствами на изобретения.

Структура диссертации: диссертация состоит из введения, четырех глав, заключения, списка литературы и трех приложений.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Гомонов, Александр Дмитриевич

Выводы.

Результаты, приведенные в таблицах 4.1 и 4.2, показывают, что двумерные В -сплайны можно применять к аппроксимации гладких функций, в том числе заданных с некоторыми погрешностями, о чем свидетельствует уменьшение норм ошибок моделирования Кп и £2>п ■ Данное обстоятельство дает возможность применять методы двумерной В -сплайн аппроксимации для построения полей аномалий уровенной поверхности Мирового океана.

Кроме того, как следует с расчетов приведенных на графиках рисунков 4.19—4.21, определен минимально-достаточный временной интервал накопления вдольтрековых альтиметрических данных необходимый для качественного восстановления уровенной поверхности океана. С помощью предлагаемой математической модели это время определено по различным акваториям океана: 5 суток для акватории 560Ы-64°М и 22^-01^ (Северная Атлантика);

4 суток для акватории 10°К-25°К и 250\\г-15°\\'г (Северо-Западная Африка);

8 суток для акватории 70°1Ч-75оК и 20°Е-50°Е (Баренцево море). Как и следовало ожидать, продолжительность данного периода зависит от координат акватории, на которой планируется восстановление уровенной поверхности океана. Вблизи экватора плотность получаемых спутниковых данных выше, а ближе к полюсам уменьшается. Это связано с количеством спутников, которые «перекрывают» данную акваторию (рисунок 2.1).

114

ЗАКЛЮЧЕНИЕ

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

Рассматривая теорию аппроксимации функций, Р. В. Хемминг [62. с. 91] писал: «Прежде чем приступать к вычислениям, необходимо дать ответ на следующие вопросы:

1. Какие узлы мы будем использовать?

2. Какой класс приближающих функций мы будем использовать?

3. Какой критерий согласия мы примем?

4. Какую точность мы хотим?»

Рассмотрение уровенной поверхности океана и ее пространственно-временных изменений как модели физического явления позволило ответить на вопросы, поставленные Р. В. Хеммингом:

1. Для аппроксимации уровенной поверхности океана наиболее целесообразным является выбор узлов регулярной сети.

2. Использование в качестве класса приближающих функций кубических ¿¿-сплайнов дает ряд преимуществ перед другими классами функций.

3. В качестве критерия согласия был выбран функционал качества на базе метода наименьших квадратов, который позволил использовать избыI точную информацию, чтобы получить некоторое сглаживание неточности данных спутниковых измерений, дополненный требованиями к моделируемой поверхности океана.

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

В проведенном исследовании были получены следующие основные результаты:

1. Разработана математическая модель уровениой поверхности Мирового океана, позволяющая осуществлять моделирование поверхности по спутниковым вдольтрековым альтиметрическим данным с использованием двумерных £>-сплайнов.

2. Предложен целевой функционал для определения параметров модели, учитывающий рассогласование моделируемой поверхности со спутниковыми вдольтрековыми альтиметрическими данными и априорную информацию о характере поведения восстанавливаемого поля.

3. Предложена методика определения временных периодов устойчивости физических процессов, влияющих на формирование уровенной поверхности. Эти периоды определены для регионов Мирового океана:

-5-10 суток для акватории 56°N-64°N и 22°W-01°W (Северная Атлантика);

- 5-8 суток для акватории 10°N-25°N и 25°W-15°W (Северо-Западная

Африка);

- 5-10 суток для акватории 45°N-55°N и 145°Е-155°Е (Дальний Восток);

- 9-14 суток для акватории 70°N-75°N и 20°Е-50°Е (Баренцево море).

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

- 5 суток для акватории 56°N-64°N и 22°W-01°W (Северная Атлантика);

- 4 суток для акватории 10°N-25°N и 25°W-15°W (Северо-Западная Африка);

- 8 суток для акватории 70°N-75°N и 20°Е-50°Е (Баренцево море).

5. Разработан программный комплекс в среде Borland Delphi 7.0 с развитым человеко-машинным интерфейсом, реализующий предлагаемую методику построения уровенной поверхности океана и их аномалий по спутниковым альтиметрическим данным. Программный комплекс обеспечивает трехмерную визуализацию и отображение результатов вычислений на электронной картографической системе.

6. Проведен экспериментальный анализ согласованности моделируемых уровенных поверхностей океана с реальными данными спутниковых измерений.

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

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

Список литературы диссертационного исследования кандидат технических наук Гомонов, Александр Дмитриевич, 2011 год

1. Алберг, Дж. Теория сплайнов и ее приложения / Дж. Алберг, Э. Ниль-сон, Дж. Уолш. М. : Мир, 1972. - 316 с.

2. Алгоритмы и структуры данных геоинформационных систем / Крас-нояр. гос. техн. ун-т ; сост.: И. В. Варфоломеев, И. Г. Ермакова, И. С. Савельев. Красноярск : Изд-во КГТУ, 2003. - 34 с.

3. Белоненко, Т. В. Вихри или волны? / Т. В. Белоненко, Е. А. Захарчук,

4. B. Р. Фукс // Вестник С.-Петерб. ун-та. Сер. 7. 1998. - Вып. 3, № 21.1. C. 37-44.

5. Белоненко, Т. В. Градиентно-вихревые волны в океане / Т. В. Белоненко, Е. А. Захарчук, В. Р. Фукс. СПб. : Изд-во С.-Петерб. гос. ун-та, 2004.-215 с.

6. Белышев, А. П. Вероятностный анализ морских течений / А. П. Бе-лышев, Ю. П. Клеванцов, В. А. Рожков. Л. : Гидрометеоиздат, 1983. -264 с.

7. Березин, И. С. Методы вычислений / И. С. Березин, Н. П. Жидков. -М. : Наука, 1966. 632 с.

8. Блаттер, К. Вейвлет-анализ. Основы теории / К. Блаттер ; пер с нем. Т. Э. Кренкеля ; под ред. А. Г. Кюркчана. М. : Техносфера, 2006. -272 с. - (Мир математики).

9. Бор, К. де Практическое руководство по сплайнам / К. де Бор. — М. : Радио и связь, 1985. 304 с.

10. Василенко, В. А. Сплайн-функции : теория, алгоритмы, программы / В. А. Василенко. Новосибирск : Наука, 1983. - 224 с.

11. Васильев, Ф. П. Лекции по методам решения экстремальных задач / Ф. П. Васильев. М. : Изд-во Моск. ун-та, 1974. - 364 с.

12. Васильев, Ф. П. Методы оптимизации / Ф. П. Васильев. — М. : Факториал Пресс, 2002. 824 с.

13. Вершинин, В. В. Экстремальные свойства сплайнов и задача сглаживания / В. В. Вершинин, Ю. С. Завьялов, Н. Н. Павлов. Новосибирск : Наука, 1988. - 102 с.

14. Воеводин, В. В. Численные методы алгебры. Теория и алгоритмы / В. В. Воеводин. М. : Наука, 1966. - 248 с.

15. Геоинформационная система Golden Software Surfer 8 : учебно-метод. пособие для вузов / К. Ю. Силкин, И. В. Варфоломеев, И. Г. Ермакова, И. С. Савельев. Воронеж : Изд-во ВГУ, 2008. - 66 с.

16. Гребенников, А. И. Метод сплайнов в численном анализе / А. И. Гребенников. М.: Изд-во МГУ, 1979. - 208 с.

17. Демьянов, Г. В. Планетарные модели гравитационного поля Земли, их роль в современных условиях развития геодезии / Г. В. Демьянов, Р. А. Сермягин // Изв. вузов. Геодезия и картография / Моск. гос. ун-т геодезии и картографии. 2009. — № 10. - С. 8-12.

18. Дьяконов, В. П. Mathematica 4.1/4.2/5.0 в математических и научно-технических расчетах / В. П. Дьяконов. М. : COJIOH-Пресс, 2004. -696 с.

19. Завьялов, Ю. С. Интерполирование кубическими многозвенками / Ю. С. Завьялов // Вычислительные системы : сб. науч. тр. / АН СССР, Сиб. отд-ние, Ин-т математики. Новосибирск, 1970. — Т. 38. - С. 2373.

20. Завьялов, Ю. С. Методы сплайн-функций / Ю. С. Завьялов, Б. И. Квасов, В. Л. Мирошниченко. — Новосибирск : Наука, 1980. — 350 с.

21. Завьялов, Ю. С. Сплайны в инженерной геометрии / Ю. С. Завьялов, В. А. Леус, В. А. Скороспелов. М. : Машиностроение, 1985. - 224 с.

22. Завьялов, Ю. С. Экстремальное свойство кубических многозвенников и задача сглаживания / Ю. С. Завьялов // Вычислительные системы : сб. науч. тр. / АН СССР, Сиб. отд-ние, Ин-т математики. Новосибирск, 1970.-Т. 42.-С. 89-108.

23. Инженерная геодезия / Е. Б. Клюшин, М. И. Кисилев, Д. Ш. Михелев,

24. В. Д. Фельдман. М. : Академия, 2004. - 480 с.

25. Карманов, В. Г. Математическое программирование / В. Г. Карманов. -М. : Наука, 1986.-288 с.

26. Корнейчук, Н. П. Сплайны в теории приближения / Н. П. Корнейчук. -М. : Наука, 1984.-352 с.

27. Ландау, Л. Д. Теория упругости / Л. Д. Ландау, Е. М. Лифшиц. М. : Наука, 1965. -204 с. - (Теоретическая физика, т. 7).

28. Латфуллин, М. Г. Технология автоматизированного картирования геофизических параметров с учетом дизъюнктивных нарушений : авто-реф. дис. . канд. техн. наук : 04.00.12 / М. Г. Латфуллин. М., 1997. -21 с.

29. Мамаев, О. И. Морские течения / О. И. Мамаев. М. : Изд-во МГУ, 1986. - 103 с.

30. Моисеев, Н. Н. Методы оптимизации / Н. Н. Моисеев, Ю. П. Ивани-лов, Е. М. Столярова. М. : Наука, 1978. - 352с.

31. Мориц, Г. Современная физическая геодезия / Г. Мориц ; пер. с англ. П. П. Медведева. М. : Недра, 1983. - 392 с.

32. Нейман, Г. Н. Океанские течения / Г. Н. Нейман. М. : ГИМИЗ, 1973.-207 с.

33. Параметры Земли 1990 года (ПЗ-90.02). Параметры общеземного эллипсоида и гравитационного поля Земли (проект) / Военно-топограф. Упр. Генерального штаба. М., 2005. - 57 с.

34. Поклад, Г. Г. Геодезия / Г. Г. Поклад, С. П. Гриднев. М. : Академ, проект, 2007. - 592 с.

35. Пшеничный, Б. Н. Численные методы в экстремальных задачах / Б. Н. Пшеничный, Ю. М. Данилин. М. : Наука, 1975. - 311 с.

36. Романов, А. А. Восстановление мезомасштабной изменчивости аномалий высоты поверхности океана по данным спутниковой альтиметрии : дис. . канд. физико-мат. наук : 25.00.29. / А. А. Романов. М., 2004. - 125 с.

37. Романов, А. А. Методика восстановления карт аномалий морской поверхности методом Вт -сплайнов Электронный ресурс. / А. А. Романов // Исследовано в России. 2004. - № 42. - С. 454^-63. - Режим доступа: http://zlшrnal.ape■relam.ru/2004/042.pdf■

38. Рябенький, В. С. Введение в вычислительную математику: учеб. пособие / В. С. Рябенький. 2-е изд. - М. : ФИЗМАТЛИТ, 2000. - 296с.

39. Скворцов, А. В. Триангуляция Делоне и ее применение / А. В. Скворцов. Томск : Изд-во Том. ун-та, 2002. - 128 с.

40. Смоленцев, Н. К. Основы теории вейвлетов. Вейвлеты в МАТЬАВ / Н. К. Смоленцев. —М. : Пресс, 2008.-448 с.

41. Смоляк, С. А. Сплайны и их применение / С. А. Смоляк // Экономика и математические методы / АН СССР, Центр, экономико-мат. ин-т. -М., 1971.-Т. 7, вып. З.-С. 419—431.

42. Старицын, Д. К. Основы использования спутниковой альтиметриче-ской информации для оценки условий промысла сайры / Д. К. Старицын, В. Н. Филатов, В. Р. Фукс // Изв. ТИНРО / Тихоокеан. науч.-исслед. рыбохозяйств. центр. 2004. - Т. 137. - С. 398-408.

43. Стечкин, С. Б. Сплайны в вычислительной математике / С. Б. Стеч-кин, Ю. Н. Субботин. М. : Наука, 1976. - 248 с.

44. Субботин, Ю. Н. Экстремальные задачи функциональной интерполяции и интерполяционные в среднем сплайны / Ю. Н. Субботин // Тр. МИАН СССР / Мат. ин-т им. В. А. Стеклова РАН. 1975. - Т. 138. - С. 118-173.

45. Фиакко, А. Нелинейное программирование. Методы последовательной безусловной минимизации / А. Фиакко, Г. Мак-Кормик. — М. : Мир, 1972.-237с.

46. Хемминг, Р. В. Численные методы : для науч. работников и инженеров / Р. В. Хемминг ; пер. с англ. В. Л. Арлазарова, Г. С. Разиной, А. В. Усакова ; под. ред. Р. С. Гутера. М. : Наука, 1968. - 400 с.

47. Шериф, Р. Е. Англо-русский энциклопедический словарь терминов разведочной геофизики / Р. Е. Шериф ; пер. с англ. А. А. Богданова ; ред. пер. И. И. Гурвич. М. : Недра, 1984. - 351 с.

48. Bretheron, F. P.A Technique for Objective Analysis and Design of Ocea-nographic Experiments Applied to MODE-73 / F. P. Bretheron, R. E. Davis,

49. C. B. Fandry // Deep Sea Res. 1976. - № 23. - P. 559-582. (42)

50. Dcrber, J. A global oceanic data assimilation system / J. Derber, A. Rosati // J. Phys. Oceanogr. 1989. -№ 19. - P. 1333-1347. (50)

51. Ducet, N. Global high resolution mapping of ocean circulation from the combination of TOPEX/POSEIDON and ERS-1/2 / N. Ducet, P. Y. Le Traon, G. Reverdin // J. of Geophysical Res. 2000. - Vol. 105, № C8. -P.477-498. (32)

52. Field, D . A. Laplacian Smoothing And Delaunay Triangulations /

53. D. A. Field // Communications in Applied Numerical Methods. 1988. -Vol. 4.-P. 709-712. (59)

54. Forbes, C. Assimilation of sea surface height data into an isypicnic ocean model / C. Forbes, O. Brown // J. Phys. Oceanogr. 1996. - Vol. 26. -P. 1189-1213.(33)

55. Gandin, L. S. Objective analysis of meteorological fields / L. S. Gandin. -Leningrad : GIMIZ, 1963. 242 p. (41)

56. Global statistical analysis of TOPEX and Poseidon data / P. Y. Le Traon et all. // J. Geophys. Res. 1994. - Vol. 99, № C12. - P. 619-631.

57. Harder, R. L. Interpolation using surface splines / R. L. Harder, R. N. Desmarais // J. of Aircraft. 1972. - Vol. 9, № 2. - P. 189-191.

58. Le Traon, P. Y. ERS-1/2 orbit improvement using Topex/Poseidon: the 2 cm challenge / Le Traon P. Y., F. Ogor // J. Geophys. Res. 1998. -№95.-P. 8045-8057.

59. Le Traon, P. Y. Spatial scales of mesoscale variability in the Nort Atlantic as deduced from Geosat data / P. Y. Le Traon, M. C. Rouquet, C. Boisser // J. Gcophys. Res 1990. - Vol. 95. - P. 267-285.

60. Le Traon, P. Y. Use of a High-Resolution Model to Analyze the Mapping Capabilities of Multiple-Altimeter Missions / P. Y Le Traon, G. Dibarboure, N. Ducet // J. of Atmospheric and Oceanic Technology. 2001. - № 18. -P. 1277-1288.

61. Rogers, F. Mathematical elements for computer graphics (McGraw-Hill Science/Engineering/Math) / F. Rogers, J. Adams. 2 ed. - New York, 1989.-512 p.

62. Schoenberg, I. J. Contributions to problem of approximation of equidistant data by analytic function / I. J. Schoenberg // Quart. Appl. Math. 1946. -№ 4. - P. 45-99.

63. Schoenberg, I. J. Spline functions and the problem of graduation / I. G. Schoenberg // Proc. Nat. Acad. Sci. U. S. A. 1964. - Vol. 52, № 4. -P. 947-949.

64. Schoenberg, I. J. Spline interpolation and higher derivatives /1. J. Schoen-berg 11 Abhandlungen aus Zahlentheorie und Analysis / Deutscher der Wissenschaften ; ed. P. Turan. Berlin, 1968. - P. 279-295.

65. Shewchuk, J. R. Delaunay Refinement Algorithms for Triangular Mesh Generation / J. R. Shewchuk // Computational Geometry: Theory and Applications. 2002. - № 22 (1-3). - P. 21-74.

66. Short-arc analysis of intersatellite tracking data in a gravity mapping mission / D. D. Rowlands et all. // J. of Geodesy. 2002. - Vol. 76, № 6-7. -P. 307-316.

67. Stammer, D. Global characteristics of ocean variability estimated from regional Topex/Poseidon Altimeter Measurements / D. Stammer // J. Phys. Oceanogr. 1997. - № 27. - P. 1743-1769.

68. Stammer, D. Global characteristics of ocean variability estimated from regional Topex/Poseidon altimeter measurements / D. Stammer // J. Phys. Oceanogr/- 1997. № 27. - P. 1743-1769.

69. The Joint Gravity Model-3 / Tapley B. D. et al. // J. Geophys. Res. -1996. Vol. 101, №B12.- P. 29-49.

70. The midlatitude resolution capability of sea level fields constructed from single and multiple satellite altimeter datasets / Greenslade et all. // J. At-mos. Oceanic. Technol. 1997. - № 14. - P. 849-870.

71. Wunsch, C. Sampling characteristics of satellite orbits / C. Wunsch // J. Atmos. Oceanic Technol. 1989. - № 6. - P. 891-907.

72. СПИСОК И ДАННЫЕ СУЩЕСТВУЮЩИХ ЭЛЛИПСОИДОВ LIST OF DATUMS AND ELLIPSOIDS

73. No DATUM NAME ELLIPSOID NAME

74. EUROPEAN DATUM 1950 International Hayford 19242. WGS-66 NWL 8D/WGS 19663. WGS-72 NWL 10D/WGS 19724. WGS-84 WGS 19845. DUTCH DATUM Bessel6. NAD 1927 Clarke 1866

75. ORDINANCE SURVEY UK 36 Airy

76. ORDINANCE SURVEY UK 70 Airy

77. NAHRWAN DATUM 1976 Clarke 1880 Modified

78. MERCURY 1960 Fisher (Mercury) 1960

79. MERCURY 1968 Modified Mercury 1968 Modified

80. ARGENTINE International Hayford 1924

81. AMERICAN SAMOA 1962 Clarke 1866

82. ASCENSION ISL. 1958 International Hayford 1924

83. CHUA ASTRO-BRAZ. GEODETIC International Hayford 1924

84. CORREGO-ALLEGRE-br map International Hayford 1924

85. EASTER ISL. 1967 ASTRO International Hayford 1924

86. PROV. S. AMERICAN 56 International Hayford 1924

87. PROV. S. CHILE 1963 International Hayford 1924

88. SOUTH AMERICAN 1969 Reference Ellipsoid 196721. SAO-C5 SOA-C522. SAO-C6 SOA-C623. SAO-C7 SOA-C724. OLD BAVARIAN Bessel25. POTSDAM Bessel26. BERNE 1898 Bessel27. TOKYO Bessel28. BERMUDA 1957 Clarke 186629. OLD HAWAIIN Clarke 1866

89. S.E. ISLAND SOLOMONS Clarke 1880.

90. YOF ASTRO 1967-DAHAR Clarke 1880

91. ADINDAN -EPHIOPIA Clarke 1880 Modified33. INDIAN Everest 1830

92. PULKOVO 1942 Krassovsky 1932

93. AUSTRALIAN GEODETIC Reference Ellipsoid 1967

94. NEW ZEALAND 1949 International Hayford 1924

95. WAKE ISL. ASTRO 1952 International Hayford 1924

96. FITILEVU 1916-FIJI Clarke 1880

97. SOUTH GEORGIA ASTRO International Hayford 1924

98. LUZON 1911 -PHILIPPINES Clarke 1866

99. ISLA SOCORRO ASTRO Clarke 1866

100. No DATUM NAME ELLIPSOID NAME42. GUAM Clarke 1866

101. GACIOSO ISL.-AZORES International Hayford 1924

102. CANTON ASTRO 1966 International Hayford 1924

103. CAMP AREA 61-62 ASTRO International Hayford 1924

104. BETIO ISLAND 1966 International Hayford 1924

105. CHRISTMAS ISLAND ASTRO 1967 International Hayford 1924

106. EFTATE International Hayford 1924

107. TRISTAN ASTRO 196B International Hayford 1924

108. No ELLIPSOID NAME SEMI-MAJOR AXIS (a) SEMI-MINOR AXIS (b)

109. Airy 6 377 563.396 6 356 256.909

110. Bessel 6 377 97.155 6 356 078.963

111. Clarke 1858 6 378 93.645 6 356 617.938

112. Clarke 1866 6 378 206.400 6 356 583.800

113. Clarke 1880 6 378 301.000 6 356 566.918

114. Clarke 1880 Modified 6 378 249.145 6 356 514.870

115. Everest 1830 6 377 276.345 6 356 075.413

116. Fisher (Mcrcury) 1960 6 378 166.000 6 356 784.283

117. Fisher 1960 (S.Asia) 6 378 155.000 6 356 773.320

118. Fisher 1968 6 378 150.000 6 356 768.955

119. Hayford 1909 6 378 388.000 6 356 911.946

120. Helmert 1906 6 378 200.000 6 356 818.170

121. International Astro Union 1968 6 378 160.000 6 356 774.719

122. International Hayford 1924 6 378 388.000 6 356 911.946

123. Krassovsky 1932 6 378 245.000 6 356 863.019

124. Mercury 1968 Modified 6 378 150.000 6 356 768.337

125. NWL8D/WGS 1966 6 378 145.000 6 356 759.769

126. NWL 10DAVGS 1972 6 378 135.000 6 356 750.520

127. WGS-84 6 378 137.000 6 356 752.314

128. Reference Ellipsoid 1967 6 378 160.000 6 356 774.719

129. SOA-C5 6 378 165.000 6 356 779.702

130. SOA-C6 6 378 155.000 6 356 769.736

131. SOA-C7 6 378 142.000 . 6 356 757.1381. Flattening (f):f = ( a b )/aшетіши фвдіращііжшSш,