Question

我有一个大型数据集：20,000 x 40,000作为numpy数组。我把它保存为腌菜文件。

我不想将这个庞大的数据集读入内存，而是一次只读取几行（比方说100行），作为一个小批量使用。

如何从pickle文件中只读取一些随机选择的（无替换）行？

Answer 1

您可以逐步将腌菜写入文件，以便加载它们也是渐进的。

采用以下示例。在这里，我们迭代列表的项目，和每个人依次腌制。

>>> import cPickle
>>> myData = [1, 2, 3]
>>> f = open('mydata.pkl', 'wb')
>>> pickler = cPickle.Pickler(f)
>>> for e in myData:
...     pickler.dump(e)
<cPickle.Pickler object at 0x7f3849818f68>
<cPickle.Pickler object at 0x7f3849818f68>
<cPickle.Pickler object at 0x7f3849818f68>
>>> f.close()

现在我们可以反向执行相同的过程并根据需要加载每个对象。对于例如，让我们说我们只想要第一个项而不是想迭代整个文件。

>>> f = open('mydata.pkl', 'rb')
>>> unpickler = cPickle.Unpickler(f)
>>> unpickler.load()
1

此时，文件流只有第一个进展宾语。其余的对象未加载，这正是您的行为想。为了证明，您可以尝试阅读文件的其余部分，然后查看其余部分还坐在那里。

>>> f.read()
'I2\n.I3\n.'

Answer 2

由于您不知道pickle的内部工作原理，您需要使用另一种存储方法。下面的脚本使用tobytes()函数将数据按行保存在原始文件中。

由于每行的长度是已知的，因此可以通过seek()和read()计算和访问文件中的偏移量。之后，它将转换回具有frombuffer()函数的数组。

然而，最大的免责声明是数据的大小未保存（这也可以添加，但需要更多的复杂性），并且此方法可能不像腌制数组那样可移植。

正如@PadraicCunningham在comment中指出的那样，memmap可能是另一种优雅的解决方案。

评论性能：阅读评论后，我做了一个简短的基准测试。在我的机器上（16GB RAM，加密的SSD）我能够在24秒内完成40000个随机行读取（当然，使用20000x40000矩阵，而不是示例中的10x10）。

from __future__ import print_function
import numpy
import random

def dumparray(a, path):
    lines, _ = a.shape
    with open(path, 'wb') as fd:
        for i in range(lines):
            fd.write(a[i,...].tobytes())

class RandomLineAccess(object):
    def __init__(self, path, cols, dtype):
        self.dtype = dtype
        self.fd = open(path, 'rb')
        self.line_length = cols*dtype.itemsize

    def read_line(self, line):
        offset = line*self.line_length
        self.fd.seek(offset)
        data = self.fd.read(self.line_length)

        return numpy.frombuffer(data, self.dtype)

    def close(self):
        self.fd.close()


def main():
    lines = 10
    cols = 10
    path = '/tmp/array'

    a = numpy.zeros((lines, cols))
    dtype = a.dtype

    for i in range(lines):
        # add some data to distinguish lines
        numpy.ndarray.fill(a[i,...], i)

    dumparray(a, path)
    rla = RandomLineAccess(path, cols, dtype)

    line_indices = list(range(lines))
    for _ in range(20):
        line_index = random.choice(line_indices)
        print(line_index, rla.read_line(line_index))

if __name__ == '__main__':
    main()

Answer 3

谢谢大家。我最终找到了一种解决方法（一台具有更多RAM的机器，所以我实际上可以将数据集加载到内存中）。

如何从pickle文件一次加载一行？

3 个答案: