forked from sixty-north/segpy
-
Notifications
You must be signed in to change notification settings - Fork 0
/
testsegy.py
56 lines (40 loc) · 925 Bytes
/
testsegy.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
#!/usr/bin/env python
#
# testsegy.py
#
import struct
import segypy
filename='ld0057_file_00095.sgy';
#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';
# Read Segy File
Data,SH,STH=segypy.readSegy(filename);
exit
filename_out='testout.segy';
SH['DataSampleFormat']=5;
SH['SegyFormatRevisionNumber']=100;
segypy.writeSegyStructure(filename_out,Data,SH,STH);
segypy.wiggle(Data,SH,2,.1,.1)
f_ieee='data_IEEE.segy';
f_ibm='data_IBM_REV1.segy';
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()