'''
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'))