R子设置netcdf文件答案

【问题标题】：R sub-setting netcdf fileR子设置netcdf文件
【发布时间】：2020-09-16 10:01:03
【问题描述】：

我正在使用来自该来源的 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 文件太大，读取文件需要很长时间。我只需要一个纬度子集，其边界框定义如下：

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

如何使用这些边界框对上述 netcdf 进行子集化？

【问题讨论】：

标签： r r-raster netcdf4

【解决方案1】：

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

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

【讨论】：