关于导入特定区域卫星影像及创建危险性栅格文件的技术问询

一、导入特定区域的卫星影像 #

在R中处理卫星影像，常用raster包（或更现代的terra包）来完成。核心步骤是定义感兴趣区域（ROI），再导入影像并裁剪到目标范围：

# 加载所需包
library(raster)

# 1. 定义特定区域的边界（示例为UTM坐标，可根据你的数据调整）
roi_extent <- extent(696583.6, 696799.6, 167579.6, 167789.6)

# 2. 导入卫星影像（替换为你的影像文件路径）
# 单波段影像用raster()，多波段用stack()
sat_image <- raster("your_satellite_image.tif")

# 3. 裁剪影像到特定区域
sat_image_roi <- crop(sat_image, roi_extent)

# 可选：查看裁剪后的影像
plot(sat_image_roi)

如果使用更高效的terra包，代码类似：

library(terra)
roi_extent <- ext(696583.6, 696799.6, 167579.6, 167789.6)
sat_image <- rast("your_satellite_image.tif")
sat_image_roi <- crop(sat_image, roi_extent)

二、基于模拟数据创建危险性栅格（危险性地图） #

假设你有模拟的落石轨迹数据（包含位置和动能等危险性指标），可以通过以下步骤创建危险性栅格：

步骤说明 #

• 创建匹配目标区域的空栅格模板（定义范围、分辨率、坐标系）；
• 将模拟轨迹数据栅格化，计算每个栅格单元的统计值（如平均动能）作为危险性评分。

完整代码（补全你提供的片段） #

library(raster)

# 1. 创建空栅格模板，匹配目标区域参数
rockfall_intensity <- raster(
  xmn = 696583.6, xmx = 696799.6, 
  ymn = 167579.6, ymx = 167789.6, 
  res = 2,  # 栅格分辨率（单位：米）
  crs = "+proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +k_0=1 +x_0=2600000 +y_0=1200000 +ellps=bessel +towgs84=674.374,15.056,405.346,0,0,0,0 +units=m +no_defs"
)

# 假设trjct是你的模拟轨迹数据框，包含x、y坐标和kinetic_energy（动能）列
# 2. 栅格化轨迹数据，计算每个栅格的平均动能作为危险性强度
rockfall_intensity <- rasterize(
  x = trjct[, c('x', 'y')],  # 轨迹点的坐标
  y = rockfall_intensity,    # 空栅格模板
  field = trjct$kinetic_energy,  # 用于计算危险性的指标（动能）
  fun = mean,  # 用平均动能代表该栅格的危险性，可替换为sum、max等
  na.rm = TRUE
)

# 查看最终的危险性地图
plot(rockfall_intensity, main = "落石危险性强度地图")

注意事项 #

• 确保trjct数据框包含有效的x/y坐标和对应的危险性指标（如动能）；
• fun参数可根据需求调整：比如用max表示栅格内的最大动能，sum表示总动能，以此匹配不同的危险性评估逻辑。

内容的提问来源于stack exchange，提问作者user1607