xumi1993 / seispy

Python module of seismology and receiver functions

Home Page:https://seispy.xumijian.me

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Moveoutcorrect module

Hannah0363 opened this issue · comments

Hi!!! I've tried the example of "moveout correction" shown in the video at 01:30:30. But, there was an error of plotting.
Code
rf_corr = rfsta.moveoutcorrect(dep_range=np.arange(0,200),ref_rayp=0.06)
correct_stack = np.mean(rf_corr,axis=0)
plt.plot(rfsta.time_axis,raw_stack)
plt.plot(rfsta.time_axis,correct_stack)
plt.xlim([-1,25])
Error
ValueError: x and y must have same first dimension, but have shapes (13001,) and (99, 13001)

Then, I checked the shape of "rf_corr"
AttributeError: 'tuple' object has no attribute 'shape'

However, the shape of "correct_stack" is (99, 13001)!!! Why??
Shouldn't "correct_stack" be like (13001,) and "rf_corr" is (99, 13001)?

In the current version, the code can be written as

rf_corr, t_corr = rfsta.moveoutcorrect(dep_range=np.arange(0,200),ref_rayp=0.06)

It works! Thanks for your timely help!!!