如何在 NCO 中按经度裁剪 NetCDF 文件
How to crop a NetCDF file by longitude in NCO
我正在尝试从 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
我正在尝试从 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