Python SEGYFile.writeSU Examples

Example #1

File: sufix.py Project: valeryTech/pygeo

elif (not os.path.isfile(args[0])):
    parser.error('File %s does not exist!' % (args[0], ))

infile = args[0]

if (len(args) > 1):
    outfile = args[1]
else:
    outfile = infile + '.su'

if (options.verbose):
    print(
        '\nSEG-Y --> SU File Converter v1.0\nBrendan Smithyman, July 2011\n\n\tConverting \'%s\' to \'%s\'...\n'
        % (infile, outfile))

sf = SEGYFile(infile, verbose=options.verbose, majorheadersonly=False)
intr = sf.readTraces()
if (options.normalize):
    intr = sf.sNormalize(intr)

if (options.verbose):
    print('Generating output file.')

if (len(intr.shape) == 1):
    intr = [intr]

sf.writeSU(outfile, intr, sf.trhead)

if (options.verbose):
    print('Done!\n')

Example #2

File: vel2stack.py Project: bsmithyman/pygeo

vntr = sfvel.ntr
vns = sfvel.ns

scalco = float(strh0['scalco'])
if (scalco < 0):
  scalco = -1./scalco

scalel = float(strh0['scalel'])
if (scalel < 0):
  scalel = -1./scalel

coordarr = np.zeros((sntr,3), dtype=np.float32)

coordarr[:,0] = np.array([float(trh['sx'])*scalco for trh in sfstack.trhead])
coordarr[:,1] = np.array([float(trh['sy'])*scalco for trh in sfstack.trhead])

newcoordarr = reduceToLocal(coordarr, angle, basis)

velarr = sfvel[:]
outarr = np.zeros((sntr,sns))

vdx = (x1 - x0) / (vntr - 1)

for i in xrange(sntr):
  newx = newcoordarr[i, 0]
  vi = np.round((newx - x0) / vdx)
  outarr[i,:] = velarr[vi,:]

sfstack.writeSU(outfile + '.su', outarr)
os.system('< %s segyhdrs | segywrite endian=0 tape=%s'%(outfile + '.su', outfile))

Example #3

File: sufix.py Project: EvanBianco/pygeo

  parser.error('Please specify a valid filename to convert.')
elif (not os.path.isfile(args[0])):
  parser.error('File %s does not exist!'%(args[0],))

infile = args[0]

if (len(args) > 1):
  outfile = args[1]
else:
  outfile = infile + '.su'

if (options.verbose):
  print('\nSEG-Y --> SU File Converter v1.0\nBrendan Smithyman, July 2011\n\n\tConverting \'%s\' to \'%s\'...\n'%(infile, outfile))

sf = SEGYFile(infile, verbose=options.verbose, majorheadersonly=False)
intr = sf.readTraces()
if (options.normalize):
  intr = sf.sNormalize(intr)

if (options.verbose):
  print('Generating output file.')

if (len(intr.shape) == 1):
    intr = [intr]

sf.writeSU(outfile, intr, sf.trhead)

if (options.verbose):
  print('Done!\n')

Example #4

scalco = float(strh0['scalco'])
if (scalco < 0):
    scalco = -1. / scalco

scalel = float(strh0['scalel'])
if (scalel < 0):
    scalel = -1. / scalel

coordarr = np.zeros((sntr, 3), dtype=np.float32)

coordarr[:,
         0] = np.array([float(trh['sx']) * scalco for trh in sfstack.trhead])
coordarr[:,
         1] = np.array([float(trh['sy']) * scalco for trh in sfstack.trhead])

newcoordarr = reduceToLocal(coordarr, angle, basis)

velarr = sfvel[:]
outarr = np.zeros((sntr, sns))

vdx = (x1 - x0) / (vntr - 1)

for i in xrange(sntr):
    newx = newcoordarr[i, 0]
    vi = np.round((newx - x0) / vdx)
    outarr[i, :] = velarr[vi, :]

sfstack.writeSU(outfile + '.su', outarr)
os.system('< %s segyhdrs | segywrite endian=0 tape=%s' %
          (outfile + '.su', outfile))