171 lines
4.7 KiB
Python
171 lines
4.7 KiB
Python
from collections import defaultdict
|
|
from datetime import datetime, time
|
|
|
|
def is_blank(l: str) -> bool:
|
|
return len(l.strip(" ")) == 0
|
|
|
|
def parse_flt(v:str) -> float | None:
|
|
try:
|
|
t = float(v)
|
|
return t
|
|
except ValueError:
|
|
return None
|
|
|
|
def parse_int(v:str) -> int | None:
|
|
try:
|
|
t = int(v)
|
|
return t
|
|
except ValueError:
|
|
return None
|
|
|
|
|
|
def parse():
|
|
fp = open("dados.txt")
|
|
data = [l for l in fp.read().split("\n")]
|
|
chunks = boundaries(data)
|
|
|
|
for c in chunks:
|
|
parse_chunk(data[c[0]:c[1]])
|
|
|
|
fp.close()
|
|
|
|
|
|
def boundaries(data: list[str]):
|
|
boundaries = []
|
|
start = None
|
|
for (idx,l) in enumerate(data):
|
|
if start is None:
|
|
if not is_blank(l):
|
|
start = idx
|
|
else:
|
|
if is_blank(l):
|
|
boundaries.append((start,idx))
|
|
start = None
|
|
return boundaries
|
|
|
|
|
|
def parse_chunk(chunk_lines: list[str]):
|
|
hIdx = None
|
|
for (idx, l) in enumerate(chunk_lines):
|
|
if l[-1] == "7":
|
|
hIdx = idx
|
|
break
|
|
headersRet = parse_header(chunk_lines[:hIdx])
|
|
# TODO: implementar o parser das fases parser_type_7
|
|
|
|
return headersRet
|
|
|
|
|
|
def parse_header(hLines: list[str]):
|
|
aux = defaultdict(list)
|
|
|
|
for line in hLines:
|
|
match line[-1]:
|
|
case "1":
|
|
aux[1].append(line)
|
|
case "3":
|
|
aux[3].append(line)
|
|
case "6":
|
|
aux[6].append(line)
|
|
case "E":
|
|
aux["E"].append(line)
|
|
case "I":
|
|
aux["I"].append(line)
|
|
case "F":
|
|
aux["F"].append(line)
|
|
case _:
|
|
raise NotImplemented
|
|
|
|
headerDict = dict()
|
|
for (k,v) in aux.items():
|
|
if len(v) != 0:
|
|
headerDict.update(FUNCS[k](v))
|
|
return headerDict
|
|
|
|
|
|
def parse_mag(line: str):
|
|
magnitudes = []
|
|
base = 55
|
|
while base < 79:
|
|
m = line[base:base+4]
|
|
mt = line[base+4]
|
|
if not is_blank(m):
|
|
magnitudes.append({"M": m, "T": mt})
|
|
base += 8
|
|
return magnitudes
|
|
|
|
|
|
def parse_type_1(data: list[str]):
|
|
aux = data[0]
|
|
y = int(aux[1:5])
|
|
mo = int(aux[6:8])
|
|
d = int(aux[8:10])
|
|
h = int(aux[11:13])
|
|
m = int(aux[13:15])
|
|
s = int(aux[16:18])
|
|
mil = int(aux[19]) * 10**5
|
|
dt = datetime(y,mo,d,h,m,s,mil)
|
|
|
|
dist_ind = aux[21]
|
|
eId = aux[22]
|
|
lat = float(aux[23:30])
|
|
long = float(aux[30:38])
|
|
depth = float(aux[38:43])
|
|
rep_ag = aux[45:48]
|
|
|
|
hypo = {"DateTime": dt.isoformat(), "Distance Indicator": dist_ind, "Event ID": eId, "Lat": lat, "Long": long, "Depth": depth, "Agency": rep_ag, "Magnitudes": list()}
|
|
|
|
for l in data:
|
|
hypo["Magnitudes"] = hypo["Magnitudes"] + parse_mag(l)
|
|
|
|
return hypo
|
|
|
|
def parse_type_3(data: list[str]):
|
|
comments = []
|
|
for line in data:
|
|
comments.append(line[:-2].strip())
|
|
return {"Comments": comments}
|
|
|
|
|
|
def parse_type_6(data: list[str]):
|
|
waves = []
|
|
for l in data:
|
|
waves.append(l.strip().split(" ")[0])
|
|
return {"Wave": waves}
|
|
|
|
|
|
def parse_type_7(data: list[str]):
|
|
phases = []
|
|
# nordic format
|
|
for l in data:
|
|
print(l)
|
|
h = int(l[18:20])
|
|
m = int(l[20:22])
|
|
sec = int(l[23:25])
|
|
mil = int(l[26:28]) * 10**4
|
|
t = time(h,m,sec,mil)
|
|
phases.append({"Stat:":l[1:5], "Com": l[6:10], "I": l[9].strip(), "Phase": l[10:15].strip(), "Polarity": l[16].strip(), "Time": t.isoformat(), "Duration": parse_flt(l[29:33]), "Amplitude": parse_flt(l[34:40]), "Period": parse_flt(l[41:45]), "Azimuth": parse_flt(l[46:51]), "Velocity":parse_int(l[52:56]), "AIN": parse_int(l[57:60]), "AR": l[61:63], "Travel Time": parse_flt(l[63:67]), "Weigth": parse_int(l[67:70]), "Distance": float(l[71:75]), "CAZ": int(l[76:79])})
|
|
|
|
return {"Phases": phases}
|
|
|
|
def parse_type_e(data: list[str]):
|
|
aux = data[0]
|
|
error = {"Gap": int(aux[5:8]), "Origin": float(aux[14:20]), "Error_lat": float(aux[24:30]), "Error_long": float(aux[32:38]), "Error_depth": float(aux[38:43]), "Cov_xy": float(aux[43:55]), "Cov_xz": float(aux[55:67]), "Cov_yz": float(aux[67:79])}
|
|
return error
|
|
|
|
|
|
def parse_type_f(data: list[str]):
|
|
return {}
|
|
|
|
|
|
def parse_type_i(data: list[str]):
|
|
aux = data[0]
|
|
dt = datetime.strptime(aux[12:26], "%y-%m-%d %H:%M")
|
|
return {"Action": aux[8:11], "Action Extra": {"Date": dt.isoformat(), "OP": aux[30:35].strip(), "Status": aux[42:57].strip(), "ID":int(aux[60:74])}}
|
|
|
|
|
|
FUNCS = {1: parse_type_1, 3: parse_type_3, 6: parse_type_6, "E": parse_type_e, "F": parse_type_f, "I": parse_type_i}
|
|
|
|
|
|
parse()
|