Python读取Fortran二进制浮点文件出现乱码问题求助

首先，你遇到乱码的核心原因非常明确：你用文本模式打开了二进制文件！

Fortran生成的二进制浮点数文件，存储的是原始的字节序列（比如双精度浮点数占8字节、单精度占4字节），而你用open(..., "r")这种文本模式打开时，Python会默认把这些字节当成UTF-8（或系统默认编码）的文本字符来解析——但二进制浮点数的字节序列根本不是合法的文本编码，自然会输出一堆乱码字符。

接下来给你具体的解决步骤：

1. 必须用二进制模式打开文件 #

先把文件打开方式改成"rb"（read binary），这是读取二进制文件的基础：

f = open(args.model_files[0], "rb")

2. 解析Fortran二进制浮点数的正确方式 #

Fortran的未格式化二进制文件通常有两种常见格式，对应不同的读取方法：

方法一：用Numpy直接读取（推荐，简单高效） #

Numpy对Fortran的二进制格式支持很好，能自动处理大部分场景：

• 如果你的Fortran程序是用直接访问方式写的文件（比如OPEN(unit, FORM='UNFORMATTED', ACCESS='DIRECT')），直接用fromfile读取：

import numpy as np

# 根据Fortran代码选择dtype：双精度用np.float64，单精度用np.float32
# Fortran默认REAL通常是双精度，选float64即可
data = np.fromfile(args.model_files[0], dtype=np.float64)
# 读取第一个值，对应你Fortran里得到的284.69
print(data[0])

• 如果你的Fortran程序是用顺序访问的未格式化文件（比如WRITE(unit) variable），这种文件会在每个数据块前后添加4字节的记录长度标记，需要跳过这些标记再读取：

import numpy as np

with open(args.model_files[0], "rb") as f:
    # 读取开头的记录长度（4字节整数）
    record_length = np.frombuffer(f.read(4), dtype=np.int32)[0]
    # 读取对应长度的浮点数，这里假设是双精度（每个8字节）
    num_values = record_length // 8
    data = np.frombuffer(f.read(record_length), dtype=np.float64)
    # 跳过结尾的记录长度标记
    f.read(4)

print(data[0])

方法二：用struct模块手动解析（适合小数据量） #

如果你不想依赖Numpy，可以用Python内置的struct模块手动解析字节：

import struct

with open(args.model_files[0], "rb") as f:
    # Fortran默认用大端字节序，格式字符串用'>d'（>表示大端，d表示双精度浮点数）
    # 单精度则用'>f'
    # 先跳过记录长度（如果是顺序访问的未格式化文件）
    f.read(4)
    # 读取第一个双精度浮点数的8字节
    byte_data = f.read(8)
    value = struct.unpack('>d', byte_data)[0]
    print(value)

额外注意点 #

• 字节序匹配：Fortran默认使用大端字节序（big-endian），而x86架构的电脑是小端字节序，所以读取时一定要指定字节序（比如Numpy里用'>f8'，struct里用'>d'），否则会得到错误的数值。
• 数据类型一致：确保Python里用的dtype和Fortran里的变量类型匹配——Fortran的REAL通常是双精度（对应Python的float64），REAL(KIND=4)是单精度（对应float32）。

内容的提问来源于stack exchange，提问作者i'mlaguiar