double
difference;
+ if ((fabs(reconstruct_moments[channel].I[i]) < MagickEpsilon) ||
+ (fabs(image_moments[channel].I[i]) < MagickEpsilon))
+ continue;
difference=log10(fabs(reconstruct_moments[channel].I[i]))-
log10(fabs(image_moments[channel].I[i]));
distortion[channel]+=difference*difference;
double
difference;
+ if ((fabs(reconstruct_moments[channel].I[i]) < MagickEpsilon) ||
+ (fabs(image_moments[channel].I[i]) < MagickEpsilon))
+ continue;
difference=log10(fabs(reconstruct_moments[channel].I[i]))-
log10(fabs(image_moments[channel].I[i]));
distortion[channel]+=difference*difference;
*/
if (fabs(M00[channel]) < MagickEpsilon)
continue;
- channel_moments[channel].I[0]=M20[channel]+M02[channel];
+ channel_moments[channel].I[0]=M20[channel]+M02[channel];
channel_moments[channel].I[1]=(M20[channel]-M02[channel])*
(M20[channel]-M02[channel])+4.0*M11[channel]*M11[channel];
channel_moments[channel].I[2]=(M30[channel]-3.0*M12[channel])*