Xarray在时间上找到每个网格点的最大值

Question

我有以下NETCDF文件，其中包含在一个地区的每日降雨（爆炸），并且我以Xarray开放。

<xarray.Dataset>
Dimensions:    (latitude: 500, longitude: 600, time: 120)
Coordinates:
  * latitude   (latitude) float32 -35.22 -35.17 -35.12 ... -10.38 -10.33 -10.28
  * longitude  (longitude) float32 10.27 10.32 10.38 10.43 ... 40.12 40.18 40.22
  * time       (time) datetime64[ns] 2022-01-01 2022-01-02 ... 2022-04-30
Data variables:
    precip     (time, latitude, longitude) float32 ...
Attributes: (12/15)
    Conventions:       CF-1.6
    title:             CHIRPS Version 2.0
    history:           created by Climate Hazards Group
    version:           Version 2.0
    date_created:      2022-05-16
    creator_name:      Pete Peterson
    ...                ...
    reference:         Funk, C.C., Peterson, P.J., Landsfeld, M.F., Pedreros,...
    comments:           time variable denotes the first day of the given day....
    acknowledgements:  The Climate Hazards Group InfraRed Precipitation with ...
    ftp_url:           ftp://chg-ftpout.geog.ucsb.edu/pub/org/chg/products/CH...
    website:           http://chg.geog.ucsb.edu/data/chirps/index.html
    faq:               http://chg-wiki.geog.ucsb.edu/wiki/CHIRPS_FAQ

我想在每个网格点（即每个LAT/LON点）拉出最大降水值。我需要最大值才能保持其纬度，经度和时间的属性，因为然后我将分析这些最大值何时何地。

我知道我可以使用以下内容在第一个网格点跨时间获得最大值。但是，我该如何查看最大发生在哪里以及在什么日期发生？

ds1a.precip[:,0,0].max()

然后，我真的很努力地为所有网格点这样做。我真的很感谢指导。

Answer 1

一般而言

• ​A>找到最大值
• 代码> 找到达到最大值的坐标。

两种方法都有可选的参数dim，以指定要找到最大值的维度。

在您的情况下，

ds1a.max(dim='time')

每个latitude和pentitude的数据集将返回一个数据集，并将

ds1a.argmax(dim='time')

返回达到此最大值的时间。

这是一个可复制的示例，其中包括示例数据：

import xarray as xr

ds = xr.tutorial.open_dataset('air_temperature').load()
ds.max(dim='time')
ds.argmax(dim='time')