double
difference;
- if ((fabs(reconstruct_moments[channel].I[i]) < MagickEpsilon) ||
- (fabs(image_moments[channel].I[i]) < MagickEpsilon))
+ if ((fabs(image_moments[channel].ellipse_intensity) < MagickEpsilon) ||
+ (fabs(reconstruct_moments[channel].ellipse_intensity) < MagickEpsilon))
continue;
difference=log10(fabs(reconstruct_moments[channel].I[i]))-
log10(fabs(image_moments[channel].I[i]));
double
difference;
- if ((fabs(reconstruct_moments[channel].I[i]) < MagickEpsilon) ||
- (fabs(image_moments[channel].I[i]) < MagickEpsilon))
+ if ((fabs(image_moments[channel].ellipse_intensity) < MagickEpsilon) ||
+ (fabs(reconstruct_moments[channel].ellipse_intensity) < MagickEpsilon))
continue;
difference=log10(fabs(reconstruct_moments[channel].I[i]))-
- log10(fabs(image_moments[channel].I[i]));
+ log10(fabs(image_moments[channel].I[i]));
distortion[channel]+=difference*difference;
distortion[CompositePixelChannel]+=difference*difference;
}