From 3b7a1ec421d28214e840f80870fd884146125daf Mon Sep 17 00:00:00 2001 From: cultpenguin Date: Mon, 9 May 2005 20:15:39 +0000 Subject: [PATCH] *** empty log message *** --- segypy.py | 37 +++++++++++++++++++++++++++++++------ testsegy.py | 24 +++++++----------------- 2 files changed, 38 insertions(+), 23 deletions(-) diff --git a/segypy.py b/segypy.py index 1095b79..9233afd 100644 --- a/segypy.py +++ b/segypy.py @@ -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' diff --git a/testsegy.py b/testsegy.py index 2913c60..73d110b 100644 --- a/testsegy.py +++ b/testsegy.py @@ -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);