MagickWand 7.1.1
Convert, Edit, Or Compose Bitmap Images
Loading...
Searching...
No Matches
compare.c
1/*
2%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3% %
4% %
5% %
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 %
11% %
12% %
13% Image Comparison Methods %
14% %
15% Software Design %
16% Cristy %
17% December 2003 %
18% %
19% %
20% Copyright @ 1999 ImageMagick Studio LLC, a non-profit organization %
21% dedicated to making software imaging solutions freely available. %
22% %
23% You may not use this file except in compliance with the License. You may %
24% obtain a copy of the License at %
25% %
26% https://imagemagick.org/script/license.php %
27% %
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. %
33% %
34%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
35%
36% Use the compare program to mathematically and visually annotate the
37% difference between an image and its reconstruction.
38%
39*/
40
41/*
42 Include declarations.
43*/
44#include "MagickWand/studio.h"
45#include "MagickWand/MagickWand.h"
46#include "MagickWand/mogrify-private.h"
47#include "MagickCore/image-private.h"
48#include "MagickCore/string-private.h"
49
50/*
51%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
52% %
53% %
54% %
55% C o m p a r e I m a g e C o m m a n d %
56% %
57% %
58% %
59%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
60%
61% CompareImagesCommand() compares two images and returns the difference between
62% them as a distortion metric and as a new image visually annotating their
63% differences.
64%
65% The format of the CompareImagesCommand method is:
66%
67% MagickBooleanType CompareImagesCommand(ImageInfo *image_info,int argc,
68% char **argv,char **metadata,ExceptionInfo *exception)
69%
70% A description of each parameter follows:
71%
72% o image_info: the image info.
73%
74% o argc: the number of elements in the argument vector.
75%
76% o argv: A text array containing the command line arguments.
77%
78% o metadata: any metadata is returned here.
79%
80% o exception: return any errors or warnings in this structure.
81%
82*/
83
84static MagickBooleanType CompareUsage(void)
85{
86 static const char
87 channel_operators[] =
88 " -separate separate an image channel into a grayscale image",
89 miscellaneous[] =
90 " -channel mask set the image channel mask\n"
91 " -debug events display copious debugging information\n"
92 " -help print program options\n"
93 " -list type print a list of supported option arguments\n"
94 " -log format format of debugging information",
95 operators[] =
96 " -auto-orient automagically orient (rotate) image\n"
97 " -brightness-contrast geometry\n"
98 " improve brightness / contrast of the image\n"
99 " -distort method args\n"
100 " distort images according to given method and args\n"
101 " -level value adjust the level of image contrast\n"
102 " -resize geometry resize the image\n"
103 " -rotate degrees apply Paeth rotation to the image\n"
104 " -sigmoidal-contrast geometry\n"
105 " increase the contrast without saturating highlights or\n"
106 " -trim trim image edges\n"
107 " -write filename write images to this file",
108 sequence_operators[] =
109 " -crop geometry cut out a rectangular region of the image",
110 settings[] =
111 " -adjoin join images into a single multi-image file\n"
112 " -alpha option on, activate, off, deactivate, set, opaque, copy\n"
113 " transparent, extract, background, or shape\n"
114 " -authenticate password\n"
115 " decipher image with this password\n"
116 " -background color background color\n"
117 " -colorspace type alternate image colorspace\n"
118 " -compose operator set image composite operator\n"
119 " -compress type type of pixel compression when writing the image\n"
120 " -decipher filename convert cipher pixels to plain pixels\n"
121 " -define format:option\n"
122 " define one or more image format options\n"
123 " -density geometry horizontal and vertical density of the image\n"
124 " -depth value image depth\n"
125 " -dissimilarity-threshold value\n"
126 " maximum distortion for (sub)image match\n"
127 " -encipher filename convert plain pixels to cipher pixels\n"
128 " -extract geometry extract area from image\n"
129 " -format \"string\" output formatted image characteristics\n"
130 " -fuzz distance colors within this distance are considered equal\n"
131 " -gravity type horizontal and vertical text placement\n"
132 " -highlight-color color\n"
133 " emphasize pixel differences with this color\n"
134 " -identify identify the format and characteristics of the image\n"
135 " -interlace type type of image interlacing scheme\n"
136 " -limit type value pixel cache resource limit\n"
137 " -lowlight-color color\n"
138 " de-emphasize pixel differences with this color\n"
139 " -metric type measure differences between images with this metric\n"
140 " -monitor monitor progress\n"
141 " -negate replace every pixel with its complementary color \n"
142 " -passphrase filename get the passphrase from this file\n"
143 " -precision value maximum number of significant digits to print\n"
144 " -profile filename add, delete, or apply an image profile\n"
145 " -quality value JPEG/MIFF/PNG compression level\n"
146 " -quiet suppress all warning messages\n"
147 " -quantize colorspace reduce colors in this colorspace\n"
148 " -read-mask filename associate a read mask with the image\n"
149 " -regard-warnings pay attention to warning messages\n"
150 " -respect-parentheses settings remain in effect until parenthesis boundary\n"
151 " -sampling-factor geometry\n"
152 " horizontal and vertical sampling factor\n"
153 " -seed value seed a new sequence of pseudo-random numbers\n"
154 " -set attribute value set an image attribute\n"
155 " -quality value JPEG/MIFF/PNG compression level\n"
156 " -repage geometry size and location of an image canvas\n"
157 " -similarity-threshold value\n"
158 " minimum distortion for (sub)image match\n"
159 " -size geometry width and height of image\n"
160 " -subimage-search search for subimage\n"
161 " -synchronize synchronize image to storage device\n"
162 " -taint declare the image as modified\n"
163 " -transparent-color color\n"
164 " transparent color\n"
165 " -type type image type\n"
166 " -verbose print detailed information about the image\n"
167 " -version print version information\n"
168 " -virtual-pixel method\n"
169 " virtual pixel access method\n"
170 " -write-mask filename associate a write mask with the image",
171 stack_operators[] =
172 " -delete indexes delete the image from the image sequence";
173
174 ListMagickVersion(stdout);
175 (void) printf("Usage: %s [options ...] image reconstruct difference\n",
176 GetClientName());
177 (void) printf("\nImage Settings:\n");
178 (void) puts(settings);
179 (void) printf("\nImage Operators:\n");
180 (void) puts(operators);
181 (void) printf("\nImage Channel Operators:\n");
182 (void) puts(channel_operators);
183 (void) printf("\nImage Sequence Operators:\n");
184 (void) puts(sequence_operators);
185 (void) printf("\nImage Stack Operators:\n");
186 (void) puts(stack_operators);
187 (void) printf("\nMiscellaneous Options:\n");
188 (void) puts(miscellaneous);
189 (void) printf(
190 "\nBy default, the image format of 'file' is determined by its magic\n");
191 (void) printf(
192 "number. To specify a particular image format, precede the filename\n");
193 (void) printf(
194 "with an image format name and a colon (i.e. ps:image) or specify the\n");
195 (void) printf(
196 "image type as the filename suffix (i.e. image.ps). Specify 'file' as\n");
197 (void) printf("'-' for standard input or output.\n");
198 return(MagickTrue);
199}
200
201WandExport MagickBooleanType CompareImagesCommand(ImageInfo *image_info,
202 int argc,char **argv,char **metadata,ExceptionInfo *exception)
203{
204#define CompareEpsilon (1.0e-06)
205#define DefaultDissimilarityThreshold (1.0/MagickPI)
206#define DefaultSimilarityThreshold (-1.0)
207#define DestroyCompare() \
208{ \
209 if (similarity_image != (Image *) NULL) \
210 similarity_image=DestroyImageList(similarity_image); \
211 if (difference_image != (Image *) NULL) \
212 difference_image=DestroyImageList(difference_image); \
213 DestroyImageStack(); \
214 for (i=0; i < (ssize_t) argc; i++) \
215 argv[i]=DestroyString(argv[i]); \
216 argv=(char **) RelinquishMagickMemory(argv); \
217}
218#define ThrowCompareException(asperity,tag,option) \
219{ \
220 if (exception->severity < (asperity)) \
221 (void) ThrowMagickException(exception,GetMagickModule(),asperity,tag, \
222 "`%s'",option); \
223 DestroyCompare(); \
224 return(MagickFalse); \
225}
226#define ThrowCompareInvalidArgumentException(option,argument) \
227{ \
228 (void) ThrowMagickException(exception,GetMagickModule(),OptionError, \
229 "InvalidArgument","'%s': %s",option,argument); \
230 DestroyCompare(); \
231 return(MagickFalse); \
232}
233
234 char
235 *filename,
236 *option;
237
238 const char
239 *format;
240
241 double
242 dissimilarity_threshold,
243 distortion,
244 similarity_metric,
245 similarity_threshold;
246
247 Image
248 *difference_image,
249 *image = (Image *) NULL,
250 *reconstruct_image,
251 *similarity_image;
252
254 image_stack[MaxImageStackDepth+1];
255
256 MagickBooleanType
257 fire,
258 pend,
259 respect_parentheses,
260 similar = MagickTrue,
261 subimage_search;
262
263 MagickStatusType
264 status;
265
266 MetricType
267 metric;
268
269 RectangleInfo
270 offset;
271
272 ssize_t
273 i;
274
275 ssize_t
276 j,
277 k;
278
279 /*
280 Set defaults.
281 */
282 assert(image_info != (ImageInfo *) NULL);
283 assert(image_info->signature == MagickCoreSignature);
284 assert(exception != (ExceptionInfo *) NULL);
285 if (IsEventLogging() != MagickFalse)
286 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
287 if (argc == 2)
288 {
289 option=argv[1];
290 if ((LocaleCompare("version",option+1) == 0) ||
291 (LocaleCompare("-version",option+1) == 0))
292 {
293 ListMagickVersion(stdout);
294 return(MagickTrue);
295 }
296 }
297 if (argc < 3)
298 {
299 (void) ThrowMagickException(exception,GetMagickModule(),OptionError,
300 "MissingArgument","%s","");
301 (void) CompareUsage();
302 return(MagickFalse);
303 }
304 difference_image=NewImageList();
305 similarity_image=NewImageList();
306 dissimilarity_threshold=DefaultDissimilarityThreshold;
307 similarity_threshold=DefaultSimilarityThreshold;
308 distortion=0.0;
309 similarity_metric=0.0;
310 format=(char *) NULL;
311 j=1;
312 k=0;
313 metric=UndefinedErrorMetric;
314 NewImageStack();
315 option=(char *) NULL;
316 pend=MagickFalse;
317 reconstruct_image=NewImageList();
318 respect_parentheses=MagickFalse;
319 status=MagickTrue;
320 subimage_search=MagickFalse;
321 /*
322 Compare an image.
323 */
324 ReadCommandlLine(argc,&argv);
325 status=ExpandFilenames(&argc,&argv);
326 if (status == MagickFalse)
327 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
328 GetExceptionMessage(errno));
329 for (i=1; i < (ssize_t) (argc-1); i++)
330 {
331 option=argv[i];
332 if (LocaleCompare(option,"(") == 0)
333 {
334 FireImageStack(MagickTrue,MagickTrue,pend);
335 if (k == MaxImageStackDepth)
336 ThrowCompareException(OptionError,"ParenthesisNestedTooDeeply",
337 option);
338 PushImageStack();
339 continue;
340 }
341 if (LocaleCompare(option,")") == 0)
342 {
343 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
344 if (k == 0)
345 ThrowCompareException(OptionError,"UnableToParseExpression",option);
346 PopImageStack();
347 continue;
348 }
349 if (IsCommandOption(option) == MagickFalse)
350 {
351 Image
352 *images;
353
354 /*
355 Read input image.
356 */
357 FireImageStack(MagickFalse,MagickFalse,pend);
358 filename=argv[i];
359 if ((LocaleCompare(filename,"--") == 0) && (i < (ssize_t) (argc-1)))
360 filename=argv[++i];
361 images=ReadImages(image_info,filename,exception);
362 status&=(MagickStatusType) (images != (Image *) NULL) &&
363 (exception->severity < ErrorException);
364 if (images == (Image *) NULL)
365 continue;
366 AppendImageStack(images);
367 continue;
368 }
369 pend=image != (Image *) NULL ? MagickTrue : MagickFalse;
370 switch (*(option+1))
371 {
372 case 'a':
373 {
374 if (LocaleCompare("adjoin",option+1) == 0)
375 break;
376 if (LocaleCompare("alpha",option+1) == 0)
377 {
378 ssize_t
379 type;
380
381 if (*option == '+')
382 break;
383 i++;
384 if (i == (ssize_t) argc)
385 ThrowCompareException(OptionError,"MissingArgument",option);
386 type=ParseCommandOption(MagickAlphaChannelOptions,MagickFalse,
387 argv[i]);
388 if (type < 0)
389 ThrowCompareException(OptionError,
390 "UnrecognizedAlphaChannelOption",argv[i]);
391 break;
392 }
393 if (LocaleCompare("auto-orient",option+1) == 0)
394 break;
395 if (LocaleCompare("authenticate",option+1) == 0)
396 {
397 if (*option == '+')
398 break;
399 i++;
400 if (i == (ssize_t) argc)
401 ThrowCompareException(OptionError,"MissingArgument",option);
402 break;
403 }
404 ThrowCompareException(OptionError,"UnrecognizedOption",option);
405 }
406 case 'b':
407 {
408 if (LocaleCompare("background",option+1) == 0)
409 {
410 if (*option == '+')
411 break;
412 i++;
413 if (i == (ssize_t) argc)
414 ThrowCompareException(OptionError,"MissingArgument",option);
415 break;
416 }
417 if (LocaleCompare("brightness-contrast",option+1) == 0)
418 {
419 i++;
420 if (i == (ssize_t) argc)
421 ThrowCompareException(OptionError,"MissingArgument",option);
422 if (IsGeometry(argv[i]) == MagickFalse)
423 ThrowCompareInvalidArgumentException(option,argv[i]);
424 break;
425 }
426 ThrowCompareException(OptionError,"UnrecognizedOption",option);
427 }
428 case 'c':
429 {
430 if (LocaleCompare("cache",option+1) == 0)
431 {
432 if (*option == '+')
433 break;
434 i++;
435 if (i == (ssize_t) argc)
436 ThrowCompareException(OptionError,"MissingArgument",option);
437 if (IsGeometry(argv[i]) == MagickFalse)
438 ThrowCompareInvalidArgumentException(option,argv[i]);
439 break;
440 }
441 if (LocaleCompare("channel",option+1) == 0)
442 {
443 ssize_t
444 channel;
445
446 if (*option == '+')
447 break;
448 i++;
449 if (i == (ssize_t) argc)
450 ThrowCompareException(OptionError,"MissingArgument",option);
451 channel=ParseChannelOption(argv[i]);
452 if (channel < 0)
453 ThrowCompareException(OptionError,"UnrecognizedChannelType",
454 argv[i]);
455 break;
456 }
457 if (LocaleCompare("colorspace",option+1) == 0)
458 {
459 ssize_t
460 colorspace;
461
462 if (*option == '+')
463 break;
464 i++;
465 if (i == (ssize_t) argc)
466 ThrowCompareException(OptionError,"MissingArgument",option);
467 colorspace=ParseCommandOption(MagickColorspaceOptions,MagickFalse,
468 argv[i]);
469 if (colorspace < 0)
470 ThrowCompareException(OptionError,"UnrecognizedColorspace",
471 argv[i]);
472 break;
473 }
474 if (LocaleCompare("compose",option+1) == 0)
475 {
476 ssize_t
477 compose;
478
479 if (*option == '+')
480 break;
481 i++;
482 if (i == (ssize_t) argc)
483 ThrowCompareException(OptionError,"MissingArgument",option);
484 compose=ParseCommandOption(MagickComposeOptions,MagickFalse,
485 argv[i]);
486 if (compose < 0)
487 ThrowCompareException(OptionError,"UnrecognizedComposeOperator",
488 argv[i]);
489 break;
490 }
491 if (LocaleCompare("compress",option+1) == 0)
492 {
493 ssize_t
494 compress;
495
496 if (*option == '+')
497 break;
498 i++;
499 if (i == (ssize_t) argc)
500 ThrowCompareException(OptionError,"MissingArgument",option);
501 compress=ParseCommandOption(MagickCompressOptions,MagickFalse,
502 argv[i]);
503 if (compress < 0)
504 ThrowCompareException(OptionError,"UnrecognizedImageCompression",
505 argv[i]);
506 break;
507 }
508 if (LocaleCompare("concurrent",option+1) == 0)
509 break;
510 if (LocaleCompare("crop",option+1) == 0)
511 {
512 if (*option == '+')
513 break;
514 i++;
515 if (i == (ssize_t) argc)
516 ThrowCompareException(OptionError,"MissingArgument",option);
517 if (IsGeometry(argv[i]) == MagickFalse)
518 ThrowCompareInvalidArgumentException(option,argv[i]);
519 break;
520 }
521 ThrowCompareException(OptionError,"UnrecognizedOption",option)
522 }
523 case 'd':
524 {
525 if (LocaleCompare("debug",option+1) == 0)
526 {
527 LogEventType
528 event_mask;
529
530 if (*option == '+')
531 break;
532 i++;
533 if (i == (ssize_t) argc)
534 ThrowCompareException(OptionError,"MissingArgument",option);
535 event_mask=SetLogEventMask(argv[i]);
536 if (event_mask == UndefinedEvents)
537 ThrowCompareException(OptionError,"UnrecognizedEventType",
538 argv[i]);
539 break;
540 }
541 if (LocaleCompare("decipher",option+1) == 0)
542 {
543 if (*option == '+')
544 break;
545 i++;
546 if (i == (ssize_t) argc)
547 ThrowCompareException(OptionError,"MissingArgument",option);
548 break;
549 }
550 if (LocaleCompare("define",option+1) == 0)
551 {
552 i++;
553 if (i == (ssize_t) argc)
554 ThrowCompareException(OptionError,"MissingArgument",option);
555 if (*option == '+')
556 {
557 const char
558 *define;
559
560 define=GetImageOption(image_info,argv[i]);
561 if (define == (const char *) NULL)
562 ThrowCompareException(OptionError,"NoSuchOption",argv[i]);
563 break;
564 }
565 break;
566 }
567 if (LocaleCompare("delete",option+1) == 0)
568 {
569 if (*option == '+')
570 break;
571 i++;
572 if (i == (ssize_t) argc)
573 ThrowCompareException(OptionError,"MissingArgument",option);
574 if (IsSceneGeometry(argv[i],MagickFalse) == MagickFalse)
575 ThrowCompareInvalidArgumentException(option,argv[i]);
576 break;
577 }
578 if (LocaleCompare("density",option+1) == 0)
579 {
580 if (*option == '+')
581 break;
582 i++;
583 if (i == (ssize_t) argc)
584 ThrowCompareException(OptionError,"MissingArgument",option);
585 if (IsGeometry(argv[i]) == MagickFalse)
586 ThrowCompareInvalidArgumentException(option,argv[i]);
587 break;
588 }
589 if (LocaleCompare("depth",option+1) == 0)
590 {
591 if (*option == '+')
592 break;
593 i++;
594 if (i == (ssize_t) argc)
595 ThrowCompareException(OptionError,"MissingArgument",option);
596 if (IsGeometry(argv[i]) == MagickFalse)
597 ThrowCompareInvalidArgumentException(option,argv[i]);
598 break;
599 }
600 if (LocaleCompare("dissimilarity-threshold",option+1) == 0)
601 {
602 if (*option == '+')
603 break;
604 i++;
605 if (i == (ssize_t) argc)
606 ThrowCompareException(OptionError,"MissingArgument",option);
607 if (IsGeometry(argv[i]) == MagickFalse)
608 ThrowCompareInvalidArgumentException(option,argv[i]);
609 if (*option == '+')
610 dissimilarity_threshold=DefaultDissimilarityThreshold;
611 else
612 dissimilarity_threshold=StringToDouble(argv[i],(char **) NULL);
613 break;
614 }
615 if (LocaleCompare("distort",option+1) == 0)
616 {
617 ssize_t
618 op;
619
620 i++;
621 if (i == (ssize_t) argc)
622 ThrowCompareException(OptionError,"MissingArgument",option);
623 op=ParseCommandOption(MagickDistortOptions,MagickFalse,argv[i]);
624 if (op < 0)
625 ThrowCompareException(OptionError,"UnrecognizedDistortMethod",
626 argv[i]);
627 i++;
628 if (i == (ssize_t) argc)
629 ThrowCompareException(OptionError,"MissingArgument",option);
630 break;
631 }
632 if (LocaleCompare("duration",option+1) == 0)
633 {
634 if (*option == '+')
635 break;
636 i++;
637 if (i == (ssize_t) argc)
638 ThrowCompareException(OptionError,"MissingArgument",option);
639 if (IsGeometry(argv[i]) == MagickFalse)
640 ThrowCompareInvalidArgumentException(option,argv[i]);
641 break;
642 }
643 ThrowCompareException(OptionError,"UnrecognizedOption",option)
644 }
645 case 'e':
646 {
647 if (LocaleCompare("encipher",option+1) == 0)
648 {
649 if (*option == '+')
650 break;
651 i++;
652 if (i == (ssize_t) argc)
653 ThrowCompareException(OptionError,"MissingArgument",option);
654 break;
655 }
656 if (LocaleCompare("extract",option+1) == 0)
657 {
658 if (*option == '+')
659 break;
660 i++;
661 if (i == (ssize_t) argc)
662 ThrowCompareException(OptionError,"MissingArgument",option);
663 if (IsGeometry(argv[i]) == MagickFalse)
664 ThrowCompareInvalidArgumentException(option,argv[i]);
665 break;
666 }
667 ThrowCompareException(OptionError,"UnrecognizedOption",option)
668 }
669 case 'f':
670 {
671 if (LocaleCompare("format",option+1) == 0)
672 {
673 if (*option == '+')
674 break;
675 i++;
676 if (i == (ssize_t) argc)
677 ThrowCompareException(OptionError,"MissingArgument",option);
678 format=argv[i];
679 break;
680 }
681 if (LocaleCompare("fuzz",option+1) == 0)
682 {
683 if (*option == '+')
684 break;
685 i++;
686 if (i == (ssize_t) argc)
687 ThrowCompareException(OptionError,"MissingArgument",option);
688 if (IsGeometry(argv[i]) == MagickFalse)
689 ThrowCompareInvalidArgumentException(option,argv[i]);
690 break;
691 }
692 ThrowCompareException(OptionError,"UnrecognizedOption",option)
693 }
694 case 'g':
695 {
696 if (LocaleCompare("gravity",option+1) == 0)
697 {
698 ssize_t
699 gravity;
700
701 if (*option == '+')
702 break;
703 i++;
704 if (i == (ssize_t) argc)
705 ThrowCompareException(OptionError,"MissingArgument",option);
706 gravity=ParseCommandOption(MagickGravityOptions,MagickFalse,
707 argv[i]);
708 if (gravity < 0)
709 ThrowCompareException(OptionError,"UnrecognizedGravityType",
710 argv[i]);
711 break;
712 }
713 ThrowCompareException(OptionError,"UnrecognizedOption",option)
714 }
715 case 'h':
716 {
717 if ((LocaleCompare("help",option+1) == 0) ||
718 (LocaleCompare("-help",option+1) == 0))
719 {
720 DestroyCompare();
721 return(CompareUsage());
722 }
723 if (LocaleCompare("highlight-color",option+1) == 0)
724 {
725 if (*option == '+')
726 break;
727 i++;
728 if (i == (ssize_t) argc)
729 ThrowCompareException(OptionError,"MissingArgument",option);
730 break;
731 }
732 ThrowCompareException(OptionError,"UnrecognizedOption",option)
733 }
734 case 'i':
735 {
736 if (LocaleCompare("identify",option+1) == 0)
737 break;
738 if (LocaleCompare("interlace",option+1) == 0)
739 {
740 ssize_t
741 interlace;
742
743 if (*option == '+')
744 break;
745 i++;
746 if (i == (ssize_t) argc)
747 ThrowCompareException(OptionError,"MissingArgument",option);
748 interlace=ParseCommandOption(MagickInterlaceOptions,MagickFalse,
749 argv[i]);
750 if (interlace < 0)
751 ThrowCompareException(OptionError,"UnrecognizedInterlaceType",
752 argv[i]);
753 break;
754 }
755 ThrowCompareException(OptionError,"UnrecognizedOption",option)
756 }
757 case 'l':
758 {
759 if (LocaleCompare("level",option+1) == 0)
760 {
761 i++;
762 if (i == (ssize_t) argc)
763 ThrowCompareException(OptionError,"MissingArgument",option);
764 if (IsGeometry(argv[i]) == MagickFalse)
765 ThrowCompareInvalidArgumentException(option,argv[i]);
766 break;
767 }
768 if (LocaleCompare("limit",option+1) == 0)
769 {
770 char
771 *p;
772
773 double
774 value;
775
776 ssize_t
777 resource;
778
779 if (*option == '+')
780 break;
781 i++;
782 if (i == (ssize_t) argc)
783 ThrowCompareException(OptionError,"MissingArgument",option);
784 resource=ParseCommandOption(MagickResourceOptions,MagickFalse,
785 argv[i]);
786 if (resource < 0)
787 ThrowCompareException(OptionError,"UnrecognizedResourceType",
788 argv[i]);
789 i++;
790 if (i == (ssize_t) argc)
791 ThrowCompareException(OptionError,"MissingArgument",option);
792 value=StringToDouble(argv[i],&p);
793 (void) value;
794 if ((p == argv[i]) && (LocaleCompare("unlimited",argv[i]) != 0))
795 ThrowCompareInvalidArgumentException(option,argv[i]);
796 break;
797 }
798 if (LocaleCompare("list",option+1) == 0)
799 {
800 ssize_t
801 list;
802
803 if (*option == '+')
804 break;
805 i++;
806 if (i == (ssize_t) argc)
807 ThrowCompareException(OptionError,"MissingArgument",option);
808 list=ParseCommandOption(MagickListOptions,MagickFalse,argv[i]);
809 if (list < 0)
810 ThrowCompareException(OptionError,"UnrecognizedListType",argv[i]);
811 status=MogrifyImageInfo(image_info,(int) (i-j+1),(const char **)
812 argv+j,exception);
813 DestroyCompare();
814 return(status == 0 ? MagickFalse : MagickTrue);
815 }
816 if (LocaleCompare("log",option+1) == 0)
817 {
818 if (*option == '+')
819 break;
820 i++;
821 if ((i == (ssize_t) argc) || (strchr(argv[i],'%') == (char *) NULL))
822 ThrowCompareException(OptionError,"MissingArgument",option);
823 break;
824 }
825 if (LocaleCompare("lowlight-color",option+1) == 0)
826 {
827 if (*option == '+')
828 break;
829 i++;
830 if (i == (ssize_t) argc)
831 ThrowCompareException(OptionError,"MissingArgument",option);
832 break;
833 }
834 ThrowCompareException(OptionError,"UnrecognizedOption",option)
835 }
836 case 'm':
837 {
838 if (LocaleCompare("matte",option+1) == 0)
839 break;
840 if (LocaleCompare("metric",option+1) == 0)
841 {
842 ssize_t
843 type;
844
845 if (*option == '+')
846 break;
847 i++;
848 if (i == (ssize_t) argc)
849 ThrowCompareException(OptionError,"MissingArgument",option);
850 type=ParseCommandOption(MagickMetricOptions,MagickTrue,argv[i]);
851 if (type < 0)
852 ThrowCompareException(OptionError,"UnrecognizedMetricType",
853 argv[i]);
854 metric=(MetricType) type;
855 break;
856 }
857 if (LocaleCompare("monitor",option+1) == 0)
858 break;
859 ThrowCompareException(OptionError,"UnrecognizedOption",option)
860 }
861 case 'n':
862 {
863 if (LocaleCompare("negate",option+1) == 0)
864 break;
865 ThrowCompareException(OptionError,"UnrecognizedOption",option)
866 }
867 case 'p':
868 {
869 if (LocaleCompare("passphrase",option+1) == 0)
870 {
871 if (*option == '+')
872 break;
873 i++;
874 if (i == (ssize_t) argc)
875 ThrowCompareException(OptionError,"MissingArgument",option);
876 break;
877 }
878 if (LocaleCompare("precision",option+1) == 0)
879 {
880 if (*option == '+')
881 break;
882 i++;
883 if (i == (ssize_t) argc)
884 ThrowCompareException(OptionError,"MissingArgument",option);
885 if (IsGeometry(argv[i]) == MagickFalse)
886 ThrowCompareInvalidArgumentException(option,argv[i]);
887 break;
888 }
889 if (LocaleCompare("profile",option+1) == 0)
890 {
891 i++;
892 if (i == (ssize_t) argc)
893 ThrowCompareException(OptionError,"MissingArgument",option);
894 break;
895 }
896 ThrowCompareException(OptionError,"UnrecognizedOption",option)
897 }
898 case 'q':
899 {
900 if (LocaleCompare("quality",option+1) == 0)
901 {
902 if (*option == '+')
903 break;
904 i++;
905 if (i == (ssize_t) argc)
906 ThrowCompareException(OptionError,"MissingArgument",option);
907 if (IsGeometry(argv[i]) == MagickFalse)
908 ThrowCompareInvalidArgumentException(option,argv[i]);
909 break;
910 }
911 if (LocaleCompare("quantize",option+1) == 0)
912 {
913 ssize_t
914 colorspace;
915
916 if (*option == '+')
917 break;
918 i++;
919 if (i == (ssize_t) argc)
920 ThrowCompareException(OptionError,"MissingArgument",option);
921 colorspace=ParseCommandOption(MagickColorspaceOptions,
922 MagickFalse,argv[i]);
923 if (colorspace < 0)
924 ThrowCompareException(OptionError,"UnrecognizedColorspace",
925 argv[i]);
926 break;
927 }
928 if (LocaleCompare("quiet",option+1) == 0)
929 break;
930 ThrowCompareException(OptionError,"UnrecognizedOption",option)
931 }
932 case 'r':
933 {
934 if (LocaleCompare("read-mask",option+1) == 0)
935 {
936 if (*option == '+')
937 break;
938 i++;
939 if (i == (ssize_t) argc)
940 ThrowCompareException(OptionError,"MissingArgument",option);
941 break;
942 }
943 if (LocaleCompare("regard-warnings",option+1) == 0)
944 break;
945 if (LocaleCompare("repage",option+1) == 0)
946 {
947 if (*option == '+')
948 break;
949 i++;
950 if (i == (ssize_t) argc)
951 ThrowCompareException(OptionError,"MissingArgument",option);
952 if (IsGeometry(argv[i]) == MagickFalse)
953 ThrowCompareInvalidArgumentException(option,argv[i]);
954 break;
955 }
956 if (LocaleCompare("resize",option+1) == 0)
957 {
958 if (*option == '+')
959 break;
960 i++;
961 if (i == (ssize_t) argc)
962 ThrowCompareException(OptionError,"MissingArgument",option);
963 if (IsGeometry(argv[i]) == MagickFalse)
964 ThrowCompareInvalidArgumentException(option,argv[i]);
965 break;
966 }
967 if (LocaleNCompare("respect-parentheses",option+1,17) == 0)
968 {
969 respect_parentheses=(*option == '-') ? MagickTrue : MagickFalse;
970 break;
971 }
972 if (LocaleCompare("rotate",option+1) == 0)
973 {
974 i++;
975 if (i == (ssize_t) argc)
976 ThrowCompareException(OptionError,"MissingArgument",option);
977 if (IsGeometry(argv[i]) == MagickFalse)
978 ThrowCompareInvalidArgumentException(option,argv[i]);
979 break;
980 }
981 ThrowCompareException(OptionError,"UnrecognizedOption",option)
982 }
983 case 's':
984 {
985 if (LocaleCompare("sampling-factor",option+1) == 0)
986 {
987 if (*option == '+')
988 break;
989 i++;
990 if (i == (ssize_t) argc)
991 ThrowCompareException(OptionError,"MissingArgument",option);
992 if (IsGeometry(argv[i]) == MagickFalse)
993 ThrowCompareInvalidArgumentException(option,argv[i]);
994 break;
995 }
996 if (LocaleCompare("seed",option+1) == 0)
997 {
998 if (*option == '+')
999 break;
1000 i++;
1001 if (i == (ssize_t) argc)
1002 ThrowCompareException(OptionError,"MissingArgument",option);
1003 if (IsGeometry(argv[i]) == MagickFalse)
1004 ThrowCompareInvalidArgumentException(option,argv[i]);
1005 break;
1006 }
1007 if (LocaleCompare("separate",option+1) == 0)
1008 break;
1009 if (LocaleCompare("set",option+1) == 0)
1010 {
1011 i++;
1012 if (i == (ssize_t) argc)
1013 ThrowCompareException(OptionError,"MissingArgument",option);
1014 if (*option == '+')
1015 break;
1016 i++;
1017 if (i == (ssize_t) argc)
1018 ThrowCompareException(OptionError,"MissingArgument",option);
1019 break;
1020 }
1021 if (LocaleCompare("sigmoidal-contrast",option+1) == 0)
1022 {
1023 i++;
1024 if (i == (ssize_t) argc)
1025 ThrowCompareException(OptionError,"MissingArgument",option);
1026 if (IsGeometry(argv[i]) == MagickFalse)
1027 ThrowCompareInvalidArgumentException(option,argv[i]);
1028 break;
1029 }
1030 if (LocaleCompare("similarity-threshold",option+1) == 0)
1031 {
1032 if (*option == '+')
1033 break;
1034 i++;
1035 if (i == (ssize_t) argc)
1036 ThrowCompareException(OptionError,"MissingArgument",option);
1037 if (IsGeometry(argv[i]) == MagickFalse)
1038 ThrowCompareInvalidArgumentException(option,argv[i]);
1039 if (*option == '+')
1040 similarity_threshold=DefaultSimilarityThreshold;
1041 else
1042 similarity_threshold=StringToDouble(argv[i],(char **) NULL);
1043 break;
1044 }
1045 if (LocaleCompare("size",option+1) == 0)
1046 {
1047 if (*option == '+')
1048 break;
1049 i++;
1050 if (i == (ssize_t) argc)
1051 ThrowCompareException(OptionError,"MissingArgument",option);
1052 if (IsGeometry(argv[i]) == MagickFalse)
1053 ThrowCompareInvalidArgumentException(option,argv[i]);
1054 break;
1055 }
1056 if (LocaleCompare("subimage-search",option+1) == 0)
1057 {
1058 if (*option == '+')
1059 {
1060 subimage_search=MagickFalse;
1061 break;
1062 }
1063 subimage_search=MagickTrue;
1064 break;
1065 }
1066 if (LocaleCompare("synchronize",option+1) == 0)
1067 break;
1068 ThrowCompareException(OptionError,"UnrecognizedOption",option)
1069 }
1070 case 't':
1071 {
1072 if (LocaleCompare("taint",option+1) == 0)
1073 break;
1074 if (LocaleCompare("transparent-color",option+1) == 0)
1075 {
1076 if (*option == '+')
1077 break;
1078 i++;
1079 if (i == (ssize_t) argc)
1080 ThrowCompareException(OptionError,"MissingArgument",option);
1081 break;
1082 }
1083 if (LocaleCompare("trim",option+1) == 0)
1084 break;
1085 if (LocaleCompare("type",option+1) == 0)
1086 {
1087 ssize_t
1088 type;
1089
1090 if (*option == '+')
1091 break;
1092 i++;
1093 if (i == (ssize_t) argc)
1094 ThrowCompareException(OptionError,"MissingArgument",option);
1095 type=ParseCommandOption(MagickTypeOptions,MagickFalse,argv[i]);
1096 if (type < 0)
1097 ThrowCompareException(OptionError,"UnrecognizedImageType",
1098 argv[i]);
1099 break;
1100 }
1101 ThrowCompareException(OptionError,"UnrecognizedOption",option)
1102 }
1103 case 'v':
1104 {
1105 if (LocaleCompare("verbose",option+1) == 0)
1106 break;
1107 if ((LocaleCompare("version",option+1) == 0) ||
1108 (LocaleCompare("-version",option+1) == 0))
1109 {
1110 ListMagickVersion(stdout);
1111 break;
1112 }
1113 if (LocaleCompare("virtual-pixel",option+1) == 0)
1114 {
1115 ssize_t
1116 method;
1117
1118 if (*option == '+')
1119 break;
1120 i++;
1121 if (i == (ssize_t) argc)
1122 ThrowCompareException(OptionError,"MissingArgument",option);
1123 method=ParseCommandOption(MagickVirtualPixelOptions,MagickFalse,
1124 argv[i]);
1125 if (method < 0)
1126 ThrowCompareException(OptionError,
1127 "UnrecognizedVirtualPixelMethod",argv[i]);
1128 break;
1129 }
1130 ThrowCompareException(OptionError,"UnrecognizedOption",option)
1131 }
1132 case 'w':
1133 {
1134 if (LocaleCompare("write",option+1) == 0)
1135 {
1136 i++;
1137 if (i == (ssize_t) argc)
1138 ThrowCompareException(OptionError,"MissingArgument",option);
1139 break;
1140 }
1141 if (LocaleCompare("write-mask",option+1) == 0)
1142 {
1143 if (*option == '+')
1144 break;
1145 i++;
1146 if (i == (ssize_t) argc)
1147 ThrowCompareException(OptionError,"MissingArgument",option);
1148 break;
1149 }
1150 ThrowCompareException(OptionError,"UnrecognizedOption",option)
1151 }
1152 case '?':
1153 break;
1154 default:
1155 ThrowCompareException(OptionError,"UnrecognizedOption",option)
1156 }
1157 fire=(GetCommandOptionFlags(MagickCommandOptions,MagickFalse,option) &
1158 FireOptionFlag) == 0 ? MagickFalse : MagickTrue;
1159 if (fire != MagickFalse)
1160 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
1161 }
1162 if (k != 0)
1163 ThrowCompareException(OptionError,"UnbalancedParenthesis",argv[i]);
1164 if (i-- != (ssize_t) (argc-1))
1165 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
1166 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
1167 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
1168 FinalizeImageSettings(image_info,image,MagickTrue);
1169 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
1170 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
1171 image=GetImageFromList(image,0);
1172 reconstruct_image=GetImageFromList(image,1);
1173 offset.x=0;
1174 offset.y=0;
1175 if (subimage_search != MagickFalse)
1176 {
1177 similarity_image=SimilarityImage(image,reconstruct_image,metric,
1178 similarity_threshold,&offset,&similarity_metric,exception);
1179 if (dissimilarity_threshold == DefaultDissimilarityThreshold)
1180 switch (metric)
1181 {
1182 case PhaseCorrelationErrorMetric:
1183 case PeakSignalToNoiseRatioErrorMetric:
1184 {
1185 dissimilarity_threshold=1.0;
1186 break;
1187 }
1188 default:
1189 break;
1190 }
1191 if (similarity_metric > (dissimilarity_threshold+MagickEpsilon))
1192 (void) ThrowMagickException(exception,GetMagickModule(),ImageWarning,
1193 "ImagesTooDissimilar","`%s'",image->filename);
1194 }
1195 if (similarity_image == (Image *) NULL)
1196 difference_image=CompareImages(image,reconstruct_image,metric,&distortion,
1197 exception);
1198 else
1199 {
1200 Image
1201 *composite_image;
1202
1203 /*
1204 Determine if reconstructed image is a subimage of the image.
1205 */
1206 composite_image=CloneImage(image,0,0,MagickTrue,exception);
1207 if (composite_image == (Image *) NULL)
1208 difference_image=CompareImages(image,reconstruct_image,metric,
1209 &distortion,exception);
1210 else
1211 {
1212 Image
1213 *distort_image;
1214
1215 RectangleInfo
1216 page;
1217
1218 (void) CompositeImage(composite_image,reconstruct_image,
1219 CopyCompositeOp,MagickTrue,offset.x,offset.y,exception);
1220 difference_image=CompareImages(image,composite_image,metric,
1221 &distortion,exception);
1222 if (difference_image != (Image *) NULL)
1223 {
1224 difference_image->page.x=offset.x;
1225 difference_image->page.y=offset.y;
1226 }
1227 composite_image=DestroyImage(composite_image);
1228 page.width=reconstruct_image->columns;
1229 page.height=reconstruct_image->rows;
1230 page.x=offset.x;
1231 page.y=offset.y;
1232 distort_image=CropImage(image,&page,exception);
1233 if (distort_image != (Image *) NULL)
1234 {
1235 Image
1236 *sans_image;
1237
1238 sans_image=CompareImages(distort_image,reconstruct_image,metric,
1239 &distortion,exception);
1240 distort_image=DestroyImage(distort_image);
1241 if (sans_image != (Image *) NULL)
1242 sans_image=DestroyImage(sans_image);
1243 }
1244 }
1245 if (difference_image != (Image *) NULL)
1246 {
1247 AppendImageToList(&difference_image,similarity_image);
1248 similarity_image=(Image *) NULL;
1249 }
1250 }
1251 switch (metric)
1252 {
1253 case DotProductCorrelationErrorMetric:
1254 case NormalizedCrossCorrelationErrorMetric:
1255 case StructuralSimilarityErrorMetric:
1256 case StructuralDissimilarityErrorMetric:
1257 case UndefinedErrorMetric:
1258 {
1259 distortion=1.0-distortion;
1260 similarity_metric=1.0-similarity_metric;
1261 break;
1262 }
1263 case PeakSignalToNoiseRatioErrorMetric:
1264 {
1265 distortion=fabs(distortion);
1266 similarity_metric=fabs(similarity_metric);
1267 break;
1268 }
1269 case PhaseCorrelationErrorMetric:
1270 {
1271 distortion=1.0-distortion;
1272 break;
1273 }
1274 default: break;
1275 }
1276 if (fabs(distortion) > CompareEpsilon)
1277 similar=MagickFalse;
1278 if (difference_image == (Image *) NULL)
1279 status=0;
1280 else
1281 {
1282 if (image_info->verbose != MagickFalse)
1283 (void) SetImageColorMetric(image,reconstruct_image,exception);
1284 if (*difference_image->magick == '\0')
1285 (void) CopyMagickString(difference_image->magick,image->magick,
1286 MagickPathExtent);
1287 if (image_info->verbose == MagickFalse)
1288 {
1289 switch (metric)
1290 {
1291 case AbsoluteErrorMetric:
1292 case DotProductCorrelationErrorMetric:
1293 case FuzzErrorMetric:
1294 case NormalizedCrossCorrelationErrorMetric:
1295 case PerceptualHashErrorMetric:
1296 case PhaseCorrelationErrorMetric:
1297 case StructuralSimilarityErrorMetric:
1298 case StructuralDissimilarityErrorMetric:
1299 case MeanAbsoluteErrorMetric:
1300 case MeanSquaredErrorMetric:
1301 case PeakAbsoluteErrorMetric:
1302 case RootMeanSquaredErrorMetric:
1303 {
1304 (void) FormatLocaleFile(stderr,"%.*g (%.*g)",GetMagickPrecision(),
1305 (double) QuantumRange*distortion,GetMagickPrecision(),
1306 distortion);
1307 break;
1308 }
1309 case PeakSignalToNoiseRatioErrorMetric:
1310 {
1311 (void) FormatLocaleFile(stderr,"%.*g (%.*g)",GetMagickPrecision(),
1312 48.1647*distortion,GetMagickPrecision(),distortion);
1313 break;
1314 }
1315 case MeanErrorPerPixelErrorMetric:
1316 {
1317 (void) FormatLocaleFile(stderr,"%.*g (%.*g, %.*g)",
1318 GetMagickPrecision(),distortion,
1319 GetMagickPrecision(),image->error.normalized_mean_error,
1320 GetMagickPrecision(),image->error.normalized_maximum_error);
1321 break;
1322 }
1323 case UndefinedErrorMetric:
1324 break;
1325 }
1326 if (subimage_search != MagickFalse)
1327 (void) FormatLocaleFile(stderr," @ %.20g,%.20g [%.*g]",
1328 (double) offset.x,(double) offset.y,GetMagickPrecision(),
1329 similarity_metric);
1330 }
1331 else
1332 {
1333 double
1334 *channel_distortion;
1335
1336 channel_distortion=GetImageDistortions(image,reconstruct_image,
1337 metric,exception);
1338 (void) FormatLocaleFile(stderr,"Image: %s\n",image->filename);
1339 if ((reconstruct_image->columns != image->columns) ||
1340 (reconstruct_image->rows != image->rows))
1341 (void) FormatLocaleFile(stderr,"Offset: %.20g,%.20g\n",(double)
1342 difference_image->page.x,(double) difference_image->page.y);
1343 (void) FormatLocaleFile(stderr," Channel distortion: %s\n",
1344 CommandOptionToMnemonic(MagickMetricOptions,(ssize_t) metric));
1345 switch (metric)
1346 {
1347 case FuzzErrorMetric:
1348 case MeanAbsoluteErrorMetric:
1349 case MeanSquaredErrorMetric:
1350 case PeakAbsoluteErrorMetric:
1351 case RootMeanSquaredErrorMetric:
1352 {
1353 switch (image->colorspace)
1354 {
1355 case RGBColorspace:
1356 default:
1357 {
1358 (void) FormatLocaleFile(stderr," red: %.*g (%.*g)\n",
1359 GetMagickPrecision(),(double) QuantumRange*
1360 channel_distortion[RedPixelChannel],GetMagickPrecision(),
1361 channel_distortion[RedPixelChannel]);
1362 (void) FormatLocaleFile(stderr," green: %.*g (%.*g)\n",
1363 GetMagickPrecision(),(double) QuantumRange*
1364 channel_distortion[GreenPixelChannel],GetMagickPrecision(),
1365 channel_distortion[GreenPixelChannel]);
1366 (void) FormatLocaleFile(stderr," blue: %.*g (%.*g)\n",
1367 GetMagickPrecision(),(double) QuantumRange*
1368 channel_distortion[BluePixelChannel],GetMagickPrecision(),
1369 channel_distortion[BluePixelChannel]);
1370 if (image->alpha_trait != UndefinedPixelTrait)
1371 (void) FormatLocaleFile(stderr," alpha: %.*g (%.*g)\n",
1372 GetMagickPrecision(),(double) QuantumRange*
1373 channel_distortion[AlphaPixelChannel],
1374 GetMagickPrecision(),
1375 channel_distortion[AlphaPixelChannel]);
1376 break;
1377 }
1378 case CMYKColorspace:
1379 {
1380 (void) FormatLocaleFile(stderr," cyan: %.*g (%.*g)\n",
1381 GetMagickPrecision(),(double) QuantumRange*
1382 channel_distortion[CyanPixelChannel],GetMagickPrecision(),
1383 channel_distortion[CyanPixelChannel]);
1384 (void) FormatLocaleFile(stderr," magenta: %.*g (%.*g)\n",
1385 GetMagickPrecision(),(double) QuantumRange*
1386 channel_distortion[MagentaPixelChannel],
1387 GetMagickPrecision(),
1388 channel_distortion[MagentaPixelChannel]);
1389 (void) FormatLocaleFile(stderr," yellow: %.*g (%.*g)\n",
1390 GetMagickPrecision(),(double) QuantumRange*
1391 channel_distortion[YellowPixelChannel],GetMagickPrecision(),
1392 channel_distortion[YellowPixelChannel]);
1393 (void) FormatLocaleFile(stderr," black: %.*g (%.*g)\n",
1394 GetMagickPrecision(),(double) QuantumRange*
1395 channel_distortion[BlackPixelChannel],GetMagickPrecision(),
1396 channel_distortion[BlackPixelChannel]);
1397 if (image->alpha_trait != UndefinedPixelTrait)
1398 (void) FormatLocaleFile(stderr," alpha: %.*g (%.*g)\n",
1399 GetMagickPrecision(),(double) QuantumRange*
1400 channel_distortion[AlphaPixelChannel],
1401 GetMagickPrecision(),
1402 channel_distortion[AlphaPixelChannel]);
1403 break;
1404 }
1405 case LinearGRAYColorspace:
1406 case GRAYColorspace:
1407 {
1408 (void) FormatLocaleFile(stderr," gray: %.*g (%.*g)\n",
1409 GetMagickPrecision(),(double) QuantumRange*
1410 channel_distortion[GrayPixelChannel],GetMagickPrecision(),
1411 channel_distortion[GrayPixelChannel]);
1412 if (image->alpha_trait != UndefinedPixelTrait)
1413 (void) FormatLocaleFile(stderr," alpha: %.*g (%.*g)\n",
1414 GetMagickPrecision(),(double) QuantumRange*
1415 channel_distortion[AlphaPixelChannel],
1416 GetMagickPrecision(),
1417 channel_distortion[AlphaPixelChannel]);
1418 break;
1419 }
1420 }
1421 (void) FormatLocaleFile(stderr," all: %.*g (%.*g)\n",
1422 GetMagickPrecision(),(double) QuantumRange*
1423 channel_distortion[MaxPixelChannels],GetMagickPrecision(),
1424 channel_distortion[MaxPixelChannels]);
1425 break;
1426 }
1427 case AbsoluteErrorMetric:
1428 case DotProductCorrelationErrorMetric:
1429 case NormalizedCrossCorrelationErrorMetric:
1430 case PeakSignalToNoiseRatioErrorMetric:
1431 case PerceptualHashErrorMetric:
1432 case PhaseCorrelationErrorMetric:
1433 case StructuralSimilarityErrorMetric:
1434 case StructuralDissimilarityErrorMetric:
1435 {
1436 switch (image->colorspace)
1437 {
1438 case RGBColorspace:
1439 default:
1440 {
1441 (void) FormatLocaleFile(stderr," red: %.*g\n",
1442 GetMagickPrecision(),channel_distortion[RedPixelChannel]);
1443 (void) FormatLocaleFile(stderr," green: %.*g\n",
1444 GetMagickPrecision(),channel_distortion[GreenPixelChannel]);
1445 (void) FormatLocaleFile(stderr," blue: %.*g\n",
1446 GetMagickPrecision(),channel_distortion[BluePixelChannel]);
1447 if (image->alpha_trait != UndefinedPixelTrait)
1448 (void) FormatLocaleFile(stderr," alpha: %.*g\n",
1449 GetMagickPrecision(),
1450 channel_distortion[AlphaPixelChannel]);
1451 break;
1452 }
1453 case CMYKColorspace:
1454 {
1455 (void) FormatLocaleFile(stderr," cyan: %.*g\n",
1456 GetMagickPrecision(),channel_distortion[CyanPixelChannel]);
1457 (void) FormatLocaleFile(stderr," magenta: %.*g\n",
1458 GetMagickPrecision(),
1459 channel_distortion[MagentaPixelChannel]);
1460 (void) FormatLocaleFile(stderr," yellow: %.*g\n",
1461 GetMagickPrecision(),
1462 channel_distortion[YellowPixelChannel]);
1463 (void) FormatLocaleFile(stderr," black: %.*g\n",
1464 GetMagickPrecision(),
1465 channel_distortion[BlackPixelChannel]);
1466 if (image->alpha_trait != UndefinedPixelTrait)
1467 (void) FormatLocaleFile(stderr," alpha: %.*g\n",
1468 GetMagickPrecision(),
1469 channel_distortion[AlphaPixelChannel]);
1470 break;
1471 }
1472 case LinearGRAYColorspace:
1473 case GRAYColorspace:
1474 {
1475 (void) FormatLocaleFile(stderr," gray: %.*g\n",
1476 GetMagickPrecision(),channel_distortion[GrayPixelChannel]);
1477 if (image->alpha_trait != UndefinedPixelTrait)
1478 (void) FormatLocaleFile(stderr," alpha: %.*g\n",
1479 GetMagickPrecision(),
1480 channel_distortion[AlphaPixelChannel]);
1481 break;
1482 }
1483 }
1484 (void) FormatLocaleFile(stderr," all: %.*g\n",
1485 GetMagickPrecision(),channel_distortion[MaxPixelChannels]);
1486 break;
1487 }
1488 case MeanErrorPerPixelErrorMetric:
1489 {
1490 (void) FormatLocaleFile(stderr," %.*g (%.*g, %.*g)\n",
1491 GetMagickPrecision(),channel_distortion[MaxPixelChannels],
1492 GetMagickPrecision(),image->error.normalized_mean_error,
1493 GetMagickPrecision(),image->error.normalized_maximum_error);
1494 break;
1495 }
1496 case UndefinedErrorMetric:
1497 break;
1498 }
1499 channel_distortion=(double *) RelinquishMagickMemory(
1500 channel_distortion);
1501 if (subimage_search != MagickFalse)
1502 {
1503 (void) FormatLocaleFile(stderr," Offset: %.20g,%.20g\n",(double)
1504 difference_image->page.x,(double) difference_image->page.y);
1505 (void) FormatLocaleFile(stderr," Similarity metric: %*g\n",
1506 GetMagickPrecision(),similarity_metric);
1507 }
1508 }
1509 (void) ResetImagePage(difference_image,"0x0+0+0");
1510 if (difference_image->next != (Image *) NULL)
1511 (void) ResetImagePage(difference_image->next,"0x0+0+0");
1512 status&=(MagickStatusType) WriteImages(image_info,difference_image,
1513 argv[argc-1],exception);
1514 if ((metadata != (char **) NULL) && (format != (char *) NULL))
1515 {
1516 char
1517 *text;
1518
1519 text=InterpretImageProperties(image_info,difference_image,format,
1520 exception);
1521 if (text == (char *) NULL)
1522 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
1523 GetExceptionMessage(errno));
1524 (void) ConcatenateString(&(*metadata),text);
1525 text=DestroyString(text);
1526 }
1527 difference_image=DestroyImageList(difference_image);
1528 }
1529 DestroyCompare();
1530 if (similar == MagickFalse)
1531 (void) SetImageOption(image_info,"compare:dissimilar","true");
1532 return(status != 0 ? MagickTrue : MagickFalse);
1533}