Python读取MICAPs气象数据文件实战指南

# 导入所需的库
import numpy
# 打开 micaps 文件
f1 = open('13052520.000', 'rt')
f2 = open('data.txt', 'wt')
# 前面4行为注释数据，没有用
for i in range(4):
 temp = f1.readline()
# 定义一个数组保存文件中的数据，29行53列
data = numpy.zeros((29, 53))
# 一行一行地读取，读取的是一个字符串
for lin in range(29):
 con = 0
 # 文件中每6行构成所需的1行
 for i in range(6):
  # 读1行
  dd = f1.readline()
  # 将读取的\n去掉，将字符按空格分开，构成了若干个字符串构成的列表
  dd = dd.split()
  # 下面是将字符串一个一个分开导入数组中
  for d in dd:
   data[lin][con] = float(d)
   con = con + 1

# 将数组转化为列表
data = data.tolist()
# 一个一个地写
# 因为写文件只能写字符串,所以需要将列表转化为字符串
for i in range(len(data)):
 for j in range(len(data[i])):
  f2.write(str(data[i][j]))
  f2.write(' ')
 # 每写完一行，换行
 f2.write('\n')
f1.close()
f2.close()