如何解决R gdalUtils::gdalwarp 函数在粗化单元分辨率时跳过单元
我将 fasterized
一个 sf
多边形转换为 30m 像元分辨率栅格。我现在试图将 30m 栅格聚合到 150m。 raster::aggregate
函数很慢,所以我尝试使用 gdalUtils::gdalwarp
。当我粗化分辨率时,gdalwarp
输出会跳过单元格。有办法解决吗?
input <- ".../globalIP/alp/footprint/pl30m.tif"
output <- ".../globalIP/alp/footprint/pl150m.tif"
gdalUtils::gdalwarp(srcfile = input,dstfile = output,tr = c(150,150),t_srs='+proj=utm +zone=11 +datum=NAD83',overwrite=TRUE,r="bilinear")
图 1. 绿色单元格 = 30m fasterize
栅格。黑色单元格 = 150m gdalwarp
栅格
解决方法
最终使用了 terra::aggregate()
。仍然很好奇如何在 gdal 中重新采样以便单元格不会跳过...
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。