Обсудить в форуме Комментариев 9Редактировать в вики
Авторы: Андрей Шихов, Ян Маракулин
Содержание |
Летом 2010 года на Европейской территории России и в Пермском крае наблюдались катастрофические лесные пожары. По территории Пермского края были получены оценки пройденных пожарами площадей от Агентства по Природопользованию Пермского края, МЧС, ИКИ РАН и Гринпис России, значительно различающиеся между собой. Этим была обусловлена постановка первой задачи работы – создание векторного слоя гарей 2010 года на территорию Пермского края, на основе космических снимков среднего (LANDSAT 5/TM) и высокого (SPOT-4) разрешения.
Вторая задача – оценка степени повреждения лесной растительности на участках, пройденных пожарами летом 2010 года. Оценка проведена по вегетационным индексам, данные полевых наблюдений не использовались. Поэтому методика оценки степени повреждения растительности будет дорабатываться.
В качестве исходных данных использовались:
Участки лесной растительности, пройденные пожарами, характеризуются пониженной спектральной яркостью в ближней инфракрасной зоне (далее по тексту - NIR). Это объясняется уменьшением содержания хлорофилла в вегетативных органах усыхающих деревьев. Также для гарей характерно повышение спектральной яркости в средней инфракрасной зоне (далее по тексту - SWIR). Оно, в свою очередь, объясняется уменьшением содержания влаги в листьях или хвое. В видимой зоне спектра для гарей характерна более высокая, чем у здоровой растительности, спектральная яркость. Это также объясняется уменьшением содержания хлорофилла, которое внешне проявляется в дефолиации и дехромации листьев при усыхании деревьев[3].
На перечисленных выше закономерностях основаны различные вегетационные индексы, а также комбинации каналов, позволяющие дешифрировать гари. Разработаны также автоматизированные алгоритмы оконтуривания гарей. Однако для их эффективного применения необходимо, чтобы гарь была полностью свободна от облачности и теней[4].
Методика дешифрирования гарей включала следующие этапы:
- Отбор снимков LANDSAT 5-TM за период до и после зафиксированных пожаров.
- Расчет вегетационных индексов:
Индекс NDVI = (NIR–RED)/(NIR+RED)
Индекс SWVI = (NIR–SWIR)/(NIR+SWIR), где SWIR - 5-й канал снимка LANDSAT 5-TM. Индекс DNDVI = NDVIpre – NDVIpost, где NDVIpre –индекс за период до пожара, NDVIpost –индекс за период после пожара.
Индекс DSWVI = SWVIpre – SWVIpost
Также для всех пар снимков рассчитана степень изменения яркости в ближнем инфракрасном канале с помощью инструмента Change Detection (метод PCA).
Дешифрирование гарей проводилось по следующим признакам:
Примеры использования дешифровочных признаков приведены ниже на рис. 1 - 6
Рис. 1. Визуализация разновременных снимков LANDSAT 5-TM в синтезе каналов SWIR, NIR, GREEN (7-4-2). Пермский край, Гайнский район, Веслянское лесничество. Выделены две гари площадью более 8000 га.
С использованием всех вышеперечисленных признаков проведено ручное оконтуривание гарей. Всего нам удалось выделить 46 гарей 2010 года на общей площади 31959 га, что существенно больше, чем по данным АПП ПК. При этом площадь отдельных выделенных нами гарей составляет от 6 га до 4245 га. Выделено 20 крупных гарей площадью свыше 200 га, причем на них приходится более 90% всех выгоревших площадей, из них 7 гарей площадью более 1000 га.
Суммарная площадь выделенных нами гарей более чем на 7 тыс. га превышает площадь гарей, выделенных Гринпис России в 2010 году[5]. Такое расхождение объясняется тем, что в Гринпис России использовали снимки за период до завершения пожароопасного сезона.
Было также выполнено сопоставление полученного нами слоя гарей с результатами детектирования выгоревших площадей по снимкам низкого разрешения MODIS, проведенного ИКИ РАН [2]. Площадь гарей, выделенных ИКИ РАН на территории Пермского края, составляет 27969 га. Однако значительная часть объектов в данном слое не совпадает по расположению с тепловыми аномалиями MODIS, то есть ошибочно классифицирована как гари. На юге Пермского края обнаружена «ложная» гарь площадью более 7300 га. Площадь гарей, выделенных ИКИ РАН и совпадающих с нашими данными, составляет менее 16 тыс. га.
По полученным данным рассчитана площадь выгоревших территорий по участковым лесничествам. Установлены даты действия пожаров, с которыми связаны выделенные гари. Также были сопоставлены площади гарей от крупных пожаров по нашим данным и по данным АПП ПК, примеры приведены в Табл. 1.
Таблица 1. Сопоставление данных о крупных пожарах, полученных по ДДЗЗ и предоставленных АПП ПК
Номер пожара (АПП ПК) | Участковое лесничество | Площадь ликвидации пожара (АПП ПК), га | Площадь оцифрованной по снимку гари, га | Период действия (FIRMS) | Период между обнаружением и ликвидацией пожара (АПП ПК) | Максимальная зарегистрированная температура в 21 канале, К (FIRMS) |
256 | Вайское | 976 | 977+140 | 04.08.2011 – 15.08.2011 | 05.08.2011 – 22.08.2011 | 320,9
|
342 | Велсовское | 100 | 336 | 10.08.2011 – 11.08.2011 | 13.08.2011 – 14.08.2011 | 342 |
347 | Вишерское | 1120 | 3212 | 07.08.2011 – 15.08.2011 | 14.08.2011 – 19.08.2011 | 356 |
339 | Кутимское | 380 | 805+360 | 10.08.2011 – 16.08.2011 | 12.08.2011 – 25.08.2011 | 371,4
|
314, 320 | Ульвинское, Камское | 1207+428 | 3168 | 25.07.2011 – 08.08.2011 | 10.08.2011 – 23.08.2011 | 346,7 |
247 | Верхнекамское, Дозовское | 2100 | 3998 | 02.08.2011 – 15.08.2011 | 04.08.2011 – 24.08.2011 | 391,4 |
Таким образом, наши оценки площадей некоторых крупных пожаров значительно (в три раза и более) превышают оценки АПП ПК. В то же время по данным АПП ПК зафиксировано 23 крупных пожара (более 200 га), нами же выделено только 20 крупных гарей.
Для выборочной валидации контуров гарей использованы панхроматические и мультиспектральные снимки SPOT4, за вегетационный период 2011 года. Более высокое разрешение этих снимков позволяет определить контур гари с большей точностью. В то же время снимки SPOT4 за период до пожароопасного сезона у нас отсутствуют, что не позволяет использовать разностные вегетационные индексы. Валидация контуров гарей проведена методом визуального дешифрирования с использованием синтезов каналов nir-red-green и swir-nir-green, а также панхроматического канала высокого разрешения (10 м).
За вегетационный период 2011 года территории, в наибольшей степени пострадавшие от пожаров (западная часть Гайнского района) ни разу не были покрыты качественными малооблачными снимками SPOT4. Поэтому валидацию контуров удалось провести только для 9 гарей, которые закрыты безоблачными фрагментами снимков.
Выборочная валидация позволяет оценить точность оконтуривания гарей. Полученные расхождения в оценке площади гарей по данным LANDSAT5-TM и SPOT4 составляют от 11 до 66 га. Относительная ошибка оценки площади минимальна для наиболее крупных гарей (ошибка не более 5%). Для гарей площадью менее 100 га ошибка оценки площади может достигать 25% и более.
Пример валидации контуров гари по данным SPOT4 приведен на рис. 5
Выделенные крупные гари являются сложными объектами, включающими в себя различные типы местности – лесные и нелесные (криволесье, горная тундра) земли, относительно свежие рубки, фрагменты болот, долины рек. В пределах контуров гарей значительную площадь занимают участки растительности, относительно слабо пострадавшие от огня (сохранившие достаточно высокую спектральную яркость в NIR канале). К ним относятся, прежде всего, фрагменты болот, долины малых рек и участки, поврежденные низовым пожаром малой интенсивности. Для их выделения необходима оценка степени повреждения растительности в пределах выделенных гарей. Она может быть реализована с применением вегетационных индексов и методов классификации изображений.
Оценка степени повреждения растительности на всех крупных гарях (площадью более 200 га) была проведена по методике, предложенной в работе )[6]. Авторы данной работы предложили ориентировочную переклассификацию значений NDVI (после пожара) в степень повреждения растительности (Табл. 2). Для более корректной переклассификации необходимо проведение полевых обследований выгоревших участков с целью получения калибровочной зависимости NDVI и других вегетационных индексов (SWVI, dSWVI) от степени повреждения древостоя пожаром.
После переклассификации были рассчитаны площади выделенных классов для каждой крупной гари. Результаты, суммированные по всем крупным гарям площадью более 200 га, приведены в таблице 2. Пример результата переклассификации для гари от пожара № 247 (нумерация АПП ПК) приведен на рис. 6.
Таблица 2. Перевод NDVI в степень повреждения растительности
Значение NDVI | Класс поражения | Площадь, га |
Менее 0,0 | 4 (отсутствие растительности, открытая почва после пожара) | 123 |
0,0 - 0,1 | 3 | 3201 |
0,1 - 0,2 | 2 | 4467 |
0,2 - 0,3 | 3 | 7324 |
Более 0,3 | 0 (неповрежденная растительность, или участки восстанавливающейся растительности) | 13797 |
К классу неповрежденной или восстанавливающейся растительности были отнесены 13797 га площади крупных гарей, то есть почти половина всей площади. Часть этой территории (прежде всего болота) по-видимому, вообще не была пройдена пожарами, но оказалась полностью окружена выгоревшими участками. Другие фрагменты, несомненно, пострадали от огня, на что указывают высокие значения разностного индекса DSWVI (0,15 – 0,2). Но скорее всего это были низовые пожары небольшой интенсивности.
Участки, относящиеся к 2, 3 и 4 классу повреждения, вероятно, были пройдены верховыми пожарами и требуют лесовосстановительных мероприятий. На ряде участков облачность, попадающая в границы гарей, сделала невозможной корректную классификацию. Так как NDVI облаков и из теней незначительно ниже NDVI гари, в отдельных случаях применялось маскирование облаков и теней. Для маскирования облаков мы использовали 1-й канал LANDSAT, для маскирования теней – 5-й канал.
Анализ значений разностных индексов показал, что в пределах площади гарей они имеют, как правило, положительные значения, в диапазоне от 0,1 до 0,5. Максимальные значения DNDVI (0,45…0,55) характерны для участков выгоревших мелколиственных лесов на старых вырубках. В отдельных случаях индекс DNDVI в пределах гари имеет отрицательные значения, что связано с наличием облаков или теней в пределах контура гари на снимке, сделанном до пожара.
Индекс DSWVI лучше выделяет гари, чем DNDVI, поскольку характеризуется большей устойчивостью. На территориях с неповрежденной растительностью его значения обычно близки к нулю, при условии отсутствия сезонных различий в состоянии растительности. Поэтому легко выделить все участки, где произошли значительные изменения растительного покрова (гари, свежие рубки, ветровалы). Участки с максимальными значениями DSWVI (0,4…0,6) в пределах гари визуально дешифрируются как наиболее поврежденные пожаром.
Таким образом, в результате данной работы нами создан векторный слой гарей 2011 года в Пермском крае, выполнено сопоставление площадей крупных пожаров, полученных по ДДЗЗ и по данным АПП ПК. Проведена оценка степени повреждения растительности на крупных гарях.
В перспективе мы планируем использовать более корректные методы оценки состояния растительности на выгоревших площадях, применяя в комплексе различные вегетационные индексы и методы классификации космических снимков.
Скачать слой можно здесь.
Формат файла – ESRI SHAPE
Система координат – WGS 1984
Структура атрибутивной таблицы
FID – идентификатор
SHAPE – геометрия
Scene1 – полное название сцены LANDSAT5-TM, использованной для расчета вегетационных индексов до пожара. Название включает дату снимка, например 2011193 – 2011 год, 193-й день года, т.е снимок сделан 12.07.2011г.
Scene2- полное название сцены LANDSAT5-TM, по которой оконтуривалась гарь
Area – площадь гари, оконтуренной по снимку LANDSAT5-TM
Обсудить в форуме Комментариев 9Редактировать в вики
Последнее обновление: 2014-05-14 21:49
Дата создания: 13.01.2012
Автор(ы): Андрей Шихов, Ян Маракулин
© GIS-Lab и авторы, 2002-2021. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов. (подробнее).