Rotations fixed

examples/scripts/plot_phase_tensor_map.py

@@ -15,11 +15,12 @@
 import mtpy.imaging.phase_tensor_maps as pptmaps
 
 # directory containing edis
+
 #edipath = r'C:\mtpywin\mtpy\examples\data\edi2'
 edipath = '/mtpy/examples/data/ET_edi'
 
 # whether or not to save the figure to file
-save = True
+save = False
 
 # full path to file to save to
 savepath = '/tmp'
@@ -28,8 +29,9 @@
 plot_freq = 1e-2
 
 # value to color ellipses by, options are phimin,phimax,skew
-colorby = 'skew'
+colorby = 'strike'
 ellipse_range = [-9, 9]
+rotation_angle = 0
 
 image_fn = 'phase_tensor_map%1is_' % (int(1. / plot_freq)) + colorby + '.png'
 
@@ -61,7 +63,7 @@
                                 background_image='/mtpy/examples/data/gravity/ET_gravity.tif'
  #                               bimg_band=1, #  Optional, set to None by default
                                 #bimg_cmap='viridis'  # Optional, set to 'viridis' by default
-                                )
+                                rotation_angle=rotation_angle)
 
 
 if save:

examples/scripts/plot_phase_tensor_section.py

@@ -28,6 +28,7 @@
                                  plot_tipper = 'yri', # plot tipper ('y') + 'ri' means real+imag
                                  font_size=5,
                                  lw=0.5,
+                                 rotation_angle=20,
                                  ellipse_dict = {'ellipse_colorby':'skew_seg',# option to colour by phimin, phimax, skew, skew_seg
                                                  'ellipse_range':[-12, 12, 3]} # set color limits - default 0,90 for phimin or max,
                                                                          # [-12,12] for skew. If plotting skew_seg need to provide
@@ -37,8 +38,7 @@
 # update ellipse size (tweak for your dataset)
 plotObj.ellipse_size = 2.5
 
-
 plotObj.plot()
 
 #plotObj.save_figure(save_fn = op.join(savepath,'PhaseTensorSection.png'),
-#                    fig_dpi=400) # change to your preferred file resolution
+#                    fig_dpi=400) # change to your preferred file resolution

examples/scripts/plot_strike_roseplot.py

@@ -11,12 +11,11 @@
 import os
 os.chdir(r'C:\mtpywin\mtpy') # change to path where mtpy is installed
 
-
 from mtpy.imaging.plotstrike import PlotStrike
 
 
 # directory containing edis
-edipath = r'C:\mtpywin\mtpy\examples\data\edi_files_2'
+edipath = r'C:\mtpywin\mtpy\data\edifiles2'
 
 
 # full path to file to save to
@@ -27,13 +26,50 @@
 elst = [op.join(edipath,f) for f in os.listdir(edipath) if (f.endswith('.edi'))]# and f.startswith('GL')
 
 
-strikeplot = PlotStrike(fn_list=elst,
-                        fold=False,
-                        show_ptphimin=False,
-                        plot_type=2 # 1 means divide into separate plots for different decades
-                                    # 2 means combine all data into one rose plot
-                                    )
-#strikeplot.save_plot(savepath,
-#                     file_format='png',
-#                     fig_dpi=400 
-#                     )
+### this will plot the estimated strike duplicated across quadrants
+#  Plot type 2 will plot all estimates of  strike into one ploe 
+strike_plot = PlotStrike(fn_list=elst,
+                         fold=False,
+                         plot_type=2)
+
+# If you want to plot the orthogonal estimation
+strike_plot.plot_orthogonal = True
+strike_plot.fig_num = 2
+strike_plot.plot()
+
+# if you want to plot Tipper strike estimates as well
+strike_plot.plot_orthogonal = False
+strike_plot.plot_tipper = True
+strike_plot.fig_num = 3
+strike_plot.plot()
+
+# if you want to rotate the data
+# note this will rotate the data N30W
+strike_plot.rotation_angle = -30
+strike_plot.fig_num = 4
+strike_plot.plot()
+
+# to rotate back
+strike_plot.rotation_angle = 30
+
+# if you want to plot per decade plots
+strike_plot.plot_type = 1
+strike_plot.fig_num = 5
+strike_plot.text_pad = 1.4
+strike_plot.plot()
+
+# if you want to only look at a few period ranges
+# not the range is given in log10 of the period
+strike_plot.plot_range = [-2, 0] 
+strike_plot.fig_num = 6
+strike_plot.plot()
+
+# if you want a vertical orientation instead of horizontal
+strike_plot.plot_orientation = 'v'
+strike_plot.plot_range = 'data'
+strike_plot.fig_num = 7
+strike_plot.plot()
+
+
+
+