使用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.改变本文件的属性
import
os
import
stat
os.chmod( filename, stat.S_IWRITE )
2.改变本目录及其子目录属性
import
os
os.system(r
'
attrib -r' + path +'\\*.* /s
'
)
3.介绍改变文件属性的dos指令
Attrib
显示、设置或删除指派给文件或目录的只读、存档、系统以及隐藏属性。如果在不含参数的� ......
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 ......
编译了一个windows下的python3连接Mysql的库
mypysql
版本是 0.5.1 ，根据官方的修改日志，这个版本修改了0.5中一个内存泄漏问题。
源代码和编译后的文件为：
http://211.147.215.55/down/mypysql-0.5.1-win.zip
mypysql的官方地址 https://sourceforge.net/projects/mypysql/
......
折腾计算机，被计算机折腾，前前后后也有十几年。
最可悲的是居然一无所成，虽然仅仅出于兴趣，并无工作之压力，但是这种挫折感确实是实打实的。
已经记不起是哪一天了，好像突然就知道了python的存在，从最初的不经意，到越来越喜爱，我一步步的向python靠近。
我喜欢Python，这是一个事实，我撰� ......
译者言：
早在 2008 年 8 月，我就曾在自己的博客发表了一篇《为什么<Dive into Python>不值得推荐》（http://blog.csdn.net/lanphaday/archive/2008/08/28/2845258.aspx
），当时引起的讨论就不多说了，不过事实上到今天仍然有许多朋友留言与我讨论，让我几次贼心不死，想写再深入批判。好运的时候在我真正攒足动� ......
