Hello, I'm a master student from National Chen Kung University, Taiwan.
I read the article "Rayleigh wave group velocities at periods of 6-23s across Brazil from ambient noise tomography", and i'm very interested in the tomographic inversion.
I have already finished the cross-correlation function processing, and it output the reference waveform with the miniseed format. I saw the code in dispersion_curves.py that the waveform format is pickle. I have also tried to get the pickle file from obspy with the following code:
import obspy.core as obs
st = obs.read('test.mseed')
st.write('test.pickle', format='PICKLE')
And it prints
AttributeError: 'Stream' object has no attribute 'FTANs'
Then i tried another way.
import pickle
from pysismo import pscrosscorr
f = open(name='test.mseed', mode='rb')
xc = pscrosscorr.CrossCorrelationCollection()
with open('test.pickle', 'wb') as f:
pickle.dump(xc, f, protocol=2)
It can work, but it prints
Exporting FTANs of 0 pairs to file ./output/FTAN/FTAN_test.pdf
and dispersion curves to file ./output/FTAN/FTAN_test.pickle
Hello, I'm a master student from National Chen Kung University, Taiwan.
I read the article "Rayleigh wave group velocities at periods of 6-23s across Brazil from ambient noise tomography", and i'm very interested in the tomographic inversion.
I have already finished the cross-correlation function processing, and it output the reference waveform with the miniseed format. I saw the code in dispersion_curves.py that the waveform format is pickle. I have also tried to get the pickle file from obspy with the following code:
And it prints
AttributeError: 'Stream' object has no attribute 'FTANs'Then i tried another way.
It can work, but it prints
Here is my reference miniseed file. https://goo.gl/VFzXuz