Python:HYPODD输入文件.pha读取
此脚本意在读取HypoDD的输入.pha文件将结果保存为字典。因为可能会碰到对.pha文件进行修改,比如修改拾取数据的权重,可能需要修改台站名或台网名,所以此脚本可能可以作为辅助帮助修改.pha文件。
.pha的文件格式如下:
# 2016 10 14 00 00 09.270 42.8085 13.2121 6.63 0 0 0 0 1
T1214 1.3196 1 P
ED03 7.7896 1 S
ED24 3.0596 1 P
ED01 5.6196 1 P
# 2016 10 14 00 01 50.170 42.7388 13.1928 5.58 0 0 0 0 2
ED03 5.2796 1 S
T1218 3.6296 1 S
ED02 8.5396 1 S
......
代码如下:
import os
import sys
def read_data_from_pha_file(file_path):
fileID = open(file_path, "r")
all_lines = fileID.readlines()
fileID.close()
results_dict = dict()
for i in range(0, len(all_lines)):
sgl_line = all_lines[i]
if sgl_line[0] == '#':
if i > 0:
results_dict[key] = dict()
results_dict[key]['hl'] = head_line
results_dict[key]['rcd'] = temp_list
temp_list = []
elements = sgl_line.split()
key = " ".join(elements[1:7])
head_line = sgl_line
else:
temp_list.append(sgl_line)
if i == (len(all_lines)-1):
results_dict[key] = dict()
results_dict[key]['hl'] = head_line
results_dict[key]['rcd'] = temp_list
return results_dict
if __name__ == "__main__":
file_path = "./test.pha"
results_dict = read_data_from_pha_file(file_path)
print(results_dict)
以上是代码分享,希望能帮到需要的人,有错误的地方望告知。
作为编程小白,希望能在这里记录自己的一步步成长。