make-rotated-copies.py

import sys
import os
import shutil
import itertools as itr
import segyio
import segyio.su as su

def product(f):
    return itr.product(range(len(f.ilines)), range(len(f.xlines)))

def pathjoin(prefix, path, directory):
    _, base = os.path.split(path)
    return os.path.join(directory, '-'.join((prefix, base)))

# this program copies the source-file and creates eight copies, each with a
# modified set of CDP-X and CDP-Y coordinates, rotating the field around the
# north (increasing CDP-Y) axis.
def main():
    if len(sys.argv) < 2:
        sys.exit("Usage: {} [source-file] [destination-file] [prefix]".format(sys.argv[0]))

    srcfile = sys.argv[1]
    dstfile = sys.argv[2] if len(sys.argv) > 2 else srcfile
    prefix  = sys.argv[3] if len(sys.argv) > 3 else os.path.split(dstfile)

    for pre in ['normal', 'acute', 'right', 'obtuse', 'straight', 'reflex', 'left', 'inv-acute']:
        fname = pathjoin(pre, dstfile, prefix)
        shutil.copyfile(srcfile, fname)

        with segyio.open(srcfile) as src, segyio.open(fname, 'r+') as dst:
            for i in range(1 + src.ext_headers):
                dst.text[i] = src.text[i]

            dst.bin = src.bin
            dst.trace = src.trace
            dst.header = src.header

    with segyio.open(pathjoin('normal', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = x
            trh[su.cdpy] = y
            trh[su.scalco] = 10

    with segyio.open(pathjoin('acute', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[y + x * len(src.ilines)]
            trh[su.cdpx] = x + y
            trh[su.cdpy] = (100 - x) + y
            trh[su.scalco] = -10

    with segyio.open(pathjoin('right', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = y
            trh[su.cdpy] = 100 - x
            trh[su.scalco] = 1

    with segyio.open(pathjoin('obtuse', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = (100 - x) + y
            trh[su.cdpy] = (100 - x) - y
            trh[su.scalco] = 2

    with segyio.open(pathjoin('straight', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = 100 - x
            trh[su.cdpy] = 100 - y
            trh[su.scalco] = -7

    with segyio.open(pathjoin('reflex', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = 100 - (x + y)
            trh[su.cdpy] = 100 + (x - y)
            trh[su.scalco] = 7

    with segyio.open(pathjoin('left', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = 100 - y
            trh[su.cdpy] = x
            trh[su.scalco] = 21

    with segyio.open(pathjoin('inv-acute', dstfile, prefix), 'r+') as dst:
        for i, (x, y) in enumerate(product(src)):
            trh = dst.header[i]
            trh[su.cdpx] = 100 + x - y
            trh[su.cdpy] = x + y
            trh[su.scalco] = 100

if __name__ == '__main__':
    main()