首先，在这个问题被标记为重复之前，我知道其他人也问过类似的问题，但似乎没有明确的解释。我正在尝试将二进制文件读入二维数组(在此处记录良好 http://nsidc.org/data/docs/daac/nsidc0051_gsfc_seaice.gd.html )。

header 是一个 300 字节的数组。

到目前为止，我有；

import struct

with open("nt_197912_n07_v1.1_n.bin",mode='rb') as file:
    filecontent = file.read()

x = struct.unpack("iiii",filecontent[:300])

抛出字符串参数长度的错误。

最佳答案

读取数据(简答题)

从标题(见下文)确定网格的大小(n_rowsxn_cols = 448x304)后，您可以使用 numpy.frombuffer 简单地读取数据。

import numpy as np

#...

#Get data from Numpy buffer
dt = np.dtype(('>u1', (n_rows, n_cols)))
x = np.frombuffer(filecontent[300:], dt) #we know the data starts from idx 300 onwards

#Remove unnecessary dimension that numpy gave us
x = x[0,:,:]

'>u1' 指定数据格式，在本例中为大小为 1 字节的 无符号 整数，即大端格式。

用 matplotlib.pyplot 绘图

import matplotlib.pyplot as plt

#...

plt.imshow(x, extent=[0,3,-3,3], aspect="auto")
plt.show()

extent= 选项仅指定轴值，例如您可以将这些更改为纬度/经度(从您的标题中解析)

.unpack() 错误的解释

从 docs for struct.unpack(fmt, string) 开始:

The string must contain exactly the amount of data required by the format (len(string) must equal calcsize(fmt))

您可以通过查看 Format Characters 部分来确定格式字符串 (fmt) 中指定的大小。

你在struct.unpack("iiii",filecontent[:300])中的fmt，指定了4个int类型(你也可以用4i = iiii 为简单起见)，每个大小为 4，需要长度为 16 的字符串。

您的字符串 (filecontent[:300]) 的长度为 300，而您的 fmt 要求的字符串长度为 16，因此出现错误。

.unpack() 的用法示例

例如，阅读您提供的文档，我提取了前 21*6 个字节，其格式为:

a 21-element array of 6-byte character strings that contain information such as polar stereographic grid characteristics

与:

x = struct.unpack("6s"*21, filecontent[:126])

这将返回一个包含 21 个元素的元组。请注意某些元素中的空格填充以满足 6 字节的要求。

>> print x
    # ('00255\x00', '  304\x00', '  448\x00', '1.799\x00', '39.43\x00', '45.00\x00', '558.4\x00', '154.0\x00', '234.0\x00', '
    # SMMR\x00', '07 cn\x00', '  336\x00', ' 0000\x00', ' 0034\x00', '  364\x00', ' 0000\x00', ' 0046\x00', ' 1979\x00', '  33
    # 6\x00', '  000\x00', '00250\x00')

注意事项:

• 第一个参数fmt，"6s"*21是6s重复21的字符串 次。每个格式字符6s代表一个6字节的字符串 (见下文)，这将匹配您在 文件。
• filecontent[:126] 中的数字126 计算为6*21 = 126。
• 请注意，对于 s(字符串)说明符，前面的数字确实 not 意味着重复格式字符 6 次(因为它会 通常用于其他格式字符)。相反，它指定大小 的字符串。 s代表一个1字节的字符串，而6s代表 一个 6 字节的字符串。

更广泛的标题读取解决方案(长)

因为必须手动指定二进制数据，所以这在源代码中可能很繁琐。您可以考虑使用一些配置文件(如 .ini 文件)

此函数将读取 header 并将其存储在字典中，其中的结构来自 .ini 文件

# user configparser for Python 3x
import ConfigParser

def read_header(data, config_file):
    """
    Read binary data specified by a INI file which specifies the structure
    """

    with open(config_file) as fd:

        #Init the config class
        conf = ConfigParser.ConfigParser()
        conf.readfp(fd)

        #preallocate dictionary to store data
        header = {}

        #Iterate over the key-value pairs under the
        #'Structure' section
        for key in conf.options('structure'):

            #determine the string properties
            start_idx, end_idx = [int(x) for x in conf.get('structure', key).split(',')]
            start_idx -= 1 #remember python is zero indexed!
            strLength = end_idx - start_idx

            #Get the data
            header[key] = struct.unpack("%is" % strLength, data[start_idx:end_idx])

            #Format the data
            header[key] = [x.strip() for x in header[key]]
            header[key] = [x.replace('\x00', '') for x in header[key]]

        #Unmap from list-type
        #use .items() for Python 3x
        header = {k:v[0] for k, v in header.iteritems()}

    return header

下面是一个示例 .ini 文件。键是存储数据时使用的名称，值是逗号分隔的一对值，第一个是起始索引，第二个是结束索引。 这些值取自您文档中的表 1。

[structure]
missing_data: 1, 6
n_cols: 7, 12
n_rows: 13, 18
latitude_enclosed: 25, 30

这个函数可以按如下方式使用:

header = read_header(filecontent, 'headerStructure.ini')
n_cols = int(header['n_cols'])

关于python - 在python中读取二进制数据，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/42365830/