One issue with this is that sh2peaks generates values between -1.5 and 1.5 so if you do -acos you get quite a few NaN values.
Could you write out how we calculate the angle (in degrees) between the z-axis and peak orientation from the sh2peaks output?