Gitiles
Code Review Sign In
gerrit-public.fairphone.software / platform / external / ImageMagick / 8cd5b3193212b4aebce08c4e7afbb66b09778029 / . / wand / compare.c
blob: d02c67c265d7e1f4ed0363e8daf4b62c76794342 [file] [log] [blame]
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]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% John Cristy %
17% December 2003 %
18% %
19% %
cristy16af1cb2009-12-11 21:38:29 +0000[diff] [blame]20% Copyright 1999-2010 ImageMagick Studio LLC, a non-profit organization %
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]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% 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 +0000[diff] [blame]37% difference between an image and its reconstruction.
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]38%
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 +0000[diff] [blame]47#include "magick/string-private.h"
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]48
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 +0000[diff] [blame]151 (void) printf("Copyright: %s\n",GetMagickCopyright());
152 (void) printf("Features: %s\n\n",GetMagickFeatures());
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]153 (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 +0000[diff] [blame]176#define DefaultDissimilarityThreshold 0.31830988618379067154
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]177#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 +0000[diff] [blame]264 (void) fprintf(stdout,"Copyright: %s\n",GetMagickCopyright());
cristy104cea82009-10-25 02:26:51 +0000[diff] [blame]265 (void) fprintf(stdout,"Features: %s\n\n",GetMagickFeatures());
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]266 return(MagickFalse);
267 }
268 }
269 if (argc < 3)
270 {
271 (void) CompareUsage();
272 return(MagickTrue);
273 }
274 channels=AllChannels;
275 difference_image=NewImageList();
276 similarity_image=NewImageList();
277 dissimilarity_threshold=DefaultDissimilarityThreshold;
278 distortion=0.0;
279 format=(char *) NULL;
280 j=1;
281 k=0;
282 metric=UndefinedMetric;
283 NewImageStack();
284 option=(char *) NULL;
285 pend=MagickFalse;
286 reconstruct_image=NewImageList();
287 status=MagickTrue;
288 /*
289 Compare an image.
290 */
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 < (long) (argc-1); i++)
297 {
298 option=argv[i];
299 if (LocaleCompare(option,"(") == 0)
300 {
301 FireImageStack(MagickTrue,MagickTrue,pend);
302 if (k == MaxImageStackDepth)
303 ThrowCompareException(OptionError,"ParenthesisNestedTooDeeply",
304 option);
305 PushImageStack();
306 continue;
307 }
308 if (LocaleCompare(option,")") == 0)
309 {
310 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
311 if (k == 0)
312 ThrowCompareException(OptionError,"UnableToParseExpression",option);
313 PopImageStack();
314 continue;
315 }
316 if (IsMagickOption(option) == MagickFalse)
317 {
318 Image
319 *images;
320
321 /*
322 Read input image.
323 */
324 FireImageStack(MagickFalse,MagickFalse,pend);
325 filename=argv[i];
326 if ((LocaleCompare(filename,"--") == 0) && (i < (argc-1)))
327 filename=argv[++i];
328 (void) CopyMagickString(image_info->filename,filename,MaxTextExtent);
329 images=ReadImages(image_info,exception);
330 status&=(images != (Image *) NULL) &&
331 (exception->severity < ErrorException);
332 if (images == (Image *) NULL)
333 continue;
334 AppendImageStack(images);
335 continue;
336 }
337 pend=image != (Image *) NULL ? MagickTrue : MagickFalse;
338 switch (*(option+1))
339 {
340 case 'a':
341 {
342 if (LocaleCompare("alpha",option+1) == 0)
343 {
344 long
345 type;
346
347 if (*option == '+')
348 break;
349 i++;
350 if (i == (long) argc)
351 ThrowCompareException(OptionError,"MissingArgument",option);
352 type=ParseMagickOption(MagickAlphaOptions,MagickFalse,argv[i]);
353 if (type < 0)
354 ThrowCompareException(OptionError,"UnrecognizedAlphaChannelType",
355 argv[i]);
356 break;
357 }
358 if (LocaleCompare("authenticate",option+1) == 0)
359 {
360 if (*option == '+')
361 break;
362 i++;
363 if (i == (long) argc)
364 ThrowCompareException(OptionError,"MissingArgument",option);
365 break;
366 }
367 ThrowCompareException(OptionError,"UnrecognizedOption",option);
368 }
369 case 'c':
370 {
371 if (LocaleCompare("cache",option+1) == 0)
372 {
373 if (*option == '+')
374 break;
375 i++;
376 if (i == (long) argc)
377 ThrowCompareException(OptionError,"MissingArgument",option);
378 if (IsGeometry(argv[i]) == MagickFalse)
379 ThrowCompareInvalidArgumentException(option,argv[i]);
380 break;
381 }
382 if (LocaleCompare("channel",option+1) == 0)
383 {
384 long
385 channel;
386
387 if (*option == '+')
388 break;
389 i++;
390 if (i == (long) (argc-1))
391 ThrowCompareException(OptionError,"MissingArgument",option);
392 channel=ParseChannelOption(argv[i]);
393 if (channel < 0)
394 ThrowCompareException(OptionError,"UnrecognizedChannelType",
395 argv[i]);
396 channels=(ChannelType) channel;
397 break;
398 }
399 if (LocaleCompare("colorspace",option+1) == 0)
400 {
401 long
402 colorspace;
403
404 if (*option == '+')
405 break;
406 i++;
407 if (i == (long) (argc-1))
408 ThrowCompareException(OptionError,"MissingArgument",option);
409 colorspace=ParseMagickOption(MagickColorspaceOptions,MagickFalse,
410 argv[i]);
411 if (colorspace < 0)
412 ThrowCompareException(OptionError,"UnrecognizedColorspace",
413 argv[i]);
414 break;
415 }
416 if (LocaleCompare("compose",option+1) == 0)
417 {
418 long
419 compose;
420
421 if (*option == '+')
422 break;
423 i++;
424 if (i == (long) argc)
425 ThrowCompareException(OptionError,"MissingArgument",option);
426 compose=ParseMagickOption(MagickComposeOptions,MagickFalse,
427 argv[i]);
428 if (compose < 0)
429 ThrowCompareException(OptionError,"UnrecognizedComposeOperator",
430 argv[i]);
431 break;
432 }
433 if (LocaleCompare("compress",option+1) == 0)
434 {
435 long
436 compress;
437
438 if (*option == '+')
439 break;
440 i++;
441 if (i == (long) (argc-1))
442 ThrowCompareException(OptionError,"MissingArgument",option);
443 compress=ParseMagickOption(MagickCompressOptions,MagickFalse,
444 argv[i]);
445 if (compress < 0)
446 ThrowCompareException(OptionError,"UnrecognizedImageCompression",
447 argv[i]);
448 break;
449 }
cristy22879752009-10-25 23:55:40 +0000[diff] [blame]450 if (LocaleCompare("concurrent",option+1) == 0)
451 break;
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]452 ThrowCompareException(OptionError,"UnrecognizedOption",option)
453 }
454 case 'd':
455 {
456 if (LocaleCompare("debug",option+1) == 0)
457 {
458 LogEventType
459 event_mask;
460
461 if (*option == '+')
462 break;
463 i++;
464 if (i == (long) argc)
465 ThrowCompareException(OptionError,"MissingArgument",option);
466 event_mask=SetLogEventMask(argv[i]);
467 if (event_mask == UndefinedEvents)
468 ThrowCompareException(OptionError,"UnrecognizedEventType",
469 argv[i]);
470 break;
471 }
472 if (LocaleCompare("decipher",option+1) == 0)
473 {
474 if (*option == '+')
475 break;
476 i++;
477 if (i == (long) (argc-1))
478 ThrowCompareException(OptionError,"MissingArgument",option);
479 break;
480 }
481 if (LocaleCompare("define",option+1) == 0)
482 {
483 i++;
484 if (i == (long) argc)
485 ThrowCompareException(OptionError,"MissingArgument",option);
486 if (*option == '+')
487 {
488 const char
489 *define;
490
491 define=GetImageOption(image_info,argv[i]);
492 if (define == (const char *) NULL)
493 ThrowCompareException(OptionError,"NoSuchOption",argv[i]);
494 break;
495 }
496 break;
497 }
498 if (LocaleCompare("density",option+1) == 0)
499 {
500 if (*option == '+')
501 break;
502 i++;
503 if (i == (long) argc)
504 ThrowCompareException(OptionError,"MissingArgument",option);
505 if (IsGeometry(argv[i]) == MagickFalse)
506 ThrowCompareInvalidArgumentException(option,argv[i]);
507 break;
508 }
509 if (LocaleCompare("depth",option+1) == 0)
510 {
511 if (*option == '+')
512 break;
513 i++;
514 if (i == (long) argc)
515 ThrowCompareException(OptionError,"MissingArgument",option);
516 if (IsGeometry(argv[i]) == MagickFalse)
517 ThrowCompareInvalidArgumentException(option,argv[i]);
518 break;
519 }
520 if (LocaleCompare("dissimilarity-threshold",option+1) == 0)
521 {
522 if (*option == '+')
523 break;
524 i++;
525 if (i == (long) argc)
526 ThrowCompareException(OptionError,"MissingArgument",option);
527 if (IsGeometry(argv[i]) == MagickFalse)
528 ThrowCompareInvalidArgumentException(option,argv[i]);
529 if (*option == '+')
530 dissimilarity_threshold=DefaultDissimilarityThreshold;
531 else
cristyf2f27272009-12-17 14:48:46 +0000[diff] [blame]532 dissimilarity_threshold=StringToDouble(argv[i]);
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]533 break;
534 }
cristy22879752009-10-25 23:55:40 +0000[diff] [blame]535 if (LocaleCompare("duration",option+1) == 0)
536 {
537 if (*option == '+')
538 break;
539 i++;
540 if (i == (long) (argc-1))
541 ThrowCompareException(OptionError,"MissingArgument",option);
542 if (IsGeometry(argv[i]) == MagickFalse)
543 ThrowCompareInvalidArgumentException(option,argv[i]);
544 break;
545 }
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]546 ThrowCompareException(OptionError,"UnrecognizedOption",option)
547 }
548 case 'e':
549 {
550 if (LocaleCompare("encipher",option+1) == 0)
551 {
552 if (*option == '+')
553 break;
554 i++;
555 if (i == (long) (argc-1))
556 ThrowCompareException(OptionError,"MissingArgument",option);
557 break;
558 }
559 if (LocaleCompare("extract",option+1) == 0)
560 {
561 if (*option == '+')
562 break;
563 i++;
564 if (i == (long) (argc-1))
565 ThrowCompareException(OptionError,"MissingArgument",option);
566 if (IsGeometry(argv[i]) == MagickFalse)
567 ThrowCompareInvalidArgumentException(option,argv[i]);
568 break;
569 }
570 ThrowCompareException(OptionError,"UnrecognizedOption",option)
571 }
572 case 'f':
573 {
574 if (LocaleCompare("format",option+1) == 0)
575 {
576 if (*option == '+')
577 break;
578 i++;
579 if (i == (long) argc)
580 ThrowCompareException(OptionError,"MissingArgument",option);
581 format=argv[i];
582 break;
583 }
584 if (LocaleCompare("fuzz",option+1) == 0)
585 {
586 if (*option == '+')
587 break;
588 i++;
589 if (i == (long) (argc-1))
590 ThrowCompareException(OptionError,"MissingArgument",option);
591 if (IsGeometry(argv[i]) == MagickFalse)
592 ThrowCompareInvalidArgumentException(option,argv[i]);
593 break;
594 }
595 ThrowCompareException(OptionError,"UnrecognizedOption",option)
596 }
597 case 'h':
598 {
599 if ((LocaleCompare("help",option+1) == 0) ||
600 (LocaleCompare("-help",option+1) == 0))
601 return(CompareUsage());
602 if (LocaleCompare("highlight-color",option+1) == 0)
603 {
604 if (*option == '+')
605 break;
606 i++;
607 if (i == (long) (argc-1))
608 ThrowCompareException(OptionError,"MissingArgument",option);
609 break;
610 }
611 ThrowCompareException(OptionError,"UnrecognizedOption",option)
612 }
613 case 'i':
614 {
615 if (LocaleCompare("identify",option+1) == 0)
616 break;
617 if (LocaleCompare("interlace",option+1) == 0)
618 {
619 long
620 interlace;
621
622 if (*option == '+')
623 break;
624 i++;
625 if (i == (long) argc)
626 ThrowCompareException(OptionError,"MissingArgument",option);
627 interlace=ParseMagickOption(MagickInterlaceOptions,MagickFalse,
628 argv[i]);
629 if (interlace < 0)
630 ThrowCompareException(OptionError,"UnrecognizedInterlaceType",
631 argv[i]);
632 break;
633 }
634 ThrowCompareException(OptionError,"UnrecognizedOption",option)
635 }
636 case 'l':
637 {
638 if (LocaleCompare("limit",option+1) == 0)
639 {
640 char
641 *p;
642
643 double
644 value;
645
646 long
647 resource;
648
649 if (*option == '+')
650 break;
651 i++;
652 if (i == (long) argc)
653 ThrowCompareException(OptionError,"MissingArgument",option);
654 resource=ParseMagickOption(MagickResourceOptions,MagickFalse,
655 argv[i]);
656 if (resource < 0)
657 ThrowCompareException(OptionError,"UnrecognizedResourceType",
658 argv[i]);
659 i++;
660 if (i == (long) argc)
661 ThrowCompareException(OptionError,"MissingArgument",option);
662 value=strtod(argv[i],&p);
663 if ((p == argv[i]) && (LocaleCompare("unlimited",argv[i]) != 0))
664 ThrowCompareInvalidArgumentException(option,argv[i]);
665 break;
666 }
667 if (LocaleCompare("list",option+1) == 0)
668 {
669 long
670 list;
671
672 if (*option == '+')
673 break;
674 i++;
675 if (i == (long) argc)
676 ThrowCompareException(OptionError,"MissingArgument",option);
677 list=ParseMagickOption(MagickListOptions,MagickFalse,argv[i]);
678 if (list < 0)
679 ThrowCompareException(OptionError,"UnrecognizedListType",argv[i]);
680 (void) MogrifyImageInfo(image_info,(int) (i-j+1),(const char **)
681 argv+j,exception);
682 DestroyCompare();
683 return(MagickTrue);
684 }
685 if (LocaleCompare("log",option+1) == 0)
686 {
687 if (*option == '+')
688 break;
689 i++;
690 if ((i == (long) argc) || (strchr(argv[i],'%') == (char *) NULL))
691 ThrowCompareException(OptionError,"MissingArgument",option);
692 break;
693 }
694 if (LocaleCompare("lowlight-color",option+1) == 0)
695 {
696 if (*option == '+')
697 break;
698 i++;
699 if (i == (long) (argc-1))
700 ThrowCompareException(OptionError,"MissingArgument",option);
701 break;
702 }
703 ThrowCompareException(OptionError,"UnrecognizedOption",option)
704 }
705 case 'm':
706 {
707 if (LocaleCompare("matte",option+1) == 0)
708 break;
709 if (LocaleCompare("metric",option+1) == 0)
710 {
711 long
712 type;
713
714 if (*option == '+')
715 break;
716 i++;
717 if (i == (long) argc)
718 ThrowCompareException(OptionError,"MissingArgument",option);
719 type=ParseMagickOption(MagickMetricOptions,MagickTrue,argv[i]);
720 if (type < 0)
721 ThrowCompareException(OptionError,"UnrecognizedMetricType",
722 argv[i]);
723 metric=(MetricType) type;
724 break;
725 }
726 if (LocaleCompare("monitor",option+1) == 0)
727 break;
728 ThrowCompareException(OptionError,"UnrecognizedOption",option)
729 }
730 case 'p':
731 {
732 if (LocaleCompare("passphrase",option+1) == 0)
733 {
734 if (*option == '+')
735 break;
736 i++;
737 if (i == (long) argc)
738 ThrowCompareException(OptionError,"MissingArgument",option);
739 break;
740 }
741 if (LocaleCompare("profile",option+1) == 0)
742 {
743 i++;
744 if (i == (long) (argc-1))
745 ThrowCompareException(OptionError,"MissingArgument",option);
746 break;
747 }
748 ThrowCompareException(OptionError,"UnrecognizedOption",option)
749 }
750 case 'q':
751 {
752 if (LocaleCompare("quality",option+1) == 0)
753 {
754 if (*option == '+')
755 break;
756 i++;
757 if (i == (long) (argc-1))
758 ThrowCompareException(OptionError,"MissingArgument",option);
759 if (IsGeometry(argv[i]) == MagickFalse)
760 ThrowCompareInvalidArgumentException(option,argv[i]);
761 break;
762 }
763 if (LocaleCompare("quantize",option+1) == 0)
764 {
765 long
766 colorspace;
767
768 if (*option == '+')
769 break;
770 i++;
771 if (i == (long) (argc-1))
772 ThrowCompareException(OptionError,"MissingArgument",option);
773 colorspace=ParseMagickOption(MagickColorspaceOptions,
774 MagickFalse,argv[i]);
775 if (colorspace < 0)
776 ThrowCompareException(OptionError,"UnrecognizedColorspace",
777 argv[i]);
778 break;
779 }
780 if (LocaleCompare("quiet",option+1) == 0)
781 break;
782 ThrowCompareException(OptionError,"UnrecognizedOption",option)
783 }
784 case 'r':
785 {
786 if (LocaleCompare("regard-warnings",option+1) == 0)
787 break;
788 if (LocaleNCompare("respect-parentheses",option+1,17) == 0)
789 {
790 respect_parenthesis=(*option == '-') ? MagickTrue : MagickFalse;
791 break;
792 }
793 ThrowCompareException(OptionError,"UnrecognizedOption",option)
794 }
795 case 's':
796 {
797 if (LocaleCompare("sampling-factor",option+1) == 0)
798 {
799 if (*option == '+')
800 break;
801 i++;
802 if (i == (long) argc)
803 ThrowCompareException(OptionError,"MissingArgument",option);
804 if (IsGeometry(argv[i]) == MagickFalse)
805 ThrowCompareInvalidArgumentException(option,argv[i]);
806 break;
807 }
808 if (LocaleCompare("seed",option+1) == 0)
809 {
810 if (*option == '+')
811 break;
812 i++;
813 if (i == (long) (argc-1))
814 ThrowCompareException(OptionError,"MissingArgument",option);
815 if (IsGeometry(argv[i]) == MagickFalse)
816 ThrowCompareInvalidArgumentException(option,argv[i]);
817 break;
818 }
819 if (LocaleCompare("set",option+1) == 0)
820 {
821 i++;
822 if (i == (long) argc)
823 ThrowCompareException(OptionError,"MissingArgument",option);
824 if (*option == '+')
825 break;
826 i++;
827 if (i == (long) argc)
828 ThrowCompareException(OptionError,"MissingArgument",option);
829 break;
830 }
831 if (LocaleCompare("size",option+1) == 0)
832 {
833 if (*option == '+')
834 break;
835 i++;
836 if (i == (long) argc)
837 ThrowCompareException(OptionError,"MissingArgument",option);
838 if (IsGeometry(argv[i]) == MagickFalse)
839 ThrowCompareInvalidArgumentException(option,argv[i]);
840 break;
841 }
842 ThrowCompareException(OptionError,"UnrecognizedOption",option)
843 }
844 case 't':
845 {
846 if (LocaleCompare("transparent-color",option+1) == 0)
847 {
848 if (*option == '+')
849 break;
850 i++;
851 if (i == (long) (argc-1))
852 ThrowCompareException(OptionError,"MissingArgument",option);
853 break;
854 }
855 if (LocaleCompare("type",option+1) == 0)
856 {
857 long
858 type;
859
860 if (*option == '+')
861 break;
862 i++;
863 if (i == (long) argc)
864 ThrowCompareException(OptionError,"MissingArgument",option);
865 type=ParseMagickOption(MagickTypeOptions,MagickFalse,argv[i]);
866 if (type < 0)
867 ThrowCompareException(OptionError,"UnrecognizedImageType",
868 argv[i]);
869 break;
870 }
871 ThrowCompareException(OptionError,"UnrecognizedOption",option)
872 }
873 case 'v':
874 {
875 if (LocaleCompare("verbose",option+1) == 0)
876 break;
877 if ((LocaleCompare("version",option+1) == 0) ||
878 (LocaleCompare("-version",option+1) == 0))
879 {
880 (void) fprintf(stdout,"Version: %s\n",
881 GetMagickVersion((unsigned long *) NULL));
cristy3980b0d2009-10-25 14:37:13 +0000[diff] [blame]882 (void) fprintf(stdout,"Copyright: %s\n",GetMagickCopyright());
cristy104cea82009-10-25 02:26:51 +0000[diff] [blame]883 (void) fprintf(stdout,"Features: %s\n\n",GetMagickFeatures());
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]884 break;
885 }
886 if (LocaleCompare("virtual-pixel",option+1) == 0)
887 {
888 long
889 method;
890
891 if (*option == '+')
892 break;
893 i++;
894 if (i == (long) (argc-1))
895 ThrowCompareException(OptionError,"MissingArgument",option);
896 method=ParseMagickOption(MagickVirtualPixelOptions,MagickFalse,
897 argv[i]);
898 if (method < 0)
899 ThrowCompareException(OptionError,
900 "UnrecognizedVirtualPixelMethod",argv[i]);
901 break;
902 }
903 ThrowCompareException(OptionError,"UnrecognizedOption",option)
904 }
905 case '?':
906 break;
907 default:
908 ThrowCompareException(OptionError,"UnrecognizedOption",option)
909 }
910 fire=ParseMagickOption(MagickImageListOptions,MagickFalse,option+1) < 0 ?
911 MagickFalse : MagickTrue;
912 if (fire != MagickFalse)
913 FireImageStack(MagickTrue,MagickTrue,MagickTrue);
914 }
915 if (k != 0)
916 ThrowCompareException(OptionError,"UnbalancedParenthesis",argv[i]);
917 if (i-- != (long) (argc-1))
918 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
919 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
920 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
921 FinalizeImageSettings(image_info,image,MagickTrue);
cristy94f20072009-09-12 02:12:36 +0000[diff] [blame]922 if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
923 ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]924 image=GetImageFromList(image,0);
925 reconstruct_image=GetImageFromList(image,1);
926 similarity_image=SimilarityImage(image,reconstruct_image,&offset,
927 &similarity_metric,exception);
928 if (similarity_metric > dissimilarity_threshold)
929 ThrowCompareException(ImageError,"ImagesTooDissimilar",image->filename);
930 if ((reconstruct_image->columns == image->columns) &&
931 (reconstruct_image->rows == image->rows))
932 difference_image=CompareImageChannels(image,reconstruct_image,channels,
933 metric,&distortion,exception);
934 else
935 if (similarity_image != (Image *) NULL)
936 {
937 Image
938 *composite_image;
939
940 /*
941 Determine if reconstructed image is a subimage of the image.
942 */
943 composite_image=CloneImage(image,0,0,MagickTrue,exception);
944 if (composite_image == (Image *) NULL)
945 difference_image=CompareImageChannels(image,reconstruct_image,
946 channels,metric,&distortion,exception);
947 else
948 {
949 (void) CompositeImage(composite_image,CopyCompositeOp,
950 reconstruct_image,offset.x,offset.y);
951 difference_image=CompareImageChannels(image,composite_image,
952 channels,metric,&distortion,exception);
953 if (difference_image != (Image *) NULL)
954 {
955 difference_image->page.x=offset.x;
956 difference_image->page.y=offset.y;
957 }
958 composite_image=DestroyImage(composite_image);
959 }
960 if (difference_image == (Image *) NULL)
961 similarity_image=DestroyImage(similarity_image);
962 else
963 {
964 AppendImageToList(&difference_image,similarity_image);
965 similarity_image=(Image *) NULL;
966 }
967 }
968 if (difference_image == (Image *) NULL)
969 status=0;
970 else
971 {
972 if (image_info->verbose != MagickFalse)
973 (void) IsImagesEqual(image,reconstruct_image);
974 if (*difference_image->magick == '\0')
975 (void) CopyMagickString(difference_image->magick,image->magick,
976 MaxTextExtent);
977 if (image_info->verbose == MagickFalse)
978 {
979 switch (metric)
980 {
981 case MeanAbsoluteErrorMetric:
982 case MeanSquaredErrorMetric:
983 case RootMeanSquaredErrorMetric:
984 case PeakAbsoluteErrorMetric:
985 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]986 (void) fprintf(stderr,"%.15g (%.15g)",QuantumRange*distortion,
987 (double) distortion);
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]988 if ((reconstruct_image->columns != image->columns) ||
989 (reconstruct_image->rows != image->rows))
990 (void) fprintf(stderr," @ %ld,%ld",difference_image->page.x,
991 difference_image->page.y);
992 (void) fprintf(stderr,"\n");
993 break;
994 }
995 case AbsoluteErrorMetric:
996 case PeakSignalToNoiseRatioMetric:
997 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]998 (void) fprintf(stderr,"%.15g",distortion);
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]999 if ((reconstruct_image->columns != image->columns) ||
1000 (reconstruct_image->rows != image->rows))
1001 (void) fprintf(stderr," @ %ld,%ld",difference_image->page.x,
1002 difference_image->page.y);
1003 (void) fprintf(stderr,"\n");
1004 break;
1005 }
1006 case MeanErrorPerPixelMetric:
1007 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1008 (void) fprintf(stderr,"%.15g (%.15g, %.15g)",distortion,
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1009 image->error.normalized_mean_error,
1010 image->error.normalized_maximum_error);
1011 if ((reconstruct_image->columns != image->columns) ||
1012 (reconstruct_image->rows != image->rows))
1013 (void) fprintf(stderr," @ %ld,%ld",difference_image->page.x,
1014 difference_image->page.y);
1015 (void) fprintf(stderr,"\n");
1016 break;
1017 }
1018 case UndefinedMetric:
1019 break;
1020 }
1021 }
1022 else
1023 {
1024 double
1025 *channel_distortion;
1026
1027 channel_distortion=GetImageChannelDistortions(image,reconstruct_image,
1028 metric,&image->exception);
1029 (void) fprintf(stderr,"Image: %s\n",image->filename);
1030 if ((reconstruct_image->columns != image->columns) ||
1031 (reconstruct_image->rows != image->rows))
1032 (void) fprintf(stderr,"Offset: %ld,%ld\n",difference_image->page.x,
1033 difference_image->page.y);
1034 (void) fprintf(stderr," Channel distortion: %s\n",
1035 MagickOptionToMnemonic(MagickMetricOptions,(long) metric));
1036 switch (metric)
1037 {
1038 case MeanAbsoluteErrorMetric:
1039 case MeanSquaredErrorMetric:
1040 case RootMeanSquaredErrorMetric:
1041 case PeakAbsoluteErrorMetric:
1042 {
1043 switch (image->colorspace)
1044 {
1045 case RGBColorspace:
1046 default:
1047 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1048 (void) fprintf(stderr," red: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1049 QuantumRange*channel_distortion[RedChannel],
1050 channel_distortion[RedChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1051 (void) fprintf(stderr," green: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1052 QuantumRange*channel_distortion[GreenChannel],
1053 channel_distortion[GreenChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1054 (void) fprintf(stderr," blue: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1055 QuantumRange*channel_distortion[BlueChannel],
1056 channel_distortion[BlueChannel]);
1057 if (image->matte != MagickFalse)
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1058 (void) fprintf(stderr," alpha: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1059 QuantumRange*channel_distortion[OpacityChannel],
1060 channel_distortion[OpacityChannel]);
1061 break;
1062 }
1063 case CMYKColorspace:
1064 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1065 (void) fprintf(stderr," cyan: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1066 QuantumRange*channel_distortion[CyanChannel],
1067 channel_distortion[CyanChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1068 (void) fprintf(stderr," magenta: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1069 QuantumRange*channel_distortion[MagentaChannel],
1070 channel_distortion[MagentaChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1071 (void) fprintf(stderr," yellow: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1072 QuantumRange*channel_distortion[YellowChannel],
1073 channel_distortion[YellowChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1074 (void) fprintf(stderr," black: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1075 QuantumRange*channel_distortion[BlackChannel],
1076 channel_distortion[BlackChannel]);
1077 if (image->matte != MagickFalse)
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1078 (void) fprintf(stderr," alpha: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1079 QuantumRange*channel_distortion[OpacityChannel],
1080 channel_distortion[OpacityChannel]);
1081 break;
1082 }
1083 case GRAYColorspace:
1084 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1085 (void) fprintf(stderr," gray: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1086 QuantumRange*channel_distortion[GrayChannel],
1087 channel_distortion[GrayChannel]);
1088 if (image->matte != MagickFalse)
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1089 (void) fprintf(stderr," alpha: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1090 QuantumRange*channel_distortion[OpacityChannel],
1091 channel_distortion[OpacityChannel]);
1092 break;
1093 }
1094 }
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1095 (void) fprintf(stderr," all: %.15g (%.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1096 QuantumRange*channel_distortion[AllChannels],
1097 channel_distortion[AllChannels]);
1098 break;
1099 }
1100 case AbsoluteErrorMetric:
1101 case PeakSignalToNoiseRatioMetric:
1102 {
1103 switch (image->colorspace)
1104 {
1105 case RGBColorspace:
1106 default:
1107 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1108 (void) fprintf(stderr," red: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1109 channel_distortion[RedChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1110 (void) fprintf(stderr," green: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1111 channel_distortion[GreenChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1112 (void) fprintf(stderr," blue: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1113 channel_distortion[BlueChannel]);
1114 if (image->matte != MagickFalse)
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1115 (void) fprintf(stderr," alpha: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1116 channel_distortion[OpacityChannel]);
1117 break;
1118 }
1119 case CMYKColorspace:
1120 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1121 (void) fprintf(stderr," cyan: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1122 channel_distortion[CyanChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1123 (void) fprintf(stderr," magenta: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1124 channel_distortion[MagentaChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1125 (void) fprintf(stderr," yellow: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1126 channel_distortion[YellowChannel]);
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1127 (void) fprintf(stderr," black: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1128 channel_distortion[BlackChannel]);
1129 if (image->matte != MagickFalse)
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1130 (void) fprintf(stderr," alpha: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1131 channel_distortion[OpacityChannel]);
1132 break;
1133 }
1134 case GRAYColorspace:
1135 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1136 (void) fprintf(stderr," gray: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1137 channel_distortion[GrayChannel]);
1138 if (image->matte != MagickFalse)
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1139 (void) fprintf(stderr," alpha: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1140 channel_distortion[OpacityChannel]);
1141 break;
1142 }
1143 }
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1144 (void) fprintf(stderr," all: %.15g\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1145 channel_distortion[AllChannels]);
1146 break;
1147 }
1148 case MeanErrorPerPixelMetric:
1149 {
cristy8cd5b312010-01-07 01:10:24 +0000[diff] [blame^]1150 (void) fprintf(stderr," %.15g (%.15g, %.15g)\n",
cristy3ed852e2009-09-05 21:47:34 +0000[diff] [blame]1151 channel_distortion[AllChannels],
1152 image->error.normalized_mean_error,
1153 image->error.normalized_maximum_error);
1154 break;
1155 }
1156 case UndefinedMetric:
1157 break;
1158 }
1159 channel_distortion=(double *) RelinquishMagickMemory(
1160 channel_distortion);
1161 }
1162 status&=WriteImages(image_info,difference_image,argv[argc-1],exception);
1163 if ((metadata != (char **) NULL) && (format != (char *) NULL))
1164 {
1165 char
1166 *text;
1167
1168 text=InterpretImageProperties(image_info,difference_image,format);
1169 if (text == (char *) NULL)
1170 ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
1171 GetExceptionMessage(errno));
1172 (void) ConcatenateString(&(*metadata),text);
1173 (void) ConcatenateString(&(*metadata),"\n");
1174 text=DestroyString(text);
1175 }
1176 difference_image=DestroyImageList(difference_image);
1177 }
1178 DestroyCompare();
1179 return(status != 0 ? MagickTrue : MagickFalse);
1180}