ДИСТАНЦИОННЫЕ МЕТОДЫ
В ЛЕСНОМ ХОЗЯЙСТВЕ
DOI 10.24419/Ш12304-3083.2020.2.05 УДК 630.5
Автоматизированное выявление низкополнотных насаждений зоны смешанных лесов с применением спутниковых данных среднего разрешения на примере Учебно-опытного лесхоза БГТУ
В.В. Вечеров
«Заплеспроект», филиал ФГБУ «Рослесинфорг», ведущий инженер, кандидат сельскохозяйственных наук, г. Брянск, Российская Федерация, vecherovvv32@gmail.com
В.И. Дзубан
Управление лесами Брянской области, начальник управления,
г. Брянск, Российская Федерация,
dzu-vladimir@mail.ru
Ю.И. Перепечина
Брянский государственный инженерно-технологический университет, профессор кафедры «Лесное дело», доктор сельскохозяйственных наук, г. Брянск, Российская Федерация, lhf_bryansk@mail.ru
В.И. Шошин
Брянский государственный инженерно-технологический университет, кандидат сельскохозяйственных наук, г. Брянск, Российская Федерация, VIShoshin44@gmail.com
Относительная полнота - один из важнейших таксационных показателей, позволяющий определять запас насаждения и характеризовать его состояние.
В работе проанализирована возможность автоматизированного выявления низкополнотных насаждений, пройденных выборочными санитарными рубками, по показателю сомкнутости лесного полога (FCD) на основании данных среднего пространственного разрешения Sentinel-2.
Для ссылок: http://dx.d0i.0rg/10.24419/LHI.2304-3083.2020.2.05
Автоматизированное выявление низкополнотных насаждений зоны смешанных лесов с применением спутниковых данных среднего разрешения на примере Учебно-опытного лесхоза БГТУ/В.В. Вечеров, В.И. Дзубан, Ю.И. Перепечина, В.И. Шошин. - DOI 10.24419/LHI.2304-3083.2020.2.05. - Текст : электронный//Лесохоз. ин-форм.: электрон. сетевой журн. - 2020. - № 2. - С. 54-65. URL: http://lhi.vniilm.ru/
В результате корреляционно-регрессионного анализа связи наземных данных и данных, полученных из карт сомкнутости лесного полога, выявлена сильная экспоненциальная зависимость (r = 0,787, R2 = 0,829 при p-value < 0,001) между показателем FCD и относительной полнотой древостоя. На основании установленной зависимости созданы карты относительной полноты для лесных участков, пройденных выборочными санитарными рубками. Ключевые слова: дистанционное зондирование Земли, низкопол-нотные насаждения, вегетационный индекс, forest canopy density, QGIS
Для оценки плотности стояния деревьев в лесном хозяйстве применяют понятие «полнота древостоя», под которым понимают показатель, характеризующий степень использования деревьями занимаемого ими пространства [1]. Относительная полнота - один из важнейших таксационных показателей, позволяющий определять запас насаждения и характеризовать его состояние. Таким образом, относительная полнота древостоя выступает в качестве характеристики фактической продуктивности того или иного насаждения [2].
Существует 3 различных подхода к определению или расчёту относительной полноты насаждений: а) непосредственные натурные измерения на пробных площадях; б) статистический метод, для которого необходимы данные о сумме площадей поперечных сечений или среднем диаметре древостоя и числе стволов на 1 га; в) применение данных дистанционного зондирования Земли (ДЗЗ).
Наземные и статистические методы при их практическом применении требуют трудоёмких вычислений и значительных затрат. Именно поэтому в настоящее время особую значимость для мониторинга состояния земель лесного фонда приобрели современные и доступные методы дистанционного зондирования. Дистанционное зондирование Земли - процесс или метод получения информации об объекте, участке поверхности или явлении путем анализа данных, собранных без контакта с изучаемым объектом [3]. Преимущества данных дистанционного зондирования Земли при изучении динамики лесных сообществ объясняются, прежде всего, их «первичностью» по сравнению картографическими материалами, а также возможностью извлекать разную, в зависимости от необходимости, информацию [4]. В данной работе для оценки относительной полноты древостоев с целью выявления низкополнотных насаждений использован метод картирования по показателю сомкнутости лесного полога (Forest Canopy Density - FCD) на основании данных среднего пространственного разрешения Sentinel-2 [5]. Показатель FCD является удобным
в применении индикатором для оценки сомкнутости древесной растительности и состояния лесов при использовании данных космической съёмки [6, 7]. Преимущество этого метода заключается в возможности моделирования густоты крон (сомкнутости) без наземных тренировочных данных. Показатель густоты полога FCD выражается в процентах для каждого пикселя изображения, что позволяет, используя полевые данные и регрессионный анализ, составлять попиксельные карты относительной полноты древостоев.
Цель исследования - проанализировать возможность автоматизированного выявления низкополнотных насаждений, пройденных выборочными санитарными рубками.
Решаемые задачи:
Объекты и методы исследований
Объект исследования - Учебно-опытный лесхоз Брянского государственного инженерно-технологического университета (БГИТУ), находящийся на территории Брянского лесного массива [8]. Для Брянского лесного массива характерна большая пестрота почвенного покрова и неоднородность лесоводственно-таксационных показателей насаждений [9]. Это осложняет применение ДЗЗ для выявления насаждений, изре-женных в результате санитарно-оздоровительных
Автоматизированное выявление низкополнотных насаждений зоны смешанных лесов с применением спутниковых данных среднего разрешения на примере Учебно-опытного лесхоза БГТУ
мероприятий (СОМ). В качестве натурных объектов на территории Учебно-опытного лесхоза БГИТУ отобрано 20 лесотаксационных выделов, отличающихся породным составом насаждений, возрастом и интенсивностью изреживания (табл. 1).
Для автоматизированного выявления низкополнотных насаждений, пройденных выборочными санитарными рубками, применяли материалы бесплатной съёмки среднего разрешения Sentinel-2, которая на данный момент является наилучшей комбинацией высокого
спектрального (13 каналов), пространственного (10-20 м/пкс) и временного разрешения (5 сут для одной и той же территории)[10].
Учитывая, что лесные насаждения являются динамической экосистемой, для оценки трендов взаимосвязи показателя относительной полноты древостоя и показателя FCD использованы 3 раз-носезонные сцены Sentinel-2 - поздневесенняя (25 мая 2019 г.), летняя (14 июня 2019 г.) и ран-неосенняя (9 сентября 2019 г.).
В ходе натурных работ в мае 2019 г. заложено 20 круговых пробных площадей постоянного
Таблица 1. Таксационная характеристика пробных площадей
№ п/п Координаты центра (WGS84) X Y Квартал Выдел Состав по ярусам Возраст, лет Высота, м Диаметр, см Класс бонитета ТЛ1 ТЛУ2 Отн. полнота
№ п/п Координаты центра (WGS84) X Y Квартал Выдел Состав по ярусам Возраст, лет Высота, м Диаметр, см Класс бонитета ТЛ1 ТЛУ2 Отн. полнота
радиуса (R=17,8 м, S = 1000 м2) в насаждениях различной полноты (0,1-0,5), пройденных выборочными санитарными рубками. Характеристика пробных площадей приведена в табл. 1.
Картирование сомкнутости лесного полога по FCD основано на нескольких индикаторах: улучшенный вегетационный индекс (AVI), индекс оголённости почвы (BI) и индекс затенённости древесного полога (SI) [7, 11] (табл. 2). Расчёт вегетационных индексов проводился попиксельно с применением растрового калькулятора в QGIS.
Перед процедурой расчёта индексов была проведена атмосферная и радиометрическая коррекция съёмки в QGIS (модуль SCP), а также нормализация данных с применением линейной трансформации [12].
Улучшенный вегетационный индекс (AVI) характеризует величину биомасс и очень чувствителен к изменениям в густоте полога растительности благодаря усилению отражения в инфракрасном диапазоне, в отличие от широко применяемого вегетационного индекса [13].
Индекс оголённости почвы (BI) позволяет разделять растительность, расположенную на различных почвах, благодаря тому, что красный канал реагирует на минеральное строение почв, а ближний инфракрасный и зелёный канал характеризуют растительность [13, 14].
Индекс затенённости древесного полога (SI) позволяет оценить степень затенённости крон древостоя, что служит признаком для разделения насаждений как по возрасту (молодые насаждения имеют низкие значения индекса), так и по сомкнутости (сомкнутые насаждения имеют низкую степень затенения); выражается в процентах [15].
Показатель сомкнутости лесного полога FCD рассчитывался попиксельно в растровом калькуляторе QGIS по формуле:
FCD=[(VD х SI + 1)/2] - 1,
где VD - показатель густоты растительности (первая компонента преобразованного по методу главных компонент (PCA) композитного изображения из растров AVI и BI), % [7].
Таблица2. Индикаторы, используемые при моделировании густоты древесного полога
Индикатор Формула расчёта*
Улучшенный вегетационный индекс (AVI) AVI=[B4x256-B3)x(B4-B3)]1/3
Индекс оголённости почвы (BI) BI={[(B4+B2)-B3]/[(B4+B2)+B3]}
Индекс затенённости древесного полога (SI) SI=[(256-B1)x(256-B2)x(256-B3)]1/3
* В1 - синий канал космоснимка, В2 - зелёный канал, В3 - красный канал, В4 - ближний инфракрасный канал.
Автоматизированное выявление низкополнотных насаждений зоны смешанных лесов с применением спутниковых данных среднего разрешения на примере Учебно-опытного лесхоза БГТУ
Результаты и обсуждение
В результате проведённых исследований составлена карта сомкнутости древесного полога насаждений Учебно-опытного лесхоза БГИТУ. Согласно полученным данным (рис. 1), наибольшую долю среди древостоев с низкой сомкнутостью составляют хвойные насаждения с преобладанием ели европейской (Picea abies L.) в составе. Это связано с проведением санитарно-оздоровительных мероприятий вследствие увеличения площади насаждений, поврежденных короедом-типографом. Предыдущее исследование, посвящённое глазомерному выявлению низкополнотных хвойных насаждений после проведения СОМ, показало, что чрезмерное изреживание при выборочных санитарных рубках (ВСР) приводит к замедлению восстановления природоохранных
функций лесов и на длительное время снижает их продуктивность [16].
После создания попиксельной карты сомкнутости лесного полога в программе QGIS пробным площадям присвоены полученные значения показателя FCD. Затем в среде RStudio проведен корреляционно-регрессионный анализ связи между показателем FCD и данными натурных обследований относительной полноты древостоя:
Результат корреляционно-регрессионного анализа связи между FCD и относительной полнотой
Уравнение регрессии P=0,0024xe(°&°881xFCD)
Коэффициент корреляции, г 0,787
Коэффициент детерминации, R2 0,829
F, при p-value<0,001 10,57
Стандартная ошибка ±0,07
Рис. 1. Карта сомкнутости древесного полога (fcd, %) на территорию объекта исследований
Из-за ограниченного набора полевых натурных данных, их неравномерного распределения по территории объекта работ, а также необходимости включения в модель высокополнотных лесных насаждений, проведена рандомизированная выборка 10-ти дополнительных выделов с относительной полнотой 0,6-1,0 по векторным данным последнего лесоустройства.
Таким образом, корреляционно-регрессионный анализ связи между показателем FCD и показателем относительной полноты древостоя выполнен по 30-ти независимым переменным (рис. 2), что многие авторы [17-19] считают достаточным условием репрезентативности выборочной совокупности.
Между показателем относительной полноты древостоя и величиной сомкнутости лесного полога имеется выраженная экспоненциальная зависимость. Значения коэффициентов корреляции (г = 0,787) и детерминации = 0,829) близки к единице, что говорит о сильной зависимости относительной полноты от показателя FCD. Коэффициенты уравнения регрессии значимы на уровне р^а1ие < 0,001.
Судя по значению коэффициента детерминации = 0,829) для всех изучаемых насаждений, около 83% изменчивости относительной полноты древостоя обусловлено изменчивостью сомкнутости древесного полога, что коррелирует с увеличением показателя FCD.
Р = 0,0024хе(0,0881хрС0) R2 = 0,829
Рис. 2. График зависимости между относительной полнотой древостоев и показателем сомкнутости древесного полога (fcd, %) для территории объекта исследований
Оценка точности полученной регрессионной модели выполнена широко применяемым методом перекрёстной проверки (k-foldcгoss-validation) в среде RStudio (пакет DAAG). При его использовании имеющиеся натурные данные разбивают на к частей (в данной работе к = 5), затем на к-1 частях (т.е. на 25 независимых переменных) проводят обучение модели, а на оставшейся части - тестирование. Процедуру повторяли к раз, в итоге получили значение средней точности регрессионной модели, которое составило 0,96. Данный показатель значим при р^а1ие < 0,001, что говорит о высокой точности полученной модели.
Результаты перекрёстной проверки регрессионной модели
Уравнение регрессии Средняя точность модели F, при р^а!ие < 0,001 Стандартная ошибка
Р = 0,0024хе(0&0881хГС0) 0,96 169,00 ±0,01
После установления уравнения регрессии между показателями FCD и относительной полноты в программе QGIS с помощью модуля «Калькулятор растров» создана тематическая карта вычисленных значений относительной полноты на 20 отобранных для исследования лесотаксационных выделах. Затем в программе QGIS с помощью модуля SCP проведена неконтролируемая классификация полученной карты относительной полноты и последующая фильтрация и векторизация классифицированной тематической карты относительной полноты древостоев, пройденных выборочными санитарными рубками (табл. 3, рис. 3).
Согласно данным табл. 5, среди насаждений, пройденных выборочными санитарными рубками, наибольшая доля приходится на древостои с полнотой 0,3 (42,9%). Значительную площадь (8,45 га, или 14,9% обследованных насаждений) занимают прогалины (относительная полнота 0,1 и менее), пригодные для создания лесных культур с целью повышения продуктивности насаждений.
Автоматизированное выявление низкополнотных насаждений зоны смешанных лесов с применением спутниковых данных среднего разрешения на примере Учебно-опытного лесхоза БГТУ
Относительная полнота Площадь
<0,1 3,41 6,0
Итого 56,71 100
Рис.3. Карта относительной полноты насаждений, пройденных выборочными санитарными рубками, на территорию объекта исследований
ДИСТАНЦИОННЫЕ МЕТОДЫ В ЛЕСНОМ ХОЗЯЙСТВЕ
Выводы
Таким образом, применение спутниковых данных среднего разрешения Sentinel-2 в зоне смешанных лесов позволяет с достаточно высокой точностью определять величину относительной полноты древостоев, т.е. осуществлять автоматизированное выявление низкополнотных (0,5-0,1) насаждений, пройденных выборочными санитарными рубками.
В ходе статистического анализа данных в среде RStudio получена экспоненциальная регрессионная модель. Значения коэффициентов корреляции и детерминации (г = 0,787;
R2 = 0,829) свидетельствуют о сильной зависимости относительной полноты от показателя FCD.
Средняя точность модели, полученная методом перекрёстной проверки на натурных данных, составляет 0,96 (при р^а1ие < 0,001).
На основании установленной зависимости созданы попиксельные карты относительной полноты для лесных участков, пройденных выборочными рубками.
Согласно результатам исследования, среди насаждений, пройденных выборочными санитарными рубками, наибольшую площадь занимают древостои с полнотой 0,3 (42,9%).
Автоматизированное выявление низкополнотных насаждений зоны смешанных лесов с применением спутниковых данных среднего разрешения на примере Учебно-опытного лесхоза БГТУ
Список использованных источников
C.К. Гош. - М. : Техносфера, 2008. - 312 с.
D. Panagiotidis, P. Surovy // Journal of forest science. - 2017. - 63(3). - Р. 106-115.
References
DOI 10.24419/LHI.2304-3083.2020.2.05
Automated Identification of Low-Density Stands of a Mixed Forest Area with Use of Мedium Resolution Satellite Data on the Example of Uchebno-Opytnyj Leschoz BSUET
V. Vecherov
"Zaplesproekt", Branch of FSBI "Roslesinforg", Lead Engineer, Candidate of Agricultural Sciences, Bryansk, Russian Federation, vecherovvv32@gmail.com
V. Dzuban
Forest Management in the Bryansk region, Head of Department, Bryansk, Russian Federation, dzu-vladimir@mail.ru
Yu. Perepechina
Bryansk State University of Engineering and Technology, Professors of the Forestry Department, Doctor of Agricultural Sciences, Bryansk, Russian Federation, lhf_bryansk@mail.ru
V. Shoshin
Bryansk State University of Engineering and Technology, Candidate of Agricultural Sciences, Bryansk, Russian Federation,VIShoshin44@gmail.com
Relativelystand density is one of the most important taxation indicators, which allows determining the stock of stands and characterizing its condition. The relativelystand density of aforest as a characteristic of the actual productivity of a particular plantation.
Terrestrial and statistical methods in their practical application require laborious calculations and significant costs. That is why at present, modern and affordable methods of remote sensing have acquired special significance for monitoring the state of the lands of the forest fund. The paper analyzes the possibility of automated detection of low-density stands, passed through selective sanitary felling, according to the index of forest canopy closure (FCD) based on the average spatial resolution of Sentinel-2. The FCD indicator is a convenient indicator to use for assessing the closeness of woody vegetation and the state of forests when using satellite imagery data. The advantage of this method lies in the possibility of modeling the density of crowns (closeness) without ground training data. After a regression analysis of ground data and data obtained from forest canopy closure maps, a strong exponential dependence (r = 0.787, R2 = 0.829 with p-value < 0.001) between the FCD and the relative fullness of the stand was revealed. Because of the established dependence, maps of relative completeness were created for forest plots covered by selective sanitary felling.