概述
numpy 提供了高效的二进制文件格式(.npy / .npz),读写速度远快于 csv,且保留 dtype 和 shape 等元信息。
| 格式 | 适用场景 | 速度 | 压缩 |
|---|---|---|---|
| .npy | 单个数组 | 极快 | 否 |
| .npz | 多个数组 | 极快 | 可压缩 |
| .csv / .txt | 与其他工具交互 | 慢 | 否 |
| memmap | 超大文件(>内存) | 按需 | 否 |
.npy— 单数组保存与加载
import numpy as np
arr = np.arange(1000000).reshape(1000, 1000)
# 保存
np.save('data.npy', arr)
# 文件大小 ≈ arr.nbytes + 128 字节(头信息)
# 加载
loaded = np.load('data.npy')
print(loaded.shape) # (1000, 1000)
print(loaded.dtype) # int32(保留原始类型)
# 也可以用 savez 的 save 格式存储
# 注意: np.load 返回的是 ndarray(对 .npy)或 npzfile(对 .npz)
优势:保留 shape、dtype、字节序等所有元信息,加载即原样。
.npz— 多数组保存与加载
未压缩版本
train_data = np.random.rand(10000, 28, 28)
train_labels = np.random.randint(0, 10, 10000)
test_data = np.random.rand(2000, 28, 28)
test_labels = np.random.randint(0, 10, 2000)
# 保存多个数组(键值对)
np.savez('dataset.npz',
x_train=train_data, y_train=train_labels,
x_test=test_data, y_test=test_labels)
# 加载
data = np.load('dataset.npz')
# data 是 npzfile 对象(类似字典,懒加载)
print(data.files) # ['x_train', 'y_train', 'x_test', 'y_test']
print(data['x_train'].shape) # (10000, 28, 28)
print(data['y_train'].shape) # (10000,)
压缩版本(推荐)
# savez_compressed: 压缩存储(zip 压缩)
np.savez_compressed('dataset_compressed.npz',
x_train=train_data, y_train=train_labels,
x_test=test_data, y_test=test_labels)
# 加载方式相同
data = np.load('dataset_compressed.npz')
文本格式(csv / txt)
np.loadtxt()— 加载简单文本
# 数据文件 data.txt:
# 1.0 2.0 3.0
# 4.0 5.0 6.0
# 7.0 8.0 9.0
data = np.loadtxt('data.txt')
print(data)
# [[1. 2. 3.]
# [4. 5. 6.]
# [7. 8. 9.]]
# 常用参数
data = np.loadtxt('data.csv',
delimiter=',', # 分隔符
skiprows=1, # 跳过表头
dtype=np.float32, # 指定类型
usecols=(0, 2), # 只读第 0、2 列
max_rows=100) # 最多读取行数
np.savetxt()— 保存文本
arr = np.random.rand(100, 5)
np.savetxt('output.csv', arr,
delimiter=',',
fmt='%.4f', # 格式: 4 位小数
header='col1,col2,col3,col4,col5',
comments='') # 表头前不加 #
np.genfromtxt()— 处理复杂文本
比 loadtxt 更强大:处理缺失值、混合类型、自动推断。
# 数据含缺失值(用 ? 或空白表示)
data = np.genfromtxt('messy.csv',
delimiter=',',
dtype=none, # 自动推断类型
names=true, # 第 1 行作为列名
missing_values='?', # 缺失值标记
filling_values=0, # 缺失值填充为 0
encoding='utf-8')
# 访问结构化数组的列
print(data['column_name'])
# 常用参数
data = np.genfromtxt('data.csv',
delimiter=',',
skip_header=1, # 跳过的表头行数
skip_footer=1, # 跳过的尾部行数
usecols=(0, 1, 3), # 读取的列索引
dtype=[('name', 'u10'), ('age', 'i4'), ('score', 'f8')])
内存映射(memory-mapped files)
对于超大文件(>内存容量),使用 memmap 按需加载,不会一次性读入内存。
# 1. 创建内存映射文件
fp = np.memmap('large.dat', dtype='float32', mode='w+', shape=(100000, 1000))
# 此时只分配了虚拟地址空间,不占物理内存
# 像普通数组一样操作
fp[0, :] = np.random.rand(1000).astype('float32')
fp[1, :] = np.random.rand(1000).astype('float32')
# 写入磁盘
fp.flush() # 强制刷到磁盘
del fp # 释放映射
# 2. 加载已有的内存映射
fp = np.memmap('large.dat', dtype='float32', mode='r', shape=(100000, 1000))
# 只加载访问到的部分,不会把整个文件读入内存
print(fp[0:10, :5]) # 硬盘只读这 10×5 的数据
# 3. 就地修改
fp = np.memmap('large.dat', dtype='float32', mode='r+')
fp[0, 0] = 3.14 # 直接写入磁盘文件
fp.flush()
memmap 模式
| mode | 说明 |
|---|---|
| 'r' | 只读 |
| 'r+' | 读写(文件必须存在) |
| 'w+' | 创建或覆盖,读写 |
memmap vs 普通加载
# 10gb 的数据集
# ❌ 普通加载:需要 10gb 内存
# arr = np.load('huge.npy')
# ✅ memmap:内存占用极小
arr = np.load('huge.npy', mmap_mode='r') # 以 memmap 方式打开
print(arr[10000:10050]) # 只加载这几行
从标准 io 读取
# 从字符串读取 import io text = "1,2,3\n4,5,6\n7,8,9" data = np.genfromtxt(io.stringio(text), delimiter=',') print(data) # 从字节读取 binary_data = b"1.0 2.0\n3.0 4.0" data = np.loadtxt(io.bytesio(binary_data))
np.savez的懒加载
# npz 文件只在访问时才加载对应的数组(省内存)
data = np.load('dataset.npz')
# 此时只有元数据加载,实际数据还在磁盘
print(data.files) # 列出键名
# 访问时才加载
x_train = data['x_train'] # 现在才从磁盘读取
# 用完后关闭(释放文件句柄)
# data.close()
hdf5 格式(大数据推荐)
对于 gb 以上的数据集,推荐使用 hdf5 格式(更灵活的分层存储):
# pip install h5py
import h5py
# 保存
with h5py.file('dataset.h5', 'w') as f:
f.create_dataset('images', data=images_array, compression='gzip')
f.create_dataset('labels', data=labels_array, compression='gzip')
# 加载(支持切片,不会全读入内存)
with h5py.file('dataset.h5', 'r') as f:
images = f['images']
print(images.shape)
batch = images[0:32] # 只加载 32 张图
格式选择建议
| 场景 | 推荐格式 |
|---|---|
| 单个 numpy 数组 | .npy |
| 多个 numpy 数组 | .npz(压缩用 .npz zipped) |
| 与其他工具(excel/r)共享 | .csv / .txt |
| 超大文件(>内存) | memmap 或 hdf5 |
| 深度学习训练 | .npy / .npz / hdf5 |
| 中等数据 + 可读性 | .csv |
速查表
| 需求 | 代码 |
|---|---|
| 保存单个数组 | np.save('a.npy', arr) |
| 加载单个数组 | np.load('a.npy') |
| 保存多个数组 | np.savez('a.npz', x=x, y=y) |
| 压缩保存 | np.savez_compressed('a.npz', x=x, y=y) |
| 简单文本加载 | np.loadtxt('a.txt') |
| 复杂文本加载 | np.genfromtxt('a.csv', delimiter=',') |
| 保存文本 | np.savetxt('a.csv', arr, delimiter=',') |
| 内存映射(创建) | np.memmap('f.dat', dtype='f4', mode='w+', shape=s) |
| 内存映射(加载) | np.load('f.npy', mmap_mode='r') |
| 强制刷盘 | fp.flush() |
| 列出 npz 键 | data.files |
总结
以上为个人经验,希望能给大家一个参考,也希望大家多多支持代码网。
发表评论