問題タブ [python-iris]

これは少し長いショットですが、自分でこれを書く前にここで質問したいと思いました。

経度、緯度、高さの座標を持つ 3D 立方体のデータがあります。経度、緯度、高さ、データのすべての点で、データの 4 つの 1D ベクトルが必要です。これは、ポイントとその位置のリストとして ASCII ファイルに書き込むことができるようにするためです。データに対してこれを行うのは、形状を変更することで簡単ですが、よりトリッキーな部分は、座標を正しいベクトルに変換することです。

誰かがすでにこれを行っていて、いくつかのヒントがありますか?

NetCDF ファイルからの降水データに基づいて、地域の降水レベルのマップを作成しました。降水量が 800mm 未満の場合は 1 つの色、800 ～ 1000mm の場合は別の色になるように、カスタム スケールを追加したいと思います。http: //www.metmalawi.com/climate/climate. php

現在、グラデーション スケールを使用していますが、必要な詳細が表示されません。これは、現時点でのプロットのコードです (「平均」は、既にフォーマットしたデータです)。

どうすればこれができるか知っている人はいますか？

ありがとうございました！

いくつかのポリゴン(カナダの州) があり、 を使用して読み取り、GeoPandasこれらを使用してマスクを作成し、2 次元の緯度経度グリッド上のグリッド データに適用したいと考えています ( を使用してnetcdfファイルから読み取りますiris)。最終的な目標は、特定の州のデータのみを残し、残りのデータをマスクすることです。したがって、マスクは、州内のグリッド ボックスの場合は 1 になり、州外のグリッド ボックスの場合は 0 または NaN になります。

ポリゴンは、こちらのシェープファイルから取得できます: https://www.dropbox.com/s/o5elu01fetwnobx/CAN_adm1.shp?dl=0

私が使用している netcdf ファイルは、 https ://www.dropbox.com/s/kxb2v2rq17m7lp7/t2m.20090815.nc?dl=0 からダウンロードできます。

ここには2つのアプローチがあると思いますが、両方に苦労しています：

1) 多角形を使用して緯度経度グリッドにマスクを作成し、これを Python 以外の多くのデータファイルに適用できるようにします (推奨)。

2) 多角形を使用して、読み込まれたデータをマスクし、関心のある州内のデータのみを抽出して、対話的に操作します。

これまでの私のコード:

助けやアドバイスをありがとう。

更新: 以下の @DPeterK の優れたソリューションに従って、元のデータをマスクして、次のようにすることができます。