python⽓象插值_Python⽓象数据处理进阶之Xarray（3）：插
值
Xarray库关于插值的操作内容不多，但是很实⽤，需要安装scipy库。
⾸先我们还是创建⼀个DataArray:
da = xr.DataArray(np.sin(0.3 * np.arange(12).reshape(4, 3)),[('time', np.arange(4)),('lat', [30, 60, 90])])
print(da)
#
#array([[ 0. , 0.29552021, 0.56464247],
# [ 0.78332691, 0.93203909, 0.99749499],
python蟒蛇编程# [ 0.97384763, 0.86320937, 0.67546318],
# [ 0.42737988, 0.14112001, -0.15774569]])
#Coordinates:
# * time (time) int64 0 1 2 3
vue三大核心
# * lat (lat) int64 30 60 90
exp1
我们先对时间维插值，⽐如说这⾥有第0天，1天，2天和3天，我们想获取第1.5天的数据：
print(da.interp(time=1.5))
#
#array([0.87858727, 0.89762423, 0.83647908])
#Coordinates:
# * lat (lat) int64 30 60 90
# time float64 1.5
socket连接是什么意思
可以看出实际上是对第1天和第2天的数据取了平均.
springcloud集成log4j2如果要往外插，⽐如说得到第3.5天的数据，得到的结果全是np.nan
也可以通过datetime格式的字符串指定插值时间：da.interp(time='2000-01-02')
也可以同时指定对多个维度插值da.interp(time=1.5，lat=45)
通过这个功能，就能实现不同⽹格的插值，甚⾄是站点数据到格点的插值。
exp2
这部分主要 是介绍⼏种插值的⽅法
import matplotlib.pyplot as plt
da = xr.DataArray(np.sin(np.linspace(0, 2 * np.pi, 10)), dims='x',coords={'x': np.linspace(0, 1, 10)})
da.plot.line('o', label='original')
da.interp(x=np.linspace(0, 1, 100)).plot.line(label='linear (default)')
da.interp(x=np.linspace(0, 1, 100), method='cubic').plot.line(label='cubic')
plt.legend()
插值效果对⽐
.interp()函数通过参数设置，可以实⽤不同的插值⽅法，⽐如linear线性，cubic⽴⽅等等
为了解决EXP1中提到的外插全是缺测的问题，可以通过如下参数的设置进⾏调整：
# ⽤0填充
da.interp(x=np.linspace(-0.5, 1.5, 10), kwargs={'fill_value': 0.0})
#
#array([ 0. , 0. , 0. , 0.814, 0.604, -0.604, -0.814, 0. , 0. , 0. ])
#Coordinates:
# * x (x) float64 -0.5 -0.2778 -0.05556 0.1667 ... 0.8333 1.056 1.278 1.5
# 外推
da.interp(x=np.linspace(-0.5, 1.5, 10), kwargs={'fill_value': 'extrapolate'})
#
#array([-2.893, -1.607, -0.321, 0.814, 0.604, -0.604, -0.814, 0.321, 1.607, 2.893]) #Coordinates:
# * x (x) float64 -0.5 -0.2778 -0.05556 0.1667 ... 0.8333 1.056 1.278 1.5
原本⼩于0和⼤于1的部分是没⽤数值的，但是通过这种⽅法，填充或外推出了对应的数值。exp3
接下来实现更⾼端⼀点的操作：
气象python零基础入门教程插值exp3
左图表⽰通过索引，选取出了⼀条斜线。
右图表⽰通过插值，选取出了⼀条不在格点上的斜线。
这种⽅法，在我们绘制不规则垂直剖⾯时⾮常需要。
da = xr.DataArray(np.sin(0.3 * np.arange(20).reshape(5, 4)),
[('x', np.arange(5)),('y', [0.1, 0.2, 0.3, 0.4])])
x = xr.DataArray([0, 2, 4], dims='z')
y = xr.DataArray([0.1, 0.2, 0.3], dims='z')
da.sel(x=x, y=y)
print(da.sel(x=x, y=y))
#
#array([ 0. , 0.42737988, -0.77276449])
#Coordinates:
# x (z) int64 0 2 4
# y (z) float64 0.1 0.2 0.3
#Dimensions without coordinates: z
以上代码实现了左图的操作。
da = xr.DataArray(np.sin(0.3 * np.arange(20).reshape(5, 4)),
[('x', np.arange(5)),('y', [0.1, 0.2, 0.3, 0.4])])
x = xr.DataArray([0.5, 1.5, 2.5], dims='z')
y = xr.DataArray([0.15, 0.25, 0.35], dims='z')
print(da.interp(x=x, y=y))
#
#array([ 0.55626357, 0.63496063, -0.46643289])
#Coordinates:
# x (z) float64 0.5 1.5 2.5
# y (z) float64 0.15 0.25 0.35
#Dimensions without coordinates: z
以上代码实现了右图的操作。
Metpy库中有类似操作的实现，本质就是对这个函数的封装。
exp4
接下来部分涉及到缺测值的问题。
da = xr.DataArray([0, 2, np.nan, 3, 3.25], dims='x',coords={'x': range(5)}) print(da.interp(x=[0.5, 1.5, 2.5]))
#
#array([ 1., nan, nan])
php货币80k是多少钱
#Coordinates:
# * x (x) float64 0.5 1.5 2.5
可以看到默认设置是⽆法处理缺测值的，可以⾸先去除缺测部分：
dropped = da.dropna('x')
print(dropped)
#
#array([0. , 2. , 3. , 3.25])
#Coordinates:
# * x (x) int64 0 1 3 4
print(dropped.interp(x=[0.5, 1.5, 2.5], method='cubic'))
#
#array([1.19010417, 2.5078125 , 2.9296875 ])
#Coordinates:
# * x (x) float64 0.5 1.5 2.5
但是这样之后数组的shape发⽣了变化⼗分不⽅便。XARRAY提供了填补缺测数据的⽅法：
filled = da.interpolate_na(dim='x')
print(filled)
#
#array([0. , 2. , 2.5 , 3. , 3.25])
#Coordinates:
# * x (x) int64 0 1 2 3 4
在填补之后，在进⾏插值便可。
总结
Xarray关于插值的⽤法还是⾮常好⽤的，可以实现包括：
1、填补缺测
2、站点插格点
3、任意剖⾯截取
4、不同分辨率格点互插
5、规则⽹格插⾮规则⽹格
6、数据延长
等等操作，这边的例⼦仅是参考官⽅⽂档的⼀些基础举例。

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。