如何在Python中为NetCDF时间序列图X轴显示年份（850-1849）

我懂你的烦恼——用这段代码绘制千年跨度的时间序列时，X轴的年份显示完全不是你想要的样子对吧？问题出在matplotlib默认的日期刻度逻辑上，它没法自动适配这么长的时间跨度，咱们手动调整几下就能搞定！

问题根源 #

你用datetime生成的日期序列是正确的，但matplotlib默认会根据数据范围自动选择刻度间隔，对于850-1849年这种千年跨度，默认的刻度要么挤成一团看不清，要么显示的不是纯年份格式，自然不符合需求。

修改后的完整代码 #

import netCDF4 as netcdf
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
import datetime

# 读取NetCDF数据
nc = netcdf.Dataset('tas_850-1849.nc')
time = nc.variables['time']
lat = nc.variables['lat'][:]
lon = nc.variables['lon'][:]
var = nc.variables['pdo'][:]

# 生成日期序列（这部分你原来的代码没问题）
start = datetime.date(850, 1, 1)
dates = [start + datetime.timedelta(days=n) for n in range(len(var))]

# 绘制折线图
plt.figure(figsize=(12, 6))  # 可选：调整图的大小，让X轴显示更宽松
plt.plot(dates, var[:,0,0], color='red', linewidth=2, label='Data1')

# 关键：自定义X轴刻度与格式
# 每隔100年显示一个主刻度（可按需改成50、200等）
years = mdates.YearLocator(100)
# 设置刻度显示为4位纯年份格式
year_fmt = mdates.DateFormatter('%Y')

ax = plt.gca()
ax.xaxis.set_major_locator(years)
ax.xaxis.set_major_formatter(year_fmt)

# 旋转X轴标签，避免文字重叠
plt.xticks(rotation=45)

# 锁定X轴范围，确保严格覆盖850-1849年
plt.xlim(datetime.date(850, 1, 1), datetime.date(1849, 12, 31))

# 添加图例与标签，让图表更清晰
plt.legend()
plt.xlabel('年份')
plt.ylabel('PDO指数')
plt.title('850-1849年PDO时间序列')

plt.tight_layout()  # 自动调整布局，防止标签被截断
plt.show()

关键调整点说明 #

• YearLocator(100)：指定每隔100年放置一个主刻度，你可以根据需求调整间隔（比如改成50，X轴会显示更多年份刻度）。
• DateFormatter('%Y')：强制刻度显示为4位年份，彻底避免matplotlib自动显示其他冗余日期格式。
• plt.xticks(rotation=45)：旋转标签角度，解决年份文字重叠的问题。
• plt.xlim()：锁定X轴起止范围，确保图表只展示850到1849年的区间，不会出现多余的刻度。

这样修改后，你的X轴就会清晰显示你想要的年份刻度啦！

内容的提问来源于stack exchange，提问作者S.Br