(M20[channel]-M02[channel])*(M20[channel]-M02[channel]))));
channel_moments[channel].ellipse_angle=RadiansToDegrees(0.5*atan(2.0*
M11[channel]/(M20[channel]-M02[channel]+MagickEpsilon)));
- if (channel_moments[channel].ellipse_angle < 0.0)
+ while (channel_moments[channel].ellipse_angle < 0.0)
channel_moments[channel].ellipse_angle+=90.0;
- if (channel_moments[channel].ellipse_angle >= 90.0)
+ while (channel_moments[channel].ellipse_angle >= 90.0)
channel_moments[channel].ellipse_angle-=90.0;
channel_moments[channel].ellipse_eccentricity=sqrt(1.0-(
channel_moments[channel].ellipse_axis.y/