更好的 python 减少内存使用的代码？

Question

我有一个大约 1900 万行的数据框，其中 4 个变量是纬度和经度。我在 python haversine 包的帮助下创建了一个计算纬度和经度距离的函数。

# function to calculate distance of 2 coordinates

def measure_distance(lat_1, long_1, lat_2, long_2):

    coordinate_start = list(zip(lat_1, long_1))
    coodrinate_end = list(zip(lat_2, long_2))
    
    distance = haversine_vector(coordinate_start, coodrinate_end, Unit.KILOMETERS)

    return distance

我使用魔术命令 %%memit 来测量内存使用情况以执行计算。平均而言，内存使用量在 8 - 10 GB 之间。我运行我在 Google Colab 上工作，它有 12GB RAM，因此，有时操作会达到运行时间的限制并重新启动。

%%memit

measure_distance(df.station_latitude_start.values, 
                 df.station_longitude_start.values, 
                 df.station_latitude_end.values, 
                 df.station_longitude_end.values)

peak memory: 7981.16 MiB, increment: 5312.66 MiB

有没有办法优化我的代码？

Answer 1

TL;DR: 使用 Numpy 并通过 chunk.

对于大输入大小，CPython 解释器占用的内存量是预期的。

事实上，CPython 使用 references 将值存储在列表中。在 64 位系统上，引用占用 8 个字节，基本类型（浮点数和小整数）通常占用 32 个字节。两个浮点数的元组是一个复杂类型，它包含元组的大小以及两个浮点数的 references（不是值本身）。它的大小应该接近 64 字节。由于您有 2 个列表包含 1900 万个（引用）浮点数对和 4 个列表包含 1900 万个（引用）浮点数，因此占用的内存应该约为 4*19e6*(8+32) + 2*19e6*(8+64) = 5.7 GB。更不用说 Haversine 可以制作一些内部副本，结果也需要一些 space。

如果您想减少内存使用，那么使用 Numpy。事实上，float Numpy 数组以更紧凑的方式存储值（没有引用，没有内部标签）。您可以用 N x 2 Numpy 二维数组替换元组列表。生成的大小应该约为 4*19e6*8 + 2*19e6*(8*2) = 1.2 GB。此外，Haversine 在内部使用 Numpy 的计算速度会快得多。这是一个例子：

import numpy as np
 
# Assume lat_1, long_1, lat_2 and long_2 are of type np.array.
# Use np.array(yourList) if you want to convert it.
def measure_distance(lat_1, long_1, lat_2, long_2):
    coordinate_start = np.column_stack((lat_1, long_1))
    coordinate_end = np.column_stack((lat_2, long_2))
    return haversine_vector(coordinate_start, coordinate_end, Unit.KILOMETERS)

上面的代码大约快 25 倍。

如果您想进一步减少内存使用，可以通过块计算坐标（例如 32K 值），然后连接输出块。如果您不太关心计算距离的准确性，您也可以使用 单精度 数字而不是双精度。

这是一个如何按块计算结果的示例：

def better_measure_distance(lat_1, long_1, lat_2, long_2):
    chunckSize = 65536
    result = np.zeros(len(lat_1))
    for i in range(0, len(lat_1), chunckSize):
        coordinate_start = np.column_stack((lat_1[i:i+chunckSize], long_1[i:i+chunckSize]))
        coordinate_end = np.column_stack((lat_2[i:i+chunckSize], long_2[i:i+chunckSize]))
        result[i:i+chunckSize] = haversine_vector(coordinate_start, coordinate_end, Unit.KILOMETERS)
    return result

在我的机器上，使用 double 精度，上面的代码大约需要 800 MB，而初始实现需要 8 GB。因此，内存减少 10 倍！它仍然23 倍！使用简单精度，上面的代码大约需要 500 MB，因此 内存减少了 16 倍，速度提高了 48 倍!

更好的 python 减少内存使用的代码？

Better python code for reduce memory usage?

python

performance

haversine

google-colaboratory