問題タブ [python-xarray]

私は次元、、および でDataArray呼び出す小さな python-xarrayを持っています。次元、および変数、および_ _TINYxytTINYDataSetHUGEx'y't'foobarbaz

座標を使用してデータセットを選択する方法を探しています (小さな座標と座標のペアはすべて巨大に存在します)。 HUGETINY

私は急いで次のイラスト を書き上げました

HUGE選択は、 のを共有するcoordsのサブセットを返しますTINY。返されるのはDataSet. これを処理する簡単な x-array メソッドはありますか?

nl4124 個の測定値を含む配列があるとします。それぞれは、測定が行われた場所を指定する( lat, ) カップルに関連付けられています。lonこれらの位置はグリッド化されていません。つまり、等間隔の値に整列していません。

とを座標としてnl指定して、 の DataArray を作成します。latlon

たとえば、これらの値を経度または緯度のビンにグループ化して、それらを操作できることを知っています。

私がやりたいのは、経度 x 緯度の 2D ビンの値をグループ化することです。これにより、結果をマップとして表示できます。groupby_bins ではそれができないと思いますが、別の解決策はありますか?

例で更新：

これは、適切なnumpyで私がやりたいことをする方法です:

xarray とダッシュの統合を維持するために numpy に頼るのは避けたいと思います。

3 次元の DataArray とPanel.ffill(axis=None, inplace=False, limit=None, downcast=None)同等のものがあるとすれば、それは 何ですか?Panel.fillna(value=None, method=None, ...)

DataArray には fillna メソッドがありますが、パラメーターとして値のみを受け取ります。DataArray.fillna(value).

例: xarray の既存のメソッドで以前の非 nan 値を使用して nan 値を埋める方法は?

たとえば、次を使用してDJF降水量の値を選択したいとします。

その期間から値を取得するためだけに、全期間にわたって平均したくありません(その.mean('time')部分は必要ありません) DJF。

これは、ND データを取得する方法です ( funcIRL はベクトル化できません)。

自分のしていることは単純化できると感じています。データを再形成せずにこれを入れたいと思いxarray.DataArrayます。しかし、これは私が今それを行う方法です：

itertools.productこれは単純な例ですが、通常、 (並列に)をマッピングして取得した ~10D データを使用します。

私の質問: データを再形成せずに、 、、およびvalsの長さを使用および使用せずにこれを行うにはどうすればよいですか?xsyszs

あなたが何をするかと同様の方法で：

編集： @hpauljの答えに触発されて、これが私がそれを解決した方法です、ありがとう！

EDIT 2 この問題を参照してください: https://github.com/pydata/xarray/issues/1914

landcover タイプの dataArray があります。リストにある特定の値をマスクしたいと思います。複数の条件で xr.where() 関数を使用することは可能ですか?

これにより、次のことが得られます。

上記は 2 つの土地被覆の値に対して機能しますが、30 のタイプに対してこれを行うのは面倒です。より良い方法はありますか？

世界中の気温を含むこのサンプル データセットと、さらに重要なことに、land陸上/非水域をマークする mask があります。

海をマスクしてプロットできるようになりました

dry_areasこのように見えます

すべての非 nan 値の座標を取得するにはどうすればよいですか?

dry_areas.coords(55, 143)バウンディングボックスが表示され、マスクを適用できるように緯度と経度を形状に入れることができません。

私が見つけた唯一の有効な回避策は dry_areas.to_dataframe().dropna().reset_index()[['lat', 'lon']].valuesです。

これは非常に単純だと思いますが、明らかにnumpy/matrixの忍者ではありません。

これまでの最善の解決策

これは私がこれまでに持ってきた最短のものです: