CC BY
4
ГЕОИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ
Лесохозяйственная информация.
2024. № 1. С. 67-80 Forestry Information. 2024. № 1. P. 67-80
Научная статья УДК 630.2 EDN HMMVJL
DOI 10.24419/LHI.2304-3083.2024.1.05
Методы лесотаксационного дешифрирования радиолокационной спутниковой съемки в горных лесах Южной Сибири
Виктор Михайлович Сидоренков1
кандидат сельскохозяйственных наук
Даниил Олегович Астапов2 Юлия Сергеевна Ачиколова3 Любава Анатольевна Бадак4
Людмила Викторовна Буряк5
доктор сельскохозяйственных наук
Елена Алексеевна Шикалова6
Олег Викторович Рябцев7
кандидат сельскохозяйственных наук
Дмитрий Юрьевич Капиталинин8
Аннотация. В работе рассматриваются методы лесотаксационного дешифрирования данных радиолокационной съемки со спутника БвпИпв-г. Применяемые алгоритмы лесотаксационного дешифрирования основаны на комплексном подходе к анализу территории с учетом особенностей роста и продуктивности лесных насаждений, а также определения степени воздействия различных факторов на рост леса и лесообразовательные процессы. Разработанные методы лесотаксационного дешифрирования основаны на математических моделях, верифицированных на значительном объеме таксационных данных сети пробных площадей.
Ключевые слова: модели лесотаксационного дешифрирования, радиолокационные спутниковые данные БепИпе-г, сеть пробных площадей, РЛИ, УЭПР, карта, верхняя граница леса, класс бонитета, полнота, запас насаждений. Для цитирования: Сидоренков В.М., Астапов Д.О., Ачиколова Ю.С., БадакЛ.А., БурякЛ.В., Шикалова Е.А., Рябцев О.В., Капиталинин Д.Ю. Методылесотаксаци-онного дешифрирования радиолокационной спутниковой съемки в горных лесах Южной Сибири. - Текст: электронный //Лесохозяйственная информация. 2024. № 1. С. 67-80. йй! 10.24419/1Н!.2304-308з.2024.1.05. https://elibrary.ru/hmmvjl.
1 Всероссийский научно-исследовательский институт лесоводства и механизации лесного хозяйства, заместитель директора (Пушкино, Московская обл., Российская Федерация), lesvn@yandex.ru
2 Всероссийский научно-исследовательский институт лесоводства и механизации лесного хозяйства, заведующий лабораторией лесной таксации и лесоустройства (Пушкино, Московская обл., Российская Федерация), astdann09@yandex.ru
3 Всероссийский научно-исследовательский институт лесоводства и механизации лесного хозяйства, ведущий инженер отдела лесоводства и лесоустройства (Пушкино, Московская обл., Российская Федерация), pipintook@yandex.ru
4 Научный центр оперативного мониторинга Земли АО «Российские космические системы», начальник сектора 2 (Москва, Российская Федерация), badakla@ntsomz.ru
5 Центр лесной пирологии, филиал Всероссийского научно-исследовательского института лесоводства и механизации лесного хозяйства, главный научный сотрудник (Красноярск, Российская Федерация), buryaklv@firescience.ru
6 Объединенная дирекция государственного природного биосферного заповедника «Саяно-Шушенский» и национального парка «Шушенский бор», заместитель директора по научной работе (Шушенское, Красноярский край, Российская Федерация), e.shikalova@mail.ru
7 Всероссийский научно-исследовательский институт лесоводства и механизации лесного хозяйства, заведующий отделом инновационных технологий (Пушкино, Московская обл., Российская Федерация), sectorles@yandex.ru
8 Департамент лесного хозяйства по Центральному федеральному округу, начальник департамента (Пушкино, Московская обл., Российская Федерация), kapitalinin.d@yandex.ru
Original аrticle EDN HMMVJL
DOI 10.24419/LHI.2304-3083.2024.1.05
Methods of Radar Satellite Imagery Interpretation for Southern Siberia Mountain Forests
Victor M. Sidorenkov1
Candidate of Agricultural Sciences
Daniil O. Astapov2 Iuliia S. Achikolova3 Lyubava A. Badak4
Lyudmila V. Buryak5
Doctor of Agricultural Sciences
Elena A. Shikalova6
Oleg V. Ryabtsev7
Candidate of Agricultural Sciences
Dmitriy Yu. Kapitalinin8
Abstract. The paper discusses methods of radar imagery data interpretation from the Sentinel-1 satellite for obtaining forest characteristics. The applied algorithms of forest taxation interpretation are based on an integrated approach to the territory analysis considering features of forest growth and productivity, and various factors' influence degree on forest growth and formation processes. The developed methods of satellite imagery interpretation include mathematical models verified on a massive pull of forest taxation data from a network of trial plots.
Key words: forest taxation interpretation models, radar imagery data, Sentinel-1, specific radar cross-section, map, treeline, forest class, forest density, stock volume. For citation: Sidorenkov V., Astapov D., Achikolova Y., Badak L., Buryak L., Shikalova E., Ryabtsev O., Kapitalinin D. Methods of Radar Satellite Imagery Interpretation for Southern Siberia Mountain Forests. - Text: electronic// Forestry Information. 2024. № 1. P. 67-80. DOI 10.24419/LHI.2304-3083.2024.1.05. https://elibrary.ru/hmmvjl.
1 Russian Research Institute for Silviculture and Mechanization of Forestry, Deputy Director (Pushkino, Moscow region, Russian Federation), lesvn@yandex.ru
2 Russian Research Institute for Silviculture and Mechanization of Forestry, Head of Laboratory for Forest Inventory and Management (Pushkino, Moscow region, Russian Federation), astdann09@yandex.ru
3 Russian Research Institute for Silviculture and Mechanization of Forestry, Lead Engineer of Department for Silviculture and Forest Management (Pushkino, Moscow region, Russian Federation), pipintook@yandex.ru
4 Scientific Center for Operational Earth Monitoring of Russian Space Systems JSC, Head of Sector 2 (Moscow, Russian Federation), badakla@ntsomz.ru
5 Center for Forest Pyrology, Branch Russian Research Institute for Silviculture and Mechanization of Forestry, Chief Researcher (Krasnoyarsk, Russian Federation), buryaklv@firescience.ru
6 Joint Directorate of the State Natural Biosphere Reserve «Sayano-Shushensky» and the National Park «Shushensky Bor», Deputy Director for Scientific Work (Shushenskoye, Krasnoyarsk Territory, Russian Federation), e.shikalova@mail.ru
7 Russian Research Institute for Silviculture and Mechanization of Forestry, Head of the Innovative Technologies Department (Pushkino, Moscow region, Russian Federation), sectorles@yandex.ru
8 Forestry Department for the Central Federal District, Head of Department (Pushkino, Moscow region, Russian Federation), kapitalinin.d@yandex.ru
Введение
Разработка методов лесотаксационного дешифрирования по данным спутниковой съемки для труднодоступных территорий горных лесов является важным аспектом общей концепции технологического развития России. Рассматриваемые в статье технологии лесотаксационного дешифрирования данных радиолокационной съемки позволяют реализовать автоматизированные процессы обработки информации с учетом современных методов анализа информации на основе систем машинного обучения.
Важность решаемых вопросов отражена в основных стратегических и отраслевых документах в сфере развития лесного хозяйства России [1, 2]. Проведенные ранее исследования показывают, что применение методов дешифрирования данных спутниковой съемки для горных лесов имеет значительное различие с аналогичными работами в условиях равнинных лесов [3-5].
К одной из основных сложностей обработки данных радиолокационной съемки в горных лесах относится необходимость устранения погрешностей и искажений, вызванных спецификой отражения радиолокационного сигнала в горной местности.
В ряде научных работ показана возможность использования данных с радиолокационного спутника Sentinel-1 для определения количественных и качественных характеристик лесов [6-8], преимущество которых состоит в том, что за короткий период времени можно получить изображения с разных углов обзора. В большинстве проанализированных исследований данные радиолокационной спутниковой съемки используют с целью определения биомассы насаждений, мониторинга изменения границ земель, занятых лесной растительностью. Лесотаксационное дешифрирование радиолокационной съемки - новое направление, для разработки методов которого требуются значительные объемы данных полевых исследований. Важным этапом при разработке методов является их апробация с целью определения точности моделей и пригодности данных.
Приведенные методы лесотаксационного дешифрирования спутниковой съемки обоснованы
их верификацией по таксационным данным сети из 169 пробных площадей.
Цель исследования - изучить возможность использования радиолокационных данных в решении задач лесотаксационного дешифрирования по определению количественных и качественных характеристик насаждений для эксплуатационных и резервных горных лесов.
Объекты и методы исследований
Наземные исследования для сбора экспериментального материала выполнены в Южно-Сибирской горной зоне: в Алтае-Саянском горно-таежном, Алтае-Саянском горно-лесостепном и Байкальском горном лесных районах. Исследование охватило все выделенные В.Н. Смагиным с соавт. [9] в горах Южной Сибири лесорастительные области: Алтае-Саянскую горную, Центральноазиатскую котловинно-горную, Восточнотувинско-Южноза-байкальскую горную и Прибайкальскую горную.
В Алтае-Саянской горной лесорастительной области пробные площади с отбором образцов были заложены в Северной Алтае-Саянской горной лесорастительной провинции пихтовых и кедровых лесов (Джебашско-Амыльском округе черневых и горно-таежных пихтовых и кедровых лесов и Осевом Западносибирском округе горно-таежных и подгольцово-таежных кедровых лесов), в Восточнокузнецко-Минусинской кот-ловинно-горной лесорастительной провинции лиственничных и сосновых лесов (Кузнецко-Ала-таусском округе подтаежно-лесостепных лиственничных, горно-таежных темнохвойных лесов) и Восточносаянской горной лесорастительной провинции кедровых лесов (Бирюсинско-Китой-ском округе подтаежных сосновых и горно-таежных кедровых лесов).
Для проведения исследований заложено 169 пробных площадей, по которым получены таксационные данные, дано описание типов леса, лесорастительных условий, живого напочвенного покрова и почв (рис. 1). Результаты проведенных за 5 лет наземных исследований позволили сформировать сеть пробных площадей, включающую
Ка|43-слеиа НК1ЫКИ проЛоыи ЛЯОШНЯ » К(**«1й*(»;»ии ир«. Рка^баиго Ти пл. Ркп)блип Хикмя, Б^-ряшм
Рис. 1. Карта пробных площадей, заложенных для проведения исследований
значительное разнообразие насаждений по специфике лесорастительных условий, породному составу, лесообразовательным процессам. Сформированная сеть пробных площадей создана с учетом возможности ее использования для целей лесотак-сационного дешифрирования различных данных спутниковой съемки.
Пробные площади закладывали по общепринятым в лесоустройстве методам [10] - круговыми площадками постоянного радиуса. Для повышения точности данных радиус площадок в приспевающих, спелых и перестойных насаждениях составлял 20 м, а в средневозрастных - 16 м. В молодняках учёт осуществлялся площадками небольшого радиуса - 1,83 м по диагоналям или параллельным линиям в пределах пробной площади 50х50 м. На пробной площади выполнялась таксация древостоя, описание живого напочвенного покрова и закладка почвенного разреза. Для каждой пробной площади составлена карточка таксационно-де-шифровочного участка. При таксации деревьев
верхнего яруса, помимо учета всех таксационных показателей, дополнительно определяли их видимость и жизнеспособность. Возраст насаждений оценивался по годичным кольцам на кернах.
Предварительная обработка данных спутниковой съемки с аппарата 8епйпе1-1 выполнялась в 2 этапа. На первом этапе обработки радиолокационных изображений (РЛИ) проводилась их геопривязка в значении удельной эффективной площади рассеивания (УЭПР) с применением процедур уточнения орбитальных параметров, радиометрической калибровки, подавления спекл-шума, ортотрансформирования. Далее осуществлялось уточнение орбитальных параметров с целью достижения более высокой точности привязки исходного снимка. Для этого использовались файлы, поставляемые спустя некоторое время после съемки и содержащие уточнённую информацию о движении космического аппарата.
Следующим шагом предварительной обработки данных РЛИ была радиометрическая
калибровка - процесс перевода «сырых» значений пикселей (DN, digital number) к значениям УЭПР. В зависимости от проекции линейного разрешения РЛИ по дальности, существует три варианта нормировки эффективной площади рассеивания (ЭПР) геометрической площадью:
а) бета-ноль - ЭПР объекта, нормированная геометрической площадью, у которой линейный размер поперек линии пути находится в плоскости наклонной дальности;
б) сигма-ноль - ЭПР объекта, нормированная геометрической площадью, у которой линейный размер поперек линии пути находится в плоскости наземной дальности;
в) гамма-ноль - ЭПР объекта, нормированная геометрической площадью, у которой линейный размер поперек линии пути находится в плоскости, перпендикулярной к линии визирования.
С учетом расположения территории исследования ниже 60° с. ш. радиометрическая калибровка осуществлялась в 2 этапа: пересчет DN значения бета-ноль и нормирование до значений гамма-ноль с коррекцией уклонов по цифровой модели рельефа (ЦМР). В последующем осуществлялась координатная привязка радиометрически откалиброванного изображения (ортотрансфор-мирование). В процессе данной процедуры устранялись искажения, связанные с рельефом местности, с помощью ЦМР. Географические координаты записывались в систему WGS84; использовалась SRTM-модель рельефа с шагом пикселя 30 м.
На следующем шаге выполнено подавление спекл-шума - снижение влияния спекл-шума, возникающего на РЛИ вследствие когерентного наложения откликов от нескольких объектов в элементе разрешения изображения. В работе использовался метод некогерентного накопления, заключающийся в усреднении изображения в окне N*M пикселей с последующим укрупнением элементов изображения. Для процедуры выбран размер «окна», равный 2*2, из-за небольшого размера пробных площадей - до 1 га.
Вторым шагом в процедуре предварительной подготовки данных РЛИ было создание и применение маски затенений и переналожений вершин и склонов с целью исключения территорий, на
которых возможно получение ошибочных данных для анализа. В горной местности в зависимости от углов визирования для объектов разной высоты и наклона рельефа могут наблюдаться разные смещения склонов гор по наклонной дальности: растяжение, сжатие склонов, переналожение вершин и оснований склонов, радиолокационные тени. Первичная корректировка геометрических искажений РЛИ из-за эффектов сжатия и растяжения склонов осуществляется на этапе ортотранс-формирования путем подбора интерполяторов ЦМР высокого порядка - бикубического и на основе sinc-функции. Однако в зонах переналожений и затенений данная коррекция невозможна. Эта ситуация требует исключения из обработки и анализа этих зон путем маскирования. Маска затенений и переналожений строилась на основе данных геореференцированного РЛИ с Sentinel-1, снятого в интерферометрическом широкозахватном режиме (Interferometric Wide Swath, IW), уровня обработки GRDH (Ground Range Detected High Resolution - амплитудное РЛИ в проекции наклонной дальности c высоким уровнем детальности) в модуле SAR-Simulation Terrain Correction программного комплекса SNAP.
Подготовка данных радарной съемки по приведенному алгоритму позволила получить первичные данные для лесотаксационного дешифрирования. Ранее проведенные исследования показывают, что использование только данных спутниковой съемки, без связей с особенностями произрастания насаждений, не позволяет получить количественные и качественные характеристики лесов [8]. Наши исследования позволили разработать методы интеграции факторов, влияющих на продуктивность лесов, в методы лесотаксационного дешифрирования.
Обсуждение
Район проведения исследований характеризуется высоким разнообразием природных факторов, оказывающих влияние на породный состав и продуктивность лесов, процессы лесооб-разования. В частности, район включает зоны
тайги, лесостепи, степи и полупустыни, сухой степи. Разнообразие природно-климатических условий осложняется циркуляцией воздушных потоков, в том числе определяющих специфику формирующихся высотно-поясных комплексов. В районе исследований, согласно районированию горных лесов по В.Н. Смагину, выделены 4 лесорастительные области, 14 провинций, 29 округов.
С продвижением в глубь материка в районе исследований уменьшается годовая сумма осадков, происходит резкое колебание температуры воздуха и ее понижение. С продвижением на юг усиливается аридизация климата и снижается количество осадков. Согласно ранее проведенным исследованиям, оптимальные условия произрастания для основных лесообразующих пород наблюдаются на большей части Алтае-Са-янской горной области. Доминирующей формой рельефа является среднегорный эрозионный (800-2 000 м над ур. моря), большая часть которого занята лесами [9]. На высотах 1 500-2 000 м часто встречается рельеф древних поверхностей выравнивания, представленный широкими плато, покрытыми лесом или россыпями камней. Поверхности выравнивания заняты в основном подгольцово-таежными, подгольцово-субаль-пийскими, реже горно-таежными природными комплексами с преобладанием в составе пихты, кедра, иногда лиственницы и ели.
Особенности рельефа и высотная зональность оказывают значительное влияние на природно-климатические факторы и, как следствие, на растительность. Среди элементов рельефа основными факторами, влияющими на лесную растительность, являются высота местности над уровнем моря, экспозиция и крутизна склона [11, 12]. Разнообразие рельефа в горной местности приводит к значительному различию природных комплексов. При этом высота границы леса зависит от широты расположения местности, специфики рельефа.
В среднем с возрастанием абсолютной высоты на 100 м происходит снижение температуры на 0,5 °С, но на северных склонах падение температуры может быть ниже, чем на южных [13].
Каждые 100 м высоты в горах эквивалентны 1° широты. При этом в горных условиях, вследствие термической неоднородности вдоль склона и разности температур в приземном слое, возникают местные циркуляции воздуха, приводящие к изменению влажности. В ночной период в горном рельефе с вершины и склонов холодный воздух стекает вниз, приводя к резкому изменению температуры и влажности.
Все перечисленные процессы имеют существенную неоднородность, так как зависят от экспозиции и крутизны склона. Экспозиция склона определяет уровень радиации, распределение осадков, в том числе снега, и скорость ветра [14, 15]. Все эти факторы оказывают влияние на распределение лесообразующих пород в зависимости от экспозиции склона и высотной зональности. Приход солнечной радиации на склоны северной и южной экспозиции сильно различается. В пределах одного климатического пояса в условиях сложного горного рельефа склоны разной экспозиции имеют разный климат, что приводит к различиям и в лесооб-разовательном процессе. Наиболее высокие температуры отмечаются на юго-западных склонах и, как следствие, для этих склонов характерна пирогенная динамика формирования лесов. Северные склоны прогреваются хуже, на них в основном преобладают темнохвойные породы. Роль экспозиции склонов резко возрастает в районах с сухим климатом, где различия растительных сообществ на горных склонах выражены наиболее отчетливо [16, 17].
Крутизна склона влияет на возможность произрастания растений, образование почвы и ее прогрев. Каждый 1° южного склона в сторону увеличения перемещает условия произрастания на 10° географической широты к югу и, наоборот, каждый градус северного склона -к северу. Даже незначительный склон оказывает влияние на лесную растительность [18-20]. От крутизны склонов зависит степень эрозии и, соответственно, возможность почвообразовательного процесса. При высоких показателях крутизны склона, более 70°, скалы становятся практически недоступными для поселения
лесной растительности: она селится в выступах, промежутках между скалами, где есть возможность для формирования почвообразовательного процесса.
Проведённые под руководством профессора В.Н. Смагина работы выявили значительные различия в породном составе и растительности горных лесов юга Сибири в зависимости от широты, долготы, высотной зональности, экспозиции склона, особенностей климата [9].
Обобщение ранее проведенных исследований и анализ результатов экспериментальных работ показывают возможность создания математических моделей роста леса с учетом особенностей рельефа. Их эффективность зависит от специфики воздействия на произрастание лесной растительности региональных факторов, в частности, высоты над уровнем моря, экспозиции и крутизны склона. Перечисленные показатели можно получить из цифровых моделей рельефа, в том числе SRTM (NASA Shuttle Radar Topography Mission) [21, 22]. Оптимальным вариантом для целей лесотаксационного дешифрирования является использование моделей рельефа, полученных по данным лидарной аэрофотосъёмки, что исключает искажение данных в области тени при радиолокационной съемке.
Для решения задач лесотаксационного дешифрирования и разработки математических моделей необходим критерий, которым можно охарактеризовать комплекс изменяющихся условий. На основе работ В.Н. Смагина и проведенных нами исследований таким критерием, комплексно учитывающим воздействие природных факторов, определена верхняя граница леса. Для создания карты границ леса использовались данные со спутника Метеор-М, аппарата комплекса многозональной спутниковой съёмки (КМСС) с полосой захвата 900 км, разрешением 60 м. Анализ спутниковых данных позволил получить уравнение связи верхней границы леса с высотой над уровнем моря, экспозицией и крутизной склона, широтой и долготой. Математическая модель позволяет осуществлять прогноз высоты верхней границы леса для Южно-Сибирской горной зоны (рис. 2).
Зависимость можно выразить уравнением полиномиальной регрессии (1):
ВГЛ = -164240,689044 + 1867,06010105хХ + 3103,25029464хУ + 215,646402577хК + 2109,23354366хЭ - 34,9410542175xXxY -
2,46454303617хХхК + (1)
0,230633676384xYxК - 31,0209047004хХхЭ + 14,0276065139xYxЭ -79,8877612692 хКхЭ + 0xXxYxК + 0xXxYxЭ + 0,949230352913 хХхКхЭ -0,127414214933 xYxКxС + 0xXxYxКxС,
где:
ВГЛ - высота верхней границы леса, м; X - восточная долгота, градус; У - северная широта, градус; К - крутизна склона, градус; Э - экспозиция склона.
Уравнение характеризуется высокой степенью связи между наблюдаемыми и прогнозными значениями, коэффициент корреляции R=0,6. Повысить уровень связи можно при увеличении точности моделей рельефа и интеграции в модель климатических факторов.
Получена модель прогноза класса бонитета насаждений в зависимости от природных факторов. Результаты работы модели представлены
Рис. 2. Фрагмент карты верхней границы леса
на рис. 3. Модель можно описать уравнением факторной регрессии (2):
В = -87,9154224999 + 0,937624218007*Э + 5,89803386852хК + 0,0914915031603хВНМ + 0,0526150661517хВГЛ -0,41589335311хЭхК -0,00274025455906 хЭхВНМ -0,00555276512657хКхВНМ -0,000460769802699хЭхВГЛ - (2)
0,00359820294032хКхВГЛ -5,07543622073е - 005хВНМхВГЛ + 0,000428870647473хЭхКхВНМ + 0,000247069532734хЭхКхВГЛ + 1,25719843157е - 006хЭхВНМх ВГЛ3,2077470722е - 006хКхВНМхВГЛ -2,32337283732е - 007хЭхКхВНМхВГЛ,
где:
ВГЛ - высота верхней границы леса, м;
К - крутизна склона, градус;
Э - экспозиция склона;
ВНМ - высота над уровнем моря, м.
Оценка наблюдаемых и прогнозных значений показывает высокий уровень корреляции модели (Я=0,96, Р = 12, р=0). На рис. 4 показано,
что значения прогнозных показателей класса бонитета не пересекаются между собой, что позволяет провести зонирование исследуемой территории с высокой точностью. Повысить точность модели стало возможно благодаря исследованию по зонированию верхней границы леса и получению математической модели, позволяющей провести классификацию исследуемой территории по данному признаку.
Полученные данные демонстрируют, что показатель продуктивности горных лесов в значительной степени зависит от рельефа, в частности, от высотной зональности, экспозиции и крутизны склона. Высотной зональностью и экспозицией склона определяется не только продуктивность лесов, но и доминирование определённой породы в составе насаждения. Такая связь обусловлена воздействием различных природных процессов и климатическими особенностями.
Результаты
Рис. з. Дешифрирование класса бонитета насаждения на основании верхней границы леса (применение маски границы леса)
С учетом результатов ранее проведенных исследований [3, 4, 8, 23, 24] были определены показатели радиолокационной съемки, имеющие связь с таксационными характеристиками насаждений. Разработка моделей определения запаса, относительной полноты, среднего возраста насаждений осуществлялась с использованием следующих показателей радиолокационной спутниковой съемки: ^ Gamma_VH - УЭПР объекта в значении
гамма-ноль на поляризации VH. ^ Gamma_VV - УЭПР объекта в значении
гамма-ноль на поляризации ^ ^ ИУ1 - модификация радарного вегетационного индекса для двух поляризаций VV и УН. Представляет собой отношение четырех УЭПР объектов в значении гамма-ноль на поляризации УН к сумме УЭПР объектов в значении гамма-ноль на VV и VH поляризациях. ^ GLCMMean - суммарное среднее. Текстурный признак. Является произведением значения пикселя с частотой его появления с соседним пикселем.
^ GLCMCorre1ation - корреляция. Текстурный признак. Мера линейной связи уровней яркостей пар пикселей.
В уравнениях зависимостей применяются обозначения:
В - класс бонитета насаждений;
ВНМ - высота над уровнем моря, м;
ВГЛ - высота верхней границы леса, м;
Э - экспозиция склона;
К - крутизна склона, градус.
В моделях оценка экспозиции осуществлялась согласно классификации: 1 - Север, 2 - Северо-Восток, 3 - Восток, 4 - Юго-Восток, 5 - Юг, 6 - Юго-Запад, 7 - Запад, 8 - Северо-Запад.
Перечисленные показатели используются во всех формулах, представленных ниже.
Реализация расчетов предусматривает использование информации о зонировании территории района исследований по классу бонитета лесных насаждений и показателю верхней границы леса. Зонирование территории по классам бонитета лесных насаждений можно осуществить с использованием материалов лесоустройства, при их отсутствии или неактуальности - по данным цифровых моделей рельефа.
Апробация разработанных методов лесотакса-ционного дешифрирования проведена на территории Саяно-Шушенского заповедника с определением показателей запаса и полноты насаждений.
Результаты исследований показывают возможность определения запаса насаждений с высокой достоверностью (3). Коэффициент корреляции между наблюдаемыми и прогнозными значениями составляет Я = 0,74 № = 18, р = 0) (рис. 5).
Связь запаса насаждений (М) с показателями радиолокационной съемки характеризуется полиномиальной зависимостью, которую можно представить формулой (3):
М = 583,639211956 -3938,63692387 хGamma_VH + 58308,9862794хGamma_VH^2 + 644,947333587хGamma_VV -1044,7974875хGamma_W^2 + 216,809476031 хЯУ1 -79,0660294597хЯУ1^2 + 11,0477411441 хGLCMMean -0,0704335968988хGLCMMean^2 - 195,947 (3) 842693 хGLCMCorre1ation -61,1594919669хGLCMCorre1ation^2 -198,2247664хВ + 16,4666729983хВ^2 -0,0325109965151хВНМ + 2,59467437601е-005хВНМ^2 -
0,0461101777147хВГЛ + 5,20766016225е-006хВГЛ^2 -34,0988187075 хЭ + 3,93852859431 хЭ^2 + 2,64464290178хК - 0,0949154300019хК^2.
Рис. 4. Оценка эффективности модели по прогнозу класса бонитета насаждения
Рис. 5. Оценка эффективности модели по прогнозу запаса насаждения
URL http://lhi.vniilm.ru/ 75
Результаты работы модели по формуле представлены на рис. 6.
Полнота насаждения (4) является наиболее трудно определяемым показателем. Уровень
Рис. 6. Дешифрирование запаса насаждений по данным радиолокационной съемки
Рис. 7. Дешифрирование полноты насаждений по данным радиолокационной съемки
корреляции между прогнозными и экспериментальными значениями не превышает 0,5 (Р = 5, р = 0). Это может быть связано с тем, что большинство исследований проведено в сложных двухъярусных насаждениях. Отражение сигнала Б-диапазона от верхнего полога насаждений, слабое его проникновение под полог ограничивает определение показателя относительной полноты насаждений (Р):
Р = 1,43606484464 -8,93928619974хGamma_VH + 114,169968515 хGamma_VH^2 +
2,56959386264хGamma_VV -6,17665299064хGamma_W^2 + 0,117400468196хЯУ1 + 0,144535500605хЯУ1^2 + 0,0160772011437хGLCMMean -0,000117121519459хGLCMMean^2 -1,51816304489 хGLCMCorre1ation + (4) 0,622054320851 хGLCMCorre1ation^2 -0,192946906605хВ + 0,0316702329601 хВ^2 -0,00029731827993хВНМ + 5,22537524734е - 008хВНМ^2 -
0,000311058319558хВГЛ + 1,55849566358е - 007хВГЛ^2 -0,0168759474307хЭ + 0,00148781112477хЭ^2 + 0,0115133245654хК0,000419238141324хК^2.
Результаты работы модели по формуле представлены на рис. 7.
Выводы
Результаты представленного исследования демонстрируют возможность использования радиолокационных данных в решении задач лесотаксационного дешифрирования по определению количественных и качественных характеристик насаждений для эксплуатационных и резервных горных лесов. Разработанные предложения позволяют осуществить весь спектр подготовительных работ по данным
спутниковой съемки и последующей ее обработки с целью получения таксационных показателей лесов. Использование комплексного подхода с возможностью бонитировки территории на основе моделей рельефа или данных лесоустройства дает возможность увеличить точность моделей лесотаксационного дешифрирования до уровня, необходимого при выполнении таксации по третьему разряду лесоустройства.
Разработанные модели позволяют существенно сократить объем полевых таксационных работ и ускорить процесс лесоустройства труднодоступных горных лесов. На основе показанного принципа проведения лесотаксационных работ возможно создание автоматизированных сервисов с учетом оптимизации различных алгоритмов машинного обучения для решения задач лесоустройства горных лесов.
ГЕОИНФОРМПЦИОННЫЕ ТЕХНОЛОГИИ
Список источников
1. Об утверждении Стратегии развития лесного комплекса Российской Федерации до 2030 года. Распоряжение Правительства Российской Федерации от 11.01.2021 № 312-р 2021.
2. Об утверждении государственной программы Российской Федерации «Развитие лесного хозяйства». Постановление Правительства Российской Федерации от 15.04.2014 № 318. - Текст электронный - Режим доступа: URL: https://www.garant.ru/products/ipo/prime/doc/70544228/ (дата обращения: 29.05.2019).
3. Возможности таксации лесов на основе данных радиолокационной спутниковой съемки / В.М. Сидоренков, О.В. Кушнырь, Л.В. Бадак, Д.О. Астапов, А.В. Букась // Исследование Земли из космоса. - 2021. - № 5. - C. 72-84.
4. Возможности таксации лесов на основе данных радиолокационной спутниковой съемки на примере ленточных боров Алтайского края / В.М. Сидоренков, Д.О. Астапов, Л.А. Бадак, Ю.С. Ачиколова // Современные проблемы дистанционного зондирования Земли из космоса : материалы 19-й Международной конференции. -Москва : Институт космических исследований РАН, 2021. - C. 379.
5. Radar satellite survey data-based forest inventory opportunities / V.M. Sidorenkov, O.V. Kusnir, L.V. Badak, D.O. Astapov, A.V. Bukas // Izvestiya, Atmospheric and Oceanic Physics. - 2021. - Vol. 57. - № 12. - P. 1549-1560.
6. Potential of Sentinel-1 Data for Monitoring Temperate Mixed Forest Phenology / P.-L. Frison, В. Fruneau, S. Kmiha, К. Soudani, Е. Dufrene, Т. Le Toan, Т. Koleck, L.Villard, Е. Mougin, J.-P. Rudant // Remote Sensing. - 2018. - Т. 10. - C. 2049.
7. Sensitivity analysis of multi-temporal Sentinel-1 SAR parameters to crop height and canopy coverage / R. Nas, F. Balik Sanli, S. Abdikan, Z. Qakir, А. Sekertekin, М. UstUner // Applied Sciences. - 2019. - Т. 9. - C. 655.
8. Возможности определения таксационных показателей смешанных насаждений по данным радиолокаци-оной съемки со спутника Sentinel-1 / В.М. Сидоренков, В.Н. Косицын, Л.А. Бадак, Д.О. Астапов, Ю.С. Ачиколова // Исследование Земли из космоса. - 2023. - № 1. - C. 53-65.
9. Смагин, В.Н. Типы лесов гор Южной Сибири / В.Н. Смагин. - Новосибирск : Наука. Сибирское отделение, 1980. - 334 с.
10. Об утверждении Лесоустроительной инструкции. Приказ Министерства природных ресурсов и экологии Российской Федерации от 29.03.2018 № 122.
11. Соколова, Г.Г. Влияние высоты местности, экспозиции и крутизны склона на особенности пространственного распределения растений / Г.Г. Соколова // Acta Biologica Sibirica. - 2016. - Т. 2. - № 3. - C. 34-45.
12. Фридланд, В.М. Структура почвенного покрова / В.М. Фриндланд. - Москва : Наука, 1983. - 196 с.
13. Берг, Л.С. Географические зоны Советского Союза / Л.С. Берг. - Москва : Географгиз, 1938. - 297 с.
14. Isard, S.A. Factor influencing soil moisture and plant community distribution on Niwot Ridge, Front Range, Colorado, USA / S.A. Isard // Arctic and Alpine Research. - 1986. - № 18. - C. 83-96.
15. Zhang, Y.Z. Physical geography in Tibet / Y.Z. Zhang, D. Zheng, Q.Y. Yang. - Beijing : Science press, 2002. - 178 с.
16. Захарова, А.Ф. Радиационный режим северных и южных склонов в зависимости от географической широты / А.Ф. Захарова. - Ученые записки ЛГУ. Сер. географические науки. - № 269, 1959. - 269 с.
17. Станюкович, К.В. Растительность высокогорий СССР / К.В. Станюкович. - Душанбе : Дониш, 1960. - 523 с.
18. Мандельштам, Н.Л. Инсоляция склонов различной экспозиции и крутизны в условиях Саратовской области / Н.Л. Мандельштам // Вопросы климата и погоды Нижнего Поволжья. - 1966. - C. 122-130.
19. Щербаков, Ю.А. Из опыта изучения роли экспозиции в ландшафтоведении. Влияние экспозиции на ландшафты / Ю.А. Щербаков. - Москва : изд-во АН СССР, 1970. - 399 с.
20. Антюфеев, В.В. Динамика продуктивности лугово-степных сообществ крымской яйлы в связи с изменчивостью солнечной радиации / Закономерности развития органического мира / В.В. Антюфеев. - Минск : Наука и техника, 1978. - С. 5.
21. The Shuttle Radar Topography Mission / T.G. Farr, Р.А. Rosen, Е. Caro, R. Crippen, R. Duren, S. Hensley [et. al.] // Reviews of Geophysics. - 2007. - Vol. 45. - № 2. - RG2004. doi:10.1029/2005RG000183.
22. An assessment of the SRTM topographic products, Technical Report JPL D-31639 / Е. Rodriguez, C.S. Morris, J.E. Belz, Е. Chapin, J. Martin, W. Daffer, S. Hensley. - JPL, NASA, 2005. - 143 с.
23. Современные методы таксации лесов на основе данных радарной спутниковой съемки / В.М. Сидоренков, Л.А. Бадак, О.В. Кушнырь, А.С. Рыбкин, Д.О. Астапов // Современные проблемы дистанционного зондирования Земли из космоса : материалы 18-й Всероссийской открытой конференции. - Москва : Институт космических исследований РАН, 2020. - C. 441.
24. Методические рекомендации по подготовительным работам к лесотаксационному дешифрированию данных радиолокационной съемки / В.М. Сидоренков, В.Н. Косицын, О.В. Кушнырь, Д.О. Астапов, А.С. Рыбкин, Л.А. Бадак, Л.А. Костюк. - Пушкино : ВНИИЛМ, 2020. - 88 с.
References
1. Ob utverzhdenii Strategii razvitiya lesnogo kompleksa Rossijskoj Federacii do 2030 goda. Rasporyazhenie Pravitel'stva Rossijskoj Federacii ot 11.01.2021 № 312-r 2021.
2. Ob utverzhdenii gosudarstvennoj programmy Rossijskoj Federacii «Razvitie lesnogo hozyajstva». Postanovlenie Pravitel'stva Rossijskoj Federacii ot 15.04.2014 № 318. - Tekst elektronnyj - Rezhim dostupa: URL: https://www. garant.ru/products/ipo/prime/doc/70544228/ (data obrashcheniya: 29.05.2019).
3. Vozmozhnosti taksacii lesov na osnove dannyh radiolokacionnoj sputnikovoj s"emki / V.M. Sidorenkov, O.V. Kushnyr', L.V. Badak, D.O. Astapov, A.V. Bukas' // Issledovanie Zemli iz kosmosa. - 2021. - № 5. - C. 72-84.
4. Vozmozhnosti taksacii lesov na osnove dannyh radiolokacionnoj sputnikovoj s"emki na primere lentochnyh borov Altajskogo kraya / V.M. Sidorenkov, D.O. Astapov, L.A. Badak, Yu.S. Achikolova // Sovremennye problemy distancionnogo zondirovaniya Zemli iz kosmosa : materialy 19-j Mezhdunarodnoj konferencii. - Moskva : Institut kosmicheskih issledovanij RAN, 2021. - C. 379.
5. Radar satellite survey data-based forest inventory opportunities / V.M. Sidorenkov, O.V. Kusnir, L.V. Badak, D.O. Astapov, A.V. Bukas // Izvestiya, Atmospheric and Oceanic Physics. - 2021. - Vol. 57. - № 12. - P. 1549-1560.
6. Potential of Sentinel-1 Data for Monitoring Temperate Mixed Forest Phenology / P.-L. Frison, V. Fruneau, S. Kmiha, K. Soudani, E. Dufrene, T. Le Toan, T. Koleck, L.Villard, E. Mougin, J.-P. Rudant // Remote Sensing. - 2018. - T. 10. - C. 2049.
7. Sensitivity analysis of multi-temporal Sentinel-1 SAR parameters to crop height and canopy coverage / R. Nas, F. Balik Sanli, S. Abdikan, Z. Qakir, A. Sekertekin, M. Üstüner // Applied Sciences. - 2019. - T. 9. - C. 655.
8. Vozmozhnosti opredeleniya taksacionnyh pokazatelej smeshannyh nasazhdenij po dannym radiolokacionoj s"emki so sputnika Sentinel-1 / V.M. Sidorenkov, V.N. Kosicyn, L.A. Badak, D.O. Astapov, Yu.S. Achikolova // Issledovanie Zemli iz kosmosa. - 2023. - № 1. - C. 53-65.
9. Smagin, V.N. Tipy lesov gor YUzhnoj Sibiri / V.N. Smagin. - Novosibirsk : Nauka. Sibirskoe otdelenie, 1980. - 334 s.
10. Ob utverzhdenii Lesoustroitel'noj instrukcii. Prikaz Ministerstva prirodnyh resursov i ekologii Rossijskoj Federacii ot 29.03.2018 № 122.
11. Sokolova, G.G. Vliyanie vysoty mestnosti, ekspozicii i krutizny sklona na osobennosti prostranstvennogo raspredeleniya rastenij / G.G. Sokolova // Acta Biologica Sibirica. - 2016. - T. 2. - № 3. - C. 34-45.
12. Fridland, V.M. Struktura pochvennogo pokrova / V.M. Frindland. - Moskva : Nauka, 1983. - 196 s.
13. Berg, L.S. Geograficheskie zony Sovetskogo Soyuza / L.S. Berg. - Moskva : Geografgiz, 1938. - 297 s.
14. Isard, S.A. Factor influencing soil moisture and plant community distribution on Niwot Ridge, Front Range, Colorado, USA / S.A. Isard // Arctic and Alpine Research. - 1986. - № 18. - C. 83-96.
15. Zhang, Y.Z. Physical geography in Tibet / Y.Z. Zhang, D. Zheng, Q.Y. Yang. - Beijing : Science press, 2002. - 178 s.
16. Zaharova, A.F. Radiacionnyj rezhim severnyh i yuzhnyh sklonov v zavisimosti ot geograficheskoj shiroty / A.F. Zaharova. - Uchenye zapiski LGU. Ser. geograficheskie nauki. - № 269, 1959. - 269 s.
ГЕОИНФОРМПЦИОННЫЕ ТЕХНОЛОГИИ
17. Stanyukovich, K.V. Rastitel'nost' vysokogorij SSSR / K.V. Stanyukovich. - Dushanbe : Donish, 1960. - 523 s.
18. Mandel'shtam, N.L. Insolyaciya sklonov razlichnoj ekspozicii i krutizny v usloviyah Saratovskoj oblasti / N.L. Mandel'shtam // Voprosy klimata i pogody Nizhnego Povolzh'ya. - 1966. - C. 122-130.
19. Shcherbakov, Yu.A. Iz opyta izucheniya roli ekspozicii v landshaftovedenii. Vliyanie ekspozicii na landshafty / Yu.A. Shcherbakov. - Moskva : izd-vo AN SSSR, 1970. - 399 s.
20. Antyufeev, V.V. Dinamika produktivnosti lugovo-stepnyh soobshchestv krymskoj yajly v svyazi s izmenchivost'yu solnechnoj radiacii / Zakonomernosti razvitiya organicheskogo mira / V.V. Antyufeev. - Minsk : Nauka i tekhnika, 1978. - S. 5.
21. The Shuttle Radar Topography Mission / T.G. Farr, R.A. Rosen, E. Caro, R. Crippen, R. Duren, S. Hensley [et. al.] // Reviews of Geophysics. - 2007. - Vol. 45. - № 2. - RG2004. doi:10.1029/2005RG000183.
22. An assessment of the SRTM topographic products, Technical Report JPL D-31639 / E. Rodriguez, C.S. Morris, J.E. Belz, E. Chapin, J. Martin, W. Daffer, S. Hensley. - JPL, NASA, 2005. - 143 s.
23. Sovremennye metody taksacii lesov na osnove dannyh radarnoj sputnikovoj s"emki / V.M. Sidorenkov, L.A. Badak, O.V. Kushnyr', A.S. Rybkin, D.O. Astapov // Sovremennye problemy distancionnogo zondirovaniya Zemli iz kosmosa : materialy 18-j Vserossijskoj otkrytoj konferencii. - Moskva : Institut kosmicheskih issledovanij RAN, 2020. - C. 441.
24. Metodicheskie rekomendacii po podgotovitel'nym rabotam k lesotaksacionnomu deshifrirovaniyu dannyh radiolokacionnoj s"emki / V.M. Sidorenkov, V.N. Kosicyn, O.V. Kushnyr', D.O. Astapov, A.S. Rybkin, L.A. Badak, L.A. Kostyuk. - Pushkino : VNIILM, 2020. - 88 s.
