blob: 674c792298db47382edf54339ab7632833d22b74 [file] [log] [blame]
cristy3ed852e2009-09-05 21:47:34 +00001/*
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% John Cristy %
17% December 2003 %
18% %
19% %
cristy16af1cb2009-12-11 21:38:29 +000020% Copyright 1999-2010 ImageMagick Studio LLC, a non-profit organization %
cristy3ed852e2009-09-05 21:47:34 +000021% 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% http://www.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
cristy6a917d92009-10-06 19:23:54 +000037% difference between an image and its reconstruction.
cristy3ed852e2009-09-05 21:47:34 +000038%
39*/
40
41/*
42 Include declarations.
43*/
44#include "wand/studio.h"
45#include "wand/MagickWand.h"
46#include "wand/mogrify-private.h"
cristyf2f27272009-12-17 14:48:46 +000047#include "magick/string-private.h"
cristy3ed852e2009-09-05 21:47:34 +000048
49/*
50%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
51% %
52% %
53% %
54% C o m p a r e I m a g e C o m m a n d %
55% %
56% %
57% %
58%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
59%
60% CompareImageCommand() compares two images and returns the difference between
61% them as a distortion metric and as a new image visually annotating their
62% differences.
63%
64% The format of the CompareImageCommand method is:
65%
66% MagickBooleanType CompareImageCommand(ImageInfo *image_info,int argc,
67% char **argv,char **metadata,ExceptionInfo *exception)
68%
69% A description of each parameter follows:
70%
71% o image_info: the image info.
72%
73% o argc: the number of elements in the argument vector.
74%
75% o argv: A text array containing the command line arguments.
76%
77% o metadata: any metadata is returned here.
78%
79% o exception: return any errors or warnings in this structure.
80%
81*/
82
83static MagickBooleanType CompareUsage(void)
84{
85 const char
86 **p;
87
88 static const char
89 *miscellaneous[]=
90 {
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",
95 (char *) NULL
96 },
97 *settings[]=
98 {
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 RMSE 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 "-passphrase filename get the passphrase from this file",
128 "-profile filename add, delete, or apply an image profile",
129 "-quality value JPEG/MIFF/PNG compression level",
130 "-quiet suppress all warning messages",
131 "-quantize colorspace reduce colors in this colorspace",
132 "-regard-warnings pay attention to warning messages",
133 "-respect-parentheses settings remain in effect until parenthesis boundary",
134 "-sampling-factor geometry",
135 " horizontal and vertical sampling factor",
136 "-seed value seed a new sequence of pseudo-random numbers",
137 "-set attribute value set an image attribute",
138 "-quality value JPEG/MIFF/PNG compression level",
139 "-size geometry width and height of image",
140 "-transparent-color color",
141 " transparent color",
142 "-type type image type",
143 "-verbose print detailed information about the image",
144 "-version print version information",
145 "-virtual-pixel method",
146 " virtual pixel access method",
147 (char *) NULL
148 };
149
150 (void) printf("Version: %s\n",GetMagickVersion((unsigned long *) NULL));
cristy610b2e22009-10-22 14:59:43 +0000151 (void) printf("Copyright: %s\n",GetMagickCopyright());
152 (void) printf("Features: %s\n\n",GetMagickFeatures());
cristy3ed852e2009-09-05 21:47:34 +0000153 (void) printf("Usage: %s [options ...] image reconstruct difference\n",
154 GetClientName());
155 (void) printf("\nImage Settings:\n");
156 for (p=settings; *p != (char *) NULL; p++)
157 (void) printf(" %s\n",*p);
158 (void) printf("\nMiscellaneous Options:\n");
159 for (p=miscellaneous; *p != (char *) NULL; p++)
160 (void) printf(" %s\n",*p);
161 (void) printf(
162 "\nBy default, the image format of `file' is determined by its magic\n");
163 (void) printf(
164 "number. To specify a particular image format, precede the filename\n");
165 (void) printf(
166 "with an image format name and a colon (i.e. ps:image) or specify the\n");
167 (void) printf(
168 "image type as the filename suffix (i.e. image.ps). Specify 'file' as\n");
169 (void) printf("'-' for standard input or output.\n");
170 return(MagickFalse);
171}
172
173WandExport MagickBooleanType CompareImageCommand(ImageInfo *image_info,
174 int argc,char **argv,char **metadata,ExceptionInfo *exception)
175{
cristy2b70e802009-11-23 19:40:25 +0000176#define DefaultDissimilarityThreshold 0.31830988618379067154
cristy3ed852e2009-09-05 21:47:34 +0000177#define DestroyCompare() \
178{ \
179 if (similarity_image != (Image *) NULL) \
180 similarity_image=DestroyImageList(similarity_image); \
181 if (difference_image != (Image *) NULL) \
182 difference_image=DestroyImageList(difference_image); \
183 DestroyImageStack(); \
184 for (i=0; i < (long) argc; i++) \
185 argv[i]=DestroyString(argv[i]); \
186 argv=(char **) RelinquishMagickMemory(argv); \
187}
188#define ThrowCompareException(asperity,tag,option) \
189{ \
190 if (exception->severity < (asperity)) \
191 (void) ThrowMagickException(exception,GetMagickModule(),asperity,tag, \
192 "`%s'",option); \
193 DestroyCompare(); \
194 return(MagickFalse); \
195}
196#define ThrowCompareInvalidArgumentException(option,argument) \
197{ \
198 (void) ThrowMagickException(exception,GetMagickModule(),OptionError, \
199 "InvalidArgument","`%s': %s",option,argument); \
200 DestroyCompare(); \
201 return(MagickFalse); \
202}
203
204 char
205 *filename,
206 *option;
207
208 const char
209 *format;
210
211 ChannelType
212 channels;
213
214 double
215 dissimilarity_threshold,
216 distortion,
217 similarity_metric;
218
219 Image
220 *difference_image,
221 *image,
222 *reconstruct_image,
223 *similarity_image;
224
225 ImageStack
226 image_stack[MaxImageStackDepth+1];
227
228 long
229 j,
230 k;
231
232 MagickBooleanType
233 fire,
234 pend;
235
236 MagickStatusType
237 status;
238
239 MetricType
240 metric;
241
242 RectangleInfo
243 offset;
244
245 register long
246 i;
247
248 /*
249 Set defaults.
250 */
251 assert(image_info != (ImageInfo *) NULL);
252 assert(image_info->signature == MagickSignature);
253 if (image_info->debug != MagickFalse)
254 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
255 assert(exception != (ExceptionInfo *) NULL);
256 if (argc == 2)
257 {
258 option=argv[1];
259 if ((LocaleCompare("version",option+1) == 0) ||
260 (LocaleCompare("-version",option+1) == 0))
261 {
262 (void) fprintf(stdout,"Version: %s\n",
263 GetMagickVersion((unsigned long *) NULL));
cristy3980b0d2009-10-25 14:37:13 +0000264 (void) fprintf(stdout,"Copyright: %s\n",GetMagickCopyright());
cristy104cea82009-10-25 02:26:51 +0000265 (void) fprintf(stdout,"Features: %s\n\n",GetMagickFeatures());
cristy3ed852e2009-09-05 21:47:34 +0000266 return(MagickFalse);
267 }
268 }
269 if (argc < 3)
cristy13e61a12010-02-04 20:19:00 +0000270 return(CompareUsage());
cristy3ed852e2009-09-05 21:47:34 +0000271 channels=AllChannels;
272 difference_image=NewImageList();
273 similarity_image=NewImageList();
274 dissimilarity_threshold=DefaultDissimilarityThreshold;
275 distortion=0.0;
276 format=(char *) NULL;
277 j=1;
278 k=0;
279 metric=UndefinedMetric;
280 NewImageStack();
281 option=(char *) NULL;
282 pend=MagickFalse;
283 reconstruct_image=NewImageList();
284 status=MagickTrue;
285 /*
286 Compare an image.
287 */
288 ReadCommandlLine(argc,&argv);
289 status=ExpandFilenames(&argc,&argv);
290 if (status == MagickFalse)
291 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
292 GetExceptionMessage(errno));
293 for (i=1; i < (long) (argc-1); i++)
294 {
295 option=argv[i];
296 if (LocaleCompare(option,"(") == 0)
297 {
298 FireImageStack(MagickTrue,MagickTrue,pend);
299 if (k == MaxImageStackDepth)
300 ThrowCompareException(OptionError,"ParenthesisNestedTooDeeply",
301 option);
302 PushImageStack();
303 continue;
304 }
305 if (LocaleCompare(option,")") == 0)
306 {
307 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
308 if (k == 0)
309 ThrowCompareException(OptionError,"UnableToParseExpression",option);
310 PopImageStack();
311 continue;
312 }
313 if (IsMagickOption(option) == MagickFalse)
314 {
315 Image
316 *images;
317
318 /*
319 Read input image.
320 */
321 FireImageStack(MagickFalse,MagickFalse,pend);
322 filename=argv[i];
323 if ((LocaleCompare(filename,"--") == 0) && (i < (argc-1)))
324 filename=argv[++i];
325 (void) CopyMagickString(image_info->filename,filename,MaxTextExtent);
326 images=ReadImages(image_info,exception);
327 status&=(images != (Image *) NULL) &&
328 (exception->severity < ErrorException);
329 if (images == (Image *) NULL)
330 continue;
331 AppendImageStack(images);
332 continue;
333 }
334 pend=image != (Image *) NULL ? MagickTrue : MagickFalse;
335 switch (*(option+1))
336 {
337 case 'a':
338 {
339 if (LocaleCompare("alpha",option+1) == 0)
340 {
341 long
342 type;
343
344 if (*option == '+')
345 break;
346 i++;
347 if (i == (long) argc)
348 ThrowCompareException(OptionError,"MissingArgument",option);
349 type=ParseMagickOption(MagickAlphaOptions,MagickFalse,argv[i]);
350 if (type < 0)
351 ThrowCompareException(OptionError,"UnrecognizedAlphaChannelType",
352 argv[i]);
353 break;
354 }
355 if (LocaleCompare("authenticate",option+1) == 0)
356 {
357 if (*option == '+')
358 break;
359 i++;
360 if (i == (long) argc)
361 ThrowCompareException(OptionError,"MissingArgument",option);
362 break;
363 }
364 ThrowCompareException(OptionError,"UnrecognizedOption",option);
365 }
366 case 'c':
367 {
368 if (LocaleCompare("cache",option+1) == 0)
369 {
370 if (*option == '+')
371 break;
372 i++;
373 if (i == (long) argc)
374 ThrowCompareException(OptionError,"MissingArgument",option);
375 if (IsGeometry(argv[i]) == MagickFalse)
376 ThrowCompareInvalidArgumentException(option,argv[i]);
377 break;
378 }
379 if (LocaleCompare("channel",option+1) == 0)
380 {
381 long
382 channel;
383
384 if (*option == '+')
385 break;
386 i++;
387 if (i == (long) (argc-1))
388 ThrowCompareException(OptionError,"MissingArgument",option);
389 channel=ParseChannelOption(argv[i]);
390 if (channel < 0)
391 ThrowCompareException(OptionError,"UnrecognizedChannelType",
392 argv[i]);
393 channels=(ChannelType) channel;
394 break;
395 }
396 if (LocaleCompare("colorspace",option+1) == 0)
397 {
398 long
399 colorspace;
400
401 if (*option == '+')
402 break;
403 i++;
404 if (i == (long) (argc-1))
405 ThrowCompareException(OptionError,"MissingArgument",option);
406 colorspace=ParseMagickOption(MagickColorspaceOptions,MagickFalse,
407 argv[i]);
408 if (colorspace < 0)
409 ThrowCompareException(OptionError,"UnrecognizedColorspace",
410 argv[i]);
411 break;
412 }
413 if (LocaleCompare("compose",option+1) == 0)
414 {
415 long
416 compose;
417
418 if (*option == '+')
419 break;
420 i++;
421 if (i == (long) argc)
422 ThrowCompareException(OptionError,"MissingArgument",option);
423 compose=ParseMagickOption(MagickComposeOptions,MagickFalse,
424 argv[i]);
425 if (compose < 0)
426 ThrowCompareException(OptionError,"UnrecognizedComposeOperator",
427 argv[i]);
428 break;
429 }
430 if (LocaleCompare("compress",option+1) == 0)
431 {
432 long
433 compress;
434
435 if (*option == '+')
436 break;
437 i++;
438 if (i == (long) (argc-1))
439 ThrowCompareException(OptionError,"MissingArgument",option);
440 compress=ParseMagickOption(MagickCompressOptions,MagickFalse,
441 argv[i]);
442 if (compress < 0)
443 ThrowCompareException(OptionError,"UnrecognizedImageCompression",
444 argv[i]);
445 break;
446 }
cristy22879752009-10-25 23:55:40 +0000447 if (LocaleCompare("concurrent",option+1) == 0)
448 break;
cristy3ed852e2009-09-05 21:47:34 +0000449 ThrowCompareException(OptionError,"UnrecognizedOption",option)
450 }
451 case 'd':
452 {
453 if (LocaleCompare("debug",option+1) == 0)
454 {
455 LogEventType
456 event_mask;
457
458 if (*option == '+')
459 break;
460 i++;
461 if (i == (long) argc)
462 ThrowCompareException(OptionError,"MissingArgument",option);
463 event_mask=SetLogEventMask(argv[i]);
464 if (event_mask == UndefinedEvents)
465 ThrowCompareException(OptionError,"UnrecognizedEventType",
466 argv[i]);
467 break;
468 }
469 if (LocaleCompare("decipher",option+1) == 0)
470 {
471 if (*option == '+')
472 break;
473 i++;
474 if (i == (long) (argc-1))
475 ThrowCompareException(OptionError,"MissingArgument",option);
476 break;
477 }
478 if (LocaleCompare("define",option+1) == 0)
479 {
480 i++;
481 if (i == (long) argc)
482 ThrowCompareException(OptionError,"MissingArgument",option);
483 if (*option == '+')
484 {
485 const char
486 *define;
487
488 define=GetImageOption(image_info,argv[i]);
489 if (define == (const char *) NULL)
490 ThrowCompareException(OptionError,"NoSuchOption",argv[i]);
491 break;
492 }
493 break;
494 }
495 if (LocaleCompare("density",option+1) == 0)
496 {
497 if (*option == '+')
498 break;
499 i++;
500 if (i == (long) argc)
501 ThrowCompareException(OptionError,"MissingArgument",option);
502 if (IsGeometry(argv[i]) == MagickFalse)
503 ThrowCompareInvalidArgumentException(option,argv[i]);
504 break;
505 }
506 if (LocaleCompare("depth",option+1) == 0)
507 {
508 if (*option == '+')
509 break;
510 i++;
511 if (i == (long) argc)
512 ThrowCompareException(OptionError,"MissingArgument",option);
513 if (IsGeometry(argv[i]) == MagickFalse)
514 ThrowCompareInvalidArgumentException(option,argv[i]);
515 break;
516 }
517 if (LocaleCompare("dissimilarity-threshold",option+1) == 0)
518 {
519 if (*option == '+')
520 break;
521 i++;
522 if (i == (long) argc)
523 ThrowCompareException(OptionError,"MissingArgument",option);
524 if (IsGeometry(argv[i]) == MagickFalse)
525 ThrowCompareInvalidArgumentException(option,argv[i]);
526 if (*option == '+')
527 dissimilarity_threshold=DefaultDissimilarityThreshold;
528 else
cristyf2f27272009-12-17 14:48:46 +0000529 dissimilarity_threshold=StringToDouble(argv[i]);
cristy3ed852e2009-09-05 21:47:34 +0000530 break;
531 }
cristy22879752009-10-25 23:55:40 +0000532 if (LocaleCompare("duration",option+1) == 0)
533 {
534 if (*option == '+')
535 break;
536 i++;
537 if (i == (long) (argc-1))
538 ThrowCompareException(OptionError,"MissingArgument",option);
539 if (IsGeometry(argv[i]) == MagickFalse)
540 ThrowCompareInvalidArgumentException(option,argv[i]);
541 break;
542 }
cristy3ed852e2009-09-05 21:47:34 +0000543 ThrowCompareException(OptionError,"UnrecognizedOption",option)
544 }
545 case 'e':
546 {
547 if (LocaleCompare("encipher",option+1) == 0)
548 {
549 if (*option == '+')
550 break;
551 i++;
552 if (i == (long) (argc-1))
553 ThrowCompareException(OptionError,"MissingArgument",option);
554 break;
555 }
556 if (LocaleCompare("extract",option+1) == 0)
557 {
558 if (*option == '+')
559 break;
560 i++;
561 if (i == (long) (argc-1))
562 ThrowCompareException(OptionError,"MissingArgument",option);
563 if (IsGeometry(argv[i]) == MagickFalse)
564 ThrowCompareInvalidArgumentException(option,argv[i]);
565 break;
566 }
567 ThrowCompareException(OptionError,"UnrecognizedOption",option)
568 }
569 case 'f':
570 {
571 if (LocaleCompare("format",option+1) == 0)
572 {
573 if (*option == '+')
574 break;
575 i++;
576 if (i == (long) argc)
577 ThrowCompareException(OptionError,"MissingArgument",option);
578 format=argv[i];
579 break;
580 }
581 if (LocaleCompare("fuzz",option+1) == 0)
582 {
583 if (*option == '+')
584 break;
585 i++;
586 if (i == (long) (argc-1))
587 ThrowCompareException(OptionError,"MissingArgument",option);
588 if (IsGeometry(argv[i]) == MagickFalse)
589 ThrowCompareInvalidArgumentException(option,argv[i]);
590 break;
591 }
592 ThrowCompareException(OptionError,"UnrecognizedOption",option)
593 }
594 case 'h':
595 {
596 if ((LocaleCompare("help",option+1) == 0) ||
597 (LocaleCompare("-help",option+1) == 0))
598 return(CompareUsage());
599 if (LocaleCompare("highlight-color",option+1) == 0)
600 {
601 if (*option == '+')
602 break;
603 i++;
604 if (i == (long) (argc-1))
605 ThrowCompareException(OptionError,"MissingArgument",option);
606 break;
607 }
608 ThrowCompareException(OptionError,"UnrecognizedOption",option)
609 }
610 case 'i':
611 {
612 if (LocaleCompare("identify",option+1) == 0)
613 break;
614 if (LocaleCompare("interlace",option+1) == 0)
615 {
616 long
617 interlace;
618
619 if (*option == '+')
620 break;
621 i++;
622 if (i == (long) argc)
623 ThrowCompareException(OptionError,"MissingArgument",option);
624 interlace=ParseMagickOption(MagickInterlaceOptions,MagickFalse,
625 argv[i]);
626 if (interlace < 0)
627 ThrowCompareException(OptionError,"UnrecognizedInterlaceType",
628 argv[i]);
629 break;
630 }
631 ThrowCompareException(OptionError,"UnrecognizedOption",option)
632 }
633 case 'l':
634 {
635 if (LocaleCompare("limit",option+1) == 0)
636 {
637 char
638 *p;
639
640 double
641 value;
642
643 long
644 resource;
645
646 if (*option == '+')
647 break;
648 i++;
649 if (i == (long) argc)
650 ThrowCompareException(OptionError,"MissingArgument",option);
651 resource=ParseMagickOption(MagickResourceOptions,MagickFalse,
652 argv[i]);
653 if (resource < 0)
654 ThrowCompareException(OptionError,"UnrecognizedResourceType",
655 argv[i]);
656 i++;
657 if (i == (long) argc)
658 ThrowCompareException(OptionError,"MissingArgument",option);
659 value=strtod(argv[i],&p);
660 if ((p == argv[i]) && (LocaleCompare("unlimited",argv[i]) != 0))
661 ThrowCompareInvalidArgumentException(option,argv[i]);
662 break;
663 }
664 if (LocaleCompare("list",option+1) == 0)
665 {
666 long
667 list;
668
669 if (*option == '+')
670 break;
671 i++;
672 if (i == (long) argc)
673 ThrowCompareException(OptionError,"MissingArgument",option);
674 list=ParseMagickOption(MagickListOptions,MagickFalse,argv[i]);
675 if (list < 0)
676 ThrowCompareException(OptionError,"UnrecognizedListType",argv[i]);
cristyaeb2cbc2010-05-07 13:28:58 +0000677 status=MogrifyImageInfo(image_info,(int) (i-j+1),(const char **)
cristy3ed852e2009-09-05 21:47:34 +0000678 argv+j,exception);
679 DestroyCompare();
cristyaeb2cbc2010-05-07 13:28:58 +0000680 return(status != 0 ? MagickFalse : MagickTrue);
cristy3ed852e2009-09-05 21:47:34 +0000681 }
682 if (LocaleCompare("log",option+1) == 0)
683 {
684 if (*option == '+')
685 break;
686 i++;
687 if ((i == (long) argc) || (strchr(argv[i],'%') == (char *) NULL))
688 ThrowCompareException(OptionError,"MissingArgument",option);
689 break;
690 }
691 if (LocaleCompare("lowlight-color",option+1) == 0)
692 {
693 if (*option == '+')
694 break;
695 i++;
696 if (i == (long) (argc-1))
697 ThrowCompareException(OptionError,"MissingArgument",option);
698 break;
699 }
700 ThrowCompareException(OptionError,"UnrecognizedOption",option)
701 }
702 case 'm':
703 {
704 if (LocaleCompare("matte",option+1) == 0)
705 break;
706 if (LocaleCompare("metric",option+1) == 0)
707 {
708 long
709 type;
710
711 if (*option == '+')
712 break;
713 i++;
714 if (i == (long) argc)
715 ThrowCompareException(OptionError,"MissingArgument",option);
716 type=ParseMagickOption(MagickMetricOptions,MagickTrue,argv[i]);
717 if (type < 0)
718 ThrowCompareException(OptionError,"UnrecognizedMetricType",
719 argv[i]);
720 metric=(MetricType) type;
721 break;
722 }
723 if (LocaleCompare("monitor",option+1) == 0)
724 break;
725 ThrowCompareException(OptionError,"UnrecognizedOption",option)
726 }
727 case 'p':
728 {
729 if (LocaleCompare("passphrase",option+1) == 0)
730 {
731 if (*option == '+')
732 break;
733 i++;
734 if (i == (long) argc)
735 ThrowCompareException(OptionError,"MissingArgument",option);
736 break;
737 }
738 if (LocaleCompare("profile",option+1) == 0)
739 {
740 i++;
741 if (i == (long) (argc-1))
742 ThrowCompareException(OptionError,"MissingArgument",option);
743 break;
744 }
745 ThrowCompareException(OptionError,"UnrecognizedOption",option)
746 }
747 case 'q':
748 {
749 if (LocaleCompare("quality",option+1) == 0)
750 {
751 if (*option == '+')
752 break;
753 i++;
754 if (i == (long) (argc-1))
755 ThrowCompareException(OptionError,"MissingArgument",option);
756 if (IsGeometry(argv[i]) == MagickFalse)
757 ThrowCompareInvalidArgumentException(option,argv[i]);
758 break;
759 }
760 if (LocaleCompare("quantize",option+1) == 0)
761 {
762 long
763 colorspace;
764
765 if (*option == '+')
766 break;
767 i++;
768 if (i == (long) (argc-1))
769 ThrowCompareException(OptionError,"MissingArgument",option);
770 colorspace=ParseMagickOption(MagickColorspaceOptions,
771 MagickFalse,argv[i]);
772 if (colorspace < 0)
773 ThrowCompareException(OptionError,"UnrecognizedColorspace",
774 argv[i]);
775 break;
776 }
777 if (LocaleCompare("quiet",option+1) == 0)
778 break;
779 ThrowCompareException(OptionError,"UnrecognizedOption",option)
780 }
781 case 'r':
782 {
783 if (LocaleCompare("regard-warnings",option+1) == 0)
784 break;
785 if (LocaleNCompare("respect-parentheses",option+1,17) == 0)
786 {
787 respect_parenthesis=(*option == '-') ? MagickTrue : MagickFalse;
788 break;
789 }
790 ThrowCompareException(OptionError,"UnrecognizedOption",option)
791 }
792 case 's':
793 {
794 if (LocaleCompare("sampling-factor",option+1) == 0)
795 {
796 if (*option == '+')
797 break;
798 i++;
799 if (i == (long) argc)
800 ThrowCompareException(OptionError,"MissingArgument",option);
801 if (IsGeometry(argv[i]) == MagickFalse)
802 ThrowCompareInvalidArgumentException(option,argv[i]);
803 break;
804 }
805 if (LocaleCompare("seed",option+1) == 0)
806 {
807 if (*option == '+')
808 break;
809 i++;
810 if (i == (long) (argc-1))
811 ThrowCompareException(OptionError,"MissingArgument",option);
812 if (IsGeometry(argv[i]) == MagickFalse)
813 ThrowCompareInvalidArgumentException(option,argv[i]);
814 break;
815 }
816 if (LocaleCompare("set",option+1) == 0)
817 {
818 i++;
819 if (i == (long) argc)
820 ThrowCompareException(OptionError,"MissingArgument",option);
821 if (*option == '+')
822 break;
823 i++;
824 if (i == (long) argc)
825 ThrowCompareException(OptionError,"MissingArgument",option);
826 break;
827 }
828 if (LocaleCompare("size",option+1) == 0)
829 {
830 if (*option == '+')
831 break;
832 i++;
833 if (i == (long) argc)
834 ThrowCompareException(OptionError,"MissingArgument",option);
835 if (IsGeometry(argv[i]) == MagickFalse)
836 ThrowCompareInvalidArgumentException(option,argv[i]);
837 break;
838 }
839 ThrowCompareException(OptionError,"UnrecognizedOption",option)
840 }
841 case 't':
842 {
843 if (LocaleCompare("transparent-color",option+1) == 0)
844 {
845 if (*option == '+')
846 break;
847 i++;
848 if (i == (long) (argc-1))
849 ThrowCompareException(OptionError,"MissingArgument",option);
850 break;
851 }
852 if (LocaleCompare("type",option+1) == 0)
853 {
854 long
855 type;
856
857 if (*option == '+')
858 break;
859 i++;
860 if (i == (long) argc)
861 ThrowCompareException(OptionError,"MissingArgument",option);
862 type=ParseMagickOption(MagickTypeOptions,MagickFalse,argv[i]);
863 if (type < 0)
864 ThrowCompareException(OptionError,"UnrecognizedImageType",
865 argv[i]);
866 break;
867 }
868 ThrowCompareException(OptionError,"UnrecognizedOption",option)
869 }
870 case 'v':
871 {
872 if (LocaleCompare("verbose",option+1) == 0)
873 break;
874 if ((LocaleCompare("version",option+1) == 0) ||
875 (LocaleCompare("-version",option+1) == 0))
876 {
877 (void) fprintf(stdout,"Version: %s\n",
878 GetMagickVersion((unsigned long *) NULL));
cristy3980b0d2009-10-25 14:37:13 +0000879 (void) fprintf(stdout,"Copyright: %s\n",GetMagickCopyright());
cristy104cea82009-10-25 02:26:51 +0000880 (void) fprintf(stdout,"Features: %s\n\n",GetMagickFeatures());
cristy3ed852e2009-09-05 21:47:34 +0000881 break;
882 }
883 if (LocaleCompare("virtual-pixel",option+1) == 0)
884 {
885 long
886 method;
887
888 if (*option == '+')
889 break;
890 i++;
891 if (i == (long) (argc-1))
892 ThrowCompareException(OptionError,"MissingArgument",option);
893 method=ParseMagickOption(MagickVirtualPixelOptions,MagickFalse,
894 argv[i]);
895 if (method < 0)
896 ThrowCompareException(OptionError,
897 "UnrecognizedVirtualPixelMethod",argv[i]);
898 break;
899 }
900 ThrowCompareException(OptionError,"UnrecognizedOption",option)
901 }
902 case '?':
903 break;
904 default:
905 ThrowCompareException(OptionError,"UnrecognizedOption",option)
906 }
907 fire=ParseMagickOption(MagickImageListOptions,MagickFalse,option+1) < 0 ?
908 MagickFalse : MagickTrue;
909 if (fire != MagickFalse)
910 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
911 }
912 if (k != 0)
913 ThrowCompareException(OptionError,"UnbalancedParenthesis",argv[i]);
914 if (i-- != (long) (argc-1))
915 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
916 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
917 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
918 FinalizeImageSettings(image_info,image,MagickTrue);
cristy94f20072009-09-12 02:12:36 +0000919 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
920 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
cristy3ed852e2009-09-05 21:47:34 +0000921 image=GetImageFromList(image,0);
922 reconstruct_image=GetImageFromList(image,1);
923 similarity_image=SimilarityImage(image,reconstruct_image,&offset,
924 &similarity_metric,exception);
925 if (similarity_metric > dissimilarity_threshold)
926 ThrowCompareException(ImageError,"ImagesTooDissimilar",image->filename);
927 if ((reconstruct_image->columns == image->columns) &&
928 (reconstruct_image->rows == image->rows))
929 difference_image=CompareImageChannels(image,reconstruct_image,channels,
930 metric,&distortion,exception);
931 else
932 if (similarity_image != (Image *) NULL)
933 {
934 Image
935 *composite_image;
936
937 /*
938 Determine if reconstructed image is a subimage of the image.
939 */
940 composite_image=CloneImage(image,0,0,MagickTrue,exception);
941 if (composite_image == (Image *) NULL)
942 difference_image=CompareImageChannels(image,reconstruct_image,
943 channels,metric,&distortion,exception);
944 else
945 {
946 (void) CompositeImage(composite_image,CopyCompositeOp,
947 reconstruct_image,offset.x,offset.y);
948 difference_image=CompareImageChannels(image,composite_image,
949 channels,metric,&distortion,exception);
950 if (difference_image != (Image *) NULL)
951 {
952 difference_image->page.x=offset.x;
953 difference_image->page.y=offset.y;
954 }
955 composite_image=DestroyImage(composite_image);
956 }
957 if (difference_image == (Image *) NULL)
958 similarity_image=DestroyImage(similarity_image);
959 else
960 {
961 AppendImageToList(&difference_image,similarity_image);
962 similarity_image=(Image *) NULL;
963 }
964 }
965 if (difference_image == (Image *) NULL)
966 status=0;
967 else
968 {
969 if (image_info->verbose != MagickFalse)
970 (void) IsImagesEqual(image,reconstruct_image);
971 if (*difference_image->magick == '\0')
972 (void) CopyMagickString(difference_image->magick,image->magick,
973 MaxTextExtent);
974 if (image_info->verbose == MagickFalse)
975 {
976 switch (metric)
977 {
978 case MeanAbsoluteErrorMetric:
979 case MeanSquaredErrorMetric:
980 case RootMeanSquaredErrorMetric:
981 case PeakAbsoluteErrorMetric:
982 {
cristye7f51092010-01-17 00:39:37 +0000983 (void) fprintf(stderr,"%g (%g)",QuantumRange*distortion,
cristy8cd5b312010-01-07 01:10:24 +0000984 (double) distortion);
cristy3ed852e2009-09-05 21:47:34 +0000985 if ((reconstruct_image->columns != image->columns) ||
986 (reconstruct_image->rows != image->rows))
987 (void) fprintf(stderr," @ %ld,%ld",difference_image->page.x,
988 difference_image->page.y);
989 (void) fprintf(stderr,"\n");
990 break;
991 }
992 case AbsoluteErrorMetric:
993 case PeakSignalToNoiseRatioMetric:
994 {
cristye7f51092010-01-17 00:39:37 +0000995 (void) fprintf(stderr,"%g",distortion);
cristy3ed852e2009-09-05 21:47:34 +0000996 if ((reconstruct_image->columns != image->columns) ||
997 (reconstruct_image->rows != image->rows))
998 (void) fprintf(stderr," @ %ld,%ld",difference_image->page.x,
999 difference_image->page.y);
1000 (void) fprintf(stderr,"\n");
1001 break;
1002 }
1003 case MeanErrorPerPixelMetric:
1004 {
cristye7f51092010-01-17 00:39:37 +00001005 (void) fprintf(stderr,"%g (%g, %g)",distortion,
cristy3ed852e2009-09-05 21:47:34 +00001006 image->error.normalized_mean_error,
1007 image->error.normalized_maximum_error);
1008 if ((reconstruct_image->columns != image->columns) ||
1009 (reconstruct_image->rows != image->rows))
1010 (void) fprintf(stderr," @ %ld,%ld",difference_image->page.x,
1011 difference_image->page.y);
1012 (void) fprintf(stderr,"\n");
1013 break;
1014 }
1015 case UndefinedMetric:
1016 break;
1017 }
1018 }
1019 else
1020 {
1021 double
1022 *channel_distortion;
1023
1024 channel_distortion=GetImageChannelDistortions(image,reconstruct_image,
1025 metric,&image->exception);
1026 (void) fprintf(stderr,"Image: %s\n",image->filename);
1027 if ((reconstruct_image->columns != image->columns) ||
1028 (reconstruct_image->rows != image->rows))
1029 (void) fprintf(stderr,"Offset: %ld,%ld\n",difference_image->page.x,
1030 difference_image->page.y);
1031 (void) fprintf(stderr," Channel distortion: %s\n",
1032 MagickOptionToMnemonic(MagickMetricOptions,(long) metric));
1033 switch (metric)
1034 {
1035 case MeanAbsoluteErrorMetric:
1036 case MeanSquaredErrorMetric:
1037 case RootMeanSquaredErrorMetric:
1038 case PeakAbsoluteErrorMetric:
1039 {
1040 switch (image->colorspace)
1041 {
1042 case RGBColorspace:
1043 default:
1044 {
cristye7f51092010-01-17 00:39:37 +00001045 (void) fprintf(stderr," red: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001046 QuantumRange*channel_distortion[RedChannel],
1047 channel_distortion[RedChannel]);
cristye7f51092010-01-17 00:39:37 +00001048 (void) fprintf(stderr," green: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001049 QuantumRange*channel_distortion[GreenChannel],
1050 channel_distortion[GreenChannel]);
cristye7f51092010-01-17 00:39:37 +00001051 (void) fprintf(stderr," blue: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001052 QuantumRange*channel_distortion[BlueChannel],
1053 channel_distortion[BlueChannel]);
1054 if (image->matte != MagickFalse)
cristye7f51092010-01-17 00:39:37 +00001055 (void) fprintf(stderr," alpha: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001056 QuantumRange*channel_distortion[OpacityChannel],
1057 channel_distortion[OpacityChannel]);
1058 break;
1059 }
1060 case CMYKColorspace:
1061 {
cristye7f51092010-01-17 00:39:37 +00001062 (void) fprintf(stderr," cyan: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001063 QuantumRange*channel_distortion[CyanChannel],
1064 channel_distortion[CyanChannel]);
cristye7f51092010-01-17 00:39:37 +00001065 (void) fprintf(stderr," magenta: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001066 QuantumRange*channel_distortion[MagentaChannel],
1067 channel_distortion[MagentaChannel]);
cristye7f51092010-01-17 00:39:37 +00001068 (void) fprintf(stderr," yellow: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001069 QuantumRange*channel_distortion[YellowChannel],
1070 channel_distortion[YellowChannel]);
cristye7f51092010-01-17 00:39:37 +00001071 (void) fprintf(stderr," black: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001072 QuantumRange*channel_distortion[BlackChannel],
1073 channel_distortion[BlackChannel]);
1074 if (image->matte != MagickFalse)
cristye7f51092010-01-17 00:39:37 +00001075 (void) fprintf(stderr," alpha: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001076 QuantumRange*channel_distortion[OpacityChannel],
1077 channel_distortion[OpacityChannel]);
1078 break;
1079 }
1080 case GRAYColorspace:
1081 {
cristye7f51092010-01-17 00:39:37 +00001082 (void) fprintf(stderr," gray: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001083 QuantumRange*channel_distortion[GrayChannel],
1084 channel_distortion[GrayChannel]);
1085 if (image->matte != MagickFalse)
cristye7f51092010-01-17 00:39:37 +00001086 (void) fprintf(stderr," alpha: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001087 QuantumRange*channel_distortion[OpacityChannel],
1088 channel_distortion[OpacityChannel]);
1089 break;
1090 }
1091 }
cristye7f51092010-01-17 00:39:37 +00001092 (void) fprintf(stderr," all: %g (%g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001093 QuantumRange*channel_distortion[AllChannels],
1094 channel_distortion[AllChannels]);
1095 break;
1096 }
1097 case AbsoluteErrorMetric:
1098 case PeakSignalToNoiseRatioMetric:
1099 {
1100 switch (image->colorspace)
1101 {
1102 case RGBColorspace:
1103 default:
1104 {
cristye7f51092010-01-17 00:39:37 +00001105 (void) fprintf(stderr," red: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001106 channel_distortion[RedChannel]);
cristye7f51092010-01-17 00:39:37 +00001107 (void) fprintf(stderr," green: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001108 channel_distortion[GreenChannel]);
cristye7f51092010-01-17 00:39:37 +00001109 (void) fprintf(stderr," blue: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001110 channel_distortion[BlueChannel]);
1111 if (image->matte != MagickFalse)
cristye7f51092010-01-17 00:39:37 +00001112 (void) fprintf(stderr," alpha: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001113 channel_distortion[OpacityChannel]);
1114 break;
1115 }
1116 case CMYKColorspace:
1117 {
cristye7f51092010-01-17 00:39:37 +00001118 (void) fprintf(stderr," cyan: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001119 channel_distortion[CyanChannel]);
cristye7f51092010-01-17 00:39:37 +00001120 (void) fprintf(stderr," magenta: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001121 channel_distortion[MagentaChannel]);
cristye7f51092010-01-17 00:39:37 +00001122 (void) fprintf(stderr," yellow: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001123 channel_distortion[YellowChannel]);
cristye7f51092010-01-17 00:39:37 +00001124 (void) fprintf(stderr," black: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001125 channel_distortion[BlackChannel]);
1126 if (image->matte != MagickFalse)
cristye7f51092010-01-17 00:39:37 +00001127 (void) fprintf(stderr," alpha: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001128 channel_distortion[OpacityChannel]);
1129 break;
1130 }
1131 case GRAYColorspace:
1132 {
cristye7f51092010-01-17 00:39:37 +00001133 (void) fprintf(stderr," gray: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001134 channel_distortion[GrayChannel]);
1135 if (image->matte != MagickFalse)
cristye7f51092010-01-17 00:39:37 +00001136 (void) fprintf(stderr," alpha: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001137 channel_distortion[OpacityChannel]);
1138 break;
1139 }
1140 }
cristye7f51092010-01-17 00:39:37 +00001141 (void) fprintf(stderr," all: %g\n",
cristy3ed852e2009-09-05 21:47:34 +00001142 channel_distortion[AllChannels]);
1143 break;
1144 }
1145 case MeanErrorPerPixelMetric:
1146 {
cristye7f51092010-01-17 00:39:37 +00001147 (void) fprintf(stderr," %g (%g, %g)\n",
cristy3ed852e2009-09-05 21:47:34 +00001148 channel_distortion[AllChannels],
1149 image->error.normalized_mean_error,
1150 image->error.normalized_maximum_error);
1151 break;
1152 }
1153 case UndefinedMetric:
1154 break;
1155 }
1156 channel_distortion=(double *) RelinquishMagickMemory(
1157 channel_distortion);
1158 }
1159 status&=WriteImages(image_info,difference_image,argv[argc-1],exception);
1160 if ((metadata != (char **) NULL) && (format != (char *) NULL))
1161 {
1162 char
1163 *text;
1164
1165 text=InterpretImageProperties(image_info,difference_image,format);
1166 if (text == (char *) NULL)
1167 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
1168 GetExceptionMessage(errno));
1169 (void) ConcatenateString(&(*metadata),text);
1170 (void) ConcatenateString(&(*metadata),"\n");
1171 text=DestroyString(text);
1172 }
1173 difference_image=DestroyImageList(difference_image);
1174 }
1175 DestroyCompare();
1176 return(status != 0 ? MagickTrue : MagickFalse);
1177}