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 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)

以上是代码分享,希望能帮到需要的人,有错误的地方望告知。
作为编程小白,希望能在这里记录自己的一步步成长。

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值