Entry Point for Making MTH5 files based on User Inputs and Various Ingestion Protocols about mth5

[timronan]

@timronan already created this class but does not have permissions to branch from the code main branch to submit a a Pull request for review.

Class could be /mth5/mth5/utils/make_mth5.py

And to start would contain the code

 from obspy import UTCDateTime
 
 from mt_metadata.timeseries.stationxml import XMLInventoryMTExperiment
 from mth5.mth5 import MTH5
 from mth5.timeseries import RunTS
 import os
 
 
 class Make_MTH5():
 
     def make_mth5_from_iris(network, station, starttime, endtime, path=None,
                             client='IRIS'):
 
         if path is None:
             path = str(os.getcwd()) + '/'
         network = network
         station = station
         start = UTCDateTime(starttime)
         end = UTCDateTime(endtime)
         file_name = path + network + '_' + station + '_' + starttime + '.h5'
         # need to know network, station, start and end times before hand
         client = fdsn.Client(client)
 
         # get the data
         streams = client.get_waveforms(network, station, None, None, start, end)
 
         # get the metadata
         inventory = client.get_stations(
             start, end, network=network, station=station, level="channel"
         )
         # translate obspy.core.Inventory to an mt_metadata.timeseries.Experiment
         translator = XMLInventoryMTExperiment()
         experiment = translator.xml_to_mt(inventory)
 
         # initiate MTH5 file
         m = MTH5()
         m.open_mth5(r"" + file_name, "w")
         # fill metadata
         m.from_experiment(experiment)
         station_group = m.get_station(station)
 
         # runs can be split into channels with similar start times and sample rates
         start_times = sorted(list(set([tr.stats.starttime.isoformat() for tr in streams])))
         end_times = sorted(list(set([tr.stats.endtime.isoformat() for tr in streams])))
 
         for index, times in enumerate(zip(start_times, end_times), 1):
             run_stream = streams.slice(UTCDateTime(times[0]), UTCDateTime(times[1]))
             run_ts_obj = RunTS()
             run_ts_obj.from_obspy_stream(run_stream)
             run_group = station_group.add_run(f"{index:03}")
             run_group.from_runts(run_ts_obj)
         print(file_name + " has been built.")
         m.close_mth5()
         return file_name

or similar.

[timronan]

Closed by merging branch enhancement_42