查找横线在各自的Y轴值上穿过的X轴数据点。

Question

我在x-axis上有一个x-axis数据点，在y-axis (比如[10,15,10,10,20] )上有相应的值。

通常，为了通过给定的y-axis数据点(如y=f(x) )找到x-axis的值点，我检查了它，我们可以通过interpolation使用numpy实现这一点。但我没有找到如何用给定的y轴值插值x轴.根据附加的屏幕，我想找到分别的x轴值，其中线12交叉..。

所以我期望得到类似于[1, 1.x, 2, 2.x, 3, 4, 4.x, 5, 5.x] on x-axis的结果

​

​

Answer 1

如果这是一条平滑的曲线，您可以使用InterpolatedUnivariateSpline

import numpy as np
from scipy import interpolate

x = np.linspace(0, 20, 100)
y = np.sin(x + 0.1)

y0 = 0.3
spline = interpolate.InterpolatedUnivariateSpline(x, y - y0)

xp = spline.roots()

情节如下：

pl.plot(x, y)
pl.axhline(0.3, color="black", linestyle="dashed")
pl.vlines(xp, 0, 0.3, color="gray", linestyle="dotted")

​

​

如果您想要线性内插：

x = np.linspace(0, 20, 20)
y = np.sin(x + 0.1)
y0 = 0.3
y_offset = y - y0
pos = np.where((y_offset[1:] * y_offset[:-1]) <= 0)[0]

x1 = x[pos]
x2 = x[pos+1]
y1 = y[pos]
y2 = y[pos+1]

xp = (y0 - y1) / (y2 - y1) * (x2 - x1) + x1

​

​

Answer 2

如果将interp1d(x,y)更改为interp1d(y,x)，则将x表示为y的函数。

请注意，如果f(x)不是唯一的，则可能会出现意外或未定义的行为。