問題タブ [netcdf4]

不規則なグリッド (ランバート コニカル) 上のデータを通常のグリッドに再グリッド化する必要があります。pyresample が私の最善の策だと思います。実際、私の元の緯度、経度は1Dではありません(basemap.interpまたはscipy.interpolate.griddataを使用するために必要なようです)。

このSOの回答は役に立ちました。ただし、空の補間データが得られます。これは、影響範囲の選択と、データがラップされている (??) ことに関係していると思います。

これは私のコードです:

lat_n と lon_n は不規則で、緯度と経度は投影された座標 x、y に対応します。

lon_n のやり方のため、次のように追加しました。

私がそれらをプロットすると、見栄えがよく、問題ないようになります。

次に、通常の座標の新しいセットを作成します。

上記の回答に続いて、次のコードを書きました。

pyresample.geometry から SwathDefinition をインポート pyresample.kd_tree から resample_nearest をインポート

データの解像度は約 30km なので、70km と入力しました。入力した fill_value はデータから取得したものですが、もちろんゼロまたは nan を入力することもできます。

ただし、空の配列を取得します。

私は何を間違っていますか？また、別の方法がある場合は、それを知りたいと思っています。Pyresample のドキュメントは少し薄いので、もう少し助けが必要です。

別の griddata 関数を使用することを提案するこの回答を見つけました:

それは大丈夫だと思われます：

しかし、非常に強力なように見えるので、pyresample についてもっと理解したいと思います。

次のように、netcdf4 で変数を作成し、属性をコピーできます。

xarray の対応バージョンは？

以前にラスター スタックから書き込まれたように、netcdf ファイルからラスター レイヤー名を抽出しようとしています。ラスター スタックを ncdf にエクスポートすると、正常に動作します。例えば：

ただし、ncdf ファイルを R に読み込むと、Z 次元 (年) は保持されません。例えば：

したがって、得られるのはバンド番号、つまり 1、2、3 です。しかし、私が必要とするのは、以下で定義されている年 (例: 2014,2015,2016) で定義されたテキストです。

これを行うことは可能ですか?? この問題は新しいものではありません。こちらを参照してください: https://gis.stackexchange.com/questions/122167/export-band-names-with-netcdf-file-in-r

必要なものを取得するためのいくつかの複雑な回避策がありますが、それらはほとんど効率が悪いようです (つまり、スタックをマトリックスに変換してから、ここから操作します)。大量の余分なコードを記述したり、不要な RAM を占有したりすることなく、よりエレガントな方法があるかどうか疑問に思っています。

xarray オブジェクトを netcdf ファイルに書き込もうとすると、次のエラーが発生しました。

データがメモリに対して大きすぎるため、チャンクする必要があります。
ルーチンは基本的に次のとおりです。

基本的には内容をいじって書き直すだけです。それにもかかわらず、チャンクは悪いようです。1 つの配列に再チャンクする場合と同じです。私もdsを書き換えることができません。エラーを追跡する方法、またはこれを解決する方法はありますか?

netCDF4 のバージョンは 1.2.4
、xarray (以前の xray) のバージョンは 0.8.2
、dask のバージョンは 0.10.1 です。

netCDF ファイルのグローバル メタデータ _FillValue を変更するための python netCDF4 コマンド/例はありますか? netCDF ファイル内のすべての -ve 値を置き換えようとしましたが、_FillValue 属性が設定されるまでは機能しません。

PythonモジュールnetCDF4を使用し、mpl_toolkits Basemapで画像を作成して、気象モデルの出力を分析しています。

気圧の等高線図をいくつか作成しましたが、これは非常に簡単でした。今、私は風にとげを追加したいのですが、最終的に海の上に風だけのプロットになりました. いじり回した後、ベースマップ関数 fillcontinents() が実際にこれらの風のバーブをオーバーレイしていることを発見しましたが、バーブ関数は後で呼び出されます。そのため、fillcontinents() 関数を削除すると、必要に応じてすべての風刺が得られました。

大陸を埋める方法はありますか?

私のコードは次のようになります (一部のグラフィカル パラメータは省略しましたが、質問には重要ではありません)。

https://code.google.com/archive/p/netcdf4-python/wikis/UbuntuInstall.wikiに従って netcdf4-python をインストールしようとして います。 、以下のエラーが発生しています

これをソースとして使用: Windows 7 で NCO または R を使用して、毎月の TRMM netCDF ファイルを単一の netCDF ファイルに連結する方法は?

rdsファイルを作成できました。ただし、ncファイルとして保存することに本当に興味があります。12 ステップ (各月に 1 つ) の時間ディメンションを持つ新しい netCDF ファイルを次のように作成してみました。

ここでの課題は、降水データを各月に追加することです。
最初のスクリプトに従って、ncvar_put関数を使用してみましたが成功しませんでした。

私が得たもの

ncvar_put(precip.nccreate, precip.ncvar, data2d, start = c(1), のエラー: オブジェクト 'precip.nccreate' が見つかりません

nc_create("precip.nccreate.nc", precip.ncvar, force_v4 = FALSE, のエラー: オブジェクト 'precip.ncvar' が見つかりません

とにかく、複数の netcdf ファイルを単一の netcdf に連結する簡単な方法を見つけようとしているだけだと思います。

ありがとう