为什么raster::extract设置buffer参数时提取栅格值会返回多个值？

为raster::extract()的点输入添加buffer参数后，函数的底层执行逻辑会发生本质切换：

• 无buffer的点提取场景：输入几何为0维点，每个点仅匹配其几何位置落在的单个栅格单元，默认返回等长数值向量，每个点对应1个栅格值，与你观察到的行为一致。
• 带buffer的提取场景：函数会首先执行几何运算，把每个临时点转换为以点为圆心、buffer值为半径的圆形矢量面，后续提取逻辑完全等同于“输入一组圆形多边形提取栅格值”——只要和该圆形面存在空间重叠的所有栅格像元，都会被纳入提取范围。

如果不指定额外的汇总规则，函数默认会把每个缓冲区覆盖的所有栅格值全部返回，因此每个点对应的返回结果是一个值集合（长度等于缓冲区内覆盖的栅格像元总数），自然就出现了“一个点对应多个栅格值”的现象，这是函数的固有设计逻辑，不属于程序报错或异常行为。

需要纠正一个普遍的认知偏差：buffer参数并不是“以点为中心的近距搜索半径”，它的本质是临时生成矢量面用于后续面域提取，不存在“半径内只取1个值”的默认逻辑。

根据实际分析需求，可选择对应的写法让每个点仅返回1个汇总结果：

• 需求为提取250米缓冲区内的栅格值统计量（均值、最大值、中位数等）
  直接通过fun参数指定汇总函数，函数会自动将每个缓冲区覆盖的所有栅格值传入该函数计算，每个点仅返回1个统计值。例如计算缓冲区内的栅格平均值：
  raster::extract(my_raster, my_points, buffer=250, fun=mean, na.rm=TRUE)
• 需求为按像元到点中心的距离做权重计算，得到更贴合点位置的代表值
  开启weights=TRUE参数，函数会自动计算每个重叠像元到缓冲区中心的距离权重，搭配自定义汇总函数即可得到加权结果。例如计算距离加权的栅格平均值：
  raster::extract(my_raster, my_points, buffer=250, weights=TRUE, fun=function(x, w) weighted.mean(x, w, na.rm=TRUE))
• 需求为保留缓冲区覆盖的所有像元明细，但整理成规整的表格结构
  传入df=TRUE参数，函数会返回带点ID的规整数据框，每个点的多个像元值会对应多行记录，不会出现零散的列表结构。

注意：如果原本需求是“仅提取点本身位置的栅格值”，完全不需要添加buffer参数；如果是为了规避点落在栅格边界导致的空值提取，更优的方案是设置method='bilinear'，直接计算点所在位置的双线性插值栅格值，既不会产生多值，也能避免边界空值问题。

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