Comments (10)
Something is wrong with your error estimation. You achieve a chi-square < 1 with a relatively large RRMS (and then the inversion stops), so your error estimate is definitely too pessimistic. How did you do that? Moreover, I noticed that your model depth is probably too small.
from gimli.
Based on this issue, I would like to know how to improve the appearance of the plotted results. In the above figure, the horizontal axis has a length of around 140, while the vertical axis has a length of around 20. The aspect ratio makes the plot results look unappealing. thanks! @halbmy
mgr.showResult()
from gimli.
As I've told you before, the 20m depth might be too low and you should change this.
But what do you mean with unappealing? Visualization is something you should invest time. pyGIMLi provides useful tools but it is up to the user. The result of the showResult
function is a matplotlib axes object with which you can improve your figure like aspect ratio, labels etc.
from gimli.
Something is wrong with your error estimation. You achieve a chi-square < 1 with a relatively large RRMS (and then the inversion stops), so your error estimate is definitely too pessimistic. How did you do that? Moreover, I noticed that your model depth is probably too small. @halbmy
-
First question about the measurement error
Thank you for the response. I checked the error data in the raw file. If I did reciprocal measurement, I used the mean of two error data. If there is only one measurement file, I directly used the err column in syscal pro machine measurement file.
Here is one plot of the err. I removed data with errors larger than 20.
And this is the apparent rhoa data and relative error:
-
Second question about the depth.
I chose the above depth because this is a measurement for a landfill which is around 150m * 20m. The area below 20m has different components and is not our interested part. Why do I need to use a larger depth? Is it better for inversion?
from gimli.
The stacking errors of instruments are often helpful for data processing and also for having an idea of your error model, but usually not good for weighting individual data. Try data.estimateError()
instead.
Even if you're not interested in structures below 20m depth, it will affect your data and you should account for that, otherwise you might get artifacts or bad data fit.
from gimli.
Good to know!
And may I know what is the recommended depth for the mesh?
Since there should be a different layer under 20m. Should I consider it in the mesh generation process?
I used data.estimateError() and increased the inversion depth to be around 1/3 of length. Here is the new result. The rrms is lower than 10 and chi-square < 1.
from gimli.
looks good. There are a few low apparent resistivity points that cannot be fitted but otherwise good. I guess there won't be much resolution below 25-30m depth. Maybe you can improve things be some zWeight
< 1.
from gimli.
Thanks a lot Thomas. I have a following question: if I use the inverted result as prior and generate a synthetic measurement data file, may I get the data back if I run inversion again by the synthetic measurement? I have this werid question because I want to validate some assumption by a synthetic experiment. Instead of manually generating the prior, I want to use the inverted result as a prior. And then I change part of the resistivity map to see if it can be reflected by the new synthetic inversion.
I find the new inversion is hard to reach a low rrms. It's around 40%. Do you think I can improve the result by changing inversion parameters? Or it's just difficult to run an inversion with such a diverse resistivity map?
from gimli.
I don't quite understand. It sounds like it is a rather simple workflow and if you're inverting synthetic data you should be able to find a model fitting the data even though the result may differ from the synthetic model. I guess you should write down some sort of a code. And maybe discuss this in a separate issue.
from gimli.
Yes, you are right. The workflow is correct. I find there was a mistake in my former script. Thank you!
from gimli.
Related Issues (20)
- 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 HOT 8
- Discussion regarding `errorVals`-`absoluteError`-`relativeError` HOT 4
- Petrophysical joint inversion: smooth the output result HOT 3
- Running pyGIMLi on GHA fails to import HOT 3
- Can the gimli module perform joint HVSR-Dispersion Curve inversion? HOT 2
- The resistivity values ββare inverted after inversion with dipole-dipole array HOT 9
- Matplotlib version >3.8.4 causes failure to load module HOT 5
- colab and /tmp HOT 3
- Marker in the Model is not readable in the final result HOT 2
- Coupled DC-Magnetic Inversion HOT 6
- Problem with abortion criteria in ERT data inversion HOT 4
- Question about GravityModelling operator HOT 2
- 3D topography and paraBoundary settings make strange topography HOT 1
- pg.show() plotting the depth values with an opposite sign HOT 3
- Question convertMesh midcell data HOT 1
- neighborCellCount incorrect for boundary cells? HOT 2
- Can't import Pygimli library after conda installation HOT 2
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.