人生的意义就是思考其意义分享 http://blog.sciencenet.cn/u/qianlivan 天体物理工作者,家乡云南昆明。

博文

python读hdf5文件

已有 568 次阅读 2020-10-23 14:14 |个人分类:知识|系统分类:科研笔记| python, hdf5

读hdf5文件和读fits差不多,需要知道类似文件头的东西,这里是一个变量的名称。这可以通过定义一个prt函数结合h5py报的工具来实现如下

def prt(name):

    print(name)

 

filename=sys.argv[1]

f = h5py.File(filename, 'r')

f.visit(prt)

 

知道了变量名以后就可以读取相应的数据了。例如

vx=np.array(f['velocity_x'])

 

之后就和通常的操作一样了。总的读取数据画图的脚本如下。

 

import h5py

import sys

import numpy as np

import matplotlib.pyplot as plt

import matplotlib

 

matplotlib.style.use('classic')

 

def prt(name):

    print(name)

 

filename=sys.argv[1]

f = h5py.File(filename, 'r')

f.visit(prt)

vx=np.array(f['velocity_x'])

integrated_vx=np.sum(vx,axis=0)

 

plt.imshow(integrated_vx)

plt.show()




http://blog.sciencenet.cn/blog-117333-1255475.html

上一篇:arXiv周记202010_4(20201019-20201023)
下一篇:思考题(六十五)如何漂洗衣服最省水?

0

评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2020-12-3 04:17

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部