about.py

import sys
from segyio import TraceField
import segyio


def list_byte_offset_names():
    print("Available offsets and their corresponding byte value:")
    for x in TraceField.enums():
        print("  %s: %d" % (str(x), x))


if __name__ == '__main__':
    if len(sys.argv) < 4:
        list_byte_offset_names()
        sys.exit("Usage: about.py [file] [inline] [crossline]")

    # we need a way to convert from run-time inline/crossline argument (as
    # text) to the internally used TraceField enum. Make a string -> TraceField
    # map and look up into that. this dictionary comprehension creates that
    fieldmap = {str(x).lower(): x for x in TraceField.enums()}

    filename = sys.argv[1]
    inline_name, crossline_name = sys.argv[2].lower(), sys.argv[3].lower()

    # exit if inline or crossline are unknown
    if inline_name not in fieldmap:
        list_byte_offset_names()
        sys.exit("Unknown inline field '%s'" % sys.argv[2])

    if crossline_name not in fieldmap:
        list_byte_offset_names()
        sys.exit("Unknown crossline field '%s'" % sys.argv[3])

    inline, crossline = fieldmap[inline_name], fieldmap[crossline_name]

    with segyio.open(filename, "r", inline, crossline) as f:
        print("About '%s':" % filename)
        print("Format type: %s" % f.format)
        print("Offset count: %d" % f.offsets)
        print("ilines: %s" % ", ".join(map(str, f.ilines)))
        print("xlines: %s" % ", ".join(map(str, f.xlines)))

    print("+------+")

    with segyio.open(filename, "r", crossline, inline) as f:
        # with swapped inline/crossline
        print("About '%s':" % filename)
        print("Format type: %s" % f.format)
        print("Offset count: %d" % f.offsets)
        print("ilines: %s" % ", ".join(map(str, f.ilines)))
        print("xlines: %s" % ", ".join(map(str, f.xlines)))