读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)
下一篇:思考题(六十五)如何漂洗衣服最省水?

3621

被折叠的 条评论
为什么被折叠?



