使用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，下载org.python.pydev.feature-1.5.0.1251989166.zip http://sourceforge.net/projects/pydev/files/
2，安装插件到eclipse
3，重启eclipse
注意：使用1.5.6版本pydev插件，创建python工程会报错，使用1.5.0版本无此问题。 ......
前言：
最近又想学习python,又想去温习一下算法，于是就想出了这个两不误的方法,^_^
堆栈：
使用python的列表结构，详情可以查看help(list)
#Filename: stack.py
shoplist=['apple','mango','carrot','banana']
shoplist.append('rice')
popitem=shoplist[-1]
del shoplist[-1]
print 'the popitem is',popitem
......
队列：
与堆栈类似，通过python的列表类型来实现，参考 help(list)
shoplist=['apple','mango','carrot','banana']
print 'I have',len(shoplist),'items to purchase'
print 'these items are:'
for item in shoplist:
print item,
shoplist.append('rice')
print 'my shopping list is now', shoplist
shoplist. ......
