Раздел IV. Математические методы искусственного интеллекта
УДК 528.854:004.932.2
ВТ. Бондур, И.А. Матвеев, АЛ. Мурынин, А.Н. Трёкин
РАСПОЗНАВАНИЕ ВЫГОРЕВШИХ ТЕРРИТОРИЙ НА МУЛЬТИСПЕКТРАЛЬНЫХ ИЗОБРАЖЕНИЯХ С АДАПТИРУЕМОЙ МАСКОЙ ОБЛАЧНОСТИ*
Важной задачей дистанционного зондирования Земли является оценка площади, пройденной огнём в результате природных пожаров. Для этого можно применить два метода: распознавание пожаров как таковых и распознавание гарей. Выделение видимых пожаров достаточно надёжно, но далеко не все пожары видны из-за облаков. Распознавание гарей может быть произведено для любого участка поверхности, в момент, когда нет облаков, но это распознавание менее точно. Предлагается сочетание двух способов, в том числе путём применения данных о гарях, полученных первым методом, для построения классификатора во втором методе. Необходимо эффективно отсеивать области изо, . использования маски облачности, автоматически сгенерированной на основании спутниковых снимков MODIS.
Мультиспектральные изображения; распознавание природных пожаров; маска об.
V.G. Bondur, I.A. Matveev, A.B. Murynin, A.N. Trekin
RECOGNITION OF BURNED AREAS IN MULTISPECTRAL IMAGES USING
ADAPTIVE CLOUD MASK
Estimation of areas burned as result of wildfires is an important problem in Earth remote sensing. Two approaches can be exploited: recognition of fire itself and recognition of burned areas. Recognition of visible fires is quite reliable however not every fire is visible due to clouds. Recognition of burned areas can be done for any point of surface in the moment when it is not covered with clouds, but the precision is low. Combination of the two methods is proposed, particularly by using data on fires obtained by the first approach to construct the classifier to be used in second method. It is necessary to effectively discriminate images areas covered with clouds. Peculiarities of using MODIS cloud mask in these tasks are studied.
Multispectral images; wildfire recognition; cloud mask.
Введение. Определение площади природных пожаров - важная задача, решение которой необходимо при учёте экономических и экологических последствий пожаров [1]. Современные системы дистанционного зондирования Земли (ДЗЗ) позволяют получать снимки низкого и среднего разрешения (250-1000 м)
1-2 раза в сутки для каждой точки поверхности. Рассматриваются два основных , ( ): -щадей детектированных активных пожаров (дапее ДАП) за исследуемый период и определение выгоревших территорий путём распознавания их на изображениях. Первый метод хорош тем, что активный пожар достаточно надёжно обнаруживается на орбитальных снимках в инфракрасном диапазоне, и объединив такие об-
ласти, можно получить выгоревшую площадь. Однако, сведения об активных пожарах обладают существенной неполнотой в связи с тем, что области АП могут быть замаскированы облаками, дымом других пожаров и т.п. Тем не менее, рано , -
, -
, .
* Работа выполнена при финансовой поддержке РФФИ (проект № 12-07-00778-а).
гарей относится малая надёжность детектирования гарей, которые могут довольно слабо отличаться от водоемов, вспаханной земли или других слабо отражающих поверхностей. Таким образом, два метода распознавания дополняют друг друга.
, -
ности, соответствующие ДАП, аккумулированным за предыдущий период (с вы, ), -пользованы для обучения классификатора, применяемого для обработки этого же изображения. Модель различения ППО и "не-ППО", построенная по самому изображению, может оказаться более удачной, чем обобщённая модель.
. -
бражения поверхности, поставляемые MODIS. Для работы используются 22 спек, ( тепловое излучение). 2 из этих каналов имеют разрешение 250 метров, 5 - разрешение 500 метров, а остальные 1 километр. Изображения более грубых масштабов
билинейно интерполируются до 250-метрового. Таким образом, каждая точка изображения, соответствующая приблизительно 250-метровая области земной по-
([)\N
верхности, представляется 22-мерным вектором. Обозначим \x j;=1 - обучающая
1 N
выборка, состоящая из векторов, ^Х) = — ^ XXT - средний вектор, X - матрица,
N i=1
составленная из нормализованных векторов-колонок XX) X) , XXT - матрица
, 22 22. , -
Лпах
ное и минимальное собственные значения этой матрицы —>> 1, т.е. векторы
Amin
.
(уменьшения размерности) был использован метод главных компонент [3]. Число используемых главных компонент выбиралось из расчета, что собственные числа
использованных компонент дают не менее 99,9 % нормы матрицы ковариаций, т.е.
к 22
выбиралось минимальное к : ^ А2 > 0,999^ А2. Типичные значения к , полу-
j=1 j=1
ченные на данных - от 7 до 10. Таким образом, размерность задачи понижается в
2-3 .
Множество точек гарей моделировалось многомерным нормальным распределением p(X) = Cexp((X — р)ТЕ—1 (X — р)), где Е и р - ковариационная матрица и среднее подмножества обучающей выборки, состоящее из точек гарей. Таким же образом моделировалось множество "не-гарь". Обучающая выборка для класса «гарь» выбирается согласно карте ДАП, зарегистрированных на предыдущие периоды вплоть до нескольких дней до момента регистрации обрабатываемого изображения, а для « » - , . -ласть на несколько порядков больше сгоревшей, точки второго класса берутся со случайным интервалом, равномерно распределенным на отрезке [0; 100000], т.е. в среднем каждая пятидесятитысячная. Для распознавания используется байесовский клас,
как g (x ) = (X — р1 )т Е—1 (X — р1)— (X — р2 )т Е—1 (X — р2)+ C - разделяющая функция, C - константа, выбираемая из желаемого соотношения ошибок первого и второ.
. -
, , -
ность не соответствует ни гарям, ни поверхности земли. Для учёта таких точек
использован продукт MOD35 Cloud Mask, который генерируется для каждого на-
MODIS
облаками каждого пикселя, построенную с помощью алгоритмов автоматической обработки изображения [2]. Классификатор облачности MODIS выдаёт 4 градации : 00 - , 01 - , 10 - ,
11 - . " -"
,
и результаты работы мало отличаются от применения только первого подхода (ак): . -
этому приходится использовать также области с частичной облачностью. Для более точного разделения был введён третий класс, условно называемый "^ман".
250 , -
рой выделены участки выгоревшей земли, области невыгоревшей земли и облака.
Маска облачности содержит для каждого пикселя 48-битную запись, из которой для наших целей извлекаются следующие данные: флаг вычисления маски
(1 ); (2 ): 00 - , 01 - -
мнительно ясно, 10 - вероятно ясно, и 11 - уверенно ясно; флаги «день/ночь», наличия солнечного блика, наличия снега/льда - (по 1 биту); флаг «земля/вода» (2 ): 00 - , 01 - , 10 - , 11 - ; помехи (плотный аэрозоль), тонких перистых облаков, облачной тени. Берутся в рассмотрение только те пиксели, для которых маска облачности показывает, что: для данного пиксела она вычислена; на снимке день, нет льда, солнечного блика; флаг «земля/вода» принимает значение «11», «10» или «01»; нет аэрозоля, перистых облаков, облачной тени; значение облачности «11», то есть «о»; индиви-250- , -
ками. Этот подход предполагает, что не рассматриваются точки, на которых маска MODIS -
ности. В данный момент описанный алгоритм проходит экспериментальную проверку на сериях изображений MODIS, полученных в лето 2010 года, в качестве обучающей выборки использована серия шейпфайлов, содержащих активные пожары за этот период.
Первоначально планировалось не проверять маску облачности на флаги тени, перистых облаков и аэрозоля, а также использовать те пиксели, значение облачности для которых равно «10», т.е. «вероятно, ясно». Но данный подход уже в ходе предварительных тестов показал свою неприемлемость, так как значительная , , -, .
. :
MODIS и классификатора гарей (представляющего собой байесовский с многомерным нормальным распределением) построена комбинированная система рас,
результатами работы второго, причём обучение второго происходит на данных работы первого детектора для текущего обрабатываемого изображения. Это позволяет адаптивно настраивать классификатор гарей и компенсировать различия в условиях съёмки. Помимо классов «гарь» и «не гарь», вводится класс «туман», обучающую выборку для которого можно строить на зонах, отбракованных мас.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. . . -
жаров в России // Вестник ОНЗ РАН, NZ11001. - 2010 (а). - Т. 2.
2. Ackerman S.A., Strabala K.I., Menzel W.P., Frey R.A., Moeller C.C., Gumley L.E., Baum B.A., Schaaf C., Riggs G. 1997: Discriminating clear-sky from cloud with MODIS algorithm theoretical basis document (MOD35). EOS ATBD web site. - 125 p.
3. Айвазян C.A., Бухштабер B.M., Енюков КС., Мешалкин Л.Д. Прикладная статистика: классификация и снижение размерности. - М.: Финансы и статистика, 1989.
4. .. . . .
применение. Вып. 2. - М.: Наука, 1989.
Статью рекомендовал к опубликованию д.ф.-м.н. В.И. Цурков.
Бондур Валерий Григорьевич - Научно-исследовательский институт аэрокосмического мониторинга <^эрокосмос»; e-mail: [email protected]; 105064, г. Москва, Гороховский пер. 4; тел.: +74956321654; директор; академик РАН; д.т.н.; профессор.
Матвеев Иван Алексеевич - Федеральное государственное бюджетное учреждение науки Вычислительный центр им. А.А. Дородницына Российской академии наук; e-mail: [email protected]; 119333, г. Москва, ул. Вавилова, 40; тел.: 84991352228; отдел сложных ; . ; . .- . .
Мурынин Александр Борисович - e-mail: [email protected]; тел.: +79266902722; старший научный сотрудник; к.т.н.
Трёкин Алексей Николаевич - Московский физико-технический институт (государственный университет); e-mail: [email protected]; 141700, Московская область, г. Долгопрудный, Институтский переулок, 9; тел.: 84954084554; отдел сложных систем; факультет управления и прикладной математики; студент.
Bondur Valery Grigoryevich - Research Institute for Aerospace Monitoring “Aerocosmos”; e-mail: office@ aerocosmos.info; 4, Gorokhovskij, Moscow, 105064, Russia; phone:
+74956321654; director, academician of Russian Academy of Sciences; dr. of eng. sc.; professor.
Matveev Ivan Alexeevich - Institution of Russian Academy of Sciences Dorodnicyn Computing Centre of RAS; e-mail: [email protected]; 40, Vavilov street, Moscow, 119333, Russia; phone: +74991352228; the department of complex systems; head of sector; cand. of phis.-math. sc.
Murynin Alexander Borisovich - e-mail: [email protected]; 40 phone: +79266902722; senior researcher; cand. of eng. sc.
Trekin Alexej Nikolaevich - Moscow Institute of Physics and Technology; e-mail: [email protected]; 9, Institutskii per., Dolgoprudny, Moscow Region, 141700, Russia; phone: +74954084554; the department of complex systems; student.
УДК 621.391.25
М.В. Стремоухое, С.В. Чистяков
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ФОРМИРУЮЩЕЙ СИСТЕМЫ, УЧИТЫВАЮЩАЯ НЕЛИНЕЙНЫЕ СТАТИСТИЧЕСКИЕ ЗАВИСИМОСТИ
ИЗОБРАЖЕНИЙ
Рассматривается проблема повышения качества при эффективном кодировании изображений с учетом нелинейных статистических зависимостей между их элементами. Предлагается способ параметрического кодирования изображений на основе математической модели нелинейной формирующей системы при заданной возбуждающей последовательности, обеспечивающей наилучшее по выбранному критерию совпадение синтезированного изображения с исходным. Авторы указывают на достоинства предлагаемого под-