...
Now calculating the ellipticity components:. We have the position angle calculated, and we calculate rotazel in the same way as before (with calculate_rotazel function). Now, we rotate the points before calculating the ellipticity components. Wondering if we can just rotate the angle?