if (FALSE) {
# These examples require external data files not included with the package
# Extract terrain variables for study sites
sites_with_terrain <- integrate_terrain_analysis(
vector_data = "study_sites.shp",
elevation_raster = "dem.tif",
terrain_vars = c("slope", "aspect", "TRI", "TPI")
)
# Use custom terrain functions
custom_functions <- list(
ruggedness = function(sf_data) {
sf_data$slope * sf_data$TRI
}
)
terrain_analysis <- integrate_terrain_analysis(
vector_data = field_boundaries,
elevation_raster = dem_raster,
custom_terrain_functions = custom_functions
)
}
Run the code above in your browser using DataLab