我在terra中使用行作为SpatVect对象,因此我更喜欢terra::rasterize解决方案。我也会对使用stars::st_rasterize的解决方案感到满意(也请参阅问题)。尽管如此,我还是尝试过使用terra::rasterize,参见下面的示例。在这里您可以看到,当通过max进行栅格化时,mean函数和terra::rasterize(... fun = "max"/"mean
在gdal 2.1中使用gdal.Rasterize时,是否可以设置数据类型(字节、Float32)?目前,我使用gdal.Translate来转换为字节,但是效率很低。/solution.tifgdal.Rasterize(tiftemp, shp, options=rasterizeOptions)
#translate to Byte data type (not supporte