使用Python读取二进制文件的实例讲解
2018-07-09 11:11
901 查看
目标:目标文件为一个float32型存储的二进制文件,按列优先方式存储。本文使用Python读取该二进制文件并使用matplotlib.pyplot相关工具画出图像
工具:Python3, matplotlib,os,struct,numpy
1. 读取二进制文件
首先使用open函数打开文件,打开模式选择二进制读取"rb"。
f = open(filename, "rb")
第二步,需要打开按照行列读取文件,由于是纯二进制文件,内部不含邮任何的数据结构信息,因此我们需要给定二进制数据的行数列数(nx和ny)来确定图像的形状。这里我们的数据类型是float32型的,对应过来是4bytes,使用for循环逐个read4个字节。
for i in range(nx): for j in range(ny): data = f.read(4)
得到的结果如下:
b'5\x9d\x82\xc3' b'\xb1\x04\x10\xc4' b'\xc1\x9eD\xc4' b'a\x86R\xc4' b'\x15\x01=\xc4'
可以看到,读入的数据为十六进制数,这个时候的数据并不能直接用来画图,需要转换为float32型。这里使用struct.unpack来转译
data_float = struct.unpack("f", data)[0]
就可以得到正确的结果了
-261.2281799316406 -576.0733032226562 -786.4805297851562 -842.0996704101562 -756.0169067382812
建立一个numpy数组,将读入的数据分别按列优先的方式放入数组,就完成了图像读入的操作。
将上述操作包装成函数,代码如下:
def xshow(filename, nx, nz):f = open(filename, "rb") pic = np.zeros((nx, nz)) for i in range(nx): for j in range(nz): data = f.read(4) elem = struct.unpack("f", data)[0] pic[i][j] = elem f.close() return pic
2. 画出图像
这里我们的目的是将矩阵输出灰度图,并保存为tiff格式的图像,效果图如下(该图为使用弹性波波动方程在Marmousi模型下的波场快照):
不需要坐标轴的信息,直接使用imsave命令,使用cmap选定颜色格式即可。
plt.imsave('output.tiff', img, format='tiff', cmap=plt.cm.gray)
如果使用PIL库读取图像,可能会得到不同的效果,可能需要做直方图均衡,目前还没有找到原因。
相关文章推荐
- python使用xmlrpc实例讲解
- 二种python发送邮件实例讲解(python发邮件附件可以使用email模块实现)
- 实例讲解Python的函数闭包使用中应注意的问题
- 实例讲解Python设计模式编程之工厂方法模式的使用
- Python使用pandas处理CSV文件的实例讲解
- 使用C#实现读取系统配置文件的代码实例讲解
- python使用xmlrpc实例讲解
- python django使用haystack:全文检索的框架(实例讲解)
- python队列通信:rabbitMQ的使用(实例讲解)
- Python的字符串模板(Template)使用操作实例讲解
- 【python】如何批量读取文件夹的所有文件数据,os模块与open函数结合使用实例
- Python使用urllib2获取网络资源实例讲解
- python爬虫 正则表达式使用技巧及爬取个人博客的实例讲解
- Python读取和处理文件后缀为.sqlite的数据文件(实例讲解)
- 实例讲解Python设计模式编程之工厂方法模式的使用
- python中文分词,使用结巴分词对python进行分词(实例讲解)
- Python读取图片EXIF信息类库介绍和使用实例
- python使用多线程实例讲解
- Python读取图片EXIF信息类库介绍和使用实例
- 实例讲解Python的函数闭包使用中应注意的问题