Comments (16)
Note that I added the line data["t"] /= 1000
because I guess the traveltime in the file is in obviously ms and not s. This needs to be changed in the import routine (if the TOM file always contains ms) but until then you will have to multiply by yourself.
from gimli.
what is data[t]/=1000
in line 5? Just saw that same comment was made by Prof @halbmy
from gimli.
the traveltimes are divided by 1000
from gimli.
Thank you all for replying, The code is now running fine, but the second step apparent velocity is not showing well. Attached is the screenshot of the code and the displace. Hope to get a response
from gimli.
Note that I added the line
data["t"] /= 1000
because I guess the traveltime in the file is in obviously ms and not s. This needs to be changed in the import routine (if the TOM file always contains ms) but until then you will have to multiply by yourself.
Thank you prof., The code is now running fine, but the second step apparent velocity is not showing well. Attached is the screenshot of the code and the displace. Hope to get a response
from gimli.
As the output says, you seem to have zero traveltime leading to infinite apparent velocity. Remove those data with data.remove(data["t"] <= 0])
.
from gimli.
As the output says, you seem to have zero travel-time leading to infinite apparent velocity. Remove those data with
data.remove(data["t"] <= 0])
.
Thanks prof., the output from the code displaced this; as in the attached file but my apparent velocity is showing range 0.26. to 5.26. I don't know if the range is okay or needs to be multiplied. moreover there are some empty space in the model which i don't know if is okay too.
from gimli.
As the output says, you seem to have zero traveltime leading to infinite apparent velocity. Remove those data with
data.remove(data["t"] <= 0])
.
Thanks prof., the output from the code displaced this; as in the attached file but my apparent velocity is showing range 0.26. to 5.26. I don't know if the range is okay or needs to be multiplied. moreover there are some empty space in the model which i don't know if is okay too.
from gimli.
I am pretty sure you DO need to multiply your data with 0.001 as discussed before. Which wave is so slow to have a traveltime of 60 seconds on a distance of 100m?? The lowest values would then be around 300m/s (the speed in air) and the highest at 5000 m/s (typical values for hard rock). The missing values are probably the ones with zero traveltime.
And please, don't send smartphone photos. You already have a notebook and can copy the output into the issue.
from gimli.
I am pretty sure you DO need to multiply your data with 0.001 as discussed before. Which wave is so slow to have a traveltime of 60 seconds on a distance of 100m?? The lowest values would then be around 300m/s (the speed in air) and the highest at 5000 m/s (typical values for hard rock). The missing values are probably the ones with zero travel-time.
And please, don't send smartphone photos. You already have a notebook and can copy the output into the issue.
Thank you prof., I multiplied 0.001 in the code and my velocity range is still showing 0.26 to 5.26. Attached is my output and script.
SRTD.txt
from gimli.
I am pretty sure you DO need to multiply your data with 0.001 as discussed before. Which wave is so slow to have a traveltime of 60 seconds on a distance of 100m?? The lowest values would then be around 300m/s (the speed in air) and the highest at 5000 m/s (typical values for hard rock). The missing values are probably the ones with zero traveltime.
And please, don't send smartphone photos. You already have a notebook and can copy the output into the issue.
Prof. the code is running fine if i run it in anaconda prompt but in jupyter-notebook is displacing 0.26 to 5.27 m/s for velocity
from gimli.
The plots look correct. But if you run the same code on the command line and in a notebook (using the same pygimli in the same environment), the results should be the same.
from gimli.
The plots look correct. But if you run the same code on the command line and in a notebook (using the same pygimli in the same environment), the results should be the same.
Yes prof. is the same now IF I JOIN all the code together but i run this code below for fit between measured (crosses) and modeled (lines) travel-times and i am having a empty displace cells saying ( Line2D.set() got an unexpected keyword argument 'orientation'). Attached is the code, output and error message.
mgr.invert(secNodes=3, paraMaxCellSize=5.0,
zWeight=0.2, vTop=500, vBottom=5000, verbose=1)
mgr.showFit(firstPicks=True)
from gimli.
Yes, that's definitely an error that need to be corrected. Thank you for point out.
from gimli.
Yes, that's definitely an error that need to be corrected. Thank you for point out.
Good evening Prof.,
Prof., the code is still not working yet. still showing as before
from gimli.
TravelTimeManager.showFit()
(firstPicks
is True by default) shows the measured data as crosses along with the forward response as lines like demonstrated in the traveltime field data example
There was apparently some bug in one of the previous versions which was fixed with the last version v1.4.6.
from gimli.
Related Issues (20)
- Running inversion multiple times give different result HOT 11
- 1. How get j of specify nodes after forward? HOT 2
- Joint inversion module missing HOT 3
- bug-typo/Cm05.cols HOT 1
- The cellMarkers is different before and after the inversion process AND inversion mesh setting HOT 2
- This may be a corrupt mesh error HOT 7
- pyvista shows no result graph when I tried to reproduce the example '3D surface ERT inversion' HOT 2
- ERT: blanked out cells when adding new regions to the mesh HOT 4
- Cannot import name 'RVector3' form 'pygimli.core.core HOT 6
- how to fix import pygimli error libcholmod.so.3 on Mybinder HOT 5
- Plot Inversion Result from Another Software use Pygimli HOT 2
- fitReciprocalErrorModel IP data HOT 4
- Seeking Help with 3D Magnetic Inversion HOT 2
- Question about geometric factor for an ERT survey with submerged and terrestrial electrodes HOT 3
- ERTManager.invert() --> error: malloc(): unsorted double linked list corrupted HOT 3
- Failure to install on mac HOT 6
- Wrong coordinates in tutorial notebook HOT 4
- pyGIMLi and NUM_THREADS HOT 7
- 3d geometry with topography and fixed thin layer HOT 5
- pg.interpolate() giving error with regionwise regularization
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from gimli.