下面的代码应该可以解决问题。首先,它打开一个文件并在lzma 中对其进行解码,然后使用struct 解压缩二进制数据。
import lzma
import struct
import pandas as pd
def bi5_to_df(filename, fmt):
chunk_size = struct.calcsize(fmt)
data = []
with lzma.open(filename) as f:
while True:
chunk = f.read(chunk_size)
if chunk:
data.append(struct.unpack(fmt, chunk))
else:
break
df = pd.DataFrame(data)
return df
最重要的是要知道正确的格式。我四处搜索并试图猜测,'>3i2f'(或>3I2f)工作得很好。 (它是大端 3 个整数 2 个浮点数。您的建议是:'i4f' 不会产生合理的浮点数 - 无论是大端还是小端。)对于 struct 和格式语法,请参阅 docs。
df = bi5_to_df('13h_ticks.bi5', '>3i2f')
df.head()
Out[177]:
0 1 2 3 4
0 210 110218 110216 1.87 1.12
1 362 110219 110216 1.00 5.85
2 875 110220 110217 1.00 1.12
3 1408 110220 110218 1.50 1.00
4 1884 110221 110219 3.94 1.00
更新
比较bi5_to_df 和https://github.com/ninety47/dukascopy 的输出,
我从那里编译并运行test_read_bi5。输出的第一行是:
time, bid, bid_vol, ask, ask_vol
2012-Dec-03 01:00:03.581000, 131.945, 1.5, 131.966, 1.5
2012-Dec-03 01:00:05.142000, 131.943, 1.5, 131.964, 1.5
2012-Dec-03 01:00:05.202000, 131.943, 1.5, 131.964, 2.25
2012-Dec-03 01:00:05.321000, 131.944, 1.5, 131.964, 1.5
2012-Dec-03 01:00:05.441000, 131.944, 1.5, 131.964, 1.5
和bi5_to_df 在同一个输入文件上给出:
bi5_to_df('01h_ticks.bi5', '>3I2f').head()
Out[295]:
0 1 2 3 4
0 3581 131966 131945 1.50 1.5
1 5142 131964 131943 1.50 1.5
2 5202 131964 131943 2.25 1.5
3 5321 131964 131944 1.50 1.5
4 5441 131964 131944 1.50 1.5
所以一切似乎都很好(ninety47 的代码重新排列了列)。
另外,使用'>3I2f' 代替'>3i2f' 可能更准确(即unsigned int 代替int)。