降雨-径流过程图

import pandas as pd
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
# 读取Excel文件
df = pd.read_excel(r'C:\Users\Admin\Desktop\东湾.xlsx', sheet_name='测试')
# 将日期列转换为日期格式
df['日期'] = pd.to_datetime(df['日期'])
# 提取实测流量和模拟流量数据
true = df['流量']
sim = df['测试3']
pre = df['雨量']
# 创建一个新的图形和轴对象
fig, ax1 = plt.subplots(figsize=(10, 6))
plt.rcParams['font.sans-serif'] = 'SimHei'  # 设置中文显示
# 绘制实测流量折线
ax1.plot(df['日期'], true, label='实测流量', marker='.', linestyle='--', color='blue')
# 绘制模拟流量折线
ax1.plot(df['日期'], sim, label='模拟流量', color='orange')
# 设置第一个y轴标签和颜色
ax1.set_ylabel('流量/($m^3·s^{-1}$)', color='black')
ax1.tick_params(axis='y', labelcolor='black', direction='in')
# 设置主y轴范围
ax1.set_ylim(0, 500)
# 创建第二个y轴对象，共享x轴
ax2 = ax1.twinx()
# 绘制雨量柱状图
ax2.bar(df['日期'], pre, label='雨量', alpha=0.5, color='blue', width=2)
# 设置第二个y轴标签和颜色
ax2.set_ylabel('雨量/(mm)', color='black')
ax2.tick_params(axis='y', labelcolor='black', direction='in')
# 设置次y轴范围并倒置
ax2.set_ylim(100, 0)
# 设置x轴日期刻度为1年为间隔
ax1.xaxis.set_major_locator(mdates.YearLocator(1))
# 设置x轴刻度线向内
ax1.tick_params(axis='x', direction='in')
# 合并图例
lines, labels = ax1.get_legend_handles_labels()
lines2, labels2 = ax2.get_legend_handles_labels()
ax2.legend(lines + lines2, labels + labels2, loc='upper center', bbox_to_anchor=(0.5, -0.06), ncol=3, frameon=False)  # 调整bbox_to_anchor参数

# 设置整个图形的xlabel
fig.text(0.5, 0.06, '年份', ha='center', va='center')
plt.title('降雨-径流图')
# 显示图形
plt.show()