如何使用 ferret 或 cdo 遮盖盆地

how to mask basin using ferret or cdo

我在 netcdf 文件中有世界地图 GRACE 数据。我正在寻找选定区域的可视化数据,例如只是可视化亚马逊或伏尔加河盆地的数据。有人建议我使用面盆面膜

如果您有盆地的 shapefile,您可以执行以下操作。

# convert shp to nc
gdal_rasterize -of netCDF -burn 1 -tr 0.01 0.01 basin.shp basin.nc

# regrid
cdo remapbil,infile.nc basin.nc basin2.nc

# mask
cdo ifthen basin2.nc infile.nc masked_infile.nc

如果您不需要完全遵循盆地的边界并且乐于简单地查看“矩形”经纬度区域,那么您也可以简单地做

cdo sellonlatbox,lon1,lon2,lat1,lat2 world_file.nc region.nc