std_profile = np.asarray([np.std(hp_hdu.data[pix]) for pix in a_pix])