Phase tensor map

[MohammadHakimRezayee]

I tried to plot the phase tensor map. However, I could plot for one station but not for all station like a grid.

1. how to create an mth5 file from list of edifies?
2. When I try to plot the phase tensor, encounter the following error as this is the script:

ptmap = PlotPhaseTensorMaps(mt_data = edi_list,
plot_freq = plot_freq , # frequency to plot
fig_size = (4,3), # x, y dimensions of figure
xpad = 0.02, ypad = 0.02, # pad around stations
plot_tipper = 'yr', # 'y' + 'r' and/or 'i' to plot
# real and/or imaginary
edgecolor='k', # a matplotlib colour or None for no borders
lw=0.5, # linewidth for the ellipses
minorticks_on=False, # whether or not to turn on minor ticks
ellipse_colorby='skew', # 'phimin', 'phimax', or 'skew'
ellipse_range = [-12,12,2], # [min,max,step]
ellipse_size=0.01, # scaling factor for the ellipses
arrow_size=0.1,
arrow_head_width=0.002, # scaling for arrows (head width)
arrow_head_length=0.002, # scaling for arrows (head length)
ellipse_cmap='bwr', # matplotlib colormap
station_dict={'id':(5,7)} ,
cb_dict = {'position':
[1.05,0.2,0.02,0.3]}, # colorbar position [x,y,dx,dy]
font_size=5
)

save the plot

ptmap.save_plot(os.path.join(savepath,image_fn))

File ~\AppData\Local\anaconda3\Lib\site-packages\mtpy\imaging\plot_phase_tensor_maps.py:109, in PlotPhaseTensorMaps.init(self, mt_data, **kwargs)
107 # --> plot if desired ------------------------
108 if self.show_plot:
--> 109 self.plot()

File ~\AppData\Local\anaconda3\Lib\site-packages\mtpy\imaging\plot_phase_tensor_maps.py:727, in PlotPhaseTensorMaps.plot(self, fig, save_path, show, raster_file, raster_kwargs)
725 self.plot_xarr = np.zeros(len(self.mt_data))
726 self.plot_yarr = np.zeros(len(self.mt_data))
--> 727 for index, tf in enumerate(self.mt_data.values()):
728 if self.pt_type == "ellipses":
729 plot_x, plot_y = self._get_patch_ellipse(tf)

AttributeError: 'list' object has no attribute 'values'

[wieweiwei]

The values attribute is a property inherent to the dictionary class. When the method is designed to accept an mt_data parameter, you have provided it with edi_list. Upon inspecting the method's implementation, it becomes evident that this is not the intended input. The method requires a dictionary that contains multiple MT (Magnetotelluric) objects. However, what you have passed is a list of full paths to EDI files. This discrepancy between the expected input and the provided argument may lead to an error or unexpected behavior in the method's execution.