2017-03-03 39 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主頁介紹如何克隆/獲取最新快照。 –