如题。我的数据是用shapefile边界地图结合一个栅格数据文件来统计。
shapefile 名字(inZoneData) ： state.shp
栅格数据(InVa1)：temp.tif
zoneField (字段)：ID
方法是：加总(sum)
找到的资料显示R有个zonal() function, 但是只应用于栅格数据(RASTER):
> zonal(x, z, fun='mean', digits=0, na.rm=TRUE, ...)
#Compute zonal statistics, that is summarized values of a Raster* object for **each "zone" defined by a RasterLayer.**
并且这个function并没有zone field的选项，所以并不适用。
我用python的code是:
outZonalStatistics = ZonalStatistics(inZoneData, zoneField, InVa1, "SUM", "NODATA")
因为课堂要求必须使用R，所以想了解清楚这个问题。至于R的reticulate包，暂时不考虑 。
有熟悉R的朋友可以帮忙解答一下吗？谢谢！！
zonal(x, z, fun='mean', digits=0, na.rm=TRUE, ...)
#Compute zonal statistics, that is summarized values of a Raster* object for each "zone" defined by a RasterLayer.