Esta es una leccion que aprendimos gracias a chatgpt. Dice que primero ha que instalar “pip install h5py numpy”
import h5py
import numpy as np
# Abre el archivo HDF5 en modo lectura
with h5py.File('eps-000000.00.h5', 'r') as f:
# Muestra las claves principales del archivo
print(f"Claves principales del archivo: {list(f.keys())}")
Claves principales del archivo: ['eps']
with h5py.File('eps-000000.00.h5', 'r') as f:
dataset = f['eps']
data_array = np.array(dataset)
con eso ya se leen los archgivos, por ejemplo, podemos graficar
import matplotlib.pyplot as plt plt.imshow(data_array)
El campo electromagetico tambien lo podemos visualizar si hacemos
# Abre el archivo HDF5 en modo lectura
with h5py.File('ez.h5', 'r') as f:
# Muestra las claves principales del archivo
print(f"Claves principales del archivo: {list(f.keys())}")
plt.imshow(data3[:,:,332])
# Abre el archivo HDF5 en modo lectura
with h5py.File('ez.h5', 'r') as f:
aux1 = f['ez']
# Convierte el dataset a un arreglo de numpy
aux2 = np.array(aux1)
plt.imshow(data3[:,:,332])
plt.plot(data3[40,40,:])
plt.imshow(data3[40,:,:].transpose())
plt.imshow(data3[:,40,:].transpose())
h5topng -t 0:332 -R -Zc dkbluered -a yarg -A eps-000000.00.h5 ez.h5 convert ez.t*.png ez.gif
