Closed JMcKercher closed 1 week ago
Sorry for the very late response. I've had a look into your scripts. There might be some further modifications necessary (filtering ranges, error estimation) but as to the technical point, you can simple access the IP model by mgr.modelIP
.
However, I suggest NOT using the ratio of the IP values but the difference. If you consider resistivity as complex number of amplitude and phase A*exp(w P), the ratio of two values exhibits the phase difference, i.e. A2 exp(w P2) / (A1 exp(w P1)) = A2/A1 exp(w(P2-P1))
Done by
lowipmgr.showResult(lowipmgr.modelIP - modipmgr.modelIP,
cMap="bwr", logScale=False)
I ended up in this image
Thank you for getting back to me. I am able to return the same response using the provided code.
I consider this issue solved. If not, reopen and describe the problem in more detail.
Problem description
Hello, I have been loving this software and starting to get a good understanding of it. My knowledge of coding is still quite limited and I am unsure if this has already been covered in an example or a tutorial, but I would like to know if it is possible to create an IP ratio (graph and/or inversion plot) from two different inversions or datasets.
Problem: I have collected two sets of data each at different frequencies from the same field transect. I would like to test or evaluate the IP difference between the two frequencies to ultimately determine which frequency is better at detecting a signal from the field. Please let me know if you require further information.
Your environment
Python 3.9.18 | packaged by conda-forge | (main, Dec 23 2023, 16:29:04) [MSCv.1929 64 bit (AMD64)] pygimli : 1.4.6 pgcore : 1.4.0 numpy : 1.21.6 matplotlib : 3.7.2 scipy : 1.11.2 IPython : 8.18.1 pyvista : 0.43.2
Steps to reproduce
Below is the basic inversion of two data sets showing the IP results. I am unsure how to go about what I would like to achieve, but I believe it to be like ratio.
Expected behavior
Create a inversion graphical output displaying the ratio between the two datasets IP models
Actual behavior
None so far, is this possible?
Data to reproduce
Rest point 03-04 cross 0.52.txt Rest point 03-04 cross 4.16.txt