Using miniseed file in dispersion_curves.py

[jimmytseng79]

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

Here is my reference miniseed file.
https://goo.gl/VFzXuz