重新投影栅格堆栈列表

Question

我有一个很大的栅格堆栈列表，我想重新投影然后在R中裁剪。我在ArcGIS中使用batch processing完成了相同的过程，这要快得多。 在R中进展不顺利。是否有任何改进流程的建议？

我的rstack.lst由19个栅格堆栈组成，如下所示：

class       : RasterStack 
dimensions  : 4800, 7200, 34560000, 46  (nrow, ncol, ncell, nlayers)
resolution  : 463.3127, 463.3127  (x, y)
extent      : 11119505, 14455357, -5559753, -3335852  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs 
names       : ...

我正在使用以下代码：

for (i in 1:length(rstack.lst)){
  rstack.lst[[i]] <- projectRaster(rstack.lst[[i]], crs = m.crs) # reproject
  rstack.lst[[i]] <- crop(rstack.lst[[i]], brdshp) # crop
  rstack.lst[[i]] <- mask(rstack.lst[[i]], brdshp) # mask
  print (i)
}

两个小时后连一个我都没有打印出来！