如何用Python读取METEOSAT卫星的BUFR文件并转为CSV格式？

核心问题说明 #

你当前用FlatTextRenderer得到的是无结构化的纯文本输出，自然难以解析。要提取可用数据，需要直接访问解码后BUFR消息的结构化对象，而非仅渲染成文本。

解决方案步骤 #

1. 先明确BUFR的基本结构 #

BUFR文件由消息(Message)组成，每个消息包含一个或多个子集(Subset)，每个子集里是带编码的气象要素（对应WMO的BUFR码表）。你需要先明确目标要素的代码（比如温度对应020004、纬度对应005001等），可以通过查看FlatTextRenderer输出的完整内容，找到要素对应的代码（通常格式为(BUFR代码) 要素名称）。

2. 结构化提取数据的代码示例 #

下面用pybufrkit结合pandas提取数据并导出CSV：

import os
import pandas as pd
from pybufrkit.decoder import Decoder
from pybufrkit.dataquery import DataQuery

# 目标BUFR文件目录
directory_path = '/data/swadhin/meteosat'
# 定义需要提取的要素（BUFR代码 -> 自定义字段名）
target_elements = {
    "005001": "latitude",
    "006001": "longitude",
    "020004": "air_temperature",
    "012004": "total_cloud_cover"
    # 可根据需求添加更多要素
}

decoder = Decoder()
# 初始化数据查询器，批量查询目标要素
query = DataQuery(",".join(target_elements.keys()))

all_data = []

for filename in os.listdir(directory_path):
    if filename.endswith('.bfr'):
        file_path = os.path.join(directory_path, filename)
        with open(file_path, 'rb') as f:
            bufr_data = f.read()
            bufr_message = decoder.process(bufr_data)
            
            # 遍历每个子集提取数据
            for subset_idx, subset in enumerate(bufr_message.subsets):
                # 查询当前子集的目标要素值
                results = query.execute(subset)
                row = {"filename": filename, "subset_id": subset_idx + 1}
                
                # 将BUFR代码映射为自定义字段名
                for bufr_code, value in zip(target_elements.keys(), results):
                    row[target_elements[bufr_code]] = value
                
                all_data.append(row)

# 转换为DataFrame并导出CSV
df = pd.DataFrame(all_data)
df.to_csv('meteosat_bufr_data.csv', index=False, encoding='utf-8')
print("数据已导出到meteosat_bufr_data.csv")

3. 关键注意事项 #

• 要素补全：如果需要提取更多要素，先通过FlatTextRenderer输出单个文件的完整内容，从中找到对应要素的BUFR代码，添加到target_elements字典即可。
• 异常处理：部分文件或子集可能缺少某个要素，可添加try-except块避免程序中断，比如：

  try:
      results = query.execute(subset)
  except Exception as e:
      print(f"处理子集{subset_idx+1}出错: {e}")
      continue
  

• 替代工具：如果pybufrkit存在兼容性问题，可使用ECMWF官方的ecCodes工具（Python绑定为eccodes-python），它自带完整的BUFR码表，处理卫星数据更稳定。

内容的提问来源于stack exchange，提问作者The Emerging Star