2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6 % CCCC OOO M M PPPP AAA RRRR EEEEE %
7 % C O O MM MM P P A A R R E %
8 % C O O M M M PPPP AAAAA RRRR EEE %
9 % C O O M M P A A R R E %
10 % CCCC OOO M M P A A R R EEEEE %
13 % Image Comparison Methods %
20 % Copyright 1999-2013 ImageMagick Studio LLC, a non-profit organization %
21 % dedicated to making software imaging solutions freely available. %
23 % You may not use this file except in compliance with the License. You may %
24 % obtain a copy of the License at %
26 % http://www.imagemagick.org/script/license.php %
28 % Unless required by applicable law or agreed to in writing, software %
29 % distributed under the License is distributed on an "AS IS" BASIS, %
30 % WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. %
31 % See the License for the specific language governing permissions and %
32 % limitations under the License. %
34 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
36 % Use the compare program to mathematically and visually annotate the
37 % difference between an image and its reconstruction.
44 #include "MagickWand/studio.h"
45 #include "MagickWand/MagickWand.h"
46 #include "MagickWand/mogrify-private.h"
47 #include "MagickCore/string-private.h"
50 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
54 % C o m p a r e I m a g e C o m m a n d %
58 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
60 % CompareImagesCommand() compares two images and returns the difference between
61 % them as a distortion metric and as a new image visually annotating their
64 % The format of the CompareImagesCommand method is:
66 % MagickBooleanType CompareImagesCommand(ImageInfo *image_info,int argc,
67 % char **argv,char **metadata,ExceptionInfo *exception)
69 % A description of each parameter follows:
71 % o image_info: the image info.
73 % o argc: the number of elements in the argument vector.
75 % o argv: A text array containing the command line arguments.
77 % o metadata: any metadata is returned here.
79 % o exception: return any errors or warnings in this structure.
83 static MagickBooleanType CompareUsage(void)
91 "-debug events display copious debugging information",
92 "-help print program options",
93 "-list type print a list of supported option arguments",
94 "-log format format of debugging information",
99 "-alpha option on, activate, off, deactivate, set, opaque, copy",
100 " transparent, extract, background, or shape",
101 "-authenticate password",
102 " decipher image with this password",
103 "-channel type apply option to select image channels",
104 "-colorspace type alternate image colorspace",
105 "-compose operator set image composite operator",
106 "-compress type type of pixel compression when writing the image",
107 "-decipher filename convert cipher pixels to plain pixels",
108 "-define format:option",
109 " define one or more image format options",
110 "-density geometry horizontal and vertical density of the image",
111 "-depth value image depth",
112 "-dissimilarity-threshold value",
113 " maximum distortion for (sub)image match",
114 "-encipher filename convert plain pixels to cipher pixels",
115 "-extract geometry extract area from image",
116 "-format \"string\" output formatted image characteristics",
117 "-fuzz distance colors within this distance are considered equal",
118 "-highlight-color color",
119 " empasize pixel differences with this color",
120 "-identify identify the format and characteristics of the image",
121 "-interlace type type of image interlacing scheme",
122 "-limit type value pixel cache resource limit",
123 "-lowlight-color color",
124 " de-emphasize pixel differences with this color",
125 "-metric type measure differences between images with this metric",
126 "-monitor monitor progress",
127 "-profile filename add, delete, or apply an image profile",
128 "-quality value JPEG/MIFF/PNG compression level",
129 "-quiet suppress all warning messages",
130 "-quantize colorspace reduce colors in this colorspace",
131 "-regard-warnings pay attention to warning messages",
132 "-respect-parentheses settings remain in effect until parenthesis boundary",
133 "-sampling-factor geometry",
134 " horizontal and vertical sampling factor",
135 "-seed value seed a new sequence of pseudo-random numbers",
136 "-set attribute value set an image attribute",
137 "-quality value JPEG/MIFF/PNG compression level",
138 "-similarity-threshold value",
139 " minimum distortion for (sub)image match",
140 "-size geometry width and height of image",
141 "-subimage-search search for subimage",
142 "-synchronize synchronize image to storage device",
143 "-taint declare the image as modified",
144 "-transparent-color color",
145 " transparent color",
146 "-type type image type",
147 "-verbose print detailed information about the image",
148 "-version print version information",
149 "-virtual-pixel method",
150 " virtual pixel access method",
154 ListMagickVersion(stdout);
155 (void) printf("Usage: %s [options ...] image reconstruct difference\n",
157 (void) printf("\nImage Settings:\n");
158 for (p=settings; *p != (char *) NULL; p++)
159 (void) printf(" %s\n",*p);
160 (void) printf("\nMiscellaneous Options:\n");
161 for (p=miscellaneous; *p != (char *) NULL; p++)
162 (void) printf(" %s\n",*p);
164 "\nBy default, the image format of 'file' is determined by its magic\n");
166 "number. To specify a particular image format, precede the filename\n");
168 "with an image format name and a colon (i.e. ps:image) or specify the\n");
170 "image type as the filename suffix (i.e. image.ps). Specify 'file' as\n");
171 (void) printf("'-' for standard input or output.\n");
175 WandExport MagickBooleanType CompareImagesCommand(ImageInfo *image_info,
176 int argc,char **argv,char **metadata,ExceptionInfo *exception)
178 #define CompareEpsilon (1.0e-06)
179 #define DefaultDissimilarityThreshold 0.31830988618379067154
180 #define DefaultSimilarityThreshold (-1.0)
181 #define DestroyCompare() \
183 if (similarity_image != (Image *) NULL) \
184 similarity_image=DestroyImageList(similarity_image); \
185 if (difference_image != (Image *) NULL) \
186 difference_image=DestroyImageList(difference_image); \
187 DestroyImageStack(); \
188 for (i=0; i < (ssize_t) argc; i++) \
189 argv[i]=DestroyString(argv[i]); \
190 argv=(char **) RelinquishMagickMemory(argv); \
192 #define ThrowCompareException(asperity,tag,option) \
194 if (exception->severity < (asperity)) \
195 (void) ThrowMagickException(exception,GetMagickModule(),asperity,tag, \
198 return(MagickFalse); \
200 #define ThrowCompareInvalidArgumentException(option,argument) \
202 (void) ThrowMagickException(exception,GetMagickModule(),OptionError, \
203 "InvalidArgument","'%s': %s",option,argument); \
205 return(MagickFalse); \
216 dissimilarity_threshold,
219 similarity_threshold;
228 image_stack[MaxImageStackDepth+1];
255 assert(image_info != (ImageInfo *) NULL);
256 assert(image_info->signature == MagickSignature);
257 if (image_info->debug != MagickFalse)
258 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
259 assert(exception != (ExceptionInfo *) NULL);
263 if ((LocaleCompare("version",option+1) == 0) ||
264 (LocaleCompare("-version",option+1) == 0))
266 ListMagickVersion(stdout);
271 return(CompareUsage());
272 difference_image=NewImageList();
273 similarity_image=NewImageList();
274 dissimilarity_threshold=DefaultDissimilarityThreshold;
275 similarity_threshold=DefaultSimilarityThreshold;
277 format=(char *) NULL;
280 metric=UndefinedErrorMetric;
282 option=(char *) NULL;
284 reconstruct_image=NewImageList();
285 respect_parenthesis=MagickFalse;
287 subimage_search=MagickFalse;
291 ReadCommandlLine(argc,&argv);
292 status=ExpandFilenames(&argc,&argv);
293 if (status == MagickFalse)
294 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
295 GetExceptionMessage(errno));
296 for (i=1; i < (ssize_t) (argc-1); i++)
299 if (LocaleCompare(option,"(") == 0)
301 FireImageStack(MagickTrue,MagickTrue,pend);
302 if (k == MaxImageStackDepth)
303 ThrowCompareException(OptionError,"ParenthesisNestedTooDeeply",
308 if (LocaleCompare(option,")") == 0)
310 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
312 ThrowCompareException(OptionError,"UnableToParseExpression",option);
316 if (IsCommandOption(option) == MagickFalse)
324 FireImageStack(MagickFalse,MagickFalse,pend);
326 if ((LocaleCompare(filename,"--") == 0) && (i < (ssize_t) (argc-1)))
328 images=ReadImages(image_info,filename,exception);
329 status&=(images != (Image *) NULL) &&
330 (exception->severity < ErrorException);
331 if (images == (Image *) NULL)
333 AppendImageStack(images);
336 pend=image != (Image *) NULL ? MagickTrue : MagickFalse;
341 if (LocaleCompare("alpha",option+1) == 0)
349 if (i == (ssize_t) argc)
350 ThrowCompareException(OptionError,"MissingArgument",option);
351 type=ParseCommandOption(MagickAlphaChannelOptions,MagickFalse,argv[i]);
353 ThrowCompareException(OptionError,"UnrecognizedAlphaChannelOption",
357 if (LocaleCompare("authenticate",option+1) == 0)
362 if (i == (ssize_t) argc)
363 ThrowCompareException(OptionError,"MissingArgument",option);
366 ThrowCompareException(OptionError,"UnrecognizedOption",option);
370 if (LocaleCompare("cache",option+1) == 0)
375 if (i == (ssize_t) argc)
376 ThrowCompareException(OptionError,"MissingArgument",option);
377 if (IsGeometry(argv[i]) == MagickFalse)
378 ThrowCompareInvalidArgumentException(option,argv[i]);
381 if (LocaleCompare("channel",option+1) == 0)
389 if (i == (ssize_t) (argc-1))
390 ThrowCompareException(OptionError,"MissingArgument",option);
391 channel=ParseChannelOption(argv[i]);
393 ThrowCompareException(OptionError,"UnrecognizedChannelType",
395 SetPixelChannelMask(image,(ChannelType) channel);
398 if (LocaleCompare("colorspace",option+1) == 0)
406 if (i == (ssize_t) (argc-1))
407 ThrowCompareException(OptionError,"MissingArgument",option);
408 colorspace=ParseCommandOption(MagickColorspaceOptions,MagickFalse,
411 ThrowCompareException(OptionError,"UnrecognizedColorspace",
415 if (LocaleCompare("compose",option+1) == 0)
423 if (i == (ssize_t) argc)
424 ThrowCompareException(OptionError,"MissingArgument",option);
425 compose=ParseCommandOption(MagickComposeOptions,MagickFalse,
428 ThrowCompareException(OptionError,"UnrecognizedComposeOperator",
432 if (LocaleCompare("compress",option+1) == 0)
440 if (i == (ssize_t) (argc-1))
441 ThrowCompareException(OptionError,"MissingArgument",option);
442 compress=ParseCommandOption(MagickCompressOptions,MagickFalse,
445 ThrowCompareException(OptionError,"UnrecognizedImageCompression",
449 if (LocaleCompare("concurrent",option+1) == 0)
451 ThrowCompareException(OptionError,"UnrecognizedOption",option)
455 if (LocaleCompare("debug",option+1) == 0)
463 if (i == (ssize_t) argc)
464 ThrowCompareException(OptionError,"MissingArgument",option);
465 event_mask=SetLogEventMask(argv[i]);
466 if (event_mask == UndefinedEvents)
467 ThrowCompareException(OptionError,"UnrecognizedEventType",
471 if (LocaleCompare("decipher",option+1) == 0)
476 if (i == (ssize_t) (argc-1))
477 ThrowCompareException(OptionError,"MissingArgument",option);
480 if (LocaleCompare("define",option+1) == 0)
483 if (i == (ssize_t) argc)
484 ThrowCompareException(OptionError,"MissingArgument",option);
490 define=GetImageOption(image_info,argv[i]);
491 if (define == (const char *) NULL)
492 ThrowCompareException(OptionError,"NoSuchOption",argv[i]);
497 if (LocaleCompare("density",option+1) == 0)
502 if (i == (ssize_t) argc)
503 ThrowCompareException(OptionError,"MissingArgument",option);
504 if (IsGeometry(argv[i]) == MagickFalse)
505 ThrowCompareInvalidArgumentException(option,argv[i]);
508 if (LocaleCompare("depth",option+1) == 0)
513 if (i == (ssize_t) argc)
514 ThrowCompareException(OptionError,"MissingArgument",option);
515 if (IsGeometry(argv[i]) == MagickFalse)
516 ThrowCompareInvalidArgumentException(option,argv[i]);
519 if (LocaleCompare("dissimilarity-threshold",option+1) == 0)
524 if (i == (ssize_t) argc)
525 ThrowCompareException(OptionError,"MissingArgument",option);
526 if (IsGeometry(argv[i]) == MagickFalse)
527 ThrowCompareInvalidArgumentException(option,argv[i]);
529 dissimilarity_threshold=DefaultDissimilarityThreshold;
531 dissimilarity_threshold=StringToDouble(argv[i],(char **) NULL);
534 if (LocaleCompare("duration",option+1) == 0)
539 if (i == (ssize_t) (argc-1))
540 ThrowCompareException(OptionError,"MissingArgument",option);
541 if (IsGeometry(argv[i]) == MagickFalse)
542 ThrowCompareInvalidArgumentException(option,argv[i]);
545 ThrowCompareException(OptionError,"UnrecognizedOption",option)
549 if (LocaleCompare("encipher",option+1) == 0)
554 if (i == (ssize_t) (argc-1))
555 ThrowCompareException(OptionError,"MissingArgument",option);
558 if (LocaleCompare("extract",option+1) == 0)
563 if (i == (ssize_t) (argc-1))
564 ThrowCompareException(OptionError,"MissingArgument",option);
565 if (IsGeometry(argv[i]) == MagickFalse)
566 ThrowCompareInvalidArgumentException(option,argv[i]);
569 ThrowCompareException(OptionError,"UnrecognizedOption",option)
573 if (LocaleCompare("format",option+1) == 0)
578 if (i == (ssize_t) argc)
579 ThrowCompareException(OptionError,"MissingArgument",option);
583 if (LocaleCompare("fuzz",option+1) == 0)
588 if (i == (ssize_t) (argc-1))
589 ThrowCompareException(OptionError,"MissingArgument",option);
590 if (IsGeometry(argv[i]) == MagickFalse)
591 ThrowCompareInvalidArgumentException(option,argv[i]);
594 ThrowCompareException(OptionError,"UnrecognizedOption",option)
598 if ((LocaleCompare("help",option+1) == 0) ||
599 (LocaleCompare("-help",option+1) == 0))
600 return(CompareUsage());
601 if (LocaleCompare("highlight-color",option+1) == 0)
606 if (i == (ssize_t) (argc-1))
607 ThrowCompareException(OptionError,"MissingArgument",option);
610 ThrowCompareException(OptionError,"UnrecognizedOption",option)
614 if (LocaleCompare("identify",option+1) == 0)
616 if (LocaleCompare("interlace",option+1) == 0)
624 if (i == (ssize_t) argc)
625 ThrowCompareException(OptionError,"MissingArgument",option);
626 interlace=ParseCommandOption(MagickInterlaceOptions,MagickFalse,
629 ThrowCompareException(OptionError,"UnrecognizedInterlaceType",
633 ThrowCompareException(OptionError,"UnrecognizedOption",option)
637 if (LocaleCompare("limit",option+1) == 0)
651 if (i == (ssize_t) argc)
652 ThrowCompareException(OptionError,"MissingArgument",option);
653 resource=ParseCommandOption(MagickResourceOptions,MagickFalse,
656 ThrowCompareException(OptionError,"UnrecognizedResourceType",
659 if (i == (ssize_t) argc)
660 ThrowCompareException(OptionError,"MissingArgument",option);
661 value=StringToDouble(argv[i],&p);
663 if ((p == argv[i]) && (LocaleCompare("unlimited",argv[i]) != 0))
664 ThrowCompareInvalidArgumentException(option,argv[i]);
667 if (LocaleCompare("list",option+1) == 0)
675 if (i == (ssize_t) argc)
676 ThrowCompareException(OptionError,"MissingArgument",option);
677 list=ParseCommandOption(MagickListOptions,MagickFalse,argv[i]);
679 ThrowCompareException(OptionError,"UnrecognizedListType",argv[i]);
680 status=MogrifyImageInfo(image_info,(int) (i-j+1),(const char **)
683 return(status != 0 ? MagickFalse : MagickTrue);
685 if (LocaleCompare("log",option+1) == 0)
690 if ((i == (ssize_t) argc) || (strchr(argv[i],'%') == (char *) NULL))
691 ThrowCompareException(OptionError,"MissingArgument",option);
694 if (LocaleCompare("lowlight-color",option+1) == 0)
699 if (i == (ssize_t) (argc-1))
700 ThrowCompareException(OptionError,"MissingArgument",option);
703 ThrowCompareException(OptionError,"UnrecognizedOption",option)
707 if (LocaleCompare("matte",option+1) == 0)
709 if (LocaleCompare("metric",option+1) == 0)
717 if (i == (ssize_t) argc)
718 ThrowCompareException(OptionError,"MissingArgument",option);
719 type=ParseCommandOption(MagickMetricOptions,MagickTrue,argv[i]);
721 ThrowCompareException(OptionError,"UnrecognizedMetricType",
723 metric=(MetricType) type;
726 if (LocaleCompare("monitor",option+1) == 0)
728 ThrowCompareException(OptionError,"UnrecognizedOption",option)
732 if (LocaleCompare("profile",option+1) == 0)
735 if (i == (ssize_t) (argc-1))
736 ThrowCompareException(OptionError,"MissingArgument",option);
739 ThrowCompareException(OptionError,"UnrecognizedOption",option)
743 if (LocaleCompare("quality",option+1) == 0)
748 if (i == (ssize_t) (argc-1))
749 ThrowCompareException(OptionError,"MissingArgument",option);
750 if (IsGeometry(argv[i]) == MagickFalse)
751 ThrowCompareInvalidArgumentException(option,argv[i]);
754 if (LocaleCompare("quantize",option+1) == 0)
762 if (i == (ssize_t) (argc-1))
763 ThrowCompareException(OptionError,"MissingArgument",option);
764 colorspace=ParseCommandOption(MagickColorspaceOptions,
765 MagickFalse,argv[i]);
767 ThrowCompareException(OptionError,"UnrecognizedColorspace",
771 if (LocaleCompare("quiet",option+1) == 0)
773 ThrowCompareException(OptionError,"UnrecognizedOption",option)
777 if (LocaleCompare("regard-warnings",option+1) == 0)
779 if (LocaleNCompare("respect-parentheses",option+1,17) == 0)
781 respect_parenthesis=(*option == '-') ? MagickTrue : MagickFalse;
784 ThrowCompareException(OptionError,"UnrecognizedOption",option)
788 if (LocaleCompare("sampling-factor",option+1) == 0)
793 if (i == (ssize_t) argc)
794 ThrowCompareException(OptionError,"MissingArgument",option);
795 if (IsGeometry(argv[i]) == MagickFalse)
796 ThrowCompareInvalidArgumentException(option,argv[i]);
799 if (LocaleCompare("seed",option+1) == 0)
804 if (i == (ssize_t) (argc-1))
805 ThrowCompareException(OptionError,"MissingArgument",option);
806 if (IsGeometry(argv[i]) == MagickFalse)
807 ThrowCompareInvalidArgumentException(option,argv[i]);
810 if (LocaleCompare("set",option+1) == 0)
813 if (i == (ssize_t) argc)
814 ThrowCompareException(OptionError,"MissingArgument",option);
818 if (i == (ssize_t) argc)
819 ThrowCompareException(OptionError,"MissingArgument",option);
822 if (LocaleCompare("similarity-threshold",option+1) == 0)
827 if (i == (ssize_t) argc)
828 ThrowCompareException(OptionError,"MissingArgument",option);
829 if (IsGeometry(argv[i]) == MagickFalse)
830 ThrowCompareInvalidArgumentException(option,argv[i]);
832 similarity_threshold=DefaultSimilarityThreshold;
834 similarity_threshold=StringToDouble(argv[i],(char **) NULL);
837 if (LocaleCompare("size",option+1) == 0)
842 if (i == (ssize_t) argc)
843 ThrowCompareException(OptionError,"MissingArgument",option);
844 if (IsGeometry(argv[i]) == MagickFalse)
845 ThrowCompareInvalidArgumentException(option,argv[i]);
848 if (LocaleCompare("subimage-search",option+1) == 0)
852 subimage_search=MagickFalse;
855 subimage_search=MagickTrue;
858 if (LocaleCompare("synchronize",option+1) == 0)
860 ThrowCompareException(OptionError,"UnrecognizedOption",option)
864 if (LocaleCompare("taint",option+1) == 0)
866 if (LocaleCompare("transparent-color",option+1) == 0)
871 if (i == (ssize_t) (argc-1))
872 ThrowCompareException(OptionError,"MissingArgument",option);
875 if (LocaleCompare("type",option+1) == 0)
883 if (i == (ssize_t) argc)
884 ThrowCompareException(OptionError,"MissingArgument",option);
885 type=ParseCommandOption(MagickTypeOptions,MagickFalse,argv[i]);
887 ThrowCompareException(OptionError,"UnrecognizedImageType",
891 ThrowCompareException(OptionError,"UnrecognizedOption",option)
895 if (LocaleCompare("verbose",option+1) == 0)
897 if ((LocaleCompare("version",option+1) == 0) ||
898 (LocaleCompare("-version",option+1) == 0))
900 ListMagickVersion(stdout);
903 if (LocaleCompare("virtual-pixel",option+1) == 0)
911 if (i == (ssize_t) (argc-1))
912 ThrowCompareException(OptionError,"MissingArgument",option);
913 method=ParseCommandOption(MagickVirtualPixelOptions,MagickFalse,
916 ThrowCompareException(OptionError,
917 "UnrecognizedVirtualPixelMethod",argv[i]);
920 ThrowCompareException(OptionError,"UnrecognizedOption",option)
925 ThrowCompareException(OptionError,"UnrecognizedOption",option)
927 fire=(GetCommandOptionFlags(MagickCommandOptions,MagickFalse,option) &
928 FireOptionFlag) == 0 ? MagickFalse : MagickTrue;
929 if (fire != MagickFalse)
930 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
933 ThrowCompareException(OptionError,"UnbalancedParenthesis",argv[i]);
934 if (i-- != (ssize_t) (argc-1))
935 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
936 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
937 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
938 FinalizeImageSettings(image_info,image,MagickTrue);
939 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
940 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
941 image=GetImageFromList(image,0);
942 reconstruct_image=GetImageFromList(image,1);
943 if (subimage_search != MagickFalse)
945 similarity_image=SimilarityImage(image,reconstruct_image,metric,
946 similarity_threshold,&offset,&similarity_metric,exception);
947 if (similarity_metric > dissimilarity_threshold)
948 ThrowCompareException(ImageError,"ImagesTooDissimilar",image->filename);
950 if ((reconstruct_image->columns == image->columns) &&
951 (reconstruct_image->rows == image->rows))
952 difference_image=CompareImages(image,reconstruct_image,metric,&distortion,
955 if (similarity_image == (Image *) NULL)
956 ThrowCompareException(OptionError,"ImageWidthsOrHeightsDiffer",
964 Determine if reconstructed image is a subimage of the image.
966 composite_image=CloneImage(image,0,0,MagickTrue,exception);
967 if (composite_image == (Image *) NULL)
968 difference_image=CompareImages(image,reconstruct_image,metric,
969 &distortion,exception);
978 (void) CompositeImage(composite_image,reconstruct_image,
979 CopyCompositeOp,MagickTrue,offset.x,offset.y,exception);
980 difference_image=CompareImages(image,composite_image,metric,
981 &distortion,exception);
982 if (difference_image != (Image *) NULL)
984 difference_image->page.x=offset.x;
985 difference_image->page.y=offset.y;
987 composite_image=DestroyImage(composite_image);
988 page.width=reconstruct_image->columns;
989 page.height=reconstruct_image->rows;
992 distort_image=CropImage(image,&page,exception);
993 if (distort_image != (Image *) NULL)
998 sans_image=CompareImages(distort_image,reconstruct_image,metric,
999 &distortion,exception);
1000 distort_image=DestroyImage(distort_image);
1001 if (sans_image != (Image *) NULL)
1002 sans_image=DestroyImage(sans_image);
1005 if (difference_image != (Image *) NULL)
1007 AppendImageToList(&difference_image,similarity_image);
1008 similarity_image=(Image *) NULL;
1011 if (difference_image == (Image *) NULL)
1015 if (image_info->verbose != MagickFalse)
1016 (void) IsImagesEqual(image,reconstruct_image,exception);
1017 if (*difference_image->magick == '\0')
1018 (void) CopyMagickString(difference_image->magick,image->magick,
1020 if (image_info->verbose == MagickFalse)
1024 case FuzzErrorMetric:
1025 case MeanAbsoluteErrorMetric:
1026 case MeanSquaredErrorMetric:
1027 case RootMeanSquaredErrorMetric:
1028 case PeakAbsoluteErrorMetric:
1030 (void) FormatLocaleFile(stderr,"%g (%g)",QuantumRange*distortion,
1031 (double) distortion);
1032 if ((reconstruct_image->columns != image->columns) ||
1033 (reconstruct_image->rows != image->rows))
1034 (void) FormatLocaleFile(stderr," @ %.20g,%.20g",(double)
1035 difference_image->page.x,(double) difference_image->page.y);
1038 case AbsoluteErrorMetric:
1039 case NormalizedCrossCorrelationErrorMetric:
1040 case PeakSignalToNoiseRatioMetric:
1042 (void) FormatLocaleFile(stderr,"%g",distortion);
1043 if ((reconstruct_image->columns != image->columns) ||
1044 (reconstruct_image->rows != image->rows))
1045 (void) FormatLocaleFile(stderr," @ %.20g,%.20g",(double)
1046 difference_image->page.x,(double) difference_image->page.y);
1049 case MeanErrorPerPixelMetric:
1051 (void) FormatLocaleFile(stderr,"%g (%g, %g)",distortion,
1052 image->error.normalized_mean_error,
1053 image->error.normalized_maximum_error);
1054 if ((reconstruct_image->columns != image->columns) ||
1055 (reconstruct_image->rows != image->rows))
1056 (void) FormatLocaleFile(stderr," @ %.20g,%.20g",(double)
1057 difference_image->page.x,(double) difference_image->page.y);
1060 case UndefinedErrorMetric:
1067 *channel_distortion;
1069 channel_distortion=GetImageDistortions(image,reconstruct_image,
1071 (void) FormatLocaleFile(stderr,"Image: %s\n",image->filename);
1072 if ((reconstruct_image->columns != image->columns) ||
1073 (reconstruct_image->rows != image->rows))
1074 (void) FormatLocaleFile(stderr,"Offset: %.20g,%.20g\n",(double)
1075 difference_image->page.x,(double) difference_image->page.y);
1076 (void) FormatLocaleFile(stderr," Channel distortion: %s\n",
1077 CommandOptionToMnemonic(MagickMetricOptions,(ssize_t) metric));
1080 case FuzzErrorMetric:
1081 case MeanAbsoluteErrorMetric:
1082 case MeanSquaredErrorMetric:
1083 case RootMeanSquaredErrorMetric:
1084 case PeakAbsoluteErrorMetric:
1086 switch (image->colorspace)
1091 (void) FormatLocaleFile(stderr," red: %g (%g)\n",
1092 QuantumRange*channel_distortion[RedPixelChannel],
1093 channel_distortion[RedPixelChannel]);
1094 (void) FormatLocaleFile(stderr," green: %g (%g)\n",
1095 QuantumRange*channel_distortion[GreenPixelChannel],
1096 channel_distortion[GreenPixelChannel]);
1097 (void) FormatLocaleFile(stderr," blue: %g (%g)\n",
1098 QuantumRange*channel_distortion[BluePixelChannel],
1099 channel_distortion[BluePixelChannel]);
1100 if (image->alpha_trait == BlendPixelTrait)
1101 (void) FormatLocaleFile(stderr," alpha: %g (%g)\n",
1102 QuantumRange*channel_distortion[AlphaPixelChannel],
1103 channel_distortion[AlphaPixelChannel]);
1106 case CMYKColorspace:
1108 (void) FormatLocaleFile(stderr," cyan: %g (%g)\n",
1109 QuantumRange*channel_distortion[CyanPixelChannel],
1110 channel_distortion[CyanPixelChannel]);
1111 (void) FormatLocaleFile(stderr," magenta: %g (%g)\n",
1112 QuantumRange*channel_distortion[MagentaPixelChannel],
1113 channel_distortion[MagentaPixelChannel]);
1114 (void) FormatLocaleFile(stderr," yellow: %g (%g)\n",
1115 QuantumRange*channel_distortion[YellowPixelChannel],
1116 channel_distortion[YellowPixelChannel]);
1117 (void) FormatLocaleFile(stderr," black: %g (%g)\n",
1118 QuantumRange*channel_distortion[BlackPixelChannel],
1119 channel_distortion[BlackPixelChannel]);
1120 if (image->alpha_trait == BlendPixelTrait)
1121 (void) FormatLocaleFile(stderr," alpha: %g (%g)\n",
1122 QuantumRange*channel_distortion[AlphaPixelChannel],
1123 channel_distortion[AlphaPixelChannel]);
1126 case GRAYColorspace:
1128 (void) FormatLocaleFile(stderr," gray: %g (%g)\n",
1129 QuantumRange*channel_distortion[GrayPixelChannel],
1130 channel_distortion[GrayPixelChannel]);
1131 if (image->alpha_trait == BlendPixelTrait)
1132 (void) FormatLocaleFile(stderr," alpha: %g (%g)\n",
1133 QuantumRange*channel_distortion[AlphaPixelChannel],
1134 channel_distortion[AlphaPixelChannel]);
1138 (void) FormatLocaleFile(stderr," all: %g (%g)\n",
1139 QuantumRange*channel_distortion[MaxPixelChannels],
1140 channel_distortion[MaxPixelChannels]);
1143 case AbsoluteErrorMetric:
1144 case NormalizedCrossCorrelationErrorMetric:
1145 case PeakSignalToNoiseRatioMetric:
1147 switch (image->colorspace)
1152 (void) FormatLocaleFile(stderr," red: %g\n",
1153 channel_distortion[RedPixelChannel]);
1154 (void) FormatLocaleFile(stderr," green: %g\n",
1155 channel_distortion[GreenPixelChannel]);
1156 (void) FormatLocaleFile(stderr," blue: %g\n",
1157 channel_distortion[BluePixelChannel]);
1158 if (image->alpha_trait == BlendPixelTrait)
1159 (void) FormatLocaleFile(stderr," alpha: %g\n",
1160 channel_distortion[AlphaPixelChannel]);
1163 case CMYKColorspace:
1165 (void) FormatLocaleFile(stderr," cyan: %g\n",
1166 channel_distortion[CyanPixelChannel]);
1167 (void) FormatLocaleFile(stderr," magenta: %g\n",
1168 channel_distortion[MagentaPixelChannel]);
1169 (void) FormatLocaleFile(stderr," yellow: %g\n",
1170 channel_distortion[YellowPixelChannel]);
1171 (void) FormatLocaleFile(stderr," black: %g\n",
1172 channel_distortion[BlackPixelChannel]);
1173 if (image->alpha_trait == BlendPixelTrait)
1174 (void) FormatLocaleFile(stderr," alpha: %g\n",
1175 channel_distortion[AlphaPixelChannel]);
1178 case GRAYColorspace:
1180 (void) FormatLocaleFile(stderr," gray: %g\n",
1181 channel_distortion[GrayPixelChannel]);
1182 if (image->alpha_trait == BlendPixelTrait)
1183 (void) FormatLocaleFile(stderr," alpha: %g\n",
1184 channel_distortion[AlphaPixelChannel]);
1188 (void) FormatLocaleFile(stderr," all: %g\n",
1189 channel_distortion[MaxPixelChannels]);
1192 case MeanErrorPerPixelMetric:
1194 (void) FormatLocaleFile(stderr," %g (%g, %g)\n",
1195 channel_distortion[MaxPixelChannels],
1196 image->error.normalized_mean_error,
1197 image->error.normalized_maximum_error);
1200 case UndefinedErrorMetric:
1203 channel_distortion=(double *) RelinquishMagickMemory(
1204 channel_distortion);
1206 status&=WriteImages(image_info,difference_image,argv[argc-1],exception);
1207 if ((metadata != (char **) NULL) && (format != (char *) NULL))
1212 text=InterpretImageProperties(image_info,difference_image,format,
1214 if (text == (char *) NULL)
1215 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
1216 GetExceptionMessage(errno));
1217 (void) ConcatenateString(&(*metadata),text);
1218 text=DestroyString(text);
1220 difference_image=DestroyImageList(difference_image);
1223 if ((metric == NormalizedCrossCorrelationErrorMetric) ||
1224 (metric == UndefinedErrorMetric))
1226 if (fabs(distortion-1.0) > CompareEpsilon)
1230 if (fabs(distortion) > CompareEpsilon)
1232 return(status != 0 ? MagickTrue : MagickFalse);