Plot Pseudosection Phase tensors

Give an ellipse representation of phase tensors in pseudo-section format.

# Author: L.Kouadio
# Licence: BSD-3-clause

Plot the phase tensors with different parameters for ellipses such as [ 'phimin' | 'phimax' | 'skew' |'skew_seg' | 'phidet' |'ellipticity' ] For demonstration, we will use the data collected from Huayuan (watex.datasets.load_huayuan()) area with 27 samples. Here, we will give three representations of phase tensors ellipsis.

Import required modules as:

import watex as wx
from watex.view import TPlot
# fetch the Huayuan data
edi_samples = wx.fetch_data ('huayuan', samples =27 , return_data = True )
# can also be
# >>> from watex.datasets import load_huayuan
# >>> edi_samples = load_huayuan ( samples =27, return_data =True )

• Ellipse of 'phimin'

tplot= TPlot (fig_size =( 5, 2 )).fit(edi_samples )
# -------------------------------------------------------------------------------------
# TO BUILD THE DOC WITH  MATPLOTLIB >3.5.3, COMMENT
# THE LINES STARTING WITH  "tplot." SEE RELEASE NOTES v0.1.6 FOR THE REASON
# -------------------------------------------------------------------------------------
#
# We can skip the ellip_dic config by using the default customization
# from the Matplotlib conventional colormap instead.
ellip_dic = {'ellipse_colorby':'phimin',
            'ellipse_range':[0 , 90], # Color limits for skew
            'ellip_size': 2,
            'ellipse_cmap': 'mt_bl2wh2rd'# or color 'bwr' defined in `cmap` parameter
            }
tplot.plot_phase_tensors (ellipse_dict= ellip_dic) # by default the color limit for 'phmin' is [0, 90]

Reading 27 stations
-------------------------
phimin min = -65.62
phimin max = 72.32
-------------------------

<'TPlot':survey_area=None, distance=50.0, prefix='S', window_size=5, component='xy', mode='same', method='slinear', out='srho', how='py', c=2>

• Ellipse of 'skew' visualization

we can setup the ellipses dictionary and passed to the parameter ellipse_dict:

ellip_dic_sk = {'ellipse_colorby':'skew',
            'ellipse_range':[-1 , 1], # Color limits for skew
            'ellip_size': 2,
            'ellipse_cmap':'PuOr' # or color defined in `cmap` parameter such as `mt_bl2wh2rd`
            }
# tplot.plot_phase_tensor (ellipse_dict=ellip_dic_sk )
# or simply turn on `tensor` parameter to ``skew`` like:
tplot.plot_phase_tensors (tensor='skew')

Reading 27 stations
-------------------------
skew min = -88.72
skew max = 86.05
-------------------------

<'TPlot':survey_area=None, distance=50.0, prefix='S', window_size=5, component='xy', mode='same', method='slinear', out='srho', how='py', c=2>

• Determinant of phase tensor 'phidet'

tplot.plot_phase_tensors (tensor='phidet')

Reading 27 stations
-------------------------
phidet min = 3.39
phidet max = 737.21
-------------------------

<'TPlot':survey_area=None, distance=50.0, prefix='S', window_size=5, component='xy', mode='same', method='slinear', out='srho', how='py', c=2>