裁剪多个光栅（ASCII）

Question

我有50多个需要裁剪的光栅文件（ASCII格式）。我已经以ASCII格式从ArcMap导出了遮罩，并将其加载到R中。如何使其适用于一行中的所有光栅，并以与之前相同的名称导出它们（当然是在不同的文件夹中，以避免覆盖）？

我知道光栅软件包中有裁剪功能，但到目前为止我从未使用过。我只是把它们堆放起来做进一步的栖息地分析。

到目前为止，我的代码：

#### Use only part of area
files2 <- list.files(path="D:/",full.names=TRUE, pattern = "\\.asc$")
files2
# Create a RasterLayer from first file
mask_raster <- raster(files2[1])
# Crop. But how??   
crop(x = , y=mask_raster)
writeRaster(...)`

Answer 1

我没有找到一个简单的解决方案来通过光栅裁剪多个光栅，而是通过一个形状文件。所以我回到ArcMap并将光栅转换为一个shapefile。然后在R中，裁剪和遮罩是关键步骤。参见下面的代码（根据Mauricio Zambrano Bigiarini的代码修改）。希望这有帮助。

# Reading the shapefile (mask to crop later)
Maskshp <- readOGR("mask.shp")

# Getting the spatial extent of the shapefile
e <- extent(Maskshp)

# for loop with (in this case) 60 runs
for (i in 1:60) {

# Reading the raster to crop
files <- list.files(path="...your path",full.names=TRUE, pattern = "\\.asc$")
Env_raster <- raster(files[i])
# Save the filename
filename <- (paste(basename(files[i]), sep=""))

# Crop the raster
Env_raster.crop <- crop(Env_raster, e, snap="out")

# Dummy raster with a spatial extension equal to the cropped raster,
# but full of NA values
crop <- setValues(Env_raster.crop, NA)

#  Rasterize the catchment boundaries, with NA outside the catchment boundaries
Maskshp.r <- rasterize(Maskshp, crop)

# Putting NA values in all the raster cells outside the shapefile boundaries
Maskshp.masked <- mask(x=Env_raster.crop, mask=Maskshp.r) 
plot(Maskshp.masked)

#Export file to working directory with original name as new name
writeRaster(Maskshp.masked, filename)
}