私はpythonが初めてで、[時間、高さ、経度、緯度]の次元を持つnetcdfファイルからnumpy配列を取得し、指定された緯度と経度に関数を補間する関数を作成しようとしています。私は scipy.interpolate を調べましたが、そこからどこへ行くべきかまだわかりません。助けはありますか？
質問する
3819 次
2 に答える
2
(a) データ セット全体を新しい緯度/経度座標のセットに補間するか、(b) 単一の緯度/経度位置で値を取得するか、質問からはわかりません。
いずれにせよ、Python でこれを行うことに結婚していない場合は、 Climate Data Operatorsで remapgrid を使用する方がはるかに簡単です (必ず netCDF サポートを使用してインストールしてください)。これには、netCDF ファイルが CF 規則に従う必要があります ( WRF ファイルを CF に準拠させる方法については、この記事を参照してください)。
remapgrid のドキュメントは、こちらにあります。特に、グリッド オプションに関する情報についてはセクション 1.3 (pp. 9-12) を、補間手法に関する情報についてはセクション 2.10 (p.106 から始まる) を参照してください。再マップ方法にはいくつかのオプションがあり、使用するコマンドはアプリケーションによって異なります。
(a)ガウス (128x64) グリッドへの双一次補間を使用した再マッピングの例を次に示します。
cdo remapbil,n32 wrfout_d01_1999-01-01-01_00:00:00.nc out.nc
(b) 1 つの緯度経度ポイントへの最近傍マッピングを使用した再マッピングの例を次に示します。
cdo remapnn,lon=-107.0_lat=34.0 wrfout_d01_1999-01-01-01_00:00:00.nc out.nc
于 2015-05-21T16:10:56.030 に答える