(re-posted question!)
I am trying to use the Obspy get_ray_paths_geo tool. I already installed geographiclib
Requirement already satisfied: geographiclib in c:\users\pahan\anaconda3\lib\site-packages (2.0)
, but when plotting what I need I get this message:
UserWarning: Not able to evaluate positions of points on path. Arrivals object will not be modified. Install the Python module 'geographiclib' to solve this issue. warnings.warn(msg)
The code that I am using is:
from obspy.taup import TauPyModel
from obspy import read_events
"python get-pip.py"
lines_latF = [None]*len(st); lines_longF = [None]*len(st)
%pip install geographiclib
for i in range(len(st)):
lines_latF[i] = float(lines_lat[i])
lines_longF[i] = float(lines_long[i])
figTT[i]=plt.figure(figsize=(1080,1080))
model = TauPyModel(model='iasp91')
arrivals = model.get_ray_paths_geo(lines_depth[i], lines_latF[i],
lines_longF[i], InjLat, InjLong, phase_list=('P', 'PP'))
arrivals.plot_rays(plot_type='spherical', phase_list=['P', 'PP'],
legend=True)
plt.savefig("Ray_paths" +str(OriginTimeStringIni[i][0:13:1])+ "_"
+str(st[i].stats.station)+ ".jpg")
i += 0
Any idea about this? TIA!
