mirror of
https://github.com/wassname/segpy.git
synced 2026-07-25 13:30:14 +08:00
57 lines
865 B
Python
57 lines
865 B
Python
#
|
|
#
|
|
# testsegy.py
|
|
#
|
|
import struct
|
|
|
|
import segypy
|
|
|
|
|
|
filename='mini.sgy'
|
|
#filename='data_IEEE.segy';
|
|
#filename='data_IBM_REV1.segy';
|
|
#filename='data_IBM_REV0.segy';
|
|
#filename='data_1byteINT.segy';
|
|
#filename='data_2byteINT.segy';
|
|
#filename='data_4byteINT.segy';
|
|
|
|
f_ieee='data_IEEE.segy';
|
|
f_ibm='data_IBM_REV1.segy';
|
|
|
|
filename='data_IEEE.segy';
|
|
|
|
# Set verbose level
|
|
segypy.verbose=1;
|
|
|
|
# Read Segy File
|
|
Data,SH,STH=segypy.readSegy(filename);
|
|
|
|
|
|
filename_out='testout.segy';
|
|
segypy.writeSegyStructure(filename_out,Data,SH,STH);
|
|
segypy.wiggle(Data,SH,1,2)
|
|
|
|
|
|
d_ieee,SH,STH=segypy.readSegy(f_ieee);
|
|
d_ibm,SH,STH=segypy.readSegy(f_ibm);
|
|
|
|
|
|
import pylab
|
|
|
|
|
|
#imshow(Data)
|
|
pylab.figure(1)
|
|
pylab.imshow(d_ieee)
|
|
pylab.title('ieee')
|
|
pylab.show()
|
|
|
|
pylab.figure(2)
|
|
pylab.imshow(d_ibm)
|
|
pylab.title('IBM')
|
|
pylab.show()
|
|
|
|
pylab.figure(3)
|
|
pylab.imshow(Data)
|
|
pylab.title('TEST')
|
|
pylab.show()
|