start_index end_index在plot绘制中的作用

时间: 2024-02-24 08:54:36 浏览: 21
在绘制图形时,`start_index` 和 `end_index` 通常用于限制绘制的数据范围,从而实现对图形的缩放和平移。例如,在绘制折线图时,可以通过限制 x 轴的数据范围来实现横向缩放和平移。 以下是一个使用 Matplotlib 绘制折线图并使用 `start_index` 和 `end_index` 限制 x 轴数据范围的示例代码: ```python import tkinter as tk from matplotlib.backends.backend_tkagg import FigureCanvasTkAgg from matplotlib.figure import Figure # 创建一个顶级窗口对象 root = tk.Tk() root.title("Plot Canvas Example") # 创建一个 Matplotlib 图形对象 fig = Figure(figsize=(5, 4), dpi=100) ax = fig.add_subplot(111) ax.set_xlabel('X Label') ax.set_ylabel('Y Label') # 创建数据 x = [1, 2, 3, 4, 5] y = [1, 4, 9, 16, 25] # 绘制折线图 line, = ax.plot(x, y, 'ro-') # 创建一个 Matplotlib 画布对象 canvas = FigureCanvasTkAgg(fig, master=root) canvas.draw() canvas.get_tk_widget().pack(side=tk.TOP, fill=tk.BOTH, expand=1) # 定义缩放和平移函数 def on_zoom_in(): global start_index, end_index start_index = max(0, start_index-1) end_index = min(len(x)-1, end_index+1) ax.set_xlim(x[start_index], x[end_index]) canvas.draw() def on_zoom_out(): global start_index, end_index start_index = min(len(x)-1, start_index+1) end_index = max(0, end_index-1) ax.set_xlim(x[start_index], x[end_index]) canvas.draw() # 初始化数据范围 start_index = 0 end_index = len(x) - 1 ax.set_xlim(x[start_index], x[end_index]) # 创建两个按钮用于缩放和平移 zoom_in_button = tk.Button(root, text="Zoom In", command=on_zoom_in) zoom_out_button = tk.Button(root, text="Zoom Out", command=on_zoom_out) zoom_in_button.pack(side=tk.LEFT) zoom_out_button.pack(side=tk.LEFT) # 进入事件循环 root.mainloop() ``` 这段代码创建了两个按钮 `zoom_in_button` 和 `zoom_out_button`,用于实现横向缩放和平移。`start_index` 和 `end_index` 变量用于记录当前 x 轴的数据范围,并在缩放和平移操作中更新。`ax.set_xlim()` 方法用于设置 x 轴数据范围,并在每次缩放和平移操作后重新绘制图形。

相关推荐

import pandas as pd import matplotlib.pyplot as plt import numpy as np import pywt file_name = 'E:/liuyuan/ceshi/zhongyao/Subject_1_0cmH20_norm_breaths.csv' data = pd.read_csv(file_name, skiprows=1, usecols=[0, 2], names=['Time', 'Flow']) x = list() y = list() for i in range(len(data)): x.append(float(data.values[i][0])) y.append(float(data.values[i][1])) start_index = 0 end_index = 5372 time = np.arange(start_index, end_index) flow = np.arange(start_index, end_index) time = data['Time'][start_index:end_index] flow = data['Flow'] def wavelet_filter(data): wavelet = 'db4' # 选择小波基函数 level = 5 # 小波变换的层数 # 小波变换 coeffs = pywt.wavedec(data, wavelet, level=level) threshold = np.std(coeffs[-level]) * np.sqrt(2 * np.log(len(data))) coeffs[1:] = (pywt.threshold(c, threshold, mode='soft') for c in coeffs[1:]) filtered_data = pywt.waverec(coeffs, wavelet) return filtered_data 对Flow进行小波变换滤波 filtered_flow = wavelet_filter(flow) fig, ax = plt.subplots(figsize=(10, 5)) plt.xlim(0, 60) ax.set_ylim(-0.7, 0.7) ax.set_xlabel('Time(s)', fontsize=10) ax.set_ylabel('Flow(L/s)', fontsize=10) ax.plot(time, filtered_flow, label='Filtered Flow') ax.legend() ax.grid(True, linewidth=0.3, alpha=0.5, color='gray') plt.tight_layout() # 自动调整子图的布局 plt.show()import pandas as pd import matplotlib.pyplot as plt import numpy as np import pywt file_name = 'E:/liuyuan/ceshi/zhongyao/Subject_1_0cmH20_norm_breaths.csv' data = pd.read_csv(file_name, skiprows=1, usecols=[0, 2], names=['Time', 'Flow']) x = list() y = list() for i in range(len(data)): x.append(float(data.values[i][0])) y.append(float(data.values[i][1])) start_index = 0 end_index = 5372 time = np.arange(start_index, end_index) flow = np.arange(start_index, end_index) time = data['Time'][start_index:end_index] flow = data['Flow'] def wavelet_filter(data): wavelet = 'db4' # 选择小波基函数 level = 5 # 小波变换的层数 coeffs = pywt.wavedec(data, wavelet, level=level) threshold = np.std(coeffs[-level]) * np.sqrt(2 * np.log(len(data))) coeffs[1:] = (pywt.threshold(c, threshold, mode='soft') for c in coeffs[1:]) # 逆小波变换 filtered_data = pywt.waverec(coeffs, wavelet) return filtered_data filtered_flow = wavelet_filter(flow) fig, ax = plt.subplots(figsize=(10, 5)) plt.xlim(0, 60) ax.set_ylim(-0.7, 0.7) ax.set_xlabel('Time(s)', fontsize=10) ax.set_ylabel('Flow(L/s)', fontsize=10) ax.plot(time, filtered_flow, label='Filtered Flow') ax.legend() ax.grid(True, linewidth=0.3, alpha=0.5, color='gray') plt.tight_layout() plt.show()在此代码上添加切分代码,并按照合适窗口大小切分完准确显示

import numpy as np import matplotlib.pyplot as plt from matplotlib.ticker import MaxNLocator # 创建画布和子图对象 fig, ax = plt.subplots(figsize=(9, 6), dpi=100) # 绘制折线图 ax.plot(x, y) # 绘制平均值线 #ax.axhline(y=-650, color='r', linestyle='--',label='流域整体物质平衡=-650mm w.e.') # 添加阴影带 start_year = 2006 end_year = 2016 mask = np.logical_and(years >= start_year, years <= end_year) years_to_plot = years[mask] ax.fill_between(years_to_plot, -680- 220, -680 + 220, alpha=0.2,color='yellow',label='Brun et al.2017') ax.axhline(-680, color='yellow', linestyle='--',xmin=0.65, xmax=0.89) start_year_2 = 2000 end_year_2 = 2014 mask_2 = np.logical_and(years >= start_year_2, years <= end_year_2) years_to_plot_2 = years[mask_2] ax.fill_between(years_to_plot_2, -790-110, -790+110, alpha=0.2, color='green',label='Wu et al.2018') ax.axhline(-790, color='green', linestyle='--',xmin=0.51, xmax=0.840) start_year_3 = 2000 end_year_3 = 2018 mask_3 = np.logical_and(years >= start_year_3, years <= end_year_3) years_to_plot_3 = years[mask_3] ax.fill_between(years_to_plot_3, -540-160, -540+160, alpha=0.2, color='blue',label='Shean et al.2020') ax.axhline(-540, color='blue', linestyle='--',xmin=0.51, xmax=0.93) start_year_4 = 2000 end_year_4 = 2019 mask_4 = np.logical_and(years >= start_year_4, years <= end_year_4) years_to_plot_4 = years[mask_4] ax.fill_between(years_to_plot_4, -580-220, -580+220, alpha=0.2, color='red',label='Hugonnet et al.2021') ax.axhline(-580, color='red', linestyle='--',xmin=0.51, xmax=0.957) # 设置 x 轴标签和标题 ax.set_xlabel('年份',fontproperties=font_prop,fontsize=14) ax.set_ylabel('物质平衡(mm w.e.)',fontproperties=font_prop,fontsize=14) ax.set_title('图8 帕隆藏布流域1980-2019物质平衡',fontproperties=font_prop,fontsize=14,y=-0.17) # 强制显示整数刻度 ax.xaxis.set_major_locator(MaxNLocator(integer=True)) # 添加网格 ax.grid(True, which='major', linestyle='--') # 将坐标轴的刻度字体大小设置为12 ax.tick_params(axis='both', which='major', labelsize=12) # 添加图例 ax.legend(fontsize=24,loc='lower left',prop=font_prop) # 设置图形的边距 plt.tight_layout() # 显示图形 plt.show()

最新推荐

recommend-type

2280.宁乡杨氏绍纶谱: 十卷.pdf

2280.宁乡杨氏绍纶谱: 十卷
recommend-type

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节 线性代数的核心问题是求解方程组。这些方程都是线性的,即未知数仅与数相乘——我们绝不会 遇见 x 乘以 y。我们的第一个线性方程组较小。接下来你来看看它引申出多远: 两个方程 两个未知数 x − 2y = 1 3x + 2y = 11 (1) 我们一次从一个行开始。第一个方程 x − 2y = 1 得出了 xy 平面的一条直线。由于点 x = 1, y = 0 解 出该方程,因此它在这条直线上。因为 3 − 2 = 1,所以点 x = 3, y = 1 也在这条直线上。若我们选择 x = 101,那我们求出 y = 50。 这条特定直线的斜率是 12,是因为当 x 变化 2 时 y 增加 1。斜率在微积分中很重要,然而这是线 性代数! 图 2.1 将展示第一条直线 x − 2y = 1。此“行图”中的第二条直线来自第二个方程 3x + 2y = 11。你 不能错过两条线的交点 x = 3, y = 1。点 (3, 1) 位于两条线上并且解出两个方程。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

识别MATLAB微分方程求解中的混沌行为:分析非线性方程混沌行为的实用技巧

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程是描述物理、工程和金融等领域中动态系统的数学模型。MATLAB提供了强大的工具来求解微分方程,包括内置求解器和自定义函数
recommend-type

physon如何做ERP系统

作为一个AI语言模型,我可以给您提供一些基本的建议和思路: 1.明确需求和目标:在开始开发ERP系统之前,需要明确您的业务需求和目标,了解您的业务流程、流程数据、员工角色和权限、财务管理等方面的要求,这样才能更好的设计和开发ERP系统。 2.选择技术框架和开发工具:选择合适的技术框架和开发工具是ERP系统开发的关键。选择一种流行的技术框架和工具可以提高开发效率和质量。 3.设计数据库:ERP系统需要一个功能强大的数据库来存储数据。设计数据库需要考虑数据的完整性、安全性和可扩展性。 4.设计系统架构:系统架构是ERP系统的骨架,需要考虑系统的可扩展性、可维护性和性能。 5.开发和测试:
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

探索MATLAB微分方程求解中的分岔分析:揭示方程动态行为的秘密

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程在科学、工程和金融等领域有着广泛的应用。MATLAB作为一种强大的数值计算软件,提供了丰富的微分方程求解工具。本章将概述
recommend-type

ic验证工作中如何在平台中加入发数的总数?

在进行IC验证工作时,可以通过以下步骤在平台中加入发数的总数: 1. 打开IC验证工具(如Cadence Virtuoso)并打开对应的设计文件。 2. 在设计文件中选择需要计算发数的部分电路或模块。 3. 在IC验证工具中打开时序分析工具(如Cadence Tempus)。 4. 在时序分析工具中设置好时钟频率、时钟周期等参数。 5. 执行时序分析,生成时序报告。 6. 在时序报告中查找发数统计信息,将其记录下来。 7. 将发数统计信息添加到平台中,以便在之后的仿真或验证中使用。 需要注意的是,发数统计信息可能因为设计文件的不同而有所差异,需要根据实际情况进行调整和计算。
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。