使用struct.unpack解压缩二进制文件VS np.frombuffer VS np.ndarray VS np.fromfile

时间:2019-02-13 21:46:17

标签: python numpy struct

我正在解压缩具有许多不同数据类型的大型二进制文件(〜1GB)。我正处于创建隐蔽每个字节的循环的早期阶段。 我一直在使用struct.unpack,但最近认为如果我使用numpy,它将运行得更快。但是,切换到numpy减慢了我的程序的速度。 我尝试过:

struct.unpack
np.fromfile
np.frombuffer
np.ndarray

注意:在np.fromfile方法中,我将文件保持打开状态,并且不将其加载到内存中并进行查找

1)

with open(file="file_loc" , mode='rb') as file: 
    RAW = file.read()
byte=0
len = len(RAW)
while( byte < len):
    header = struct.unpack(">HHIH", RAW[byte:(byte+10)])
    size = header[1]
    loc  = str(header[3])
    data[loc] = struct.unpack(">B", RAW[byte+10:byte+size-10)
    byte+=size

2)

dt=('>u2,>u2,>u4,>u2')
with open(file="file_loc" , mode='rb') as RAW:
    same loop as above:
        header = np.fromfile(RAW[byte:byte+10], dtype=dt, count=1)[0]
        data   = np.fromfile(RAW[byte+10:byte+size-10], dtype=">u1", count=size-10)

3)

dt=('>u2,>u2,>u4,>u2')
with open(file="file_loc" , mode='rb') as file:
    RAW = file.read()
same loop:
    header = np.ndarray(buffer=RAW[byte:byte+10], dtype=dt_header, shape= 1)[0]
    data   = np.ndarray(buffer=RAW[byte+10:byte+size-10], dtype=">u1", shape=size-10)

4)pretty much the same as 3 except using np.frombuffer()

所有numpy实现的速度大约是struct.unpack方法的一半,这不是我期望的。

让我知道我有什么办法可以改善性能。

此外,我只是从内存中键入了此内容,因此可能会有一些错误。

1 个答案:

答案 0 :(得分:0)

我并没有使用struct,但是在您的代码和文档之间,我可以将其用于存储整数数组的缓冲区。

numpy数组中获取字节数组/字符串。

In [81]: arr = np.arange(1000)
In [82]: barr = arr.tobytes()
In [83]: type(barr)
Out[83]: bytes
In [84]: len(barr)
Out[84]: 8000

相反的是tobytes

In [85]: x = np.frombuffer(barr, dtype=int)
In [86]: x[:10]
Out[86]: array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
In [87]: np.allclose(x,arr)
Out[87]: True

ndarray也可以使用,尽管通常不建议直接使用此构造函数:

In [88]: x = np.ndarray(buffer=barr, dtype=int, shape=(1000,))
In [89]: np.allclose(x,arr)
Out[89]: True

要使用struct,我需要创建一种包含长度“ 1000 long”的格式:

In [90]: tup = struct.unpack('1000l', barr)
In [91]: len(tup)
Out[91]: 1000
In [92]: tup[:10]
Out[92]: (0, 1, 2, 3, 4, 5, 6, 7, 8, 9)
In [93]: np.allclose(np.array(tup),arr)
Out[93]: True

因此,既然我们已经建立了读取缓冲区的等效方法,请执行一些计时:

In [94]: timeit x = np.frombuffer(barr, dtype=int)
617 ns ± 0.806 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)

In [95]: timeit x = np.ndarray(buffer=barr, dtype=int, shape=(1000,))
1.11 µs ± 1.76 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)

In [96]: timeit tup = struct.unpack('1000l', barr)
19 µs ± 38.4 ns per loop (mean ± std. dev. of 7 runs, 100000 loops each)

In [97]: timeit tup = np.array(struct.unpack('1000l', barr))
87.5 µs ± 25.1 ns per loop (mean ± std. dev. of 7 runs, 10000 loops each)

frombuffer看起来不错。

您的struct.unpack循环使我感到困惑。我认为它的作用与frombuffer相同。但是就像开始时所说的,我并没有使用struct