import xarray as xr
# 用xarray读取nc文件数据
data = xr.open_dataset(r"D:\论文\论文数据\nc文件\2017——2020降水数据.nc")
# 读取经纬度
lat = data.latitude
lon = data.longitude
# 选取开始时间
start_date='2017-12-01 00:00:00'
# 选取结束时间
end_date='2021-02-28 23:00:00'
# 选择气压场
levels=500 # 500hPa
# 选择数据,经纬度,时间,压强
rain =data['tp'].sel(longitude=(105.0,105.8), latitude=( 23, 22.75),time=slice(start_date, end_date)).loc[dict(level=levels)]
"""
另一种裁剪文件的方法
time=d.loc['1991-01-01':'2020-12-01']
# 选取特定变量
ds1 = xr.Dataset({"pr_wtr ":time})
# 导出文件
ds1.to_netcdf(r"D:\论文\论文数据\nc文件\newpr_wtr.mon.mean.nc")
"""
有时候下载数据会下载很多时间点很多气压场的数据,通过筛选特定经纬度范围,时间和特定气压场的数据进行数据处理。
最后一句代码,可以选择删除一些,例如只对纬度进行筛选,则可以去除关于经度的代码。