raster包 设置高程栅格数据，计算坡度坡向

# 安装并加载raster包
if (!require(raster)) install.packages("raster")
library(raster)

# 设置高程栅格数据的路径
elevation_path <- "F:/D/qing/wc2.1_2.5m_elev.tif"

# 读取高程栅格数据
elevation_raster <- raster(elevation_path)

# 使用terrain函数计算坡度
slope_raster <- terrain(elevation_raster, opt = "slope", unit = "degrees")

# 使用terrain函数计算坡向
# 坡向的默认单位是度（degrees），也可以使用'rad'获取弧度值
aspect_raster <- terrain(elevation_raster, opt = "aspect")

# （可选）保存坡度结果
writeRaster(slope_raster, "F:/D/qing/slope.tif", format = "GTiff")

# （可选）保存坡向结果
writeRaster(aspect_raster, "F:/D/qing/aspect.tif", format = "GTiff")

# 如果需要查看结果，可以使用plot函数
plot(slope_raster, main = "Slope in degrees")
plot(aspect_raster, main = "Aspect")