MTPy User Guide.zip_大地电磁_大地电磁反演

<html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta charset="utf-8"> <meta name="generator" content="pdf2htmlEX"> <meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1"> <link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css"> <link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css"> <link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/622bb0643d2fbb0007fbbdef/raw.css"> <script src="https://static.pudn.com/base/js/compatibility.min.js"></script> <script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script> <script> try{ pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({}); }catch(e){} </script> <title></title> </head> <body> <div id="sidebar" style="display: none"> <div id="outline"> </div> </div> <div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622bb0643d2fbb0007fbbdef/bg1.jpg"><div class="t m0 x1 h2 y1 ff1 fs0 fc0 sc0 ls0 ws0">MTPy User Guide</div><div class="t m0 x2 h3 y2 ff2 fs1 fc0 sc0 ls0 ws0">October 10, 2018</div><div class="t m0 x3 h4 y3 ff3 fs2 fc0 sc0 ls0 ws0">1 Introduction</div><div class="t m0 x3 h5 y4 ff2 fs3 fc0 sc0 ls0 ws0">This workbook contains some examples for reading, analysing and plotting processed MT data.</div><div class="t m0 x3 h5 y5 ff2 fs3 fc0 sc0 ls0 ws0">It covers most of the steps available in MTPy. For more details on speciﬁc input parameters and</div><div class="t m0 x3 h5 y6 ff2 fs3 fc0 sc0 ls0 ws0">other functionality, we recommend looking at the mtpy documentation, which can be found at:</div><div class="t m0 x3 h5 y7 ff2 fs3 fc0 sc0 ls0 ws0">https://mtpy2.readthedocs.io/en/develop/master_doc.html.</div><div class="t m0 x4 h5 y8 ff2 fs3 fc0 sc0 ls0 ws0">This workbook is structured according to some of the key modules in MTPy: Core, Analysis,</div><div class="t m0 x3 h5 y9 ff2 fs3 fc0 sc0 ls0 ws0">Imaging, and Modeling.</div><div class="t m0 x3 h6 ya ff3 fs1 fc0 sc0 ls0 ws0">1.1 Getting Started</div><div class="t m0 x3 h5 yb ff2 fs3 fc0 sc0 ls0 ws0">To start with, you will need to make sure MTPy is installed and is working correctly. Please see the</div><div class="t m0 x3 h5 yc ff2 fs3 fc0 sc0 ls0 ws0">installation guide (https://github.com/MTgeophysics/mtpy/wiki/MTPy-installation-guide-on-</div><div class="t m0 x3 h5 yd ff2 fs3 fc0 sc0 ls0 ws0">Windows-PC) for details.</div><div class="t m0 x4 h5 ye ff2 fs3 fc0 sc0 ls0 ws0">Before you begin these examples, we suggest you make a temporary folder (e.g. C:/tmp) to</div><div class="t m0 x3 h5 yf ff2 fs3 fc0 sc0 ls0 ws0">save all example outputs.</div><div class="t m0 x3 h6 y10 ff3 fs1 fc0 sc0 ls0 ws0">1.2 Useful tricks and tips</div><div class="t m0 x3 h5 y11 ff2 fs3 fc0 sc0 ls0 ws0">This workbook exists as a Jupyter notebook and a pdf. If you are running the Jupyter notebook,</div><div class="t m0 x3 h5 y12 ff2 fs3 fc0 sc0 ls0 ws0">you can run each of the cells, modifying the inputs to suit your requirements. Most of these</div><div class="t m0 x3 h5 y13 ff2 fs3 fc0 sc0 ls0 ws0">examples have been written to be self contained.</div><div class="t m0 x4 h5 y14 ff2 fs3 fc0 sc0 ls0 ws0">In Jupyter, you can add the following line to the top of any cell and it will write the contents</div><div class="t m0 x3 h5 y15 ff2 fs3 fc0 sc0 ls0 ws0">of that cell to a python script: %%writeﬁle example.py</div><div class="t m0 x4 h5 y16 ff2 fs3 fc0 sc0 ls0 ws0">Many of the examples below make use of the matplotlib colour maps. Please see</div><div class="t m0 x3 h5 y17 ff2 fs3 fc0 sc0 ls0 ws0">https://matplotlib.org/examples/color/colormaps_reference.html for colour map options.</div><div class="t m0 x3 h4 y18 ff3 fs2 fc0 sc0 ls0 ws0">2 Core</div><div class="t m0 x3 h5 y19 ff2 fs3 fc0 sc0 ls0 ws0">These ﬁrst few examples cover some of the basic functions and tools that can be used to look at</div><div class="t m0 x3 h5 y1a ff2 fs3 fc0 sc0 ls0 ws0">data contained in an edi ﬁle, plot it, and make changes (e.g. sample onto different frequencies).</div><div class="t m0 x3 h6 y1b ff3 fs1 fc0 sc0 ls0 ws0">2.1 Read an edi ﬁle into an MT object</div><div class="t m0 x3 h7 y1c ff4 fs3 fc1 sc0 ls0 ws0">In [1]: # import required modules</div><div class="t m0 x5 h7 y1d ff4 fs3 fc3 sc0 ls0 ws0">from mtpy.core.mt import MT</div><div class="t m0 x6 h5 y1e ff2 fs3 fc0 sc0 ls0 ws0">1</div></div><div class="pi" data-data='{"ctm":[1.568627,0.000000,0.000000,1.568627,0.000000,0.000000]}'></div></div> </body> </html>

<div id="pf2" class="pf w0 h0" data-page-no="2"><div class="pc pc2 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622bb0643d2fbb0007fbbdef/bg2.jpg"><div class="t m0 x5 h8 y1f ff5 fs4 fc2 sc0 ls0 ws0"># Define the path to your edi file</div><div class="t m0 x5 h8 y20 ff4 fs4 fc0 sc0 ls0 ws0">edi_file = "C:/mtpywin/mtpy/examples/data/edi_files_2/Synth00.edi"</div><div class="t m0 x5 h8 y21 ff5 fs4 fc2 sc0 ls0 ws0"># Create an MT object</div><div class="t m0 x5 h8 y22 ff4 fs4 fc0 sc0 ls0 ws0">mt_obj = MT(edi_file)</div><div class="t m0 x4 h9 y23 ff2 fs4 fc0 sc0 ls0 ws0">The mt_obj contains all the data from the edi ﬁle, e.g. impedance, tipper, frequency as well as</div><div class="t m0 x3 h9 y24 ff2 fs4 fc0 sc0 ls0 ws0">station information (lat/long). To look at any of these parameters you can type, for example:</div><div class="t m0 x3 h8 y25 ff4 fs4 fc1 sc0 ls0 ws0">In [2]: # To see the latitude and longitude</div><div class="t m0 x5 h8 y26 ff4 fs4 fc3 sc0 ls0 ws0">print mt_obj.lat, mt_obj.lon</div><div class="t m0 x3 h8 y27 ff4 fs4 fc0 sc0 ls0 ws0">-19.01 136.01</div><div class="t m0 x3 h8 y28 ff4 fs4 fc1 sc0 ls0 ws0">In [3]: # To see the easting, northing, and elevation</div><div class="t m0 x5 h8 y29 ff4 fs4 fc3 sc0 ls0 ws0">print mt_obj.east, mt_obj.north, mt_obj.elev</div><div class="t m0 x3 h8 y2a ff4 fs4 fc0 sc0 ls0 ws0">606300.40602 7897760.86059 95.0</div><div class="t m0 x4 h9 y2b ff2 fs4 fc0 sc0 ls0 ws0">There are many other parameters you can look at in the mt_obj. Just type mt_obj.[TAB] to see</div><div class="t m0 x3 h9 y2c ff2 fs4 fc0 sc0 ls0 ws0">what is available. In the MT object are the Z and Tipper objects (mt_obj.Z; mt_obj.Tipper). These</div><div class="t m0 x3 h9 y2d ff2 fs4 fc0 sc0 ls0 ws0">contain all information related to, respectively, the impedance tensor and the tipper.</div><div class="t m0 x3 h8 y2e ff4 fs4 fc1 sc0 ls0 ws0">In [4]: # for example, to see the frequency values represented in the impedance tensor:</div><div class="t m0 x5 h8 y2f ff4 fs4 fc3 sc0 ls0 ws0">print mt_obj.Z.freq</div><div class="t m0 x3 h8 y30 ff4 fs4 fc0 sc0 ls0 ws0">[1.256500e+04 9.751601e+03 7.876300e+03 6.188500e+03 5.250801e+03</div><div class="t m0 x7 h8 y31 ff4 fs4 fc0 sc0 ls0 ws0">4.265799e+03 3.515799e+03 8.437800e+02 6.562798e+02 4.922399e+02</div><div class="t m0 x7 h8 y32 ff4 fs4 fc0 sc0 ls0 ws0">3.867599e+02 3.164400e+02 2.578400e+02 2.109600e+02 1.728900e+02</div><div class="t m0 x7 h8 y33 ff4 fs4 fc0 sc0 ls0 ws0">1.367200e+02 1.015600e+02 7.421900e+01 5.761700e+01 4.882800e+01</div><div class="t m0 x7 h8 y34 ff4 fs4 fc0 sc0 ls0 ws0">4.101600e+01 3.222700e+01 2.636700e+01 2.148400e+01 1.757800e+01</div><div class="t m0 x7 h8 y35 ff4 fs4 fc0 sc0 ls0 ws0">1.440400e+01 1.147500e+01 8.593800e+00 6.591801e+00 5.371100e+00</div><div class="t m0 x7 h8 y36 ff4 fs4 fc0 sc0 ls0 ws0">4.394500e+00 3.601100e+00 2.868700e+00 2.304700e+00 1.914100e+00</div><div class="t m0 x7 h8 y37 ff4 fs4 fc0 sc0 ls0 ws0">1.601600e+00 1.328100e+00 1.074200e+00 8.789100e-01 6.835900e-01</div><div class="t m0 x7 h8 y38 ff4 fs4 fc0 sc0 ls0 ws0">5.078100e-01 3.710900e-01 2.880900e-01 2.050800e-01 1.318400e-01</div><div class="t m0 x7 h8 y39 ff4 fs4 fc0 sc0 ls0 ws0">8.789098e-02 6.835900e-02 5.127000e-02 4.028299e-02 3.295900e-02</div><div class="t m0 x7 h8 y3a ff4 fs4 fc0 sc0 ls0 ws0">2.685500e-02 2.197300e-02 1.709000e-02 1.281700e-02 1.007100e-02</div><div class="t m0 x7 h8 y3b ff4 fs4 fc0 sc0 ls0 ws0">8.239700e-03 6.713900e-03 5.493201e-03 4.272499e-03 2.822900e-03</div><div class="t m0 x7 h8 y3c ff4 fs4 fc0 sc0 ls0 ws0">2.059900e-03 1.678500e-03 1.373300e-03 1.068100e-03 7.629400e-04]</div><div class="t m0 x3 h8 y3d ff4 fs4 fc1 sc0 ls0 ws0">In [5]: # or to see the impedance tensor (first 4 elements)</div><div class="t m0 x5 h8 y3e ff4 fs4 fc3 sc0 ls0 ws0">print mt_obj.Z.z[:4]</div><div class="t m0 x3 h8 y3f ff4 fs4 fc0 sc0 ls0 ws0">[[[ 26.58566 -4.302123j 482.4492 +604.7747j ]</div><div class="t m0 x8 h8 y40 ff4 fs4 fc0 sc0 ls0 ws0">[-410.0502 -800.4257j 8.994784 +44.07396j ]]</div><div class="t m0 x6 h9 y1e ff2 fs4 fc0 sc0 ls0 ws0">2</div></div><div class="pi" data-data='{"ctm":[1.568627,0.000000,0.000000,1.568627,0.000000,0.000000]}'></div></div>

<div id="pf3" class="pf w0 h0" data-page-no="3"><div class="pc pc3 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622bb0643d2fbb0007fbbdef/bg3.jpg"><div class="t m0 x7 h8 y1f ff4 fs4 fc0 sc0 ls0 ws0">[[ 12.43271 +7.519158j 434.8246 +514.6176j ]</div><div class="t m0 x8 h8 y20 ff4 fs4 fc0 sc0 ls0 ws0">[-372.7205 -666.402j 17.64062 +36.09528j ]]</div><div class="t m0 x7 h8 y21 ff4 fs4 fc0 sc0 ls0 ws0">[[ 7.652151 +6.28703j 398.3996 +460.0998j ]</div><div class="t m0 x8 h8 y22 ff4 fs4 fc0 sc0 ls0 ws0">[-349.9875 -580.3959j 21.57495 +33.98854j ]]</div><div class="t m0 x7 h8 y41 ff4 fs4 fc0 sc0 ls0 ws0">[[ 3.59474 +1.225811j 362.5121 +413.2823j ]</div><div class="t m0 x8 h8 y42 ff4 fs4 fc0 sc0 ls0 ws0">[-328.0029 -501.5329j 25.02421 +33.02813j ]]]</div><div class="t m0 x3 h8 y43 ff4 fs4 fc1 sc0 ls0 ws0">In [6]: # or the resistivity or phase (first 4 values)</div><div class="t m0 x5 h8 y44 ff4 fs4 fc3 sc0 ls0 ws0">print mt_obj.Z.resistivity[:4]</div><div class="t m0 x5 h8 y45 ff4 fs4 fc3 sc0 ls0 ws0">print mt_obj.Z.phase[:4]</div><div class="t m0 x3 h8 y46 ff4 fs4 fc0 sc0 ls0 ws0">[[[1.15448560e-02 9.52661629e+00]</div><div class="t m0 x8 h8 y47 ff4 fs4 fc0 sc0 ls0 ws0">[1.28742136e+01 3.22072438e-02]]</div><div class="t m0 x7 h8 y48 ff4 fs4 fc0 sc0 ls0 ws0">[[4.32975088e-03 9.30931663e+00]</div><div class="t m0 x8 h8 y49 ff4 fs4 fc0 sc0 ls0 ws0">[1.19572611e+01 3.31035019e-02]]</div><div class="t m0 x7 h8 y4a ff4 fs4 fc0 sc0 ls0 ws0">[[2.49056438e-03 9.40578869e+00]</div><div class="t m0 x8 h8 y4b ff4 fs4 fc0 sc0 ls0 ws0">[1.16641228e+01 4.11538240e-02]]</div><div class="t m0 x7 h8 y4c ff4 fs4 fc0 sc0 ls0 ws0">[[4.66179794e-04 9.76706091e+00]</div><div class="t m0 x8 h8 y4d ff4 fs4 fc0 sc0 ls0 ws0">[1.16060807e+01 5.54922342e-02]]]</div><div class="t m0 x3 h8 y4e ff4 fs4 fc0 sc0 ls0 ws0">[[[ -9.19198953 51.41945343]</div><div class="t m0 x8 h8 y4f ff4 fs4 fc0 sc0 ls0 ws0">[-117.1256163 78.46525668]]</div><div class="t m0 x7 h8 y50 ff4 fs4 fc0 sc0 ls0 ws0">[[ 31.16505948 49.80396162]</div><div class="t m0 x8 h8 y51 ff4 fs4 fc0 sc0 ls0 ws0">[-119.21840697 63.95414249]]</div><div class="t m0 x7 h8 y52 ff4 fs4 fc0 sc0 ls0 ws0">[[ 39.40662667 49.11076951]</div><div class="t m0 x8 h8 y53 ff4 fs4 fc0 sc0 ls0 ws0">[-121.09061209 57.59379146]]</div><div class="t m0 x7 h8 y54 ff4 fs4 fc0 sc0 ls0 ws0">[[ 18.82944281 48.74426019]</div><div class="t m0 x8 h8 y55 ff4 fs4 fc0 sc0 ls0 ws0">[-123.18471678 52.85011776]]]</div><div class="t m0 x4 h9 y56 ff2 fs4 fc0 sc0 ls0 ws0">As with the MT object, you can explore the object by typing mt_obj.Z.[TAB] to see the available</div><div class="t m0 x3 h9 y57 ff2 fs4 fc0 sc0 ls0 ws0">attributes.</div><div class="t m0 x3 ha y58 ff3 fs5 fc0 sc0 ls0 ws0">2.2 Plot an edi ﬁle</div><div class="t m0 x3 h9 y59 ff2 fs4 fc0 sc0 ls0 ws0">In this example we plot MT data from an edi ﬁle.</div><div class="t m0 x3 h8 y5a ff4 fs4 fc1 sc0 ls0 ws0">In [7]: # import required modules</div><div class="t m0 x5 h8 y1c ff4 fs4 fc3 sc0 ls0 ws0">from mtpy.core.mt import MT</div><div class="t m0 x5 h8 y1d ff4 fs4 fc3 sc0 ls0 ws0">import os</div><div class="t m0 x6 h9 y1e ff2 fs4 fc0 sc0 ls0 ws0">3</div></div><div class="pi" data-data='{"ctm":[1.568627,0.000000,0.000000,1.568627,0.000000,0.000000]}'></div></div>

<div id="pf4" class="pf w0 h0" data-page-no="4"><div class="pc pc4 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622bb0643d2fbb0007fbbdef/bg4.jpg"><div class="t m0 x5 h8 y1f ff5 fs4 fc2 sc0 ls0 ws0"># Define the path to your edi file and save path</div><div class="t m0 x5 h8 y20 ff4 fs4 fc0 sc0 ls0 ws0">edi_file = "C:/mtpywin/mtpy/examples/data/edi_files_2/Synth00.edi"</div><div class="t m0 x5 h8 y5b ff4 fs4 fc0 sc0 ls0 ws0">savepath = r"C:/tmp"</div><div class="t m0 x5 h8 y22 ff5 fs4 fc2 sc0 ls0 ws0"># Create an MT object</div><div class="t m0 x5 h8 y5c ff4 fs4 fc0 sc0 ls0 ws0">mt_obj = MT(edi_file)</div><div class="t m0 x5 h8 y42 ff5 fs4 fc2 sc0 ls0 ws0"># To plot the edi file we read in in Part 1 & save to file:</div><div class="t m0 x5 h8 y5d ff4 fs4 fc0 sc0 ls0 ws0">pt_obj = mt_obj.plot_mt_response(plot_num=1, # 1 = yx and xy; 2 = all 4 components</div><div class="t m0 x9 h8 y5e ff5 fs4 fc2 sc0 ls0 ws0"># 3 = off diagonal + determinant</div><div class="t m0 x6 h8 y5f ff4 fs4 fc0 sc0 ls0 ws0">plot_tipper = 'yri',</div><div class="t m0 x6 h8 y60 ff4 fs4 fc0 sc0 ls0 ws0">plot_pt = 'y' # plot phase tensor 'y' or 'n'</div><div class="t m0 xa h8 y61 ff4 fs4 fc0 sc0 ls0 ws0">)</div><div class="t m0 x5 h8 y62 ff5 fs4 fc2 sc0 ls0 ws0">#pt_obj.save_plot(os.path.join(savepath,"Synth00.png"), fig_dpi=400)</div><div class="t m0 x3 h8 y63 ff4 fs4 fc0 sc0 ls0 ws0"><Figure size 960x720 with 5 Axes></div><div class="t m0 x3 ha y64 ff3 fs5 fc0 sc0 ls0 ws0">2.3 Make some change to the data and save to a new ﬁle</div><div class="t m0 x3 h9 y65 ff2 fs4 fc0 sc0 ls0 ws0">This example demonstrates how to resample the data onto new frequency values and write to a</div><div class="t m0 x3 h9 y66 ff2 fs4 fc0 sc0 ls0 ws0">new edi ﬁle. In the example below, you can either choose every second frequency or resample</div><div class="t m0 x3 h9 y67 ff2 fs4 fc0 sc0 ls0 ws0">onto ﬁve periods per decade. To do this we need to make a new Z object, and save it to a ﬁle.</div><div class="t m0 x3 h8 y68 ff4 fs4 fc1 sc0 ls0 ws0">In [8]: # import required modules</div><div class="t m0 x5 h8 y69 ff4 fs4 fc3 sc0 ls0 ws0">from mtpy.core.mt import MT</div><div class="t m0 x5 h8 y6a ff4 fs4 fc3 sc0 ls0 ws0">import os</div><div class="t m0 x5 h8 y6b ff5 fs4 fc2 sc0 ls0 ws0"># Define the path to your edi file and save path</div><div class="t m0 x5 h8 y6c ff4 fs4 fc0 sc0 ls0 ws0">edi_file = r"C:/mtpywin/mtpy/examples/data/edi_files_2/Synth00.edi"</div><div class="t m0 x5 h8 y6d ff4 fs4 fc0 sc0 ls0 ws0">savepath = r"C:/tmp"</div><div class="t m0 x5 h8 y6e ff5 fs4 fc2 sc0 ls0 ws0"># Create an MT object</div><div class="t m0 x5 h8 y6f ff4 fs4 fc0 sc0 ls0 ws0">mt_obj = MT(edi_file)</div><div class="t m0 x5 h8 y70 ff5 fs4 fc2 sc0 ls0 ws0"># First, define a frequency array:</div><div class="t m0 x5 h8 y71 ff5 fs4 fc2 sc0 ls0 ws0"># Every second frequency:</div><div class="t m0 x5 h8 y72 ff4 fs4 fc0 sc0 ls0 ws0">new_freq_list = mt_obj.Z.freq[::2]</div><div class="t m0 x5 h8 y73 ff5 fs4 fc2 sc0 ls0 ws0"># OR 5 periods per decade from 10^-4 to 10^3 seconds</div><div class="t m0 x5 h8 y74 ff4 fs4 fc3 sc0 ls0 ws0">from mtpy.utils.calculator import get_period_list</div><div class="t m0 x5 h8 y75 ff4 fs4 fc0 sc0 ls0 ws0">new_freq_list = 1./get_period_list(1e-4,1e3,5)</div><div class="t m0 x5 h8 y76 ff5 fs4 fc2 sc0 ls0 ws0"># Create new Z and Tipper objects containing interpolated data</div><div class="t m0 x5 h8 y77 ff4 fs4 fc0 sc0 ls0 ws0">new_Z_obj, new_Tipper_obj = mt_obj.interpolate(new_freq_list)</div><div class="t m0 x5 h8 y78 ff5 fs4 fc2 sc0 ls0 ws0"># Write a new edi file using the new data</div><div class="t m0 x5 h8 y79 ff4 fs4 fc0 sc0 ls0 ws0">mt_obj.write_mt_file(</div><div class="t m0 x6 h9 y1e ff2 fs4 fc0 sc0 ls0 ws0">4</div></div><div class="pi" data-data='{"ctm":[1.568627,0.000000,0.000000,1.568627,0.000000,0.000000]}'></div></div>

<div id="pf5" class="pf w0 h0" data-page-no="5"><div class="pc pc5 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622bb0643d2fbb0007fbbdef/bg5.jpg"><div class="t m0 xb h8 y1f ff4 fs4 fc0 sc0 ls0 ws0">save_dir=savepath,</div><div class="t m0 xb h8 y20 ff4 fs4 fc0 sc0 ls0 ws0">fn_basename='Synth00_5ppd',</div><div class="t m0 xb h8 y5b ff4 fs4 fc0 sc0 ls0 ws0">file_type='edi',</div><div class="t m0 xb h8 y21 ff4 fs4 fc0 sc0 ls0 ws0">new_Z_obj=new_Z_obj, # provide a z object to update the data</div><div class="t m0 xb h8 y22 ff4 fs4 fc0 sc0 ls0 ws0">new_Tipper_obj=new_Tipper_obj, # provide a tipper object</div><div class="t m0 xb h8 y5c ff4 fs4 fc0 sc0 ls0 ws0">longitude_format='LONG', # write longitudes as 'LONG' not ‘LON’</div><div class="t m0 xb h8 y41 ff4 fs4 fc0 sc0 ls0 ws0">latlon_format='dd'# write as decimal degrees (any other input</div><div class="t m0 xc h8 y42 ff5 fs4 fc2 sc0 ls0 ws0"># will write as degrees:minutes:seconds</div><div class="t m0 xb h8 y5d ff4 fs4 fc0 sc0 ls0 ws0">)</div><div class="t m0 x3 h8 y7a ff4 fs4 fc0 sc0 ls0 ws0">C:\mtpywin\mtpy\mtpy\utils\calculator.py:322: RuntimeWarning: invalid value encou ntered in double_scalars</div><div class="t m0 x8 h8 y7b ff4 fs4 fc0 sc0 ls0 ws0">z_rel_err = error/z_amp</div><div class="t m0 x3 h8 y7c ff4 fs4 fc7 sc0 ls0 ws0">Out[8]: 'C:\\tmp\\Synth00_5ppd_1.edi'</div><div class="t m0 x3 hb y7d ff3 fs6 fc0 sc0 ls0 ws0">3 Analysis</div><div class="t m0 x3 h9 y7e ff2 fs4 fc0 sc0 ls0 ws0">The analysis modules contain functions to analyse MT data, for example analyse dimensionality</div><div class="t m0 x3 h9 y7f ff2 fs4 fc0 sc0 ls0 ws0">or strike angle for either a single edi ﬁle or a collection of edi ﬁles.</div><div class="t m0 x3 ha y80 ff3 fs5 fc0 sc0 ls0 ws0">3.1 Look at the Phase Tensor</div><div class="t m0 x3 h9 y81 ff2 fs4 fc0 sc0 ls0 ws0">The phase tensor is an object within the MT object. So when you read in an edi ﬁle, it will auto-</div><div class="t m0 x3 h9 y82 ff2 fs4 fc0 sc0 ls0 ws0">matically contain the phase tensor. To look at the attributes of the phase tensor, you can (as before)</div><div class="t m0 x3 h9 y83 ff2 fs4 fc0 sc0 ls0 ws0">type mt_obj.[TAB] to see what is available. For example, to look at the skew angle (beta):</div><div class="t m0 x3 h8 y84 ff4 fs4 fc1 sc0 ls0 ws0">In [1]: from mtpy.core.mt import MT</div><div class="t m0 x5 h8 y85 ff4 fs4 fc3 sc0 ls0 ws0">import matplotlib.pyplot as plt</div><div class="t m0 x5 h8 y86 ff5 fs4 fc2 sc0 ls0 ws0"># Define the path to your edi file and save path</div><div class="t m0 x5 h8 y87 ff4 fs4 fc0 sc0 ls0 ws0">edi_file = r"C:/mtpywin/mtpy/examples/data/edi_files_2/Synth00.edi"</div><div class="t m0 x5 h8 y88 ff4 fs4 fc0 sc0 ls0 ws0">mt_obj = MT(edi_file)</div><div class="t m0 x5 h8 y89 ff5 fs4 fc2 sc0 ls0 ws0"># look at the skew values as a histogram</div><div class="t m0 x5 h8 y8a ff4 fs4 fc0 sc0 ls0 ws0">plt.hist(mt_obj.pt.beta,bins=50)</div><div class="t m0 x5 h8 y8b ff4 fs4 fc0 sc0 ls0 ws0">plt.xlabel('Skew angle (degree)')</div><div class="t m0 x5 h8 y8c ff4 fs4 fc0 sc0 ls0 ws0">plt.ylabel('Number of values')</div><div class="t m0 x3 h8 y8d ff4 fs4 fc7 sc0 ls0 ws0">Out[1]: Text(0,0.5,'Number of values')</div><div class="t m0 x3 ha y8e ff3 fs5 fc0 sc0 ls0 ws0">3.2 Look at dimensionality</div><div class="t m0 x3 h9 y8f ff2 fs4 fc0 sc0 ls0 ws0">Here, we use the phase tensor to determine which parts of the impedance tensor are 1D, 2D or 3D.</div><div class="t m0 x3 h9 y90 ff2 fs4 fc0 sc0 ls0 ws0">The code returns an array with a value (1, 2 or 3) representing the dimensionality</div><div class="t m0 x6 h9 y1e ff2 fs4 fc0 sc0 ls0 ws0">5</div></div><div class="pi" data-data='{"ctm":[1.568627,0.000000,0.000000,1.568627,0.000000,0.000000]}'></div></div>