Question

我的文件夹名称的日期时间格式为 2018-08-21至2018-10-16

每个文件夹中都有一个zip文件，其中包含时间值（线性比例），时间线性增加。

我试图为每一天绘制图表，其中包含很多.bz2文件时间序列数据，即该日期的时间值。

现在我正在尝试这样做：

 timearr = np.asarray(data1['time'])
            ax.plot(np.asarray(timeStamps), timearr)

            ax.set_title('title')
            ax.set_ylabel('date vs time ')
            ax.grid(True)

            # Format the x-axis for dates (label formatting, rotation)
            fig.autofmt_xdate(rotation=45)
            fig.tight_layout()

            plt.show()

但是我收到一条错误消息，两个尺寸都不匹配。

dateStamps are list[2018-08-21
2018-08-22
2018-08-23
2018-08-24
2018-08-25]

data1['time'] = list of EPOC values.

Answer 1

不幸的是，我不知道文件的详细结构，所以我不得不猜测问题出在哪里

以下是一些代码，用于使用通用bz2文件生成某些文件夹：

import bz2
import numpy as np
import datetime
import os

startDate = datetime.datetime(2000,5,2,10,15,0,0)

for day in range(5):
    theDate = startDate +  datetime.timedelta(days=day)
    folder = "{}".format( theDate.replace( microsecond = 0 ).strftime("%Y-%m-%d") ) 
    os.mkdir( folder )
    data = ""
    for k in range(100):
        zzz = theDate + datetime.timedelta(seconds=137*k)
        data += "{} ".format( zzz.replace( microsecond = 0 ).strftime("%H:%M:%S") ) 
        d = zzz.day
        m = zzz.minute
        data += " {}\n".format( .17 * d + .003 * m**2 -.001 * m )
    myZip = bz2.BZ2File(os.path.join( folder, 'dat.bz2' ), 'w' )
    myZip.write( data )
    myZip.close()

这些文件夹和文件的处理方式为：

import bz2
import numpy as np
import datetime
import os
import matplotlib.pyplot as plt
import matplotlib.dates as mdates

"""
SE posts I used 
https://stackoverflow.com/questions/1574088/plotting-time-in-python-with-matplotlib
https://stackoverflow.com/questions/11264521/date-ticks-and-rotation-in-matplotlib
"""

def split_data( inData ):
    rows=data.strip().split('\n')
    rowcol = [x.split() for x in rows ]
    x,y = zip(*rowcol)
    y = [float(z) for z in  y ]
    x = [ datetime.datetime.strptime(z, '%H:%M:%S') for z in x]
    return x,y

dataDict = dict()
for root, dirs, files in os.walk("."):
   for name in files:
       if name.split('.')[-1]=='bz2':
           base = os.path.basename( root )
           myPath = (os.path.join(root, name))
           bz = bz2.BZ2File( myPath, 'r' )
           data = bz.read()
           dataDict[ base ] = split_data( data )

myFmt = mdates.DateFormatter('%H:%M')

fig = plt.figure()
ax = fig.add_subplot( 1, 1, 1 )

for key, dt in dataDict.iteritems():
    ax.plot( *dt , label=key )
ax.xaxis.set_major_formatter(myFmt)
for label in ax.get_xmajorticklabels():
    label.set_rotation(30)
ax.set_ylabel('data (arb. u.)')
ax.set_xlabel('time')
ax.legend( loc=0 )
plt.tight_layout()
plt.show()

提供：

希望我做对了。

从文件绘制日期与时间

1 个答案: