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

Министерство образования и науки Российской Федерации ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ

ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ "РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ"

(РГГМУ)

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

005538330

Воинов Никита Евгеньевич

МЕТОДЫ ЦИФРОВОМ ОБРАБОТКИ МНОГОЗОНАЛЬНЫХ

СПУТНИКОВЫХ СНИМКОВ В ЗАДАЧАХ АНАЛИЗА МЕЗОМАСШТАБНЫХ КОНВЕКТИВНЫХ АТМОСФЕРНЫХ

ДВИЖЕНИЙ

25.00.30 - Метеорология, климатология, агрометеорология

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

14 НОЯ ¿013

Санкт-Петербург 2013

005538330

Министерство образования и науки Российской Федерации ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ

ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ "РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ"

(РГГМУ)

МЕТОДЫ ЦИФРОВОЙ ОБРАБОТКИ МНОГОЗОНАЛЬНЫХ

СПУТНИКОВЫХ СНИМКОВ В ЗАДАЧАХ АНАЛИЗА МЕЗОМАСШТАБНЫХ КОНВЕКТИВНЫХ АТМОСФЕРНЫХ

ДВИЖЕНИЙ

25.00.30 - Метеорология, климатология, агрометеорология

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

Войнов Никита Евгеньевич

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

РГГМУ

Санкт-Петербург

2013

Работа выполнена в ФГБОУ ВПО «Российский Государственный гидрометеорологический университет».

Научный руководитель Научный консультант Официальные оппоненты

Ведущая организация

доктор физ.-мат. наук, профессор А.Д. Кузнецов кандидат физ.-мат. наук, доцент Т.Е. Симакина

доктор физ.-мат. наук, профессор Л.С. Ивлев, кандидат физ.-мат. наук, доцент В.Д. Еникеева

Международный центр по окружающей среде и дистанционному зондированию имени Нансена

-30

Защита состоится «5» декабря 2013 г. в 15 " на заседании диссертационного совета Д212.197.01 при Российском Государственном гидрометеорологическом университете.

Адрес: 195196, г. Санкт-Петербург, Малоохтинский пр., 98.

С диссертацией можно ознакомиться в библиотеке Российского Государственного гидрометеорологического университета.

Автореферат разослан « 1 » ноября 2013 г.

Ученый секретарь

диссертационного совета Д212.197.01, кандидат физ.-мат. наук

Кашлева Л.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы

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

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

Предмет исследования

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

Цель исследования

Целью настоящей работы являются:

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

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

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

3

О

Для достижения поставленной цели автором:

1. Рассмотрены спутниковые средства мониторинга мезомасштабных конвективных систем.

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

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

4. Разработана методика расчета текстурной сложности области, занятой разномасштабными облачными системами, на основе вычисления ее фрактальной размерности.

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

Методы исследования

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

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

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

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

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

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

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

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

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

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

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

Апробация результатов

Основные результаты работы докладывались на Итоговой сессии Ученого совета РГГМУ (Санкт-Петербург, январь 2013 г.), на научных семинарах кафедры экспериментальной физики атмосферы РГГМУ (2013 г.).

Основные результаты диссертации опубликованы в двух печатных работах, в том числе две - в издании, рекомендованном ВАК.

Структура и объем

Материал диссертации состоит из введения, трех глав, заключения, списка литературы. Текст содержит 154 страницы, включая 76 рисунков, 20 таблиц. Список использованных источников содержит 101 наименование.

Содержание диссертационной работы

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

Первая глава содержит информацию как об объекте исследования — атмосферных мезоконвективных процессах, - так и о спутниковых средствах мониторинга этих процессов.

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

вективных процессов требует соответствующего диапазона пространственного и временного разрешения спутниковых данных.

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

Данный раздел содержит обзор существующих классификаций мезомас-штабных конвективных систем, выявлены ограничения таких классификаций. Составлена классификация систем мелкой конвекции по аналогии с наиболее популярной и удобной как с генетической, так и с морфологической точки зрения классификацией систем глубокой конвекции Мэддокса. Однако отмечается, что основное внимание в существующих классификациях обращается на горизонтальные размеры и морфологическое описание рассматриваемых процессов. Для целей же сверхкраткосрочного прогноза и пошсав!^ необходимым условием является также информация о времени полного цикла "жизни" объекта и о стадии цикла, на которой данный объект в конкретный момент находится. Разработка методов идентификации стадии развития мезоконвективной системы в режиме реального времени является важной задачей.

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

Системы мелкой конвекции не отражаются на синоптических картах и редко продуцируют опасные явления, однако анализ пространственной и временной изменчивости мелкой конвекции необходим как дополнительная, а над районами, малоосвещеннными данными наблюдений, как основная информация для оценки синоптической обстановки и возможных условий погоды. Тип конвективных ячеек может служить довольно надежным индикатором воздушных масс. Распространение открытых облачных ячеек позволяет следить за проникновением масс холодного воздуха, закрытых - теплого. Учет влияния конвективных ячеек позволяет локализовать появление экстремальных концентраций примеси в приземном слое. Образование крупных открытых конвективных ячеек диаметром более 30 км часто указывает на штормовой характер погоды в океане: проходят кратковременные, но интенсивные осадки, сопровождающиеся шквалистым усилением ветра. Причем форма открытых облачных ячеек связана с высотой волн в районах океана, над которыми прослеживаются облачные образования: симметричные ячейки соответствуют волнам, высота которых не превышает 4 м, ячейки неправильной (серповидной) формы указывают на значительное волнение в 5 - 6 м. С мелкими закрытыми конвективны-

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

Раздел 1.3.2 содержит результаты исследования мезомасштабных конвективных комплексов над территорией Европы, проведенного в рамках настоящей работы по данным 2008 - 2011 гг. геостационарного спутника Ме1ео5а1-8 - рис.1. МКК определялись по признакам Мэддокса. В данном разделе отмечено, что с точки зрения сверхкраткосрочного прогноза и по\уса51:т§ идентификация МКК по признакам Мэддокса лишена смысла, поскольку необходимое условие сохранения геометрических размеров в течение 6 часов, входящее в перечень идентифицируемых признаков, означает, что от момента появления МКК до его окончательной идентификации проходит столько времени, что большая часть опасных явлений к этому моменту уже произошла. Другим признаком является наличие округлой области перистых облаков определенного размера, однако такая область зачастую маскирует различные по форме скопления кучево-дождевой облачности, в том числе и мезомасштабные линии шквалов.

Для всех идентифицированных МКК измерялось время их существования в поле облачности, а также направление перемещения путем построения вектора, проведенного между центрами тяжести облачного массива МКК в стадиях зарождения и диссипации, нанесенного на физическую карту Европы масштаба 1:250000.

В результате обработки спутниковой информации было установлено, что в среднем за летний период на территории Европы образуется около 20 мезомасштабных конвективных комплексов. В отличие от территории Северной Америки, где наибольшая повторяемость и интенсивность возникновения МКК наблюдается в весенние месяцы, над Европой в переходные сезоны были зафиксированы лишь единичные случаи возникновения МКК: в апреле, мае и сентябре - рис.1 (а). Связаны они с очень активными атмосферными фронтами. Регулярное формирование МКК на территории Европы начинается с десятых чисел июня и отмечается практически до конца августа, пик приходится на июль. В холодное полугодие над территорией Европы мезомасштабные конвективные комплексы не формируются.

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

По результатам исследования наиболее часто МКК образуются в центральной и восточной Европе, причем большинство из них (около 80% случаев) формируется под воздействием орографии, т.е. в районе Кавказского хребта, Карпат и Альпийских гор — рис. 1(6). Над поверхностью моря случаев возникновения МКК не зафиксировано. Также можно отметить отсутствие случаев МКК на широтах выше 60°с.ш.

юв ююв

Рис. 1. Результаты исследования МКК над территорией Европы: а) гистограмма случаев возникновения МКК; 6) районы наибольшей повторяемости МКК; в) диаграмма направлений

смещения МКК

Большинство МКК смещается согласно ведущему потоку, господствующему над Европой в теплый период года, а именно в северо-восточном и северном направлении - рис. 1(в).

Характерным отличием является длительность жизненного цикла мезо-масштабных конвективных комплексов, которая над США составляет 15—16 часов, в то время как над Европой облачная система, связанная с МКК, наблюдается почти вдвое меньше по времени - около 8-9 часов.

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

Высота верхней границы облачности, ее фазовый состав, направление и интенсивность вертикальных токов и другие факторы влияют на интенсивность излучения и отражение радиации, что в свою очередь формирует яркость пикселей многозонального спутникового снимка. Сравнение яркостей объектов и их разностей в различных каналах дает возможность получить несоизмеримо большую информацию о процессах конвективного облакообразования, о стадии жизни отдельных ячеек конвективной системы, чем анализ одиночных снимков. Так, синтез трех изображений в одно, осуществляемый в рамках той или иной колориметрической модели, позволяет получить мультипликативный эффект. Для 12-тиканального сканера SEVIRI объединение по три исходных снимка в модели RGB дает 12x11x10 = 1320 комбинаций. Количество комбинаций возрастает, если использовать разности каналов, которые, как уже упоминалось, несут большую информацию.

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

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

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

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

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

В подразделе 2.2.1 отмечается, что цветное изображение передает существенно больший объем информации об объекте съемки, чем черно-белое, что способствует более точному решению задач дешифрирования (опознавания, классификации). Методы анализа цветных изображений зависят от цветового пространства. В подразделе 2.2.2 рассмотрены основные принципы следующих цветовых моделей: RGB, ХУ2, Lab, HSB (HSL, HSV). Приводятся достоинства и недостатки рассмотренных моделей, даны формулы перехода из пространства RGB в Lab через промежуточное пространство XYZ и из RGB в HS В.

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

Результат идентификации объектов на композитном изображении, составленном по схеме RGB из каналов VIS 0.8, IR 3.9 и IR 10.8 мкм, представлен на рис. 2.1. Три ряда цветных прямоугольников отражают цвет верхней грани-

цы облаков, находящихся на трех уровнях по высоте. Трехстолбчатые диаграммы показывают вклад в результирующий цвет яркостей исходных снимков. Яркости в канале 0,8 мкм, чувствительном к плотности облачности и к количеству капель и ледяных кристаллов, модулирующих отражаемость солнечного излучения, управляют красным столбиком диаграммы. Присутствие среднеин-фракрасного канала 3,9 мкм, управляющего зеленым столбиком, делает такую комбинацию чувствительной к размеру и фазе облачных частиц. Канал 10,8 мкм, модулируемый температурой верхней границы, управляет синим столбиком. Данная комбинация VIS 0,8; IR 3,9; IR 0,8 доступна только в дневное время из-за наличия канала VIS 0,8 мкм.

а)

6)

J.

д)

ш

и)

I

он

I. К

ж)

II

Л)

Растительность

Рис. 2.1. Цвета облачности и объектов подстилающей поверхности в схеме RGB 0.8, 3.9, 10.8: Верхний ярус: а) тонкая перистообразная мелкокристаллическая; б) тонкая перисто-образная крупнокристаллическая; в) конвективная облачность с мощным вертикальным подъемом и осадками; г) конвективная облачность.

Средний ярус: д) тонкий слой с переохлажденными мелкими каплями; е) тонкий слой с переохлажденными крупными каплями; ж) мощный слой с переохлажденными мелкими каплями; з) мощный слой с переохлажденными крупными каплями.

Нижний ярус: и) тонкий слой с мелкими каплями; к) тонкий слой с крупными каплями; л) мощный слой с мелкими каплями без осадков; м) мощный слой с крупными каплями с осадками

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

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

Далее рассмотрение процесса конвекции проводится в цветовой системе Lab. Компоненты этой системы - Lightness, а и Ь, - рассчитывались по формулам подраздела 2.2.2.

Вследствие разделимости информации о цвете и яркости в системе Lab появляется возможность проводить анализ развития конвективной облачности в прямоугольной плоскости хроматических компонент а — Ь. Начало координат находится в центре в точке белого цвета. Пример анализа значений интенсивности в трех каналах К/5 0.8, IR 3.9 и IR 10.8 мкм (см. рис.2.1) приведен на рис. 2.2. Цвета облачности в начальной фазе развития (мелкие частицы, малая оптическая толщина слоя) и конечной (крупные частицы, большая оптическая толщина) на каждом ярусе представлены точками.

Римские цифры обозначают ярус верхней границы облачности. Цифры расположены у точки наименьших толщины облачного слоя и размера частиц на соответствующей высоте. Рост размера частиц на любой высоте отражается увеличением координаты а и уменьшением координаты b (синие линии на рисунке). Возрастание толщины облачного слоя без укрупнения частиц сопровождается ростом координаты b (серые линии). Красный пунктир указывает на одновременный рост размера частиц и толщины облачного слоя. Заметно возрастание амплитуды изменения (длина красного пунктирного отрезка) с ростом высоты верхней границы.

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

Рис. 2.2. Этапы развития облачности в координатах a- b

таблицы для схемы, составленной из каналов VIS 0.6; NIR 1.6; IR 10.8, приведен далее.

Таблица 2.1

Цвета верхней границы СЬ на разных стадиях в комбинации VIS 0.6; NIR 1.6;

IR 10.8

Стадии R G В Н S Вг Цвет

I 255 255 200 60 22 100

II 255 200 0 47 100 100

III 255 100 0 24 100 100

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

В таблицах также представлены значения компонент моделей RGB и HSB каждого цвета.

Проведен анализ развития конвективной облачности в полярной плоскости компонент цветовой тон (Я) - насыщенность (S), где Я является аналогом длины волны излучения, S - аналогом интенсивности волны. Значения цветового тона откладывают на периферии цветового круга и характеризуют величиной угла в диапазоне от 0 до 360°. Насыщенность цвета - радиус цветового круга. Отображение цветов верхней границы СЬ на трех стадиях развития в плоскости Н- S, представлено на рис. 2.3.

120° 60°

Рис. 2.3. Цветовая эволюция в плоскости тон-насыщенность в комбинации каналов

Ж 0.6; Л®? 1.6; Ш 10.8

Анализ всех рассмотренных в данном параграфе диаграмм позволил сделать следующие выводы. Направление векторов, соединяющих точки I, II и III, соответствующих номерам стадий, стремится в сторону увеличения длины волны, т.е. к красному тону со значением 0°. Значение насыщенности при этом всегда возрастает, что особенно характерно для вектора, соединяющего I и II стадии. По мере развития конвекции цвет стремится к максимально чистой тоновой волне, поэтому локализация цвета на диаграмме может позволить оценить степень эволюции СЬ, а, следовательно, косвенно и угрозу опасных явлений ими вызываемых.

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

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

Раздел 3.1 содержит как теорию дискретного (быстрого) вейвлет-преобразования (БВП), так и практические примеры разложения спутникового снимка разномасштабных облачных систем.

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

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

= 2(рср^ х - к,у - /), =2'Г щ/^ х-к,2-1 у-1), = 23 у^2-7 * - У ~ 1\

= 2У у(х)\|/(л: ^ х-к,2-1 у-¡],

где к - целые числа.

Рекурсивное использование процедуры свертки происходит с уменьшением количества отсчетов в 2 раза при переходе от одного уровня («масштаба») к другому.

В подразделе 3.1.2 отмечается, что присутствие на исходном снимке разномасштабных облачных объектов, как в примере на рис.3.1, - циклона синоптического масштаба, линии шквалов масштаба мезо-а, скопления кучево-дождевой облачности масштаба мезо-0 и отдельных облаков масштаба мезо-у, составляющих поля конвективных ячеек, - не позволяет применять ко всему изображению одну простую статистическую модель. Результат четырех уровней вейвлет-декомпозиции представлен на рис. 3.1(6). Исходное изображение

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

Рис.3.1. Вейвлет-преобразование спутникового снимка: а) исходный снимок видимого диапазона, МОШБ/Тегга, 29.04.2009, 600x860 пиксел, пространственное разрешение 250 м.; б) четыре уровня вейвлет-декомпозиции

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

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

Для описания текстуры изображений облачности использовался статистический подход, в основе которого лежит формирование матриц пространственной смежности яркости пикселей. Матрица смежности содержит относительные частоты р{ ] наличия на изображении соседних элементов, расположенных на расстоянии с/ друг от друга, с яркостями г,у. Средняя по направлениям матрица р^ = 1/4(/?10 + р^ +¿>1,90 +Р1,135) использовалась для расчета некоррелированных текстурных признаков.

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

Однородность является мерой гомогенности изображения:

N N

Энтропия выражает неравномерность распределения яркостных свойств элементов изображения:

Обратный момент отражает степень разброса элементов матрицы градиентов вокруг главной диагонали:

Анализ гистограмм текстурных признаков, вычисленных для каждого типа систем мелкой конвекции и представленных на рис. 3.2, позволил сделать следующие выводы. Увереннее всего выделены закрытые ячейки по трем текстурным признакам: контраст, энтропия и обратный момент, причем меньшую дисперсию и наибольшее расстояние от других текстур имеет обратный момент. Для открытых ячеек и облачных гряд наиболее компактным диапазоном значений обладает текстурный признак контраст. Наименьшую разделимость представленных мезоструктур показал признак однородность. Совместное использование признаков обратный момент, контраст и энтропия повысит возможность идентификации типов конвекции и позволит локализовать границу между ними. Изображения данных текстурных признаков приведены на рис. 3.3. Представив текстурные признаки в кодировке RGB: обратный момент красным цветом, контраст - зеленым и энтропию - синим, получим псевдоцветное текстурное изображение — рис. З.З(г). Высокие значения обратного момента и контраста присущи открытым ячейкам, поэтому такие объекты на текстурном ^GiJ-изображении будут иметь оттенки смеси двух цветов зеленого и красного, т. е. оттенки желтого цвета. Закрытые ячейки обладают высокими значениями энтропии (синий цвет), контраста (зеленый цвет) и низкими значениями обратного момента (красный цвет). Т.е. для закрытых ячеек характерным цветом будут оттенки синего и голубого цветов. Текстурное цветное изображение позволяет сегментировать различные виды мелкой конвекции и локализовать границу воздушных масс с разными свойствами.

Гистограмма значений контраста

ilk_ні ні і ill пі

Значения контраста

Однородности

Значения однородности

Гистограм м а значений энтропии

Количество значений ill

£ -if ^ ^ С^СГ^

Значения энтропии

Гистограмма значений обратного момента

я я я й я й з "9 ? g в s i

Значения обратного момента

закрытые ячейки ■ открытые ячейки облачные гряды О

Рис. 3.2. Гистограммы текстурных признаков: контраст, энтропия, однородность и обратный момент

а) б) в) г)

Рис.3.3. Текстурные признаки и псевдоцветное текстурное изображение: а), б) и в) обратный момент, контраст и энтропия соответственно; г) RGB — синтез текстурных признаков: обратный момент - (R), контраст - (G) и энтропия - (5)

В разделе 3.3 рассматривается разработанный алгоритм определения высоты конвективных ячеек открытого и закрытого типа с помощью текстурных признаков. Алгоритм позволяет в случае полного отсутствия инструментальных измерений в районе с однородной мелкоконвективной облачностью грубо оценить мощность слоя конвекции по спутниковому снимку. Алгоритм опирается на известную статистическую связь между горизонтальными и вертикальными размерами элементов систем мелкой конвекции и использует эмпирические соотношения для каждого типа мезоструктуры: h = D/30 для открытых ячеек, h = D/21 для закрытых ячеек, где h — высота слоя конвекции, D - диаметр ячеек.

Построены графические зависимости текстурных признаков и мощности конвективного слоя на основе анализа шестидесяти фрагментов спутниковых

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

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

1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0

у = -6,Ох + 2,6

г =-0,74

0.25 0.3

Значение обратного момента

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

0.55 0.6 0.65

Значение обратного момента

Рис. 3.4. Связь мощности конвективного слоя закрытых (слева) и открытых ячеек (справа) и текстурного признака обратный момент

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

Поверхность облачного поля представляет из себя фрактал, размерность которого находится в диапазоне от 2 до 3. Способы расчета фрактальной размерности приводятся в подразделе 3.4.1.

Алгоритм вычисления фрактальной размерности текстуры поверхности, модернизированный и используемый в данной работе, является развитием алгоритма оценки фрактальной размерности линии. Предварительно, с помощью программы АБСАЫ, разработанной ранее на кафедре ЭФА, РГГМУ, по ИК-изображению производится формирование массивов бинарных изображений, содержащих области облачности, охваченные различными изотермами выше температурного порога. Для каждого из таких изображений производится вычисление фрактальной размерности, среднее из полученных значений используется в качестве оценки фрактальной размерности исходного изображения.

Процедура вычисления фрактальной размерности контура бинарного изображения строится следующим образом. Для серии размеров элемента разложения (шагов) Я,- производится оценка длин контуров единичной области Ь путем подсчета числа переходов уровней яркости бинарного изображения через ноль. Увеличение шага эквивалентно интерпретации анализируемого изображения с уменьшенным разрешением. Фрактальная размерность рассчитывается из наклона |3, аппроксимированной методом наименьших квадратов линии на графике в двойном логарифмическом масштабе как £> = 3 - (3.

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

= cov(logL,logs) = E[[\ogL - M\ogi \\ogS - MlojS )j ri°gi'l0gS /-ClogX x D\ogS fDiogL x DlogS

где D - дисперсия, М - математическое ожидание соответствующих случайных величин log L и log S.

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

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

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

Поле фрактальных размерностей, построенное по спутниковым снимкам мезомасштабной облачности, позволяет оценить текстурную сложность района, занятого разными типами систем мелкой конвекции. Пример построения такого поля для территории рис. 3.5 представлен в градациях серого на рис. 3.16: чем ярче - тем выше фрактальная размерность. Диапазон изменения фрактальной размерности ячеек открытого типа составил 0,25, в то время как среднее значение фрактальной размерности приблизительно равно 2,8. Диапазон изменения фрактальной размерности закрытых ячеек составил 0,45, со средним значением 2,5. Строчный анализ поля фрактальных размерностей может использоваться для оценки степени неоднородности поверхности облачных ячеек и позволяет локализовать границу открытые-закрытые конвективные ячейки, которая соответствует переходу от холодных воздушных масс к теплым. В качестве иллюстрации приведен рис. 3.7, где изображено изменение фрактальной размерности по направлению, обозначенному на исходном изображении рис. 3.5 белым цветом. Четко виден скачок значений фрактальной размерности на границе открытые-закрытые ячейки, в данном случае этому состоянию соответствует 415 точка.

Рис.3.5. Поле конвективных ячеек. Рис.3.6. Поле фрактальных

MSG, 23.08.2001, 11:00 UTC размерностей

Рис.3.7. Профиль фрактальной размерности по направлению, указанному белой линией на рис. 3.5 (слева-направо)

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

Таблица 3.1

Фрактальная размерность исходного изображения рис.3.1 и четырех уровней вейвлет-_преобразования___

Разрешение, м 250 500 1000 2000 4000

Фрактальная размерность 2,788 2,6634 2,4894 2,3012 2,1834

СКО 0,1802 0,2204 0,2294 0,205 0,1596

Диапазон фраталь-ной размерности 2-2,9831 2-2,9565 2,0003 -2,9185 2-2,8428 2,0001 -0,6762

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

Произведено сопоставление вычисленных значений фрактальной размерности со спектром цветных шумов. Размерность близкая к 2 для плоскости определяется как «черный шум» и связывается обычно с процессами в очень вязкой среде. Такая размерность характерна для выровненных поверхностей, для которых величины оптической плотности изменяются незначительно. Тип размерности 2,5 обычно называют «коричневым шумом», и он отражает формы поверхности, порождаемые процессами, связанными с теплопереносом, или диффузией, в основу которых может быть положена модель случайного блуждания. На снимке такая поверхность отражается через резкие изменения оптической плотности с определенной степенью регулярности, для которой характерны высокочастотные (низкопериодические) составляющие. Тип размерности 2,9 называется «розовым шумом» и связывается с турбулентными процессами в среде малой вязкости.

На изображении поля фрактальных размерностей (см. рис. 3.6 и 3.8) светлому тону соответствуют территории, фрактальная размерность которых близка к «розовому шуму», серому тону - к «коричневому шуму», а наиболее темному тону - территории с фрактальной размерностью близкой к «черному шуму».

На основе спектра «цветности шума» полученные оценки фрактальной размерности были связаны с атмосферными процессами, вследствие чего сделаны следующие выводы. При самом высоком разрешении (250 м) к «розовому шуму» близка фрактальная размерность открытых ячеек. Так как текстура территории есть функция мощности действия факторов, определивших ее генезис, то можно сказать, что на объекты с текстурой, соответствующей «розовому шуму», действует множество возмущающих факторов. Действительно, характеристики образующихся в пограничном слое атмосферы над океаном облачных ячеек, состоящих из облаков масштаба мезо-у, определяются совместным воздействием нескольких механизмов: нагревом на нижней границе, радиационными процессами на верхней границе, высвобождением латентного тепла в облаках, крупномасштабными вертикальными движениями, температурой и влажностью воздуха в слое инверсии, вовлечением и др. процессами.

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

При разрешении 500 м и особенно 1 км, что соответствует снимкам с типичных метеорологических спутников (ЫОАА), средняя фрактальная размерность соответствует "коричневому шуму". На изображениях первого и второго уровня вейвлет-декомпозиции просматриваются системы масштаба мезо-а и -р с соответствующим временем жизни порядка суток.

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

неизменные во времени, что характерно для синоптических объектов (их фрактальная размерность на всех иерархических уровнях близка к 2). Доля территорий с "черным шумом" по визуальной оценке существенно больше доли территорий с "розовым шумом". Разрешение последнего уровня (4 км) типично для геостационарных спутников.

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

"розовый шум"

0 12 3 4

<г- (3 -» <-а->

Разрешение, км

Рис.3.9. Фрактальная размерность на разных масштабах

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

2. Выполнено исследование мезомасштабных конвективных комплексов над территорией Европы по данным геостационарного спутника Meteosat-8 (2008 - 2011 гг.). На основе сравнительного анализа полученных данных с данными США выявлены особенности их статистических характеристик.

3. Выполнен анализ многозональных снимков в различных колориметрических моделях. Для анализа стадий развития конвекции предложено использование прямоугольной плоскости хроматических компонент колориметрической модели Lab и полярной плоскости Тон - Насыщенность колориметрической модели HSB.

4. Предложено и обосновано применение вейвлет-преобразования для анализа разномасштабных облачных процессов при разном разрешении.

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

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

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

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

8. Модернизирована методика определения фрактальной размерности текстур мезомасштабной облачности по ИК снимкам.

Основные положения диссертации представлены в следующих публикациях:

1. Воинов Н.Е, Кузнецов А.Д., Симакина Т.Е., Сероухова О.С. Измерение фрактальной размерности мезомасштабной конвективной облачности на спутниковых снимках.// Ученые записки РГГМУ. 2013. № 30. С.93-100.

2. Войнов Н.Е., Кузнецов А.Д., Симакина Т.Е., Сероухова О.С. Кратно-масштабный анализ спутниковых снимков облачности. //Вестник Тверского государственного университета. Серия: "Прикладная математика", выпуск 1 (28). 2013. № 15. С. 81-90.

3. Отчет о научно-исследовательской работе «Исследование дистанционного обнаружения опасных явлений погоды с помощью автоматизированного импульсного когерентного доплеровского метеорологического радиолокатора с двойной поляризацией Meteor 50DX» (Шифр: «Метеор») под научной редакцией проф. Кузнецова А.Д. № государственной регистрации НИР: 01201251681, 2012.- 137 с.

4. Воинов Н.Е., Кузнецов А.Д. д.ф-м.н., проф., зав. каф., Симакина Т.Е. к.ф-м.н., доц. Текстурный анализ спутниковых изображений мезомасштабной облачности. // Тезисы доклада на Итоговой сессии Ученого совета РГГМУ. 2223 января 2013.

5. Войнов Н.Е., Восканян K.JL, Кузнецов А.Д., Сероухова О.С., Симакина Т.Е. Методика текущего прогнозирования эволюции зон осадков по радиолокационным данным // Тезисы доклада на заседании постоянно действующего семинара «Дистанционные методы зондирования природной среды». СПб. Военно-космическая академия имени А. Ф. Можайского. 2012

Воинов Никита Евгеньевич

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

МЕТОДЫ ЦИФРОВОЙ ОБРАБОТКИ МНОГОЗОНАЛЬНЫХ

СПУТНИКОВЫХ СНИМКОВ В ЗАДАЧАХ АНАЛИЗА МЕЗОМАСШТАБНЫХ КОНВЕКТИВНЫХ АТМОСФЕРНЫХ

ДВИЖЕНИЙ

ЛР№ 020309 от 30.19.96.

Подписано в печать 30.10.13. Формат 60x90 1/16. Гарнитура Times New Roman. Печать цифровая. Усл. печ. л. 1,5. Тираж 100 экз. Зак. № 222. РГГМУ, 195196, Санкт-Петербург, Малоохтинский пр., 98. _Отпечатано в ЦОП РГГМУ_

Министерство образования и науки Российской Федерации ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ "РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ

УНИВЕРСИТЕТ" (РГГМУ)

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

МЕЗОМАСШТАБНЫХ КОНВЕКТИВНЫХ АТМОСФЕРНЫХ

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

04201450393

ВОЙНОВ НИКИТА ЕВГЕНЬЕВИЧ

ДВИЖЕНИИ

25.00.30 - метеорология, климатология, агрометеорология

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель: доктор физ.-мат. наук, проф. Кузнецов А.Д.

Научный консультант:

кандидат физ.-мат. наук, доцент Т.Е. Симакина

РГГМУ

Санкт-Петербург 2013

ОГЛАВЛЕНИЕ

Стр.

Введение 4

1. АТМОСФЕРНЫЕ МЕЗОКОНВЕКТИВНЫЕ ПРОЦЕССЫ ПО

ДАННЫМ ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ 9

1.1. Спутниковые средства мониторинга мезомасштабных конвективных систем 9

1.2. Классификации мезомасштабных конвективных 23 систем

1.3. Конвективные системы мелкой и глубокой конвекции 30

1.3.1. Мезомасштабные системы мелкой конвекции 3 0

1.3.2. Мезомасштабные конвективные комплексы 46

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

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

аппаратуры дистанционного зондирования Земли 56

2.2. Методы колориметрических преобразований спутниковых изображений 72

2.2.1. Цвет как важная характеристика восприятия объекта 72

2.2.2. Колориметрические модели 76

2.3. Синтез псевдоцветных изображений по данным многоспектрального дистанционного зондирования 85

2.3.1. Идентификация облачности на основе синтезирования снимков с помощью различных колориметрических моделей 85

2.3.2. Идентификация стадии развития конвективной облачности 97

3. МЕТОДЫ ИССЛЕДОВАНИЯ ТЕКСТУРЫ МЕЗОКОНВЕКЦИИ НА СПУТНИКОВЫХ СНИМКАХ 104

3.1. Разномасштабный анализ на основе

вейвлет- преобразования 104

3.1.1. Теория дискретного вейвлет- преобразования 104 изображений

3.1.2. Пример использования вейвлет-преобразования для формирования пирамиды изображений 108

3.2. Текстурный анализ снимков облачности мелкой 112 конвекции

3.3. Алгоритмы определения высоты слоя конвективных

ячеек открытого и закрытого типа 118

3.4. Методика расчета текстурной сложности области,

занятой разномасштабными облачными системами, на

основе вычисления ее фрактальной размерности 122

3.4.1. Фрактальный подход к описанию текстур 122

3.4.2. Исследование мезоконвективных систем на основе вычисления фрактальной размерности 131

Заключение 142

Список использованных источников 145

ВВЕДЕНИЕ

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

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

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

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

Оперативная идентификация таких систем, а также определение стадии их «жизненного» цикла позволяет осуществлять диагноз и сверхкраткосрочный прогноз. Изучение положения, структуры и типа облачных образований позволяет проанализировать пространственно-временную иерархию мезомас-штабных систем /2-4/. Так, изменение пространственной организации и интенсивности конвективных ячеек масштабов мезо-у и -р определяет жизненный цикл более крупных долгоживущих мезомасштабных конвективных систем масштаба мезо-а/4-6/.

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

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

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

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

Цель исследования. Целью настоящей работы являются:

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

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

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

Для достижения поставленной цели автором:

1. Рассмотрены спутниковые средства мониторинга мезомасштабных конвективных систем.

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

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

4. Разработана методика расчета текстурной сложности области, занятой разномасштабными облачными системами, на основе вычисления ее фрактальной размерности.

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

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

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

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

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

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

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

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

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

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

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

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

Апробация результатов. Основные результаты работы докладывались на Итоговой сессии Ученого совета РГГМУ (Санкт-Петербург, январь 2013 г.), на научных семинарах кафедры экспериментальной физики атмосферы РГГМУ (2013 г.).

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

Структура и объем. Материал диссертации состоит из введения, трех глав, заключения, списка литературы. Текст содержит 154 страницы, включая 76 рисунков, 20 таблиц. Список использованных источников содержит 101 наименование.

1. АТМОСФЕРНЫЕ МЕЗОКОНВЕКТИВНЫЕ ПРОЦЕССЫ ПО ДАННЫМ ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ

1.1. Спутниковые средства мониторинга мезомасштабных

конвективных систем

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

Метеорологические спутники

Первыми странами, начавшими метеорологические программы, были Россия (тогда СССР) и США. Началом систематического изучения поверхности Земли из космоса можно считать запуск в 1960 г. американского метеорологического спутника Tiros-1. В 1969 г. Советским Союзом был запущен спутник Метеор серии "Космос". До 1990-х годов спутников серии "Космос" было запущено более ста. В сентябре 2009 г. Россия возобновила метеорологическую программу, осуществив запуск спутника Метеор-М.

На сегодняшний день безусловным лидером среди данных дистанционного зондирования благодаря бесплатному распространению и очень высокому радиометрическому разрешению являются данные с американского спутника NOAA (National Oceanic and Atmospheric Administration). Очередной спутник NOAA-19 был запущен в 2009 г. и сейчас в активной эксплуатации 5 космических аппаратов NOAA, что позволяет осуществлять съемку с высокой периодичностью.

В 2006 году запущен европейский метеоспутник MetOp-A, в 2008 г. -китайский спутник Fengyun 3.

В группу метеорологических спутников можно отнести космические аппараты (КА) широкого назначения, запущенные в рамках американской экологической программы ЕОС (Earth Observing System) - глобального мониторинга поверхности Земли, проводимой NASA /18/. В конце 1999 г. выведен на орбиту первый спутник программы EOS АМ-1, который известен под названием Terra. Два других спутника этой серии - Aqua и Aura - запущены в 2002 и 2004 гг. соответственно.

К спутникам метеорологического назначения относятся геостационарные спутники, международная глобальная система которых формируется с 70-х гг. прошлого века. Одними из первых были запущены два геостационарных спутника GOES (США), координаты подспутниковых точек которых 75° з.д. и 135° з.д. Запуск очередного геостационарного метеоспутника GOES-13 состоялся в 2006 г. Европейский спутник Meteosat расположен над точкой с координатой 0°. Последний спутник этой серии MSG (Meteosat второго поколения) несет на своем борту 12-тиканальный сканер SEVIRI. Раздел 2.1 настоящей работы посвящен информации спектральных каналов данной аппаратуры. Японский GMS MTSAT и индийский Insat KALPANA находятся над точками с координатами 140° в.д. и 74° в.д. соответственно. Российский GOMS-2 «Электро-JI № 1» заменил в орбитальной позиции 76° в. д. космический аппарат GOMS-1 «Электро», прекративший работу в 1998 г. Вслед за Электро-Л № 1 будут запущены аналогичные спутники GOMS-3 Электро-Л № 2 в настоящем году и GOMS-4 Электро-Л № 3 - в 2015 г. /19/.

Характерными особенностями метеорологических спутников является большой охват территории съемкой (около 2800 км и больше для геостационарных спутников) и низкое пространственное разрешение (около 1 км).

Основным средством ДЗЗ на борту метеорологических спутников является съемочная камера типа оптико-механического сканера.

Оптико-механические сканеры, схема которых представлена на рис. 1.1 /14/, включают в себя зеркальный объектив с очень малым углом зрения, устройство, отклоняющее сканирующий луч, - зеркало, систему фильтров и приемник сигналов. В каждый момент на сканирующее зеркало проецируется малый участок местности, составляющий единичный элемент будущего изображения. Зеркало поворачивается вокруг некоторой оси, обеспечивая сканирование местности в направлении, перпендикулярном движению спутника. Отраженное от сканирующего зеркала излучение попадает на систему фильтров, которая разделяет излучение на спектральные диапазоны. После этого система фотоприемников регистрирует среднюю интенсивность излучения с поверхности элементарного участка.

Рисунок 1.1. Принцип действия оптико-механического сканера

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

качающееся зеркало

фотоприемники

Ширина полосы обзора земной поверхности сканером на рис. 1.1 обозначена буквой Ь и определяется следующей формулой

L = 2H-tgam2LX,

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

Существует два класса твердотельных фоточувствительных приборов: квантовые и тепловые приемники излучения. Квантовые фотоприемники основаны на изменении электрических свойств полупроводника при п