*** empty log message ***

This commit is contained in:
cultpenguin
2005-05-09 20:15:39 +00:00
parent a645b0ea28
commit 3b7a1ec421
2 changed files with 38 additions and 23 deletions
+31 -6
View File
@@ -300,6 +300,7 @@ STH_def["SourceMeasurementUnit"]["descr"] = {1: {
6: "Kilograms (kg)"}}
STH_def["UnassignedInt1"]={"pos":232 ,"type":"int32"} #'int32'); %232
STH_def["UnassignedInt2"]={"pos":236 ,"type":"int32"} #'int32'); %236
##############
# FUNCTIONS
@@ -335,7 +336,7 @@ def getSegyTraceHeader(SH,THN='cdp',data='none'):
pos=THpos+3600+(SH["ns"]*4+240)*(itrace-1);
txt="Reading trace header " + str(itrace) + " of " + str(ntraces) + " " +str(pos)
txt="Reading trace header " + THN + " " + str(itrace) + " of " + str(ntraces) + " " +str(pos)
printverbose(txt,10);
thv[itrace-1],index = getValue(data,pos,THformat,'>',1)
@@ -617,9 +618,33 @@ def printverbose(txt,level=1):
print 'SegyPY',version,': ',txt
##############
# segy class
class SegyTraceheaderClass:
def __init__(self):
self.cdp=0
#def read(filename):
# read_segy(filename)
#
#def write(filename):
# write_segy(filename)
class SegyHeaderClass:
def __str__(self):
return "SegyHeaderClass "
def __init__(self):
self.filename=0
self.Trace = version
def cdp(self):
return "Getting CDP trace header"
def InlineX(self):
return "Getting CDP trace header"
class SegyClass:
STH_def=STH_def
SH_def=SH_def
STH=SegyTraceheaderClass()
SH=SegyHeaderClass()
def __init__(self):
self.THOMAS='Thomas'
+7 -17
View File
@@ -7,26 +7,16 @@ import segypy
filename='mini.sgy'
data = open(filename).read()
# reading SEG-Y header
SH=segypy.getSegyHeader(filename)
print SH
#index=3200;
#Job,size = segypy.getValue(data,index,'l','>')
#Job,size = segypy.getValue(data,index,'long','>')
#index=index+size
#Line,size = segypy.getValue(data,index,'l','>')
#index =index+size
#Reel,size = segypy.getValue(data,index,'l','>')
#index=index+size
#DataTracePerEnsemble, size = segypy.getValue(data,index,'h','>')
#DataTracePerEnsemble, size = segypy.getValue(data,index,'short','>')
#print "GET SEGY GEADER"
#SH=segypy.getSegyHeader(filename)
#print SH
#print SH['Reel']
print SH['Reel']
# Set verbose level
segypy.verbose=1;
Data,SH,STH=segypy.readSegyFast(filename);
# Read Segy File
#Data,SH,STH=segypy.readSegyFast(filename);