R语言地图截取：更快替代mask函数的crop和subset

在R语言中，截取地图可以使用raster包中的crop函数或者terra包中的subset函数代替mask函数。这两个函数的运行速度通常比mask函数快。

例如，使用raster包中的crop函数可以这样实现：

library(raster)
# 导入地图
map <- raster('map.tif')
# 定义要截取的区域
xmin <- -10
xmax <- 10
ymin <- -10
ymax <- 10
# 使用crop函数截取地图
map_cropped <- crop(map, extent(xmin, xmax, ymin, ymax))

使用terra包中的subset函数可以这样实现：

library(terra)
# 导入地图
map <- rast('map.tif')
# 定义要截取的区域
xmin <- -10
xmax <- 10
ymin <- -10
ymax <- 10
# 使用subset函数截取地图
map_cropped <- subset(map, xmin, xmax, ymin, ymax)

需要注意的是，crop函数和subset函数的参数顺序略有不同，需要根据实际情况进行调整。