Dear all
I would like to share my code that converts from seisan to Ascii
TSPAIR format and finally to ASCII 2 colums format.
Any comments please to : tonino.gafm@gmail.com
I know it is not the pythonic way but it can be improves if you want
## De Seisan a Ascii 2 coumnas
## --------------byG0D&+0n1nO-----------------
from obspy import read
import os
import csv
import string
import numpy as np
patho=raw_input('Ruta del Archivo-->')
os.chdir(patho)
listaa=[f for f in os.listdir(patho)]
#print listaa
no_ar=raw_input('Nombre del Archivo-->')
data_in=read(no_ar)
data_in.write('salida.ascii',format='TSPAIR')
data_inn=open('salida.ascii')
data_out=open('evento_T_A.ascii','w')
for i in range(1):
data_inn.next()
for line in data_inn:
data_out.writelines(line)
data_inn.close()
data_out.close()
ff=open('evento_T_A.ascii','r')
tiem=[]
amp=[]
for line in ff:
line=line.strip()
columns=line.split()
Tiempo=columns[0]
Amplitud=columns[1]
Tiemp=str(Tiempo)
if "T" in Tiemp:
tiem.append(Tiemp[11:])
amp.append(Amplitud)
# print Tiemp[11:]
#else:
# print "None"
ff.close()
outo=open('ASCII_XY.txt','w')
for t, a in zip(tiem,amp):
outo.write('{} {}\n'.format(t, a))
outo.close()
os.remove('evento_T_A.ascii')
Hope it can help someone else
best regards.
Gonzalo A Fernandez M