使用zlib在C中压缩数据后Python读取报错：头部校验错误

你的问题核心在于C端用gz系列函数生成的是gzip格式文件，但Python直接用zlib.decompress处理的是原始zlib流，两者格式不匹配导致头部校验失败。下面是具体的调试方向和解决方法：

1. 最直接的解决：改用Python的gzip模块读取 #

gzopen/gzprintf生成的是标准gzip格式的压缩文件（包含gzip头部、校验信息等），Python的gzip模块可以直接识别并处理这种格式，完美适配你原来的np.loadtxt读取逻辑：

import numpy as np
import gzip

# 用gzip.open以文本模式打开文件，直接传给np.loadtxt
with gzip.open("output_file.dat", 'rt') as f:
    data = np.loadtxt(f)

这个方法和你原来读取未压缩文件的逻辑几乎一致，不需要额外处理字符串转数组，最省心。

2. 若坚持用zlib.decompress：指定gzip格式兼容参数 #

如果一定要用zlib.decompress，需要告诉它处理gzip格式（默认只处理原始zlib压缩流），通过设置wbits参数实现：

import numpy as np
import zlib

with open("output_file.dat", 'rb') as f:
    compressed_data = f.read()

# 16 + zlib.MAX_WBITS 表示识别gzip格式
decompressor = zlib.decompressobj(16 + zlib.MAX_WBITS)
decompressed_str = decompressor.decompress(compressed_data) + decompressor.flush()

# 将解压后的文本字符串转为numpy数组
data = np.fromstring(decompressed_str, sep=' ')
# 记得替换成你实际的Nx和Ny值
data = data.reshape(Nx, Ny)

这种方式需要手动处理数组的重塑，适合需要对解压后的字符串做额外处理的场景。

3. 验证C端压缩文件的正确性 #

先排除C端的问题：用系统自带的gzip命令解压你的输出文件：

gzip -d output_file.dat

如果能成功解压出正常的文本数据文件，说明C端的压缩逻辑是对的，问题完全在Python的读取方式；如果解压失败，那需要检查C端代码：

• 确认gzclose(fp1)是否被正确调用（如果程序提前退出没执行gzclose，可能导致文件损坏）
• 检查循环中gzprintf的返回值，是否有写入失败的情况

内容的提问来源于stack exchange，提问作者Tropilio