search

xumi1993 / seispy

Python module of seismology and receiver functions
https://seispy.xumijian.me
GNU General Public License v3.0
93 stars 29 forks source link

Time-to-depth conversion of SRF #74

Open Wufaxing opened 10 months ago

Wufaxing commented 10 months ago

ls the Time-to-depth conversion of SRF available now?

xumi1993 commented 10 months ago

Is there any errors in the S-wave time-to-depth conversion? Theoretically, Both Ps and Sp phases have the same procedure on this function.

Wufaxing commented 10 months ago

yes, I use LQT to calculate SRFs, when I run rf2depth ccp.cfg, it only search files with comp of R wyb@Wufatian:/mnt/c/Users/22247/Desktop/all_SRF/SRF_1.0$ rf2depth ccp.cfg Traceback (most recent call last): File "/home/wyb/.local/bin/rf2depth", line 8, in sys.exit(rf2depth()) File "/home/wyb/.local/lib/python3.8/site-packages/seispy/rf2depth_makedata.py", line 182, in rf2depth makedata(cpara) File "/home/wyb/.local/lib/python3.8/site-packages/seispy/rf2depth_makedata.py", line 77, in makedata stadatar = RFStation(rfpath, only_r=True) File "/home/wyb/.local/lib/python3.8/site-packages/seispy/rfcorrect.py", line 56, in init self.read_sample(data_path) File "/home/wyb/.local/lib/python3.8/site-packages/seispy/rfcorrect.py", line 76, in read_sample raise FileNotFoundError('No such files with comp of {} in {}'.format(self.comp, data_path)) FileNotFoundError: No such files with comp of R in ./RFresult/TB01

xumi1993 commented 10 months ago

Thanks for reporting this. I'm in holiday now, I'll debug this after back to work at January.

JouSF-1409 commented 10 months ago

I'm developing S receiver function ccp stacking via seispy during this time and may know this one.

This is partially due to class RFStation and it's api. It needs a sample trace to init RFStation and collect neccessary information like stla, stlo .etc. Unfortunately, rf2depth didn't leave space for different comp. So by default, seispy is searching xxx_Rxxx.sac file rather than xxx_Lxxx.sac.

It's a little bit awkward time point, since I am working on ccp moduel right now to add S receiver function support and it requires weeks to write, run and go through a set of test.

Here's a couple of suggestions:

  1. Take a deeper look into rf2depth function
  2. write a py script to do rf2depth rather than rf2depth in commnad.
xumi1993 commented 10 months ago

Hi all In current version, I have written an API for read L component in RFStation as below 

RFStation('path/to/SRF', prime_comp='L')

But this function hasn't been embedded into rf2depth. In the future, the 'rf2depth' will search components order by R, Q, L and Z.