diff --git a/MagickWand/compare.c b/MagickWand/compare.c
new file mode 100644
index 0000000..ad328ca
--- /dev/null
+++ b/MagickWand/compare.c
@@ -0,0 +1,1210 @@
+/*
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% %
+% %
+% %
+% CCCC OOO M M PPPP AAA RRRR EEEEE %
+% C O O MM MM P P A A R R E %
+% C O O M M M PPPP AAAAA RRRR EEE %
+% C O O M M P A A R R E %
+% CCCC OOO M M P A A R R EEEEE %
+% %
+% %
+% Image Comparison Methods %
+% %
+% Software Design %
+% John Cristy %
+% December 2003 %
+% %
+% %
+% Copyright 1999-2011 ImageMagick Studio LLC, a non-profit organization %
+% dedicated to making software imaging solutions freely available. %
+% %
+% You may not use this file except in compliance with the License. You may %
+% obtain a copy of the License at %
+% %
+% http://www.imagemagick.org/script/license.php %
+% %
+% Unless required by applicable law or agreed to in writing, software %
+% distributed under the License is distributed on an "AS IS" BASIS, %
+% WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. %
+% See the License for the specific language governing permissions and %
+% limitations under the License. %
+% %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% Use the compare program to mathematically and visually annotate the
+% difference between an image and its reconstruction.
+%
+*/
+
+/*
+ Include declarations.
+*/
+#include "MagickWand/studio.h"
+#include "MagickWand/MagickWand.h"
+#include "MagickWand/mogrify-private.h"
+#include "MagickCore/string-private.h"
+
+/*
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% %
+% %
+% %
+% C o m p a r e I m a g e C o m m a n d %
+% %
+% %
+% %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% CompareImageCommand() compares two images and returns the difference between
+% them as a distortion metric and as a new image visually annotating their
+% differences.
+%
+% The format of the CompareImageCommand method is:
+%
+% MagickBooleanType CompareImageCommand(ImageInfo *image_info,int argc,
+% char **argv,char **metadata,ExceptionInfo *exception)
+%
+% A description of each parameter follows:
+%
+% o image_info: the image info.
+%
+% o argc: the number of elements in the argument vector.
+%
+% o argv: A text array containing the command line arguments.
+%
+% o metadata: any metadata is returned here.
+%
+% o exception: return any errors or warnings in this structure.
+%
+*/
+
+static MagickBooleanType CompareUsage(void)
+{
+ const char
+ **p;
+
+ static const char
+ *miscellaneous[]=
+ {
+ "-debug events display copious debugging information",
+ "-help print program options",
+ "-list type print a list of supported option arguments",
+ "-log format format of debugging information",
+ (char *) NULL
+ },
+ *settings[]=
+ {
+ "-alpha option on, activate, off, deactivate, set, opaque, copy",
+ " transparent, extract, background, or shape",
+ "-authenticate password",
+ " decipher image with this password",
+ "-channel type apply option to select image channels",
+ "-colorspace type alternate image colorspace",
+ "-compose operator set image composite operator",
+ "-compress type type of pixel compression when writing the image",
+ "-decipher filename convert cipher pixels to plain pixels",
+ "-define format:option",
+ " define one or more image format options",
+ "-density geometry horizontal and vertical density of the image",
+ "-depth value image depth",
+ "-dissimilarity-threshold value",
+ " maximum distortion for (sub)image match",
+ "-encipher filename convert plain pixels to cipher pixels",
+ "-extract geometry extract area from image",
+ "-format \"string\" output formatted image characteristics",
+ "-fuzz distance colors within this distance are considered equal",
+ "-highlight-color color",
+ " empasize pixel differences with this color",
+ "-identify identify the format and characteristics of the image",
+ "-interlace type type of image interlacing scheme",
+ "-limit type value pixel cache resource limit",
+ "-lowlight-color color",
+ " de-emphasize pixel differences with this color",
+ "-metric type measure differences between images with this metric",
+ "-monitor monitor progress",
+ "-passphrase filename get the passphrase from this file",
+ "-profile filename add, delete, or apply an image profile",
+ "-quality value JPEG/MIFF/PNG compression level",
+ "-quiet suppress all warning messages",
+ "-quantize colorspace reduce colors in this colorspace",
+ "-regard-warnings pay attention to warning messages",
+ "-respect-parentheses settings remain in effect until parenthesis boundary",
+ "-sampling-factor geometry",
+ " horizontal and vertical sampling factor",
+ "-seed value seed a new sequence of pseudo-random numbers",
+ "-set attribute value set an image attribute",
+ "-quality value JPEG/MIFF/PNG compression level",
+ "-size geometry width and height of image",
+ "-subimage-search search for subimage",
+ "-transparent-color color",
+ " transparent color",
+ "-type type image type",
+ "-verbose print detailed information about the image",
+ "-version print version information",
+ "-virtual-pixel method",
+ " virtual pixel access method",
+ (char *) NULL
+ };
+
+ (void) printf("Version: %s\n",GetMagickVersion((size_t *) NULL));
+ (void) printf("Copyright: %s\n",GetMagickCopyright());
+ (void) printf("Features: %s\n\n",GetMagickFeatures());
+ (void) printf("Usage: %s [options ...] image reconstruct difference\n",
+ GetClientName());
+ (void) printf("\nImage Settings:\n");
+ for (p=settings; *p != (char *) NULL; p++)
+ (void) printf(" %s\n",*p);
+ (void) printf("\nMiscellaneous Options:\n");
+ for (p=miscellaneous; *p != (char *) NULL; p++)
+ (void) printf(" %s\n",*p);
+ (void) printf(
+ "\nBy default, the image format of `file' is determined by its magic\n");
+ (void) printf(
+ "number. To specify a particular image format, precede the filename\n");
+ (void) printf(
+ "with an image format name and a colon (i.e. ps:image) or specify the\n");
+ (void) printf(
+ "image type as the filename suffix (i.e. image.ps). Specify 'file' as\n");
+ (void) printf("'-' for standard input or output.\n");
+ return(MagickFalse);
+}
+
+WandExport MagickBooleanType CompareImageCommand(ImageInfo *image_info,
+ int argc,char **argv,char **metadata,ExceptionInfo *exception)
+{
+#define DefaultDissimilarityThreshold 0.31830988618379067154
+#define DestroyCompare() \
+{ \
+ if (similarity_image != (Image *) NULL) \
+ similarity_image=DestroyImageList(similarity_image); \
+ if (difference_image != (Image *) NULL) \
+ difference_image=DestroyImageList(difference_image); \
+ DestroyImageStack(); \
+ for (i=0; i < (ssize_t) argc; i++) \
+ argv[i]=DestroyString(argv[i]); \
+ argv=(char **) RelinquishMagickMemory(argv); \
+}
+#define ThrowCompareException(asperity,tag,option) \
+{ \
+ if (exception->severity < (asperity)) \
+ (void) ThrowMagickException(exception,GetMagickModule(),asperity,tag, \
+ "`%s'",option); \
+ DestroyCompare(); \
+ return(MagickFalse); \
+}
+#define ThrowCompareInvalidArgumentException(option,argument) \
+{ \
+ (void) ThrowMagickException(exception,GetMagickModule(),OptionError, \
+ "InvalidArgument","`%s': %s",option,argument); \
+ DestroyCompare(); \
+ return(MagickFalse); \
+}
+
+ char
+ *filename,
+ *option;
+
+ const char
+ *format;
+
+ ChannelType
+ channels;
+
+ double
+ dissimilarity_threshold,
+ distortion,
+ similarity_metric;
+
+ Image
+ *difference_image,
+ *image,
+ *reconstruct_image,
+ *similarity_image;
+
+ ImageStack
+ image_stack[MaxImageStackDepth+1];
+
+ MagickBooleanType
+ fire,
+ pend,
+ respect_parenthesis,
+ subimage_search;
+
+ MagickStatusType
+ status;
+
+ MetricType
+ metric;
+
+ RectangleInfo
+ offset;
+
+ register ssize_t
+ i;
+
+ ssize_t
+ j,
+ k;
+
+ /*
+ Set defaults.
+ */
+ assert(image_info != (ImageInfo *) NULL);
+ assert(image_info->signature == MagickSignature);
+ if (image_info->debug != MagickFalse)
+ (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
+ assert(exception != (ExceptionInfo *) NULL);
+ if (argc == 2)
+ {
+ option=argv[1];
+ if ((LocaleCompare("version",option+1) == 0) ||
+ (LocaleCompare("-version",option+1) == 0))
+ {
+ (void) FormatLocaleFile(stdout,"Version: %s\n",
+ GetMagickVersion((size_t *) NULL));
+ (void) FormatLocaleFile(stdout,"Copyright: %s\n",
+ GetMagickCopyright());
+ (void) FormatLocaleFile(stdout,"Features: %s\n\n",
+ GetMagickFeatures());
+ return(MagickFalse);
+ }
+ }
+ if (argc < 3)
+ return(CompareUsage());
+ channels=CompositeChannels;
+ difference_image=NewImageList();
+ similarity_image=NewImageList();
+ dissimilarity_threshold=DefaultDissimilarityThreshold;
+ distortion=0.0;
+ format=(char *) NULL;
+ j=1;
+ k=0;
+ metric=UndefinedMetric;
+ NewImageStack();
+ option=(char *) NULL;
+ pend=MagickFalse;
+ reconstruct_image=NewImageList();
+ respect_parenthesis=MagickFalse;
+ status=MagickTrue;
+ subimage_search=MagickFalse;
+ /*
+ Compare an image.
+ */
+ ReadCommandlLine(argc,&argv);
+ status=ExpandFilenames(&argc,&argv);
+ if (status == MagickFalse)
+ ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
+ GetExceptionMessage(errno));
+ for (i=1; i < (ssize_t) (argc-1); i++)
+ {
+ option=argv[i];
+ if (LocaleCompare(option,"(") == 0)
+ {
+ FireImageStack(MagickTrue,MagickTrue,pend);
+ if (k == MaxImageStackDepth)
+ ThrowCompareException(OptionError,"ParenthesisNestedTooDeeply",
+ option);
+ PushImageStack();
+ continue;
+ }
+ if (LocaleCompare(option,")") == 0)
+ {
+ FireImageStack(MagickTrue,MagickTrue,MagickTrue);
+ if (k == 0)
+ ThrowCompareException(OptionError,"UnableToParseExpression",option);
+ PopImageStack();
+ continue;
+ }
+ if (IsCommandOption(option) == MagickFalse)
+ {
+ Image
+ *images;
+
+ /*
+ Read input image.
+ */
+ FireImageStack(MagickFalse,MagickFalse,pend);
+ filename=argv[i];
+ if ((LocaleCompare(filename,"--") == 0) && (i < (ssize_t) (argc-1)))
+ filename=argv[++i];
+ (void) CopyMagickString(image_info->filename,filename,MaxTextExtent);
+ images=ReadImages(image_info,exception);
+ status&=(images != (Image *) NULL) &&
+ (exception->severity < ErrorException);
+ if (images == (Image *) NULL)
+ continue;
+ AppendImageStack(images);
+ continue;
+ }
+ pend=image != (Image *) NULL ? MagickTrue : MagickFalse;
+ switch (*(option+1))
+ {
+ case 'a':
+ {
+ if (LocaleCompare("alpha",option+1) == 0)
+ {
+ ssize_t
+ type;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ type=ParseCommandOption(MagickAlphaOptions,MagickFalse,argv[i]);
+ if (type < 0)
+ ThrowCompareException(OptionError,"UnrecognizedAlphaChannelType",
+ argv[i]);
+ break;
+ }
+ if (LocaleCompare("authenticate",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option);
+ }
+ case 'c':
+ {
+ if (LocaleCompare("cache",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("channel",option+1) == 0)
+ {
+ ssize_t
+ channel;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ channel=ParseChannelOption(argv[i]);
+ if (channel < 0)
+ ThrowCompareException(OptionError,"UnrecognizedChannelType",
+ argv[i]);
+ channels=(ChannelType) channel;
+ break;
+ }
+ if (LocaleCompare("colorspace",option+1) == 0)
+ {
+ ssize_t
+ colorspace;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ colorspace=ParseCommandOption(MagickColorspaceOptions,MagickFalse,
+ argv[i]);
+ if (colorspace < 0)
+ ThrowCompareException(OptionError,"UnrecognizedColorspace",
+ argv[i]);
+ break;
+ }
+ if (LocaleCompare("compose",option+1) == 0)
+ {
+ ssize_t
+ compose;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ compose=ParseCommandOption(MagickComposeOptions,MagickFalse,
+ argv[i]);
+ if (compose < 0)
+ ThrowCompareException(OptionError,"UnrecognizedComposeOperator",
+ argv[i]);
+ break;
+ }
+ if (LocaleCompare("compress",option+1) == 0)
+ {
+ ssize_t
+ compress;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ compress=ParseCommandOption(MagickCompressOptions,MagickFalse,
+ argv[i]);
+ if (compress < 0)
+ ThrowCompareException(OptionError,"UnrecognizedImageCompression",
+ argv[i]);
+ break;
+ }
+ if (LocaleCompare("concurrent",option+1) == 0)
+ break;
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'd':
+ {
+ if (LocaleCompare("debug",option+1) == 0)
+ {
+ LogEventType
+ event_mask;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ event_mask=SetLogEventMask(argv[i]);
+ if (event_mask == UndefinedEvents)
+ ThrowCompareException(OptionError,"UnrecognizedEventType",
+ argv[i]);
+ break;
+ }
+ if (LocaleCompare("decipher",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ if (LocaleCompare("define",option+1) == 0)
+ {
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (*option == '+')
+ {
+ const char
+ *define;
+
+ define=GetImageOption(image_info,argv[i]);
+ if (define == (const char *) NULL)
+ ThrowCompareException(OptionError,"NoSuchOption",argv[i]);
+ break;
+ }
+ break;
+ }
+ if (LocaleCompare("density",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("depth",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("dissimilarity-threshold",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ if (*option == '+')
+ dissimilarity_threshold=DefaultDissimilarityThreshold;
+ else
+ dissimilarity_threshold=InterpretLocaleValue(argv[i],
+ (char **) NULL);
+ break;
+ }
+ if (LocaleCompare("duration",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'e':
+ {
+ if (LocaleCompare("encipher",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ if (LocaleCompare("extract",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'f':
+ {
+ if (LocaleCompare("format",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ format=argv[i];
+ break;
+ }
+ if (LocaleCompare("fuzz",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'h':
+ {
+ if ((LocaleCompare("help",option+1) == 0) ||
+ (LocaleCompare("-help",option+1) == 0))
+ return(CompareUsage());
+ if (LocaleCompare("highlight-color",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'i':
+ {
+ if (LocaleCompare("identify",option+1) == 0)
+ break;
+ if (LocaleCompare("interlace",option+1) == 0)
+ {
+ ssize_t
+ interlace;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ interlace=ParseCommandOption(MagickInterlaceOptions,MagickFalse,
+ argv[i]);
+ if (interlace < 0)
+ ThrowCompareException(OptionError,"UnrecognizedInterlaceType",
+ argv[i]);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'l':
+ {
+ if (LocaleCompare("limit",option+1) == 0)
+ {
+ char
+ *p;
+
+ double
+ value;
+
+ ssize_t
+ resource;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ resource=ParseCommandOption(MagickResourceOptions,MagickFalse,
+ argv[i]);
+ if (resource < 0)
+ ThrowCompareException(OptionError,"UnrecognizedResourceType",
+ argv[i]);
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ value=InterpretLocaleValue(argv[i],&p);
+ (void) value;
+ if ((p == argv[i]) && (LocaleCompare("unlimited",argv[i]) != 0))
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("list",option+1) == 0)
+ {
+ ssize_t
+ list;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ list=ParseCommandOption(MagickListOptions,MagickFalse,argv[i]);
+ if (list < 0)
+ ThrowCompareException(OptionError,"UnrecognizedListType",argv[i]);
+ status=MogrifyImageInfo(image_info,(int) (i-j+1),(const char **)
+ argv+j,exception);
+ DestroyCompare();
+ return(status != 0 ? MagickFalse : MagickTrue);
+ }
+ if (LocaleCompare("log",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if ((i == (ssize_t) argc) || (strchr(argv[i],'%') == (char *) NULL))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ if (LocaleCompare("lowlight-color",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'm':
+ {
+ if (LocaleCompare("matte",option+1) == 0)
+ break;
+ if (LocaleCompare("metric",option+1) == 0)
+ {
+ ssize_t
+ type;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ type=ParseCommandOption(MagickMetricOptions,MagickTrue,argv[i]);
+ if (type < 0)
+ ThrowCompareException(OptionError,"UnrecognizedMetricType",
+ argv[i]);
+ metric=(MetricType) type;
+ break;
+ }
+ if (LocaleCompare("monitor",option+1) == 0)
+ break;
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'p':
+ {
+ if (LocaleCompare("passphrase",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ if (LocaleCompare("profile",option+1) == 0)
+ {
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'q':
+ {
+ if (LocaleCompare("quality",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("quantize",option+1) == 0)
+ {
+ ssize_t
+ colorspace;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ colorspace=ParseCommandOption(MagickColorspaceOptions,
+ MagickFalse,argv[i]);
+ if (colorspace < 0)
+ ThrowCompareException(OptionError,"UnrecognizedColorspace",
+ argv[i]);
+ break;
+ }
+ if (LocaleCompare("quiet",option+1) == 0)
+ break;
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'r':
+ {
+ if (LocaleCompare("regard-warnings",option+1) == 0)
+ break;
+ if (LocaleNCompare("respect-parentheses",option+1,17) == 0)
+ {
+ respect_parenthesis=(*option == '-') ? MagickTrue : MagickFalse;
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 's':
+ {
+ if (LocaleCompare("sampling-factor",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("seed",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("set",option+1) == 0)
+ {
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ if (LocaleCompare("size",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ if (IsGeometry(argv[i]) == MagickFalse)
+ ThrowCompareInvalidArgumentException(option,argv[i]);
+ break;
+ }
+ if (LocaleCompare("subimage-search",option+1) == 0)
+ {
+ if (*option == '+')
+ {
+ subimage_search=MagickFalse;
+ break;
+ }
+ subimage_search=MagickTrue;
+ break;
+ }
+ if (LocaleCompare("synchronize",option+1) == 0)
+ break;
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 't':
+ {
+ if (LocaleCompare("taint",option+1) == 0)
+ break;
+ if (LocaleCompare("transparent-color",option+1) == 0)
+ {
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ break;
+ }
+ if (LocaleCompare("type",option+1) == 0)
+ {
+ ssize_t
+ type;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) argc)
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ type=ParseCommandOption(MagickTypeOptions,MagickFalse,argv[i]);
+ if (type < 0)
+ ThrowCompareException(OptionError,"UnrecognizedImageType",
+ argv[i]);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case 'v':
+ {
+ if (LocaleCompare("verbose",option+1) == 0)
+ break;
+ if ((LocaleCompare("version",option+1) == 0) ||
+ (LocaleCompare("-version",option+1) == 0))
+ {
+ (void) FormatLocaleFile(stdout,"Version: %s\n",
+ GetMagickVersion((size_t *) NULL));
+ (void) FormatLocaleFile(stdout,"Copyright: %s\n",
+ GetMagickCopyright());
+ (void) FormatLocaleFile(stdout,"Features: %s\n\n",
+ GetMagickFeatures());
+ break;
+ }
+ if (LocaleCompare("virtual-pixel",option+1) == 0)
+ {
+ ssize_t
+ method;
+
+ if (*option == '+')
+ break;
+ i++;
+ if (i == (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingArgument",option);
+ method=ParseCommandOption(MagickVirtualPixelOptions,MagickFalse,
+ argv[i]);
+ if (method < 0)
+ ThrowCompareException(OptionError,
+ "UnrecognizedVirtualPixelMethod",argv[i]);
+ break;
+ }
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ case '?':
+ break;
+ default:
+ ThrowCompareException(OptionError,"UnrecognizedOption",option)
+ }
+ fire=(GetCommandOptionFlags(MagickCommandOptions,MagickFalse,option) &
+ FireOptionFlag) == 0 ? MagickFalse : MagickTrue;
+ if (fire != MagickFalse)
+ FireImageStack(MagickTrue,MagickTrue,MagickTrue);
+ }
+ if (k != 0)
+ ThrowCompareException(OptionError,"UnbalancedParenthesis",argv[i]);
+ if (i-- != (ssize_t) (argc-1))
+ ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
+ if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
+ ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
+ FinalizeImageSettings(image_info,image,MagickTrue);
+ if ((image == (Image *) NULL) || (GetImageListLength(image) < 2))
+ ThrowCompareException(OptionError,"MissingAnImageFilename",argv[i]);
+ image=GetImageFromList(image,0);
+ reconstruct_image=GetImageFromList(image,1);
+ if (subimage_search != MagickFalse)
+ {
+ similarity_image=SimilarityImage(image,reconstruct_image,&offset,
+ &similarity_metric,exception);
+ if (similarity_metric > dissimilarity_threshold)
+ ThrowCompareException(ImageError,"ImagesTooDissimilar",image->filename);
+ }
+ if ((reconstruct_image->columns == image->columns) &&
+ (reconstruct_image->rows == image->rows))
+ difference_image=CompareImageChannels(image,reconstruct_image,channels,
+ metric,&distortion,exception);
+ else
+ if (similarity_image == (Image *) NULL)
+ ThrowCompareException(OptionError,"ImageWidthsOrHeightsDiffer",
+ image->filename)
+ else
+ {
+ Image
+ *composite_image;
+
+ /*
+ Determine if reconstructed image is a subimage of the image.
+ */
+ composite_image=CloneImage(image,0,0,MagickTrue,exception);
+ if (composite_image == (Image *) NULL)
+ difference_image=CompareImageChannels(image,reconstruct_image,
+ channels,metric,&distortion,exception);
+ else
+ {
+ (void) CompositeImage(composite_image,CopyCompositeOp,
+ reconstruct_image,offset.x,offset.y);
+ difference_image=CompareImageChannels(image,composite_image,
+ channels,metric,&distortion,exception);
+ if (difference_image != (Image *) NULL)
+ {
+ difference_image->page.x=offset.x;
+ difference_image->page.y=offset.y;
+ }
+ composite_image=DestroyImage(composite_image);
+ }
+ if (difference_image != (Image *) NULL)
+ {
+ AppendImageToList(&difference_image,similarity_image);
+ similarity_image=(Image *) NULL;
+ }
+ }
+ if (difference_image == (Image *) NULL)
+ status=0;
+ else
+ {
+ if (image_info->verbose != MagickFalse)
+ (void) IsImagesEqual(image,reconstruct_image);
+ if (*difference_image->magick == '\0')
+ (void) CopyMagickString(difference_image->magick,image->magick,
+ MaxTextExtent);
+ if (image_info->verbose == MagickFalse)
+ {
+ switch (metric)
+ {
+ case FuzzErrorMetric:
+ case MeanAbsoluteErrorMetric:
+ case MeanSquaredErrorMetric:
+ case RootMeanSquaredErrorMetric:
+ case PeakAbsoluteErrorMetric:
+ {
+ (void) FormatLocaleFile(stderr,"%g (%g)",QuantumRange*distortion,
+ (double) distortion);
+ if ((reconstruct_image->columns != image->columns) ||
+ (reconstruct_image->rows != image->rows))
+ (void) FormatLocaleFile(stderr," @ %.20g,%.20g",(double)
+ difference_image->page.x,(double) difference_image->page.y);
+ (void) FormatLocaleFile(stderr,"\n");
+ break;
+ }
+ case AbsoluteErrorMetric:
+ case NormalizedCrossCorrelationErrorMetric:
+ case PeakSignalToNoiseRatioMetric:
+ {
+ (void) FormatLocaleFile(stderr,"%g",distortion);
+ if ((reconstruct_image->columns != image->columns) ||
+ (reconstruct_image->rows != image->rows))
+ (void) FormatLocaleFile(stderr," @ %.20g,%.20g",(double)
+ difference_image->page.x,(double) difference_image->page.y);
+ (void) FormatLocaleFile(stderr,"\n");
+ break;
+ }
+ case MeanErrorPerPixelMetric:
+ {
+ (void) FormatLocaleFile(stderr,"%g (%g, %g)",distortion,
+ image->error.normalized_mean_error,
+ image->error.normalized_maximum_error);
+ if ((reconstruct_image->columns != image->columns) ||
+ (reconstruct_image->rows != image->rows))
+ (void) FormatLocaleFile(stderr," @ %.20g,%.20g",(double)
+ difference_image->page.x,(double) difference_image->page.y);
+ (void) FormatLocaleFile(stderr,"\n");
+ break;
+ }
+ case UndefinedMetric:
+ break;
+ }
+ }
+ else
+ {
+ double
+ *channel_distortion;
+
+ channel_distortion=GetImageChannelDistortions(image,reconstruct_image,
+ metric,&image->exception);
+ (void) FormatLocaleFile(stderr,"Image: %s\n",image->filename);
+ if ((reconstruct_image->columns != image->columns) ||
+ (reconstruct_image->rows != image->rows))
+ (void) FormatLocaleFile(stderr,"Offset: %.20g,%.20g\n",(double)
+ difference_image->page.x,(double) difference_image->page.y);
+ (void) FormatLocaleFile(stderr," Channel distortion: %s\n",
+ CommandOptionToMnemonic(MagickMetricOptions,(ssize_t) metric));
+ switch (metric)
+ {
+ case FuzzErrorMetric:
+ case MeanAbsoluteErrorMetric:
+ case MeanSquaredErrorMetric:
+ case RootMeanSquaredErrorMetric:
+ case PeakAbsoluteErrorMetric:
+ {
+ switch (image->colorspace)
+ {
+ case RGBColorspace:
+ default:
+ {
+ (void) FormatLocaleFile(stderr," red: %g (%g)\n",
+ QuantumRange*channel_distortion[RedChannel],
+ channel_distortion[RedChannel]);
+ (void) FormatLocaleFile(stderr," green: %g (%g)\n",
+ QuantumRange*channel_distortion[GreenChannel],
+ channel_distortion[GreenChannel]);
+ (void) FormatLocaleFile(stderr," blue: %g (%g)\n",
+ QuantumRange*channel_distortion[BlueChannel],
+ channel_distortion[BlueChannel]);
+ if (image->matte != MagickFalse)
+ (void) FormatLocaleFile(stderr," alpha: %g (%g)\n",
+ QuantumRange*channel_distortion[OpacityChannel],
+ channel_distortion[OpacityChannel]);
+ break;
+ }
+ case CMYKColorspace:
+ {
+ (void) FormatLocaleFile(stderr," cyan: %g (%g)\n",
+ QuantumRange*channel_distortion[CyanChannel],
+ channel_distortion[CyanChannel]);
+ (void) FormatLocaleFile(stderr," magenta: %g (%g)\n",
+ QuantumRange*channel_distortion[MagentaChannel],
+ channel_distortion[MagentaChannel]);
+ (void) FormatLocaleFile(stderr," yellow: %g (%g)\n",
+ QuantumRange*channel_distortion[YellowChannel],
+ channel_distortion[YellowChannel]);
+ (void) FormatLocaleFile(stderr," black: %g (%g)\n",
+ QuantumRange*channel_distortion[BlackChannel],
+ channel_distortion[BlackChannel]);
+ if (image->matte != MagickFalse)
+ (void) FormatLocaleFile(stderr," alpha: %g (%g)\n",
+ QuantumRange*channel_distortion[OpacityChannel],
+ channel_distortion[OpacityChannel]);
+ break;
+ }
+ case GRAYColorspace:
+ {
+ (void) FormatLocaleFile(stderr," gray: %g (%g)\n",
+ QuantumRange*channel_distortion[GrayChannel],
+ channel_distortion[GrayChannel]);
+ if (image->matte != MagickFalse)
+ (void) FormatLocaleFile(stderr," alpha: %g (%g)\n",
+ QuantumRange*channel_distortion[OpacityChannel],
+ channel_distortion[OpacityChannel]);
+ break;
+ }
+ }
+ (void) FormatLocaleFile(stderr," all: %g (%g)\n",
+ QuantumRange*channel_distortion[CompositeChannels],
+ channel_distortion[CompositeChannels]);
+ break;
+ }
+ case AbsoluteErrorMetric:
+ case NormalizedCrossCorrelationErrorMetric:
+ case PeakSignalToNoiseRatioMetric:
+ {
+ switch (image->colorspace)
+ {
+ case RGBColorspace:
+ default:
+ {
+ (void) FormatLocaleFile(stderr," red: %g\n",
+ channel_distortion[RedChannel]);
+ (void) FormatLocaleFile(stderr," green: %g\n",
+ channel_distortion[GreenChannel]);
+ (void) FormatLocaleFile(stderr," blue: %g\n",
+ channel_distortion[BlueChannel]);
+ if (image->matte != MagickFalse)
+ (void) FormatLocaleFile(stderr," alpha: %g\n",
+ channel_distortion[OpacityChannel]);
+ break;
+ }
+ case CMYKColorspace:
+ {
+ (void) FormatLocaleFile(stderr," cyan: %g\n",
+ channel_distortion[CyanChannel]);
+ (void) FormatLocaleFile(stderr," magenta: %g\n",
+ channel_distortion[MagentaChannel]);
+ (void) FormatLocaleFile(stderr," yellow: %g\n",
+ channel_distortion[YellowChannel]);
+ (void) FormatLocaleFile(stderr," black: %g\n",
+ channel_distortion[BlackChannel]);
+ if (image->matte != MagickFalse)
+ (void) FormatLocaleFile(stderr," alpha: %g\n",
+ channel_distortion[OpacityChannel]);
+ break;
+ }
+ case GRAYColorspace:
+ {
+ (void) FormatLocaleFile(stderr," gray: %g\n",
+ channel_distortion[GrayChannel]);
+ if (image->matte != MagickFalse)
+ (void) FormatLocaleFile(stderr," alpha: %g\n",
+ channel_distortion[OpacityChannel]);
+ break;
+ }
+ }
+ (void) FormatLocaleFile(stderr," all: %g\n",
+ channel_distortion[CompositeChannels]);
+ break;
+ }
+ case MeanErrorPerPixelMetric:
+ {
+ (void) FormatLocaleFile(stderr," %g (%g, %g)\n",
+ channel_distortion[CompositeChannels],
+ image->error.normalized_mean_error,
+ image->error.normalized_maximum_error);
+ break;
+ }
+ case UndefinedMetric:
+ break;
+ }
+ channel_distortion=(double *) RelinquishMagickMemory(
+ channel_distortion);
+ }
+ status&=WriteImages(image_info,difference_image,argv[argc-1],exception);
+ if ((metadata != (char **) NULL) && (format != (char *) NULL))
+ {
+ char
+ *text;
+
+ text=InterpretImageProperties(image_info,difference_image,format);
+ if (text == (char *) NULL)
+ ThrowCompareException(ResourceLimitError,"MemoryAllocationFailed",
+ GetExceptionMessage(errno));
+ (void) ConcatenateString(&(*metadata),text);
+ (void) ConcatenateString(&(*metadata),"\n");
+ text=DestroyString(text);
+ }
+ difference_image=DestroyImageList(difference_image);
+ }
+ DestroyCompare();
+ return(status != 0 ? MagickTrue : MagickFalse);
+}