我试图获取关于我的图在哪一天截获阈值的信息(在本例中,某些hline y=-1)。我的一些结果是:
(11.805201846946158, -1.0)
(12.99215619183464, -1.0)
(21.33191612908414, -1.0)
(22.633241076423914, -1.0)
(37.00391584924397, -1.0)
我想换成11.8052。。。值一个特定的日期和时间。图形以日期时间绘制。我需要在脚本中添加什么才能及时打印它们
spi_avg_1['spi3'].plot(label='SPI 3- monthly')
plt.axhline(y=-1, color='r', linestyle='-', linewidth=0.5)
plt.title('SPI 3-monthly')
yys = [-1]
xx, yy = [], []
xo,yo = [k for k in range(200)],spi_avg_1['spi3']
for i in range(1,len(spi_avg_1['spi3'])):
for k in yys:
p1 = np.array([xo[i-1],yo[i-1]],dtype='float')
p2 = np.array([xo[i],yo[i]],dtype='float')
k1 = np.array([xo[i-1],k],dtype='float')
k2 = np.array([xo[i],k],dtype='float')
if crosses((p2,p1),(k1,k2)):
seg = line_intersection((p2,p1),(k1,k2))
print(seg)
数据集:
脚本如下:
相关问题 更多 >
编程相关推荐