Python读取VAX格式文件失败，求替代Matlab旧版读取方法

问题分析 #

旧版Matlab中fopen(..., 'vaxd')指定的是VAX字节序，配合fread(..., 'float32')读取的是VAX F格式单精度浮点数，你之前误用了pyvax处理双精度的from_vax_d8函数，这是数据异常的核心原因。

修正方案 #

使用pyvax中对应VAX单精度浮点数的转换函数from_vax_f4，代码如下：

import pyvax
import numpy as np

file_path = "testfile.THK"

with open(file_path, 'rb') as f:
    raw_bytes = f.read()

# 转换VAX F格式单精度浮点数为IEEE标准格式
ieee_bytes = pyvax.from_vax_f4(raw_bytes)
# 解析为numpy float32数组
data = np.frombuffer(ieee_bytes, dtype=np.float32)

# 验证：打印前10个值，确认是否接近700
print(data[:10])

备选方案（手动转换） #

如果pyvax仍有问题，可以手动实现VAX F格式到IEEE float32的转换，逻辑基于VAX浮点数的格式定义：

import struct
import numpy as np

def vaxf_to_ieee32(vax_bytes):
    # 按4字节为单位解析VAX F格式的无符号整数
    vax_vals = struct.unpack(f'<{len(vax_bytes)//4}I', vax_bytes)
    ieee_vals = []
    for val in vax_vals:
        sign = (val >> 31) & 0x1
        vax_exp = (val >> 23) & 0xff
        mantissa = val & 0x7fffff
        
        # 转换指数偏移量：VAX偏移128，IEEE偏移127
        ieee_exp = vax_exp - 1 if vax_exp != 0 else 0
        # 组合成IEEE float32的位结构
        ieee_val = (sign << 31) | (ieee_exp << 23) | mantissa
        ieee_vals.append(ieee_val)
    
    # 转换为字节后解析为float32数组
    ieee_bytes = struct.pack(f'<{len(ieee_vals)}I', *ieee_vals)
    return np.frombuffer(ieee_bytes, dtype=np.float32)

# 使用示例
with open("testfile.THK", 'rb') as f:
    raw_data = f.read()
data = vaxf_to_ieee32(raw_data)
print(data[:10])

验证步骤 #

1. 运行代码后查看前几个数值是否接近700
2. 检查数组长度是否约为8000（执行len(data)）

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