Convert quakeML to hypoinverse format for hyp2000

Junli_Zhang · October 1, 2021, 4:53pm

Hello All,

I would like to convert quakeML to hypoinverse format in order to run hyp2000 for earthquake relocation. Is there a way in obspy to do that?

Thank you very much,
Junli

trichter · October 5, 2021, 11:29am

Hi Junli,

welcome to the forum.

ObsPy does not have hypoinverse write support.

You find a list of supported formats for Catalog.write() in the documentation.

https://docs.obspy.org/packages/autogen/obspy.core.event.catalog.Catalog.write.html

Junli_Zhang · October 7, 2021, 9:40pm

Hello @trichter,

Thank you very much for your reply!
Got it.

Best,
Junli

Miners · August 24, 2022, 4:51am

Hello @Junli_Zhang
Are you able to figure out the Quakeml to Hypo2000 conversion? Any support will be of great help.
Regards
MP

megies · August 24, 2022, 1:03pm

I’m pretty sure this is writing out data in the format you are talking about:

github.com

megies/obspyck/blob/b1d8a9b76291333d88093c84b01b1d0077977b5b/obspyck/obspyck.py#L3583


      
                  # if sensor is buried or downhole, account for the specified sensor
                  # depth
                  depth = stats.coordinates.get('local_depth')
                  if depth:
                      ele -= depth
                  hypo71_string += fmt % (sta_map[sta], lat_deg, lat_min, hem_NS,
                                          lon_deg, lon_min, hem_EW, ele)
          
          
    return hypo71_string
          
          
def dicts2hypo71Phases(self):
              """
              Returns the pick information in hypo71 phase file format
              as a string. This string can then be written to a file.
          
          
    Information on the file formats can be found at:
              http://geopubs.wr.usgs.gov/open-file/of02-171/of02-171.pdf p.30
          
          
    Quote:
              The traditional USGS phase data input format (not Y2000 compatible)
              Some fields were added after the original HYPO71 phase format

megies · August 24, 2022, 1:06pm

If it is, just a heads up on two issues:

that file format only allows 4 characters for station codes, so you will have to do some kind of mapping back and forth when writing input for hypo inverse and when loading back its output
since S phase is written as “time after P” I think you can’t use stations with only an S wave set (at least thats what I figured back then)

Miners · August 24, 2022, 5:22pm

Thank you. Will definitely try this out.

Miners · August 26, 2022, 5:56am

I tried using it but had some issue running it. Amy help?

megies · August 29, 2022, 9:05am

Well, you can’t just run it since it’s embedded in a workflow, you’ll have to copy/paste and adjust the parts that create the output from the event data