R语言提取全球栅格图中每个国家的范围

要从一张全球尺度的栅格图上提取每个国家的范围，可以使用以下步骤来实现：

1. 安装并加载必要的R包：raster、sp、rgeos和rgdal。

install.packages('raster')
install.packages('sp')
install.packages('rgeos')
install.packages('rgdal')

library(raster)
library(sp)
library(rgeos)
library(rgdal)

2. 读取栅格图像文件，并将其转换为一个RasterLayer对象。

raster_file <- 'path_to_raster_file.tif'
raster_layer <- raster(raster_file)

3. 获取国家边界数据，可以使用raster::getData函数从raster包中获取。

countries <- getData('GADM', country = 'all', level = 0)

4. 将国家边界数据转换为SpatialPolygonsDataFrame对象。

countries_sp <- as(countries, 'SpatialPolygonsDataFrame')

5. 将栅格图的坐标系转换为国家边界数据的坐标系。

raster_layer <- projectRaster(raster_layer, crs = projection(countries_sp))

6. 使用raster::extract函数提取每个国家的范围。

country_ranges <- extract(raster_layer, countries_sp, fun = 'sum', na.rm = TRUE)

在这个例子中，我们使用了fun = 'sum'来计算每个国家在栅格图上的值的总和。你可以根据你的需求选择其他的提取函数。

最后，country_ranges将包含每个国家在栅格图上的范围。