问题描述
class : SpatRaster
dimensions : 45000,45000,1 (nrow,ncol,nlyr)
resolution : 0.0008888889,0.0008888889 (x,y)
extent : 20,60,-40,0 (xmin,xmax,ymin,ymax)
coord. ref. : +proj=longlat +datum=wgs84 +no_defs
data source : N00E020_agb.tif
names : N00E020_agb
我想将所有 29 个图块合并为一个栅格,但是,在这种高分辨率下需要很长时间。所以我首先使用聚合降低了这些栅格的分辨率。我按因子 113 聚合(大致导致分辨率为 0.1):
pathz <- "C:/Users/rastertiles/" # directory where raster tiles are located
filez <- list.files(path = paste(pathz,sep = ""),pattern = ".tif") # list all files in directory
r.list_113 <- list() # make an empty list
for(i in 1:length(filez)){
tmp1 <- rast(paste(pathz,filez[i],sep = ""))
tmp1 <- terra::aggregate(tmp1,fact=113,fun="mean",na.rm=TRUE) # aggregate by a factor 113.
r.list_113[[i]] = tmp1
}
然而,虽然原始切片都具有相同的原点 ([0,0]),但在使用“聚合”后,新栅格的原点略有不同:
> lapply(r.list_113,origin)
[[1]]
[1] 0.01155556 0.00000000
[[2]]
[1] 0.03466667 0.00000000
[[3]]
[1] -0.04266667 0.00000000
[[4]]
[1] -0.01955556 0.00000000
(...)
然后当我尝试将新栅格与 'merge' 合并时,出现错误:
> m_113X <- do.call(terra::merge,r.list_113X)
Error: [merge] origin of SpatRaster 2 does not match the prevIoUs SpatRaster(s)
欢迎任何我应该做的建议。
解决方法
您使用的是什么版本的 terra
?我问是因为我相信当前版本会发出警告但继续。此外,您当然可以使用 100 而不是 113 进行聚合,这样来源就不会改变。
他们最简单的方法可能是制作这样的虚拟栅格:
v <- vrt(files,"my.vrt")
然后继续
r <- rast("my.vrt")
也许
a <- aggregate(r,100)
或者任何你想做的事情。