Traveltime Tomography - fixing values for reconstruction. #836
Unanswered
carmsignes
asked this question in
Q&A
Replies: 1 comment 1 reply
-
|
Beta Was this translation helpful? Give feedback.
1 reply
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
Uh oh!
There was an error while loading. Please reload this page.
-
Hello!
I am using pyGIMLi to simulate an imaging of a big tank of water. The idea is to obtain a map of the sound-speed (and therefore, the density and temperature) of the water using a circular array of ultrasonic transducers. I built up my code from one of the useful provided examples. However, there are a few things that I am still trying to implement:
Inside the tank of water, there will be certain solid structures with a known sound-speed and location. Is there a way to hard-fix this values for the traveltime inversion? I could not find a way to do this. One of the things I tried was to flag such obstacles with the boundaryMarker=-1 so they would be excluded, but it did not have the intended result. What would be the best way to proceed? At this moment, if I include such obstacles, I get a good measurement of them, because their sound-speed is very high and therefore first arrivals often go through them. However, I lose a lot of information of the surrounding media (which is the one I am interested in)
Is there any way to incorporate reflections in such structures into the traveltime inversion? For what I found on the documentation, the answer is no, but just wanted to check if there are any known workarounds.
Would it be possible to use a full wave propagation for the forward calculation using other software, for example Kwave , and then inputting the result (after treatment) on pyGIMLi for the inversion? I wondered if someone was experienced on this.
Here is a copy of my code so far:
Initial media

Initial media with obstacle --> Velocities are [1,1482],[2,1550],[3,1528], [4,4000]

Inversion result

Inversion result (adjusted scale)

Thank you kindly!
Beta Was this translation helpful? Give feedback.
All reactions