Экспериментальные исследования собственного излучения пространственных структур облачности в диапазоне 8-13 мкм для их объективного распознавания

Артюхов, Александр Викторович

Бесплатный автореферат и диссертация по наукам о земле на тему
Экспериментальные исследования собственного излучения пространственных структур облачности в диапазоне 8-13 мкм для их объективного распознавания
ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Экспериментальные исследования собственного излучения пространственных структур облачности в диапазоне 8-13 мкм для их объективного распознавания"

005012486

АРТЮХОВ Александр Викторович

ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ СОБСТВЕННОГО ИЗЛУЧЕНИЯ ПРОСТРАНСТВЕННЫХ СТРУКТУР ОБЛАЧНОСТИ В ДИАПАЗОНЕ 8-13 МКМ ДЛЯ ИХ ОБЪЕКТИВНОГО РАСПОЗНАВАНИЯ

Специальность 25.00.29 - Физика атмосферы и гидросферы

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

1 2 УД?

Москва-2012 г.

005012486

Работа выполнена в Институте экспериментальной метеорологии Федерального государственного бюджетного учреждения НПО «Тайфун» Федеральной службы России по гидрометеорологии и мониторингу окружающей среды.

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

доцент Третьяков Николай Дмитриевич

Научный консультант: Заслуженный метеоролог РФ, доктор технических наук, профессор Алленов Михаил Иванович.

Официальные оппоненты:

кандидат технических наук, профессор, Соловьев Виктор Александрович,

доктор физико-математических наук, Дмитриева Лидия Романовна.

Ведущая организация: Федеральное государственное бюджетное учреждение «Научно-исследовательский центр космической гидрометеорологии «Планета».

Защита состоится 21 марта 2012 г. в 14 часов на заседании совета по защите докторских и кандидатских диссертаций Д 327.003.01 при Федеральном государственном бюджетном учреждении "Гидрометеорологический научно-исследовательский центр Российской Федерации" по адресу: 123242, Россия, Москва, Большой Предтеченский переулок, д. 11-13.

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

Автореферат разослан 20 февраля 2012 г.

Ученый секретарь диссертационного совета,

доктор географических наук Нестеров Е.С.

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

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

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

Оценки формы и балла облачности являются чрезвычайно важными и для целого ряда прикладных задач - в сельском хозяйстве, где урожайность зависит от освещенности, или в авиации для аэродромных служб, для которых форма и балл облачности являются значимыми характеристиками. Однако до сих пор определение этих характеристик производится субъективно, а в ночное время - вообще не производится. Предварительно проведенные в ФГБУ «НПО «Тайфун» натурные

"Ч

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

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

Для достижения поставленной цели требуется решить следующие задачи:

- провести анализ современного состояния методов и средств определения форм и балла облачности и методов классификации;

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

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

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

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

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

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

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

Практическое значение. Разработаны методы определения формы и балла облачности. Автоматически получаемые данные могут дополнить или заменить существующую практику определения балла и формы облаков на синоптических станциях наблюдателем. Разработан способ определения направления и скорости движения облачности. На данный способ определения направления движения облачного поля получен патент на изобретение №2414728 от 20.03.2011.

Реализация результатов. Основные результаты работы внедрены в научно-производственном объединении «Тайфун», Военной академии войсковой противовоздушной обороны Вооруженных сил Российской Федерации им. А. М. Василевского, Всероссийском научно-исследовательском институте сельскохозяйственной метеорологии. Работа выполнялась по федеральным целевым программам и программам Росгидромета «Целевая научно-техническая программа Росгидромета «Научные исследования и разработки в области гидрометеорологии и мониторинга окружающих сред», программе Союзного государства «Совершенствование системы обеспечения населения и отраслей экономики Российской Федерации и Республики Беларусь» 2007-2011 гг., Федеральной Целевой программе Российской Федерации «Мировой океан».

Апробация работы. Основные результаты докладывались на конференции молодых ученых, посвященной 70-летию дрейфа «СП-1» (Москва, 2008 г.); на

восьмой Международной конференции «Прикладная оптика-2008» (Санкт-Петербург, 2008 г.); на шестнадцатой военно-научной конференции Военной академии войсковой ПВО ВС РФ (Смоленск, 2008 г.); на Всероссийской научной конференции «Исследование процессов в нижней атмосфере при помощи высотных сооружений» (Обнинск, 2009 г.); на Международном симпозиуме стран СНГ «Атмосферная радиация и динамика» (Санкт-Петербург, 2009 г.); на конференции "Проблемы теории и практики развития войсковой ПВО в современных условиях" (Смоленск, 2010 г.); на девятой Международной конференции «Прикладная оптика-2010» (Санкт-Петербург, 2010 г.).

Защищаемые положения. На защиту выносятся:

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

2. Способ определения балла облачности по уровню энергетической яркости облаков.

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

4. Способ определения направления и скорости движения облачности на основе совместного анализа соседних по времени кадров яркости.

Личный вклад.

Личный вклад соискателя заключается:

1. В непосредственном участии при разработке и проектировании автоматизированной системы для параметризации и распознавания форм облачности (АСПРФО).

2. В выполнении работ по метрологической аттестации АСПРФО и проведении измерений пространственно-временных структур собственного излучения облачности.

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

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

5. В подготовке материалов и написании статей и совместном обсуждении результатов исследований по целевым программам («Целевая научно-техническая программа Росгидромета «Научные исследования и разработки в области гидрометеорологии и мониторинга окружающих сред», программе Союзного государства «Совершенствование системы обеспечения населения и отраслей экономики Российской Федерации и Республики Беларусь» 2007-2011 гг., Федеральной Целевой программе Российской Федерации «Мировой океан»).

6. В проведении патентных исследований и участии в разработке формулы и описании изобретения (патент на изобретение №2414728 от 20.03.2011).

7. В научно-методических и организационных работ по внедрению результатов диссертации в Военной академии войсковой противовоздушной обороны Вооруженных сил Российской Федерации им. А. М. Василевского, в подразделениях ФГБУ «НПО «Тайфун», Всероссийском научно-исследовательском институте сельскохозяйственной метеорологии.

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

Структура и объем работы. Диссертационная работа состоит из введения, четырех глав, заключения, списка литературы и приложений. Основная часть диссертационной работы изложена на 124 страницах текста, содержит 9 таблиц, 44 рисунка. Список литературы включает 97 наименований.

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

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

С помощью сканирующей аппаратуры проводились систематические исследования пространственно-временных и спектральных структур излучения облачного фона в диапазонах 3 - 5 и 8 - 13 мкм. В результате, были разработаны модели фона в диапазонах 3 - 5 и 8 - 13 мкм и получены значения ряда статистических характеристик для различных форм облачности: средние значения, дисперсии, спектральные плотности и плотности распределения флуктуаций энергетической яркости.

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

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

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

многоканальный радиометр для исследования радиационной структуры полей природных образований в спектральном интервале 1,5-2,86 мкм; низкотемпературный радиометр на область спектра 8-13 мкм; сканирующий радиометр высокого пространственного разрешения на область спектра 1,4-13 мкм;

автоматизированную систему для параметризации и распознавания форм и количества облачности (АСПРФО) созданную коллективом лаборатории оптических методов исследований природных сред ФГБУ «НПО «Тайфун» (при непосредственном участии автора диссертации) под руководством доктора техзических наук, профессора М.И. Алленова и доктора технических наук Н.Д. Третьякова и запатентованную в России (№2331853 от 20.06.2008 г).

В качестве основного средства для решения поставленной проблемы была выбрана автоматизированная система для объективной параметризации форм и

балла облачности (АСПРФО). В отличие от использованных ранее приборов, эта система позволяет оперативно получать объективную информацию о форме и количестве облачности по полусферическим изображениям поля яркости.

Система представляет собой сканирующий радиометр, обеспечивающий сканирование облачного поля по азимуту от 0 до 360° с пространственным разрешением 10 минут дуги, по зенитному углу 0 с разрешением от десятков минут до нескольких градусов. В большинстве случаев, интервал сканирования по зенитному углу составлял от 30° до 62°. Полное круговое сканирование осуществляется за 1 секунду. За это время регистрируется 360 значений энергетической яркости облачного поля через каждый градус. После завершения записи данных поступает сигнал на шаговый двигатель, и зеркало меняет угол наклона на 1°, цикл повторяется. После следующего оборота угол меняется еще на один градус. Через заданное программой количество шагов (строк) например 17, шаговый двигатель возвращает сканирующее зеркало в исходное начальное положение, цикл повторяется и записывается следующий кадр. В результате работы системы, мы получаем набор матриц, где по горизонтали 360 значений, а по вертикали — 16 (или другое, заданное количество строк). Каждое из 6120 значений представляет собой конкретную область - изображение неба в ИК области на небесной сфере.

Вольт-ваттная характеристика радиометра:

где В(Т1), В(Г2) - энергетические яркости исследуемого объекта и опорного источника излучения (внутренней полости радиометра);

А - постоянный для данного прибора коэффициент - коэффициент передачи; Т], Т2 - абсолютная температура источника излучения и внутренней полости радиометра соответственно;

Излучение объекта определяется выражением

и(Т)=А[В(Т,)-В(ТЛ

где В (Т) - энергетическая яркость объекта, Вт см'2 -ср'1;

г(Х, Т) - функция Планка;

Т- абсолютная температура источника излучения;

5' (I.) - относительная спектральная характеристика приемника;

г ф (Я) - спектральная характеристика интерференционного светофильтра.

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

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

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

- средние значения энергетической яркости (ЭЯ) по зенитным углам;

- дисперсии по зенитным углам;

- нормированные автокорреляционные функции К,(к)

где 1 - азимутальный угол,) - зенитный угол, В/ - значение энергетической яркости, ц, - средние значения по зенитным углам / а) - дисперсии по зенитным углам у,

- нормированные спектральные плотности й0

Л=| Г

где / = 0,1,..., Г

в формуле (1) ы(к) обозначает корреляционное окно с точкой отсечения I (использовалось окно Тьюки (Бендат Дж. Измерение и анализ случайных процессов,

У: Р- шаг по частоте для спектральной плотности (/•" в 2-3 раза больше Ь). - нормированные взаимные корреляционные функции р„т(к) между зенитными углами

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

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

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

1971)):

Р„,Лк) = —Ц--цп){ВшМ

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

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

В результате анализа полученных данных был сделан вывод о необходимости расширения системы классификации за счет ряда дополнительных классификационных признаков. Дополнительные признаки позволят эффективней использовать вероятностный классификатор, показывающий насколько вероятна принадлежность конкретного набора признаков (полученных по конкретной ситуации) к какому-либо выбранному классу облачности. В качестве дополнительных признаков были разработаны аналоги признаков, применявшихся ранее в похожей задаче классификации облачности (Garand L., Automated Recognition of Oceanic Cloud Patterns, 1987). Эти признаки включают: общий процент облаков, процент облаков на нижнем, среднем и верхнем уровнях, индекс многоярусности, количество отдельных облаков и просветов, связность облаков и просветов, общая направленность облаков.

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

09:01:40 - 16:37:11

О 200 400 600 600 1 ИЛ 1 200 1 400 Кадры

Рис. 1. Сравнение балла облачности. Кучевая и высококучевая облачность. 17.09.09. 09.00- 16.30

3:55:43 - 17:01:41

Рис. 2. Сравнение балла облачности. Кучевая и слоисто-кучевая облачность. 06.10.09.09.00 - 17.00 На рисунках 1 и 2 представлены примеры результатов определения балла облачности для различных дней и различных типов облачности во времени (по

кадрам). На каждом графике два набора точек: значение балла облачности, определяемое штатным в начале каждого часа метеорологом (на графиках обозначен как «оператор») и вычисленный балл облачности.

В таблице 1 приведены статистические данные правильности определения балла облачности для различных форм. В качестве исходных данных использовались кадры, полученные весной и осенью 2009 года и в августе 2010 года. Значения балла облачности вычислялись для интервала +/- 5 минут от начала каждого часа, и усреднялись. Полученный балл сравнивался со значением, указанным метеорологом для этого момента времени. Среднее абсолютное отличие величины балла облачности показывает усредненную (по количеству отсчетов) сумму отличий балла облачности от величин указанных метеорологом. Средняя разность величины балла облачности суммирует отличия балла облачности с учетом знака.

Таблица 1. Статистика определения величины балла облачности

Тип облачности Общее количество отсчетов с данной облачностью Кол-во отсчетов с только с данной облачностью Среднее абсолютное отличие величины балла облачности Средняя разность величины балла облачности

Си 119 32 2.3 0

Ас 134 17 2.1 1

Se 108 16 1.8 0

St 24 23 0.2 0

c¡ 154 38 2.7 -1

Сс 42 0 2.4 1.2

Cb 27 0 2.2 0.5

Из таблицы видно, что лучше всего балл облачности определяется для кучевых (Си), слоистых (81) и слоисто-кучевых (8с) облаков. Для перистых облаков (С1) характерно занижение общего балла облачности, а для остальных типов облачности (Ас, Сс, СЬ) - завышение балла.

Помимо общей характеристики количества облаков, для классификации также

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

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

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

Определим 1к(1) и ИЬф как число ячеек принадлежащих соответственно облаку / или просвету Также, упорядочим все области облаков и просветов по размеру - от наиболее мелких к более крупным областям. Тогда связность облаков (СС) определяется как:

(Иф)/(тСР), СИ*О (1, СР = 0 '

где т - количество облаков, СГ- коэффициент покрытия облаками (от 0 до 1), к это такой номер облака, что неравенство

^с(1)>тСР72, / = 1,к удовлетворено при суммировании размеров участков облаков, упорядоченных по размеру. Аналогично, связность просветов (ВС) определяется как:

Гад/М'-с/?, СР*1

и Ср = 1

где т - количество просветов, к соответствует просвету, такому, что неравенство

ХЬЬ(,-)>т(1-СГ)/2, у = Скверно при суммировании размеров участков просветов упорядоченных по размеру.

Таким образом, в приведенных выше равенствах Ис(к) и ИЬ(к') являются средними размерами облаков и просветов, соответственно.

Рис. 3. Пример бинарного изображения кадра. Связность облаков наиболее полезна при разделении кучевых и слоисто-кучевых облаков. Значение параметра связности меньше для кучевых облаков - сильно разрозненные элементы, и выше для слоисто-кучевых облаков, которые связаны в более крупные структуры. Признак связности просветов является хорошим средством определения ячеистых структур.

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

1 ¡^ (1ги /V ^

ЛА и=о /=о V Л Ь ]

и = й,\,...К-\, V = 0,1,... ¿-1 Такой спектр симметричен относительно начала координат, причем конфигурация поля облачности влияет на конфигурацию спектра мощности, а в случаях безоблачного неба и неба, полностью закрытого облаками, вся энергия спектра яркости обычно сосредоточена около начала координат. Исходя из этого, можно определить признак направленности как абсолютное значение коэффициента корреляции, взвешенной мощности спектра излучения (8Т): вТ = |5ху| / (5х8у)"2

где

и V и V

и = 0,1,...8-1, V = 0,1,...8 — 1, (и,\') * (0,0), РБ(к, у) - мощность спектра: Р8(г/,у) = |Р(«,у)|2

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

классификатор вида:

К = ДМ S,= £,[(x-M)(x-Al)r]

где/>,(х) - вероятность принадлежности вектора признаков х (принадлежащего конкретному кадру), классу i, Е - оператор среднего значения, X, - ансамбль кадров принадлежащих классу /; d - количество признаков или размерность, //, - вектор средних значений размерности d, и симметричная матрица ковариаций размера d х d. Правило выбора класса / для вектора х есть максимизация значения вероятности р,(х). Такой классификатор полностью определяется значениями средних векторов признаков и матриц ковариаций, причем в отличие от простейшего классификатора, определяющего метрическое расстояние между вектором признаков и средним вектором каждого класса, этот классификатор учитывает корреляционные связи между различными признаками с помощью матрицы ковариаций.

Для проверки работы такого классификатора было создано два набора кадров. Первый набор предназначался для вычисления средних значений признаков каждого класса, а второй - для оценки успешности полученных значений. Каждый набор содержал 500 экземпляров, разбитых на 9 классов, представляющих собой различные формы облачности, в том числе, и несколько смешанных форм. Общее количество правильно распознанных кадров составило 70.4% (352 кадра из 500). Если ограничиться только одноярусными формами облачности, то количество правильно классифицированных кадров составит 85.4%

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

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

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

||1>; с '">§! В Л*

* № ^ " ю

Рис. 4. Изменение соседних, по времени, кадров.

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

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

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

21.09.09 14:00-15:00 Си - б баллов, вычисленная скорость 8.1 м/с

21.06.07 12:20-14:00 Си - 5 баллов, вычисленная скорость 10 м/с

Рис. 5. Сравнение результатов определения направления движения облачности.

--линия векторов вероятного движения облачности.

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

22.09.09 15:00- 16:00

21.10.09 16:00- 17:00

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

--линия векторов вероятного движения облачности.

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

На рисунках 5 и 6 приведены примеры определения направления движения облачности в сравнении с метеорологическими данными о скорости ветра. Ломаная линия указывает направление ветра, согласно вышеописанному методу. Прямая линия - направление ветра на высоте 300 м, по данным высотной метеомачты г.Обнинска в рассматриваемые интервалы времени.

Зная интервал времени Д1, в течение которого производились измерения, и координаты вектора (х, у), в декартовой системе определяются составляющие скорости Ух и Уу, и затем скорость V и направление движения нижней границы облачности ср:

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

Несмотря на заметное визуальное отличие от кадров, полученных в дневное время

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

В 3

15.10.09 00:00-01:00

16.10.09 00:00-01:00

Рис. 7. Примеры определения направления движения в ночное время суток.

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

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

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

Отмечено, что для кучевой облачности средние значения флуктуаций яркости для всех отмеченных баллов облачности увеличиваются с ростом угла визирования, причем скорость увеличения среднего значения энергетической яркости к горизонту обратная увеличению балла облачности: она больше у меньшего балла (1-3). Дисперсии флуктуаций энергетической яркости для облачности 4-6 и 7-9 баллов

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

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

Структура поля излучения высококучевых облаков более мелкая, чем у кучевых и слоисто-кучевых облаков. Взаимные корреляционные связи тесны по всему полю Ас. Пространственные спектры мощности флуктуаций яркости излучения спадают медленнее, чем спектры для Си и Sc, показатель SAl(6) даже для углов визирования 0=40...50° не превышает 1,8, тогда как в случае Си и Sc он больше 2.

Вычислены средние значения энергетической яркости слоистых (St) облаков в весенне-осеннее время. Они изменяются в интервале энергетической яркости (3,1-3,7)-10"3 Вт • см"2 • ср"1, радиационная температура Трад всегда положительна и близка к абсолютной температуре Та5с, т.к. коэффициент излучения этой формы облачности близок к единице. В структурах полей просвечивающихся As trans облаков наблюдаются оптические неоднородности поля яркости размером от 3 до 15 угловых градусов, причем неоднородности иногда имеют волнистую структуру, но чаще стохастически распределенные структуры.

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

составил 70%. Если ограничиться только одноярусными формами облачности, то количество правильно классифицированных кадров составит 85%

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

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

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

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

(* - работы, опубликованные в изданиях, рекомендованных ВАК Министерства образования и науки РФ)

1. Алленов М.И., Артюхов A.B., Иванов В.Н., Третьяков Н.Д. Структура излучения неба при смешанных формах облачности в диапазоне 8-13 мкм // Прикладная оптика - 2008: Сб. материалов Международной конференции, 20-24 октября 2008 г. - С.-Пб., 2008. - С. 263-267.

2. Алленов М.И., Артюхов AB., Иванова Н.П., Третьяков Н.Д. Состояние фоновых полей при кучевой (Си) и смешанной (Ac+Cu) формах облачности в диапазоне 8-13 мкм // Шестнадцатая военно-научная конференция Военной академии войсковой ПВО ВС РФ: Сб. научных материалов, ч. 3. - Смоленск, 2008. -С.110-114.

3. Артюхов A.B., Алленов М.И., Гирнык Е.Б., Иванова Н.П., Новиков H.H., Третьяков Н.Д. «Об угловом (по зенитному углу) распределении флуктуаций собственного излучения мощно-кучевых облаков, полученном по трехмерным изображениям // Исследование процессов в нижней атмосфере при помощи высотных сооружений: Тез. докладов Всероссийской научной конференции, 8-10 октября 2008 г. - Обнинск, 2008. - С. 134-136.

4. Артюхов A.B., Алленов М.И., Гирнык Е.Б., Иванов В.Н., Иванова Н.П., Новиков H.H., Третьяков Д.Н., Третьяков Н.Д. Экспериментальные исследования пространственно-временных структур собственного излучения смешанных форм облачности в диапазоне 8-13 мкм // Исследование процессов в нижней атмосфере при помощи высотных сооружений: Тез. докладов Всероссийской научной конференции, 8-10 октября 2008 г.-Обнинск, 2008.-С. 137-141.

5. Артюхов A.B., Алленов М.И., Новиков H.H., Третьяков Н.Д. Об угловых (по зенитному углу) распределениях флуктуаций радиационной температуры мощных

кучевых (Cb) облаков в диапазоне 8-13 мкм, полученных по трехмерным изображениям // Атмосферная радиация и динамика: Тез. докладов Международного симпозиума стран СНГ, 22 - 26 июня 2009 г. - С.-Пб., 2009. - С. 96.

6. Артюхов A.B., Алленов М.И., Кондратюк В.И., Третьяков Н.Д., Третьяков Д.Н. Противоизлучение неба при смешанных формах облачности в диапазоне 8-13 мкм // Атмосферная радиация и динамика: Тез. докладов Международного симпозиума стран СНГ, 22 - 26 июня 2009 г. - С.-Пб., 2009. - С. 97.

7. Алленов М.И., Артюхов A.B.,, Иванов В.Н., Третьяков Н.Д.. Официальный бюллетень федеральной службы по интеллектуальной собственности, патентам и товарным знакам №8 от 20.03.2011. Способ определения направления и скорости движения нижней границы облачности: Патент №2414728 от 20.03.2011.

8. Алленов М.И., Артюхов A.B., Третьяков Д.Н. «Угловые зависимости излучения облачности и морского горизонта». «Проблемы теории и практики развития войсковой ПВО в современных условиях. Материалы XVII военно-научной конференции.», Смоленск: ВА ВПВО ВС РФ, 2010. Ч. 3, стр. 191-195.

9. Алленов М.И., Артюхов A.B., Третьяков Д.Н., Третьяков Н.Д. «Экспериментальные исследования фонового излучения горизонта в диапазоне 8-13 мкм». - «Проблемы теории и практики развития войсковой ПВО в современных условиях. Материалы XVII военно-научной конференции.», Смоленск: ВА ВПВО ВС РФ, 2010. Ч. 3, стр. 195-198.

10. Алленов М.И., Артюхов A.B., Иванов В.Н., Третьяков Д.Н., Третьяков Н.Д. «Фоновые помехи при слоистых формах облачности». - «Проблемы теории и практики развития войсковой ПВО в современных условиях. Материалы XVII военно-научной конференции.», Смоленск: В А ВПВО ВС РФ, 2010. Ч. 3, стр. 181186.

11. Алленов М.И., Артюхов A.B., Иванов В.Н., Третьяков Н.Д. «Способ определения направления и скорости движения фонового поля собственного излучения в диапазоне 8-13 мкм». - «Проблемы теории и практики развития

войсковой ПВО в современных условиях. Материалы XVII военно-научной конференции.», Смоленск: ВА ВПВО ВС РФ, 2010. Ч. 3, стр. 186-191.

12. Артюхов A.B., Третьяков Н.Д., Якименко И.В. "Определение балла и формы облачности на основе векторов признаков" // Математическая морфология. Электронный математический и медико-биологический журнал. Т. 9. - Вып. 2, 2010.

13. Алленов М.И., Алленов A.M., Артюхов A.B., Иванова Н.П., Новиков H.H., Третьяков Н.Д., Якименко И.В. "Радиационные характеристики пространственных структур собственного излучения форм облачностей в диапазоне 8-13 мкм." // Математическая морфология. Электронный математический и медико-биологический журнал, 2010. Т. 9. - Вып. 3.

14. Алленов М. И., Артюхов А. В., Иванова Н. П., Новиков Н. Н., Третьяков Н. Д. «Радиационные характеристики пространственных структур собственного излучения форм облачностей в диапазоне 8-13 мкм». «Проблемы гидрометеорологии и мониторинга окружающей среды». Том 1. Сборник статей/Под ред. А.Ф. Нерушева/Юбнинск 2010. Стр. 123-137.

15. Алленов М. И., Артюхов А. В., Васильев A.C., Давлетшина P.A., Иванов В. Н., Третьяков Д. Н., Третьяков Н. Д.. «Нанотехнологии интерференционных систем ГУ «НПО «Тайфун» для спектрорадиометрии природных сред». «Проблемы гидрометеорологии и мониторинга окружающей среды». Том 1. Сборник статей/Под ред. А.Ф. Нерушева/Юбнинск 2010. Стр. 269 - 287.

16. Алленов М.И., Алленов A.M., Артюхов A.B., Иванова Н.П., Новиков H.H., Третьяков Д.Н., Третьяков Н.Д. «Диагноз и прогноз противоизлучения неба при слоистых формах облачности в диапазоне 8-13 мкм». Труды оптического общества. X международная конференция «Прикладная оптика - 2010». С.Петербург 18-22 октября 2010г. Том 1, стр. 255-259.

17. Алленов М.И., Алленов A.M., Артюхов A.B., Иванова Н.П., Новиков H.H., Третьяков Д.Н., Третьяков Н.Д. «Корреляционные пространственные связи между флуктуациями энергетической яркости кучевых (Си), слоисто-кучевых (Sc) и смешанных форм облачности в диапазоне 8-13 мкм». Труды оптического общества.

X международная конференция «Прикладная оптика - 2010». С.Петербург 18-22 октября 2010г. Том 1, стр. 260-264.

18. Алленов М.И., Артюхов A.B., Иванов В.Н., Третьяков Н.Д., Якименко И.В. «Способ определения направления и скорости движения облачности по структуре собственного излучения в диапазоне 8-13 мкм» // Математическая морфология. Электронный математический и медико-биологический журнал, 2010. Т. 9. - Вып. 4.

19. Артюхов A.B. Определение формы и балла облачности в дневное и ночное время при помощи сканирующего полусферического радиометра // Тезисы докладов конференции молодых специалистов, Обнинск, ГУ «НПО «Тайфун», 24-26 ноября 2010 г. С. 16-19.

20. Алленов М.И., Артюхов A.B., Иванов В.Н., Третьяков Н.Д. Способ определения направления и скорости движения облачностей. // Тезисы докладов конференции молодых специалистов, Обнинск, ГУ «НПО «Тайфун», 24-26 ноября 2010 г. С. 271-274.

21*. Алленов М.И., Артюхов A.B., Третьяков Д.Н., Третьяков Н.Д. Исследование противоизлучения облачного неба в дневное и ночное время в диапазоне 8-13 мкм. "Оптический журнал", 78, 9, 2011 г. С. 20-24.

22*. Алленов М.И., Артюхов A.B., Третьяков Д.Н. Диагноз и прогноз структур собственного излучения форм облачности для мониторинга их загрязнений. Вестник ТГТУ, 2011 г. Том 17. № 2 С. 327-332.

23*. Артюхов A.B., Третьяков Д.Н., Якименко И.В. Распознавание форм облачности на основе векторов признаков. Журнал «Наукоемкие технологии», №6, 2011, т.12. С. 74-76.

24. Алленов М.И., Артюхов A.B., Иванов В.Н., Третьяков Н.Д. Способ определения направления и скорости движения облачностей. Журнал «Метеорология и гидрология», 2011 (в печати)...

Подписано к печати 08.02.2012 г. Формат 60x84/16. Усл.п.л. 1,6.

Тираж 100 экз. Заказ № 2 Отпечатано в ФГБУ «ВНИИГМИ-МЦ Д», г. Обнинск, Королева, 6

Текст научной работыДиссертация по наукам о земле, кандидата физико-математических наук, Артюхов, Александр Викторович, Москва

61 12-1/606

Федеральное государственное бюджетное учреждение «Научно-производственное объединение «ТАЙФУН»

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

Артюхов Александр Викторович

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

Специальность 25.00.29 - Физика атмосферы и гидросферы

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

Третьяков Н.Д.

Научный консультант Заслуженный метеоролог РФ доктор технических наук, профессор

Алленов М.И.

Москва-2012 г.

ОГЛАВЛЕНИЕ

Введение................................................................................................................................4

1. Обзор научных результатов по исследованию пространственно-временных структур облачности в ИК-диапазоне спектра...........................................................16

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

2.1. Средства измерения пространственно-временных структур полей излучения (отражения) природных сред.........................................................................................27

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

2.1.2 Малогабаритный низкотемпературный радиометр на область.....................30

спектра 8-13 мкм..........................................................................................................30

2.1.3 Сканирующий радиометр высокого пространственного разрешения на область 1,4-13 мкм.......................................................................................................34

2.2. Автоматизированная система для параметризации и распознавания форм и балла облачности (АСПРФО)........................................................................................37

2.3. Методика измерений................................................................................................55

2.4. Методика обработки................................................................................................57

2.5 Выводы по второй главе...........................................................................................64

3. Классификатор различных форм облачности на основе векторов прнзнаков.65 3.1 Исследования пространственно-временных структур собственного излучения облачных полей в диапазоне 8-13 мкм.......................................................................65

3.1.1. Радиационная модель неба при кучевой (Си) облачности............................65

3.1.2. Радиационная модель неба при слоисто-кучевых (8с) облаках...................67

3.1.3 Радиационная модель неба при высококучевой (Ас) облачности.................69

3.1.4. Радиационные характеристики слоистых (81:), высоко-слоистых (Аб) и перисто-слоистых (Сб) облаков в диапазоне от 8 до 13 мкм..................................71

3.1.5. Пространственно-временная структура высокослоистой Аб облачности...72 3.1.5. Пространственно-временная структура слоистой облачности ................74 У

3.1.7. Характерные пространственные структуры полей собственного излучения

облачности....................................................................................................................78

3.1.8. Угловое (по зенитному углу) распределение флуктуаций собственного

излучения мощно-кучевых облаков, полученном по трехмерным изображениям. 79

3.2. Вектор признаков и метод классификации облачности.......................................82

3.2.1 Коэффициент покрытия (балл облачности).....................................................83

3.2.2 Индекс многоярусности.....................................................................................89

3.2.3 Признаки связности............................................................................................89

3.2.4 Спектральные признаки.....................................................................................90

3.3. Методология классификации и результаты..........................................................91

3.4. Выводы по третьей главе.........................................................................................94

4. Методы исследований направления и скорости движения облачности по ее пространственно-временной структуре собственного излучения в диапазоне 813 мкм..................................................................................................................................96

4.1. Существующие методы и средства измерений направления и скорости

движения облачности.....................................................................................................96

4.2. Метод определения направления и скорости движения облачности на основе

попарного анализа соседних, по времени, кадров.....................................................102

4.3. Выводы по четвертой главе..................................................................................110

Основные результаты и выводы по работе..............................................................111

Список использованной литературы:........................................................................115

Приложение.....................................................................................................................127

Введение

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

Для достижения поставленной цели требуется решить следующие задачи:

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

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

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

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

Международных конференциях «Прикладная оптика-2008» и «Прикладная оптика-2010» (Санкт-Петербург, 2008 г., 2010 г.); на шестнадцатой военно-научной конференции Военной академии войсковой ПВО ВС РФ (Смоленск, 2008 г.); на Всероссийской научной конференции «Исследование процессов в нижней атмосфере при помощи высотных сооружений» (Обнинск, 2009 г.); на Международном симпозиуме стран СНГ «Атмосферная радиация и динамика» (Санкт-Петербург, 2009 г.); на конференции "Проблемы теории и практики развития войсковой ПВО в современных условиях" (Смоленск, 2010 г.)*; на девятой Международной конференции «Прикладная оптика-2010» (Санкт-Петербург, 2010 г.).

Защищаемые положения. На защиту выносятся:

2. Способ определения балла облачности по уровню энергетической яркости облаков.

Личный вклад.

Личный вклад соискателя заключается:

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

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

/■УШ'- Л./ V"

С 'Г5 ■') -

Содержание.

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

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

Систематические исследования стохастических пространственно-временных и спектральных структур излучения облачного фона в диапазонах 3 - 5 и 8 - 13 мкм проводились с помощью сканирующей аппаратуры. Были разработаны модели фона в диапазонах 3 - 5 и 8 - 13 мкм. Основной задачей исследования являлось выявление маломасштабных флуктуаций, порождаемых мелкими облачными неоднородностями. Были получены значения ряда статистических характеристик для различных форм облачности: средние значения, дисперсии, спектральные плотности и плотности распределения флуктуаций энергетической яркости.

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

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

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

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

Информация о работе

Артюхов, Александр Викторович
кандидата физико-математических наук
Москва, 2012
ВАК 25.00.29

Диссертация

Экспериментальные исследования собственного излучения пространственных структур облачности в диапазоне 8-13 мкм для их объективного распознавания - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы