Бесплатный автореферат и диссертация по наукам о земле на тему

Разработка и исследование методики локального структурно-спектрального анализа оптических фотоизображений морской поверхности

ВАК РФ 25.00.34, Аэрокосмические исследования земли, фотограмметрия

Автореферат диссертации по теме "Разработка и исследование методики локального структурно-спектрального анализа оптических фотоизображений морской поверхности"

На правах рукописи

'' Сухорученко Александр Николаевич

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

25.00.34 — Аэрокосмические исследования Земли, фотограмметрия

Автореферат диссертации на соискание ученой степени кандидат технических наук

Работа выполнена на кафедре прикладной экологии Московского государственного университета геодезии и картографии.

Научный руководитель: доктор технических наук,

профессор Малинников В.А,

Официальные оппоненты: доктор технических наук,

профессор Ильин Ю.А.

доктор технических наук, Аржененко Н.И.

Ведущая организация: Всероссийский Научно-исследовательский

Институт по проблемам гражданской обороны и чрезвычайных ситуаций

Защита диссертации состоится «20» «декабря» 2006 г. в «~12» час, на заседании диссертационного совета 212.143.01 в Московском Государственном Университете геодезии и картографии (МИИГАиК), по адресу: 105064, Москва, К-64, Гороховский переулок, 4 (ауд. 321).

С диссертацией можно ознакомиться в библиотеке МИИГАиК.

Автореферат разослан «ноября» 2006 г. Ученый секретарь

диссертационного совета

Краснопевцев Б. В.

Общая характеристика работы

Актуальность проблемы. Многообразие форм и • сильная пространственно-временная изменчивость морской поверхности в значительной степени определяют динамику процессов взаимодействия в системе «атмосфера — морская поверхность». Знание пространственно-временных структур морской поверхности при различных метеоусловиях требуется при разработке схем параметризации процессов взаимодействия в численных динамических моделях общей циркуляции атмосферы и океана, долгосрочного прогноза погоды, решения практических задач надводного и подводного мореплавания, океанского хозяйства, предотвращения загрязнения океана.

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

Кроме того, следует отметить, что в настоящее время отсутствует методика локального структурно-спектрального анализа оптических

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

Актуальность диссертационной работы, таким образом, обусловлена:

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

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

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

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

1. Проведение натурных экспериментов по крупномасштабной фотосъемке различных явлений на морской поверхности при различных метеоусловиях.

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

3. Исследование фрактальных и мультифрактальных характеристик основных типов природных явлений на морской поверхности.

4. Исследование пространственной структуры пригребневых зон и зон обрушения морских гравитационных волн.

4

Научная новизна работы заключается в следующем:

1. Рассмотрены теоретико-методические аспекты применения спектрального, линеаментного, вейвлет и мультифрактального анализа изображений и впервые сформулирована концепция локального структурно-спектрального анализа оптических изображений морской поверхности;

2. Разработана методика и программное обеспечение локального структурно-спектрального анализа оптических изображений морской поверхности;

3. Разработан мультифрактальный метод оценки функции распределения уклонов морской поверхности;

4. Впервые исследованы спектральные и мультифрактальные характеристики пригребневых зон и зон обрушения гравитационных волн.

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

Достоверность результатов подтверждается:

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

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

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

На защиту выносятся следующие разработки и результаты:

1. Методика локального структурно-спектрального (ЛСС) анализа оптических изображений морской поверхности.

2. Алгоритмы и программное обеспечение, реализующие предложенную методику ЛСС-анализа.

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

4. Результаты натурных экспериментов и практического использования методики ЛСС-анализа для изучения особенностей пригребневых зон и зон обрушения гравитационных волн.

Апробация работы. Основные положения диссертационной работы докладывались и обсуждались: на 3 международной конференции «Современные проблемы оптики естественных вод» (Санкт-Петербург, сентябрь, 2005), на 9 международной конференции «Методы дистанционного зондирования и ГИС-технологии для оценки состояния окружающей среды» (Италия, май, 2005), на 59 (апрель, 2004), и 60 (апрель, 2005) научно-технических конференциях студентов, аспирантов и молодых ученых, проводившихся в Московском Государственном университете геодезии и картографии (МИИГАиК).

Публикации; по результатам проведенных исследований и разработок, выполненных в процессе работы, опубликовано 5 научных работ.

Структура II объем работы: диссертация состоит из введения, четырех глав, заключения, списка литературы и приложения. Ее объем составляет 147 страниц текста, включая 22 рисунка и 5 таблиц (не входящих в приложение). Приложение содержит 2 таблицыи 20 рисунков. Список литературы содержит 98 наименований, из них 25 на иностранных языках.

СОДЕРЖАНИЕ РАБОТЫ Во введении определена актуальность проблемы корректного учета стохастической геометрии взволнованной морской поверхности в динамических моделях океана и процесах взаимодействия в системе океан-атмосфера и показана необходимость проведения локального структурно-спектрального анализа оптических изображений морской поверхности, сформулированы цель и задачи диссертационной работы.

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

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

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

Во второй главе «Методика локального структурно-спектрального анализа оптических изображений морской поверхности», изложена разработанная автором, автоматизированная методика локального структурно-спектрального (ЛСС) анализа оптических фотоизображений морской поверхности.

В начале дается краткое теоретико-методическое описание линеаментного, Фурье и ввейвлет, фрактального и мультифрактального анализа

изображений и временных рядов. Затем приводится подробное описание разработанной автором методики ЛСС-анализа.

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

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

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

оцифровка исходных фотоснимков морской поверхности', оцифровка волновых регистрограмм, выполняется с помощью специального программного модуля, разработанного в УФЦ ФПК МИИГАнК при участии автора;

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

бинаризация исходного полутонового изображения по заданному уровню серого тона.

3 этап. Спектральный анализ цифровых изображений. Хотя, существует множество программных пакетов, в которых можно выполнить преобразование Фурье цифрового изображения, нами был использован специальный программный модуль «Fourier», разработанный в МИИГАиК. Данный программный модуль позволяет не только выполнить стандартное одно- и двумерное преобразования Фурье, но и вычислить его следующие информативные признаки: геометрические; топологические; энергетические;

8

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

4 этап. Линеаментный анализ цифрового изображения.

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

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

При выполнении линеаментного анализа мы использовать стандартный программный пакет «ЬЕББА», разработанный в ПГС «Аэрогеология».

5 этап. Вейвлет — анализ. Для выполнения вейвлет анализа нами использовались математические пакеты Матлаб и Маткад.

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

1. Определение интервалов самоподобия исследуемой пространственной структуры.

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

массового радиуса) равен 10Д. Максимальный размер покрывающего квадрата - не более 1/5 Nmax *Л, где Nmax - наибольший размер изображения, для покрывающего круга - не более Nmin/2, Nmi„- нименьший размер изображения, 2. Определение фрактальной размерностей исследуемого цифрового изображения и параметра Херста (Н) одномерных пространственных сечений изображения различного азимута простирания. Оценку фрактальных свойств временных и пространственных последовательностей мы получаем с помощью показателя Н. Для вычисления Н используется зависимость "дисперсия — пространство (время)". Для этого исходная последовательность {Хк} разбивается на последовательные блоки размера т. Для полученной последовательности Х(ш) вычисляется значение s2(m)=M{(X(m)k- ш )2}и строится зависимость Log (s2(m)) как функция от Log (m). Усредненный угловой наклон b этой зависимости связан с показателем Херста соотношением b =2(Н- 1).

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

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

Бинарное изображение. Исходное изображение разбивается равномерной квадратной сеткой на элементарные ячейки. Ячейкам сетки, покрывающим исследуемый объект, присваиваются значения 1, а ячейкам, не покрывающим объект - значения 0, и формируется двумерная матрица нулей и единиц. Затем проводят ряд последующих разбиений на более крупные ячейки размером г j * г j с построением для каждого разбиения характеристической меры в виде равноячеечного распределения единиц pi . Здесь pi = N, / Е Nj , где N; -количество единиц в i-ой крупной ячейке, £ N, - общее количество единиц в матрице крупных ячеек, i = 1,2,3.....N.

Полутоновое изображение. Уровни серого тона принимают значения от 0 до 255. Исходное изображение разбивается равномерной квадратной сеткой на элементарные ячейки. Каждой элементарной ячейки — пикселу соответствует число - определенный уровень серого тона. Если каждое такое число разделить на сумму всех чисел на изображении, то получим меру для каждой элементарной ячейки. Затем на основе этой меры можно сгенерировать меры для разбиения (покрытия) изображения ячейками более крупного размера. Меру для каждой ячейки более крупного размера получим просто складывая меры отдельных элементарных ячеек в укрупненных ячейках. Таким образом, получаем набор равноячеечных разбиений (покрытий) с размерами ячеек от г0

до г п и мерами {pi , i=l......п}, определенными по заданным мерам

элементарных ячеек.

2.Расчет обобгценных корреляционных функций. Для каждого q из интервала [-Q, Q] расчет обобщенных корреляционных функций осуществляется по формуле x(q) = (pi)q~rT(4>, г—>0, здесь г - размер ячейки, т - набор непрерывных вещественных параметров, -со <т< оо.

3. Вычисление статистических характеристик мулътифрактальных спектров.. Алгоритм расчета спектров состоит в вычислении сумм:

Ak (q) = Sj [р 1 ik (q)( In pik (q))], Fk(q) = ZL [Pük (q)( In Piik (q))], put (q) =(Pi04 / Xik(q), Xik(q)= £i[(Pik)q]

для каждого разбиения на ячейки размером г k . Затем методом наименьших квадратов определяются сингулярные размерности a(q) и мультифрактальный спектр f(q): a(q) = ЭАк (q)/öln rk , f(q) = SFk/31n rk. По рассчитанным величинам a(q) и f(q) определяется спектр размерностей Репьи: t(q) = qa(q) - f(q), D4 = t(q)/(q-l),D,=a(q=l).

4. Анализ корректности вычисленных мулътифрактальных спектров.

Корректность мультифрактальных спектров определяется следующими критериями:

Оя, £ Оч2, при < ча, ((а(я=0)) = тах = О0, ((а(я=1)) = а(я=1) = О,, ((а(я,)) ^ ((а(я2 )), а(яО < а(я2), при 0<я,<я2 ,и другие. Заметим, что диапазон изменения а всегда конечен, а его величину К = атах -«ты = О.» — Би можно принять в качестве положительной меры хаоса в изучаемой системе (К=0 в случае белого шума при генерации меры по Фурье спектру).

Для выполнения мультифрактального анализа может быть использован программный модуль Ргас1аЬ, входящий в программный пакет МАТЪАВ, а также разработанный в МИИГАиК программный пакет «Фрактал-ПК».

8. Комплексированиерезультатов обработки изображений.

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

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

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

12

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

Линеаменты на водной поверхности, отражают объективную информацию о разномасштабных волновых процессов. Закономерная ориентированность линейных структур, проявляющихся на изображениях морской поверхности в виде линеаментов, отражает гребни морских волн или другие линейные структуры (например, сликовые полосы). Совокупность линеаментов, их транзитность, разноранговость и организованность в зоны используются в качестве доказательства не случайной, а физически обусловленной их природы. Это обстоятельство позволяет рассматривать сеть линеаментов как отражение определенной структурной организации динамических явлений на морской поверхности. Отмеченные закономерности могут использоваться при изучении пространственного распределения гребней гравитационных и капиллярно-гравитационных волн, при изучении межволновых взаимодействий и др. После обработки схем линеаментов строятся карты плотности развития линеаментной сети по выделенным направлениям. Сопоставление этих карт с результатами вейвлет анализа позволяет уточнить физическую природу линеаментов и более корректно интерпретировать роз-диаграммы линеаментных полей.

Разработанная методика ЛСС-анализа изображений была в дальнейшем использована при изучении особенностей пригребневых зон гравитационных волн.

В третьей главе «Исследование фрактальных и мультифрактальных характеристик изображений взволнованной морской поверхности»

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

морской поверхности представлена рис.:/ ( 2 столбец). В 3 столбце рис.1 представлены результаты статистического анализа наклонов фрактальной кривой при различных интервалах значений размера покрывающего элемента, использованных при линейной аппроксимации фрактальной кривой.

бинаризированное изображение фрактальной регрессионная кривая Статистика наклонов

5 50000 >20000 • юооо ' 500D 1 2000 I. 1000 500200 100 50 20 , 10 5 »..-' 2 1 Slope[Boxl: ...... - ■ . -- _ I 1 ¿1 < : '.y.G 077« К X X -у 1 * * - N. ■ ••• \ \. . •• Л % 4 N "ч 1 .2 5 10 20 50.100 200 . 500 1 Pious

i»« IW || :ад. «ю Г. !*« гз п7 j s ; tmltoCrfnrl , 1 ! m I, 1 1m 14 1 ID.». ' j 1 »» j M ДО ! 1 : Ы» S ш { Cw Ш.т j c, l ■ • i I "K LUlli^V—— 4 < 4 1 J Л »1 ---------------, P.MfJ , »1---------— ОШЫтШф ¡¡¡t IW» <!_ С» (¡ом 'Sm [«мм*»

Рис.1. Пример фрактального анализа бинаризированного изображени.

Как следует из результатов фрактального анализа для всех исследуемых изображений можно выделить два четко выраженных интервала значений

величины г покрытий исследуемых структур (г = 2,4,6,8,........160), для которых

характерно наличие линейной зависимости. Первый интервал: 4 — 20, второй 28 -64. Хотя для различных изображений границы этих интервалов немного смещаются. Следует отметить, этим интервалам можно поставить в соответствие определенные диапазоны пространственных частот и тем самым определить частотные интервалы, в которых исследуемые цифровые изображения обладают свойством самоподобия. Так, если сопоставить полученные результаты фрактального анализа с результатами спектральной обработки цифровых изображений, то легко увидеть, что выделенные в результате фрактального анализа частотные интервалы самоподобия, хорошо согласуются с частотными интервалами для которых зависимость (lnS(k), Ink) линейна. Наблюдающиеся расхождения между частотными интервалами зон самоподобия, определенных спектральным и фрактальным методами

14

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

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

1) отрезок прямой параллельный оси абсцисс,

2) линейную функцию, возрастающую с увеличением масштаба,

3) возрастающую гладкую кривую п-го порядка.

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

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

1. С увеличением скорости ветра в приводном слое атмосферы фрактальная размерность бликового поля уменьшается.

2. Число интервалов самоподобня, при увеличении скорости ветра увеличивается.

Мультифрактальпая параметризация изображений пространственных структур на морской поверхности состоит в определении некоторых

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

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

В этом методе разбиение исследуемой структуры квадратами размером е проводится таким же образом, как и в стандартной методике для бинаризированных изображений, но мера строится как отношение локальной дисперсии к глобальной. Для так или иначе построенной меры строится обобщенная корреляционная функция Z. В случае сингулярного поведения меры lnZ(q) ~ x(q) Ine. Зависимость (InZ, Ine) должна иметь линейный характер. Отсутствие линейной зависимости ставит под сомнение правомочность использования мультифрактального подхода для параметризации исследуемых структур.

Предварительно, для изучаемых изображений мы провели исследование зависимости обобщенной фрактальной меры от размера покрытия. Анализ полученных данных показал, что в целом для размеров покрытия 4, 8, 16, 32, 64, 128 и 256 пикселов имеет место линейная зависимость. Однако для малых масштабов (размеры покрытия 2 и 4 пикселя) линейность нарушается.

Таблица 1. Мультифрактальные параметры рассчитанные для изучаемых изображений морской поверхности.

Параметры Гравитационные волны Гравитационно- капиллярные волны Бликовые образования 11епные образования

Do 1.999 1.999 1.981 2.0

D, 1.993 1.979 1.926 1.984

D2 1.989 1.971 1.918 1.977

D-40 2.882 3.346 3.69 2.236

D4O 1.926 1.911 1.865 1.877

Д = D| - D4O 0.067 0.088 0.116 0.123

f-40 -0.37 -0.331 0.003 0.347

Gmax 2.961 3.587 4.221 2.287

Ctmin 1.895 1.893 1.852 1.841

5 — Ömax-Cimm 1.066 1.694 2.369 0.446

Поэтому, в дальнейшем все вычисления проводились для масштабов покрытия: 8*8, 16*16, 32*32, 64*64, 128*128, 256*256. На рис.2 представлены

17

типичные мультифрактальные характеристики для изучаемых изображений морской поверхности. Из представленных на рис.2 результатов мультифрактального анализа видно, мультифрактальные спектры различных пространственных образований на морской поверхности существенно отличаются друг от друга. Так, кривая мультифрактального спектра для изображения блика достигает оси абсцисс, что свидетельствует о значительно меньшем разнообразии в организации пространственных элементов изображения (как по размерам, так и по плотности фототона), чем для других изображений. Наиболее заметно отличия проявляются не в форме фрактальных кривых, а в величинах мультифрактальных характеристик (табл.1). Из данных, представленных в табл.1 видно, исследуемые изображения значительно отличаются друг от друга по однородности и упорядоченности. Известно, что количественными характеристиками однородности и упорядоченности являются мультифрактальные параметры Ею и А = D| - D10 ■ Показатель А отражает степень упорядоченности и нарушения симметрии для общей конфигурации исследуемой структуры. Уменьшение для исследуемой серии структур показывает, что в структуре становиться больше периодической составляющей. Это хорошо согласуется с полученными нами данными. Так для изображения G (система гравитационных волн) значение параметра А максимально, а для изображения В (пенное образование) — минимально. Значения параметра S = amax - amin иногда принимают в качестве количественной меры хаоса в изучаемой пространственной структуре (5 = 0 в случае белого шума при генерации меры по. Фурье спектру). Рассчитанные значения параметра 5 для рассматриваемых изображений (см. табл.1) вполне адекватно отражают их хаотическую составляющую.

Функция f(a) играет весьма важную роль в физике моря. С одной стороны, она представляет спектр фрактальных размерностей мультифрактала, каковым является взволнованная морская поверхность. С другой стороны, эта величина представляет статистику значений уклонов взволнованной морской поверхности в неупорядоченной системе. Действительно, примем для простоты

размер ячейки е=1.Тогда вероятность нахождения частицы в ячейке с номером i определяется формулой p¡ ~ (1/L)"', где L — размер образца в единицах е. Эту же формулу можно переписать и по-другому a¡ ~ ln(p¡)/ln(L). Отождествим вероятность p¡, используемую для расчета обобщенной статистической суммы Z(q,l), со значением квадрата модуля наклона морской поверхности |lI'(r¡)|2 (в случае использования дисперсионного метода расчетов фрактальной меры) в некоторой точке r¡. Вероятность заполнения ячейки p¡ представляет собой случайную величину. Такой же случайной величиной является, очевидно, и показатель степени a¡. Однако распределение значений а для мультифрактала известно. Оно определяется формулой n(a) ~ (1/L) f(a) =exp[f(a)lriL]. Как мы знаем, вблизи своего максимума при некотором значении ао функция f(ct) может быть аппроксимирована параболой f(a) =d-^(a-a0)2, где С, = f(a0)"/2. Поэтому при больших L распределение п(а) принимает вид n(a) ~ ехр [-^lnL(a-ao)2). Учитывая, что а — — lnp/lnL, получаем отсюда функцию распределения вероятностей p¡ : Ñ(p) ~ ехр [-¿¡lnL(lnp/lnL -t- а0)2]. Это есть так называемое лог-нормальное распределение. Поскольку вероятности p¡ характеризуют величину |lF(r¡)|2, мы приходим к результату, что квадрат модуля уклона распределен лог-нормально по отображаемому на изображении участку морской поверхности.

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

В четвертой главе «Исследование пространственной структуры прнгребневых зон и зон обрушения морских гравитационных волн» приводятся результаты экспериментальных исследований особенностей пространствепно-частотных характеристик прнгребневых зон и зон обрушения гравитационных

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

Рис. 3. Поле линеаментов изображения морской поверхности.

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

вытянутым лепестком диаграммы.

Проведенные исследования показали, что поле изолиний и поле градиентов плотности линеаментов четко локализуют системы ГК-волн, которые при малых скоростях ветра генерируются на подветренном сколне, а

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

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

Рис. 4. Вейвлет-спектрограмма сечения изображения морской поверхности.

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

преобразования как плавное изменение яркости вдоль пространственной оси одновременно (но с разной интенсивностью) на всех масштабах. Если тренд нарастающий, то яркость увеличивается, если убывающий - уменьшается. Наличие нескольких последовательных трендов представляется в виде поочередного нарастания и убывания яркости точек на спектрограмме, при этом точка изменения тренда имеет наибольшую или наименьшую яркость. Имеющие место локальные особенности спектрограммы (например, обусловленные пенными образованиями) представляются как линии резкого перепада яркости, исходящие из точки, соответствующей времени возникновения скачка яркости.

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

Нами был разработан специальный программный модуль, позволяющий по полю вейвлет-коэффициентов рассчитать мультифрактальную меру, по которой оцениваются различные мультифрактальные характеристики исследуемых структур. Для распознавания различных явлений на морской поверхности мы предлагаем использовать: информационную фрактальную размерность - D(l); ширину фрактального спектра - 5= amax - amin; размах

фрактального спектра — 5Г= — скейлинговый показатель спектральной плотности — р = 2+т2; показатель однородности структуры - показатель упорядоченности структуры - А = О] - Ц,,.

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

Заключение

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

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

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

2. Разработана автоматизированная технология оцифровки бумажных регистрограмм возвышений морской поверхности.

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

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

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

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

на морской поверхности..

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

По теме диссертации опубликованы следующие работы:

1. Малинников В.А., Сухорученко А.Н. Оценка фрактальной размерности оптических изображений морской поверхности. Труды Международной научно-технической конференции, посвященной 225-летию МИИГАиК, М„ 2004, с. 252-259.

2. Малинников В.А.,Сухорученко А.Н. Мультифрактальная параметризация оптических изображений морской поверхности. Сб. трудов IX Международной научно-практической конференции «Методы дистанционного зондирования и ГИС-технологии для оценки состояния окружающей среды», Италия, 2005, с. 19-23.

3. Сухорученко А.Н. Исследование фрактальных свойств поля блика на морской поверхности. Сб. научных трудов VIII межвузовского научно-практического семинара «Экологическая безопасность и рациональное природопользование», М., 2004, с. 139-142.

4. Сухорученко А.Н. Методика локального структурно-спектрального анализа оптических изображений морской поверхности. Изв. Вузов. Геодезия и аэрофотосъемка. Специальный выпуск, 2006, с. 131-140.

5. Сухорученко А.Н. Фрактальные свойства зон обрушений морских гравитационных волн. Сб. научных трудов VIII межвузовского научно-практического семинара «Экологическая безопасность и рациональное природопользование», М., 2004, с. 143-147.

МГУГиК

105064, Москва К-64, Гороховский пер., 4

Подп. к печати 16.11.2006 Формат 60x90

Бумага офсетная Печ. л. 5 Уч.-изд. л. 5 Тираж 80 экз. Заказ №253 Цена договорная

Содержание диссертации, кандидата технических наук, Сухорученко, Александр Николаевич

Введение.

Глава 1. Тематическая обработка оптических изображений морской поверхности.

1.1. Явления на морской поверхности. Общее описание.

1.2. . Методы контроля состояния морской поверхности.

1.3. Тематическая обработка оптических изображений морской поверхности.

1.4. Обзор современных исследований по изучению фрактальных свойств явлений на морской поверхности.

1.5. Постановка цели и задач исследования.

Глава 2. Методика локального структурно-спектрального анализа оптических изображений морской поверхности.

2.1. Теоретические основы вейвлет и фрактального анализа.

2.1.1. Непрерывное вейвлет-преобразование.

2.1.2. Мультифрактальный анализ.

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

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

Глава 3. Исследование фрактальных и мультифрактальных характеристик изображений взволнованной морской поверхности.

3.1. Фрактальные свойства изображений взволнованной морской поверхности.

3.1.1. Фрактальные свойства изображений взволнованной морской поверхности вне зоны блика.

3.1.2. Исследование фрактальных свойств взволнованной морской поверхности в зоне блика.

3.2. Мультифрактальная параметризация явлений на морской поверхности.

3.3. Взаимосвязь спектральных и фрактальных характеристик изображений морской поверхности.

3.4. Мультифрактальная параметризация шкалы Бофорта.

Глава 4. Исследование пространственной структуры пригребневых зон и зон обрушения морских гравитационных волн.

4.2. Фурье и вейвлет анализ волнового поля.

4.3. Применение вейвлет анализа для оценки фрактальных характеристик пригребневых зон.

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

Многообразие форм и сильная пространственно-временная изменчивость морской поверхности в значительной степени определяют динамику процессов взаимодействия в системе «атмосфера - морская поверхность». Знание пространственно-временных структуры морской поверхности при различных метеоусловиях требуется при разработке схем параметризации процессов взаимодействия в численных динамических моделях общей циркуляции атмосферы и океана, долгосрочного прогноза погоды, решения практических задач надводного и подводного мореплавания, океанского хозяйства, предотвращения загрязнения океана[52,91,93,94].

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

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

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

Актуальность диссертационной работы, таким образом, обусловлена:

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

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

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

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

1. Проведение натурных экспериментов по крупномасштабной фотосъемке различных явлений на морской поверхности при различных метеоусловиях.

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

3. Исследование фрактальных и мультифрактальных характеристик основных типов природных явлений на морской поверхности.

4. Исследование пространственной структуры пригребневых зон и зон обрушения морских гравитационных волн.

Научная новизна работы заключается в следующем:

1. Рассмотрены теоретико-методические аспекты применения спектрального, линеаментного, вейвлет- и мультифрактального анализа изображений и впервые сформулирована концепция локального структурно-спектрального анализа оптических фотоизображений морской поверхности;

2. Разработана методика локального структурно-спектрального анализа оптических фотоизображений морской поверхности;

3. Разработан мультифрактальный метод оценки функции распределения уклонов морской поверхности;

4. Впервые исследованы спектральные и мультифрактальные характеристики пригребневых зон и зон обрушения гравитационных волн.

Апробация работы. Основные положения диссертационной работы докладывались и обсуждались: на 3 международной конференции «Современные проблемы оптики естественных вод» (Санкт-Петербург, сентябрь, 2005), на 9 международной конференции «Методы дистанционного зондирования и ГИС-технологии для оценки состояния окружающей среды»

Италия, май, 2005), на 59 (апрель, 2004) и 60 (апрель, 2005) научно-технических конференциях студентов, аспирантов и молодых ученых, проводившихся в Московском Государственном университете геодезии и картографии (МИИГАиК).

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

Структура и объем работы: диссертация состоит из введения, четырех глав, заключения, списка литературы и приложения. Ее объехм составляет страниц текста, включая рисунков и таблиц (не входящих в

Заключение Диссертация по теме "Аэрокосмические исследования земли, фотограмметрия", Сухорученко, Александр Николаевич

Основные результаты и выводы работы заключаются в следующем.

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

2. Разработана автоматизированная технология оцифровки бумажных регистрограмм возвышений морской поверхности.

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

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

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

6. Исследованы вейвлет-спектрограммы изображений взволнованной морской поверхностию Показано, что они могут быть использованы для локализации и описания различных динамических явлений на морской поверхности.

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

Заключение

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

Библиография Диссертация по наукам о земле, кандидата технических наук, Сухорученко, Александр Николаевич, Москва

1. Автоматизированный линеаментный анализ природных линеаментных систем,-Л.: ВСЕГЕИ, 1988. 131 с.

2. Автоматизированный линеаментный анализ при структурно-геологических и металлогенических исследованиях. М., Недра, 1988, 86 с.

3. Алмазов И.В и др. Выбор оптимального варианта выполнения крупномасштабной аэрофотосъемки. Известия высших учебных заведений. Геодезия и аэрофотосъемка. М., 1995, №3, с. 86-92.

4. Алтынов А.Е., Манченко Д.А. Выбор оптимального масштаба фотографирования для крупномасштабной аэрофотосъемки. Известия вузов. Геодезия и аэрофотосъемка. М., 1987, №3, с. 71-76.

5. Астафьева Н.М. Вейвлет анализ: основы теории и примеры применения. УФНД996, т 166, №11.

6. Бадалян С. Обработка данных дистанционного зондирования с применением программы ER Mapper. Геоинформационные системы. №5, 1999.

7. Баранов ЮБ, Королев ЮК, Миллер СА. Программное обеспечение для обработки данных дистанционного зондирования. Информационный бюллетень. ГИС Ассоциация. №2 (9) с.42-45; №4 (11), 1997, с.40-47.

8. Божокин С.В., Паршин Д.А. Фракталы и Мультифракталы. Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001.

9. Бондур В.Г. Дистанционная оценка состояния границы раздела атмосфера-океан по пространственным спектрам изображений. В книге «Оптико-метеорологические исследования земной атмосферы». Новосибирск, Наука, 1987, с. 195-207.

10. Ю.Бондур В.Г. Оперативная дистанционная оценка состояния границы раздела атмосфера океан по пространственным спектрам изображений. Оптико-метеорологические исследования земной атмосферы, Новосибирск: Наука. Сиб. отделение, 1987. С. 217-230.

11. П.Бондур В.Г., Борисов Б.Д., Гении В.Н. и др. Поле яркости морской поверхности при искусственном импульсном освещении. Перенос изображений в земной атмосфере. Томск, 1988. С. 42-45.

12. Бондур В.Г., Мурынин А.Б, Восстановление пространственных спектров морской поверхности по оптическим изображениям с учетом нелинейной модуляции поля яркости. Оптика атмосферы. 1991. Т. 4. № 4. С. 387-395.

13. Бондур В.Г., Старченков С.А. Методы и программы обработки и классификации аэрокосмических изображений. Изв. Вузов. Серия Геодезия и аэрофотосъемка. 2001, № 1, с. 118-146.

14. Н.Бондур В.Г., Шарков Е.А. Статистические характеристики линейной геометрии пенных структур на поверхности моря по данным оптического зондирования. Исследования Земли из космоса. 1986. № 4. с. 21-31.

15. Бортковский Р.С. Пространственно-временные характеристики барашков и пятен пены, образующихся при обрушении. Метрология и гидрология. 1987. № 5. с. 68-75.

16. Волков Ю. А. Анализ спектров морского волнения, развивающегося под действием турбулентного ветра. Изв. АН СССР. Физика атмосферы и океана, 1968, т. 4, № 9, .с. 968-987.

17. Встовский Г. В., Колмаков А. Г., Бунин И. Ж. Введение в мультифрактальную параметризацию структур материалов. Москва-Ижевск: Научно-издательский центр «Регулярная и хаотическая динамика», 2001, 116 с.

18. Гиляров B.JI. и др. Применение вейвлет-преобразования при изучении изменения фрактальных свойств поверхностей аморфных металлов подвоздействием механической нагрузки. Физика твердого тела, 2004, т. 46, вып. 10, с.1806-1810.

19. Гришин Г.А. и др. Фурье-анализ изображений морской поверхности для моделирования эволюции поверхностного волнения. В сб. "Проблемы исследования океана из космоса ", 1984, с.64-68.

20. Давидан И.Н., Лопатухин Л.И„ Рожков В.А. Ветровое волнение в Мировом океане. Л.: Гидрометеоиздат, 1985. 256 с.

21. Дженкинс Г., Ватте Д. Спектральный анализ и его приложения, том 1 и 2, М., Изд. Мир, 1971.

22. Дюран Б. Оден П. Кластерный анализ. М.: Статистика, 1977. 127 с.

23. Ермаков С.А., Гольдблат В.Ю., Зуйкова Э.М. и др. Натурные оптические и радиолокационные измерения ветровых гравитационно-капиллярных волн и их гашения поверхностно-активными пленками. Предварительные результаты. Нижний Новгород, 1996. -23 е.: ил.

24. Ефимов В.В. Динамика волновых процессов в пограничном слое атмосферы и океана. Киев, 1981, 253 с.

25. Заславский Г.М., Шарков Е.А. Фрактальные свойства зон обрушающихся волн на поверхности моря. Докл. АН СССР. Механика жидкостей и газов, 1987, т. 294, с. 1362-1366.

26. Захаров В. Е., Заславский М.М. Форма спектра энергонесущих компонент водной поверхности в слабо турбулентной теории ветровых волн. Изв. АН СССР. Физика атмосферы и океана. 1983, Т. 19. №3. С, 282-291.

27. Зуйкова Э.М., Титов В.И. Исследование модуляции ряби длинными поверхностными волнами с помощью оптических методов/. -Нижний Новгород, 1997. -11 с.

28. Ильин Ю.А., И.И. Стрижкин, Малинников В.А. Комплексный метод исследования высокочастотных ветровых волн. Тез. доклада. Региональная конференция «Ветровое волнение».Севастополь,26-27 февраля, 1990, с. 1617.

29. Короленко П.В. и др. Новационные методы анализа стохастических процессов и структур в оптике. Москва, изд. НИИЯФ МГУ, 2004, 82 с.

30. Киселев А. Непрерывное вейвлет-преобразование в анализе безнес-информации. Интернет-сайт BaseGroup Labs.htm, 2003.

31. Кронберг П. Дистанционное изучение Земли (Основы и методы дистанционных исследований в геологии). Перевод с немецкого. М.: Мир, 1988.

32. Кудрявцев В.Н. и др. Экспериментальные исследования механизмов модуляции радиолокационного сигнала на масштабах морских поверхностных волн. Исследование Земли из космоса, 2001, №4, с. 13-30.

33. Лупян Е.А., Мурынин А.Б. Возможности фрактального анализа оптических изображений морской поверхности. М.: Ротапринт ИКИ АН СССР, 1989.

34. Малинников В.А и др. Особенности пространственно-временной структуры морского волнения по данным фотосъемки и контактных измерений. Изв. вузов. Геодезия и аэрофотосъемка. М, 1990, № 1, с. 99-105.

35. Малинников В.А. и др. Калибровка дистанционных измерений пространственных спектров волнения по оптическим изображениям. Исследование Земли из космоса, М., № 2, 1992.

36. Малинников В.А. и др. Метод регистрации высокочастотного спектра ветрового волнения в реальном масштабе времени. Тезисы докл. IV Всесоюзной конф. «Мировой океан», ч. 2 Владивосток, 1983, с. 72—73.

37. Малинников В.А. и др. Пространственная структура высокочастотного ветрового морского волнения при различных метеоусловиях. Изв. АН СССР, ФАО, 1985, т.21, N4, с.440-442.

38. Малинников В.А. и др. Синхронные оптические и контактные исследования пространственно-спектральных характеристик морского волнения. Исследования Земли из космоса, 1986, № 2, с. 57—67.

39. Малинников В.А., Корчинский Е.В. Исследование фрактальных свойств бассейнов речных систем. Труды Международной научно-технической конференции, посвященной 225-летию МИИГАиК, 2004.

40. Малинников В.А., Савиных В.П., Малинникова Е.В. Мультифрактальный анализ и параметризация пространственных спектров космических изображений облаков. Тез. Докладов 4-го Международного аэрокосмического конгресса, Москва, изд. СИПРИА, 2003.

41. Малинников В.А., Стрижкин И.И. Экспериментальные исследования особенностей структуры спектров возвышений и уклонов ветрового морского волнения. Тез. доклада на 3 съезде советских океанологов.Л., 1987.

42. Малинников В.А., Стрижкин И.И., Тихомиров О.А. Анализ методов определения уклонов волн по фотоснимкам. Геодезия и аэрофотосъемка. Известия вузов высших учебных заведений. Москва, 1993, №3

43. Малинников В.А., Сухорученко А.Н. Оценка фрактальной размерности оптических изображений морской поверхности. Труды Международной научно-технической конференции, посвященной 225-летию МИИГАиК, 2004, с. 252-259.

44. Малинников В.А., Учаев Д.В. Мультифрактальная параметризация геопространственных структур. Труды Международной научно-технической конференции, посвященной 225-летию МИИГАиК, 2004.

45. Малинецкий Г.Г., Потапов А.Б. Современные проблемы нелинейной динамики. Изд.УРСС, Москва, 2002, 358 с.

46. Мандельброт Б.Б. Фрактальная геометрия природы. М., Ижевск: РХД, 2002.

47. Методы компьютерной обработки изображений / Под ред. В.А. Сойфера.-2-е изд., испр.-М.: ФИЗМАТЛИТ, 2003. 784с.

48. Методы обработки и преобразования дистанционных данных для целей региональных ГИС. /Научно-технический сборник. ЦНИИГАиК, М., 1996; 32с.

49. Монин А.С., Красицкий В.П. Явления на поверхности океана. Л.: Гидрометеоиздат,1985, 376 с.

50. Мун Ф. Хаотические колебания: Пер. с англ. — М.: Мир, 1990.

51. Пайтген Х.О., Рихтер X. Красота фракталов. — М.: Мир, 1993.5 5.Потапов А. Фракталы в дистанционном зондировании // Успехи современной радиоэлектроники.- 2000. № 6.

52. Применение аэрометодов для исследования моря. Под ред. Здановича В. Г. М.- Л.: Изд-во АН СССР, 1963. 544 с.

53. Смоленцев Н.К. Основы теории вейвлетов.Вейвлеты в MATLAB. Москва, издДМК, 2005,303 с.

54. Сухорученко А.Н. Исследование фрактальных свойств поля блика на морской поверхности. Сб. научных трудов VIII межвузовского научно-практического семинара «Экологическая безопасность и рациональное природопользование», М., 2004, с. 139-142.

55. Сухорученко А.Н. Методика локального структурно-спектрального анализа оптических изображений морской поверхности. Изв. Вузов. Геодезия и аэрофотосъемка. Специальный выпуск, 2006, с. 131-140.

56. Сухорученко А.Н. Фрактальные свойства зон обрушений морских гравитационных волн. Сб. научных трудов VIII межвузовского научно-практического семинара «Экологическая безопасность и рациональное природопользование», М., 2004, с. 143-147.

57. Тимофеев П.В., Шарков Е.А. Натурные исследования дисперсной фазы зоны обрушения морской волны. М., 1992. 33 с.

58. Титов В. И. Определение спектра волнения путем спектрального анализа аэрофотоснимков.— Изв. АН СССР. Физика атмосферы и океана, 1982, т. 18, №2, с. 215-216.

59. Федер Е. Фракталы: Пер. с англ. М.: Мир, 1991, 254 с

60. Филлипс О. М. Динамика верхнего слоя океана. Л.: Гидрометеоиздат, 1980. 108 с.

61. Фракталы в физике. Ред. Л. Пьетронеро, Л. М. Тозатти. М.: Мир, 1988, 205с.бб.Черный И.В.,Шарков Е.А. Радиодистанционные исследования процесса обрушения морской волны. Исследования Земли из космоса, №2, 1988.

62. Шарков Е.А. Оптические исследования временной динамики пенных структур на взволнованной морской поверхности. Исследования Земли из космоса, №1, 1994, с. 68-78.

63. Шарков Е.А., Покровская И.В. Оптико-дистанционные исследования полей обрушений гравитационных волн развивающегося морского волнения. Исследования Земли из космоса, №3, 1998, с. 11-22.

64. Шарков Е.А., Покровская И.В. Пенная активность на морской поверхности как марковский случайный процесс. Докл. АНСССР, 1987, т.293, №5, с. 1108-1111.

65. Шулейкин В.В. Физика моря. М. Наука, 1968.

66. Шустер Г. Детерминированный хаос. Пер. с англ. — М.: Мир, 1988

67. Юргенс X., Пайтген Х.-О., Заупе Д. Язык фракталов. В мире науки (Scentific American), 1990, № 10, 36-44.

68. Alparone L., Barni M., Betti M. "Fuzzy clustering of textured SAR images based on a fractal dimension feature", invited paper, IEEE IGARSS'97, Singapore, August 4-8, 1997, pp. 1184-1186

69. Andreoli G., Berizzi F., Dalle Mese E." A two-dimensional fractal model of the sea surface and sea spectrum evaluation", in Proc. IEE International Radar Conference, Edinburgh, Scotland, October 1997, pp. 189-193

70. Benelii G., Garzelli A. "Oil-spills detection in SAR images by fractal dimension estimation", Proc. IEEE IGARSS'99, Hamburg, Germany, 28 June 02 July, 1999

71. Bergamasco, L., Serio, M., 1996, "Chaotic indicators and Gaussian random processes: Some surprising results", Fractals, vol. 4, no. 1, pp. 73-90.

72. Berizzi F., Dalle Mese E. "Fractal analysis of the signal scattered from the sea surface", IEEE Transaction on Antennas and Propagation, Vol.2, No.24, pp.324338, Feebruary 1999

73. Berizzi F., Dalle Mese E., Pinelli G. "One-dimensional fractal model of sea surface", IEE Proceedings Radar, Sonar and Navigation, Vol. 146, No.l, pp.55-64,1999

74. Сох С. H., Munk W. Statistics of the sea surface derived from sea glitter J.Mar. Res., vol. 13, № 2, pp. 198—227.

75. Crawford D.R. and etc. Stability of weakey nonlinear deep-water in two and three dimensions. J.Fluid Mech.,1981, vol.105.

76. Dunckel D.E.JEwing J.A. Direction Wave spectra Observed During Jonswap из жури.: Journal of Physical Oceanography. -1980. Vol.10. - P.1264-1280.

77. ERDAS Field Guide. Fourth Edition, Revised and Expanded. /Atlanta, Georgia.: ERDAS, Inc., 656p.

78. ERDAS Imagine software //ERDAS Inc. 1995. - Vol. 1-7. - 1458 p.

79. Kudryavtsev V., Makin V., Chapron B. Coupled sea surface atmosphere model. 2. Spectrum of short wind waves. J. Geophys.Res 1999, vol.104, NC4, p.7625-7639.

80. Leykin I.A., Rosenberg A.D. Sea-Tower measurements of wind-wave spectra in the Caspain Sea // J. Phys. Oceanogr. 1984. V. 14. P. 168-176.

81. Longuet-Higgins, M.S., 1985, "Bifurcation in gravity waves", J. Fluid Mech., vol. 151, pp. 457-475.

82. Longuet-Higgins, M.S., 1994, "A fractal approach to breaking waves", J. Phys. Oceanogr, vol. 24, no. 8, pp. 1834-1838.

83. Marazzi, P. Gamba, E. Costamagna, N. Piredda: "Texture classification with fractal dimension and wavelet transform", Proc. of the 2nd International Conf. on Chaos, Fractals and Models, Pavia, Italy, Oct. 1996.

84. Morrison, A.I., Srokosz, M.A., 1993, "Estimating the fractal dimension of the sea surface: A first attempt", Ann. Geophys., vol. 11, no. 7, pp. 648-658.

85. Plass G.N.,Kattavar G.W. Humphreys T.J. Influence of the oceanic scattering phase function on the radiance // J. Geophys. Research. 1985. V. 20. №102. P. 3347-3351.

86. Stiassnie, M., Agnon, Y., Shemer, L., 1991, "Fractal dimensions of random water surfaces", Physica D: Nonlinear Phenomena, vol. 47, no. 3, pp. 341-352.

87. Stive M.J.F. Cross-shore flow in waves breaking on a beach: Diss./. -Delft, 1988. -Pag.var.: ill.Delft hydraulics publications;N395, ISSN 0166-7521).

88. Troickay Yu. I. Wind excitation of surface waves in the coupled air-water turbulent shear flow. A simple model of visco-elastic turbulence. -Nizny Novgorod, 1997. 43 p.

89. Zlatopolsky A. Description of texture orientation in remote sensing data using computer program LESSA.

90. Computers & Geosciences. 1997. - V. 23. - №. 1. p. 45-62.

91. Zlatopolsky A. Program LESSA (Lineament Extraction and Stripe Statistical Analysis). Automated linear image features analysis -experimental results. /Computers & Geosciences, v. 18, N9, 1992. p. 1121-1126.