УДК 551.521.3
Ю.Я. Матющенко, В.Е. Павлов
Программный комплекс селекции данных AERONET для определения радиационных параметров облаков
Y.I. Matyushchenko, V.E. Pavlov
The Program Package of the AERONET Data Selection to Determine the Cloud’s Radioactive Parameters
Представлено программное обеспечение для селекции наблюдательных данных сети ЛБЯОКБТ по яркости неба в альмукантарате Солнца с целью отбора ситуаций разрывной облачности. Рассмотрены критерии отбора, приведены блок-схемы алгоритмов, интерфейс и описание работы компонентов комплекса.
Ключевые слова: атмосфера, селекция, аэрозоль, яркость, облачность.
The paper presents software package for selection of the AERONET observed data of sky brightness in almancantar of the Sun to pick out situations of fractured cloudiness. Selection criteria are considered, the block diagrams of algorithms, interface and operation specification of complex components are presented.
Key words: atmosphere, selection, aerosol, brightness, cloud.
Введение. Изучение влияния разрывной облачности на характеристики радиационных полей - актуальная проблема в связи с исследованиями в области общей циркуляции атмосферы, теории климата, метеорологии, а также при решении широкого круга прикладных задач атмосферной оптики [1-4].
Глобальная автоматизированная сеть мониторинга атмосферы ЛБЯОКБТ развернута для получения и обработки данных с целью создания карты распределения аэрозоля по земному шару [5]. Одной из фундаментальных геофизических проблем, при решении которой могут быть использованы эти данные, является построение региональных аэрозольных моделей атмосферы. Специалистами МЛ8Л проводится селекция с целью исключения измерений яркости при сплошной и разрывной облачности. Однако анализ данных на сайте ЛБЯОКБТ свидетельствует о наличии облачных образований.
Ранее авторами были предложены методические разработки для отбора безоблачных ситуаций, пригодных для построения аэрозольных моделей атмосферы и корректного восстановления оптических параметров аэрозольных частиц [6, 7]. Написано и апробировано соответствующее программное обеспечение [7-9]. Разработаны методы селекции ситуаций разрывной облачности, результаты которой могут быть использованы для строго количественного анализа оптических свойств облаков [10].
В настоящей работе представлен программный комплекс [11], который обеспечивает селекцию данных ЛБЯОКБТ для последующего определения радиационных параметров облаков.
Критерии безоблачного неба в альмукантарате Солнца. Отбор безоблачных ситуаций выполняется по следующим критериям:
1. «Мягкий» критерий, суть которого в том, что необходимым условием отсутствия облаков является систематическое убывание яркости в интервале углов рассеяния от 2 до 90° и ее возрастание при углах более 1200.
2. «Жесткий» критерий заключается в проверке гладкости угловой зависимости градиента яркости и позволяет обнаруживать малоконтрастные облачные образования (с яркостью выше яркости безоблачного неба на 1-2%).
3. Суть «дополнительного» критерия в том, что яркость в симметричных относительно солнечного вертикала точках при отсутствии облака должна быть одинакова.
Критерии отбора облачных ситуаций. Селекция ситуаций разрывной облачности включает несколько этапов.
1. Изначально обработке подвергаются данные, не прошедшие селекцию безоблачного неба, т.е. на первом этапе для отбора облачных ситуаций используются описанные выше критерии безоблачности.
2. С помощью «мягкого критерия» выбираются такие кривые яркостей, для которых наблюдается неравномерность убывания яркости с увеличением угла рассеяния в какой-либо одной полуплоскости углов азимута, что говорит о наличии разрывной облачности.
3. Отбираются кривые яркости, для которых величина контраста не менее 0,3 хотя бы в одном угле азимута данной кривой. Контрасты отображают различие в величинах симметричных яркостей и вычисляются по формуле
к(Т) = |(В(Т) - В(360° - Т)| /Вт1П, где Вт1П - меньшее значение яркости В(Т) либо В(360° - Т); Т - угол азимута.
4. Отыскиваются кривые яркостей по длинам волн, принадлежащие одной серии измерений. Наличие полной серии в массиве облачности, плюс тот факт, что облачность зафиксирована в синей и красной областях спектра, значительно повышают достоверность обнаружения облачных образований.
5. Проверка соответствия величины контраста значению не менее 0,3 хотя бы для одного угла азимута найденной полной серии измерений.
Выбор средства разработки. В качестве инструмента для реализации программного комплекса было выбрано приложение Microsoft Excel со встроенным языком программирования Visual Basic for Application (VBA). Причины данного выбора следующие:
- использование распространенного, по сути, стандартного приложения;
- наличие простого встроенного языка программирования;
- некритичность выполняемой задачи селекции к скорости обработки данных;
- удобство визуального контроля результатов в табличной форме;
- возможность оперативного применения функций электронной таблицы;
- простота организации интерфейса пользователя.
Состав и описание компонентов программного комплекса. В состав программного комплекса входят программы:
- чтения исходных данных;
- селекции безоблачных данных;
- селекции ситуаций явной облачности;
- вычисления промежуточных и выходных массивов;
- вычисления окончательных результатов селекции.
Программа чтения исходных данных. Основная задача этого компонента - выбор и открытие файла альмукантарата и вывод его в табличной форме на лист Microsoft Excel. Оперативность и удобство при выборе файлов обеспечены использованием диалогового окна открытия файлов, что значительно облегчает и ускоряет поиск. Для настройки диалогового окна с фильтрацией файлов *.alm. использован метод GetOpenFileName объектной модели Microsoft Excel [12]. Для того чтобы команда «Открыть» вызывала не процедуру открытия файла, а процедуру преобразования и вывода текстового файла на лист Excel, был задействован метод OpenText класса Workbooks [13, 14].
Программа селекции безоблачных данных AERONET подробно описана в [7, 9]; отличие данной версии в том, что она написана на языке VBA.
Программа селекции ситуаций явной облачности. Входными данными для этой программы являются два массива: исходных кривых яркостей и кривых безоблачного неба. Селекция выполняется отдельно для каждой длины волны, так как интерес представляют прежде всего цветовые характеристики исследуемых рядов яркости.
При отборе кривых, не совпадающих по дате и времени измерения с безоблачными кривыми, из исходного массива выбираются значения яркостей, соответствующие углам азимута массива безоблачного неба, исключая область ореола. Диапазон углов азимута массива контрастов определяется от начального значения до минимума яркости кривой, который лежит в пределах 90 -г- 120 градусов, соответственно, диапазон углов контрастов 10 -г- 90 градусов.
Следует отметить, что, согласно представленному алгоритму, отбраковке подлежат и кривые яркостей, для которых систематическое убывание имеет место в обеих полуплоскостях, т. е. соответствующих «мягкому» критерию безоблачности. Наличие таких кривых объясняется тем, что они не прошли проверку в программе селекции безоблачных данных по одному из оставшихся критериев - «жесткому» либо «дополнительному».
Блок-схема алгоритма программы приведена на рисунке 1.
Программа вычисления промежуточных и выходных массивов. Блок-схема алгоритма программы представлена на рисунке 2.
В данной блок-схеме полная серия измерений означает, что в серии присутствуют все четыре кривых (по длинам волн). Как показывает анализ достаточно большого числа файлов AERONET, фотометр проходит полный круг в альмукантарате Солнца за время чуть больше минуты (время сканирования небосвода на одной длине волны), а полная серия из четырех проходов занимает около 5 минут, зенитный угол при этом изменяется не более чем на 1 градус. По этим данным выполняется сведение массивов соответствующих длин волн в единый с полными сериями измерений. В процессе вычислений результаты в виде промежуточных и выходных массивов выводятся на листы книги MS Excel.
При анализе контрастов программа оставляет в массиве только те полные группы контрастов (т.е. контрасты для 4-х кривых полной серии измерений), значения которых не менее 0,3 для всей группы хотя бы в одном азимутальном угле.
Программа вычисления окончательных результатов селекции. Исходными данными для программы являются выходные массивы разностей симметричных яркостей, чистого неба, разрывной облачности, полученные программой вычисления промежуточных и выходных массивов, предварительно отфильтрованные для одной длины волны.
После объединения данных в массивы по группам азимутальных углов 10° -г- 30°, 40° -г- 60°, 70° -г- 90° для каждого из указанных массивов вычисляются среднее геометрическое яркостей и среднеквадратичные отклонения среднего геометрического яркостей в четырех длинах волн. Окончательные результаты выводятся в табличной форме и в виде диаграмм на лист MS Excel.
Рис. 1. Блок-схема программы селекции ситуаций явной облачности
■И
Конец массива
Нет
Рис. 2. Блок-схема программы вычисления промежуточных и выходных массивов
Интерфейс пользователя программного ком- - фильтрация массивов по длине волны;
плекса. Основные требования к интерфейсу: - возможность просмотра результатов в таблич-
- оперативный поиск и выбор исходных файлов ном и графическом виде;
AERONET; - сохранение результатов в файле Microsoft
- удобный запуск программ нажатием кнопки на Excel в формате *.xls; панели инструментов;
Таблица 1
Окончательные результаты вычислений
Разности яркостей
длина волны G,44 G,68 G,87 1,G2
группы углов азим ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и -
10-30 11,11994 6,972971 5,37G818 G, 854333 G,793378 G,33G59 G, 315626 G,3235G8 G,3G5129 6,66233 2,797347 4, GG1797 G,923353 G,256973 G, 422945 G,81G96 G,235353 G,382517 4,2G48G7 G,447235 G,4G4239 3,275494 2,2G9785 2,392137 G,379G36 G,339724 G,19G647 G,1755G5 G,224446 G,2G5194
40-60 2,457599 G,2G5433 G,189585
70-90 2,963536 G,3G8963 G,279793
Чистое небо
длина волны G,44 G,68 G,87 1,G2
группы углов азим ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и -
10-30 4,335721 2,627421 4,169789 G,79797 G, 673935 G,279981 G,253G19 G,37574 G, 344681 2,837298 1,472145 1,485318 G,478742 G,4G9625 G,1G7512 G,G78297 3,168717 G,325G32 G,294794 2,2GGG54 G,661165 G,379577 G,25G866 G,225188 G,G33632 G,G17326
40-60 G,115982 1 ,G81351 G,G48684 G,G46586 G,G35435
70-90 G,G82654 G,663713 G,G264G7 G,G25396 G,G18155
Облака
длина волны G,44 G,68 G,87 1,G2
группы углов азим ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и - ср. геом. ср.кв.откл. + и -
10-30 18,18824 11,24877 1G,8G582 1,353284 1,259567 G,422174 G, 4G69G2 G,461396 G, 4425G2 1G, 8G261 4,91G36 6,145128 1,337785 G,335181 G,4721G1 1, 19G371 G,313763 G,438419 7,884855 G,7492G2 G,684191 5,94G6G3 3,156G19 2,987512 G,6G69G4 G,55G648 G,22G222 G,2G5858 G,231313 G,21469
40-60 3,9GGG19 G,236595 G,223G63
70-90 4,GG6168 G,317854 G,294489
B і і .1 .■
0.4 0.5 0.6 0.7 0.8 О.Э 1 1.1
Рис. 3. Цветовые характеристики массива разностей симметричных яркостей
Рис. 4. Цветовые характеристики массива чистого неба
Рис. 5. Цветовые характеристики массива облаков
Загружаемый файл программы Microsoft Excel содержит базовый набор именованных листов, на которые выводятся при открытии либо после обработки соответствующие массивы данных. Оперативность в работе обеспечивается пользовательской панелью инструментов, которая содержит многокаскадное меню для запуска программ - компонентов комплекса, очистки указанных листов, сохранения рабочей книги Excel.
Селекция данных AERONET с использованием программного комплекса. Работа комплекса была проверена на реальных данных глобальной сети AERONET. Проведена селекция облачности для нескольких географических точек планеты. В таблице 1 и на рисунках 3-5 показаны окончательные результаты вычислений в табличном и графическом виде для географической точки Tahiti (Тихий океан).
Библиографический список
1. Лившиц Г.Ш. Рассеяние света в атмосфере. - Алма-Ата, 1973.
2. Аэрозоль и климат / под ред. К.Я. Кондратьева. -Л., 1991.
3. Изменение климата. Научные аспекты [Электронный ресурс]. - URL: http://junkscience.com/ipcc_tar/index-
2.html
4. Физика атмосферных аэрозольных систем [Электронный ресурс]. - URL: http://window.edu.ru/window/ library?p_rid=30057
5. Goddard Space Flight Center [Электронный ресурс]. - URL: http://aeronet.gsfc.nasa.gov. - Яз. англ.
6. Матющенко Ю.Я., Ошлаков В.К., Павлов В.Е. О селекции данных AERONET. Ч. 1: Обоснования методик // Оптика атмосферы и океана. - 200б. - Т. 19, №4.
7. Матющенко Ю.Я. Методы селекции и коррекции экспериментальных данных AERONET по яркости неба: автореф. дис. ... канд. физ.-мат. наук. - Барнаул,
2007.
S. Матющенко Ю.Я., Павлов В.Е. Программный комплекс селекции данных яркости при мониторинговых
наблюдениях: св-во о гос. регистрации программы для ЭВМ №2008611921 от 17.04.2008.
9. Матющенко Ю.Я., Павлов В.Е. Программный комплекс селекции данных AERONET // Известия АлтГУ. -
2008. - №1.
10. Матющенко Ю.Я., Павлов В.Е. Селекция данных AERONET для исследования оптических параметров облаков // Аэрозоли Сибири. XVII Рабочая группа. - Томск, 2010.
11. Матющенко Ю.Я., Павлов В.Е. Программный комплекс определения радиационных параметров облаков при мониторинговых наблюдениях: св-во о гос. регистрации программы для ЭВМ №2009616558 от 25.11.2009.
12. Интернет Университет информационных технологий. Биллиг В.А. Основы офисного программирования и язык VBA [Электронный ресурс]. - URL: http://www.intuit.ru/ department/office/vba2000/12/1.html.
13. Колесов А., Павлова О. Интеграция VBA в бизнес-приложения независимых разработчиков [Электронный ресурс]. - URL: http://www.compress.ru.
14. Лекции для изучения VBA (Visual Basic for Applications) [Электронный ресурс]. - URL: http://www.vbahelp.ru.