R子设置netcdf文件

Question

我正在使用来自此来源的2018年的lulc封面数据

https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=form

数据格式为netcdf，并显示每个时期的土地覆盖分类。为了从netcdf读取和创建栅格，我这样做：

library(ncdf4)
temp <- nc_open(file.name)
lon <- ncvar_get(temp, "lon")
lat <- ncvar_get(temp, "lat")
lccs <- ncvar_get(temp, "lccs_class", signedbyte = FALSE)

但是，netcdf文件太大，需要很长时间才能读取该文件。我只需要lat lat的一个子集，它的边界框定义如下：

min lat: 8.125
max lat: 37.125
min lon: 68.125
max lon: 97.375 

如何使用这些边界框设置上述netcdf的子集？

Answer 1

在大多数情况下，你可以这样做

library(raster)
b <- brick("filename.nc")
e <- extent(8.125, 37.125, 68.125, 97.375)
x <- crop(b, e)