Цель chmloader — загрузить данные модели высоты купола (CHM) из недавней работы Толана и др. (2024). Краткое изложение этой работы можно найти здесь. Данные загружаются из хранилища AWS s3. Более подробную информацию о сегменте можно найти здесь.
Вы можете установить chmloader следующим образом:
# install.packages("pak")
pak :: pkg_install( " TESS-Laboratory/chmloader " ) Это базовый пример, показывающий, как загрузить некоторые данные. Функция download_chm использует gdalwarp (через sf::gdal_utils ) для эффективного получения только необходимых данных из нескольких тайлов — разрешение по умолчанию составляет 1 м, но его можно перепроецировать по мере необходимости, используя аргумент res .
library( chmloader )
parana_cuiana <- sf :: st_point(c( - 61.89 , - 4.12 )) | >
sf :: st_sfc( crs = 4326 ) | >
sf :: st_buffer( 3000 )
pc_chm <- download_chm(
parana_cuiana ,
filename = tempfile( fileext = " .tif " )
)
terra :: plot( pc_chm , col = hcl.colors( 256 , " viridis " ))Этот пакет также предоставляет простую функцию для создания графиков для сравнения различных CHM. Цель этой функции — обеспечить простую и надежную оценку Tolan et al. (2024) Данные CHM с моделями на основе LiDAR и другими продуктами ML. Пакет chmloader поставляется с небольшим набором примеров наборов данных CHM на основе LiDAR, полученных из набора данных объектной модели растительности Английского агентства по охране окружающей среды. Ниже приведен пример использования одного из этих примеров наборов данных из Фингл-Вудс, Девон, Великобритания:
fingle_woods <- reference_data( " fingle_woods " )
compare_models( fingle_woods , aggregate = 10 , drop_zeros = TRUE )
# > ℹ meta/WRI CHM not provided, downloading now...
# > ✔ CHM downloaded successfully! Обратите внимание, что в этом примере aggregate аргумент используется для уменьшения разрешения как эталонного, так и Meta/WRI CHM в 10 раз (что приводит к модели 10 м) и проверки обеих моделей более крупного масштаба в дополнение к исходной модели 1 м. модель. Эта функция может помочь выяснить, каково истинное разрешение Meta/WRI CHM и как оно соотносится с моделью на основе LiDAR в разных масштабах.
Кроме того, аргумент drop_zeros используется для удаления нулевых значений как из графика 2d плотности, так и из производной статистики, где значения из справочных/эталонных данных и Meta/WRI CHM оба равны нулю. Это особенно полезно, когда основной целью является оценка кроны деревьев, а не отсутствие деревьев и/или когда древесный покров редок; однако значение drop_zeros по умолчанию равно FALSE .
