2017-03-03 10 views
0

これは私の前の質問からの続きです。これは例とデータによってはるかに具体的です。私は、304×448の氷データ(極座標投影)を持っています。そのうちのnetcdfはhereです。各グリッドボックスの対応する緯度/経度点はfileにあります。私はこのデータを通常の緯度/経度グリッド(180x360)に変換したいと思います。これを行うためのPythonの方法はありますか? Basemapを使用する方法、またはCDO(気候データ演算子)を使用する方法通常の緯度/経度グリッドへのポロ立体投影

事前に感謝します。 ジェームズ

答えて

1

あなたが最初に次にあなたが例えば、とのこと格子変えするNCOのncremapを使用することができます

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

、例えば、とのデータファイルにグリッドを追加した場合、

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

残念ながら、この意志入力ファイルにディメンションが含まれているためNCOの現在のスナップショット(リリース版ではなく)が必要です。 regridded出力のサンプルはhereです。

私は彼のファイルを再編集してアップロードしました(here)。私は1x1のデータセットが(ncremap -gオプションで)説明書に示されるように、私が最初にそれと再グリッド、1x1のグリッドファイルを作成したので、周りに座っていませんでした:返信用

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のホームページでは、最新のスナップショットを複製/取得する方法について説明しています。 –