上一篇文章《xarray指南:插值》中介绍使用 xarray 实现插值,本文介绍如何将该插值应用到 GRIB 2 文件中的要素场。
本文使用 nwpc-oper/nwpc-data 库从本地 GRIB 2 文件中加载要素场。nwpc-data 库介绍请浏览《nwpc-data 库简介》
准备
加载需要的库import numpy as npimport pandas as pdimport matplotlib.pyplot as pltimport xarray as xrfrom nwpc_data.data_finder import find_local_filefrom nwpc_data.grib.eccodes import load_field_from_file
站点插值
获取 850hPa 温度场data_file = find_local_file( "grapes_gfs_gmf/grib2/orig", start_time=pd.to_datetime("2020-04-26"), forecast_time="24h",)data_filePosixPath("/g1/COMMONDATA/OPER/NWPC/GRAPES_GFS_GMF/Prod-grib/2020042521/ORIG/gmf.gra.2020042600024.grb2")t850 = load_field_from_file( data_file, parameter="t", level_type="pl", level=500,)t850array([[262.26335937, 262.28335938, 262.29335938, ..., 262.26335937, 262.27335937, 262.30335938], ..., [243.63335937, 243.48335937, 243.64335938,