如何使用带有R stars软件包的多边形从光栅中提取值?

如何使用带有R stars软件包的多边形从光栅中提取值?,r,r-stars,R,R Stars,使用stars软件包,可以使用st_extract()函数从定义位置的光栅中提取值 库(星型) #>正在加载所需的包:abind #>装载所需包裹:sf #>链接到GEOS 3.9.0、GDAL 3.2.1、项目7.2.1 tif尺寸:XY #>bbox:xmin:288937.2 ymin:9112173 xmax:298589.9 ymax:9120349 #>预计CR:南半球UTM 25区 #>L7_ETMs.tif几何图形 #>164分(294613.49117565) #>2 72分(

使用
stars
软件包,可以使用
st_extract()
函数从定义位置的光栅中提取值

库(星型)
#>正在加载所需的包:abind
#>装载所需包裹:sf
#>链接到GEOS 3.9.0、GDAL 3.2.1、项目7.2.1
tif尺寸:XY
#>bbox:xmin:288937.2 ymin:9112173 xmax:298589.9 ymax:9120349
#>预计CR:南半球UTM 25区
#>L7_ETMs.tif几何图形
#>164分(294613.49117565)
#>2 72分(295130 9117225)
#>394分(298589.9 9116806)
#>4 86分(296430.2 9112864)
#>587分(297481.9115176)
#>6110分(288937.2912173)
#>7 63分(290966.6 9116890)
#>8 84分(295595.5 9116938)
#>9 73分(291047.1 9120349)
#>1065分(294525.2917110)
我想做的是在这些点周围使用一个缓冲区,然后提取缓冲区中的
平均值。
创建缓冲区

10个特征的多边形几何体集 #>几何体类型:多边形 #>尺寸:XY #>bbox:xmin:288837.2 ymin:9112073 xmax:298689.9 ymax:9120449 #>预计CR:南半球UTM 25区 #>前5个几何图形: #>多边形((294713.49117565,294713.39117560,2。。。 #>多边形((295230 9117225,295229.8 9117220,295。。。 #>多边形((298689.9 9116806,298689.8 9116800,2。。。 #>多边形((296530.29112864,296530.19112859,2。。。 #>多边形((297581.9 9115176,297581.8 9115171,2。。。
这里的问题是,
st_extract()
函数只使用点,而不使用多边形

st_提取(r[,,1],多边形)
#>st_extract.stars(r[,1],多边形)中出错:all(st_维度(pts)==0)不正确
有没有办法提取多边形下的信息


由(v1.0.0)创建于2021-02-19这可以通过聚合来完成:

库(星型)
tif=system.file(“tif/L7_ETMs.tif”,package=“stars”)
r=读数(tif)
pnt=st_样本(st_as_sfc(st_bbox(r)),10)
多边形=st_缓冲区(pnt,距离=100)
#提取每个多边形的平均值
x=骨料(r、多边形、平均值)
圣阿斯福(x)
请记住,如果多边形之间存在重叠(与本例不同),则每个光栅值仅在其所在的第一个多边形中“计数”一次(以符合聚合的一般行为)

## Simple feature collection with 10 features and 6 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: 289038 ymin: 9111186 xmax: 298491.2 ymax: 9120605
## projected CRS:  UTM Zone 25, Southern Hemisphere
##    L7_ETMs.tif.V1 L7_ETMs.tif.V2 L7_ETMs.tif.V3 L7_ETMs.tif.V4 L7_ETMs.tif.V5
## 1        87.80556       78.38889       87.44444       69.13889      124.05556
## 2        59.31579       43.94737       33.34211       70.76316       63.65789
## 3        78.33333       64.25641       62.56410       57.00000       70.79487
## 4        70.87179       57.89744       55.35897       63.94872       88.87179
## 5        89.51282       78.12821       86.00000       64.28205      111.48718
## 6        83.28205       67.46154       67.38462       51.38462       86.12821
## 7        80.27027       70.81081       72.59459       77.51351      103.78378
## 8        74.91892       60.75676       54.05405       85.86486       90.00000
## 9        68.56410       59.74359       55.10256       78.23077       94.41026
## 10       74.86486       60.91892       62.35135       58.91892      102.29730
##    L7_ETMs.tif.V6                       geometry
## 1        98.41667 POLYGON ((295003.7 9116093,...
## 2        31.55263 POLYGON ((290092.1 9119590,...
## 3        50.64103 POLYGON ((294767 9112633, 2...
## 4        61.38462 POLYGON ((289238 9114301, 2...
## 5        90.94872 POLYGON ((298491.2 9120505,...
## 6        69.41026 POLYGON ((289770 9111286, 2...
## 7        73.64865 POLYGON ((294775.7 9117676,...
## 8        57.78378 POLYGON ((294266.6 9113127,...
## 9        56.92308 POLYGON ((293838.6 9118091,...
## 10       77.51351 POLYGON ((290557.6 9114384,...