Question

我想使用'cdo'来使用南美的另一个NetCDF从降水NetCDF数据集中提取数据。我尝试了多个过程，但总是遇到一些错误（例如网格大小不同，不支持的通用坐标等）。

我尝试过的代码：

cdo mul chirps_2000-2015_annual_SA.nc Extract feature.nc output.nc
# Got grid error

cdo -f nc4 setctomiss,0 -gtc,0 -remapcon,r1440x720 Chirps_2000-2015_annual_SA.nc CHIRPS_era5_pev_2000-2015_annual_SA_masked.nc
# Got unsupported generic error

可以从以下link下载数据。

Answer 1

我确信您可以做出/找到更优雅的解决方案，但是我只是结合了Python和shell可执行文件cdo来完成任务（有时/在某处调用子进程可被视为一种坏习惯）

#!/usr/bin/env ipython
import numpy as np
from netCDF4 import Dataset
import subprocess
# -------------------------------------------------
def nc_varget(filein,varname):
    ncin=Dataset(filein);
    vardata=ncin.variables[varname][:];
    ncin.close()
    return vardata
# -------------------------------------------------
gridfile='extract_feature.nc'
inputfile='precipitation_2000-2015_annual_SA.nc'
outputfile='selected_region.nc'
# -------------------------------------------------
# Detect the start/end based on gridfile:
poutlon=nc_varget(gridfile,'lon')
poutlat=nc_varget(gridfile,'lat')

pinlon=nc_varget(inputfile,'lon')
pinlat=nc_varget(inputfile,'lat')

kkx=np.where((pinlon>=np.min(poutlon)) & (pinlon<=np.max(poutlon)))
kky=np.where((pinlat>=np.min(poutlat)) & (pinlat<=np.max(poutlat)))
# -------------------------------------------------
# -------------------------------------------------
commandstr='cdo selindexbox,'+str(np.min(kkx))+','+str(np.max(kkx))+','+str(np.min(kky))+','+str(np.max(kky))+' '+inputfile+' '+outputfile
subprocess.call(commandstr,shell=True)

数据中的问题是文件“ precipitation_2000-2015_annual_SA.nc”目前未指定网格-变量lon，lat是通用的，因此网格是通用的。否则，您可以使用其他运算符代替selindexbox。文件extract_feature.nc更接近标准，因为变量lon，lat也具有名称和单位属性。

使用CDO提取仅针对特定区域的数据集

1 个答案: