極坐標投影到規則的緯度/經度網格

Question

這是我上一個問題的后續，因為示例和數據更加具體。 我有一個 304 x 448 的冰數據陣列（在極地立體投影中），其中的 netcdf 可以在這里獲得。 每個網格框對應的緯度/經度點可以在這個文件中找到。 我想將此數據轉換為常規的緯度/經度網格（例如 180x360），其中每個網格框都是一個度數。 有沒有python的方法來做到這一點？ 使用 Basemap 等，還是使用 CDO（氣候數據運營商）？

提前謝謝了。 詹姆士

Answer 1

如果您首先將網格添加到數據文件中，例如，

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

然后你可以使用 NCO 的ncremap來重新網格，例如，

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

不幸的是，這將需要 NCO 的當前快照（而不是已發布的版本），因為輸入文件包含的維度順序不同尋常，我剛剛修補了 NCO 以進行處理。 重新網格化輸出的示例是here 。

在 OP 的要求下，我重新整理了他的文件，然后上傳了它（這里）。 我沒有一個 1x1 數據集，所以我首先創建了一個 1x1 網格文件，如手冊中所示，然后用它重新網格化（使用 ncremap -g 選項）：

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

Answer 2

通常在 CDO 中執行此操作的命令是

cdo remapcon,r360x180 in.nc out.nc 

代替 con=conservative 重映射，您可以使用 bil=bilinear, nn=nearest neighbor, con2=2nd order 保守重映射。

根據這篇文章， https ://code.mpimet.mpg.de/boards/1/topics/8302?r = 8326 最新版本的 cdo 應該能夠處理極地立體投影。