分割栅格，定义新栅格的名称答案

【问题标题】：Split raster, define names of new rasters分割栅格，定义新栅格的名称
【发布时间】：2016-04-01 08:44:22
【问题描述】：

我这几天一直在寻找解决方案，我如何在 ArcMap 中将大型地理参考栅格拆分为较小的栅格。当然，还有“分割栅格”工具，它可以将栅格按多边形中的要素进行分割，但栅格名称是随机分配的。相反，我需要定义新栅格的名称 - 我使用从坐标网格剪切的 shapefile 拆分栅格，并且栅格的名称应该与网格页码相同。我根本没有用过python，所以有没有可能在那里完成？

【问题讨论】：

标签： split raster

【解决方案1】：

如果您使用：Data Management Tools -> Raster -> Raster Processing -> Clip 会发生什么？结果是否不包含生成的栅格名称中的要素 ID 或类似名称？我想这对提示很有用：https://gis.stackexchange.com/questions/84969/clip-raster-using-shapefile-into-multiple-rasters-using-arcgis-desktop 不记得在 shp 文件中使用渔网裁剪光栅时在文件名中提供特征 ID 的方法。

【讨论】：

谢谢-您共享的链接中描述的使用模型构建器和“剪辑栅格”工具的方法似乎解决了我的问题！
很高兴知道。虽然有趣的是，大多数人发现生成渔网比剪裁本身更具挑战性。

【解决方案2】：

这是一个示例，说明如何使用 R 做到这一点

library(raster)
# example data
r <- raster()
values(r) <- 1:ncell(r)

# create blocks
a <- aggregate(raster(r), 60)
b <- as(a, 'SpatialPolygons')

f <- function(i) {
   fname <- paste0('raster_', i, '.tif')
   crop(r, extent(b[i,]), filename=fname)
}

x <- sapply(1:length(b), f)

【讨论】：