(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)));
- while (channel_moments[channel].ellipse_angle < 0.0)
- 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/
(channel_moments[channel].ellipse_axis.x+MagickEpsilon)));