granicus.if.org Git - imagemagick/blob - MagickCore/draw.c

   1 /*
   2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   3 %                                                                             %
   4 %                                                                             %
   5 %                                                                             %
   6 %                        DDDD   RRRR    AAA   W   W                           %
   7 %                        D   D  R   R  A   A  W   W                           %
   8 %                        D   D  RRRR   AAAAA  W W W                           %
   9 %                        D   D  R RN   A   A  WW WW                           %
  10 %                        DDDD   R  R   A   A  W   W                           %
  11 %                                                                             %
  12 %                                                                             %
  13 %                     MagickCore Image Drawing Methods                        %
  14 %                                                                             %
  15 %                                                                             %
  16 %                              Software Design                                %
  17 %                                John Cristy                                  %
  18 %                                 July 1998                                   %
  19 %                                                                             %
  20 %                                                                             %
  21 %  Copyright 1999-2012 ImageMagick Studio LLC, a non-profit organization      %
  22 %  dedicated to making software imaging solutions freely available.           %
  23 %                                                                             %
  24 %  You may not use this file except in compliance with the License.  You may  %
  25 %  obtain a copy of the License at                                            %
  26 %                                                                             %
  27 %    http://www.imagemagick.org/script/license.php                            %
  28 %                                                                             %
  29 %  Unless required by applicable law or agreed to in writing, software        %
  30 %  distributed under the License is distributed on an "AS IS" BASIS,          %
  31 %  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.   %
  32 %  See the License for the specific language governing permissions and        %
  33 %  limitations under the License.                                             %
  34 %                                                                             %
  35 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  36 %
  37 % Bill Radcliffe of Corbis (www.corbis.com) contributed the polygon
  38 % rendering code based on Paul Heckbert's "Concave Polygon Scan Conversion",
  39 % Graphics Gems, 1990.  Leonard Rosenthal and David Harr of Appligent
  40 % (www.appligent.com) contributed the dash pattern, linecap stroking
  41 % algorithm, and minor rendering improvements.
  42 %
  43 */
  44 \f
  45 /*
  46   Include declarations.
  47 */
  48 #include "MagickCore/studio.h"
  49 #include "MagickCore/annotate.h"
  50 #include "MagickCore/artifact.h"
  51 #include "MagickCore/blob.h"
  52 #include "MagickCore/cache.h"
  53 #include "MagickCore/cache-view.h"
  54 #include "MagickCore/color.h"
  55 #include "MagickCore/composite.h"
  56 #include "MagickCore/composite-private.h"
  57 #include "MagickCore/constitute.h"
  58 #include "MagickCore/draw.h"
  59 #include "MagickCore/draw-private.h"
  60 #include "MagickCore/enhance.h"
  61 #include "MagickCore/exception.h"
  62 #include "MagickCore/exception-private.h"
  63 #include "MagickCore/gem.h"
  64 #include "MagickCore/geometry.h"
  65 #include "MagickCore/image-private.h"
  66 #include "MagickCore/list.h"
  67 #include "MagickCore/log.h"
  68 #include "MagickCore/monitor.h"
  69 #include "MagickCore/monitor-private.h"
  70 #include "MagickCore/option.h"
  71 #include "MagickCore/paint.h"
  72 #include "MagickCore/pixel-accessor.h"
  73 #include "MagickCore/property.h"
  74 #include "MagickCore/resample.h"
  75 #include "MagickCore/resample-private.h"
  76 #include "MagickCore/string_.h"
  77 #include "MagickCore/string-private.h"
  78 #include "MagickCore/thread-private.h"
  79 #include "MagickCore/token.h"
  80 #include "MagickCore/transform.h"
  81 #include "MagickCore/utility.h"
  82 \f
  83 /*
  84   Define declarations.
  85 */
  86 #define BezierQuantum  200
  87 \f
  88 /*
  89   Typedef declarations.
  90 */
  91 typedef struct _EdgeInfo
  92 {
  93   SegmentInfo
  94     bounds;
  95
  96   MagickRealType
  97     scanline;
  98
  99   PointInfo
 100     *points;
 101
 102   size_t
 103     number_points;
 104
 105   ssize_t
 106     direction;
 107
 108   MagickBooleanType
 109     ghostline;
 110
 111   size_t
 112     highwater;
 113 } EdgeInfo;
 114
 115 typedef struct _ElementInfo
 116 {
 117   MagickRealType
 118     cx,
 119     cy,
 120     major,
 121     minor,
 122     angle;
 123 } ElementInfo;
 124
 125 typedef struct _PolygonInfo
 126 {
 127   EdgeInfo
 128     *edges;
 129
 130   size_t
 131     number_edges;
 132 } PolygonInfo;
 133
 134 typedef enum
 135 {
 136   MoveToCode,
 137   OpenCode,
 138   GhostlineCode,
 139   LineToCode,
 140   EndCode
 141 } PathInfoCode;
 142
 143 typedef struct _PathInfo
 144 {
 145   PointInfo
 146     point;
 147
 148   PathInfoCode
 149     code;
 150 } PathInfo;
 151 \f
 152 /*
 153   Forward declarations.
 154 */
 155 static MagickBooleanType
 156   DrawStrokePolygon(Image *,const DrawInfo *,const PrimitiveInfo *,
 157     ExceptionInfo *);
 158
 159 static PrimitiveInfo
 160   *TraceStrokePolygon(const DrawInfo *,const PrimitiveInfo *);
 161
 162 static size_t
 163   TracePath(PrimitiveInfo *,const char *);
 164
 165 static void
 166   TraceArc(PrimitiveInfo *,const PointInfo,const PointInfo,const PointInfo),
 167   TraceArcPath(PrimitiveInfo *,const PointInfo,const PointInfo,const PointInfo,
 168     const MagickRealType,const MagickBooleanType,const MagickBooleanType),
 169   TraceBezier(PrimitiveInfo *,const size_t),
 170   TraceCircle(PrimitiveInfo *,const PointInfo,const PointInfo),
 171   TraceEllipse(PrimitiveInfo *,const PointInfo,const PointInfo,const PointInfo),
 172   TraceLine(PrimitiveInfo *,const PointInfo,const PointInfo),
 173   TraceRectangle(PrimitiveInfo *,const PointInfo,const PointInfo),
 174   TraceRoundRectangle(PrimitiveInfo *,const PointInfo,const PointInfo,
 175     PointInfo),
 176   TraceSquareLinecap(PrimitiveInfo *,const size_t,const MagickRealType);
 177 \f
 178 /*
 179 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 180 %                                                                             %
 181 %                                                                             %
 182 %                                                                             %
 183 %   A c q u i r e D r a w I n f o                                             %
 184 %                                                                             %
 185 %                                                                             %
 186 %                                                                             %
 187 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 188 %
 189 %  AcquireDrawInfo() returns a DrawInfo structure properly initialized.
 190 %
 191 %  The format of the AcquireDrawInfo method is:
 192 %
 193 %      DrawInfo *AcquireDrawInfo(void)
 194 %
 195 */
 196 MagickExport DrawInfo *AcquireDrawInfo(void)
 197 {
 198   DrawInfo
 199     *draw_info;
 200
 201   draw_info=(DrawInfo *) AcquireMagickMemory(sizeof(*draw_info));
 202   if (draw_info == (DrawInfo *) NULL)
 203     ThrowFatalException(ResourceLimitFatalError,"MemoryAllocationFailed");
 204   GetDrawInfo((ImageInfo *) NULL,draw_info);
 205   return(draw_info);
 206 }
 207 \f
 208 /*
 209 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 210 %                                                                             %
 211 %                                                                             %
 212 %                                                                             %
 213 %   C l o n e D r a w I n f o                                                 %
 214 %                                                                             %
 215 %                                                                             %
 216 %                                                                             %
 217 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 218 %
 219 %  CloneDrawInfo() makes a copy of the given draw_info structure.  If NULL
 220 %  is specified, a new draw_info structure is created initialized to
 221 %  default values, according to the given image_info.
 222 %
 223 %  The format of the CloneDrawInfo method is:
 224 %
 225 %      DrawInfo *CloneDrawInfo(const ImageInfo *image_info,
 226 %        const DrawInfo *draw_info)
 227 %
 228 %  A description of each parameter follows:
 229 %
 230 %    o image_info: the image info.
 231 %
 232 %    o draw_info: the draw info.
 233 %
 234 */
 235 MagickExport DrawInfo *CloneDrawInfo(const ImageInfo *image_info,
 236   const DrawInfo *draw_info)
 237 {
 238   DrawInfo
 239     *clone_info;
 240
 241   ExceptionInfo
 242     *exception;
 243
 244   clone_info=(DrawInfo *) AcquireMagickMemory(sizeof(*clone_info));
 245   if (clone_info == (DrawInfo *) NULL)
 246     ThrowFatalException(ResourceLimitFatalError,"MemoryAllocationFailed");
 247   GetDrawInfo(image_info,clone_info);
 248   if (draw_info == (DrawInfo *) NULL)
 249     return(clone_info);
 250   exception=AcquireExceptionInfo();
 251   (void) CloneString(&clone_info->primitive,draw_info->primitive);
 252   (void) CloneString(&clone_info->geometry,draw_info->geometry);
 253   clone_info->viewbox=draw_info->viewbox;
 254   clone_info->affine=draw_info->affine;
 255   clone_info->gravity=draw_info->gravity;
 256   clone_info->fill=draw_info->fill;
 257   clone_info->stroke=draw_info->stroke;
 258   clone_info->stroke_width=draw_info->stroke_width;
 259   if (draw_info->fill_pattern != (Image *) NULL)
 260     clone_info->fill_pattern=CloneImage(draw_info->fill_pattern,0,0,MagickTrue,
 261       exception);
 262   if (draw_info->stroke_pattern != (Image *) NULL)
 263     clone_info->stroke_pattern=CloneImage(draw_info->stroke_pattern,0,0,
 264       MagickTrue,exception);
 265   clone_info->stroke_antialias=draw_info->stroke_antialias;
 266   clone_info->text_antialias=draw_info->text_antialias;
 267   clone_info->fill_rule=draw_info->fill_rule;
 268   clone_info->linecap=draw_info->linecap;
 269   clone_info->linejoin=draw_info->linejoin;
 270   clone_info->miterlimit=draw_info->miterlimit;
 271   clone_info->dash_offset=draw_info->dash_offset;
 272   clone_info->decorate=draw_info->decorate;
 273   clone_info->compose=draw_info->compose;
 274   (void) CloneString(&clone_info->text,draw_info->text);
 275   (void) CloneString(&clone_info->font,draw_info->font);
 276   (void) CloneString(&clone_info->metrics,draw_info->metrics);
 277   (void) CloneString(&clone_info->family,draw_info->family);
 278   clone_info->style=draw_info->style;
 279   clone_info->stretch=draw_info->stretch;
 280   clone_info->weight=draw_info->weight;
 281   (void) CloneString(&clone_info->encoding,draw_info->encoding);
 282   clone_info->pointsize=draw_info->pointsize;
 283   clone_info->kerning=draw_info->kerning;
 284   clone_info->interline_spacing=draw_info->interline_spacing;
 285   clone_info->interword_spacing=draw_info->interword_spacing;
 286   clone_info->direction=draw_info->direction;
 287   (void) CloneString(&clone_info->density,draw_info->density);
 288   clone_info->align=draw_info->align;
 289   clone_info->undercolor=draw_info->undercolor;
 290   clone_info->border_color=draw_info->border_color;
 291   (void) CloneString(&clone_info->server_name,draw_info->server_name);
 292   if (draw_info->dash_pattern != (double *) NULL)
 293     {
 294       register ssize_t
 295         x;
 296
 297       for (x=0; draw_info->dash_pattern[x] != 0.0; x++) ;
 298       clone_info->dash_pattern=(double *) AcquireQuantumMemory((size_t) x+1UL,
 299         sizeof(*clone_info->dash_pattern));
 300       if (clone_info->dash_pattern == (double *) NULL)
 301         ThrowFatalException(ResourceLimitFatalError,
 302           "UnableToAllocateDashPattern");
 303       (void) CopyMagickMemory(clone_info->dash_pattern,draw_info->dash_pattern,
 304         (size_t) (x+1)*sizeof(*clone_info->dash_pattern));
 305     }
 306   clone_info->gradient=draw_info->gradient;
 307   if (draw_info->gradient.stops != (StopInfo *) NULL)
 308     {
 309       size_t
 310         number_stops;
 311
 312       number_stops=clone_info->gradient.number_stops;
 313       clone_info->gradient.stops=(StopInfo *) AcquireQuantumMemory((size_t)
 314         number_stops,sizeof(*clone_info->gradient.stops));
 315       if (clone_info->gradient.stops == (StopInfo *) NULL)
 316         ThrowFatalException(ResourceLimitFatalError,
 317           "UnableToAllocateDashPattern");
 318       (void) CopyMagickMemory(clone_info->gradient.stops,
 319         draw_info->gradient.stops,(size_t) number_stops*
 320         sizeof(*clone_info->gradient.stops));
 321     }
 322   (void) CloneString(&clone_info->clip_mask,draw_info->clip_mask);
 323   clone_info->bounds=draw_info->bounds;
 324   clone_info->clip_units=draw_info->clip_units;
 325   clone_info->render=draw_info->render;
 326   clone_info->alpha=draw_info->alpha;
 327   clone_info->element_reference=draw_info->element_reference;
 328   clone_info->debug=IsEventLogging();
 329   exception=DestroyExceptionInfo(exception);
 330   return(clone_info);
 331 }
 332 \f
 333 /*
 334 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 335 %                                                                             %
 336 %                                                                             %
 337 %                                                                             %
 338 +   C o n v e r t P a t h T o P o l y g o n                                   %
 339 %                                                                             %
 340 %                                                                             %
 341 %                                                                             %
 342 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 343 %
 344 %  ConvertPathToPolygon() converts a path to the more efficient sorted
 345 %  rendering form.
 346 %
 347 %  The format of the ConvertPathToPolygon method is:
 348 %
 349 %      PolygonInfo *ConvertPathToPolygon(const DrawInfo *draw_info,
 350 %        const PathInfo *path_info)
 351 %
 352 %  A description of each parameter follows:
 353 %
 354 %    o Method ConvertPathToPolygon returns the path in a more efficient sorted
 355 %      rendering form of type PolygonInfo.
 356 %
 357 %    o draw_info: Specifies a pointer to an DrawInfo structure.
 358 %
 359 %    o path_info: Specifies a pointer to an PathInfo structure.
 360 %
 361 %
 362 */
 363
 364 #if defined(__cplusplus) || defined(c_plusplus)
 365 extern "C" {
 366 #endif
 367
 368 static int CompareEdges(const void *x,const void *y)
 369 {
 370   register const EdgeInfo
 371     *p,
 372     *q;
 373
 374   /*
 375     Compare two edges.
 376   */
 377   p=(const EdgeInfo *) x;
 378   q=(const EdgeInfo *) y;
 379   if ((p->points[0].y-MagickEpsilon) > q->points[0].y)
 380     return(1);
 381   if ((p->points[0].y+MagickEpsilon) < q->points[0].y)
 382     return(-1);
 383   if ((p->points[0].x-MagickEpsilon) > q->points[0].x)
 384     return(1);
 385   if ((p->points[0].x+MagickEpsilon) < q->points[0].x)
 386     return(-1);
 387   if (((p->points[1].x-p->points[0].x)*(q->points[1].y-q->points[0].y)-
 388        (p->points[1].y-p->points[0].y)*(q->points[1].x-q->points[0].x)) > 0.0)
 389     return(1);
 390   return(-1);
 391 }
 392
 393 #if defined(__cplusplus) || defined(c_plusplus)
 394 }
 395 #endif
 396
 397 static void LogPolygonInfo(const PolygonInfo *polygon_info)
 398 {
 399   register EdgeInfo
 400     *p;
 401
 402   register ssize_t
 403     i,
 404     j;
 405
 406   (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin active-edge");
 407   p=polygon_info->edges;
 408   for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
 409   {
 410     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"      edge %.20g:",
 411       (double) i);
 412     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"      direction: %s",
 413       p->direction != MagickFalse ? "down" : "up");
 414     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"      ghostline: %s",
 415       p->ghostline != MagickFalse ? "transparent" : "opaque");
 416     (void) LogMagickEvent(DrawEvent,GetMagickModule(),
 417       "      bounds: %g %g - %g %g",p->bounds.x1,p->bounds.y1,
 418       p->bounds.x2,p->bounds.y2);
 419     for (j=0; j < (ssize_t) p->number_points; j++)
 420       (void) LogMagickEvent(DrawEvent,GetMagickModule(),"        %g %g",
 421         p->points[j].x,p->points[j].y);
 422     p++;
 423   }
 424   (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end active-edge");
 425 }
 426
 427 static void ReversePoints(PointInfo *points,const size_t number_points)
 428 {
 429   PointInfo
 430     point;
 431
 432   register ssize_t
 433     i;
 434
 435   for (i=0; i < (ssize_t) (number_points >> 1); i++)
 436   {
 437     point=points[i];
 438     points[i]=points[number_points-(i+1)];
 439     points[number_points-(i+1)]=point;
 440   }
 441 }
 442
 443 static PolygonInfo *ConvertPathToPolygon(
 444   const DrawInfo *magick_unused(draw_info),const PathInfo *path_info)
 445 {
 446   long
 447     direction,
 448     next_direction;
 449
 450   PointInfo
 451     point,
 452     *points;
 453
 454   PolygonInfo
 455     *polygon_info;
 456
 457   SegmentInfo
 458     bounds;
 459
 460   register ssize_t
 461     i,
 462     n;
 463
 464   MagickBooleanType
 465     ghostline;
 466
 467   size_t
 468     edge,
 469     number_edges,
 470     number_points;
 471
 472   /*
 473     Convert a path to the more efficient sorted rendering form.
 474   */
 475   polygon_info=(PolygonInfo *) AcquireMagickMemory(sizeof(*polygon_info));
 476   if (polygon_info == (PolygonInfo *) NULL)
 477     return((PolygonInfo *) NULL);
 478   number_edges=16;
 479   polygon_info->edges=(EdgeInfo *) AcquireQuantumMemory((size_t) number_edges,
 480     sizeof(*polygon_info->edges));
 481   if (polygon_info->edges == (EdgeInfo *) NULL)
 482     return((PolygonInfo *) NULL);
 483   direction=0;
 484   edge=0;
 485   ghostline=MagickFalse;
 486   n=0;
 487   number_points=0;
 488   points=(PointInfo *) NULL;
 489   (void) ResetMagickMemory(&point,0,sizeof(point));
 490   (void) ResetMagickMemory(&bounds,0,sizeof(bounds));
 491   for (i=0; path_info[i].code != EndCode; i++)
 492   {
 493     if ((path_info[i].code == MoveToCode) || (path_info[i].code == OpenCode) ||
 494         (path_info[i].code == GhostlineCode))
 495       {
 496         /*
 497           Move to.
 498         */
 499         if ((points != (PointInfo *) NULL) && (n >= 2))
 500           {
 501             if (edge == number_edges)
 502               {
 503                 number_edges<<=1;
 504                 polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
 505                   polygon_info->edges,(size_t) number_edges,
 506                   sizeof(*polygon_info->edges));
 507                 if (polygon_info->edges == (EdgeInfo *) NULL)
 508                   return((PolygonInfo *) NULL);
 509               }
 510             polygon_info->edges[edge].number_points=(size_t) n;
 511             polygon_info->edges[edge].scanline=(-1.0);
 512             polygon_info->edges[edge].highwater=0;
 513             polygon_info->edges[edge].ghostline=ghostline;
 514             polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
 515             if (direction < 0)
 516               ReversePoints(points,(size_t) n);
 517             polygon_info->edges[edge].points=points;
 518             polygon_info->edges[edge].bounds=bounds;
 519             polygon_info->edges[edge].bounds.y1=points[0].y;
 520             polygon_info->edges[edge].bounds.y2=points[n-1].y;
 521             points=(PointInfo *) NULL;
 522             ghostline=MagickFalse;
 523             edge++;
 524           }
 525         if (points == (PointInfo *) NULL)
 526           {
 527             number_points=16;
 528             points=(PointInfo *) AcquireQuantumMemory((size_t) number_points,
 529               sizeof(*points));
 530             if (points == (PointInfo *) NULL)
 531               return((PolygonInfo *) NULL);
 532           }
 533         ghostline=path_info[i].code == GhostlineCode ? MagickTrue : MagickFalse;
 534         point=path_info[i].point;
 535         points[0]=point;
 536         bounds.x1=point.x;
 537         bounds.x2=point.x;
 538         direction=0;
 539         n=1;
 540         continue;
 541       }
 542     /*
 543       Line to.
 544     */
 545     next_direction=((path_info[i].point.y > point.y) ||
 546       ((path_info[i].point.y == point.y) &&
 547        (path_info[i].point.x > point.x))) ? 1 : -1;
 548     if ((direction != 0) && (direction != next_direction))
 549       {
 550         /*
 551           New edge.
 552         */
 553         point=points[n-1];
 554         if (edge == number_edges)
 555           {
 556             number_edges<<=1;
 557             polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
 558               polygon_info->edges,(size_t) number_edges,
 559               sizeof(*polygon_info->edges));
 560             if (polygon_info->edges == (EdgeInfo *) NULL)
 561               return((PolygonInfo *) NULL);
 562           }
 563         polygon_info->edges[edge].number_points=(size_t) n;
 564         polygon_info->edges[edge].scanline=(-1.0);
 565         polygon_info->edges[edge].highwater=0;
 566         polygon_info->edges[edge].ghostline=ghostline;
 567         polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
 568         if (direction < 0)
 569           ReversePoints(points,(size_t) n);
 570         polygon_info->edges[edge].points=points;
 571         polygon_info->edges[edge].bounds=bounds;
 572         polygon_info->edges[edge].bounds.y1=points[0].y;
 573         polygon_info->edges[edge].bounds.y2=points[n-1].y;
 574         number_points=16;
 575         points=(PointInfo *) AcquireQuantumMemory((size_t) number_points,
 576           sizeof(*points));
 577         if (points == (PointInfo *) NULL)
 578           return((PolygonInfo *) NULL);
 579         n=1;
 580         ghostline=MagickFalse;
 581         points[0]=point;
 582         bounds.x1=point.x;
 583         bounds.x2=point.x;
 584         edge++;
 585       }
 586     direction=next_direction;
 587     if (points == (PointInfo *) NULL)
 588       continue;
 589     if (n == (ssize_t) number_points)
 590       {
 591         number_points<<=1;
 592         points=(PointInfo *) ResizeQuantumMemory(points,(size_t) number_points,
 593           sizeof(*points));
 594         if (points == (PointInfo *) NULL)
 595           return((PolygonInfo *) NULL);
 596       }
 597     point=path_info[i].point;
 598     points[n]=point;
 599     if (point.x < bounds.x1)
 600       bounds.x1=point.x;
 601     if (point.x > bounds.x2)
 602       bounds.x2=point.x;
 603     n++;
 604   }
 605   if (points != (PointInfo *) NULL)
 606     {
 607       if (n < 2)
 608         points=(PointInfo *) RelinquishMagickMemory(points);
 609       else
 610         {
 611           if (edge == number_edges)
 612             {
 613               number_edges<<=1;
 614               polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
 615                 polygon_info->edges,(size_t) number_edges,
 616                 sizeof(*polygon_info->edges));
 617               if (polygon_info->edges == (EdgeInfo *) NULL)
 618                 return((PolygonInfo *) NULL);
 619             }
 620           polygon_info->edges[edge].number_points=(size_t) n;
 621           polygon_info->edges[edge].scanline=(-1.0);
 622           polygon_info->edges[edge].highwater=0;
 623           polygon_info->edges[edge].ghostline=ghostline;
 624           polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
 625           if (direction < 0)
 626             ReversePoints(points,(size_t) n);
 627           polygon_info->edges[edge].points=points;
 628           polygon_info->edges[edge].bounds=bounds;
 629           polygon_info->edges[edge].bounds.y1=points[0].y;
 630           polygon_info->edges[edge].bounds.y2=points[n-1].y;
 631           ghostline=MagickFalse;
 632           edge++;
 633         }
 634     }
 635   polygon_info->number_edges=edge;
 636   qsort(polygon_info->edges,(size_t) polygon_info->number_edges,
 637     sizeof(*polygon_info->edges),CompareEdges);
 638   if (IsEventLogging() != MagickFalse)
 639     LogPolygonInfo(polygon_info);
 640   return(polygon_info);
 641 }
 642 \f
 643 /*
 644 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 645 %                                                                             %
 646 %                                                                             %
 647 %                                                                             %
 648 +   C o n v e r t P r i m i t i v e T o P a t h                               %
 649 %                                                                             %
 650 %                                                                             %
 651 %                                                                             %
 652 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 653 %
 654 %  ConvertPrimitiveToPath() converts a PrimitiveInfo structure into a vector
 655 %  path structure.
 656 %
 657 %  The format of the ConvertPrimitiveToPath method is:
 658 %
 659 %      PathInfo *ConvertPrimitiveToPath(const DrawInfo *draw_info,
 660 %        const PrimitiveInfo *primitive_info)
 661 %
 662 %  A description of each parameter follows:
 663 %
 664 %    o Method ConvertPrimitiveToPath returns a vector path structure of type
 665 %      PathInfo.
 666 %
 667 %    o draw_info: a structure of type DrawInfo.
 668 %
 669 %    o primitive_info: Specifies a pointer to an PrimitiveInfo structure.
 670 %
 671 %
 672 */
 673
 674 static void LogPathInfo(const PathInfo *path_info)
 675 {
 676   register const PathInfo
 677     *p;
 678
 679   (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin vector-path");
 680   for (p=path_info; p->code != EndCode; p++)
 681     (void) LogMagickEvent(DrawEvent,GetMagickModule(),
 682       "      %g %g %s",p->point.x,p->point.y,p->code == GhostlineCode ?
 683       "moveto ghostline" : p->code == OpenCode ? "moveto open" :
 684       p->code == MoveToCode ? "moveto" : p->code == LineToCode ? "lineto" :
 685       "?");
 686   (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end vector-path");
 687 }
 688
 689 static PathInfo *ConvertPrimitiveToPath(
 690   const DrawInfo *magick_unused(draw_info),const PrimitiveInfo *primitive_info)
 691 {
 692   PathInfo
 693     *path_info;
 694
 695   PathInfoCode
 696     code;
 697
 698   PointInfo
 699     p,
 700     q;
 701
 702   register ssize_t
 703     i,
 704     n;
 705
 706   ssize_t
 707     coordinates,
 708     start;
 709
 710   /*
 711     Converts a PrimitiveInfo structure into a vector path structure.
 712   */
 713   switch (primitive_info->primitive)
 714   {
 715     case PointPrimitive:
 716     case ColorPrimitive:
 717     case MattePrimitive:
 718     case TextPrimitive:
 719     case ImagePrimitive:
 720       return((PathInfo *) NULL);
 721     default:
 722       break;
 723   }
 724   for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
 725   path_info=(PathInfo *) AcquireQuantumMemory((size_t) (2UL*i+3UL),
 726     sizeof(*path_info));
 727   if (path_info == (PathInfo *) NULL)
 728     return((PathInfo *) NULL);
 729   coordinates=0;
 730   n=0;
 731   p.x=(-1.0);
 732   p.y=(-1.0);
 733   q.x=(-1.0);
 734   q.y=(-1.0);
 735   start=0;
 736   for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
 737   {
 738     code=LineToCode;
 739     if (coordinates <= 0)
 740       {
 741         coordinates=(ssize_t) primitive_info[i].coordinates;
 742         p=primitive_info[i].point;
 743         start=n;
 744         code=MoveToCode;
 745       }
 746     coordinates--;
 747     /*
 748       Eliminate duplicate points.
 749     */
 750     if ((i == 0) || (fabs(q.x-primitive_info[i].point.x) > MagickEpsilon) ||
 751         (fabs(q.y-primitive_info[i].point.y) > MagickEpsilon))
 752       {
 753         path_info[n].code=code;
 754         path_info[n].point=primitive_info[i].point;
 755         q=primitive_info[i].point;
 756         n++;
 757       }
 758     if (coordinates > 0)
 759       continue;
 760     if ((fabs(p.x-primitive_info[i].point.x) <= MagickEpsilon) &&
 761         (fabs(p.y-primitive_info[i].point.y) <= MagickEpsilon))
 762       continue;
 763     /*
 764       Mark the p point as open if it does not match the q.
 765     */
 766     path_info[start].code=OpenCode;
 767     path_info[n].code=GhostlineCode;
 768     path_info[n].point=primitive_info[i].point;
 769     n++;
 770     path_info[n].code=LineToCode;
 771     path_info[n].point=p;
 772     n++;
 773   }
 774   path_info[n].code=EndCode;
 775   path_info[n].point.x=0.0;
 776   path_info[n].point.y=0.0;
 777   if (IsEventLogging() != MagickFalse)
 778     LogPathInfo(path_info);
 779   return(path_info);
 780 }
 781 \f
 782 /*
 783 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 784 %                                                                             %
 785 %                                                                             %
 786 %                                                                             %
 787 %   D e s t r o y D r a w I n f o                                             %
 788 %                                                                             %
 789 %                                                                             %
 790 %                                                                             %
 791 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 792 %
 793 %  DestroyDrawInfo() deallocates memory associated with an DrawInfo
 794 %  structure.
 795 %
 796 %  The format of the DestroyDrawInfo method is:
 797 %
 798 %      DrawInfo *DestroyDrawInfo(DrawInfo *draw_info)
 799 %
 800 %  A description of each parameter follows:
 801 %
 802 %    o draw_info: the draw info.
 803 %
 804 */
 805 MagickExport DrawInfo *DestroyDrawInfo(DrawInfo *draw_info)
 806 {
 807   if (draw_info->debug != MagickFalse)
 808     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
 809   assert(draw_info != (DrawInfo *) NULL);
 810   assert(draw_info->signature == MagickSignature);
 811   if (draw_info->primitive != (char *) NULL)
 812     draw_info->primitive=DestroyString(draw_info->primitive);
 813   if (draw_info->text != (char *) NULL)
 814     draw_info->text=DestroyString(draw_info->text);
 815   if (draw_info->geometry != (char *) NULL)
 816     draw_info->geometry=DestroyString(draw_info->geometry);
 817   if (draw_info->fill_pattern != (Image *) NULL)
 818     draw_info->fill_pattern=DestroyImage(draw_info->fill_pattern);
 819   if (draw_info->stroke_pattern != (Image *) NULL)
 820     draw_info->stroke_pattern=DestroyImage(draw_info->stroke_pattern);
 821   if (draw_info->font != (char *) NULL)
 822     draw_info->font=DestroyString(draw_info->font);
 823   if (draw_info->metrics != (char *) NULL)
 824     draw_info->metrics=DestroyString(draw_info->metrics);
 825   if (draw_info->family != (char *) NULL)
 826     draw_info->family=DestroyString(draw_info->family);
 827   if (draw_info->encoding != (char *) NULL)
 828     draw_info->encoding=DestroyString(draw_info->encoding);
 829   if (draw_info->density != (char *) NULL)
 830     draw_info->density=DestroyString(draw_info->density);
 831   if (draw_info->server_name != (char *) NULL)
 832     draw_info->server_name=(char *)
 833      RelinquishMagickMemory(draw_info->server_name);
 834   if (draw_info->dash_pattern != (double *) NULL)
 835     draw_info->dash_pattern=(double *) RelinquishMagickMemory(
 836       draw_info->dash_pattern);
 837   if (draw_info->gradient.stops != (StopInfo *) NULL)
 838     draw_info->gradient.stops=(StopInfo *) RelinquishMagickMemory(
 839       draw_info->gradient.stops);
 840   if (draw_info->clip_mask != (char *) NULL)
 841     draw_info->clip_mask=DestroyString(draw_info->clip_mask);
 842   draw_info->signature=(~MagickSignature);
 843   draw_info=(DrawInfo *) RelinquishMagickMemory(draw_info);
 844   return(draw_info);
 845 }
 846 \f
 847 /*
 848 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 849 %                                                                             %
 850 %                                                                             %
 851 %                                                                             %
 852 +   D e s t r o y E d g e                                                     %
 853 %                                                                             %
 854 %                                                                             %
 855 %                                                                             %
 856 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 857 %
 858 %  DestroyEdge() destroys the specified polygon edge.
 859 %
 860 %  The format of the DestroyEdge method is:
 861 %
 862 %      ssize_t DestroyEdge(PolygonInfo *polygon_info,const int edge)
 863 %
 864 %  A description of each parameter follows:
 865 %
 866 %    o polygon_info: Specifies a pointer to an PolygonInfo structure.
 867 %
 868 %    o edge: the polygon edge number to destroy.
 869 %
 870 */
 871 static size_t DestroyEdge(PolygonInfo *polygon_info,
 872   const size_t edge)
 873 {
 874   assert(edge < polygon_info->number_edges);
 875   polygon_info->edges[edge].points=(PointInfo *) RelinquishMagickMemory(
 876     polygon_info->edges[edge].points);
 877   polygon_info->number_edges--;
 878   if (edge < polygon_info->number_edges)
 879     (void) CopyMagickMemory(polygon_info->edges+edge,polygon_info->edges+edge+1,
 880       (size_t) (polygon_info->number_edges-edge)*sizeof(*polygon_info->edges));
 881   return(polygon_info->number_edges);
 882 }
 883 \f
 884 /*
 885 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 886 %                                                                             %
 887 %                                                                             %
 888 %                                                                             %
 889 +   D e s t r o y P o l y g o n I n f o                                       %
 890 %                                                                             %
 891 %                                                                             %
 892 %                                                                             %
 893 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 894 %
 895 %  DestroyPolygonInfo() destroys the PolygonInfo data structure.
 896 %
 897 %  The format of the DestroyPolygonInfo method is:
 898 %
 899 %      PolygonInfo *DestroyPolygonInfo(PolygonInfo *polygon_info)
 900 %
 901 %  A description of each parameter follows:
 902 %
 903 %    o polygon_info: Specifies a pointer to an PolygonInfo structure.
 904 %
 905 */
 906 static PolygonInfo *DestroyPolygonInfo(PolygonInfo *polygon_info)
 907 {
 908   register ssize_t
 909     i;
 910
 911   for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
 912     polygon_info->edges[i].points=(PointInfo *)
 913       RelinquishMagickMemory(polygon_info->edges[i].points);
 914   polygon_info->edges=(EdgeInfo *) RelinquishMagickMemory(polygon_info->edges);
 915   return((PolygonInfo *) RelinquishMagickMemory(polygon_info));
 916 }
 917 \f
 918 /*
 919 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 920 %                                                                             %
 921 %                                                                             %
 922 %                                                                             %
 923 %     D r a w A f f i n e I m a g e                                           %
 924 %                                                                             %
 925 %                                                                             %
 926 %                                                                             %
 927 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 928 %
 929 %  DrawAffineImage() composites the source over the destination image as
 930 %  dictated by the affine transform.
 931 %
 932 %  The format of the DrawAffineImage method is:
 933 %
 934 %      MagickBooleanType DrawAffineImage(Image *image,const Image *source,
 935 %        const AffineMatrix *affine,ExceptionInfo *exception)
 936 %
 937 %  A description of each parameter follows:
 938 %
 939 %    o image: the image.
 940 %
 941 %    o source: the source image.
 942 %
 943 %    o affine: the affine transform.
 944 %
 945 %    o exception: return any errors or warnings in this structure.
 946 %
 947 */
 948
 949 static SegmentInfo AffineEdge(const Image *image,const AffineMatrix *affine,
 950   const double y,const SegmentInfo *edge)
 951 {
 952   double
 953     intercept,
 954     z;
 955
 956   register double
 957     x;
 958
 959   SegmentInfo
 960     inverse_edge;
 961
 962   /*
 963     Determine left and right edges.
 964   */
 965   inverse_edge.x1=edge->x1;
 966   inverse_edge.y1=edge->y1;
 967   inverse_edge.x2=edge->x2;
 968   inverse_edge.y2=edge->y2;
 969   z=affine->ry*y+affine->tx;
 970   if (affine->sx > MagickEpsilon)
 971     {
 972       intercept=(-z/affine->sx);
 973       x=intercept;
 974       if (x > inverse_edge.x1)
 975         inverse_edge.x1=x;
 976       intercept=(-z+(double) image->columns)/affine->sx;
 977       x=intercept;
 978       if (x < inverse_edge.x2)
 979         inverse_edge.x2=x;
 980     }
 981   else
 982     if (affine->sx < -MagickEpsilon)
 983       {
 984         intercept=(-z+(double) image->columns)/affine->sx;
 985         x=intercept;
 986         if (x > inverse_edge.x1)
 987           inverse_edge.x1=x;
 988         intercept=(-z/affine->sx);
 989         x=intercept;
 990         if (x < inverse_edge.x2)
 991           inverse_edge.x2=x;
 992       }
 993     else
 994       if ((z < 0.0) || ((size_t) floor(z+0.5) >= image->columns))
 995         {
 996           inverse_edge.x2=edge->x1;
 997           return(inverse_edge);
 998         }
 999   /*
1000     Determine top and bottom edges.
1001   */
1002   z=affine->sy*y+affine->ty;
1003   if (affine->rx > MagickEpsilon)
1004     {
1005       intercept=(-z/affine->rx);
1006       x=intercept;
1007       if (x > inverse_edge.x1)
1008         inverse_edge.x1=x;
1009       intercept=(-z+(double) image->rows)/affine->rx;
1010       x=intercept;
1011       if (x < inverse_edge.x2)
1012         inverse_edge.x2=x;
1013     }
1014   else
1015     if (affine->rx < -MagickEpsilon)
1016       {
1017         intercept=(-z+(double) image->rows)/affine->rx;
1018         x=intercept;
1019         if (x > inverse_edge.x1)
1020           inverse_edge.x1=x;
1021         intercept=(-z/affine->rx);
1022         x=intercept;
1023         if (x < inverse_edge.x2)
1024           inverse_edge.x2=x;
1025       }
1026     else
1027       if ((z < 0.0) || ((size_t) floor(z+0.5) >= image->rows))
1028         {
1029           inverse_edge.x2=edge->x2;
1030           return(inverse_edge);
1031         }
1032   return(inverse_edge);
1033 }
1034
1035 static AffineMatrix InverseAffineMatrix(const AffineMatrix *affine)
1036 {
1037   AffineMatrix
1038     inverse_affine;
1039
1040   double
1041     determinant;
1042
1043   determinant=1.0/(affine->sx*affine->sy-affine->rx*affine->ry);
1044   inverse_affine.sx=determinant*affine->sy;
1045   inverse_affine.rx=determinant*(-affine->rx);
1046   inverse_affine.ry=determinant*(-affine->ry);
1047   inverse_affine.sy=determinant*affine->sx;
1048   inverse_affine.tx=(-affine->tx)*inverse_affine.sx-affine->ty*
1049     inverse_affine.ry;
1050   inverse_affine.ty=(-affine->tx)*inverse_affine.rx-affine->ty*
1051     inverse_affine.sy;
1052   return(inverse_affine);
1053 }
1054
1055 static inline ssize_t MagickAbsoluteValue(const ssize_t x)
1056 {
1057   if (x < 0)
1058     return(-x);
1059   return(x);
1060 }
1061
1062 static inline double MagickMax(const double x,const double y)
1063 {
1064   if (x > y)
1065     return(x);
1066   return(y);
1067 }
1068
1069 static inline double MagickMin(const double x,const double y)
1070 {
1071   if (x < y)
1072     return(x);
1073   return(y);
1074 }
1075
1076 MagickExport MagickBooleanType DrawAffineImage(Image *image,
1077   const Image *source,const AffineMatrix *affine,ExceptionInfo *exception)
1078 {
1079   AffineMatrix
1080     inverse_affine;
1081
1082   CacheView
1083     *image_view,
1084     *source_view;
1085
1086   MagickBooleanType
1087     status;
1088
1089   PixelInfo
1090     zero;
1091
1092   PointInfo
1093     extent[4],
1094     min,
1095     max,
1096     point;
1097
1098   register ssize_t
1099     i;
1100
1101   SegmentInfo
1102     edge;
1103
1104   ssize_t
1105     start,
1106     stop,
1107     y;
1108
1109   /*
1110     Determine bounding box.
1111   */
1112   assert(image != (Image *) NULL);
1113   assert(image->signature == MagickSignature);
1114   if (image->debug != MagickFalse)
1115     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1116   assert(source != (const Image *) NULL);
1117   assert(source->signature == MagickSignature);
1118   assert(affine != (AffineMatrix *) NULL);
1119   extent[0].x=0.0;
1120   extent[0].y=0.0;
1121   extent[1].x=(double) source->columns-1.0;
1122   extent[1].y=0.0;
1123   extent[2].x=(double) source->columns-1.0;
1124   extent[2].y=(double) source->rows-1.0;
1125   extent[3].x=0.0;
1126   extent[3].y=(double) source->rows-1.0;
1127   for (i=0; i < 4; i++)
1128   {
1129     point=extent[i];
1130     extent[i].x=point.x*affine->sx+point.y*affine->ry+affine->tx;
1131     extent[i].y=point.x*affine->rx+point.y*affine->sy+affine->ty;
1132   }
1133   min=extent[0];
1134   max=extent[0];
1135   for (i=1; i < 4; i++)
1136   {
1137     if (min.x > extent[i].x)
1138       min.x=extent[i].x;
1139     if (min.y > extent[i].y)
1140       min.y=extent[i].y;
1141     if (max.x < extent[i].x)
1142       max.x=extent[i].x;
1143     if (max.y < extent[i].y)
1144       max.y=extent[i].y;
1145   }
1146   /*
1147     Affine transform image.
1148   */
1149   if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
1150     return(MagickFalse);
1151   status=MagickTrue;
1152   edge.x1=MagickMax(min.x,0.0);
1153   edge.y1=MagickMax(min.y,0.0);
1154   edge.x2=MagickMin(max.x,(double) image->columns-1.0);
1155   edge.y2=MagickMin(max.y,(double) image->rows-1.0);
1156   inverse_affine=InverseAffineMatrix(affine);
1157   GetPixelInfo(image,&zero);
1158   start=(ssize_t) ceil(edge.y1-0.5);
1159   stop=(ssize_t) ceil(edge.y2-0.5);
1160   image_view=AcquireCacheView(image);
1161   source_view=AcquireCacheView(source);
1162 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1163   #pragma omp parallel for schedule(static) shared(status)
1164 #endif
1165   for (y=start; y <= stop; y++)
1166   {
1167     PixelInfo
1168       composite,
1169       pixel;
1170
1171     PointInfo
1172       point;
1173
1174     register ssize_t
1175       x;
1176
1177     register Quantum
1178       *restrict q;
1179
1180     SegmentInfo
1181       inverse_edge;
1182
1183     ssize_t
1184       x_offset;
1185
1186     inverse_edge=AffineEdge(source,&inverse_affine,(double) y,&edge);
1187     if (inverse_edge.x2 < inverse_edge.x1)
1188       continue;
1189     q=GetCacheViewAuthenticPixels(image_view,(ssize_t) ceil(inverse_edge.x1-
1190       0.5),y,(size_t) ((ssize_t) floor(inverse_edge.x2+0.5)-(ssize_t) floor(
1191       inverse_edge.x1+0.5)+1),1,exception);
1192     if (q == (Quantum *) NULL)
1193       continue;
1194     pixel=zero;
1195     composite=zero;
1196     x_offset=0;
1197     for (x=(ssize_t) ceil(inverse_edge.x1-0.5); x <= (ssize_t) floor(inverse_edge.x2+0.5); x++)
1198     {
1199       point.x=(double) x*inverse_affine.sx+y*inverse_affine.ry+
1200         inverse_affine.tx;
1201       point.y=(double) x*inverse_affine.rx+y*inverse_affine.sy+
1202         inverse_affine.ty;
1203       (void) InterpolatePixelInfo(source,source_view,UndefinedInterpolatePixel,
1204         point.x,point.y,&pixel,exception);
1205       GetPixelInfoPixel(image,q,&composite);
1206       CompositePixelInfoOver(&pixel,pixel.alpha,&composite,composite.alpha,
1207         &composite);
1208       SetPixelInfoPixel(image,&composite,q);
1209       x_offset++;
1210       q+=GetPixelChannels(image);
1211     }
1212     if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1213       status=MagickFalse;
1214   }
1215   source_view=DestroyCacheView(source_view);
1216   image_view=DestroyCacheView(image_view);
1217   return(status);
1218 }
1219 \f
1220 /*
1221 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1222 %                                                                             %
1223 %                                                                             %
1224 %                                                                             %
1225 +   D r a w B o u n d i n g R e c t a n g l e s                               %
1226 %                                                                             %
1227 %                                                                             %
1228 %                                                                             %
1229 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1230 %
1231 %  DrawBoundingRectangles() draws the bounding rectangles on the image.  This
1232 %  is only useful for developers debugging the rendering algorithm.
1233 %
1234 %  The format of the DrawBoundingRectangles method is:
1235 %
1236 %      void DrawBoundingRectangles(Image *image,const DrawInfo *draw_info,
1237 %        PolygonInfo *polygon_info,ExceptionInfo *exception)
1238 %
1239 %  A description of each parameter follows:
1240 %
1241 %    o image: the image.
1242 %
1243 %    o draw_info: the draw info.
1244 %
1245 %    o polygon_info: Specifies a pointer to a PolygonInfo structure.
1246 %
1247 %    o exception: return any errors or warnings in this structure.
1248 %
1249 */
1250 static void DrawBoundingRectangles(Image *image,const DrawInfo *draw_info,
1251   const PolygonInfo *polygon_info,ExceptionInfo *exception)
1252 {
1253   DrawInfo
1254     *clone_info;
1255
1256   MagickRealType
1257     mid;
1258
1259   PointInfo
1260     end,
1261     resolution,
1262     start;
1263
1264   PrimitiveInfo
1265     primitive_info[6];
1266
1267   register ssize_t
1268     i;
1269
1270   SegmentInfo
1271     bounds;
1272
1273   ssize_t
1274     coordinates;
1275
1276   clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1277   (void) QueryColorCompliance("#0000",AllCompliance,&clone_info->fill,
1278     exception);
1279   resolution.x=DefaultResolution;
1280   resolution.y=DefaultResolution;
1281   if (clone_info->density != (char *) NULL)
1282     {
1283       GeometryInfo
1284         geometry_info;
1285
1286       MagickStatusType
1287         flags;
1288
1289       flags=ParseGeometry(clone_info->density,&geometry_info);
1290       resolution.x=geometry_info.rho;
1291       resolution.y=geometry_info.sigma;
1292       if ((flags & SigmaValue) == MagickFalse)
1293         resolution.y=resolution.x;
1294     }
1295   mid=(resolution.x/72.0)*ExpandAffine(&clone_info->affine)*
1296     clone_info->stroke_width/2.0;
1297   bounds.x1=0.0;
1298   bounds.y1=0.0;
1299   bounds.x2=0.0;
1300   bounds.y2=0.0;
1301   if (polygon_info != (PolygonInfo *) NULL)
1302     {
1303       bounds=polygon_info->edges[0].bounds;
1304       for (i=1; i < (ssize_t) polygon_info->number_edges; i++)
1305       {
1306         if (polygon_info->edges[i].bounds.x1 < (double) bounds.x1)
1307           bounds.x1=polygon_info->edges[i].bounds.x1;
1308         if (polygon_info->edges[i].bounds.y1 < (double) bounds.y1)
1309           bounds.y1=polygon_info->edges[i].bounds.y1;
1310         if (polygon_info->edges[i].bounds.x2 > (double) bounds.x2)
1311           bounds.x2=polygon_info->edges[i].bounds.x2;
1312         if (polygon_info->edges[i].bounds.y2 > (double) bounds.y2)
1313           bounds.y2=polygon_info->edges[i].bounds.y2;
1314       }
1315       bounds.x1-=mid;
1316       bounds.x1=bounds.x1 < 0.0 ? 0.0 : bounds.x1 >= (double)
1317         image->columns ? (double) image->columns-1 : bounds.x1;
1318       bounds.y1-=mid;
1319       bounds.y1=bounds.y1 < 0.0 ? 0.0 : bounds.y1 >= (double)
1320         image->rows ? (double) image->rows-1 : bounds.y1;
1321       bounds.x2+=mid;
1322       bounds.x2=bounds.x2 < 0.0 ? 0.0 : bounds.x2 >= (double)
1323         image->columns ? (double) image->columns-1 : bounds.x2;
1324       bounds.y2+=mid;
1325       bounds.y2=bounds.y2 < 0.0 ? 0.0 : bounds.y2 >= (double)
1326         image->rows ? (double) image->rows-1 : bounds.y2;
1327       for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
1328       {
1329         if (polygon_info->edges[i].direction != 0)
1330           (void) QueryColorCompliance("red",AllCompliance,&clone_info->stroke,
1331             exception);
1332         else
1333           (void) QueryColorCompliance("green",AllCompliance,&clone_info->stroke,
1334             exception);
1335         start.x=(double) (polygon_info->edges[i].bounds.x1-mid);
1336         start.y=(double) (polygon_info->edges[i].bounds.y1-mid);
1337         end.x=(double) (polygon_info->edges[i].bounds.x2+mid);
1338         end.y=(double) (polygon_info->edges[i].bounds.y2+mid);
1339         primitive_info[0].primitive=RectanglePrimitive;
1340         TraceRectangle(primitive_info,start,end);
1341         primitive_info[0].method=ReplaceMethod;
1342         coordinates=(ssize_t) primitive_info[0].coordinates;
1343         primitive_info[coordinates].primitive=UndefinedPrimitive;
1344         (void) DrawPrimitive(image,clone_info,primitive_info,exception);
1345       }
1346     }
1347   (void) QueryColorCompliance("blue",AllCompliance,&clone_info->stroke,
1348     exception);
1349   start.x=(double) (bounds.x1-mid);
1350   start.y=(double) (bounds.y1-mid);
1351   end.x=(double) (bounds.x2+mid);
1352   end.y=(double) (bounds.y2+mid);
1353   primitive_info[0].primitive=RectanglePrimitive;
1354   TraceRectangle(primitive_info,start,end);
1355   primitive_info[0].method=ReplaceMethod;
1356   coordinates=(ssize_t) primitive_info[0].coordinates;
1357   primitive_info[coordinates].primitive=UndefinedPrimitive;
1358   (void) DrawPrimitive(image,clone_info,primitive_info,exception);
1359   clone_info=DestroyDrawInfo(clone_info);
1360 }
1361 \f
1362 /*
1363 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1364 %                                                                             %
1365 %                                                                             %
1366 %                                                                             %
1367 %   D r a w C l i p P a t h                                                   %
1368 %                                                                             %
1369 %                                                                             %
1370 %                                                                             %
1371 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1372 %
1373 %  DrawClipPath() draws the clip path on the image mask.
1374 %
1375 %  The format of the DrawClipPath method is:
1376 %
1377 %      MagickBooleanType DrawClipPath(Image *image,const DrawInfo *draw_info,
1378 %        const char *name,ExceptionInfo *exception)
1379 %
1380 %  A description of each parameter follows:
1381 %
1382 %    o image: the image.
1383 %
1384 %    o draw_info: the draw info.
1385 %
1386 %    o name: the name of the clip path.
1387 %
1388 %    o exception: return any errors or warnings in this structure.
1389 %
1390 */
1391 MagickExport MagickBooleanType DrawClipPath(Image *image,
1392   const DrawInfo *draw_info,const char *name,ExceptionInfo *exception)
1393 {
1394   char
1395     filename[MaxTextExtent];
1396
1397   Image
1398     *clip_mask;
1399
1400   const char
1401     *value;
1402
1403   DrawInfo
1404     *clone_info;
1405
1406   MagickStatusType
1407     status;
1408
1409   assert(image != (Image *) NULL);
1410   assert(image->signature == MagickSignature);
1411   if (image->debug != MagickFalse)
1412     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1413   assert(draw_info != (const DrawInfo *) NULL);
1414   (void) FormatLocaleString(filename,MaxTextExtent,"%s",name);
1415   value=GetImageArtifact(image,filename);
1416   if (value == (const char *) NULL)
1417     return(MagickFalse);
1418   clip_mask=CloneImage(image,image->columns,image->rows,MagickTrue,exception);
1419   if (clip_mask == (Image *) NULL)
1420     return(MagickFalse);
1421   (void) QueryColorCompliance("#0000",AllCompliance,
1422     &clip_mask->background_color,exception);
1423   clip_mask->background_color.alpha=(Quantum) TransparentAlpha;
1424   (void) SetImageBackgroundColor(clip_mask,exception);
1425   if (image->debug != MagickFalse)
1426     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"\nbegin clip-path %s",
1427       draw_info->clip_mask);
1428   clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1429   (void) CloneString(&clone_info->primitive,value);
1430   (void) QueryColorCompliance("#ffffff",AllCompliance,&clone_info->fill,
1431     exception);
1432   clone_info->clip_mask=(char *) NULL;
1433   status=NegateImage(clip_mask,MagickFalse,exception);
1434   (void) SetImageMask(image,clip_mask,exception);
1435   clip_mask=DestroyImage(clip_mask);
1436   clone_info=DestroyDrawInfo(clone_info);
1437   status=DrawImage(image,clone_info,exception);
1438   if (image->debug != MagickFalse)
1439     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end clip-path");
1440   return(status != 0 ? MagickTrue : MagickFalse);
1441 }
1442 \f
1443 /*
1444 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1445 %                                                                             %
1446 %                                                                             %
1447 %                                                                             %
1448 +   D r a w D a s h P o l y g o n                                             %
1449 %                                                                             %
1450 %                                                                             %
1451 %                                                                             %
1452 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1453 %
1454 %  DrawDashPolygon() draws a dashed polygon (line, rectangle, ellipse) on the
1455 %  image while respecting the dash offset and dash pattern attributes.
1456 %
1457 %  The format of the DrawDashPolygon method is:
1458 %
1459 %      MagickBooleanType DrawDashPolygon(const DrawInfo *draw_info,
1460 %        const PrimitiveInfo *primitive_info,Image *image,
1461 %        ExceptionInfo *exception)
1462 %
1463 %  A description of each parameter follows:
1464 %
1465 %    o draw_info: the draw info.
1466 %
1467 %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
1468 %
1469 %    o image: the image.
1470 %
1471 %    o exception: return any errors or warnings in this structure.
1472 %
1473 */
1474 static MagickBooleanType DrawDashPolygon(const DrawInfo *draw_info,
1475   const PrimitiveInfo *primitive_info,Image *image,ExceptionInfo *exception)
1476 {
1477   DrawInfo
1478     *clone_info;
1479
1480   MagickRealType
1481     length,
1482     maximum_length,
1483     offset,
1484     scale,
1485     total_length;
1486
1487   MagickStatusType
1488     status;
1489
1490   PrimitiveInfo
1491     *dash_polygon;
1492
1493   register ssize_t
1494     i;
1495
1496   register MagickRealType
1497     dx,
1498     dy;
1499
1500   size_t
1501     number_vertices;
1502
1503   ssize_t
1504     j,
1505     n;
1506
1507   assert(draw_info != (const DrawInfo *) NULL);
1508   if (image->debug != MagickFalse)
1509     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin draw-dash");
1510   clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1511   clone_info->miterlimit=0;
1512   for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
1513   number_vertices=(size_t) i;
1514   dash_polygon=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
1515     (2UL*number_vertices+1UL),sizeof(*dash_polygon));
1516   if (dash_polygon == (PrimitiveInfo *) NULL)
1517     return(MagickFalse);
1518   dash_polygon[0]=primitive_info[0];
1519   scale=ExpandAffine(&draw_info->affine);
1520   length=scale*(draw_info->dash_pattern[0]-0.5);
1521   offset=draw_info->dash_offset != 0.0 ? scale*draw_info->dash_offset : 0.0;
1522   j=1;
1523   for (n=0; offset > 0.0; j=0)
1524   {
1525     if (draw_info->dash_pattern[n] <= 0.0)
1526       break;
1527     length=scale*(draw_info->dash_pattern[n]+(n == 0 ? -0.5 : 0.5));
1528     if (offset > length)
1529       {
1530         offset-=length;
1531         n++;
1532         length=scale*(draw_info->dash_pattern[n]+(n == 0 ? -0.5 : 0.5));
1533         continue;
1534       }
1535     if (offset < length)
1536       {
1537         length-=offset;
1538         offset=0.0;
1539         break;
1540       }
1541     offset=0.0;
1542     n++;
1543   }
1544   status=MagickTrue;
1545   maximum_length=0.0;
1546   total_length=0.0;
1547   for (i=1; i < (ssize_t) number_vertices; i++)
1548   {
1549     dx=primitive_info[i].point.x-primitive_info[i-1].point.x;
1550     dy=primitive_info[i].point.y-primitive_info[i-1].point.y;
1551     maximum_length=hypot((double) dx,dy);
1552     if (length == 0.0)
1553       {
1554         n++;
1555         if (draw_info->dash_pattern[n] == 0.0)
1556           n=0;
1557         length=scale*(draw_info->dash_pattern[n]+(n == 0 ? -0.5 : 0.5));
1558       }
1559     for (total_length=0.0; (total_length+length) < maximum_length; )
1560     {
1561       total_length+=length;
1562       if ((n & 0x01) != 0)
1563         {
1564           dash_polygon[0]=primitive_info[0];
1565           dash_polygon[0].point.x=(double) (primitive_info[i-1].point.x+dx*
1566             total_length/maximum_length);
1567           dash_polygon[0].point.y=(double) (primitive_info[i-1].point.y+dy*
1568             total_length/maximum_length);
1569           j=1;
1570         }
1571       else
1572         {
1573           if ((j+1) > (ssize_t) (2*number_vertices))
1574             break;
1575           dash_polygon[j]=primitive_info[i-1];
1576           dash_polygon[j].point.x=(double) (primitive_info[i-1].point.x+dx*
1577             total_length/maximum_length);
1578           dash_polygon[j].point.y=(double) (primitive_info[i-1].point.y+dy*
1579             total_length/maximum_length);
1580           dash_polygon[j].coordinates=1;
1581           j++;
1582           dash_polygon[0].coordinates=(size_t) j;
1583           dash_polygon[j].primitive=UndefinedPrimitive;
1584           status|=DrawStrokePolygon(image,clone_info,dash_polygon,exception);
1585         }
1586       n++;
1587       if (draw_info->dash_pattern[n] == 0.0)
1588         n=0;
1589       length=scale*(draw_info->dash_pattern[n]+(n == 0 ? -0.5 : 0.5));
1590     }
1591     length-=(maximum_length-total_length);
1592     if ((n & 0x01) != 0)
1593       continue;
1594     dash_polygon[j]=primitive_info[i];
1595     dash_polygon[j].coordinates=1;
1596     j++;
1597   }
1598   if ((total_length < maximum_length) && ((n & 0x01) == 0) && (j > 1))
1599     {
1600       dash_polygon[j]=primitive_info[i-1];
1601       dash_polygon[j].point.x+=MagickEpsilon;
1602       dash_polygon[j].point.y+=MagickEpsilon;
1603       dash_polygon[j].coordinates=1;
1604       j++;
1605       dash_polygon[0].coordinates=(size_t) j;
1606       dash_polygon[j].primitive=UndefinedPrimitive;
1607       status|=DrawStrokePolygon(image,clone_info,dash_polygon,exception);
1608     }
1609   dash_polygon=(PrimitiveInfo *) RelinquishMagickMemory(dash_polygon);
1610   clone_info=DestroyDrawInfo(clone_info);
1611   if (image->debug != MagickFalse)
1612     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end draw-dash");
1613   return(status != 0 ? MagickTrue : MagickFalse);
1614 }
1615 \f
1616 /*
1617 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1618 %                                                                             %
1619 %                                                                             %
1620 %                                                                             %
1621 %   D r a w I m a g e                                                         %
1622 %                                                                             %
1623 %                                                                             %
1624 %                                                                             %
1625 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1626 %
1627 %  DrawImage() draws a graphic primitive on your image.  The primitive
1628 %  may be represented as a string or filename.  Precede the filename with an
1629 %  "at" sign (@) and the contents of the file are drawn on the image.  You
1630 %  can affect how text is drawn by setting one or more members of the draw
1631 %  info structure.
1632 %
1633 %  The format of the DrawImage method is:
1634 %
1635 %      MagickBooleanType DrawImage(Image *image,const DrawInfo *draw_info,
1636 %        ExceptionInfo *exception)
1637 %
1638 %  A description of each parameter follows:
1639 %
1640 %    o image: the image.
1641 %
1642 %    o draw_info: the draw info.
1643 %
1644 %    o exception: return any errors or warnings in this structure.
1645 %
1646 */
1647
1648 static inline MagickBooleanType IsPoint(const char *point)
1649 {
1650   char
1651     *p;
1652
1653   double
1654     value;
1655
1656   value=StringToDouble(point,&p);
1657   return((value == 0.0) && (p == point) ? MagickFalse : MagickTrue);
1658 }
1659
1660 static inline void TracePoint(PrimitiveInfo *primitive_info,
1661   const PointInfo point)
1662 {
1663   primitive_info->coordinates=1;
1664   primitive_info->point=point;
1665 }
1666
1667 MagickExport MagickBooleanType DrawImage(Image *image,const DrawInfo *draw_info,
1668   ExceptionInfo *exception)
1669 {
1670 #define RenderImageTag  "Render/Image"
1671
1672   AffineMatrix
1673     affine,
1674     current;
1675
1676   char
1677     key[2*MaxTextExtent],
1678     keyword[MaxTextExtent],
1679     geometry[MaxTextExtent],
1680     name[MaxTextExtent],
1681     pattern[MaxTextExtent],
1682     *primitive,
1683     *token;
1684
1685   const char
1686     *q;
1687
1688   DrawInfo
1689     **graphic_context;
1690
1691   MagickBooleanType
1692     proceed,
1693     status;
1694
1695   MagickRealType
1696     angle,
1697     factor,
1698     primitive_extent;
1699
1700   PointInfo
1701     point;
1702
1703   PixelInfo
1704     start_color;
1705
1706   PrimitiveInfo
1707     *primitive_info;
1708
1709   PrimitiveType
1710     primitive_type;
1711
1712   register const char
1713     *p;
1714
1715   register ssize_t
1716     i,
1717     x;
1718
1719   SegmentInfo
1720     bounds;
1721
1722   size_t
1723     length,
1724     number_points;
1725
1726   ssize_t
1727     j,
1728     k,
1729     n;
1730
1731   /*
1732     Ensure the annotation info is valid.
1733   */
1734   assert(image != (Image *) NULL);
1735   assert(image->signature == MagickSignature);
1736   if (image->debug != MagickFalse)
1737     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1738   assert(draw_info != (DrawInfo *) NULL);
1739   assert(draw_info->signature == MagickSignature);
1740   if (image->debug != MagickFalse)
1741     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
1742   if ((draw_info->primitive == (char *) NULL) ||
1743       (*draw_info->primitive == '\0'))
1744     return(MagickFalse);
1745   if (image->debug != MagickFalse)
1746     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"begin draw-image");
1747   if (*draw_info->primitive != '@')
1748     primitive=AcquireString(draw_info->primitive);
1749   else
1750     primitive=FileToString(draw_info->primitive+1,~0,exception);
1751   if (primitive == (char *) NULL)
1752     return(MagickFalse);
1753   primitive_extent=(MagickRealType) strlen(primitive);
1754   (void) SetImageArtifact(image,"MVG",primitive);
1755   n=0;
1756   /*
1757     Allocate primitive info memory.
1758   */
1759   graphic_context=(DrawInfo **) AcquireMagickMemory(
1760     sizeof(*graphic_context));
1761   if (graphic_context == (DrawInfo **) NULL)
1762     {
1763       primitive=DestroyString(primitive);
1764       ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
1765         image->filename);
1766     }
1767   number_points=2047;
1768   primitive_info=(PrimitiveInfo *) AcquireQuantumMemory((size_t) number_points,
1769     sizeof(*primitive_info));
1770   if (primitive_info == (PrimitiveInfo *) NULL)
1771     {
1772       primitive=DestroyString(primitive);
1773       for ( ; n >= 0; n--)
1774         graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
1775       graphic_context=(DrawInfo **) RelinquishMagickMemory(graphic_context);
1776       ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
1777         image->filename);
1778     }
1779   graphic_context[n]=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1780   graphic_context[n]->viewbox=image->page;
1781   if ((image->page.width == 0) || (image->page.height == 0))
1782     {
1783       graphic_context[n]->viewbox.width=image->columns;
1784       graphic_context[n]->viewbox.height=image->rows;
1785     }
1786   token=AcquireString(primitive);
1787   (void) QueryColorCompliance("#000000",AllCompliance,&start_color,
1788     exception);
1789   if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
1790     return(MagickFalse);
1791   status=MagickTrue;
1792   for (q=primitive; *q != '\0'; )
1793   {
1794     /*
1795       Interpret graphic primitive.
1796     */
1797     GetMagickToken(q,&q,keyword);
1798     if (*keyword == '\0')
1799       break;
1800     if (*keyword == '#')
1801       {
1802         /*
1803           Comment.
1804         */
1805         while ((*q != '\n') && (*q != '\0'))
1806           q++;
1807         continue;
1808       }
1809     p=q-strlen(keyword)-1;
1810     primitive_type=UndefinedPrimitive;
1811     current=graphic_context[n]->affine;
1812     GetAffineMatrix(&affine);
1813     switch (*keyword)
1814     {
1815       case ';':
1816         break;
1817       case 'a':
1818       case 'A':
1819       {
1820         if (LocaleCompare("affine",keyword) == 0)
1821           {
1822             GetMagickToken(q,&q,token);
1823             affine.sx=StringToDouble(token,(char **) NULL);
1824             GetMagickToken(q,&q,token);
1825             if (*token == ',')
1826               GetMagickToken(q,&q,token);
1827             affine.rx=StringToDouble(token,(char **) NULL);
1828             GetMagickToken(q,&q,token);
1829             if (*token == ',')
1830               GetMagickToken(q,&q,token);
1831             affine.ry=StringToDouble(token,(char **) NULL);
1832             GetMagickToken(q,&q,token);
1833             if (*token == ',')
1834               GetMagickToken(q,&q,token);
1835             affine.sy=StringToDouble(token,(char **) NULL);
1836             GetMagickToken(q,&q,token);
1837             if (*token == ',')
1838               GetMagickToken(q,&q,token);
1839             affine.tx=StringToDouble(token,(char **) NULL);
1840             GetMagickToken(q,&q,token);
1841             if (*token == ',')
1842               GetMagickToken(q,&q,token);
1843             affine.ty=StringToDouble(token,(char **) NULL);
1844             break;
1845           }
1846         if (LocaleCompare("arc",keyword) == 0)
1847           {
1848             primitive_type=ArcPrimitive;
1849             break;
1850           }
1851         status=MagickFalse;
1852         break;
1853       }
1854       case 'b':
1855       case 'B':
1856       {
1857         if (LocaleCompare("bezier",keyword) == 0)
1858           {
1859             primitive_type=BezierPrimitive;
1860             break;
1861           }
1862         if (LocaleCompare("border-color",keyword) == 0)
1863           {
1864             GetMagickToken(q,&q,token);
1865             (void) QueryColorCompliance(token,AllCompliance,
1866               &graphic_context[n]->border_color,exception);
1867             break;
1868           }
1869         status=MagickFalse;
1870         break;
1871       }
1872       case 'c':
1873       case 'C':
1874       {
1875         if (LocaleCompare("clip-path",keyword) == 0)
1876           {
1877             /*
1878               Create clip mask.
1879             */
1880             GetMagickToken(q,&q,token);
1881             (void) CloneString(&graphic_context[n]->clip_mask,token);
1882             (void) DrawClipPath(image,graphic_context[n],
1883               graphic_context[n]->clip_mask,exception);
1884             break;
1885           }
1886         if (LocaleCompare("clip-rule",keyword) == 0)
1887           {
1888             ssize_t
1889               fill_rule;
1890
1891             GetMagickToken(q,&q,token);
1892             fill_rule=ParseCommandOption(MagickFillRuleOptions,MagickFalse,
1893               token);
1894             if (fill_rule == -1)
1895               status=MagickFalse;
1896             else
1897               graphic_context[n]->fill_rule=(FillRule) fill_rule;
1898             break;
1899           }
1900         if (LocaleCompare("clip-units",keyword) == 0)
1901           {
1902             ssize_t
1903               clip_units;
1904
1905             GetMagickToken(q,&q,token);
1906             clip_units=ParseCommandOption(MagickClipPathOptions,MagickFalse,
1907               token);
1908             if (clip_units == -1)
1909               {
1910                 status=MagickFalse;
1911                 break;
1912               }
1913             graphic_context[n]->clip_units=(ClipPathUnits) clip_units;
1914             if (clip_units == ObjectBoundingBox)
1915               {
1916                 GetAffineMatrix(&current);
1917                 affine.sx=draw_info->bounds.x2;
1918                 affine.sy=draw_info->bounds.y2;
1919                 affine.tx=draw_info->bounds.x1;
1920                 affine.ty=draw_info->bounds.y1;
1921                 break;
1922               }
1923             break;
1924           }
1925         if (LocaleCompare("circle",keyword) == 0)
1926           {
1927             primitive_type=CirclePrimitive;
1928             break;
1929           }
1930         if (LocaleCompare("color",keyword) == 0)
1931           {
1932             primitive_type=ColorPrimitive;
1933             break;
1934           }
1935         status=MagickFalse;
1936         break;
1937       }
1938       case 'd':
1939       case 'D':
1940       {
1941         if (LocaleCompare("decorate",keyword) == 0)
1942           {
1943             ssize_t
1944               decorate;
1945
1946             GetMagickToken(q,&q,token);
1947             decorate=ParseCommandOption(MagickDecorateOptions,MagickFalse,
1948               token);
1949             if (decorate == -1)
1950               status=MagickFalse;
1951             else
1952               graphic_context[n]->decorate=(DecorationType) decorate;
1953             break;
1954           }
1955         status=MagickFalse;
1956         break;
1957       }
1958       case 'e':
1959       case 'E':
1960       {
1961         if (LocaleCompare("ellipse",keyword) == 0)
1962           {
1963             primitive_type=EllipsePrimitive;
1964             break;
1965           }
1966         if (LocaleCompare("encoding",keyword) == 0)
1967           {
1968             GetMagickToken(q,&q,token);
1969             (void) CloneString(&graphic_context[n]->encoding,token);
1970             break;
1971           }
1972         status=MagickFalse;
1973         break;
1974       }
1975       case 'f':
1976       case 'F':
1977       {
1978         if (LocaleCompare("fill",keyword) == 0)
1979           {
1980             GetMagickToken(q,&q,token);
1981             (void) FormatLocaleString(pattern,MaxTextExtent,"%s",token);
1982             if (GetImageArtifact(image,pattern) != (const char *) NULL)
1983               (void) DrawPatternPath(image,draw_info,token,
1984                 &graphic_context[n]->fill_pattern,exception);
1985             else
1986               {
1987                 status=QueryColorCompliance(token,AllCompliance,
1988                   &graphic_context[n]->fill,exception);
1989                 if (status == MagickFalse)
1990                   {
1991                     ImageInfo
1992                       *pattern_info;
1993
1994                     pattern_info=AcquireImageInfo();
1995                     (void) CopyMagickString(pattern_info->filename,token,
1996                       MaxTextExtent);
1997                     graphic_context[n]->fill_pattern=ReadImage(pattern_info,
1998                       exception);
1999                     CatchException(exception);
2000                     pattern_info=DestroyImageInfo(pattern_info);
2001                   }
2002               }
2003             break;
2004           }
2005         if (LocaleCompare("fill-opacity",keyword) == 0)
2006           {
2007             GetMagickToken(q,&q,token);
2008             factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
2009             graphic_context[n]->fill.alpha=(MagickRealType) QuantumRange*
2010               factor*StringToDouble(token,(char **) NULL);
2011             break;
2012           }
2013         if (LocaleCompare("fill-rule",keyword) == 0)
2014           {
2015             ssize_t
2016               fill_rule;
2017
2018             GetMagickToken(q,&q,token);
2019             fill_rule=ParseCommandOption(MagickFillRuleOptions,MagickFalse,
2020               token);
2021             if (fill_rule == -1)
2022               status=MagickFalse;
2023             else
2024               graphic_context[n]->fill_rule=(FillRule) fill_rule;
2025             break;
2026           }
2027         if (LocaleCompare("font",keyword) == 0)
2028           {
2029             GetMagickToken(q,&q,token);
2030             (void) CloneString(&graphic_context[n]->font,token);
2031             if (LocaleCompare("none",token) == 0)
2032               graphic_context[n]->font=(char *)
2033                 RelinquishMagickMemory(graphic_context[n]->font);
2034             break;
2035           }
2036         if (LocaleCompare("font-family",keyword) == 0)
2037           {
2038             GetMagickToken(q,&q,token);
2039             (void) CloneString(&graphic_context[n]->family,token);
2040             break;
2041           }
2042         if (LocaleCompare("font-size",keyword) == 0)
2043           {
2044             GetMagickToken(q,&q,token);
2045             graphic_context[n]->pointsize=StringToDouble(token,(char **) NULL);
2046             break;
2047           }
2048         if (LocaleCompare("font-stretch",keyword) == 0)
2049           {
2050             ssize_t
2051               stretch;
2052
2053             GetMagickToken(q,&q,token);
2054             stretch=ParseCommandOption(MagickStretchOptions,MagickFalse,token);
2055             if (stretch == -1)
2056               status=MagickFalse;
2057             else
2058               graphic_context[n]->stretch=(StretchType) stretch;
2059             break;
2060           }
2061         if (LocaleCompare("font-style",keyword) == 0)
2062           {
2063             ssize_t
2064               style;
2065
2066             GetMagickToken(q,&q,token);
2067             style=ParseCommandOption(MagickStyleOptions,MagickFalse,token);
2068             if (style == -1)
2069               status=MagickFalse;
2070             else
2071               graphic_context[n]->style=(StyleType) style;
2072             break;
2073           }
2074         if (LocaleCompare("font-weight",keyword) == 0)
2075           {
2076             GetMagickToken(q,&q,token);
2077             graphic_context[n]->weight=StringToUnsignedLong(token);
2078             if (LocaleCompare(token,"all") == 0)
2079               graphic_context[n]->weight=0;
2080             if (LocaleCompare(token,"bold") == 0)
2081               graphic_context[n]->weight=700;
2082             if (LocaleCompare(token,"bolder") == 0)
2083               if (graphic_context[n]->weight <= 800)
2084                 graphic_context[n]->weight+=100;
2085             if (LocaleCompare(token,"lighter") == 0)
2086               if (graphic_context[n]->weight >= 100)
2087                 graphic_context[n]->weight-=100;
2088             if (LocaleCompare(token,"normal") == 0)
2089               graphic_context[n]->weight=400;
2090             break;
2091           }
2092         status=MagickFalse;
2093         break;
2094       }
2095       case 'g':
2096       case 'G':
2097       {
2098         if (LocaleCompare("gradient-units",keyword) == 0)
2099           {
2100             GetMagickToken(q,&q,token);
2101             break;
2102           }
2103         if (LocaleCompare("gravity",keyword) == 0)
2104           {
2105             ssize_t
2106               gravity;
2107
2108             GetMagickToken(q,&q,token);
2109             gravity=ParseCommandOption(MagickGravityOptions,MagickFalse,token);
2110             if (gravity == -1)
2111               status=MagickFalse;
2112             else
2113               graphic_context[n]->gravity=(GravityType) gravity;
2114             break;
2115           }
2116         status=MagickFalse;
2117         break;
2118       }
2119       case 'i':
2120       case 'I':
2121       {
2122         if (LocaleCompare("image",keyword) == 0)
2123           {
2124             ssize_t
2125               compose;
2126
2127             primitive_type=ImagePrimitive;
2128             GetMagickToken(q,&q,token);
2129             compose=ParseCommandOption(MagickComposeOptions,MagickFalse,token);
2130             if (compose == -1)
2131               status=MagickFalse;
2132             else
2133               graphic_context[n]->compose=(CompositeOperator) compose;
2134             break;
2135           }
2136         if (LocaleCompare("interline-spacing",keyword) == 0)
2137           {
2138             GetMagickToken(q,&q,token);
2139             graphic_context[n]->interline_spacing=StringToDouble(token,
2140               (char **) NULL);
2141             break;
2142           }
2143         if (LocaleCompare("interword-spacing",keyword) == 0)
2144           {
2145             GetMagickToken(q,&q,token);
2146             graphic_context[n]->interword_spacing=StringToDouble(token,
2147               (char **) NULL);
2148             break;
2149           }
2150         status=MagickFalse;
2151         break;
2152       }
2153       case 'k':
2154       case 'K':
2155       {
2156         if (LocaleCompare("kerning",keyword) == 0)
2157           {
2158             GetMagickToken(q,&q,token);
2159             graphic_context[n]->kerning=StringToDouble(token,(char **) NULL);
2160             break;
2161           }
2162         status=MagickFalse;
2163         break;
2164       }
2165       case 'l':
2166       case 'L':
2167       {
2168         if (LocaleCompare("line",keyword) == 0)
2169           primitive_type=LinePrimitive;
2170         else
2171           status=MagickFalse;
2172         break;
2173       }
2174       case 'm':
2175       case 'M':
2176       {
2177         if (LocaleCompare("matte",keyword) == 0)
2178           primitive_type=MattePrimitive;
2179         else
2180           status=MagickFalse;
2181         break;
2182       }
2183       case 'o':
2184       case 'O':
2185       {
2186         if (LocaleCompare("offset",keyword) == 0)
2187           {
2188             GetMagickToken(q,&q,token);
2189             break;
2190           }
2191         if (LocaleCompare("opacity",keyword) == 0)
2192           {
2193             GetMagickToken(q,&q,token);
2194             factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
2195             graphic_context[n]->alpha=ClampToQuantum((MagickRealType)
2196               QuantumRange*(1.0-((1.0-QuantumScale*graphic_context[n]->alpha)*
2197               factor*StringToDouble(token,(char **) NULL))));
2198             graphic_context[n]->fill.alpha=(double) graphic_context[n]->alpha;
2199             graphic_context[n]->stroke.alpha=(double) graphic_context[n]->alpha;
2200             break;
2201           }
2202         status=MagickFalse;
2203         break;
2204       }
2205       case 'p':
2206       case 'P':
2207       {
2208         if (LocaleCompare("path",keyword) == 0)
2209           {
2210             primitive_type=PathPrimitive;
2211             break;
2212           }
2213         if (LocaleCompare("point",keyword) == 0)
2214           {
2215             primitive_type=PointPrimitive;
2216             break;
2217           }
2218         if (LocaleCompare("polyline",keyword) == 0)
2219           {
2220             primitive_type=PolylinePrimitive;
2221             break;
2222           }
2223         if (LocaleCompare("polygon",keyword) == 0)
2224           {
2225             primitive_type=PolygonPrimitive;
2226             break;
2227           }
2228         if (LocaleCompare("pop",keyword) == 0)
2229           {
2230             GetMagickToken(q,&q,token);
2231             if (LocaleCompare("clip-path",token) == 0)
2232               break;
2233             if (LocaleCompare("defs",token) == 0)
2234               break;
2235             if (LocaleCompare("gradient",token) == 0)
2236               break;
2237             if (LocaleCompare("graphic-context",token) == 0)
2238               {
2239                 if (n <= 0)
2240                   {
2241                     (void) ThrowMagickException(exception,GetMagickModule(),
2242                       DrawError,"UnbalancedGraphicContextPushPop","`%s'",token);
2243                     n=0;
2244                     break;
2245                   }
2246                 if (graphic_context[n]->clip_mask != (char *) NULL)
2247                   if (LocaleCompare(graphic_context[n]->clip_mask,
2248                       graphic_context[n-1]->clip_mask) != 0)
2249                     (void) SetImageMask(image,(Image *) NULL,exception);
2250                 graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
2251                 n--;
2252                 break;
2253               }
2254             if (LocaleCompare("pattern",token) == 0)
2255               break;
2256             status=MagickFalse;
2257             break;
2258           }
2259         if (LocaleCompare("push",keyword) == 0)
2260           {
2261             GetMagickToken(q,&q,token);
2262             if (LocaleCompare("clip-path",token) == 0)
2263               {
2264                 char
2265                   name[MaxTextExtent];
2266
2267                 GetMagickToken(q,&q,token);
2268                 (void) FormatLocaleString(name,MaxTextExtent,"%s",token);
2269                 for (p=q; *q != '\0'; )
2270                 {
2271                   GetMagickToken(q,&q,token);
2272                   if (LocaleCompare(token,"pop") != 0)
2273                     continue;
2274                   GetMagickToken(q,(const char **) NULL,token);
2275                   if (LocaleCompare(token,"clip-path") != 0)
2276                     continue;
2277                   break;
2278                 }
2279                 (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
2280                 (void) SetImageArtifact(image,name,token);
2281                 GetMagickToken(q,&q,token);
2282                 break;
2283               }
2284             if (LocaleCompare("gradient",token) == 0)
2285               {
2286                 char
2287                   key[2*MaxTextExtent],
2288                   name[MaxTextExtent],
2289                   type[MaxTextExtent];
2290
2291                 SegmentInfo
2292                   segment;
2293
2294                 GetMagickToken(q,&q,token);
2295                 (void) CopyMagickString(name,token,MaxTextExtent);
2296                 GetMagickToken(q,&q,token);
2297                 (void) CopyMagickString(type,token,MaxTextExtent);
2298                 GetMagickToken(q,&q,token);
2299                 segment.x1=StringToDouble(token,(char **) NULL);
2300                 GetMagickToken(q,&q,token);
2301                 if (*token == ',')
2302                   GetMagickToken(q,&q,token);
2303                 segment.y1=StringToDouble(token,(char **) NULL);
2304                 GetMagickToken(q,&q,token);
2305                 if (*token == ',')
2306                   GetMagickToken(q,&q,token);
2307                 segment.x2=StringToDouble(token,(char **) NULL);
2308                 GetMagickToken(q,&q,token);
2309                 if (*token == ',')
2310                   GetMagickToken(q,&q,token);
2311                 segment.y2=StringToDouble(token,(char **) NULL);
2312                 if (LocaleCompare(type,"radial") == 0)
2313                   {
2314                     GetMagickToken(q,&q,token);
2315                     if (*token == ',')
2316                       GetMagickToken(q,&q,token);
2317                   }
2318                 for (p=q; *q != '\0'; )
2319                 {
2320                   GetMagickToken(q,&q,token);
2321                   if (LocaleCompare(token,"pop") != 0)
2322                     continue;
2323                   GetMagickToken(q,(const char **) NULL,token);
2324                   if (LocaleCompare(token,"gradient") != 0)
2325                     continue;
2326                   break;
2327                 }
2328                 (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
2329                 bounds.x1=graphic_context[n]->affine.sx*segment.x1+
2330                   graphic_context[n]->affine.ry*segment.y1+
2331                   graphic_context[n]->affine.tx;
2332                 bounds.y1=graphic_context[n]->affine.rx*segment.x1+
2333                   graphic_context[n]->affine.sy*segment.y1+
2334                   graphic_context[n]->affine.ty;
2335                 bounds.x2=graphic_context[n]->affine.sx*segment.x2+
2336                   graphic_context[n]->affine.ry*segment.y2+
2337                   graphic_context[n]->affine.tx;
2338                 bounds.y2=graphic_context[n]->affine.rx*segment.x2+
2339                   graphic_context[n]->affine.sy*segment.y2+
2340                   graphic_context[n]->affine.ty;
2341                 (void) FormatLocaleString(key,MaxTextExtent,"%s",name);
2342                 (void) SetImageArtifact(image,key,token);
2343                 (void) FormatLocaleString(key,MaxTextExtent,"%s-geometry",name);
2344                 (void) FormatLocaleString(geometry,MaxTextExtent,
2345                   "%gx%g%+.15g%+.15g",
2346                   MagickMax(fabs(bounds.x2-bounds.x1+1.0),1.0),
2347                   MagickMax(fabs(bounds.y2-bounds.y1+1.0),1.0),
2348                   bounds.x1,bounds.y1);
2349                 (void) SetImageArtifact(image,key,geometry);
2350                 GetMagickToken(q,&q,token);
2351                 break;
2352               }
2353             if (LocaleCompare("pattern",token) == 0)
2354               {
2355                 RectangleInfo
2356                   bounds;
2357
2358                 GetMagickToken(q,&q,token);
2359                 (void) CopyMagickString(name,token,MaxTextExtent);
2360                 GetMagickToken(q,&q,token);
2361                 bounds.x=(ssize_t) ceil(StringToDouble(token,(char **) NULL)-
2362                   0.5);
2363                 GetMagickToken(q,&q,token);
2364                 if (*token == ',')
2365                   GetMagickToken(q,&q,token);
2366                 bounds.y=(ssize_t) ceil(StringToDouble(token,(char **) NULL)-
2367                   0.5);
2368                 GetMagickToken(q,&q,token);
2369                 if (*token == ',')
2370                   GetMagickToken(q,&q,token);
2371                 bounds.width=(size_t) floor(StringToDouble(token,
2372                   (char **) NULL)+0.5);
2373                 GetMagickToken(q,&q,token);
2374                 if (*token == ',')
2375                   GetMagickToken(q,&q,token);
2376                 bounds.height=(size_t) floor(StringToDouble(token,
2377                   (char **) NULL)+0.5);
2378                 for (p=q; *q != '\0'; )
2379                 {
2380                   GetMagickToken(q,&q,token);
2381                   if (LocaleCompare(token,"pop") != 0)
2382                     continue;
2383                   GetMagickToken(q,(const char **) NULL,token);
2384                   if (LocaleCompare(token,"pattern") != 0)
2385                     continue;
2386                   break;
2387                 }
2388                 (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
2389                 (void) FormatLocaleString(key,MaxTextExtent,"%s",name);
2390                 (void) SetImageArtifact(image,key,token);
2391                 (void) FormatLocaleString(key,MaxTextExtent,"%s-geometry",name);
2392                 (void) FormatLocaleString(geometry,MaxTextExtent,
2393                   "%.20gx%.20g%+.20g%+.20g",(double) bounds.width,(double)
2394                   bounds.height,(double) bounds.x,(double) bounds.y);
2395                 (void) SetImageArtifact(image,key,geometry);
2396                 GetMagickToken(q,&q,token);
2397                 break;
2398               }
2399             if (LocaleCompare("graphic-context",token) == 0)
2400               {
2401                 n++;
2402                 graphic_context=(DrawInfo **) ResizeQuantumMemory(
2403                   graphic_context,(size_t) (n+1),sizeof(*graphic_context));
2404                 if (graphic_context == (DrawInfo **) NULL)
2405                   {
2406                     (void) ThrowMagickException(exception,GetMagickModule(),
2407                       ResourceLimitError,"MemoryAllocationFailed","`%s'",
2408                       image->filename);
2409                     break;
2410                   }
2411                 graphic_context[n]=CloneDrawInfo((ImageInfo *) NULL,
2412                   graphic_context[n-1]);
2413                 break;
2414               }
2415             if (LocaleCompare("defs",token) == 0)
2416               break;
2417             status=MagickFalse;
2418             break;
2419           }
2420         status=MagickFalse;
2421         break;
2422       }
2423       case 'r':
2424       case 'R':
2425       {
2426         if (LocaleCompare("rectangle",keyword) == 0)
2427           {
2428             primitive_type=RectanglePrimitive;
2429             break;
2430           }
2431         if (LocaleCompare("rotate",keyword) == 0)
2432           {
2433             GetMagickToken(q,&q,token);
2434             angle=StringToDouble(token,(char **) NULL);
2435             affine.sx=cos(DegreesToRadians(fmod((double) angle,360.0)));
2436             affine.rx=sin(DegreesToRadians(fmod((double) angle,360.0)));
2437             affine.ry=(-sin(DegreesToRadians(fmod((double) angle,360.0))));
2438             affine.sy=cos(DegreesToRadians(fmod((double) angle,360.0)));
2439             break;
2440           }
2441         if (LocaleCompare("roundRectangle",keyword) == 0)
2442           {
2443             primitive_type=RoundRectanglePrimitive;
2444             break;
2445           }
2446         status=MagickFalse;
2447         break;
2448       }
2449       case 's':
2450       case 'S':
2451       {
2452         if (LocaleCompare("scale",keyword) == 0)
2453           {
2454             GetMagickToken(q,&q,token);
2455             affine.sx=StringToDouble(token,(char **) NULL);
2456             GetMagickToken(q,&q,token);
2457             if (*token == ',')
2458               GetMagickToken(q,&q,token);
2459             affine.sy=StringToDouble(token,(char **) NULL);
2460             break;
2461           }
2462         if (LocaleCompare("skewX",keyword) == 0)
2463           {
2464             GetMagickToken(q,&q,token);
2465             angle=StringToDouble(token,(char **) NULL);
2466             affine.ry=sin(DegreesToRadians(angle));
2467             break;
2468           }
2469         if (LocaleCompare("skewY",keyword) == 0)
2470           {
2471             GetMagickToken(q,&q,token);
2472             angle=StringToDouble(token,(char **) NULL);
2473             affine.rx=(-tan(DegreesToRadians(angle)/2.0));
2474             break;
2475           }
2476         if (LocaleCompare("stop-color",keyword) == 0)
2477           {
2478             PixelInfo
2479               stop_color;
2480
2481             GetMagickToken(q,&q,token);
2482             (void) QueryColorCompliance(token,AllCompliance,&stop_color,
2483               exception);
2484             (void) GradientImage(image,LinearGradient,ReflectSpread,
2485               &start_color,&stop_color,exception);
2486             start_color=stop_color;
2487             GetMagickToken(q,&q,token);
2488             break;
2489           }
2490         if (LocaleCompare("stroke",keyword) == 0)
2491           {
2492             GetMagickToken(q,&q,token);
2493             (void) FormatLocaleString(pattern,MaxTextExtent,"%s",token);
2494             if (GetImageArtifact(image,pattern) != (const char *) NULL)
2495               (void) DrawPatternPath(image,draw_info,token,
2496                 &graphic_context[n]->stroke_pattern,exception);
2497             else
2498               {
2499                 status=QueryColorCompliance(token,AllCompliance,
2500                   &graphic_context[n]->stroke,exception);
2501                 if (status == MagickFalse)
2502                   {
2503                     ImageInfo
2504                       *pattern_info;
2505
2506                     pattern_info=AcquireImageInfo();
2507                     (void) CopyMagickString(pattern_info->filename,token,
2508                       MaxTextExtent);
2509                     graphic_context[n]->stroke_pattern=ReadImage(pattern_info,
2510                       exception);
2511                     CatchException(exception);
2512                     pattern_info=DestroyImageInfo(pattern_info);
2513                   }
2514               }
2515             break;
2516           }
2517         if (LocaleCompare("stroke-antialias",keyword) == 0)
2518           {
2519             GetMagickToken(q,&q,token);
2520             graphic_context[n]->stroke_antialias=
2521               StringToLong(token) != 0 ? MagickTrue : MagickFalse;
2522             break;
2523           }
2524         if (LocaleCompare("stroke-dasharray",keyword) == 0)
2525           {
2526             if (graphic_context[n]->dash_pattern != (double *) NULL)
2527               graphic_context[n]->dash_pattern=(double *)
2528                 RelinquishMagickMemory(graphic_context[n]->dash_pattern);
2529             if (IsPoint(q) != MagickFalse)
2530               {
2531                 const char
2532                   *p;
2533
2534                 p=q;
2535                 GetMagickToken(p,&p,token);
2536                 if (*token == ',')
2537                   GetMagickToken(p,&p,token);
2538                 for (x=0; IsPoint(token) != MagickFalse; x++)
2539                 {
2540                   GetMagickToken(p,&p,token);
2541                   if (*token == ',')
2542                     GetMagickToken(p,&p,token);
2543                 }
2544                 graphic_context[n]->dash_pattern=(double *)
2545                   AcquireQuantumMemory((size_t) (2UL*x+1UL),
2546                   sizeof(*graphic_context[n]->dash_pattern));
2547                 if (graphic_context[n]->dash_pattern == (double *) NULL)
2548                   {
2549                     (void) ThrowMagickException(exception,GetMagickModule(),
2550                       ResourceLimitError,"MemoryAllocationFailed","`%s'",
2551                       image->filename);
2552                     break;
2553                   }
2554                 for (j=0; j < x; j++)
2555                 {
2556                   GetMagickToken(q,&q,token);
2557                   if (*token == ',')
2558                     GetMagickToken(q,&q,token);
2559                   graphic_context[n]->dash_pattern[j]=StringToDouble(token,
2560                     (char **) NULL);
2561                 }
2562                 if ((x & 0x01) != 0)
2563                   for ( ; j < (2*x); j++)
2564                     graphic_context[n]->dash_pattern[j]=
2565                       graphic_context[n]->dash_pattern[j-x];
2566                 graphic_context[n]->dash_pattern[j]=0.0;
2567                 break;
2568               }
2569             GetMagickToken(q,&q,token);
2570             break;
2571           }
2572         if (LocaleCompare("stroke-dashoffset",keyword) == 0)
2573           {
2574             GetMagickToken(q,&q,token);
2575             graphic_context[n]->dash_offset=StringToDouble(token,
2576               (char **) NULL);
2577             break;
2578           }
2579         if (LocaleCompare("stroke-linecap",keyword) == 0)
2580           {
2581             ssize_t
2582               linecap;
2583
2584             GetMagickToken(q,&q,token);
2585             linecap=ParseCommandOption(MagickLineCapOptions,MagickFalse,token);
2586             if (linecap == -1)
2587               status=MagickFalse;
2588             else
2589               graphic_context[n]->linecap=(LineCap) linecap;
2590             break;
2591           }
2592         if (LocaleCompare("stroke-linejoin",keyword) == 0)
2593           {
2594             ssize_t
2595               linejoin;
2596
2597             GetMagickToken(q,&q,token);
2598             linejoin=ParseCommandOption(MagickLineJoinOptions,MagickFalse,
2599               token);
2600             if (linejoin == -1)
2601               status=MagickFalse;
2602             else
2603               graphic_context[n]->linejoin=(LineJoin) linejoin;
2604             break;
2605           }
2606         if (LocaleCompare("stroke-miterlimit",keyword) == 0)
2607           {
2608             GetMagickToken(q,&q,token);
2609             graphic_context[n]->miterlimit=StringToUnsignedLong(token);
2610             break;
2611           }
2612         if (LocaleCompare("stroke-opacity",keyword) == 0)
2613           {
2614             GetMagickToken(q,&q,token);
2615             factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
2616             graphic_context[n]->stroke.alpha=(MagickRealType) QuantumRange*
2617               factor*StringToDouble(token,(char **) NULL);
2618             break;
2619           }
2620         if (LocaleCompare("stroke-width",keyword) == 0)
2621           {
2622             GetMagickToken(q,&q,token);
2623             graphic_context[n]->stroke_width=StringToDouble(token,
2624               (char **) NULL);
2625             break;
2626           }
2627         status=MagickFalse;
2628         break;
2629       }
2630       case 't':
2631       case 'T':
2632       {
2633         if (LocaleCompare("text",keyword) == 0)
2634           {
2635             primitive_type=TextPrimitive;
2636             break;
2637           }
2638         if (LocaleCompare("text-align",keyword) == 0)
2639           {
2640             ssize_t
2641               align;
2642
2643             GetMagickToken(q,&q,token);
2644             align=ParseCommandOption(MagickAlignOptions,MagickFalse,token);
2645             if (align == -1)
2646               status=MagickFalse;
2647             else
2648               graphic_context[n]->align=(AlignType) align;
2649             break;
2650           }
2651         if (LocaleCompare("text-anchor",keyword) == 0)
2652           {
2653             ssize_t
2654               align;
2655
2656             GetMagickToken(q,&q,token);
2657             align=ParseCommandOption(MagickAlignOptions,MagickFalse,token);
2658             if (align == -1)
2659               status=MagickFalse;
2660             else
2661               graphic_context[n]->align=(AlignType) align;
2662             break;
2663           }
2664         if (LocaleCompare("text-antialias",keyword) == 0)
2665           {
2666             GetMagickToken(q,&q,token);
2667             graphic_context[n]->text_antialias=
2668               StringToLong(token) != 0 ? MagickTrue : MagickFalse;
2669             break;
2670           }
2671         if (LocaleCompare("text-undercolor",keyword) == 0)
2672           {
2673             GetMagickToken(q,&q,token);
2674             (void) QueryColorCompliance(token,AllCompliance,
2675               &graphic_context[n]->undercolor,exception);
2676             break;
2677           }
2678         if (LocaleCompare("translate",keyword) == 0)
2679           {
2680             GetMagickToken(q,&q,token);
2681             affine.tx=StringToDouble(token,(char **) NULL);
2682             GetMagickToken(q,&q,token);
2683             if (*token == ',')
2684               GetMagickToken(q,&q,token);
2685             affine.ty=StringToDouble(token,(char **) NULL);
2686             break;
2687           }
2688         status=MagickFalse;
2689         break;
2690       }
2691       case 'v':
2692       case 'V':
2693       {
2694         if (LocaleCompare("viewbox",keyword) == 0)
2695           {
2696             GetMagickToken(q,&q,token);
2697             graphic_context[n]->viewbox.x=(ssize_t) ceil(StringToDouble(token,
2698               (char **) NULL)-0.5);
2699             GetMagickToken(q,&q,token);
2700             if (*token == ',')
2701               GetMagickToken(q,&q,token);
2702             graphic_context[n]->viewbox.y=(ssize_t) ceil(StringToDouble(token,
2703               (char **) NULL)-0.5);
2704             GetMagickToken(q,&q,token);
2705             if (*token == ',')
2706               GetMagickToken(q,&q,token);
2707             graphic_context[n]->viewbox.width=(size_t) floor(StringToDouble(
2708               token,(char **) NULL)+0.5);
2709             GetMagickToken(q,&q,token);
2710             if (*token == ',')
2711               GetMagickToken(q,&q,token);
2712             graphic_context[n]->viewbox.height=(size_t) floor(StringToDouble(
2713               token,(char **) NULL)+0.5);
2714             break;
2715           }
2716         status=MagickFalse;
2717         break;
2718       }
2719       default:
2720       {
2721         status=MagickFalse;
2722         break;
2723       }
2724     }
2725     if (status == MagickFalse)
2726       break;
2727     if ((affine.sx != 1.0) || (affine.rx != 0.0) || (affine.ry != 0.0) ||
2728         (affine.sy != 1.0) || (affine.tx != 0.0) || (affine.ty != 0.0))
2729       {
2730         graphic_context[n]->affine.sx=current.sx*affine.sx+current.ry*affine.rx;
2731         graphic_context[n]->affine.rx=current.rx*affine.sx+current.sy*affine.rx;
2732         graphic_context[n]->affine.ry=current.sx*affine.ry+current.ry*affine.sy;
2733         graphic_context[n]->affine.sy=current.rx*affine.ry+current.sy*affine.sy;
2734         graphic_context[n]->affine.tx=current.sx*affine.tx+current.ry*affine.ty+
2735           current.tx;
2736         graphic_context[n]->affine.ty=current.rx*affine.tx+current.sy*affine.ty+
2737           current.ty;
2738       }
2739     if (primitive_type == UndefinedPrimitive)
2740       {
2741         if (image->debug != MagickFalse)
2742           (void) LogMagickEvent(DrawEvent,GetMagickModule(),"  %.*s",
2743             (int) (q-p),p);
2744         continue;
2745       }
2746     /*
2747       Parse the primitive attributes.
2748     */
2749     i=0;
2750     j=0;
2751     primitive_info[0].point.x=0.0;
2752     primitive_info[0].point.y=0.0;
2753     for (x=0; *q != '\0'; x++)
2754     {
2755       /*
2756         Define points.
2757       */
2758       if (IsPoint(q) == MagickFalse)
2759         break;
2760       GetMagickToken(q,&q,token);
2761       point.x=StringToDouble(token,(char **) NULL);
2762       GetMagickToken(q,&q,token);
2763       if (*token == ',')
2764         GetMagickToken(q,&q,token);
2765       point.y=StringToDouble(token,(char **) NULL);
2766       GetMagickToken(q,(const char **) NULL,token);
2767       if (*token == ',')
2768         GetMagickToken(q,&q,token);
2769       primitive_info[i].primitive=primitive_type;
2770       primitive_info[i].point=point;
2771       primitive_info[i].coordinates=0;
2772       primitive_info[i].method=FloodfillMethod;
2773       i++;
2774       if (i < (ssize_t) number_points)
2775         continue;
2776       number_points<<=1;
2777       primitive_info=(PrimitiveInfo *) ResizeQuantumMemory(primitive_info,
2778         (size_t) number_points,sizeof(*primitive_info));
2779       if (primitive_info == (PrimitiveInfo *) NULL)
2780         {
2781           (void) ThrowMagickException(exception,GetMagickModule(),
2782             ResourceLimitError,"MemoryAllocationFailed","`%s'",image->filename);
2783           break;
2784         }
2785     }
2786     primitive_info[j].primitive=primitive_type;
2787     primitive_info[j].coordinates=(size_t) x;
2788     primitive_info[j].method=FloodfillMethod;
2789     primitive_info[j].text=(char *) NULL;
2790     /*
2791       Circumscribe primitive within a circle.
2792     */
2793     bounds.x1=primitive_info[j].point.x;
2794     bounds.y1=primitive_info[j].point.y;
2795     bounds.x2=primitive_info[j].point.x;
2796     bounds.y2=primitive_info[j].point.y;
2797     for (k=1; k < (ssize_t) primitive_info[j].coordinates; k++)
2798     {
2799       point=primitive_info[j+k].point;
2800       if (point.x < bounds.x1)
2801         bounds.x1=point.x;
2802       if (point.y < bounds.y1)
2803         bounds.y1=point.y;
2804       if (point.x > bounds.x2)
2805         bounds.x2=point.x;
2806       if (point.y > bounds.y2)
2807         bounds.y2=point.y;
2808     }
2809     /*
2810       Speculate how many points our primitive might consume.
2811     */
2812     length=primitive_info[j].coordinates;
2813     switch (primitive_type)
2814     {
2815       case RectanglePrimitive:
2816       {
2817         length*=5;
2818         break;
2819       }
2820       case RoundRectanglePrimitive:
2821       {
2822         length*=5+8*BezierQuantum;
2823         break;
2824       }
2825       case BezierPrimitive:
2826       {
2827         if (primitive_info[j].coordinates > 107)
2828           (void) ThrowMagickException(exception,GetMagickModule(),DrawError,
2829             "TooManyBezierCoordinates","`%s'",token);
2830         length=BezierQuantum*primitive_info[j].coordinates;
2831         break;
2832       }
2833       case PathPrimitive:
2834       {
2835         char
2836           *s,
2837           *t;
2838
2839         GetMagickToken(q,&q,token);
2840         length=1;
2841         t=token;
2842         for (s=token; *s != '\0'; s=t)
2843         {
2844           double
2845             value;
2846
2847           value=StringToDouble(s,&t);
2848           (void) value;
2849           if (s == t)
2850             {
2851               t++;
2852               continue;
2853             }
2854           length++;
2855         }
2856         length=length*BezierQuantum/2;
2857         break;
2858       }
2859       case CirclePrimitive:
2860       case ArcPrimitive:
2861       case EllipsePrimitive:
2862       {
2863         MagickRealType
2864           alpha,
2865           beta,
2866           radius;
2867
2868         alpha=bounds.x2-bounds.x1;
2869         beta=bounds.y2-bounds.y1;
2870         radius=hypot((double) alpha,(double) beta);
2871         length=2*((size_t) ceil((double) MagickPI*radius))+6*BezierQuantum+360;
2872         break;
2873       }
2874       default:
2875         break;
2876     }
2877     if ((size_t) (i+length) >= number_points)
2878       {
2879         /*
2880           Resize based on speculative points required by primitive.
2881         */
2882         number_points+=length+1;
2883         primitive_info=(PrimitiveInfo *) ResizeQuantumMemory(primitive_info,
2884           (size_t) number_points,sizeof(*primitive_info));
2885         if (primitive_info == (PrimitiveInfo *) NULL)
2886           {
2887             (void) ThrowMagickException(exception,GetMagickModule(),
2888               ResourceLimitError,"MemoryAllocationFailed","`%s'",
2889               image->filename);
2890             break;
2891           }
2892       }
2893     switch (primitive_type)
2894     {
2895       case PointPrimitive:
2896       default:
2897       {
2898         if (primitive_info[j].coordinates != 1)
2899           {
2900             status=MagickFalse;
2901             break;
2902           }
2903         TracePoint(primitive_info+j,primitive_info[j].point);
2904         i=(ssize_t) (j+primitive_info[j].coordinates);
2905         break;
2906       }
2907       case LinePrimitive:
2908       {
2909         if (primitive_info[j].coordinates != 2)
2910           {
2911             status=MagickFalse;
2912             break;
2913           }
2914         TraceLine(primitive_info+j,primitive_info[j].point,
2915           primitive_info[j+1].point);
2916         i=(ssize_t) (j+primitive_info[j].coordinates);
2917         break;
2918       }
2919       case RectanglePrimitive:
2920       {
2921         if (primitive_info[j].coordinates != 2)
2922           {
2923             status=MagickFalse;
2924             break;
2925           }
2926         TraceRectangle(primitive_info+j,primitive_info[j].point,
2927           primitive_info[j+1].point);
2928         i=(ssize_t) (j+primitive_info[j].coordinates);
2929         break;
2930       }
2931       case RoundRectanglePrimitive:
2932       {
2933         if (primitive_info[j].coordinates != 3)
2934           {
2935             status=MagickFalse;
2936             break;
2937           }
2938         TraceRoundRectangle(primitive_info+j,primitive_info[j].point,
2939           primitive_info[j+1].point,primitive_info[j+2].point);
2940         i=(ssize_t) (j+primitive_info[j].coordinates);
2941         break;
2942       }
2943       case ArcPrimitive:
2944       {
2945         if (primitive_info[j].coordinates != 3)
2946           {
2947             primitive_type=UndefinedPrimitive;
2948             break;
2949           }
2950         TraceArc(primitive_info+j,primitive_info[j].point,
2951           primitive_info[j+1].point,primitive_info[j+2].point);
2952         i=(ssize_t) (j+primitive_info[j].coordinates);
2953         break;
2954       }
2955       case EllipsePrimitive:
2956       {
2957         if (primitive_info[j].coordinates != 3)
2958           {
2959             status=MagickFalse;
2960             break;
2961           }
2962         TraceEllipse(primitive_info+j,primitive_info[j].point,
2963           primitive_info[j+1].point,primitive_info[j+2].point);
2964         i=(ssize_t) (j+primitive_info[j].coordinates);
2965         break;
2966       }
2967       case CirclePrimitive:
2968       {
2969         if (primitive_info[j].coordinates != 2)
2970           {
2971             status=MagickFalse;
2972             break;
2973           }
2974         TraceCircle(primitive_info+j,primitive_info[j].point,
2975           primitive_info[j+1].point);
2976         i=(ssize_t) (j+primitive_info[j].coordinates);
2977         break;
2978       }
2979       case PolylinePrimitive:
2980         break;
2981       case PolygonPrimitive:
2982       {
2983         primitive_info[i]=primitive_info[j];
2984         primitive_info[i].coordinates=0;
2985         primitive_info[j].coordinates++;
2986         i++;
2987         break;
2988       }
2989       case BezierPrimitive:
2990       {
2991         if (primitive_info[j].coordinates < 3)
2992           {
2993             status=MagickFalse;
2994             break;
2995           }
2996         TraceBezier(primitive_info+j,primitive_info[j].coordinates);
2997         i=(ssize_t) (j+primitive_info[j].coordinates);
2998         break;
2999       }
3000       case PathPrimitive:
3001       {
3002         i=(ssize_t) (j+TracePath(primitive_info+j,token));
3003         break;
3004       }
3005       case ColorPrimitive:
3006       case MattePrimitive:
3007       {
3008         ssize_t
3009           method;
3010
3011         if (primitive_info[j].coordinates != 1)
3012           {
3013             status=MagickFalse;
3014             break;
3015           }
3016         GetMagickToken(q,&q,token);
3017         method=ParseCommandOption(MagickMethodOptions,MagickFalse,token);
3018         if (method == -1)
3019           status=MagickFalse;
3020         else
3021           primitive_info[j].method=(PaintMethod) method;
3022         break;
3023       }
3024       case TextPrimitive:
3025       {
3026         if (primitive_info[j].coordinates != 1)
3027           {
3028             status=MagickFalse;
3029             break;
3030           }
3031         if (*token != ',')
3032           GetMagickToken(q,&q,token);
3033         primitive_info[j].text=AcquireString(token);
3034         break;
3035       }
3036       case ImagePrimitive:
3037       {
3038         if (primitive_info[j].coordinates != 2)
3039           {
3040             status=MagickFalse;
3041             break;
3042           }
3043         GetMagickToken(q,&q,token);
3044         primitive_info[j].text=AcquireString(token);
3045         break;
3046       }
3047     }
3048     if (primitive_info == (PrimitiveInfo *) NULL)
3049       break;
3050     if (image->debug != MagickFalse)
3051       (void) LogMagickEvent(DrawEvent,GetMagickModule(),"  %.*s",(int) (q-p),p);
3052     if (status == MagickFalse)
3053       break;
3054     primitive_info[i].primitive=UndefinedPrimitive;
3055     if (i == 0)
3056       continue;
3057     /*
3058       Transform points.
3059     */
3060     for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
3061     {
3062       point=primitive_info[i].point;
3063       primitive_info[i].point.x=graphic_context[n]->affine.sx*point.x+
3064         graphic_context[n]->affine.ry*point.y+graphic_context[n]->affine.tx;
3065       primitive_info[i].point.y=graphic_context[n]->affine.rx*point.x+
3066         graphic_context[n]->affine.sy*point.y+graphic_context[n]->affine.ty;
3067       point=primitive_info[i].point;
3068       if (point.x < graphic_context[n]->bounds.x1)
3069         graphic_context[n]->bounds.x1=point.x;
3070       if (point.y < graphic_context[n]->bounds.y1)
3071         graphic_context[n]->bounds.y1=point.y;
3072       if (point.x > graphic_context[n]->bounds.x2)
3073         graphic_context[n]->bounds.x2=point.x;
3074       if (point.y > graphic_context[n]->bounds.y2)
3075         graphic_context[n]->bounds.y2=point.y;
3076       if (primitive_info[i].primitive == ImagePrimitive)
3077         break;
3078       if (i >= (ssize_t) number_points)
3079         ThrowFatalException(ResourceLimitFatalError,"MemoryAllocationFailed");
3080     }
3081     if (graphic_context[n]->render != MagickFalse)
3082       {
3083         if ((n != 0) && (graphic_context[n]->clip_mask != (char *) NULL) &&
3084             (LocaleCompare(graphic_context[n]->clip_mask,
3085              graphic_context[n-1]->clip_mask) != 0))
3086           (void) DrawClipPath(image,graphic_context[n],
3087             graphic_context[n]->clip_mask,exception);
3088         (void) DrawPrimitive(image,graphic_context[n],primitive_info,exception);
3089       }
3090     if (primitive_info->text != (char *) NULL)
3091       primitive_info->text=(char *) RelinquishMagickMemory(
3092         primitive_info->text);
3093     proceed=SetImageProgress(image,RenderImageTag,q-primitive,(MagickSizeType)
3094       primitive_extent);
3095     if (proceed == MagickFalse)
3096       break;
3097   }
3098   if (image->debug != MagickFalse)
3099     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end draw-image");
3100   /*
3101     Relinquish resources.
3102   */
3103   token=DestroyString(token);
3104   if (primitive_info != (PrimitiveInfo *) NULL)
3105     primitive_info=(PrimitiveInfo *) RelinquishMagickMemory(primitive_info);
3106   primitive=DestroyString(primitive);
3107   for ( ; n >= 0; n--)
3108     graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
3109   graphic_context=(DrawInfo **) RelinquishMagickMemory(graphic_context);
3110   if (status == MagickFalse)
3111     ThrowBinaryException(DrawError,"NonconformingDrawingPrimitiveDefinition",
3112       keyword);
3113   return(status);
3114 }
3115 \f
3116 /*
3117 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3118 %                                                                             %
3119 %                                                                             %
3120 %                                                                             %
3121 %     D r a w G r a d i e n t I m a g e                                       %
3122 %                                                                             %
3123 %                                                                             %
3124 %                                                                             %
3125 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3126 %
3127 %  DrawGradientImage() draws a linear gradient on the image.
3128 %
3129 %  The format of the DrawGradientImage method is:
3130 %
3131 %      MagickBooleanType DrawGradientImage(Image *image,
3132 %        const DrawInfo *draw_info,ExceptionInfo *exception)
3133 %
3134 %  A description of each parameter follows:
3135 %
3136 %    o image: the image.
3137 %
3138 %    o draw_info: the draw info.
3139 %
3140 %    o exception: return any errors or warnings in this structure.
3141 %
3142 */
3143
3144 static inline MagickRealType GetStopColorOffset(const GradientInfo *gradient,
3145   const ssize_t x,const ssize_t y)
3146 {
3147   switch (gradient->type)
3148   {
3149     case UndefinedGradient:
3150     case LinearGradient:
3151     {
3152       MagickRealType
3153         gamma,
3154         length,
3155         offset,
3156         scale;
3157
3158       PointInfo
3159         p,
3160         q;
3161
3162       const SegmentInfo
3163         *gradient_vector;
3164
3165       gradient_vector=(&gradient->gradient_vector);
3166       p.x=gradient_vector->x2-gradient_vector->x1;
3167       p.y=gradient_vector->y2-gradient_vector->y1;
3168       q.x=(double) x-gradient_vector->x1;
3169       q.y=(double) y-gradient_vector->y1;
3170       length=sqrt(q.x*q.x+q.y*q.y);
3171       gamma=sqrt(p.x*p.x+p.y*p.y)*length;
3172       gamma=1.0/(gamma <= MagickEpsilon ? 1.0 : gamma);
3173       scale=p.x*q.x+p.y*q.y;
3174       offset=gamma*scale*length;
3175       return(offset);
3176     }
3177     case RadialGradient:
3178     {
3179       MagickRealType
3180         length,
3181         offset;
3182
3183       PointInfo
3184         v;
3185
3186       v.x=(double) x-gradient->center.x;
3187       v.y=(double) y-gradient->center.y;
3188       length=sqrt(v.x*v.x+v.y*v.y);
3189       if (gradient->spread == RepeatSpread)
3190         return(length);
3191       offset=length/gradient->radius;
3192       return(offset);
3193     }
3194   }
3195   return(0.0);
3196 }
3197
3198 MagickExport MagickBooleanType DrawGradientImage(Image *image,
3199   const DrawInfo *draw_info,ExceptionInfo *exception)
3200 {
3201   CacheView
3202     *image_view;
3203
3204   const GradientInfo
3205     *gradient;
3206
3207   const SegmentInfo
3208     *gradient_vector;
3209
3210   MagickBooleanType
3211     status;
3212
3213   PixelInfo
3214     zero;
3215
3216   MagickRealType
3217     length;
3218
3219   PointInfo
3220     point;
3221
3222   RectangleInfo
3223     bounding_box;
3224
3225   ssize_t
3226     y;
3227
3228   /*
3229     Draw linear or radial gradient on image.
3230   */
3231   assert(image != (Image *) NULL);
3232   assert(image->signature == MagickSignature);
3233   if (image->debug != MagickFalse)
3234     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
3235   assert(draw_info != (const DrawInfo *) NULL);
3236   gradient=(&draw_info->gradient);
3237   gradient_vector=(&gradient->gradient_vector);
3238   point.x=gradient_vector->x2-gradient_vector->x1;
3239   point.y=gradient_vector->y2-gradient_vector->y1;
3240   length=sqrt(point.x*point.x+point.y*point.y);
3241   bounding_box=gradient->bounding_box;
3242   status=MagickTrue;
3243   GetPixelInfo(image,&zero);
3244   image_view=AcquireCacheView(image);
3245 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3246   #pragma omp parallel for schedule(static) shared(status)
3247 #endif
3248   for (y=bounding_box.y; y < (ssize_t) bounding_box.height; y++)
3249   {
3250     PixelInfo
3251       composite,
3252       pixel;
3253
3254     MagickRealType
3255       alpha,
3256       offset;
3257
3258     register Quantum
3259       *restrict q;
3260
3261     register ssize_t
3262       i,
3263       x;
3264
3265     ssize_t
3266       j;
3267
3268     if (status == MagickFalse)
3269       continue;
3270     q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
3271     if (q == (Quantum *) NULL)
3272       {
3273         status=MagickFalse;
3274         continue;
3275       }
3276     pixel=zero;
3277     composite=zero;
3278     offset=GetStopColorOffset(gradient,0,y);
3279     if (gradient->type != RadialGradient)
3280       offset/=length;
3281     for (x=bounding_box.x; x < (ssize_t) bounding_box.width; x++)
3282     {
3283       GetPixelInfoPixel(image,q,&pixel);
3284       switch (gradient->spread)
3285       {
3286         case UndefinedSpread:
3287         case PadSpread:
3288         {
3289           if ((x != (ssize_t) ceil(gradient_vector->x1-0.5)) ||
3290               (y != (ssize_t) ceil(gradient_vector->y1-0.5)))
3291             {
3292               offset=GetStopColorOffset(gradient,x,y);
3293               if (gradient->type != RadialGradient)
3294                 offset/=length;
3295             }
3296           for (i=0; i < (ssize_t) gradient->number_stops; i++)
3297             if (offset < gradient->stops[i].offset)
3298               break;
3299           if ((offset < 0.0) || (i == 0))
3300             composite=gradient->stops[0].color;
3301           else
3302             if ((offset > 1.0) || (i == (ssize_t) gradient->number_stops))
3303               composite=gradient->stops[gradient->number_stops-1].color;
3304             else
3305               {
3306                 j=i;
3307                 i--;
3308                 alpha=(offset-gradient->stops[i].offset)/
3309                   (gradient->stops[j].offset-gradient->stops[i].offset);
3310                 CompositePixelInfoBlend(&gradient->stops[i].color,1.0-alpha,
3311                   &gradient->stops[j].color,alpha,&composite);
3312               }
3313           break;
3314         }
3315         case ReflectSpread:
3316         {
3317           if ((x != (ssize_t) ceil(gradient_vector->x1-0.5)) ||
3318               (y != (ssize_t) ceil(gradient_vector->y1-0.5)))
3319             {
3320               offset=GetStopColorOffset(gradient,x,y);
3321               if (gradient->type != RadialGradient)
3322                 offset/=length;
3323             }
3324           if (offset < 0.0)
3325             offset=(-offset);
3326           if ((ssize_t) fmod(offset,2.0) == 0)
3327             offset=fmod(offset,1.0);
3328           else
3329             offset=1.0-fmod(offset,1.0);
3330           for (i=0; i < (ssize_t) gradient->number_stops; i++)
3331             if (offset < gradient->stops[i].offset)
3332               break;
3333           if (i == 0)
3334             composite=gradient->stops[0].color;
3335           else
3336             if (i == (ssize_t) gradient->number_stops)
3337               composite=gradient->stops[gradient->number_stops-1].color;
3338             else
3339               {
3340                 j=i;
3341                 i--;
3342                 alpha=(offset-gradient->stops[i].offset)/
3343                   (gradient->stops[j].offset-gradient->stops[i].offset);
3344                 CompositePixelInfoBlend(&gradient->stops[i].color,1.0-alpha,
3345                   &gradient->stops[j].color,alpha,&composite);
3346               }
3347           break;
3348         }
3349         case RepeatSpread:
3350         {
3351           MagickBooleanType
3352             antialias;
3353
3354           MagickRealType
3355             repeat;
3356
3357           antialias=MagickFalse;
3358           repeat=0.0;
3359           if ((x != (ssize_t) ceil(gradient_vector->x1-0.5)) ||
3360               (y != (ssize_t) ceil(gradient_vector->y1-0.5)))
3361             {
3362               offset=GetStopColorOffset(gradient,x,y);
3363               if (gradient->type == LinearGradient)
3364                 {
3365                   repeat=fmod(offset,length);
3366                   if (repeat < 0.0)
3367                     repeat=length-fmod(-repeat,length);
3368                   else
3369                     repeat=fmod(offset,length);
3370                   antialias=(repeat < length) && ((repeat+1.0) > length) ?
3371                     MagickTrue : MagickFalse;
3372                   offset=repeat/length;
3373                 }
3374               else
3375                 {
3376                   repeat=fmod(offset,gradient->radius);
3377                   if (repeat < 0.0)
3378                     repeat=gradient->radius-fmod(-repeat,gradient->radius);
3379                   else
3380                     repeat=fmod(offset,gradient->radius);
3381                   antialias=repeat+1.0 > gradient->radius ? MagickTrue :
3382                     MagickFalse;
3383                   offset=repeat/gradient->radius;
3384                 }
3385             }
3386           for (i=0; i < (ssize_t) gradient->number_stops; i++)
3387             if (offset < gradient->stops[i].offset)
3388               break;
3389           if (i == 0)
3390             composite=gradient->stops[0].color;
3391           else
3392             if (i == (ssize_t) gradient->number_stops)
3393               composite=gradient->stops[gradient->number_stops-1].color;
3394             else
3395               {
3396                 j=i;
3397                 i--;
3398                 alpha=(offset-gradient->stops[i].offset)/
3399                   (gradient->stops[j].offset-gradient->stops[i].offset);
3400                 if (antialias != MagickFalse)
3401                   {
3402                     if (gradient->type == LinearGradient)
3403                       alpha=length-repeat;
3404                     else
3405                       alpha=gradient->radius-repeat;
3406                     i=0;
3407                     j=(ssize_t) gradient->number_stops-1L;
3408                   }
3409                 CompositePixelInfoBlend(&gradient->stops[i].color,1.0-alpha,
3410                   &gradient->stops[j].color,alpha,&composite);
3411               }
3412           break;
3413         }
3414       }
3415       CompositePixelInfoOver(&composite,composite.alpha,&pixel,pixel.alpha,
3416         &pixel);
3417       SetPixelInfoPixel(image,&pixel,q);
3418       q+=GetPixelChannels(image);
3419     }
3420     if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
3421       status=MagickFalse;
3422   }
3423   image_view=DestroyCacheView(image_view);
3424   return(status);
3425 }
3426 \f
3427 /*
3428 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3429 %                                                                             %
3430 %                                                                             %
3431 %                                                                             %
3432 %   D r a w P a t t e r n P a t h                                             %
3433 %                                                                             %
3434 %                                                                             %
3435 %                                                                             %
3436 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3437 %
3438 %  DrawPatternPath() draws a pattern.
3439 %
3440 %  The format of the DrawPatternPath method is:
3441 %
3442 %      MagickBooleanType DrawPatternPath(Image *image,const DrawInfo *draw_info,
3443 %        const char *name,Image **pattern,ExceptionInfo *exception)
3444 %
3445 %  A description of each parameter follows:
3446 %
3447 %    o image: the image.
3448 %
3449 %    o draw_info: the draw info.
3450 %
3451 %    o name: the pattern name.
3452 %
3453 %    o image: the image.
3454 %
3455 %    o exception: return any errors or warnings in this structure.
3456 %
3457 */
3458 MagickExport MagickBooleanType DrawPatternPath(Image *image,
3459   const DrawInfo *draw_info,const char *name,Image **pattern,
3460   ExceptionInfo *exception)
3461 {
3462   char
3463     property[MaxTextExtent];
3464
3465   const char
3466     *geometry,
3467     *path;
3468
3469   DrawInfo
3470     *clone_info;
3471
3472   ImageInfo
3473     *image_info;
3474
3475   MagickBooleanType
3476     status;
3477
3478   assert(image != (Image *) NULL);
3479   assert(image->signature == MagickSignature);
3480   if (image->debug != MagickFalse)
3481     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
3482   assert(draw_info != (const DrawInfo *) NULL);
3483   assert(name != (const char *) NULL);
3484   (void) FormatLocaleString(property,MaxTextExtent,"%s",name);
3485   path=GetImageArtifact(image,property);
3486   if (path == (const char *) NULL)
3487     return(MagickFalse);
3488   (void) FormatLocaleString(property,MaxTextExtent,"%s-geometry",name);
3489   geometry=GetImageArtifact(image,property);
3490   if (geometry == (const char *) NULL)
3491     return(MagickFalse);
3492   if ((*pattern) != (Image *) NULL)
3493     *pattern=DestroyImage(*pattern);
3494   image_info=AcquireImageInfo();
3495   image_info->size=AcquireString(geometry);
3496   *pattern=AcquireImage(image_info,exception);
3497   image_info=DestroyImageInfo(image_info);
3498   (void) QueryColorCompliance("#00000000",AllCompliance,
3499     &(*pattern)->background_color,exception);
3500   (void) SetImageBackgroundColor(*pattern,exception);
3501   if (image->debug != MagickFalse)
3502     (void) LogMagickEvent(DrawEvent,GetMagickModule(),
3503       "begin pattern-path %s %s",name,geometry);
3504   clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
3505   clone_info->fill_pattern=NewImageList();
3506   clone_info->stroke_pattern=NewImageList();
3507   (void) CloneString(&clone_info->primitive,path);
3508   status=DrawImage(*pattern,clone_info,exception);
3509   clone_info=DestroyDrawInfo(clone_info);
3510   if (image->debug != MagickFalse)
3511     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end pattern-path");
3512   return(status);
3513 }
3514 \f
3515 /*
3516 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3517 %                                                                             %
3518 %                                                                             %
3519 %                                                                             %
3520 +   D r a w P o l y g o n P r i m i t i v e                                   %
3521 %                                                                             %
3522 %                                                                             %
3523 %                                                                             %
3524 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3525 %
3526 %  DrawPolygonPrimitive() draws a polygon on the image.
3527 %
3528 %  The format of the DrawPolygonPrimitive method is:
3529 %
3530 %      MagickBooleanType DrawPolygonPrimitive(Image *image,
3531 %        const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
3532 %        ExceptionInfo *exception)
3533 %
3534 %  A description of each parameter follows:
3535 %
3536 %    o image: the image.
3537 %
3538 %    o draw_info: the draw info.
3539 %
3540 %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
3541 %
3542 %    o exception: return any errors or warnings in this structure.
3543 %
3544 */
3545
3546 static PolygonInfo **DestroyPolygonThreadSet(PolygonInfo **polygon_info)
3547 {
3548   register ssize_t
3549     i;
3550
3551   assert(polygon_info != (PolygonInfo **) NULL);
3552   for (i=0; i < (ssize_t) GetOpenMPMaximumThreads(); i++)
3553     if (polygon_info[i] != (PolygonInfo *) NULL)
3554       polygon_info[i]=DestroyPolygonInfo(polygon_info[i]);
3555   polygon_info=(PolygonInfo **) RelinquishMagickMemory(polygon_info);
3556   return(polygon_info);
3557 }
3558
3559 static PolygonInfo **AcquirePolygonThreadSet(const DrawInfo *draw_info,
3560   const PrimitiveInfo *primitive_info)
3561 {
3562   PathInfo
3563     *restrict path_info;
3564
3565   PolygonInfo
3566     **polygon_info;
3567
3568   register ssize_t
3569     i;
3570
3571   size_t
3572     number_threads;
3573
3574   number_threads=GetOpenMPMaximumThreads();
3575   polygon_info=(PolygonInfo **) AcquireQuantumMemory(number_threads,
3576     sizeof(*polygon_info));
3577   if (polygon_info == (PolygonInfo **) NULL)
3578     return((PolygonInfo **) NULL);
3579   (void) ResetMagickMemory(polygon_info,0,GetOpenMPMaximumThreads()*
3580     sizeof(*polygon_info));
3581   path_info=ConvertPrimitiveToPath(draw_info,primitive_info);
3582   if (path_info == (PathInfo *) NULL)
3583     return(DestroyPolygonThreadSet(polygon_info));
3584   for (i=0; i < (ssize_t) number_threads; i++)
3585   {
3586     polygon_info[i]=ConvertPathToPolygon(draw_info,path_info);
3587     if (polygon_info[i] == (PolygonInfo *) NULL)
3588       return(DestroyPolygonThreadSet(polygon_info));
3589   }
3590   path_info=(PathInfo *) RelinquishMagickMemory(path_info);
3591   return(polygon_info);
3592 }
3593
3594 static MagickRealType GetFillAlpha(PolygonInfo *polygon_info,
3595   const MagickRealType mid,const MagickBooleanType fill,
3596   const FillRule fill_rule,const double x,const double y,
3597   MagickRealType *stroke_alpha)
3598 {
3599   MagickRealType
3600     alpha,
3601     beta,
3602     distance,
3603     subpath_alpha;
3604
3605   PointInfo
3606     delta;
3607
3608   register EdgeInfo
3609     *p;
3610
3611   register const PointInfo
3612     *q;
3613
3614   register ssize_t
3615     i;
3616
3617   ssize_t
3618     j,
3619     winding_number;
3620
3621   /*
3622     Compute fill & stroke opacity for this (x,y) point.
3623   */
3624   *stroke_alpha=0.0;
3625   subpath_alpha=0.0;
3626   p=polygon_info->edges;
3627   for (j=0; j < (ssize_t) polygon_info->number_edges; j++, p++)
3628   {
3629     if (y <= (p->bounds.y1-mid-0.5))
3630       break;
3631     if (y > (p->bounds.y2+mid+0.5))
3632       {
3633         (void) DestroyEdge(polygon_info,(size_t) j);
3634         continue;
3635       }
3636     if ((x <= (p->bounds.x1-mid-0.5)) || (x > (p->bounds.x2+mid+0.5)))
3637       continue;
3638     i=(ssize_t) MagickMax((double) p->highwater,1.0);
3639     for ( ; i < (ssize_t) p->number_points; i++)
3640     {
3641       if (y <= (p->points[i-1].y-mid-0.5))
3642         break;
3643       if (y > (p->points[i].y+mid+0.5))
3644         continue;
3645       if (p->scanline != y)
3646         {
3647           p->scanline=y;
3648           p->highwater=(size_t) i;
3649         }
3650       /*
3651         Compute distance between a point and an edge.
3652       */
3653       q=p->points+i-1;
3654       delta.x=(q+1)->x-q->x;
3655       delta.y=(q+1)->y-q->y;
3656       beta=delta.x*(x-q->x)+delta.y*(y-q->y);
3657       if (beta < 0.0)
3658         {
3659           delta.x=x-q->x;
3660           delta.y=y-q->y;
3661           distance=delta.x*delta.x+delta.y*delta.y;
3662         }
3663       else
3664         {
3665           alpha=delta.x*delta.x+delta.y*delta.y;
3666           if (beta > alpha)
3667             {
3668               delta.x=x-(q+1)->x;
3669               delta.y=y-(q+1)->y;
3670               distance=delta.x*delta.x+delta.y*delta.y;
3671             }
3672           else
3673             {
3674               alpha=1.0/alpha;
3675               beta=delta.x*(y-q->y)-delta.y*(x-q->x);
3676               distance=alpha*beta*beta;
3677             }
3678         }
3679       /*
3680         Compute stroke & subpath opacity.
3681       */
3682       beta=0.0;
3683       if (p->ghostline == MagickFalse)
3684         {
3685           alpha=mid+0.5;
3686           if ((*stroke_alpha < 1.0) &&
3687               (distance <= ((alpha+0.25)*(alpha+0.25))))
3688             {
3689               alpha=mid-0.5;
3690               if (distance <= ((alpha+0.25)*(alpha+0.25)))
3691                 *stroke_alpha=1.0;
3692               else
3693                 {
3694                   beta=1.0;
3695                   if (distance != 1.0)
3696                     beta=sqrt((double) distance);
3697                   alpha=beta-mid-0.5;
3698                   if (*stroke_alpha < ((alpha-0.25)*(alpha-0.25)))
3699                     *stroke_alpha=(alpha-0.25)*(alpha-0.25);
3700                 }
3701             }
3702         }
3703       if ((fill == MagickFalse) || (distance > 1.0) || (subpath_alpha >= 1.0))
3704         continue;
3705       if (distance <= 0.0)
3706         {
3707           subpath_alpha=1.0;
3708           continue;
3709         }
3710       if (distance > 1.0)
3711         continue;
3712       if (beta == 0.0)
3713         {
3714           beta=1.0;
3715           if (distance != 1.0)
3716             beta=sqrt(distance);
3717         }
3718       alpha=beta-1.0;
3719       if (subpath_alpha < (alpha*alpha))
3720         subpath_alpha=alpha*alpha;
3721     }
3722   }
3723   /*
3724     Compute fill opacity.
3725   */
3726   if (fill == MagickFalse)
3727     return(0.0);
3728   if (subpath_alpha >= 1.0)
3729     return(1.0);
3730   /*
3731     Determine winding number.
3732   */
3733   winding_number=0;
3734   p=polygon_info->edges;
3735   for (j=0; j < (ssize_t) polygon_info->number_edges; j++, p++)
3736   {
3737     if (y <= p->bounds.y1)
3738       break;
3739     if ((y > p->bounds.y2) || (x <= p->bounds.x1))
3740       continue;
3741     if (x > p->bounds.x2)
3742       {
3743         winding_number+=p->direction ? 1 : -1;
3744         continue;
3745       }
3746     i=(ssize_t) MagickMax((double) p->highwater,1.0);
3747     for ( ; i < (ssize_t) p->number_points; i++)
3748       if (y <= p->points[i].y)
3749         break;
3750     q=p->points+i-1;
3751     if ((((q+1)->x-q->x)*(y-q->y)) <= (((q+1)->y-q->y)*(x-q->x)))
3752       winding_number+=p->direction ? 1 : -1;
3753   }
3754   if (fill_rule != NonZeroRule)
3755     {
3756       if ((MagickAbsoluteValue(winding_number) & 0x01) != 0)
3757         return(1.0);
3758     }
3759   else
3760     if (MagickAbsoluteValue(winding_number) != 0)
3761       return(1.0);
3762   return(subpath_alpha);
3763 }
3764
3765 static MagickBooleanType DrawPolygonPrimitive(Image *image,
3766   const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
3767   ExceptionInfo *exception)
3768 {
3769   CacheView
3770     *image_view;
3771
3772   MagickBooleanType
3773     fill,
3774     status;
3775
3776   MagickRealType
3777     mid;
3778
3779   PolygonInfo
3780     **restrict polygon_info;
3781
3782   register EdgeInfo
3783     *p;
3784
3785   register ssize_t
3786     i;
3787
3788   SegmentInfo
3789     bounds;
3790
3791   ssize_t
3792     start,
3793     stop,
3794     y;
3795
3796   /*
3797     Compute bounding box.
3798   */
3799   assert(image != (Image *) NULL);
3800   assert(image->signature == MagickSignature);
3801   if (image->debug != MagickFalse)
3802     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
3803   assert(draw_info != (DrawInfo *) NULL);
3804   assert(draw_info->signature == MagickSignature);
3805   assert(primitive_info != (PrimitiveInfo *) NULL);
3806   if (primitive_info->coordinates == 0)
3807     return(MagickTrue);
3808   polygon_info=AcquirePolygonThreadSet(draw_info,primitive_info);
3809   if (polygon_info == (PolygonInfo **) NULL)
3810     return(MagickFalse);
3811   if (0)
3812     DrawBoundingRectangles(image,draw_info,polygon_info[0],exception);
3813   if (image->debug != MagickFalse)
3814     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin draw-polygon");
3815   fill=(primitive_info->method == FillToBorderMethod) ||
3816     (primitive_info->method == FloodfillMethod) ? MagickTrue : MagickFalse;
3817   mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
3818   bounds=polygon_info[0]->edges[0].bounds;
3819   for (i=1; i < (ssize_t) polygon_info[0]->number_edges; i++)
3820   {
3821     p=polygon_info[0]->edges+i;
3822     if (p->bounds.x1 < bounds.x1)
3823       bounds.x1=p->bounds.x1;
3824     if (p->bounds.y1 < bounds.y1)
3825       bounds.y1=p->bounds.y1;
3826     if (p->bounds.x2 > bounds.x2)
3827       bounds.x2=p->bounds.x2;
3828     if (p->bounds.y2 > bounds.y2)
3829       bounds.y2=p->bounds.y2;
3830   }
3831   bounds.x1-=(mid+1.0);
3832   bounds.x1=bounds.x1 < 0.0 ? 0.0 : (size_t) ceil(bounds.x1-0.5) >=
3833     image->columns ? (double) image->columns-1.0 : bounds.x1;
3834   bounds.y1-=(mid+1.0);
3835   bounds.y1=bounds.y1 < 0.0 ? 0.0 : (size_t) ceil(bounds.y1-0.5) >=
3836     image->rows ? (double) image->rows-1.0 : bounds.y1;
3837   bounds.x2+=(mid+1.0);
3838   bounds.x2=bounds.x2 < 0.0 ? 0.0 : (size_t) floor(bounds.x2+0.5) >=
3839     image->columns ? (double) image->columns-1.0 : bounds.x2;
3840   bounds.y2+=(mid+1.0);
3841   bounds.y2=bounds.y2 < 0.0 ? 0.0 : (size_t) floor(bounds.y2+0.5) >=
3842     image->rows ? (double) image->rows-1.0 : bounds.y2;
3843   status=MagickTrue;
3844   image_view=AcquireCacheView(image);
3845   if (primitive_info->coordinates == 1)
3846     {
3847       /*
3848         Draw point.
3849       */
3850       start=(ssize_t) ceil(bounds.y1-0.5);
3851       stop=(ssize_t) floor(bounds.y2+0.5);
3852 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3853       #pragma omp parallel for schedule(static) shared(status)
3854 #endif
3855       for (y=start; y <= stop; y++)
3856       {
3857         MagickBooleanType
3858           sync;
3859
3860         PixelInfo
3861           pixel;
3862
3863         register ssize_t
3864           x;
3865
3866         register Quantum
3867           *restrict q;
3868
3869         ssize_t
3870           start,
3871           stop;
3872
3873         if (status == MagickFalse)
3874           continue;
3875         start=(ssize_t) ceil(bounds.x1-0.5);
3876         stop=(ssize_t) floor(bounds.x2+0.5);
3877         x=start;
3878         q=GetCacheViewAuthenticPixels(image_view,x,y,(size_t) (stop-x+1),
3879           1,exception);
3880         if (q == (Quantum *) NULL)
3881           {
3882             status=MagickFalse;
3883             continue;
3884           }
3885         GetPixelInfo(image,&pixel);
3886         for ( ; x <= stop; x++)
3887         {
3888           if ((x == (ssize_t) ceil(primitive_info->point.x-0.5)) &&
3889               (y == (ssize_t) ceil(primitive_info->point.y-0.5)))
3890             {
3891               (void) GetStrokeColor(draw_info,x,y,&pixel,exception);
3892               SetPixelInfoPixel(image,&pixel,q);
3893             }
3894           q+=GetPixelChannels(image);
3895         }
3896         sync=SyncCacheViewAuthenticPixels(image_view,exception);
3897         if (sync == MagickFalse)
3898           status=MagickFalse;
3899       }
3900       image_view=DestroyCacheView(image_view);
3901       polygon_info=DestroyPolygonThreadSet(polygon_info);
3902       if (image->debug != MagickFalse)
3903         (void) LogMagickEvent(DrawEvent,GetMagickModule(),
3904           "    end draw-polygon");
3905       return(status);
3906     }
3907   /*
3908     Draw polygon or line.
3909   */
3910   if (image->matte == MagickFalse)
3911     (void) SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
3912   start=(ssize_t) ceil(bounds.y1-0.5);
3913   stop=(ssize_t) floor(bounds.y2+0.5);
3914 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3915   #pragma omp parallel for schedule(static) shared(status)
3916 #endif
3917   for (y=start; y <= stop; y++)
3918   {
3919     const int
3920       id = GetOpenMPThreadId();
3921
3922     MagickRealType
3923       fill_alpha,
3924       stroke_alpha;
3925
3926     PixelInfo
3927       fill_color,
3928       stroke_color;
3929
3930     register Quantum
3931       *restrict q;
3932
3933     register ssize_t
3934       x;
3935
3936     ssize_t
3937       start,
3938       stop;
3939
3940     if (status == MagickFalse)
3941       continue;
3942     start=(ssize_t) ceil(bounds.x1-0.5);
3943     stop=(ssize_t) floor(bounds.x2+0.5);
3944     q=GetCacheViewAuthenticPixels(image_view,start,y,(size_t) (stop-start+1),1,
3945       exception);
3946     if (q == (Quantum *) NULL)
3947       {
3948         status=MagickFalse;
3949         continue;
3950       }
3951     for (x=start; x <= stop; x++)
3952     {
3953       /*
3954         Fill and/or stroke.
3955       */
3956       fill_alpha=GetFillAlpha(polygon_info[id],mid,fill,draw_info->fill_rule,
3957         (double) x,(double) y,&stroke_alpha);
3958       if (draw_info->stroke_antialias == MagickFalse)
3959         {
3960           fill_alpha=fill_alpha > 0.25 ? 1.0 : 0.0;
3961           stroke_alpha=stroke_alpha > 0.25 ? 1.0 : 0.0;
3962         }
3963       (void) GetFillColor(draw_info,x,y,&fill_color,exception);
3964       fill_alpha=fill_alpha*fill_color.alpha;
3965       CompositePixelOver(image,&fill_color,fill_alpha,q,(MagickRealType)
3966         GetPixelAlpha(image,q),q);
3967       (void) GetStrokeColor(draw_info,x,y,&stroke_color,exception);
3968       stroke_alpha=stroke_alpha*stroke_color.alpha;
3969       CompositePixelOver(image,&stroke_color,stroke_alpha,q,(MagickRealType)
3970         GetPixelAlpha(image,q),q);
3971       q+=GetPixelChannels(image);
3972     }
3973     if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
3974       status=MagickFalse;
3975   }
3976   image_view=DestroyCacheView(image_view);
3977   polygon_info=DestroyPolygonThreadSet(polygon_info);
3978   if (image->debug != MagickFalse)
3979     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end draw-polygon");
3980   return(status);
3981 }
3982 \f
3983 /*
3984 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3985 %                                                                             %
3986 %                                                                             %
3987 %                                                                             %
3988 %   D r a w P r i m i t i v e                                                 %
3989 %                                                                             %
3990 %                                                                             %
3991 %                                                                             %
3992 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3993 %
3994 %  DrawPrimitive() draws a primitive (line, rectangle, ellipse) on the image.
3995 %
3996 %  The format of the DrawPrimitive method is:
3997 %
3998 %      MagickBooleanType DrawPrimitive(Image *image,const DrawInfo *draw_info,
3999 %        PrimitiveInfo *primitive_info,ExceptionInfo *exception)
4000 %
4001 %  A description of each parameter follows:
4002 %
4003 %    o image: the image.
4004 %
4005 %    o draw_info: the draw info.
4006 %
4007 %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
4008 %
4009 %    o exception: return any errors or warnings in this structure.
4010 %
4011 */
4012
4013 static void LogPrimitiveInfo(const PrimitiveInfo *primitive_info)
4014 {
4015   const char
4016     *methods[] =
4017     {
4018       "point",
4019       "replace",
4020       "floodfill",
4021       "filltoborder",
4022       "reset",
4023       "?"
4024     };
4025
4026   PointInfo
4027     p,
4028     q,
4029     point;
4030
4031   register ssize_t
4032     i,
4033     x;
4034
4035   ssize_t
4036     coordinates,
4037     y;
4038
4039   x=(ssize_t) ceil(primitive_info->point.x-0.5);
4040   y=(ssize_t) ceil(primitive_info->point.y-0.5);
4041   switch (primitive_info->primitive)
4042   {
4043     case PointPrimitive:
4044     {
4045       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4046         "PointPrimitive %.20g,%.20g %s",(double) x,(double) y,
4047         methods[primitive_info->method]);
4048       return;
4049     }
4050     case ColorPrimitive:
4051     {
4052       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4053         "ColorPrimitive %.20g,%.20g %s",(double) x,(double) y,
4054         methods[primitive_info->method]);
4055       return;
4056     }
4057     case MattePrimitive:
4058     {
4059       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4060         "MattePrimitive %.20g,%.20g %s",(double) x,(double) y,
4061         methods[primitive_info->method]);
4062       return;
4063     }
4064     case TextPrimitive:
4065     {
4066       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4067         "TextPrimitive %.20g,%.20g",(double) x,(double) y);
4068       return;
4069     }
4070     case ImagePrimitive:
4071     {
4072       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4073         "ImagePrimitive %.20g,%.20g",(double) x,(double) y);
4074       return;
4075     }
4076     default:
4077       break;
4078   }
4079   coordinates=0;
4080   p=primitive_info[0].point;
4081   q.x=(-1.0);
4082   q.y=(-1.0);
4083   for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
4084   {
4085     point=primitive_info[i].point;
4086     if (coordinates <= 0)
4087       {
4088         coordinates=(ssize_t) primitive_info[i].coordinates;
4089         (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4090           "    begin open (%.20g)",(double) coordinates);
4091         p=point;
4092       }
4093     point=primitive_info[i].point;
4094     if ((fabs(q.x-point.x) > MagickEpsilon) ||
4095         (fabs(q.y-point.y) > MagickEpsilon))
4096       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4097         "      %.20g: %.18g,%.18g",(double) coordinates,point.x,point.y);
4098     else
4099       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4100         "      %.20g: %g %g (duplicate)",(double) coordinates,point.x,point.y);
4101     q=point;
4102     coordinates--;
4103     if (coordinates > 0)
4104       continue;
4105     if ((fabs(p.x-point.x) > MagickEpsilon) ||
4106         (fabs(p.y-point.y) > MagickEpsilon))
4107       (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end last (%.20g)",
4108         (double) coordinates);
4109     else
4110       (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end open (%.20g)",
4111         (double) coordinates);
4112   }
4113 }
4114
4115 MagickExport MagickBooleanType DrawPrimitive(Image *image,
4116   const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
4117   ExceptionInfo *exception)
4118 {
4119   CacheView
4120     *image_view;
4121
4122   MagickStatusType
4123     status;
4124
4125   register ssize_t
4126     i,
4127     x;
4128
4129   ssize_t
4130     y;
4131
4132   if (image->debug != MagickFalse)
4133     {
4134       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4135         "  begin draw-primitive");
4136       (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4137         "    affine: %g %g %g %g %g %g",draw_info->affine.sx,
4138         draw_info->affine.rx,draw_info->affine.ry,draw_info->affine.sy,
4139         draw_info->affine.tx,draw_info->affine.ty);
4140     }
4141   status=MagickTrue;
4142   x=(ssize_t) ceil(primitive_info->point.x-0.5);
4143   y=(ssize_t) ceil(primitive_info->point.y-0.5);
4144   image_view=AcquireCacheView(image);
4145   switch (primitive_info->primitive)
4146   {
4147     case PointPrimitive:
4148     {
4149       PixelInfo
4150         fill_color;
4151
4152       register Quantum
4153         *q;
4154
4155       if ((y < 0) || (y >= (ssize_t) image->rows))
4156         break;
4157       if ((x < 0) || (x >= (ssize_t) image->columns))
4158         break;
4159       q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
4160       if (q == (Quantum *) NULL)
4161         break;
4162       (void) GetFillColor(draw_info,x,y,&fill_color,exception);
4163       CompositePixelOver(image,&fill_color,(MagickRealType) fill_color.alpha,q,
4164         (MagickRealType) GetPixelAlpha(image,q),q);
4165       (void) SyncCacheViewAuthenticPixels(image_view,exception);
4166       break;
4167     }
4168     case ColorPrimitive:
4169     {
4170       switch (primitive_info->method)
4171       {
4172         case PointMethod:
4173         default:
4174         {
4175           PixelInfo
4176             pixel;
4177
4178           register Quantum
4179             *q;
4180
4181           q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
4182           if (q == (Quantum *) NULL)
4183             break;
4184           GetPixelInfo(image,&pixel);
4185           (void) GetFillColor(draw_info,x,y,&pixel,exception);
4186           SetPixelInfoPixel(image,&pixel,q);
4187           (void) SyncCacheViewAuthenticPixels(image_view,exception);
4188           break;
4189         }
4190         case ReplaceMethod:
4191         {
4192           MagickBooleanType
4193             sync;
4194
4195           PixelInfo
4196             pixel,
4197             target;
4198
4199           (void) GetOneCacheViewVirtualPixelInfo(image_view,x,y,&target,
4200             exception);
4201           for (y=0; y < (ssize_t) image->rows; y++)
4202           {
4203             register Quantum
4204               *restrict q;
4205
4206             q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
4207               exception);
4208             if (q == (Quantum *) NULL)
4209               break;
4210             for (x=0; x < (ssize_t) image->columns; x++)
4211             {
4212               GetPixelInfoPixel(image,q,&pixel);
4213               if (IsFuzzyEquivalencePixelInfo(&pixel,&target) == MagickFalse)
4214                 {
4215                   q+=GetPixelChannels(image);
4216                   continue;
4217                 }
4218               (void) GetFillColor(draw_info,x,y,&pixel,exception);
4219               SetPixelInfoPixel(image,&pixel,q);
4220               q+=GetPixelChannels(image);
4221             }
4222             sync=SyncCacheViewAuthenticPixels(image_view,exception);
4223             if (sync == MagickFalse)
4224               break;
4225           }
4226           break;
4227         }
4228         case FloodfillMethod:
4229         case FillToBorderMethod:
4230         {
4231           PixelInfo
4232             target;
4233
4234           (void) GetOneVirtualPixelInfo(image,TileVirtualPixelMethod,x,y,
4235             &target,exception);
4236           if (primitive_info->method == FillToBorderMethod)
4237             {
4238               target.red=(MagickRealType) draw_info->border_color.red;
4239               target.green=(MagickRealType) draw_info->border_color.green;
4240               target.blue=(MagickRealType) draw_info->border_color.blue;
4241             }
4242           (void) FloodfillPaintImage(image,draw_info,&target,x,y,
4243             primitive_info->method == FloodfillMethod ? MagickFalse :
4244             MagickTrue,exception);
4245           break;
4246         }
4247         case ResetMethod:
4248         {
4249           MagickBooleanType
4250             sync;
4251
4252           PixelInfo
4253             pixel;
4254
4255           GetPixelInfo(image,&pixel);
4256           for (y=0; y < (ssize_t) image->rows; y++)
4257           {
4258             register Quantum
4259               *restrict q;
4260
4261             register ssize_t
4262               x;
4263
4264             q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
4265               exception);
4266             if (q == (Quantum *) NULL)
4267               break;
4268             for (x=0; x < (ssize_t) image->columns; x++)
4269             {
4270               (void) GetFillColor(draw_info,x,y,&pixel,exception);
4271               SetPixelInfoPixel(image,&pixel,q);
4272               q+=GetPixelChannels(image);
4273             }
4274             sync=SyncCacheViewAuthenticPixels(image_view,exception);
4275             if (sync == MagickFalse)
4276               break;
4277           }
4278           break;
4279         }
4280       }
4281       break;
4282     }
4283     case MattePrimitive:
4284     {
4285       if (image->matte == MagickFalse)
4286         (void) SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
4287       switch (primitive_info->method)
4288       {
4289         case PointMethod:
4290         default:
4291         {
4292           PixelInfo
4293             pixel;
4294
4295           register Quantum
4296             *q;
4297
4298           q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
4299           if (q == (Quantum *) NULL)
4300             break;
4301           (void) GetFillColor(draw_info,x,y,&pixel,exception);
4302           SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
4303           (void) SyncCacheViewAuthenticPixels(image_view,exception);
4304           break;
4305         }
4306         case ReplaceMethod:
4307         {
4308           MagickBooleanType
4309             sync;
4310
4311           PixelInfo
4312             pixel,
4313             target;
4314
4315           (void) GetOneCacheViewVirtualPixelInfo(image_view,x,y,&target,
4316             exception);
4317           for (y=0; y < (ssize_t) image->rows; y++)
4318           {
4319             register Quantum
4320               *restrict q;
4321
4322             register ssize_t
4323               x;
4324
4325             q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
4326               exception);
4327             if (q == (Quantum *) NULL)
4328               break;
4329             for (x=0; x < (ssize_t) image->columns; x++)
4330             {
4331               GetPixelInfoPixel(image,q,&pixel);
4332               if (IsFuzzyEquivalencePixelInfo(&pixel,&target) == MagickFalse)
4333                 {
4334                   q+=GetPixelChannels(image);
4335                   continue;
4336                 }
4337               (void) GetFillColor(draw_info,x,y,&pixel,exception);
4338               SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
4339               q+=GetPixelChannels(image);
4340             }
4341             sync=SyncCacheViewAuthenticPixels(image_view,exception);
4342             if (sync == MagickFalse)
4343               break;
4344           }
4345           break;
4346         }
4347         case FloodfillMethod:
4348         case FillToBorderMethod:
4349         {
4350           ChannelType
4351             channel_mask;
4352
4353           PixelInfo
4354             target;
4355
4356           (void) GetOneVirtualPixelInfo(image,TileVirtualPixelMethod,x,y,
4357             &target,exception);
4358           if (primitive_info->method == FillToBorderMethod)
4359             {
4360               target.red=(MagickRealType) draw_info->border_color.red;
4361               target.green=(MagickRealType) draw_info->border_color.green;
4362               target.blue=(MagickRealType) draw_info->border_color.blue;
4363             }
4364           channel_mask=SetPixelChannelMask(image,AlphaChannel);
4365           (void) FloodfillPaintImage(image,draw_info,&target,x,y,
4366             primitive_info->method == FloodfillMethod ? MagickFalse :
4367             MagickTrue,exception);
4368           (void) SetPixelChannelMask(image,channel_mask);
4369           break;
4370         }
4371         case ResetMethod:
4372         {
4373           MagickBooleanType
4374             sync;
4375
4376           PixelInfo
4377             pixel;
4378
4379           for (y=0; y < (ssize_t) image->rows; y++)
4380           {
4381             register Quantum
4382               *restrict q;
4383
4384             register ssize_t
4385               x;
4386
4387             q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
4388               exception);
4389             if (q == (Quantum *) NULL)
4390               break;
4391             for (x=0; x < (ssize_t) image->columns; x++)
4392             {
4393               (void) GetFillColor(draw_info,x,y,&pixel,exception);
4394               SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
4395               q+=GetPixelChannels(image);
4396             }
4397             sync=SyncCacheViewAuthenticPixels(image_view,exception);
4398             if (sync == MagickFalse)
4399               break;
4400           }
4401           break;
4402         }
4403       }
4404       break;
4405     }
4406     case TextPrimitive:
4407     {
4408       char
4409         geometry[MaxTextExtent];
4410
4411       DrawInfo
4412         *clone_info;
4413
4414       if (primitive_info->text == (char *) NULL)
4415         break;
4416       clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
4417       (void) CloneString(&clone_info->text,primitive_info->text);
4418       (void) FormatLocaleString(geometry,MaxTextExtent,"%+f%+f",
4419         primitive_info->point.x,primitive_info->point.y);
4420       (void) CloneString(&clone_info->geometry,geometry);
4421       status=AnnotateImage(image,clone_info,exception);
4422       clone_info=DestroyDrawInfo(clone_info);
4423       break;
4424     }
4425     case ImagePrimitive:
4426     {
4427       AffineMatrix
4428         affine;
4429
4430       char
4431         composite_geometry[MaxTextExtent];
4432
4433       Image
4434         *composite_image;
4435
4436       ImageInfo
4437         *clone_info;
4438
4439       RectangleInfo
4440         geometry;
4441
4442       ssize_t
4443         x1,
4444         y1;
4445
4446       if (primitive_info->text == (char *) NULL)
4447         break;
4448       clone_info=AcquireImageInfo();
4449       if (LocaleNCompare(primitive_info->text,"data:",5) == 0)
4450         composite_image=ReadInlineImage(clone_info,primitive_info->text,
4451           exception);
4452       else
4453         {
4454           (void) CopyMagickString(clone_info->filename,primitive_info->text,
4455             MaxTextExtent);
4456           composite_image=ReadImage(clone_info,exception);
4457         }
4458       clone_info=DestroyImageInfo(clone_info);
4459       if (composite_image == (Image *) NULL)
4460         break;
4461       (void) SetImageProgressMonitor(composite_image,(MagickProgressMonitor)
4462         NULL,(void *) NULL);
4463       x1=(ssize_t) ceil(primitive_info[1].point.x-0.5);
4464       y1=(ssize_t) ceil(primitive_info[1].point.y-0.5);
4465       if (((x1 != 0L) && (x1 != (ssize_t) composite_image->columns)) ||
4466           ((y1 != 0L) && (y1 != (ssize_t) composite_image->rows)))
4467         {
4468           char
4469             geometry[MaxTextExtent];
4470
4471           /*
4472             Resize image.
4473           */
4474           (void) FormatLocaleString(geometry,MaxTextExtent,"%gx%g!",
4475             primitive_info[1].point.x,primitive_info[1].point.y);
4476           composite_image->filter=image->filter;
4477           (void) TransformImage(&composite_image,(char *) NULL,geometry,
4478             exception);
4479         }
4480       if (composite_image->matte == MagickFalse)
4481         (void) SetImageAlphaChannel(composite_image,OpaqueAlphaChannel,
4482           exception);
4483       if (draw_info->alpha != OpaqueAlpha)
4484         (void) SetImageAlpha(composite_image,draw_info->alpha,exception);
4485       SetGeometry(image,&geometry);
4486       image->gravity=draw_info->gravity;
4487       geometry.x=x;
4488       geometry.y=y;
4489       (void) FormatLocaleString(composite_geometry,MaxTextExtent,
4490         "%.20gx%.20g%+.20g%+.20g",(double) composite_image->columns,(double)
4491         composite_image->rows,(double) geometry.x,(double) geometry.y);
4492       (void) ParseGravityGeometry(image,composite_geometry,&geometry,exception);
4493       affine=draw_info->affine;
4494       affine.tx=(double) geometry.x;
4495       affine.ty=(double) geometry.y;
4496       composite_image->interpolate=image->interpolate;
4497       if (draw_info->compose == OverCompositeOp)
4498         (void) DrawAffineImage(image,composite_image,&affine,exception);
4499       else
4500         (void) CompositeImage(image,composite_image,draw_info->compose,
4501           MagickTrue,geometry.x,geometry.y,exception);
4502       composite_image=DestroyImage(composite_image);
4503       break;
4504     }
4505     default:
4506     {
4507       MagickRealType
4508         mid,
4509         scale;
4510
4511       DrawInfo
4512         *clone_info;
4513
4514       if (IsEventLogging() != MagickFalse)
4515         LogPrimitiveInfo(primitive_info);
4516       scale=ExpandAffine(&draw_info->affine);
4517       if ((draw_info->dash_pattern != (double *) NULL) &&
4518           (draw_info->dash_pattern[0] != 0.0) &&
4519           ((scale*draw_info->stroke_width) > MagickEpsilon) &&
4520           (draw_info->stroke.alpha != (Quantum) TransparentAlpha))
4521         {
4522           /*
4523             Draw dash polygon.
4524           */
4525           clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
4526           clone_info->stroke_width=0.0;
4527           clone_info->stroke.alpha=(Quantum) TransparentAlpha;
4528           status=DrawPolygonPrimitive(image,clone_info,primitive_info,
4529             exception);
4530           clone_info=DestroyDrawInfo(clone_info);
4531           (void) DrawDashPolygon(draw_info,primitive_info,image,exception);
4532           break;
4533         }
4534       mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
4535       if ((mid > 1.0) &&
4536           (draw_info->stroke.alpha != (Quantum) TransparentAlpha))
4537         {
4538           MagickBooleanType
4539             closed_path;
4540
4541           /*
4542             Draw strokes while respecting line cap/join attributes.
4543           */
4544           for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
4545           closed_path=
4546             (primitive_info[i-1].point.x == primitive_info[0].point.x) &&
4547             (primitive_info[i-1].point.y == primitive_info[0].point.y) ?
4548             MagickTrue : MagickFalse;
4549           i=(ssize_t) primitive_info[0].coordinates;
4550           if ((((draw_info->linecap == RoundCap) ||
4551                 (closed_path != MagickFalse)) &&
4552                (draw_info->linejoin == RoundJoin)) ||
4553                (primitive_info[i].primitive != UndefinedPrimitive))
4554             {
4555               (void) DrawPolygonPrimitive(image,draw_info,primitive_info,
4556                 exception);
4557               break;
4558             }
4559           clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
4560           clone_info->stroke_width=0.0;
4561           clone_info->stroke.alpha=(Quantum) TransparentAlpha;
4562           status=DrawPolygonPrimitive(image,clone_info,primitive_info,
4563             exception);
4564           clone_info=DestroyDrawInfo(clone_info);
4565           status|=DrawStrokePolygon(image,draw_info,primitive_info,exception);
4566           break;
4567         }
4568       status=DrawPolygonPrimitive(image,draw_info,primitive_info,exception);
4569       break;
4570     }
4571   }
4572   image_view=DestroyCacheView(image_view);
4573   if (image->debug != MagickFalse)
4574     (void) LogMagickEvent(DrawEvent,GetMagickModule(),"  end draw-primitive");
4575   return(status != 0 ? MagickTrue : MagickFalse);
4576 }
4577 \f
4578 /*
4579 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4580 %                                                                             %
4581 %                                                                             %
4582 %                                                                             %
4583 +   D r a w S t r o k e P o l y g o n                                         %
4584 %                                                                             %
4585 %                                                                             %
4586 %                                                                             %
4587 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4588 %
4589 %  DrawStrokePolygon() draws a stroked polygon (line, rectangle, ellipse) on
4590 %  the image while respecting the line cap and join attributes.
4591 %
4592 %  The format of the DrawStrokePolygon method is:
4593 %
4594 %      MagickBooleanType DrawStrokePolygon(Image *image,
4595 %        const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
4596 %
4597 %  A description of each parameter follows:
4598 %
4599 %    o image: the image.
4600 %
4601 %    o draw_info: the draw info.
4602 %
4603 %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
4604 %
4605 %
4606 */
4607
4608 static void DrawRoundLinecap(Image *image,const DrawInfo *draw_info,
4609   const PrimitiveInfo *primitive_info,ExceptionInfo *exception)
4610 {
4611   PrimitiveInfo
4612     linecap[5];
4613
4614   register ssize_t
4615     i;
4616
4617   for (i=0; i < 4; i++)
4618     linecap[i]=(*primitive_info);
4619   linecap[0].coordinates=4;
4620   linecap[1].point.x+=(double) (10.0*MagickEpsilon);
4621   linecap[2].point.x+=(double) (10.0*MagickEpsilon);
4622   linecap[2].point.y+=(double) (10.0*MagickEpsilon);
4623   linecap[3].point.y+=(double) (10.0*MagickEpsilon);
4624   linecap[4].primitive=UndefinedPrimitive;
4625   (void) DrawPolygonPrimitive(image,draw_info,linecap,exception);
4626 }
4627
4628 static MagickBooleanType DrawStrokePolygon(Image *image,
4629   const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
4630   ExceptionInfo *exception)
4631 {
4632   DrawInfo
4633     *clone_info;
4634
4635   MagickBooleanType
4636     closed_path,
4637     status;
4638
4639   PrimitiveInfo
4640     *stroke_polygon;
4641
4642   register const PrimitiveInfo
4643     *p,
4644     *q;
4645
4646   /*
4647     Draw stroked polygon.
4648   */
4649   if (image->debug != MagickFalse)
4650     (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4651       "    begin draw-stroke-polygon");
4652   clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
4653   clone_info->fill=draw_info->stroke;
4654   if (clone_info->fill_pattern != (Image *) NULL)
4655     clone_info->fill_pattern=DestroyImage(clone_info->fill_pattern);
4656   if (clone_info->stroke_pattern != (Image *) NULL)
4657     clone_info->fill_pattern=CloneImage(clone_info->stroke_pattern,0,0,
4658       MagickTrue,exception);
4659   clone_info->stroke.alpha=(Quantum) TransparentAlpha;
4660   clone_info->stroke_width=0.0;
4661   clone_info->fill_rule=NonZeroRule;
4662   status=MagickTrue;
4663   for (p=primitive_info; p->primitive != UndefinedPrimitive; p+=p->coordinates)
4664   {
4665     stroke_polygon=TraceStrokePolygon(draw_info,p);
4666     status=DrawPolygonPrimitive(image,clone_info,stroke_polygon,exception);
4667     stroke_polygon=(PrimitiveInfo *) RelinquishMagickMemory(stroke_polygon);
4668     q=p+p->coordinates-1;
4669     closed_path=(q->point.x == p->point.x) && (q->point.y == p->point.y) ?
4670       MagickTrue : MagickFalse;
4671     if ((draw_info->linecap == RoundCap) && (closed_path == MagickFalse))
4672       {
4673         DrawRoundLinecap(image,draw_info,p,exception);
4674         DrawRoundLinecap(image,draw_info,q,exception);
4675       }
4676   }
4677   clone_info=DestroyDrawInfo(clone_info);
4678   if (image->debug != MagickFalse)
4679     (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4680       "    end draw-stroke-polygon");
4681   return(status);
4682 }
4683 \f
4684 /*
4685 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4686 %                                                                             %
4687 %                                                                             %
4688 %                                                                             %
4689 %   G e t A f f i n e M a t r i x                                             %
4690 %                                                                             %
4691 %                                                                             %
4692 %                                                                             %
4693 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4694 %
4695 %  GetAffineMatrix() returns an AffineMatrix initialized to the identity
4696 %  matrix.
4697 %
4698 %  The format of the GetAffineMatrix method is:
4699 %
4700 %      void GetAffineMatrix(AffineMatrix *affine_matrix)
4701 %
4702 %  A description of each parameter follows:
4703 %
4704 %    o affine_matrix: the affine matrix.
4705 %
4706 */
4707 MagickExport void GetAffineMatrix(AffineMatrix *affine_matrix)
4708 {
4709   (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
4710   assert(affine_matrix != (AffineMatrix *) NULL);
4711   (void) ResetMagickMemory(affine_matrix,0,sizeof(*affine_matrix));
4712   affine_matrix->sx=1.0;
4713   affine_matrix->sy=1.0;
4714 }
4715 \f
4716 /*
4717 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4718 %                                                                             %
4719 %                                                                             %
4720 %                                                                             %
4721 +   G e t D r a w I n f o                                                     %
4722 %                                                                             %
4723 %                                                                             %
4724 %                                                                             %
4725 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4726 %
4727 %  GetDrawInfo() initializes draw_info to default values from image_info.
4728 %
4729 %  The format of the GetDrawInfo method is:
4730 %
4731 %      void GetDrawInfo(const ImageInfo *image_info,DrawInfo *draw_info)
4732 %
4733 %  A description of each parameter follows:
4734 %
4735 %    o image_info: the image info..
4736 %
4737 %    o draw_info: the draw info.
4738 %
4739 */
4740 MagickExport void GetDrawInfo(const ImageInfo *image_info,DrawInfo *draw_info)
4741 {
4742   const char
4743     *option;
4744
4745   ExceptionInfo
4746     *exception;
4747
4748   /*
4749     Initialize draw attributes.
4750   */
4751   (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
4752   assert(draw_info != (DrawInfo *) NULL);
4753   (void) ResetMagickMemory(draw_info,0,sizeof(*draw_info));
4754   GetAffineMatrix(&draw_info->affine);
4755   exception=AcquireExceptionInfo();
4756   (void) QueryColorCompliance("#000F",AllCompliance,&draw_info->fill,
4757     exception);
4758   (void) QueryColorCompliance("#FFF0",AllCompliance,&draw_info->stroke,
4759     exception);
4760   draw_info->stroke_width=1.0;
4761   draw_info->alpha=OpaqueAlpha;
4762   draw_info->fill_rule=EvenOddRule;
4763   draw_info->linecap=ButtCap;
4764   draw_info->linejoin=MiterJoin;
4765   draw_info->miterlimit=10;
4766   draw_info->decorate=NoDecoration;
4767   draw_info->pointsize=12.0;
4768   draw_info->undercolor.alpha=(Quantum) TransparentAlpha;
4769   draw_info->compose=OverCompositeOp;
4770   draw_info->render=MagickTrue;
4771   draw_info->debug=IsEventLogging();
4772   if (image_info != (ImageInfo *) NULL)
4773     {
4774       draw_info->stroke_antialias=image_info->antialias;
4775       if (image_info->font != (char *) NULL)
4776         draw_info->font=AcquireString(image_info->font);
4777       if (image_info->density != (char *) NULL)
4778         draw_info->density=AcquireString(image_info->density);
4779       draw_info->text_antialias=image_info->antialias;
4780       if (image_info->pointsize != 0.0)
4781         draw_info->pointsize=image_info->pointsize;
4782       draw_info->border_color=image_info->border_color;
4783       if (image_info->server_name != (char *) NULL)
4784         draw_info->server_name=AcquireString(image_info->server_name);
4785       option=GetImageOption(image_info,"encoding");
4786       if (option != (const char *) NULL)
4787         (void) CloneString(&draw_info->encoding,option);
4788       option=GetImageOption(image_info,"kerning");
4789       if (option != (const char *) NULL)
4790         draw_info->kerning=StringToDouble(option,(char **) NULL);
4791       option=GetImageOption(image_info,"interline-spacing");
4792       if (option != (const char *) NULL)
4793         draw_info->interline_spacing=StringToDouble(option,(char **) NULL);
4794       option=GetImageOption(image_info,"interword-spacing");
4795       if (option != (const char *) NULL)
4796         draw_info->interword_spacing=StringToDouble(option,(char **) NULL);
4797       option=GetImageOption(image_info,"direction");
4798       if (option != (const char *) NULL)
4799         draw_info->direction=(DirectionType) ParseCommandOption(
4800           MagickDirectionOptions,MagickFalse,option);
4801       else
4802         draw_info->direction=UndefinedDirection;
4803       option=GetImageOption(image_info,"fill");
4804       if (option != (const char *) NULL)
4805         (void) QueryColorCompliance(option,AllCompliance,&draw_info->fill,
4806           exception);
4807       option=GetImageOption(image_info,"stroke");
4808       if (option != (const char *) NULL)
4809         (void) QueryColorCompliance(option,AllCompliance,&draw_info->stroke,
4810           exception);
4811       option=GetImageOption(image_info,"strokewidth");
4812       if (option != (const char *) NULL)
4813         draw_info->stroke_width=StringToDouble(option,(char **) NULL);
4814       option=GetImageOption(image_info,"undercolor");
4815       if (option != (const char *) NULL)
4816         (void) QueryColorCompliance(option,AllCompliance,&draw_info->undercolor,
4817           exception);
4818       option=GetImageOption(image_info,"gravity");
4819       if (option != (const char *) NULL)
4820         draw_info->gravity=(GravityType) ParseCommandOption(
4821           MagickGravityOptions,MagickFalse,option);
4822     }
4823   exception=DestroyExceptionInfo(exception);
4824   draw_info->signature=MagickSignature;
4825 }
4826 \f
4827 /*
4828 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4829 %                                                                             %
4830 %                                                                             %
4831 %                                                                             %
4832 +   P e r m u t a t e                                                         %
4833 %                                                                             %
4834 %                                                                             %
4835 %                                                                             %
4836 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4837 %
4838 %  Permutate() returns the permuation of the (n,k).
4839 %
4840 %  The format of the Permutate method is:
4841 %
4842 %      void Permutate(ssize_t n,ssize_t k)
4843 %
4844 %  A description of each parameter follows:
4845 %
4846 %    o n:
4847 %
4848 %    o k:
4849 %
4850 %
4851 */
4852 static inline MagickRealType Permutate(const ssize_t n,const ssize_t k)
4853 {
4854   MagickRealType
4855     r;
4856
4857   register ssize_t
4858     i;
4859
4860   r=1.0;
4861   for (i=k+1; i <= n; i++)
4862     r*=i;
4863   for (i=1; i <= (n-k); i++)
4864     r/=i;
4865   return(r);
4866 }
4867 \f
4868 /*
4869 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4870 %                                                                             %
4871 %                                                                             %
4872 %                                                                             %
4873 +   T r a c e P r i m i t i v e                                               %
4874 %                                                                             %
4875 %                                                                             %
4876 %                                                                             %
4877 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4878 %
4879 %  TracePrimitive is a collection of methods for generating graphic
4880 %  primitives such as arcs, ellipses, paths, etc.
4881 %
4882 */
4883
4884 static void TraceArc(PrimitiveInfo *primitive_info,const PointInfo start,
4885   const PointInfo end,const PointInfo degrees)
4886 {
4887   PointInfo
4888     center,
4889     radii;
4890
4891   center.x=0.5*(end.x+start.x);
4892   center.y=0.5*(end.y+start.y);
4893   radii.x=fabs(center.x-start.x);
4894   radii.y=fabs(center.y-start.y);
4895   TraceEllipse(primitive_info,center,radii,degrees);
4896 }
4897
4898 static void TraceArcPath(PrimitiveInfo *primitive_info,const PointInfo start,
4899   const PointInfo end,const PointInfo arc,const MagickRealType angle,
4900   const MagickBooleanType large_arc,const MagickBooleanType sweep)
4901 {
4902   MagickRealType
4903     alpha,
4904     beta,
4905     delta,
4906     factor,
4907     gamma,
4908     theta;
4909
4910   PointInfo
4911     center,
4912     points[3],
4913     radii;
4914
4915   register MagickRealType
4916     cosine,
4917     sine;
4918
4919   register PrimitiveInfo
4920     *p;
4921
4922   register ssize_t
4923     i;
4924
4925   size_t
4926     arc_segments;
4927
4928   if ((start.x == end.x) && (start.y == end.y))
4929     {
4930       TracePoint(primitive_info,end);
4931       return;
4932     }
4933   radii.x=fabs(arc.x);
4934   radii.y=fabs(arc.y);
4935   if ((radii.x == 0.0) || (radii.y == 0.0))
4936     {
4937       TraceLine(primitive_info,start,end);
4938       return;
4939     }
4940   cosine=cos(DegreesToRadians(fmod((double) angle,360.0)));
4941   sine=sin(DegreesToRadians(fmod((double) angle,360.0)));
4942   center.x=(double) (cosine*(end.x-start.x)/2+sine*(end.y-start.y)/2);
4943   center.y=(double) (cosine*(end.y-start.y)/2-sine*(end.x-start.x)/2);
4944   delta=(center.x*center.x)/(radii.x*radii.x)+(center.y*center.y)/
4945     (radii.y*radii.y);
4946   if (delta < MagickEpsilon)
4947     {
4948       TraceLine(primitive_info,start,end);
4949       return;
4950     }
4951   if (delta > 1.0)
4952     {
4953       radii.x*=sqrt((double) delta);
4954       radii.y*=sqrt((double) delta);
4955     }
4956   points[0].x=(double) (cosine*start.x/radii.x+sine*start.y/radii.x);
4957   points[0].y=(double) (cosine*start.y/radii.y-sine*start.x/radii.y);
4958   points[1].x=(double) (cosine*end.x/radii.x+sine*end.y/radii.x);
4959   points[1].y=(double) (cosine*end.y/radii.y-sine*end.x/radii.y);
4960   alpha=points[1].x-points[0].x;
4961   beta=points[1].y-points[0].y;
4962   factor=1.0/(alpha*alpha+beta*beta)-0.25;
4963   if (factor <= 0.0)
4964     factor=0.0;
4965   else
4966     {
4967       factor=sqrt((double) factor);
4968       if (sweep == large_arc)
4969         factor=(-factor);
4970     }
4971   center.x=(double) ((points[0].x+points[1].x)/2-factor*beta);
4972   center.y=(double) ((points[0].y+points[1].y)/2+factor*alpha);
4973   alpha=atan2(points[0].y-center.y,points[0].x-center.x);
4974   theta=atan2(points[1].y-center.y,points[1].x-center.x)-alpha;
4975   if ((theta < 0.0) && (sweep != MagickFalse))
4976     theta+=(MagickRealType) (2.0*MagickPI);
4977   else
4978     if ((theta > 0.0) && (sweep == MagickFalse))
4979       theta-=(MagickRealType) (2.0*MagickPI);
4980   arc_segments=(size_t) ceil(fabs((double) (theta/(0.5*MagickPI+
4981     MagickEpsilon))));
4982   p=primitive_info;
4983   for (i=0; i < (ssize_t) arc_segments; i++)
4984   {
4985     beta=0.5*((alpha+(i+1)*theta/arc_segments)-(alpha+i*theta/arc_segments));
4986     gamma=(8.0/3.0)*sin(fmod((double) (0.5*beta),DegreesToRadians(360.0)))*
4987       sin(fmod((double) (0.5*beta),DegreesToRadians(360.0)))/
4988       sin(fmod((double) beta,DegreesToRadians(360.0)));
4989     points[0].x=(double) (center.x+cos(fmod((double) (alpha+(double) i*theta/
4990       arc_segments),DegreesToRadians(360.0)))-gamma*sin(fmod((double) (alpha+
4991       (double) i*theta/arc_segments),DegreesToRadians(360.0))));
4992     points[0].y=(double) (center.y+sin(fmod((double) (alpha+(double) i*theta/
4993       arc_segments),DegreesToRadians(360.0)))+gamma*cos(fmod((double) (alpha+
4994       (double) i*theta/arc_segments),DegreesToRadians(360.0))));
4995     points[2].x=(double) (center.x+cos(fmod((double) (alpha+(double) (i+1)*
4996       theta/arc_segments),DegreesToRadians(360.0))));
4997     points[2].y=(double) (center.y+sin(fmod((double) (alpha+(double) (i+1)*
4998       theta/arc_segments),DegreesToRadians(360.0))));
4999     points[1].x=(double) (points[2].x+gamma*sin(fmod((double) (alpha+(double)
5000       (i+1)*theta/arc_segments),DegreesToRadians(360.0))));
5001     points[1].y=(double) (points[2].y-gamma*cos(fmod((double) (alpha+(double)
5002       (i+1)*theta/arc_segments),DegreesToRadians(360.0))));
5003     p->point.x=(p == primitive_info) ? start.x : (p-1)->point.x;
5004     p->point.y=(p == primitive_info) ? start.y : (p-1)->point.y;
5005     (p+1)->point.x=(double) (cosine*radii.x*points[0].x-sine*radii.y*
5006       points[0].y);
5007     (p+1)->point.y=(double) (sine*radii.x*points[0].x+cosine*radii.y*
5008       points[0].y);
5009     (p+2)->point.x=(double) (cosine*radii.x*points[1].x-sine*radii.y*
5010       points[1].y);
5011     (p+2)->point.y=(double) (sine*radii.x*points[1].x+cosine*radii.y*
5012       points[1].y);
5013     (p+3)->point.x=(double) (cosine*radii.x*points[2].x-sine*radii.y*
5014       points[2].y);
5015     (p+3)->point.y=(double) (sine*radii.x*points[2].x+cosine*radii.y*
5016       points[2].y);
5017     if (i == (ssize_t) (arc_segments-1))
5018       (p+3)->point=end;
5019     TraceBezier(p,4);
5020     p+=p->coordinates;
5021   }
5022   primitive_info->coordinates=(size_t) (p-primitive_info);
5023   for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
5024   {
5025     p->primitive=primitive_info->primitive;
5026     p--;
5027   }
5028 }
5029
5030 static void TraceBezier(PrimitiveInfo *primitive_info,
5031   const size_t number_coordinates)
5032 {
5033   MagickRealType
5034     alpha,
5035     *coefficients,
5036     weight;
5037
5038   PointInfo
5039     end,
5040     point,
5041     *points;
5042
5043   register PrimitiveInfo
5044     *p;
5045
5046   register ssize_t
5047     i,
5048     j;
5049
5050   size_t
5051     control_points,
5052     quantum;
5053
5054   /*
5055     Allocate coeficients.
5056   */
5057   quantum=number_coordinates;
5058   for (i=0; i < (ssize_t) number_coordinates; i++)
5059   {
5060     for (j=i+1; j < (ssize_t) number_coordinates; j++)
5061     {
5062       alpha=fabs(primitive_info[j].point.x-primitive_info[i].point.x);
5063       if (alpha > (MagickRealType) quantum)
5064         quantum=(size_t) alpha;
5065       alpha=fabs(primitive_info[j].point.y-primitive_info[i].point.y);
5066       if (alpha > (MagickRealType) quantum)
5067         quantum=(size_t) alpha;
5068     }
5069   }
5070   quantum=(size_t) MagickMin((double) quantum/number_coordinates,
5071     (double) BezierQuantum);
5072   control_points=quantum*number_coordinates;
5073   coefficients=(MagickRealType *) AcquireQuantumMemory((size_t)
5074     number_coordinates,sizeof(*coefficients));
5075   points=(PointInfo *) AcquireQuantumMemory((size_t) control_points,
5076     sizeof(*points));
5077   if ((coefficients == (MagickRealType *) NULL) ||
5078       (points == (PointInfo *) NULL))
5079     ThrowFatalException(ResourceLimitFatalError,"MemoryAllocationFailed");
5080   /*
5081     Compute bezier points.
5082   */
5083   end=primitive_info[number_coordinates-1].point;
5084   for (i=0; i < (ssize_t) number_coordinates; i++)
5085     coefficients[i]=Permutate((ssize_t) number_coordinates-1,i);
5086   weight=0.0;
5087   for (i=0; i < (ssize_t) control_points; i++)
5088   {
5089     p=primitive_info;
5090     point.x=0.0;
5091     point.y=0.0;
5092     alpha=pow((double) (1.0-weight),(double) number_coordinates-1.0);
5093     for (j=0; j < (ssize_t) number_coordinates; j++)
5094     {
5095       point.x+=alpha*coefficients[j]*p->point.x;
5096       point.y+=alpha*coefficients[j]*p->point.y;
5097       alpha*=weight/(1.0-weight);
5098       p++;
5099     }
5100     points[i]=point;
5101     weight+=1.0/control_points;
5102   }
5103   /*
5104     Bezier curves are just short segmented polys.
5105   */
5106   p=primitive_info;
5107   for (i=0; i < (ssize_t) control_points; i++)
5108   {
5109     TracePoint(p,points[i]);
5110     p+=p->coordinates;
5111   }
5112   TracePoint(p,end);
5113   p+=p->coordinates;
5114   primitive_info->coordinates=(size_t) (p-primitive_info);
5115   for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
5116   {
5117     p->primitive=primitive_info->primitive;
5118     p--;
5119   }
5120   points=(PointInfo *) RelinquishMagickMemory(points);
5121   coefficients=(MagickRealType *) RelinquishMagickMemory(coefficients);
5122 }
5123
5124 static void TraceCircle(PrimitiveInfo *primitive_info,const PointInfo start,
5125   const PointInfo end)
5126 {
5127   MagickRealType
5128     alpha,
5129     beta,
5130     radius;
5131
5132   PointInfo
5133     offset,
5134     degrees;
5135
5136   alpha=end.x-start.x;
5137   beta=end.y-start.y;
5138   radius=hypot((double) alpha,(double) beta);
5139   offset.x=(double) radius;
5140   offset.y=(double) radius;
5141   degrees.x=0.0;
5142   degrees.y=360.0;
5143   TraceEllipse(primitive_info,start,offset,degrees);
5144 }
5145
5146 static void TraceEllipse(PrimitiveInfo *primitive_info,const PointInfo start,
5147   const PointInfo stop,const PointInfo degrees)
5148 {
5149   MagickRealType
5150     delta,
5151     step,
5152     y;
5153
5154   PointInfo
5155     angle,
5156     point;
5157
5158   register PrimitiveInfo
5159     *p;
5160
5161   register ssize_t
5162     i;
5163
5164   /*
5165     Ellipses are just short segmented polys.
5166   */
5167   if ((stop.x == 0.0) && (stop.y == 0.0))
5168     {
5169       TracePoint(primitive_info,start);
5170       return;
5171     }
5172   delta=2.0/MagickMax(stop.x,stop.y);
5173   step=(MagickRealType) (MagickPI/8.0);
5174   if ((delta >= 0.0) && (delta < (MagickRealType) (MagickPI/8.0)))
5175     step=(MagickRealType) (MagickPI/(4*(MagickPI/delta/2+0.5)));
5176   angle.x=DegreesToRadians(degrees.x);
5177   y=degrees.y;
5178   while (y < degrees.x)
5179     y+=360.0;
5180   angle.y=(double) (DegreesToRadians(y)-MagickEpsilon);
5181   for (p=primitive_info; angle.x < angle.y; angle.x+=step)
5182   {
5183     point.x=cos(fmod(angle.x,DegreesToRadians(360.0)))*stop.x+start.x;
5184     point.y=sin(fmod(angle.x,DegreesToRadians(360.0)))*stop.y+start.y;
5185     TracePoint(p,point);
5186     p+=p->coordinates;
5187   }
5188   point.x=cos(fmod(angle.y,DegreesToRadians(360.0)))*stop.x+start.x;
5189   point.y=sin(fmod(angle.y,DegreesToRadians(360.0)))*stop.y+start.y;
5190   TracePoint(p,point);
5191   p+=p->coordinates;
5192   primitive_info->coordinates=(size_t) (p-primitive_info);
5193   for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
5194   {
5195     p->primitive=primitive_info->primitive;
5196     p--;
5197   }
5198 }
5199
5200 static void TraceLine(PrimitiveInfo *primitive_info,const PointInfo start,
5201   const PointInfo end)
5202 {
5203   TracePoint(primitive_info,start);
5204   if ((fabs(start.x-end.x) <= MagickEpsilon) &&
5205       (fabs(start.y-end.y) <= MagickEpsilon))
5206     {
5207       primitive_info->primitive=PointPrimitive;
5208       primitive_info->coordinates=1;
5209       return;
5210     }
5211   TracePoint(primitive_info+1,end);
5212   (primitive_info+1)->primitive=primitive_info->primitive;
5213   primitive_info->coordinates=2;
5214 }
5215
5216 static size_t TracePath(PrimitiveInfo *primitive_info,const char *path)
5217 {
5218   char
5219     token[MaxTextExtent];
5220
5221   const char
5222     *p;
5223
5224   int
5225     attribute,
5226     last_attribute;
5227
5228   MagickRealType
5229     x,
5230     y;
5231
5232   PointInfo
5233     end,
5234     points[4],
5235     point,
5236     start;
5237
5238   PrimitiveType
5239     primitive_type;
5240
5241   register PrimitiveInfo
5242     *q;
5243
5244   register ssize_t
5245     i;
5246
5247   size_t
5248     number_coordinates,
5249     z_count;
5250
5251   attribute=0;
5252   point.x=0.0;
5253   point.y=0.0;
5254   start.x=0.0;
5255   start.y=0.0;
5256   number_coordinates=0;
5257   z_count=0;
5258   primitive_type=primitive_info->primitive;
5259   q=primitive_info;
5260   for (p=path; *p != '\0'; )
5261   {
5262     while (isspace((int) ((unsigned char) *p)) != 0)
5263       p++;
5264     if (*p == '\0')
5265       break;
5266     last_attribute=attribute;
5267     attribute=(int) (*p++);
5268     switch (attribute)
5269     {
5270       case 'a':
5271       case 'A':
5272       {
5273         MagickBooleanType
5274           large_arc,
5275           sweep;
5276
5277         MagickRealType
5278           angle;
5279
5280         PointInfo
5281           arc;
5282
5283         /*
5284           Compute arc points.
5285         */
5286         do
5287         {
5288           GetMagickToken(p,&p,token);
5289           if (*token == ',')
5290             GetMagickToken(p,&p,token);
5291           arc.x=StringToDouble(token,(char **) NULL);
5292           GetMagickToken(p,&p,token);
5293           if (*token == ',')
5294             GetMagickToken(p,&p,token);
5295           arc.y=StringToDouble(token,(char **) NULL);
5296           GetMagickToken(p,&p,token);
5297           if (*token == ',')
5298             GetMagickToken(p,&p,token);
5299           angle=StringToDouble(token,(char **) NULL);
5300           GetMagickToken(p,&p,token);
5301           if (*token == ',')
5302             GetMagickToken(p,&p,token);
5303           large_arc=StringToLong(token) != 0 ? MagickTrue : MagickFalse;
5304           GetMagickToken(p,&p,token);
5305           if (*token == ',')
5306             GetMagickToken(p,&p,token);
5307           sweep=StringToLong(token) != 0 ? MagickTrue : MagickFalse;
5308           GetMagickToken(p,&p,token);
5309           if (*token == ',')
5310             GetMagickToken(p,&p,token);
5311           x=StringToDouble(token,(char **) NULL);
5312           GetMagickToken(p,&p,token);
5313           if (*token == ',')
5314             GetMagickToken(p,&p,token);
5315           y=StringToDouble(token,(char **) NULL);
5316           end.x=(double) (attribute == (int) 'A' ? x : point.x+x);
5317           end.y=(double) (attribute == (int) 'A' ? y : point.y+y);
5318           TraceArcPath(q,point,end,arc,angle,large_arc,sweep);
5319           q+=q->coordinates;
5320           point=end;
5321           while (isspace((int) ((unsigned char) *p)) != 0)
5322             p++;
5323           if (*p == ',')
5324             p++;
5325         } while (IsPoint(p) != MagickFalse);
5326         break;
5327       }
5328       case 'c':
5329       case 'C':
5330       {
5331         /*
5332           Compute bezier points.
5333         */
5334         do
5335         {
5336           points[0]=point;
5337           for (i=1; i < 4; i++)
5338           {
5339             GetMagickToken(p,&p,token);
5340             if (*token == ',')
5341               GetMagickToken(p,&p,token);
5342             x=StringToDouble(token,(char **) NULL);
5343             GetMagickToken(p,&p,token);
5344             if (*token == ',')
5345               GetMagickToken(p,&p,token);
5346             y=StringToDouble(token,(char **) NULL);
5347             end.x=(double) (attribute == (int) 'C' ? x : point.x+x);
5348             end.y=(double) (attribute == (int) 'C' ? y : point.y+y);
5349             points[i]=end;
5350           }
5351           for (i=0; i < 4; i++)
5352             (q+i)->point=points[i];
5353           TraceBezier(q,4);
5354           q+=q->coordinates;
5355           point=end;
5356         } while (IsPoint(p) != MagickFalse);
5357         break;
5358       }
5359       case 'H':
5360       case 'h':
5361       {
5362         do
5363         {
5364           GetMagickToken(p,&p,token);
5365           if (*token == ',')
5366             GetMagickToken(p,&p,token);
5367           x=StringToDouble(token,(char **) NULL);
5368           point.x=(double) (attribute == (int) 'H' ? x: point.x+x);
5369           TracePoint(q,point);
5370           q+=q->coordinates;
5371         } while (IsPoint(p) != MagickFalse);
5372         break;
5373       }
5374       case 'l':
5375       case 'L':
5376       {
5377         do
5378         {
5379           GetMagickToken(p,&p,token);
5380           if (*token == ',')
5381             GetMagickToken(p,&p,token);
5382           x=StringToDouble(token,(char **) NULL);
5383           GetMagickToken(p,&p,token);
5384           if (*token == ',')
5385             GetMagickToken(p,&p,token);
5386           y=StringToDouble(token,(char **) NULL);
5387           point.x=(double) (attribute == (int) 'L' ? x : point.x+x);
5388           point.y=(double) (attribute == (int) 'L' ? y : point.y+y);
5389           TracePoint(q,point);
5390           q+=q->coordinates;
5391         } while (IsPoint(p) != MagickFalse);
5392         break;
5393       }
5394       case 'M':
5395       case 'm':
5396       {
5397         if (q != primitive_info)
5398           {
5399             primitive_info->coordinates=(size_t) (q-primitive_info);
5400             number_coordinates+=primitive_info->coordinates;
5401             primitive_info=q;
5402           }
5403         i=0;
5404         do
5405         {
5406           GetMagickToken(p,&p,token);
5407           if (*token == ',')
5408             GetMagickToken(p,&p,token);
5409           x=StringToDouble(token,(char **) NULL);
5410           GetMagickToken(p,&p,token);
5411           if (*token == ',')
5412             GetMagickToken(p,&p,token);
5413           y=StringToDouble(token,(char **) NULL);
5414           point.x=(double) (attribute == (int) 'M' ? x : point.x+x);
5415           point.y=(double) (attribute == (int) 'M' ? y : point.y+y);
5416           if (i == 0)
5417             start=point;
5418           i++;
5419           TracePoint(q,point);
5420           q+=q->coordinates;
5421           if ((i != 0) && (attribute == (int) 'M'))
5422             {
5423               TracePoint(q,point);
5424               q+=q->coordinates;
5425             }
5426         } while (IsPoint(p) != MagickFalse);
5427         break;
5428       }
5429       case 'q':
5430       case 'Q':
5431       {
5432         /*
5433           Compute bezier points.
5434         */
5435         do
5436         {
5437           points[0]=point;
5438           for (i=1; i < 3; i++)
5439           {
5440             GetMagickToken(p,&p,token);
5441             if (*token == ',')
5442               GetMagickToken(p,&p,token);
5443             x=StringToDouble(token,(char **) NULL);
5444             GetMagickToken(p,&p,token);
5445             if (*token == ',')
5446               GetMagickToken(p,&p,token);
5447             y=StringToDouble(token,(char **) NULL);
5448             if (*p == ',')
5449               p++;
5450             end.x=(double) (attribute == (int) 'Q' ? x : point.x+x);
5451             end.y=(double) (attribute == (int) 'Q' ? y : point.y+y);
5452             points[i]=end;
5453           }
5454           for (i=0; i < 3; i++)
5455             (q+i)->point=points[i];
5456           TraceBezier(q,3);
5457           q+=q->coordinates;
5458           point=end;
5459         } while (IsPoint(p) != MagickFalse);
5460         break;
5461       }
5462       case 's':
5463       case 'S':
5464       {
5465         /*
5466           Compute bezier points.
5467         */
5468         do
5469         {
5470           points[0]=points[3];
5471           points[1].x=2.0*points[3].x-points[2].x;
5472           points[1].y=2.0*points[3].y-points[2].y;
5473           for (i=2; i < 4; i++)
5474           {
5475             GetMagickToken(p,&p,token);
5476             if (*token == ',')
5477               GetMagickToken(p,&p,token);
5478             x=StringToDouble(token,(char **) NULL);
5479             GetMagickToken(p,&p,token);
5480             if (*token == ',')
5481               GetMagickToken(p,&p,token);
5482             y=StringToDouble(token,(char **) NULL);
5483             if (*p == ',')
5484               p++;
5485             end.x=(double) (attribute == (int) 'S' ? x : point.x+x);
5486             end.y=(double) (attribute == (int) 'S' ? y : point.y+y);
5487             points[i]=end;
5488           }
5489           if (strchr("CcSs",last_attribute) == (char *) NULL)
5490             {
5491               points[0]=points[2];
5492               points[1]=points[3];
5493             }
5494           for (i=0; i < 4; i++)
5495             (q+i)->point=points[i];
5496           TraceBezier(q,4);
5497           q+=q->coordinates;
5498           point=end;
5499         } while (IsPoint(p) != MagickFalse);
5500         break;
5501       }
5502       case 't':
5503       case 'T':
5504       {
5505         /*
5506           Compute bezier points.
5507         */
5508         do
5509         {
5510           points[0]=points[2];
5511           points[1].x=2.0*points[2].x-points[1].x;
5512           points[1].y=2.0*points[2].y-points[1].y;
5513           for (i=2; i < 3; i++)
5514           {
5515             GetMagickToken(p,&p,token);
5516             if (*token == ',')
5517               GetMagickToken(p,&p,token);
5518             x=StringToDouble(token,(char **) NULL);
5519             GetMagickToken(p,&p,token);
5520             if (*token == ',')
5521               GetMagickToken(p,&p,token);
5522             y=StringToDouble(token,(char **) NULL);
5523             end.x=(double) (attribute == (int) 'T' ? x : point.x+x);
5524             end.y=(double) (attribute == (int) 'T' ? y : point.y+y);
5525             points[i]=end;
5526           }
5527           if (strchr("QqTt",last_attribute) == (char *) NULL)
5528             {
5529               points[0]=points[2];
5530               points[1]=points[3];
5531             }
5532           for (i=0; i < 3; i++)
5533             (q+i)->point=points[i];
5534           TraceBezier(q,3);
5535           q+=q->coordinates;
5536           point=end;
5537         } while (IsPoint(p) != MagickFalse);
5538         break;
5539       }
5540       case 'v':
5541       case 'V':
5542       {
5543         do
5544         {
5545           GetMagickToken(p,&p,token);
5546           if (*token == ',')
5547             GetMagickToken(p,&p,token);
5548           y=StringToDouble(token,(char **) NULL);
5549           point.y=(double) (attribute == (int) 'V' ? y : point.y+y);
5550           TracePoint(q,point);
5551           q+=q->coordinates;
5552         } while (IsPoint(p) != MagickFalse);
5553         break;
5554       }
5555       case 'z':
5556       case 'Z':
5557       {
5558         point=start;
5559         TracePoint(q,point);
5560         q+=q->coordinates;
5561         primitive_info->coordinates=(size_t) (q-primitive_info);
5562         number_coordinates+=primitive_info->coordinates;
5563         primitive_info=q;
5564         z_count++;
5565         break;
5566       }
5567       default:
5568       {
5569         if (isalpha((int) ((unsigned char) attribute)) != 0)
5570           (void) FormatLocaleFile(stderr,"attribute not recognized: %c\n",
5571             attribute);
5572         break;
5573       }
5574     }
5575   }
5576   primitive_info->coordinates=(size_t) (q-primitive_info);
5577   number_coordinates+=primitive_info->coordinates;
5578   for (i=0; i < (ssize_t) number_coordinates; i++)
5579   {
5580     q--;
5581     q->primitive=primitive_type;
5582     if (z_count > 1)
5583       q->method=FillToBorderMethod;
5584   }
5585   q=primitive_info;
5586   return(number_coordinates);
5587 }
5588
5589 static void TraceRectangle(PrimitiveInfo *primitive_info,const PointInfo start,
5590   const PointInfo end)
5591 {
5592   PointInfo
5593     point;
5594
5595   register PrimitiveInfo
5596     *p;
5597
5598   register ssize_t
5599     i;
5600
5601   p=primitive_info;
5602   TracePoint(p,start);
5603   p+=p->coordinates;
5604   point.x=start.x;
5605   point.y=end.y;
5606   TracePoint(p,point);
5607   p+=p->coordinates;
5608   TracePoint(p,end);
5609   p+=p->coordinates;
5610   point.x=end.x;
5611   point.y=start.y;
5612   TracePoint(p,point);
5613   p+=p->coordinates;
5614   TracePoint(p,start);
5615   p+=p->coordinates;
5616   primitive_info->coordinates=(size_t) (p-primitive_info);
5617   for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
5618   {
5619     p->primitive=primitive_info->primitive;
5620     p--;
5621   }
5622 }
5623
5624 static void TraceRoundRectangle(PrimitiveInfo *primitive_info,
5625   const PointInfo start,const PointInfo end,PointInfo arc)
5626 {
5627   PointInfo
5628     degrees,
5629     offset,
5630     point;
5631
5632   register PrimitiveInfo
5633     *p;
5634
5635   register ssize_t
5636     i;
5637
5638   p=primitive_info;
5639   offset.x=fabs(end.x-start.x);
5640   offset.y=fabs(end.y-start.y);
5641   if (arc.x > (0.5*offset.x))
5642     arc.x=0.5*offset.x;
5643   if (arc.y > (0.5*offset.y))
5644     arc.y=0.5*offset.y;
5645   point.x=start.x+offset.x-arc.x;
5646   point.y=start.y+arc.y;
5647   degrees.x=270.0;
5648   degrees.y=360.0;
5649   TraceEllipse(p,point,arc,degrees);
5650   p+=p->coordinates;
5651   point.x=start.x+offset.x-arc.x;
5652   point.y=start.y+offset.y-arc.y;
5653   degrees.x=0.0;
5654   degrees.y=90.0;
5655   TraceEllipse(p,point,arc,degrees);
5656   p+=p->coordinates;
5657   point.x=start.x+arc.x;
5658   point.y=start.y+offset.y-arc.y;
5659   degrees.x=90.0;
5660   degrees.y=180.0;
5661   TraceEllipse(p,point,arc,degrees);
5662   p+=p->coordinates;
5663   point.x=start.x+arc.x;
5664   point.y=start.y+arc.y;
5665   degrees.x=180.0;
5666   degrees.y=270.0;
5667   TraceEllipse(p,point,arc,degrees);
5668   p+=p->coordinates;
5669   TracePoint(p,primitive_info->point);
5670   p+=p->coordinates;
5671   primitive_info->coordinates=(size_t) (p-primitive_info);
5672   for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
5673   {
5674     p->primitive=primitive_info->primitive;
5675     p--;
5676   }
5677 }
5678
5679 static void TraceSquareLinecap(PrimitiveInfo *primitive_info,
5680   const size_t number_vertices,const MagickRealType offset)
5681 {
5682   MagickRealType
5683     distance;
5684
5685   register MagickRealType
5686     dx,
5687     dy;
5688
5689   register ssize_t
5690     i;
5691
5692   ssize_t
5693     j;
5694
5695   dx=0.0;
5696   dy=0.0;
5697   for (i=1; i < (ssize_t) number_vertices; i++)
5698   {
5699     dx=primitive_info[0].point.x-primitive_info[i].point.x;
5700     dy=primitive_info[0].point.y-primitive_info[i].point.y;
5701     if ((fabs((double) dx) >= MagickEpsilon) ||
5702         (fabs((double) dy) >= MagickEpsilon))
5703       break;
5704   }
5705   if (i == (ssize_t) number_vertices)
5706     i=(ssize_t) number_vertices-1L;
5707   distance=hypot((double) dx,(double) dy);
5708   primitive_info[0].point.x=(double) (primitive_info[i].point.x+
5709     dx*(distance+offset)/distance);
5710   primitive_info[0].point.y=(double) (primitive_info[i].point.y+
5711     dy*(distance+offset)/distance);
5712   for (j=(ssize_t) number_vertices-2; j >= 0;  j--)
5713   {
5714     dx=primitive_info[number_vertices-1].point.x-primitive_info[j].point.x;
5715     dy=primitive_info[number_vertices-1].point.y-primitive_info[j].point.y;
5716     if ((fabs((double) dx) >= MagickEpsilon) ||
5717         (fabs((double) dy) >= MagickEpsilon))
5718       break;
5719   }
5720   distance=hypot((double) dx,(double) dy);
5721   primitive_info[number_vertices-1].point.x=(double) (primitive_info[j].point.x+
5722     dx*(distance+offset)/distance);
5723   primitive_info[number_vertices-1].point.y=(double) (primitive_info[j].point.y+
5724     dy*(distance+offset)/distance);
5725 }
5726
5727 static PrimitiveInfo *TraceStrokePolygon(const DrawInfo *draw_info,
5728   const PrimitiveInfo *primitive_info)
5729 {
5730   typedef struct _LineSegment
5731   {
5732     double
5733       p,
5734       q;
5735   } LineSegment;
5736
5737   LineSegment
5738     dx,
5739     dy,
5740     inverse_slope,
5741     slope,
5742     theta;
5743
5744   MagickBooleanType
5745     closed_path;
5746
5747   MagickRealType
5748     delta_theta,
5749     dot_product,
5750     mid,
5751     miterlimit;
5752
5753   PointInfo
5754     box_p[5],
5755     box_q[5],
5756     center,
5757     offset,
5758     *path_p,
5759     *path_q;
5760
5761   PrimitiveInfo
5762     *polygon_primitive,
5763     *stroke_polygon;
5764
5765   register ssize_t
5766     i;
5767
5768   size_t
5769     arc_segments,
5770     max_strokes,
5771     number_vertices;
5772
5773   ssize_t
5774     j,
5775     n,
5776     p,
5777     q;
5778
5779   /*
5780     Allocate paths.
5781   */
5782   number_vertices=primitive_info->coordinates;
5783   max_strokes=2*number_vertices+6*BezierQuantum+360;
5784   path_p=(PointInfo *) AcquireQuantumMemory((size_t) max_strokes,
5785     sizeof(*path_p));
5786   path_q=(PointInfo *) AcquireQuantumMemory((size_t) max_strokes,
5787     sizeof(*path_q));
5788   polygon_primitive=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
5789     number_vertices+2UL,sizeof(*polygon_primitive));
5790   if ((path_p == (PointInfo *) NULL) || (path_q == (PointInfo *) NULL) ||
5791       (polygon_primitive == (PrimitiveInfo *) NULL))
5792     return((PrimitiveInfo *) NULL);
5793   (void) CopyMagickMemory(polygon_primitive,primitive_info,(size_t)
5794     number_vertices*sizeof(*polygon_primitive));
5795   closed_path=
5796     (primitive_info[number_vertices-1].point.x == primitive_info[0].point.x) &&
5797     (primitive_info[number_vertices-1].point.y == primitive_info[0].point.y) ?
5798     MagickTrue : MagickFalse;
5799   if ((draw_info->linejoin == RoundJoin) ||
5800       ((draw_info->linejoin == MiterJoin) && (closed_path != MagickFalse)))
5801     {
5802       polygon_primitive[number_vertices]=primitive_info[1];
5803       number_vertices++;
5804     }
5805   polygon_primitive[number_vertices].primitive=UndefinedPrimitive;
5806   /*
5807     Compute the slope for the first line segment, p.
5808   */
5809   dx.p=0.0;
5810   dy.p=0.0;
5811   for (n=1; n < (ssize_t) number_vertices; n++)
5812   {
5813     dx.p=polygon_primitive[n].point.x-polygon_primitive[0].point.x;
5814     dy.p=polygon_primitive[n].point.y-polygon_primitive[0].point.y;
5815     if ((fabs(dx.p) >= MagickEpsilon) || (fabs(dy.p) >= MagickEpsilon))
5816       break;
5817   }
5818   if (n == (ssize_t) number_vertices)
5819     n=(ssize_t) number_vertices-1L;
5820   slope.p=0.0;
5821   inverse_slope.p=0.0;
5822   if (fabs(dx.p) <= MagickEpsilon)
5823     {
5824       if (dx.p >= 0.0)
5825         slope.p=dy.p < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
5826       else
5827         slope.p=dy.p < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
5828     }
5829   else
5830     if (fabs(dy.p) <= MagickEpsilon)
5831       {
5832         if (dy.p >= 0.0)
5833           inverse_slope.p=dx.p < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
5834         else
5835           inverse_slope.p=dx.p < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
5836       }
5837     else
5838       {
5839         slope.p=dy.p/dx.p;
5840         inverse_slope.p=(-1.0/slope.p);
5841       }
5842   mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
5843   miterlimit=(MagickRealType) (draw_info->miterlimit*draw_info->miterlimit*
5844     mid*mid);
5845   if ((draw_info->linecap == SquareCap) && (closed_path == MagickFalse))
5846     TraceSquareLinecap(polygon_primitive,number_vertices,mid);
5847   offset.x=sqrt((double) (mid*mid/(inverse_slope.p*inverse_slope.p+1.0)));
5848   offset.y=(double) (offset.x*inverse_slope.p);
5849   if ((dy.p*offset.x-dx.p*offset.y) > 0.0)
5850     {
5851       box_p[0].x=polygon_primitive[0].point.x-offset.x;
5852       box_p[0].y=polygon_primitive[0].point.y-offset.x*inverse_slope.p;
5853       box_p[1].x=polygon_primitive[n].point.x-offset.x;
5854       box_p[1].y=polygon_primitive[n].point.y-offset.x*inverse_slope.p;
5855       box_q[0].x=polygon_primitive[0].point.x+offset.x;
5856       box_q[0].y=polygon_primitive[0].point.y+offset.x*inverse_slope.p;
5857       box_q[1].x=polygon_primitive[n].point.x+offset.x;
5858       box_q[1].y=polygon_primitive[n].point.y+offset.x*inverse_slope.p;
5859     }
5860   else
5861     {
5862       box_p[0].x=polygon_primitive[0].point.x+offset.x;
5863       box_p[0].y=polygon_primitive[0].point.y+offset.y;
5864       box_p[1].x=polygon_primitive[n].point.x+offset.x;
5865       box_p[1].y=polygon_primitive[n].point.y+offset.y;
5866       box_q[0].x=polygon_primitive[0].point.x-offset.x;
5867       box_q[0].y=polygon_primitive[0].point.y-offset.y;
5868       box_q[1].x=polygon_primitive[n].point.x-offset.x;
5869       box_q[1].y=polygon_primitive[n].point.y-offset.y;
5870     }
5871   /*
5872     Create strokes for the line join attribute: bevel, miter, round.
5873   */
5874   p=0;
5875   q=0;
5876   path_q[p++]=box_q[0];
5877   path_p[q++]=box_p[0];
5878   for (i=(ssize_t) n+1; i < (ssize_t) number_vertices; i++)
5879   {
5880     /*
5881       Compute the slope for this line segment, q.
5882     */
5883     dx.q=polygon_primitive[i].point.x-polygon_primitive[n].point.x;
5884     dy.q=polygon_primitive[i].point.y-polygon_primitive[n].point.y;
5885     dot_product=dx.q*dx.q+dy.q*dy.q;
5886     if (dot_product < 0.25)
5887       continue;
5888     slope.q=0.0;
5889     inverse_slope.q=0.0;
5890     if (fabs(dx.q) < MagickEpsilon)
5891       {
5892         if (dx.q >= 0.0)
5893           slope.q=dy.q < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
5894         else
5895           slope.q=dy.q < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
5896       }
5897     else
5898       if (fabs(dy.q) <= MagickEpsilon)
5899         {
5900           if (dy.q >= 0.0)
5901             inverse_slope.q=dx.q < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
5902           else
5903             inverse_slope.q=dx.q < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
5904         }
5905       else
5906         {
5907           slope.q=dy.q/dx.q;
5908           inverse_slope.q=(-1.0/slope.q);
5909         }
5910     offset.x=sqrt((double) (mid*mid/(inverse_slope.q*inverse_slope.q+1.0)));
5911     offset.y=(double) (offset.x*inverse_slope.q);
5912     dot_product=dy.q*offset.x-dx.q*offset.y;
5913     if (dot_product > 0.0)
5914       {
5915         box_p[2].x=polygon_primitive[n].point.x-offset.x;
5916         box_p[2].y=polygon_primitive[n].point.y-offset.y;
5917         box_p[3].x=polygon_primitive[i].point.x-offset.x;
5918         box_p[3].y=polygon_primitive[i].point.y-offset.y;
5919         box_q[2].x=polygon_primitive[n].point.x+offset.x;
5920         box_q[2].y=polygon_primitive[n].point.y+offset.y;
5921         box_q[3].x=polygon_primitive[i].point.x+offset.x;
5922         box_q[3].y=polygon_primitive[i].point.y+offset.y;
5923       }
5924     else
5925       {
5926         box_p[2].x=polygon_primitive[n].point.x+offset.x;
5927         box_p[2].y=polygon_primitive[n].point.y+offset.y;
5928         box_p[3].x=polygon_primitive[i].point.x+offset.x;
5929         box_p[3].y=polygon_primitive[i].point.y+offset.y;
5930         box_q[2].x=polygon_primitive[n].point.x-offset.x;
5931         box_q[2].y=polygon_primitive[n].point.y-offset.y;
5932         box_q[3].x=polygon_primitive[i].point.x-offset.x;
5933         box_q[3].y=polygon_primitive[i].point.y-offset.y;
5934       }
5935     if (fabs((double) (slope.p-slope.q)) <= MagickEpsilon)
5936       {
5937         box_p[4]=box_p[1];
5938         box_q[4]=box_q[1];
5939       }
5940     else
5941       {
5942         box_p[4].x=(double) ((slope.p*box_p[0].x-box_p[0].y-slope.q*box_p[3].x+
5943           box_p[3].y)/(slope.p-slope.q));
5944         box_p[4].y=(double) (slope.p*(box_p[4].x-box_p[0].x)+box_p[0].y);
5945         box_q[4].x=(double) ((slope.p*box_q[0].x-box_q[0].y-slope.q*box_q[3].x+
5946           box_q[3].y)/(slope.p-slope.q));
5947         box_q[4].y=(double) (slope.p*(box_q[4].x-box_q[0].x)+box_q[0].y);
5948       }
5949     if (q >= (ssize_t) (max_strokes-6*BezierQuantum-360))
5950       {
5951          max_strokes+=6*BezierQuantum+360;
5952          path_p=(PointInfo *) ResizeQuantumMemory(path_p,(size_t) max_strokes,
5953            sizeof(*path_p));
5954          path_q=(PointInfo *) ResizeQuantumMemory(path_q,(size_t) max_strokes,
5955            sizeof(*path_q));
5956          if ((path_p == (PointInfo *) NULL) || (path_q == (PointInfo *) NULL))
5957            {
5958              polygon_primitive=(PrimitiveInfo *)
5959                RelinquishMagickMemory(polygon_primitive);
5960              return((PrimitiveInfo *) NULL);
5961            }
5962       }
5963     dot_product=dx.q*dy.p-dx.p*dy.q;
5964     if (dot_product <= 0.0)
5965       switch (draw_info->linejoin)
5966       {
5967         case BevelJoin:
5968         {
5969           path_q[q++]=box_q[1];
5970           path_q[q++]=box_q[2];
5971           dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
5972             (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
5973           if (dot_product <= miterlimit)
5974             path_p[p++]=box_p[4];
5975           else
5976             {
5977               path_p[p++]=box_p[1];
5978               path_p[p++]=box_p[2];
5979             }
5980           break;
5981         }
5982         case MiterJoin:
5983         {
5984           dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
5985             (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
5986           if (dot_product <= miterlimit)
5987             {
5988               path_q[q++]=box_q[4];
5989               path_p[p++]=box_p[4];
5990             }
5991           else
5992             {
5993               path_q[q++]=box_q[1];
5994               path_q[q++]=box_q[2];
5995               path_p[p++]=box_p[1];
5996               path_p[p++]=box_p[2];
5997             }
5998           break;
5999         }
6000         case RoundJoin:
6001         {
6002           dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
6003             (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
6004           if (dot_product <= miterlimit)
6005             path_p[p++]=box_p[4];
6006           else
6007             {
6008               path_p[p++]=box_p[1];
6009               path_p[p++]=box_p[2];
6010             }
6011           center=polygon_primitive[n].point;
6012           theta.p=atan2(box_q[1].y-center.y,box_q[1].x-center.x);
6013           theta.q=atan2(box_q[2].y-center.y,box_q[2].x-center.x);
6014           if (theta.q < theta.p)
6015             theta.q+=(MagickRealType) (2.0*MagickPI);
6016           arc_segments=(size_t) ceil((double) ((theta.q-theta.p)/
6017             (2.0*sqrt((double) (1.0/mid)))));
6018           path_q[q].x=box_q[1].x;
6019           path_q[q].y=box_q[1].y;
6020           q++;
6021           for (j=1; j < (ssize_t) arc_segments; j++)
6022           {
6023             delta_theta=(MagickRealType) (j*(theta.q-theta.p)/arc_segments);
6024             path_q[q].x=(double) (center.x+mid*cos(fmod((double)
6025               (theta.p+delta_theta),DegreesToRadians(360.0))));
6026             path_q[q].y=(double) (center.y+mid*sin(fmod((double)
6027               (theta.p+delta_theta),DegreesToRadians(360.0))));
6028             q++;
6029           }
6030           path_q[q++]=box_q[2];
6031           break;
6032         }
6033         default:
6034           break;
6035       }
6036     else
6037       switch (draw_info->linejoin)
6038       {
6039         case BevelJoin:
6040         {
6041           path_p[p++]=box_p[1];
6042           path_p[p++]=box_p[2];
6043           dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
6044             (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
6045           if (dot_product <= miterlimit)
6046             path_q[q++]=box_q[4];
6047           else
6048             {
6049               path_q[q++]=box_q[1];
6050               path_q[q++]=box_q[2];
6051             }
6052           break;
6053         }
6054         case MiterJoin:
6055         {
6056           dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
6057             (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
6058           if (dot_product <= miterlimit)
6059             {
6060               path_q[q++]=box_q[4];
6061               path_p[p++]=box_p[4];
6062             }
6063           else
6064             {
6065               path_q[q++]=box_q[1];
6066               path_q[q++]=box_q[2];
6067               path_p[p++]=box_p[1];
6068               path_p[p++]=box_p[2];
6069             }
6070           break;
6071         }
6072         case RoundJoin:
6073         {
6074           dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
6075             (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
6076           if (dot_product <= miterlimit)
6077             path_q[q++]=box_q[4];
6078           else
6079             {
6080               path_q[q++]=box_q[1];
6081               path_q[q++]=box_q[2];
6082             }
6083           center=polygon_primitive[n].point;
6084           theta.p=atan2(box_p[1].y-center.y,box_p[1].x-center.x);
6085           theta.q=atan2(box_p[2].y-center.y,box_p[2].x-center.x);
6086           if (theta.p < theta.q)
6087             theta.p+=(MagickRealType) (2.0*MagickPI);
6088           arc_segments=(size_t) ceil((double) ((theta.p-theta.q)/
6089             (2.0*sqrt((double) (1.0/mid)))));
6090           path_p[p++]=box_p[1];
6091           for (j=1; j < (ssize_t) arc_segments; j++)
6092           {
6093             delta_theta=(MagickRealType) (j*(theta.q-theta.p)/arc_segments);
6094             path_p[p].x=(double) (center.x+mid*cos(fmod((double)
6095               (theta.p+delta_theta),DegreesToRadians(360.0))));
6096             path_p[p].y=(double) (center.y+mid*sin(fmod((double)
6097               (theta.p+delta_theta),DegreesToRadians(360.0))));
6098             p++;
6099           }
6100           path_p[p++]=box_p[2];
6101           break;
6102         }
6103         default:
6104           break;
6105       }
6106     slope.p=slope.q;
6107     inverse_slope.p=inverse_slope.q;
6108     box_p[0]=box_p[2];
6109     box_p[1]=box_p[3];
6110     box_q[0]=box_q[2];
6111     box_q[1]=box_q[3];
6112     dx.p=dx.q;
6113     dy.p=dy.q;
6114     n=i;
6115   }
6116   path_p[p++]=box_p[1];
6117   path_q[q++]=box_q[1];
6118   /*
6119     Trace stroked polygon.
6120   */
6121   stroke_polygon=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
6122     (p+q+2UL*closed_path+2UL),sizeof(*stroke_polygon));
6123   if (stroke_polygon != (PrimitiveInfo *) NULL)
6124     {
6125       for (i=0; i < (ssize_t) p; i++)
6126       {
6127         stroke_polygon[i]=polygon_primitive[0];
6128         stroke_polygon[i].point=path_p[i];
6129       }
6130       if (closed_path != MagickFalse)
6131         {
6132           stroke_polygon[i]=polygon_primitive[0];
6133           stroke_polygon[i].point=stroke_polygon[0].point;
6134           i++;
6135         }
6136       for ( ; i < (ssize_t) (p+q+closed_path); i++)
6137       {
6138         stroke_polygon[i]=polygon_primitive[0];
6139         stroke_polygon[i].point=path_q[p+q+closed_path-(i+1)];
6140       }
6141       if (closed_path != MagickFalse)
6142         {
6143           stroke_polygon[i]=polygon_primitive[0];
6144           stroke_polygon[i].point=stroke_polygon[p+closed_path].point;
6145           i++;
6146         }
6147       stroke_polygon[i]=polygon_primitive[0];
6148       stroke_polygon[i].point=stroke_polygon[0].point;
6149       i++;
6150       stroke_polygon[i].primitive=UndefinedPrimitive;
6151       stroke_polygon[0].coordinates=(size_t) (p+q+2*closed_path+1);
6152     }
6153   path_p=(PointInfo *) RelinquishMagickMemory(path_p);
6154   path_q=(PointInfo *) RelinquishMagickMemory(path_q);
6155   polygon_primitive=(PrimitiveInfo *) RelinquishMagickMemory(polygon_primitive);
6156   return(stroke_polygon);
6157 }