问题描述
我正在尝试从可通过 opendap 获得的全球 NetCDF 文件中提取地理子区域。所需区域从西经 13 度到东经 10 度:经度 -13 到 10。
NCO 通常可以这样做:
ncks -d lat,40.,70.-d lon,-10.,13. https://dataserver.nccs.nasa.gov/thredds/dodsC/CMIP5/NASA/GISS/rcp85/E2-H_rcp85_r2i1p1_day/tos_day_GISS-E2-H_rcp85_r2i1p1_20510101-20751231.nc out.nc
但是,生成的 NetCDF 文件中的经度从 0 到 13 度,而不是 -10 到 13。这是因为经度被记录为“东经度”,因此 NCO 似乎不处理负值。有没有人知道调整上述内容以接受否定的方法?
注意:CDO 可以处理负经度。然而,在这种情况下,CDO 不是一个合适的解决方案,因为它 a) 在裁剪之前必须先下载全局数据集,并且 b) 无法处理这个特定的 opendap 文件。
解决方法
NCO 如here 所述处理这些“包裹坐标”。请重新尝试使用 [0,360] 中的边界框限制,例如
ncks -d lat,40.,70. -d lon,347.,10. https://dataserver.nccs.nasa.gov/thredds/dodsC/CMIP5/NASA/GISS/rcp85/E2-H_rcp85_r2i1p1_day/tos_day_GISS-E2-H_rcp85_r2i1p1_20510101-20751231.nc out.nc