Example #1
0
'''
Created on 26. nov. 2010

@author: karstm
'''
import os

from processgpr.ext.hubra import HUBRA


#wdir = r'/mn/terra/hubradata/Hubra_Svalbard/Kongsvegen_2008'
wdir = r'Y:\Hubra_Svalbard\Kongsvegen_2008'
file = 'data Thu Apr 24 13-15-42 2008 gated  sweep 1 MHz0'
pro = HUBRA()

pro.convertHDF(os.path.join(wdir, file), wdir,
               1200, combine=True, stack=1, filerange="all-all",
               envelope=False, gpsfile=os.path.join(wdir, file + '.cor'),
               projection='UTM', utmzone=33)

print "Finished conversion for file: %s" % os.path.abspath(file)
print "Total Distance: %f (%i Traces)" % (pro.coords.total_distance, pro.traces)
'''
Converts Hubra data to time-domain.

Traces 500-1594 of the 20-fold stacked data are used in Paper II.

Created on 26. nov. 2010

@author: karstm
'''
import os

from processgpr.ext.hubra import HUBRA


wdir = r'/mn/terra/hubradata/Data2008_09/Wisdom/1502'
file = 'data1'
pro = HUBRA()

pro.convertHDF_prestack(os.path.join(wdir, file), wdir,
                        1000, combine=True, stack=20, filerange="all-all",
                        envelope=False)

print "Finished conversion for file: %s" % os.path.abspath(file)
try:
    print "Total Distance: %f (%i Traces)" % (pro.coords.total_distance, pro.traces)
except AttributeError:
    pass

os.rename(os.path.join(wdir, file + '.h5'),
          os.path.join(wdir, 'AT2009-02-15_pha_stk20.h5'))