2017-03-03 326 views
0

这是我上一个问题的后续内容,因为这个例子和数据更具体。我有一个304乘448阵冰数据(在极坐标投影中),其中netcdf可以得到here。每个网格框的相应经/纬度点可以在file中找到。我想将这些数据转换成一个常规的纬度/经度网格(比如180x360),其中每个网格框是一个度数。有没有python的方式来做到这一点?等使用底图,或者使用CDO(气候数据运算符)?极地立体投影到经常纬度/经度网格

非常感谢提前。 詹姆斯

回答

1

如果第一网格中添加与数据文件,例如,

ncks -A -v lat,lon NSIDC_LatLon.nc sample.nc 

然后你可以使用NCO的ncremap到REGRID与,例如,

ncremap -i sample.nc -d 1x1.nc -o out.nc 

不幸的是这将需要NCO(而不是发布版本)的当前快照,因为输入文件包含的维度与我刚刚修补NCO处理的异常顺序相关。重新输出的样本是here

根据OP的要求,我重新登录了他的文件然后上传了它(here)。我没有放置1x1数据集,因此我首先创建了1x1网格文件(如手册中所示),然后使用ncremap -g选项对其进行了修改:

ncremap -i ~/NSIDC.nc -g ${DATA}/grids/180x360_SCRIP.20150901.nc -o ~/NSIDC_1x1.nc 
+0

非常感谢您的答复。有没有一种方法可以获得NCO的最新快照(作为tar.gz文件),以便我可以尽早实施此解决方案?谢谢。 –

+0

另外,什么是1x1.nc文件? –

+0

正如文档解释,1x1.nc是已经在所需的1x1网格上的任何文件。它会被读取,但不会被写入。 NCO在GitHub上。 NCO主页介绍如何克隆/获取最新快照。 –