Read HDRs to Numpy arrays