使用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()
相关文档：
需要先安装libxml2-devel libxslt-devel这两个rpm包，如果使用非root用户安装，可以下载libxml2和libxslt的源代码进行安装。 libxml2-devel、libxslt-devel装好后，解压lxml的包，切换到这个包的路径。
加入CFLAGS进行编译和安装，在shell下依次输入如下命令： CFLAGS=-I/usr/include/libxml2:/usr/include/libxslt/ ......
前言：
最近又想学习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的存在，从最初的不经意，到越来越喜爱，我一步步的向python靠近。
我喜欢Python，这是一个事实，我撰� ......
小白是个微软迷，他的偶像是比尔大叔，原因当然是地球人都知道啦。大二的时候，他的“辍学计划”曾一度得逞，原因是他挂科太多。当然，大三新学期开始的时候，面对公众质疑，小白站在椅子上，像极了《大内密探零零发》里的西门吹雪：“世界首富比不一定有多出色，这只不过是你们这些星斗市民一厢情愿的想法� ......
1. 手动制作python的exe可执行程序
转载---------------
Python没有内建一个编 ......
