如何在給定的邊界框中（lon1，lon2，lat1，lat2）內獲得時間平均值的平均值？

Question

我使用下面的代碼繪制了EDIR變量隨時間變化的面積平均值。 但是，我想知道如何像邊界框一樣設置最小和最大緯度/經度，以便在設置緯度/經度框中繪制EDIR變量的面積平均值。 謝謝！

 ld0 = xr.open_mfdataset('path/*LDASOUT*', condcat_dim='time') t2=ld0.EDIR.mean(dim='x','y')) t2.plot.line('-',color='black',linewidth=3,label='control_run') plt.legend() plt.show 

Answer 1

如果lat和lon是數據集中的坐標變量，則可以使用sel方法。

da = ld0['EDIR'].sel(lon=slice(min_lon, max_lon), lat=slice(min_lat, max_lat))

可以在xarray文檔中找到使用xarray基於坐標的選擇的更多示例： http ://xarray.pydata.org/en/stable/indexing.html

Answer 2

從這里的文檔中，您似乎可以使用t2.where(t2.x < lon2, drop=True)選擇坐標邊界。 我在用手機打電話，所以現在無法測試，但類似：

temp = ld0.where( ( (t2.x > lon1)  &  (t2.x < lon2) & (t2.y > lat1)  &  (t2.y < lat2)  ), drop=True)

應該給你合適的口罩。 然后，您可以計算並繪制此數組的空間平均值。