我在用这个来源的2018年的lulc封面数据。
https:/cds.climate.copernicus.eucdsapp#!datasetsatellit-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 lon的子集，其边界框定义如下。
min lat: 8.125
max lat: 37.125
min lon: 68.125
max lon: 97.375
我怎样才能用这些边界框对上面的netcdf进行子集？
在大多数情况下，你可以像这样做
library(raster)
b <- brick("filename.nc")
e <- extent(8.125, 37.125, 68.125, 97.375)
x <- crop(b, e)