print("please provide the path to edi folder")
sys.exit(1)
edifiles = glob.glob(os.path.join(edidir, "*.edi"))
print(
"***************** Stats analysis on the EDI files [periods list] *******************************"
)
edis_obj = EdiCollection(edilist=edifiles)
edis_obj.select_periods(percentage=5.0)
print(
"***************** the periods with high occurance in EDI files **************************"
)
per_freq = edis_obj.get_periods_by_stats(
percentage=99.0) # occur in 99% of EDI files
print(per_freq)
for aper in per_freq:
print(aper, edis_obj.get_period_occurance(aper))
print(
"***************** For each uniqu period, show it's occurance in the EDI files ************"
)
for aper in edis_obj.all_unique_periods:
print(aper, edis_obj.get_period_occurance(aper))
print(
"****** To do the pen3d plotting, pick a value from the above periods, which have high occurance (~100%) ******************"
)
# pen3d plot had issues introduced in penetration.py https://github.com/MTgeophysics/mtpy/commit/817c9f4a6384460d57974fb7a6f80e04a8a4ce97
# change the variable below according to your edi files folder
edidir = r'c:\Users\jpeacock\OneDrive - DOI\Geysers\EDI_Files_birrp\Edited'
savepath = r'c:\Users\jpeacock\OneDrive - DOI\Geysers\jvgr\figures'
if not os.path.isdir(edidir):
print("please provide the path to edi folder")
sys.exit(1)
edifiles = glob.glob(os.path.join(edidir, "*.edi"))
# Create plot for a period index number 1,2,3 ..., 10 for determinant. This may not make sense sometimes.
# change to your preferred file resolution
#pen3d.plot_latlon_depth_profile(edidir, 4, 'det', showfig=True, savefig=True,
# savepath=savepath, fig_dpi=400)
# The recommended way is to use a float value for the period.
#############################################################
# More Stats analysis on the EDI files periods
edis_obj = EdiCollection(edilist=edifiles)
edis_obj.select_periods(percentage=5.0)
per_freq = edis_obj.get_periods_by_stats(percentage=80.0)
for aper in range(-15, 0, -1):
# pen3d.plot_bar3d_depth(edidir, aper)
pen3d.plot_latlon_depth_profile(edidir, aper, showfig=True, savefig=True,
savepath=savepath,
fig_dpi=600)