我的 NetCDF 文件有 4 个维度（经度、纬度、时间、水平压力）和几个变量。
我想将经度数据 [340 342 344 ... 360 0 2 4 ... 18 20] 替换为 [-20 -18 -16 ... -2 0 2 4 ... 18 20]。事实上，经度是从 0° 到 360° 计算的，但我需要从 -180° 到 180° 计算。
我在互联网上找到了一些想法，但我不知道如何应用它们。
这里有一个解释：
“一些 netcdf 文件包含严格的正值 [0,360]，这与 gdal/QGis 处理经度值的方式不兼容 - 它们不会环绕日期线。这也会导致 gdalwarp 以及与其他数据集结合时出现错误（与 [-180,180] ] 值），由于 180 度经度偏移，它们未正确对齐。
简单的修复方法是测试经度变量，并通过在 IReadBlock 和 SRS 检测中减去 180，将 [180,360] 区间内的任何值转换为 [-180,180]。”
我还读过有关“gdal-translate”的内容，但我不知道如何使用它。
编辑：我也遇到了纬度问题，即 N -> S （我想要它 S -> N），感谢“cdo invertlat”，我将其反转
这是我发现的：
ncap2 -O -s 'where(lon>180) lon=lon-360' ifile ofile
我还不确定它不会改变任何其他东西，但它似乎有效。
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
