我有一个Python程序,显示了温度随时间下降的曲线图。沿下降,温度在一段时间内保持恒定,几乎为0的斜率,然后继续下降。我希望程序在温度恒定时在曲线上的这个区域自动检测并显示y值。稍后将这个值放到方程中。我正在尝试找出方法。我尝试过并失败了,我的最后一次尝试是:
import numpy as np
import matplotlib.pyplot as plt
list_of_files=[('logfile.txt', 'temp')]
datalist = [ ( np.loadtxt(filename), label ) for filename, label in list_of_files]
for data, label in datalist:
plt.plot( data[:0], data[:,1], label=label )
plt.ginput(n=1, timeout=30, show_clicks=True, mouse_add=1, mouse_pops=3, mouse_stop=2)
plt.show()
我希望在高原上单击鼠标会显示并保存y坐标,只是将我引向正确的编程方向。但是,当我单击该图时,所有这些只是一个简短的红色标记。我不想点击鼠标。。谢谢Rico。
参考方案
遍历数据的小块,确定块的斜率,返回符合条件的点
def zero_slope(data, chunksize = 3, max_slope = .001):
"""return the 'first' data point with zero slope
data --> numpy ndarray - 2d [[x0,y0],[x1,y1],...]
chunksize --> odd int
returns numpy ndarray
"""
midindex = chunksize / 2
for index in xrange(len(data) - chunksize):
chunk = data[index : index + chunksize, :]
# subtract the endpoints of the chunk
# if not sufficient, maybe use a linear fit
dx, dy = abs(chunk[0] - chunk[-1])
print dy, dx, dy / dx
if 0 <= dy / dx < max_slope:
return chunk[midindex]
Python sqlite3数据库已锁定 - python我在Windows上使用Python 3和sqlite3。我正在开发一个使用数据库存储联系人的小型应用程序。我注意到,如果应用程序被强制关闭(通过错误或通过任务管理器结束),则会收到sqlite3错误(sqlite3.OperationalError:数据库已锁定)。我想这是因为在应用程序关闭之前,我没有正确关闭数据库连接。我已经试过了: connectio…
Matplotlib-固定x轴缩放比例和自动缩放y轴 - python我只想绘制部分数组,固定x部分,但让y部分自动缩放。我尝试如下所示,但是它不起作用。有什么建议么?import numpy as np import matplotlib.pyplot as plt data=[np.arange(0,101,1),300-0.1*np.arange(0,101,1)] plt.figure() plt.scatter(da…
python-docx应该在空单元格已满时返回空单元格 - python我试图遍历文档中的所有表并从中提取文本。作为中间步骤,我只是尝试将文本打印到控制台。我在类似的帖子中已经看过scanny提供的其他代码,但是由于某种原因,它并没有提供我正在解析的文档的预期输出可以在https://www.ontario.ca/laws/regulation/140300中找到该文档from docx import Document from…
Python:集群作业管理 - python我在具有两个阶段的计算群集(Slurm)上运行python脚本,它们是顺序的。我编写了两个python脚本,一个用于阶段1,另一个用于阶段2。每天早上,我检查所有第1阶段的工作是否都以视觉方式完成。只有这样,我才开始第二阶段。通过在单个python脚本中组合所有阶段和作业管理,是否有一种更优雅/自动化的方法?我如何知道工作是否完成?工作流程类似于以下内容:w…
Python-Excel导出 - python我有以下代码:import pandas as pd import requests from bs4 import BeautifulSoup res = requests.get("https://www.bankier.pl/gielda/notowania/akcje") soup = BeautifulSoup(res.cont…