如何使用带有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,...