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