使用Python

Question

我试图从一个NetCDF文件中的几个1-D变量绘制时间序列。 NetCDF文件输出类似于下面的内容，变量的名称与var [1..4]不同。

dimensions:
        time = UNLIMITED ; // (12 currently)
variables:
        float var1(time) ;
        float var2(time) ;
        float var3(time) ;
        float var4(time) ;
        double time(time) ;

以下情况适用于绘制一个变量。

import netCDF4
import pandas as pd
import matplotlib.pyplot as plt

# Data file
data = netCDF4.Dataset('./ts.nc', 'r')

# Print variables
par_list = data.variables.keys()
print par_list

# Read the variables
times = data.variables['time']
var1 = data.variables['var1']

# Convert times
t = netCDF4.num2date(times[:],times.units)

hs = pd.Series(var1[:],index=dtime)
fig = plt.figure(figsize=(12,4))
ax = fig.add_subplot(111)
hs.plot(ax=ax,title='Test plot')
ax.set_ylabel(var1.units)
plt.savefig("test.png")

除了绘制多行之外，我还想格式化这些行。变量1与变量2相关，变量3与变量4相关。因此变量1和2具有相同的线型和不同的颜色，变量3和4的不同线型将是有用的。另外，我正在寻找情节之外的传奇。