使用Python读SEGY道头
自己遇到的一个问题, WestGeco的3D地震数据, 提取导航数据:
#!/bin/env python
import sys
import struct
try:
f=open(sys.argv[1],'rb')
except (IOError,Exception):
print '''usage:
scriptname segyfilename
'''
sys.exit(1)
#skip EBCDIC header
try:
f.seek(3200)
except Exception:
print 'Oops! your file is broken..'
#read binary header
binhead = f.read(400)
ns = struct.unpack('>h',binhead[20:22])[0]
if ns < 0:
print 'file read error'
sys.exit(1)
#read trace header
while True:
trchead = f.read(240)
if trchead == '':
break
nav = trchead[224:228]+\
trchead[228:232]+\
trchead[200:204]+\
trchead[204:208]+\
trchead[60:64]
#define output format
nl = struct.unpack('>5i',nav)
print 'QL%-4d%9s%10.2f%10.2f%5s' % (nl[0], nl[1], nl[2]/100.0, nl[3]/100.0, nl[4])
f.seek(ns*4,1)
f.close()
相关文档：
1． 首先就是在编译器中把python安装目录include/与libs/加入，对于这点我在vc6中可以，但是在dev c++中即使加入了编译也会出错，说找不到python头文件，这点比较郁闷，不过考虑到一般windows编程都用的是vc，所以并没有什么影响吧！！！
然后用#include <Python.h>就可以把python的主头文件包含进来了。
但� ......
1，编写Server.py
import socket
s=socket.socket(socket.AF_INET,socket.SOCK_DGRAM)
s.bind(('localhost',8081))
while True:
data,addr=s.recvfrom(1024)
print 'receive:',data,'from',addr
2，编写Client.py
import socket
s=socket.socket(socket.AF_INET,socket.SOC ......
需要先安装libxml2-devel libxslt-devel这两个rpm包，如果使用非root用户安装，可以下载libxml2和libxslt的源代码进行安装。 libxml2-devel、libxslt-devel装好后，解压lxml的包，切换到这个包的路径。
加入CFLAGS进行编译和安装，在shell下依次输入如下命令： CFLAGS=-I/usr/include/libxml2:/usr/include/libxslt/ ......
python的变参
*args和**dargs是Python的两个可
变参数，两者有所不同的是*args是个tuple，**dargs是个dict。
*args
和**dargs并用时，*args必须放在**dargs的前面。
例如：
def func(a,b, *c):
pass
函数func至少有两个参数变参数放在tuple c中
def func(*c): 或者 def func(**d ......
Install Python Eric IDE
1 Download following things
1) Python3.1
2) PyQt for python 3.1
(http://www.riverbankcomputing.co.uk/software/pyqt/download) I am using
PyQt-Py3.1-gpl-4.7.3-2.exe
3) Eric5 IDE
(http://eric-ide.python-projects.org/eric-download.html)
2 ......
