Оценка высоты и плотности лесного полога
Высота и плотность лесного полога используются как переменные при различных исследованиях окружающей среды. К таким исследованиям можно отнести оценку биомассы, экстента и состояния леса, а также изучение биоразнообразия. Плотность лесного полога - это отношение растительного полога к открытой поверхности земли при просмотре с воздуха. Высота лесного полога означает, на какой высоте находятся кроны деревьев. Лидар может быть использован для определения обоих этих переменных.
Ниже приведены шаги, которые необходимы для вычисления плотности и высоты покрова по данным лидара. В первую очередь требуются данные лазерного сканирования с разделением на точки, находящиеся на поверхности земли, и точки, находящиеся над этой поверхностью. Такая классификация точек обычно выполняется поставщиком данных. Во-вторых, необходимо узнать, когда были получены эти данные, а также установить тип растительности в изучаемом районе. Если в изучаемом районе много лиственных деревьев, а данные были собраны осенью (во время листопада), то плотность вычислить не удастся.
Создание набора данных и слоя LAS
Набор данных LAS может использоваться для вычисления плотности и высоты полога леса. Примените инструмент геообработки Создать набор данных LAS (Create LAS Dataset) для создания ссылок на данные лазерного сканирования в формате LAS. Если вы имеете дело с данными в текстовом формате ASCII или, если вы желаете работать с точками формата LAS в базе геоданных, обратите внимание на альтернативный ход работы в конце данного раздела. После создания набора данных LAS определите слой, указав свойства фильтра, который будет определять, какие точки из файлов LAS задействованы в работе. Для создания слоя добавьте набор данных к карте, либо запустите инструмент геообработки Создать слой набора данных LAS (Make LAS Dataset Layer). Задавая слой, следует выбрать все возвращаемые значения, а затем обратить внимание на коды классификации. В основном мы будем работать с двумя наборами точек: с точками на земной поверхности и с точками над ней. В следующей таблице содержатся коды классификации формата LAS, определенные стандартом LAS 1.1:
Код классификации | Тип классификации |
|---|---|
0 | Классификация не выполнялась |
1 | Не назначено |
2 | Земля |
3 | Низкорослая растительность |
4 | Растительность средней высоты |
5 | Высокая растительность |
6 | Здание |
7 | Шум |
8 | Ключ модели |
9 | Вода |
10 | Зарезервировано ASPRS |
11 | Зарезервировано ASPRS |
12 | Наложение |
13–31 | Зарезервировано ASPRS |
Примечание:Любой файл LAS, созданный в последние несколько лет, должен использовать эти коды, если его точки были классифицированы. К сожалению, до сих пор в классификации имеется некоторая неопределенность. Например, класс 2 означает точки на земле, но класс 8 - также означает точки земли. Точки класса 8, или ключ модели, являются специальным набором точек, расположенных на земле, использующимся для построения изолиний или для других задач, где требуется уменьшенный набор точек земли. Наличие этих точек зависит от способа обработки данных. Если вы не уверены, выберите оба класса. Если окажется, что точек модели нет, это не помешает. Обозначения растительности имеют схожие проблемы. Иногда поставщики данных относят всё, что находится над землей, к классу 1, поскольку они не проводили детальной классификации. Поэтому, если вы не знаете особенностей классификации ваших данных, выберите все точки над земной поверхностью, т. е. из классов 1, 3, 4 и 5. Такое универсальное решение позволит захватить все точки растительности. Обратите внимание, если в классе 1 находятся здания или другие искусственные объекты, они также будут захвачены и могут исказить, в некоторой степени, конечные результаты.
Вычисление плотности
Наиболее эффективный способ вычислить плотность покрова - разделить изучаемый район на несколько равных участков меньшего размера, используя растеризацию. В каждой ячейке растра необходимо сравнить число точек, находящихся над поверхностью, и общее число точек.
Очень важно определить наиболее подходящий размер ячейки для данного анализа. Размер должен быть как минимум в четыре раза больше, чем средний интервал между точками. Размер ячейки можно сделать больше, но не меньше.
- Примените к точкам, лежащим над поверхностью, инструмент геообработки Статистика точек LAS в виде растра (LAS Point Statistics As Raster) с опцией Число возвращенных точек (RETURN_COUNT).
- Преобразуйте все полученные ячейки со значением NoData в 0, так чтобы в последующих операциях ячейки, не содержащие точек, рассматривались как 0. Для этого примените последовательно инструменты геообработки Равно нулю (IsNull) и Условие (Con).
- Повторите шаги 1 и 2 для точек с поверхности земли.
- Сложите вместе растры надземных и наземных точек, используя инструмент геообработки Сложить (Plus), чтобы определить общее число точек на ячейку.
- Все полученные растры содержат тип данных Long (Длинное целое). Необходимо, чтобы один растр имел данные с плавающей точкой, чтобы получить выходные данные такого же типа в результате работы инструмента геообработки Разделить (Divide), который будет использован на этапе 6. Чтобы получить растр с плавающей точкой, используйте результат работы инструмента Сложить (Plus) в качестве входных данных для инструмента Плавающая точка (Float).
- Теперь используйте инструмент геообработки Разделить (Divide), чтобы выполнить сравнение растра, соответствующего числу точек над поверхностью, с растром общего количества точек в формате с плавающей точкой. Вы получите диапазон от 0.0 до 1.0, где 0.0 означает отсутствие покрова, а 1.0 - очень плотный покров.
Показанное изображение представляет плотность полога. Самые светлые области имеют редкий покров или не имеют его совсем. Это те области, в которых высокий процент лазерных импульсов достигает земной поверхности. Области темно-зеленого цвета, где лидар не может проникнуть до земли, обозначают очень плотный покров.
Вычисление высоты
Чтобы определить высоту растительного покрова, необходимо вычесть поверхность земли (DEM) из поверхности верхних точек (DSM). Инструкция для построения этих поверхностей приведена в разделе Создание растровых ЦМР (DEM) и ЦММ (DSM) на основе обширных наборов точек лазерного сканирования.
- После того как вы получили два растра (растр первых возвращаемых значений и растр земной поверхности), используйте инструмент геообработки Вычесть (Minus) для определения разности между этими наборами данных. Полученная разность равна высоте лесного полога.
На рисунке ниже показана высота полога над поверхностью земли. Его цвет варьирует от голубого (минимальная или нулевая) до оранжевого (максимальная высота).
Лидар может использоваться для вычисления плотности и высоты растительного покрова. Это необходимо для многих задач, включая оценку объема биомассы или управление лесными ресурсами.
Альтернативный ход работы
Если точки лазерного сканирования содержатся в файлах формата ASCII, необходимо иметь два набора файлов: одни файлы должны содержать только точки земной поверхности, а другие — только точки, расположенные выше этой поверхности. Имея подобную структуру файлов, следует загрузить точки в класс пространственных объектов-мультиточек в базе геоданных при помощи инструмента геообработки ASCII 3D в класс объектов (ASCII 3D To Feature Class). Запустите инструмент дважды, отдельно для каждого набора файлов. Затем запустите инструмент геообработки Точки в растр (Point To Raster), который здесь выполняет функцию инструмента Статистика точек LAS (LAS Point Statistics), применяемого в описанном выше ходе работы с набором данных LAS.
Если данные имеют формат LAS, но вы желаете загрузить в базу геоданных только некоторые из точек и работать только с ними, примените инструмент LAS в мультиточки (LAS To Multipoint) два раза. Один раз — для загрузки точек земной поверхности, второй раз — для загрузки точек, лежащих над ней. Затем запустите инструмент геообработки Точки в растр (Point To Raster), который здесь выполняет функцию инструмента Статистика точек LAS (LAS Point Statistics), применяемого в описанном выше ходе работы с набором данных LAS.