metpy利用cross_section插值一半是NAN怎么解决

liangbb520c

sp_time = datetime(2022,7,11,00)#sp时间
lon_range = [102, 105]  # 设置经度范围
lat_range = [33, 29]  # 设置纬度范围
spfnm = r'F:/Sichuan Meteorological Bureau work/era5/2d_11_download.nc'# need only one times data
spf = xr.open_dataset(spfnm).sel(latitude = slice(lat_range[0], lat_range[1]),
longitude = slice(lon_range[0], lon_range[1]),time=sp_time,level=850)
sp_f = spf.metpy.parse_cf().squeeze()
start = (34, 105)
end = (32, 107)
cross_sp = cross_section(sp_f, start, end,steps=50)
print(cross_sp['u'])
<xarray.DataArray 'u' (index: 50)>
array([-1.07999706, -1.07790156, -1.06590956, -1.04402936, -1.01226928,
       -0.97063761, -0.91914262, -0.82433094, -0.72911539, -0.63190238,
       -0.5326935 , -0.43149031, -0.32829438, -0.24774613, -0.16921265,
       -0.08941462, -0.00835303,  0.07397116,  0.15755697,  0.09138544,
       -0.0897057 , -0.30673002, -0.55965731, -0.84845745, -1.17310034,
               nan,         nan,         nan,         nan,         nan,
               nan,         nan,         nan,         nan,         nan,
               nan,         nan,         nan,         nan,         nan,
               nan,         nan,         nan,         nan,         nan,
               nan,         nan,         nan,         nan,         nan])

小鱼糯糯

请问你解决了吗？我也遇到了这个问题，南边数据全为nan