提取矢量点位对应5×5窗口栅格像元值

文章预览

有学员自己写了一部分代码，运行不出结果，让我改改，下面是学员代码和修改过程，供大家学习参考。 任务要求 我有一个excle表格，里边包含一组带有编号、X、Y地理坐标（WGS84）的坐标点，另外给定一个含有空间参考坐标系和空间分辨率的栅格数据，想编写一个R语言代码，提取以表格中的坐标点为中心所对应的栅格数据5x5窗口中的像元值，然后存放到表格中相应的坐标点后。 学员代码 # 加载所需的包 library(readxl) library(sf) library(terra) library(writexl) # 读取Excel表格中的坐标点 coords < - read_excel( "I:/2024年课题/地面样地/KT-1.xlsx" )   # 替换为实际文件路径 # 将坐标点转换为sf对象（WGS84坐标系） pts_sf < - st_as_sf(coords, coords = c( "X" ,  "Y" ), crs = 4326) # 读取栅格数据 raster < - rast( "I:/2024年课题/无人机光谱指数/SAVI_KT1.tif" )   # 替换为实际文件路径 # 将点转换到 ………………………………