gimli-org / gimli

Geophysical Inversion and Modeling Library :earth_africa:

Home Page:https://www.pygimli.org

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Question: How to verify the inverted and observed results accuracy in traveltime inversion?

makeabhishek opened this issue · comments

Hi,
I was trying crosshole traveltime tomography example, I would like to know how can I check the difference between the observed data and predicted data? Can we plot the two data in a graph showing the observed traveltime and inverted traveltime to see the difference?
I tried to extract data as

pickTT = np.array(mgr.data('t'))
modTT = np.array(mgr.inv.response)
plt.plot(pickTT)
plt.plot(modTT)

Is there other better way to check the accuracy of the model?

Also, Why the pygimli showing
When I'm trying to allocate error to the data, error = mgr.estimateError(data=data,absoluteError=1e-6) Its giving me error TypeError: MethodManager.estimateError() got an unexpected keyword argument 'absoluteError'

pyGIMLi - ERROR - <class 'pygimli.physics.traveltime.TravelTimeManager.TravelTimeManager'>.checkError(TravelTimeManager.py:97)

Please read the example
https://www.pygimli.org/_examples_auto/2_seismics/plot_04_koenigsee.html
using the following code

ax, cbar = mgr.showData(firstPicks=True, linewidth=0)
tt.drawFirstPicks(ax, data, mgr.inv.response, marker=None)

which could actually be mgr.showFit(firstPicks=True) (otherwise the apparent velocities are shown).

just implemented (e860be2) mgr.showFit(firstPicks=True) and used it in the mentioned example

Thanks, Do i need to update pygimli for this change?
What I was thinking to draw a line plot for a particular shot (x-axis) and all receivers (y-axis) with observed data and inverted data.

you would have to use the latest dev version: https://www.pygimli.org/installation.html#staying-up-to-date

but in your case this is not necessary, as you can just use these two lines