本文主要是介绍HYPODD输入文件.pha读取,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
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 sysdef 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_lineresults_dict[key]['rcd'] = temp_listtemp_list = []elements = sgl_line.split()key = " ".join(elements[1:7])head_line = sgl_lineelse:temp_list.append(sgl_line)if i == (len(all_lines)-1):results_dict[key] = dict()results_dict[key]['hl'] = head_lineresults_dict[key]['rcd'] = temp_listreturn results_dictif __name__ == "__main__":file_path = "./test.pha"results_dict = read_data_from_pha_file(file_path)print(results_dict)
以上是代码分享,希望能帮到需要的人,有错误的地方望告知。
作为编程小白,希望能在这里记录自己的一步步成长。
这篇关于HYPODD输入文件.pha读取的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!