MagickCore  6.9.12-93
Convert, Edit, Or Compose Bitmap Images
 All Data Structures
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 % Cristy %
18 % July 1998 %
19 % %
20 % %
21 % Copyright 1999 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 % https://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 
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-private.h"
54 #include "magick/cache-view.h"
55 #include "magick/channel.h"
56 #include "magick/color.h"
57 #include "magick/color-private.h"
58 #include "magick/colorspace.h"
59 #include "magick/colorspace-private.h"
60 #include "magick/composite.h"
61 #include "magick/composite-private.h"
62 #include "magick/constitute.h"
63 #include "magick/draw.h"
64 #include "magick/draw-private.h"
65 #include "magick/enhance.h"
66 #include "magick/exception.h"
67 #include "magick/exception-private.h"
68 #include "magick/gem.h"
69 #include "magick/geometry.h"
70 #include "magick/image-private.h"
71 #include "magick/list.h"
72 #include "magick/log.h"
73 #include "magick/magick.h"
74 #include "magick/memory-private.h"
75 #include "magick/monitor.h"
76 #include "magick/monitor-private.h"
77 #include "magick/option.h"
78 #include "magick/paint.h"
79 #include "magick/pixel-accessor.h"
80 #include "magick/pixel-private.h"
81 #include "magick/property.h"
82 #include "magick/resample.h"
83 #include "magick/resample-private.h"
84 #include "magick/resource_.h"
85 #include "magick/splay-tree.h"
86 #include "magick/string_.h"
87 #include "magick/string-private.h"
88 #include "magick/thread-private.h"
89 #include "magick/token.h"
90 #include "magick/transform.h"
91 #include "magick/utility.h"
92 
93 /*
94  Define declarations.
95 */
96 #define AntialiasThreshold (1.0/3.0)
97 #define BezierQuantum 200
98 #define PrimitiveExtentPad 4296.0
99 #define MaxBezierCoordinates 67108864
100 #define ThrowPointExpectedException(image,token) \
101 { \
102  (void) ThrowMagickException(&(image)->exception,GetMagickModule(),DrawError, \
103  "NonconformingDrawingPrimitiveDefinition","`%s'",token); \
104  status=MagickFalse; \
105  break; \
106 }
107 
108 /*
109  Typedef declarations.
110 */
111 typedef struct _EdgeInfo
112 {
114  bounds;
115 
116  double
117  scanline;
118 
119  PointInfo
120  *points;
121 
122  size_t
123  number_points;
124 
125  ssize_t
126  direction;
127 
128  MagickBooleanType
129  ghostline;
130 
131  size_t
132  highwater;
133 } EdgeInfo;
134 
135 typedef struct _ElementInfo
136 {
137  double
138  cx,
139  cy,
140  major,
141  minor,
142  angle;
143 } ElementInfo;
144 
145 typedef struct _MVGInfo
146 {
148  **primitive_info;
149 
150  size_t
151  *extent;
152 
153  ssize_t
154  offset;
155 
156  PointInfo
157  point;
158 
160  *exception;
161 } MVGInfo;
162 
163 typedef struct _PolygonInfo
164 {
165  EdgeInfo
166  *edges;
167 
168  size_t
169  number_edges;
170 } PolygonInfo;
171 
172 typedef enum
173 {
174  MoveToCode,
175  OpenCode,
176  GhostlineCode,
177  LineToCode,
178  EndCode
179 } PathInfoCode;
180 
181 typedef struct _PathInfo
182 {
183  PointInfo
184  point;
185 
186  PathInfoCode
187  code;
188 } PathInfo;
189 
190 /*
191  Forward declarations.
192 */
193 static Image
194  *DrawClippingMask(Image *,const DrawInfo *,const char *,const char *,
195  ExceptionInfo *);
196 
197 static MagickBooleanType
198  DrawStrokePolygon(Image *,const DrawInfo *,const PrimitiveInfo *),
199  RenderMVGContent(Image *,const DrawInfo *,const size_t),
200  TraceArc(MVGInfo *,const PointInfo,const PointInfo,const PointInfo),
201  TraceArcPath(MVGInfo *,const PointInfo,const PointInfo,const PointInfo,
202  const double,const MagickBooleanType,const MagickBooleanType),
203  TraceBezier(MVGInfo *,const size_t),
204  TraceCircle(MVGInfo *,const PointInfo,const PointInfo),
205  TraceEllipse(MVGInfo *,const PointInfo,const PointInfo,const PointInfo),
206  TraceLine(PrimitiveInfo *,const PointInfo,const PointInfo),
207  TraceRectangle(PrimitiveInfo *,const PointInfo,const PointInfo),
208  TraceRoundRectangle(MVGInfo *,const PointInfo,const PointInfo,PointInfo),
209  TraceSquareLinecap(PrimitiveInfo *,const size_t,const double);
210 
211 static PrimitiveInfo
212  *TraceStrokePolygon(const DrawInfo *,const PrimitiveInfo *,ExceptionInfo *);
213 
214 static ssize_t
215  TracePath(Image *,MVGInfo *,const char *);
216 
217 /*
218 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
219 % %
220 % %
221 % %
222 % A c q u i r e D r a w I n f o %
223 % %
224 % %
225 % %
226 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
227 %
228 % AcquireDrawInfo() returns a DrawInfo structure properly initialized.
229 %
230 % The format of the AcquireDrawInfo method is:
231 %
232 % DrawInfo *AcquireDrawInfo(void)
233 %
234 */
235 MagickExport DrawInfo *AcquireDrawInfo(void)
236 {
237  DrawInfo
238  *draw_info;
239 
240  draw_info=(DrawInfo *) AcquireCriticalMemory(sizeof(*draw_info));
241  GetDrawInfo((ImageInfo *) NULL,draw_info);
242  return(draw_info);
243 }
244 
245 /*
246 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
247 % %
248 % %
249 % %
250 % C l o n e D r a w I n f o %
251 % %
252 % %
253 % %
254 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
255 %
256 % CloneDrawInfo() makes a copy of the given draw_info structure. If NULL
257 % is specified, a new DrawInfo structure is created initialized to default
258 % values.
259 %
260 % The format of the CloneDrawInfo method is:
261 %
262 % DrawInfo *CloneDrawInfo(const ImageInfo *image_info,
263 % const DrawInfo *draw_info)
264 %
265 % A description of each parameter follows:
266 %
267 % o image_info: the image info.
268 %
269 % o draw_info: the draw info.
270 %
271 */
272 MagickExport DrawInfo *CloneDrawInfo(const ImageInfo *image_info,
273  const DrawInfo *draw_info)
274 {
275  DrawInfo
276  *clone_info;
277 
278  clone_info=(DrawInfo *) AcquireCriticalMemory(sizeof(*clone_info));
279  GetDrawInfo(image_info,clone_info);
280  if (draw_info == (DrawInfo *) NULL)
281  return(clone_info);
282  if (draw_info->id != (char *) NULL)
283  (void) CloneString(&clone_info->id,draw_info->id);
284  if (draw_info->primitive != (char *) NULL)
285  (void) CloneString(&clone_info->primitive,draw_info->primitive);
286  if (draw_info->geometry != (char *) NULL)
287  (void) CloneString(&clone_info->geometry,draw_info->geometry);
288  clone_info->compliance=draw_info->compliance;
289  clone_info->viewbox=draw_info->viewbox;
290  clone_info->affine=draw_info->affine;
291  clone_info->gravity=draw_info->gravity;
292  clone_info->fill=draw_info->fill;
293  clone_info->stroke=draw_info->stroke;
294  clone_info->stroke_width=draw_info->stroke_width;
295  if (draw_info->fill_pattern != (Image *) NULL)
296  clone_info->fill_pattern=CloneImage(draw_info->fill_pattern,0,0,MagickTrue,
297  &draw_info->fill_pattern->exception);
298  else
299  if (draw_info->tile != (Image *) NULL)
300  clone_info->fill_pattern=CloneImage(draw_info->tile,0,0,MagickTrue,
301  &draw_info->tile->exception);
302  clone_info->tile=NewImageList(); /* tile is deprecated */
303  if (draw_info->stroke_pattern != (Image *) NULL)
304  clone_info->stroke_pattern=CloneImage(draw_info->stroke_pattern,0,0,
305  MagickTrue,&draw_info->stroke_pattern->exception);
306  clone_info->stroke_antialias=draw_info->stroke_antialias;
307  clone_info->text_antialias=draw_info->text_antialias;
308  clone_info->fill_rule=draw_info->fill_rule;
309  clone_info->linecap=draw_info->linecap;
310  clone_info->linejoin=draw_info->linejoin;
311  clone_info->miterlimit=draw_info->miterlimit;
312  clone_info->dash_offset=draw_info->dash_offset;
313  clone_info->decorate=draw_info->decorate;
314  clone_info->compose=draw_info->compose;
315  if (draw_info->text != (char *) NULL)
316  (void) CloneString(&clone_info->text,draw_info->text);
317  if (draw_info->font != (char *) NULL)
318  (void) CloneString(&clone_info->font,draw_info->font);
319  if (draw_info->metrics != (char *) NULL)
320  (void) CloneString(&clone_info->metrics,draw_info->metrics);
321  if (draw_info->family != (char *) NULL)
322  (void) CloneString(&clone_info->family,draw_info->family);
323  clone_info->style=draw_info->style;
324  clone_info->stretch=draw_info->stretch;
325  clone_info->weight=draw_info->weight;
326  if (draw_info->encoding != (char *) NULL)
327  (void) CloneString(&clone_info->encoding,draw_info->encoding);
328  clone_info->pointsize=draw_info->pointsize;
329  clone_info->kerning=draw_info->kerning;
330  clone_info->interline_spacing=draw_info->interline_spacing;
331  clone_info->interword_spacing=draw_info->interword_spacing;
332  clone_info->direction=draw_info->direction;
333  if (draw_info->density != (char *) NULL)
334  (void) CloneString(&clone_info->density,draw_info->density);
335  clone_info->align=draw_info->align;
336  clone_info->undercolor=draw_info->undercolor;
337  clone_info->border_color=draw_info->border_color;
338  if (draw_info->server_name != (char *) NULL)
339  (void) CloneString(&clone_info->server_name,draw_info->server_name);
340  if (draw_info->dash_pattern != (double *) NULL)
341  {
342  ssize_t
343  x;
344 
345  for (x=0; fabs(draw_info->dash_pattern[x]) >= MagickEpsilon; x++) ;
346  clone_info->dash_pattern=(double *) AcquireQuantumMemory((size_t) (2*x+2),
347  sizeof(*clone_info->dash_pattern));
348  if (clone_info->dash_pattern == (double *) NULL)
349  ThrowFatalException(ResourceLimitFatalError,
350  "UnableToAllocateDashPattern");
351  (void) memset(clone_info->dash_pattern,0,(size_t) (2*x+2)*
352  sizeof(*clone_info->dash_pattern));
353  (void) memcpy(clone_info->dash_pattern,draw_info->dash_pattern,(size_t)
354  (x+1)*sizeof(*clone_info->dash_pattern));
355  }
356  clone_info->gradient=draw_info->gradient;
357  if (draw_info->gradient.stops != (StopInfo *) NULL)
358  {
359  size_t
360  number_stops;
361 
362  number_stops=clone_info->gradient.number_stops;
363  clone_info->gradient.stops=(StopInfo *) AcquireQuantumMemory((size_t)
364  number_stops,sizeof(*clone_info->gradient.stops));
365  if (clone_info->gradient.stops == (StopInfo *) NULL)
366  ThrowFatalException(ResourceLimitFatalError,
367  "UnableToAllocateDashPattern");
368  (void) memcpy(clone_info->gradient.stops,draw_info->gradient.stops,
369  (size_t) number_stops*sizeof(*clone_info->gradient.stops));
370  }
371  clone_info->bounds=draw_info->bounds;
372  clone_info->fill_opacity=draw_info->fill_opacity;
373  clone_info->stroke_opacity=draw_info->stroke_opacity;
374  clone_info->element_reference=draw_info->element_reference;
375  clone_info->clip_path=draw_info->clip_path;
376  clone_info->clip_units=draw_info->clip_units;
377  if (draw_info->clip_mask != (char *) NULL)
378  (void) CloneString(&clone_info->clip_mask,draw_info->clip_mask);
379  if (draw_info->clipping_mask != (Image *) NULL)
380  clone_info->clipping_mask=CloneImage(draw_info->clipping_mask,0,0,
381  MagickTrue,&draw_info->clipping_mask->exception);
382  if (draw_info->composite_mask != (Image *) NULL)
383  clone_info->composite_mask=CloneImage(draw_info->composite_mask,0,0,
384  MagickTrue,&draw_info->composite_mask->exception);
385  clone_info->render=draw_info->render;
386  clone_info->debug=draw_info->debug;
387  return(clone_info);
388 }
389 
390 /*
391 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
392 % %
393 % %
394 % %
395 + C o n v e r t P a t h T o P o l y g o n %
396 % %
397 % %
398 % %
399 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
400 %
401 % ConvertPathToPolygon() converts a path to the more efficient sorted
402 % rendering form.
403 %
404 % The format of the ConvertPathToPolygon method is:
405 %
406 % PolygonInfo *ConvertPathToPolygon(const PathInfo *path_info,
407 % ExceptionInfo *exception)
408 %
409 % A description of each parameter follows:
410 %
411 % o ConvertPathToPolygon() returns the path in a more efficient sorted
412 % rendering form of type PolygonInfo.
413 %
414 % o draw_info: Specifies a pointer to an DrawInfo structure.
415 %
416 % o path_info: Specifies a pointer to an PathInfo structure.
417 %
418 % o exception: return any errors or warnings in this structure.
419 %
420 */
421 
422 static PolygonInfo *DestroyPolygonInfo(PolygonInfo *polygon_info)
423 {
424  ssize_t
425  i;
426 
427  if (polygon_info->edges != (EdgeInfo *) NULL)
428  {
429  for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
430  if (polygon_info->edges[i].points != (PointInfo *) NULL)
431  polygon_info->edges[i].points=(PointInfo *)
432  RelinquishMagickMemory(polygon_info->edges[i].points);
433  polygon_info->edges=(EdgeInfo *) RelinquishMagickMemory(
434  polygon_info->edges);
435  }
436  return((PolygonInfo *) RelinquishMagickMemory(polygon_info));
437 }
438 
439 #if defined(__cplusplus) || defined(c_plusplus)
440 extern "C" {
441 #endif
442 
443 static int DrawCompareEdges(const void *p_edge,const void *q_edge)
444 {
445 #define DrawCompareEdge(p,q) \
446 { \
447  if (((p)-(q)) < 0.0) \
448  return(-1); \
449  if (((p)-(q)) > 0.0) \
450  return(1); \
451 }
452 
453  const PointInfo
454  *p,
455  *q;
456 
457  /*
458  Edge sorting for right-handed coordinate system.
459  */
460  p=((const EdgeInfo *) p_edge)->points;
461  q=((const EdgeInfo *) q_edge)->points;
462  DrawCompareEdge(p[0].y,q[0].y);
463  DrawCompareEdge(p[0].x,q[0].x);
464  DrawCompareEdge((p[1].x-p[0].x)*(q[1].y-q[0].y),(p[1].y-p[0].y)*
465  (q[1].x-q[0].x));
466  DrawCompareEdge(p[1].y,q[1].y);
467  DrawCompareEdge(p[1].x,q[1].x);
468  return(0);
469 }
470 
471 #if defined(__cplusplus) || defined(c_plusplus)
472 }
473 #endif
474 
475 static void LogPolygonInfo(const PolygonInfo *polygon_info)
476 {
477  EdgeInfo
478  *p;
479 
480  ssize_t
481  i,
482  j;
483 
484  (void) LogMagickEvent(DrawEvent,GetMagickModule()," begin active-edge");
485  p=polygon_info->edges;
486  for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
487  {
488  (void) LogMagickEvent(DrawEvent,GetMagickModule()," edge %.20g:",
489  (double) i);
490  (void) LogMagickEvent(DrawEvent,GetMagickModule()," direction: %s",
491  p->direction != MagickFalse ? "down" : "up");
492  (void) LogMagickEvent(DrawEvent,GetMagickModule()," ghostline: %s",
493  p->ghostline != MagickFalse ? "transparent" : "opaque");
494  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
495  " bounds: %g,%g - %g,%g",p->bounds.x1,p->bounds.y1,
496  p->bounds.x2,p->bounds.y2);
497  for (j=0; j < (ssize_t) p->number_points; j++)
498  (void) LogMagickEvent(DrawEvent,GetMagickModule()," %g,%g",
499  p->points[j].x,p->points[j].y);
500  p++;
501  }
502  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end active-edge");
503 }
504 
505 static void ReversePoints(PointInfo *points,const size_t number_points)
506 {
507  PointInfo
508  point;
509 
510  ssize_t
511  i;
512 
513  for (i=0; i < (ssize_t) (number_points >> 1); i++)
514  {
515  point=points[i];
516  points[i]=points[number_points-(i+1)];
517  points[number_points-(i+1)]=point;
518  }
519 }
520 
521 static PolygonInfo *ConvertPathToPolygon(const PathInfo *path_info,
522  ExceptionInfo *exception)
523 {
524  long
525  direction,
526  next_direction;
527 
528  PointInfo
529  point,
530  *points;
531 
533  *polygon_info;
534 
536  bounds;
537 
538  ssize_t
539  i,
540  n;
541 
542  MagickBooleanType
543  ghostline;
544 
545  size_t
546  edge,
547  number_edges,
548  number_points;
549 
550  /*
551  Convert a path to the more efficient sorted rendering form.
552  */
553  polygon_info=(PolygonInfo *) AcquireMagickMemory(sizeof(*polygon_info));
554  if (polygon_info == (PolygonInfo *) NULL)
555  {
556  (void) ThrowMagickException(exception,GetMagickModule(),
557  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
558  return((PolygonInfo *) NULL);
559  }
560  number_edges=16;
561  polygon_info->edges=(EdgeInfo *) AcquireQuantumMemory(number_edges,
562  sizeof(*polygon_info->edges));
563  if (polygon_info->edges == (EdgeInfo *) NULL)
564  {
565  (void) ThrowMagickException(exception,GetMagickModule(),
566  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
567  return(DestroyPolygonInfo(polygon_info));
568  }
569  (void) memset(polygon_info->edges,0,number_edges*
570  sizeof(*polygon_info->edges));
571  direction=0;
572  edge=0;
573  ghostline=MagickFalse;
574  n=0;
575  number_points=0;
576  points=(PointInfo *) NULL;
577  (void) memset(&point,0,sizeof(point));
578  (void) memset(&bounds,0,sizeof(bounds));
579  polygon_info->edges[edge].number_points=(size_t) n;
580  polygon_info->edges[edge].scanline=0.0;
581  polygon_info->edges[edge].highwater=0;
582  polygon_info->edges[edge].ghostline=ghostline;
583  polygon_info->edges[edge].direction=(ssize_t) direction;
584  polygon_info->edges[edge].points=points;
585  polygon_info->edges[edge].bounds=bounds;
586  polygon_info->number_edges=0;
587  for (i=0; path_info[i].code != EndCode; i++)
588  {
589  if ((path_info[i].code == MoveToCode) || (path_info[i].code == OpenCode) ||
590  (path_info[i].code == GhostlineCode))
591  {
592  /*
593  Move to.
594  */
595  if ((points != (PointInfo *) NULL) && (n >= 2))
596  {
597  if (edge == number_edges)
598  {
599  number_edges<<=1;
600  polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
601  polygon_info->edges,(size_t) number_edges,
602  sizeof(*polygon_info->edges));
603  if (polygon_info->edges == (EdgeInfo *) NULL)
604  {
605  (void) ThrowMagickException(exception,GetMagickModule(),
606  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
607  points=(PointInfo *) RelinquishMagickMemory(points);
608  return(DestroyPolygonInfo(polygon_info));
609  }
610  }
611  polygon_info->edges[edge].number_points=(size_t) n;
612  polygon_info->edges[edge].scanline=(-1.0);
613  polygon_info->edges[edge].highwater=0;
614  polygon_info->edges[edge].ghostline=ghostline;
615  polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
616  if (direction < 0)
617  ReversePoints(points,(size_t) n);
618  polygon_info->edges[edge].points=points;
619  polygon_info->edges[edge].bounds=bounds;
620  polygon_info->edges[edge].bounds.y1=points[0].y;
621  polygon_info->edges[edge].bounds.y2=points[n-1].y;
622  points=(PointInfo *) NULL;
623  ghostline=MagickFalse;
624  edge++;
625  polygon_info->number_edges=edge;
626  }
627  if (points == (PointInfo *) NULL)
628  {
629  number_points=16;
630  points=(PointInfo *) AcquireQuantumMemory((size_t) number_points,
631  sizeof(*points));
632  if (points == (PointInfo *) NULL)
633  {
634  (void) ThrowMagickException(exception,GetMagickModule(),
635  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
636  return(DestroyPolygonInfo(polygon_info));
637  }
638  }
639  ghostline=path_info[i].code == GhostlineCode ? MagickTrue : MagickFalse;
640  point=path_info[i].point;
641  points[0]=point;
642  bounds.x1=point.x;
643  bounds.x2=point.x;
644  direction=0;
645  n=1;
646  continue;
647  }
648  /*
649  Line to.
650  */
651  next_direction=((path_info[i].point.y > point.y) ||
652  ((fabs(path_info[i].point.y-point.y) < MagickEpsilon) &&
653  (path_info[i].point.x > point.x))) ? 1 : -1;
654  if ((points != (PointInfo *) NULL) && (direction != 0) &&
655  (direction != next_direction))
656  {
657  /*
658  New edge.
659  */
660  point=points[n-1];
661  if (edge == number_edges)
662  {
663  number_edges<<=1;
664  polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
665  polygon_info->edges,(size_t) number_edges,
666  sizeof(*polygon_info->edges));
667  if (polygon_info->edges == (EdgeInfo *) NULL)
668  {
669  (void) ThrowMagickException(exception,GetMagickModule(),
670  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
671  points=(PointInfo *) RelinquishMagickMemory(points);
672  return(DestroyPolygonInfo(polygon_info));
673  }
674  }
675  polygon_info->edges[edge].number_points=(size_t) n;
676  polygon_info->edges[edge].scanline=(-1.0);
677  polygon_info->edges[edge].highwater=0;
678  polygon_info->edges[edge].ghostline=ghostline;
679  polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
680  if (direction < 0)
681  ReversePoints(points,(size_t) n);
682  polygon_info->edges[edge].points=points;
683  polygon_info->edges[edge].bounds=bounds;
684  polygon_info->edges[edge].bounds.y1=points[0].y;
685  polygon_info->edges[edge].bounds.y2=points[n-1].y;
686  polygon_info->number_edges=edge+1;
687  points=(PointInfo *) NULL;
688  number_points=16;
689  points=(PointInfo *) AcquireQuantumMemory((size_t) number_points,
690  sizeof(*points));
691  if (points == (PointInfo *) NULL)
692  {
693  (void) ThrowMagickException(exception,GetMagickModule(),
694  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
695  return(DestroyPolygonInfo(polygon_info));
696  }
697  n=1;
698  ghostline=MagickFalse;
699  points[0]=point;
700  bounds.x1=point.x;
701  bounds.x2=point.x;
702  edge++;
703  }
704  direction=next_direction;
705  if (points == (PointInfo *) NULL)
706  continue;
707  if (n == (ssize_t) number_points)
708  {
709  number_points<<=1;
710  points=(PointInfo *) ResizeQuantumMemory(points,(size_t) number_points,
711  sizeof(*points));
712  if (points == (PointInfo *) NULL)
713  {
714  (void) ThrowMagickException(exception,GetMagickModule(),
715  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
716  return(DestroyPolygonInfo(polygon_info));
717  }
718  }
719  point=path_info[i].point;
720  points[n]=point;
721  if (point.x < bounds.x1)
722  bounds.x1=point.x;
723  if (point.x > bounds.x2)
724  bounds.x2=point.x;
725  n++;
726  }
727  if (points != (PointInfo *) NULL)
728  {
729  if (n < 2)
730  points=(PointInfo *) RelinquishMagickMemory(points);
731  else
732  {
733  if (edge == number_edges)
734  {
735  number_edges<<=1;
736  polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
737  polygon_info->edges,(size_t) number_edges,
738  sizeof(*polygon_info->edges));
739  if (polygon_info->edges == (EdgeInfo *) NULL)
740  {
741  (void) ThrowMagickException(exception,GetMagickModule(),
742  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
743  return(DestroyPolygonInfo(polygon_info));
744  }
745  }
746  polygon_info->edges[edge].number_points=(size_t) n;
747  polygon_info->edges[edge].scanline=(-1.0);
748  polygon_info->edges[edge].highwater=0;
749  polygon_info->edges[edge].ghostline=ghostline;
750  polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
751  if (direction < 0)
752  ReversePoints(points,(size_t) n);
753  polygon_info->edges[edge].points=points;
754  polygon_info->edges[edge].bounds=bounds;
755  polygon_info->edges[edge].bounds.y1=points[0].y;
756  polygon_info->edges[edge].bounds.y2=points[n-1].y;
757  points=(PointInfo *) NULL;
758  ghostline=MagickFalse;
759  edge++;
760  polygon_info->number_edges=edge;
761  }
762  }
763  polygon_info->number_edges=edge;
764  polygon_info->number_edges=edge;
765  polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(polygon_info->edges,
766  polygon_info->number_edges,sizeof(*polygon_info->edges));
767  if (polygon_info->edges == (EdgeInfo *) NULL)
768  {
769  (void) ThrowMagickException(exception,GetMagickModule(),
770  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
771  return(DestroyPolygonInfo(polygon_info));
772  }
773  for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
774  {
775  EdgeInfo
776  *edge_info;
777 
778  edge_info=polygon_info->edges+i;
779  edge_info->points=(PointInfo *) ResizeQuantumMemory(edge_info->points,
780  edge_info->number_points,sizeof(*edge_info->points));
781  if (edge_info->points == (PointInfo *) NULL)
782  {
783  (void) ThrowMagickException(exception,GetMagickModule(),
784  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
785  return(DestroyPolygonInfo(polygon_info));
786  }
787  }
788  qsort(polygon_info->edges,(size_t) polygon_info->number_edges,
789  sizeof(*polygon_info->edges),DrawCompareEdges);
790  if ((GetLogEventMask() & DrawEvent) != 0)
791  LogPolygonInfo(polygon_info);
792  return(polygon_info);
793 }
794 
795 /*
796 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
797 % %
798 % %
799 % %
800 + C o n v e r t P r i m i t i v e T o P a t h %
801 % %
802 % %
803 % %
804 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
805 %
806 % ConvertPrimitiveToPath() converts a PrimitiveInfo structure into a vector
807 % path structure.
808 %
809 % The format of the ConvertPrimitiveToPath method is:
810 %
811 % PathInfo *ConvertPrimitiveToPath(const DrawInfo *draw_info,
812 % const PrimitiveInfo *primitive_info,ExceptionInfo *exception)
813 %
814 % A description of each parameter follows:
815 %
816 % o ConvertPrimitiveToPath() returns a vector path structure of type
817 % PathInfo.
818 %
819 % o draw_info: a structure of type DrawInfo.
820 %
821 % o primitive_info: Specifies a pointer to an PrimitiveInfo structure.
822 %
823 %
824 */
825 
826 static void LogPathInfo(const PathInfo *path_info)
827 {
828  const PathInfo
829  *p;
830 
831  (void) LogMagickEvent(DrawEvent,GetMagickModule()," begin vector-path");
832  for (p=path_info; p->code != EndCode; p++)
833  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
834  " %g,%g %s",p->point.x,p->point.y,p->code == GhostlineCode ?
835  "moveto ghostline" : p->code == OpenCode ? "moveto open" :
836  p->code == MoveToCode ? "moveto" : p->code == LineToCode ? "lineto" :
837  "?");
838  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end vector-path");
839 }
840 
841 static PathInfo *ConvertPrimitiveToPath(
842  const DrawInfo *magick_unused(draw_info),const PrimitiveInfo *primitive_info,
843  ExceptionInfo *exception)
844 {
845  MagickBooleanType
846  closed_subpath;
847 
848  PathInfo
849  *path_info;
850 
851  PathInfoCode
852  code;
853 
854  PointInfo
855  p,
856  q;
857 
858  ssize_t
859  i,
860  n;
861 
862  ssize_t
863  coordinates,
864  start;
865 
866  magick_unreferenced(draw_info);
867 
868  /*
869  Converts a PrimitiveInfo structure into a vector path structure.
870  */
871  switch (primitive_info->primitive)
872  {
873  case PointPrimitive:
874  case ColorPrimitive:
875  case MattePrimitive:
876  case TextPrimitive:
877  case ImagePrimitive:
878  return((PathInfo *) NULL);
879  default:
880  break;
881  }
882  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
883  path_info=(PathInfo *) AcquireQuantumMemory((size_t) (3UL*i+1UL),
884  sizeof(*path_info));
885  if (path_info == (PathInfo *) NULL)
886  {
887  (void) ThrowMagickException(exception,GetMagickModule(),
888  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
889  return((PathInfo *) NULL);
890  }
891  coordinates=0;
892  closed_subpath=MagickFalse;
893  n=0;
894  p.x=(-1.0);
895  p.y=(-1.0);
896  q.x=(-1.0);
897  q.y=(-1.0);
898  start=0;
899  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
900  {
901  code=LineToCode;
902  if (coordinates <= 0)
903  {
904  /*
905  New subpath.
906  */
907  coordinates=(ssize_t) primitive_info[i].coordinates;
908  p=primitive_info[i].point;
909  start=n;
910  code=MoveToCode;
911  closed_subpath=primitive_info[i].closed_subpath;
912  }
913  coordinates--;
914  if ((code == MoveToCode) || (coordinates <= 0) ||
915  (fabs(q.x-primitive_info[i].point.x) >= MagickEpsilon) ||
916  (fabs(q.y-primitive_info[i].point.y) >= MagickEpsilon))
917  {
918  /*
919  Eliminate duplicate points.
920  */
921  path_info[n].code=code;
922  path_info[n].point=primitive_info[i].point;
923  q=primitive_info[i].point;
924  n++;
925  }
926  if (coordinates > 0)
927  continue; /* next point in current subpath */
928  if (closed_subpath != MagickFalse)
929  {
930  closed_subpath=MagickFalse;
931  continue;
932  }
933  /*
934  Mark the p point as open if the subpath is not closed.
935  */
936  path_info[start].code=OpenCode;
937  path_info[n].code=GhostlineCode;
938  path_info[n].point=primitive_info[i].point;
939  n++;
940  path_info[n].code=LineToCode;
941  path_info[n].point=p;
942  n++;
943  }
944  path_info[n].code=EndCode;
945  path_info[n].point.x=0.0;
946  path_info[n].point.y=0.0;
947  if (IsEventLogging() != MagickFalse)
948  LogPathInfo(path_info);
949  path_info=(PathInfo *) ResizeQuantumMemory(path_info,(size_t) (n+1),
950  sizeof(*path_info));
951  return(path_info);
952 }
953 
954 /*
955 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
956 % %
957 % %
958 % %
959 % D e s t r o y D r a w I n f o %
960 % %
961 % %
962 % %
963 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
964 %
965 % DestroyDrawInfo() deallocates memory associated with an DrawInfo structure.
966 %
967 % The format of the DestroyDrawInfo method is:
968 %
969 % DrawInfo *DestroyDrawInfo(DrawInfo *draw_info)
970 %
971 % A description of each parameter follows:
972 %
973 % o draw_info: the draw info.
974 %
975 */
976 MagickExport DrawInfo *DestroyDrawInfo(DrawInfo *draw_info)
977 {
978  assert(draw_info != (DrawInfo *) NULL);
979  assert(draw_info->signature == MagickCoreSignature);
980  if (IsEventLogging() != MagickFalse)
981  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
982  if (draw_info->id != (char *) NULL)
983  draw_info->id=DestroyString(draw_info->id);
984  if (draw_info->primitive != (char *) NULL)
985  draw_info->primitive=DestroyString(draw_info->primitive);
986  if (draw_info->text != (char *) NULL)
987  draw_info->text=DestroyString(draw_info->text);
988  if (draw_info->geometry != (char *) NULL)
989  draw_info->geometry=DestroyString(draw_info->geometry);
990  if (draw_info->tile != (Image *) NULL)
991  draw_info->tile=DestroyImage(draw_info->tile);
992  if (draw_info->fill_pattern != (Image *) NULL)
993  draw_info->fill_pattern=DestroyImage(draw_info->fill_pattern);
994  if (draw_info->stroke_pattern != (Image *) NULL)
995  draw_info->stroke_pattern=DestroyImage(draw_info->stroke_pattern);
996  if (draw_info->font != (char *) NULL)
997  draw_info->font=DestroyString(draw_info->font);
998  if (draw_info->metrics != (char *) NULL)
999  draw_info->metrics=DestroyString(draw_info->metrics);
1000  if (draw_info->family != (char *) NULL)
1001  draw_info->family=DestroyString(draw_info->family);
1002  if (draw_info->encoding != (char *) NULL)
1003  draw_info->encoding=DestroyString(draw_info->encoding);
1004  if (draw_info->density != (char *) NULL)
1005  draw_info->density=DestroyString(draw_info->density);
1006  if (draw_info->server_name != (char *) NULL)
1007  draw_info->server_name=(char *)
1008  RelinquishMagickMemory(draw_info->server_name);
1009  if (draw_info->dash_pattern != (double *) NULL)
1010  draw_info->dash_pattern=(double *) RelinquishMagickMemory(
1011  draw_info->dash_pattern);
1012  if (draw_info->gradient.stops != (StopInfo *) NULL)
1013  draw_info->gradient.stops=(StopInfo *) RelinquishMagickMemory(
1014  draw_info->gradient.stops);
1015  if (draw_info->clip_mask != (char *) NULL)
1016  draw_info->clip_mask=DestroyString(draw_info->clip_mask);
1017  if (draw_info->clipping_mask != (Image *) NULL)
1018  draw_info->clipping_mask=DestroyImage(draw_info->clipping_mask);
1019  if (draw_info->composite_mask != (Image *) NULL)
1020  draw_info->composite_mask=DestroyImage(draw_info->composite_mask);
1021  if (draw_info->image_info != (ImageInfo *) NULL)
1022  draw_info->image_info=DestroyImageInfo(draw_info->image_info);
1023  draw_info->signature=(~MagickCoreSignature);
1024  draw_info=(DrawInfo *) RelinquishMagickMemory(draw_info);
1025  return(draw_info);
1026 }
1027 
1028 /*
1029 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1030 % %
1031 % %
1032 % %
1033 % D r a w A f f i n e I m a g e %
1034 % %
1035 % %
1036 % %
1037 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1038 %
1039 % DrawAffineImage() composites the source over the destination image as
1040 % dictated by the affine transform.
1041 %
1042 % The format of the DrawAffineImage method is:
1043 %
1044 % MagickBooleanType DrawAffineImage(Image *image,const Image *source,
1045 % const AffineMatrix *affine)
1046 %
1047 % A description of each parameter follows:
1048 %
1049 % o image: the image.
1050 %
1051 % o source: the source image.
1052 %
1053 % o affine: the affine transform.
1054 %
1055 */
1056 
1057 static SegmentInfo AffineEdge(const Image *image,const AffineMatrix *affine,
1058  const double y,const SegmentInfo *edge)
1059 {
1060  double
1061  intercept,
1062  z;
1063 
1064  double
1065  x;
1066 
1067  SegmentInfo
1068  inverse_edge;
1069 
1070  /*
1071  Determine left and right edges.
1072  */
1073  inverse_edge.x1=edge->x1;
1074  inverse_edge.y1=edge->y1;
1075  inverse_edge.x2=edge->x2;
1076  inverse_edge.y2=edge->y2;
1077  z=affine->ry*y+affine->tx;
1078  if (affine->sx >= MagickEpsilon)
1079  {
1080  intercept=(-z/affine->sx);
1081  x=intercept;
1082  if (x > inverse_edge.x1)
1083  inverse_edge.x1=x;
1084  intercept=(-z+(double) image->columns)/affine->sx;
1085  x=intercept;
1086  if (x < inverse_edge.x2)
1087  inverse_edge.x2=x;
1088  }
1089  else
1090  if (affine->sx < -MagickEpsilon)
1091  {
1092  intercept=(-z+(double) image->columns)/affine->sx;
1093  x=intercept;
1094  if (x > inverse_edge.x1)
1095  inverse_edge.x1=x;
1096  intercept=(-z/affine->sx);
1097  x=intercept;
1098  if (x < inverse_edge.x2)
1099  inverse_edge.x2=x;
1100  }
1101  else
1102  if ((z < 0.0) || ((size_t) floor(z+0.5) >= image->columns))
1103  {
1104  inverse_edge.x2=edge->x1;
1105  return(inverse_edge);
1106  }
1107  /*
1108  Determine top and bottom edges.
1109  */
1110  z=affine->sy*y+affine->ty;
1111  if (affine->rx >= MagickEpsilon)
1112  {
1113  intercept=(-z/affine->rx);
1114  x=intercept;
1115  if (x > inverse_edge.x1)
1116  inverse_edge.x1=x;
1117  intercept=(-z+(double) image->rows)/affine->rx;
1118  x=intercept;
1119  if (x < inverse_edge.x2)
1120  inverse_edge.x2=x;
1121  }
1122  else
1123  if (affine->rx < -MagickEpsilon)
1124  {
1125  intercept=(-z+(double) image->rows)/affine->rx;
1126  x=intercept;
1127  if (x > inverse_edge.x1)
1128  inverse_edge.x1=x;
1129  intercept=(-z/affine->rx);
1130  x=intercept;
1131  if (x < inverse_edge.x2)
1132  inverse_edge.x2=x;
1133  }
1134  else
1135  if ((z < 0.0) || ((size_t) floor(z+0.5) >= image->rows))
1136  {
1137  inverse_edge.x2=edge->x2;
1138  return(inverse_edge);
1139  }
1140  return(inverse_edge);
1141 }
1142 
1143 static AffineMatrix InverseAffineMatrix(const AffineMatrix *affine)
1144 {
1145  AffineMatrix
1146  inverse_affine;
1147 
1148  double
1149  determinant;
1150 
1151  determinant=PerceptibleReciprocal(affine->sx*affine->sy-affine->rx*
1152  affine->ry);
1153  inverse_affine.sx=determinant*affine->sy;
1154  inverse_affine.rx=determinant*(-affine->rx);
1155  inverse_affine.ry=determinant*(-affine->ry);
1156  inverse_affine.sy=determinant*affine->sx;
1157  inverse_affine.tx=(-affine->tx)*inverse_affine.sx-affine->ty*
1158  inverse_affine.ry;
1159  inverse_affine.ty=(-affine->tx)*inverse_affine.rx-affine->ty*
1160  inverse_affine.sy;
1161  return(inverse_affine);
1162 }
1163 
1164 MagickExport MagickBooleanType DrawAffineImage(Image *image,
1165  const Image *source,const AffineMatrix *affine)
1166 {
1167  AffineMatrix
1168  inverse_affine;
1169 
1170  CacheView
1171  *image_view,
1172  *source_view;
1173 
1175  *exception;
1176 
1177  MagickBooleanType
1178  status;
1179 
1181  zero;
1182 
1183  PointInfo
1184  extent[4],
1185  min,
1186  max,
1187  point;
1188 
1189  ssize_t
1190  i;
1191 
1192  SegmentInfo
1193  edge;
1194 
1195  ssize_t
1196  start,
1197  stop,
1198  y;
1199 
1200  /*
1201  Determine bounding box.
1202  */
1203  assert(image != (Image *) NULL);
1204  assert(image->signature == MagickCoreSignature);
1205  assert(source != (const Image *) NULL);
1206  assert(source->signature == MagickCoreSignature);
1207  if (IsEventLogging() != MagickFalse)
1208  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1209  assert(affine != (AffineMatrix *) NULL);
1210  extent[0].x=0.0;
1211  extent[0].y=0.0;
1212  extent[1].x=(double) source->columns;
1213  extent[1].y=0.0;
1214  extent[2].x=(double) source->columns;
1215  extent[2].y=(double) source->rows;
1216  extent[3].x=0.0;
1217  extent[3].y=(double) source->rows;
1218  for (i=0; i < 4; i++)
1219  {
1220  point=extent[i];
1221  extent[i].x=point.x*affine->sx+point.y*affine->ry+affine->tx;
1222  extent[i].y=point.x*affine->rx+point.y*affine->sy+affine->ty;
1223  }
1224  min=extent[0];
1225  max=extent[0];
1226  for (i=1; i < 4; i++)
1227  {
1228  if (min.x > extent[i].x)
1229  min.x=extent[i].x;
1230  if (min.y > extent[i].y)
1231  min.y=extent[i].y;
1232  if (max.x < extent[i].x)
1233  max.x=extent[i].x;
1234  if (max.y < extent[i].y)
1235  max.y=extent[i].y;
1236  }
1237  /*
1238  Affine transform image.
1239  */
1240  if (SetImageStorageClass(image,DirectClass) == MagickFalse)
1241  return(MagickFalse);
1242  status=MagickTrue;
1243  edge.x1=min.x;
1244  edge.y1=min.y;
1245  edge.x2=max.x;
1246  edge.y2=max.y;
1247  inverse_affine=InverseAffineMatrix(affine);
1248  if (edge.y1 < 0.0)
1249  edge.y1=0.0;
1250  if (edge.y2 > (image->rows-1.0))
1251  edge.y2=image->rows-1.0;
1252  GetMagickPixelPacket(image,&zero);
1253  exception=(&image->exception);
1254  start=CastDoubleToLong(ceil(edge.y1-0.5));
1255  stop=CastDoubleToLong(floor(edge.y2+0.5));
1256  source_view=AcquireVirtualCacheView(source,exception);
1257  image_view=AcquireAuthenticCacheView(image,exception);
1258 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1259  #pragma omp parallel for schedule(static) shared(status) \
1260  magick_number_threads(source,image,stop-start,1)
1261 #endif
1262  for (y=start; y <= stop; y++)
1263  {
1264  IndexPacket
1265  *magick_restrict indexes;
1266 
1268  composite,
1269  pixel;
1270 
1271  PointInfo
1272  point;
1273 
1274  PixelPacket
1275  *magick_restrict q;
1276 
1277  SegmentInfo
1278  inverse_edge;
1279 
1280  ssize_t
1281  x,
1282  x_offset;
1283 
1284  if (status == MagickFalse)
1285  continue;
1286  inverse_edge=AffineEdge(source,&inverse_affine,(double) y,&edge);
1287  if (inverse_edge.x2 < inverse_edge.x1)
1288  continue;
1289  if (inverse_edge.x1 < 0.0)
1290  inverse_edge.x1=0.0;
1291  if (inverse_edge.x2 > image->columns-1.0)
1292  inverse_edge.x2=image->columns-1.0;
1293  q=GetCacheViewAuthenticPixels(image_view,CastDoubleToLong(
1294  ceil(inverse_edge.x1-0.5)),y,(size_t) CastDoubleToLong(floor(
1295  inverse_edge.x2+0.5)-ceil(inverse_edge.x1-0.5)+1),1,exception);
1296  if (q == (PixelPacket *) NULL)
1297  continue;
1298  indexes=GetCacheViewAuthenticIndexQueue(image_view);
1299  pixel=zero;
1300  composite=zero;
1301  x_offset=0;
1302  for (x=CastDoubleToLong(ceil(inverse_edge.x1-0.5));
1303  x <= CastDoubleToLong(floor(inverse_edge.x2+0.5)); x++)
1304  {
1305  point.x=(double) x*inverse_affine.sx+y*inverse_affine.ry+
1306  inverse_affine.tx;
1307  point.y=(double) x*inverse_affine.rx+y*inverse_affine.sy+
1308  inverse_affine.ty;
1309  status=InterpolateMagickPixelPacket(source,source_view,
1310  UndefinedInterpolatePixel,point.x,point.y,&pixel,exception);
1311  if (status == MagickFalse)
1312  break;
1313  SetMagickPixelPacket(image,q,indexes+x_offset,&composite);
1314  MagickPixelCompositeOver(&pixel,pixel.opacity,&composite,
1315  composite.opacity,&composite);
1316  SetPixelPacket(image,&composite,q,indexes+x_offset);
1317  x_offset++;
1318  q++;
1319  }
1320  if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1321  status=MagickFalse;
1322  }
1323  source_view=DestroyCacheView(source_view);
1324  image_view=DestroyCacheView(image_view);
1325  return(status);
1326 }
1327 
1328 /*
1329 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1330 % %
1331 % %
1332 % %
1333 + D r a w B o u n d i n g R e c t a n g l e s %
1334 % %
1335 % %
1336 % %
1337 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1338 %
1339 % DrawBoundingRectangles() draws the bounding rectangles on the image. This
1340 % is only useful for developers debugging the rendering algorithm.
1341 %
1342 % The format of the DrawBoundingRectangles method is:
1343 %
1344 % MagickBooleanType DrawBoundingRectangles(Image *image,
1345 % const DrawInfo *draw_info,PolygonInfo *polygon_info)
1346 %
1347 % A description of each parameter follows:
1348 %
1349 % o image: the image.
1350 %
1351 % o draw_info: the draw info.
1352 %
1353 % o polygon_info: Specifies a pointer to a PolygonInfo structure.
1354 %
1355 */
1356 
1357 static MagickBooleanType DrawBoundingRectangles(Image *image,
1358  const DrawInfo *draw_info,const PolygonInfo *polygon_info)
1359 {
1360  double
1361  mid;
1362 
1363  DrawInfo
1364  *clone_info;
1365 
1366  MagickStatusType
1367  status;
1368 
1369  PointInfo
1370  end,
1371  resolution,
1372  start;
1373 
1375  primitive_info[6];
1376 
1377  ssize_t
1378  i;
1379 
1380  SegmentInfo
1381  bounds;
1382 
1383  ssize_t
1384  coordinates;
1385 
1386  (void) memset(primitive_info,0,sizeof(primitive_info));
1387  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1388  status=QueryColorDatabase("#0000",&clone_info->fill,&image->exception);
1389  if (status == MagickFalse)
1390  {
1391  clone_info=DestroyDrawInfo(clone_info);
1392  return(MagickFalse);
1393  }
1394  resolution.x=96.0;
1395  resolution.y=96.0;
1396  if (clone_info->density != (char *) NULL)
1397  {
1398  GeometryInfo
1399  geometry_info;
1400 
1401  MagickStatusType
1402  flags;
1403 
1404  flags=ParseGeometry(clone_info->density,&geometry_info);
1405  if ((flags & RhoValue) != 0)
1406  resolution.x=geometry_info.rho;
1407  resolution.y=resolution.x;
1408  if ((flags & SigmaValue) != 0)
1409  resolution.y=geometry_info.sigma;
1410  }
1411  mid=(resolution.x/96.0)*ExpandAffine(&clone_info->affine)*
1412  clone_info->stroke_width/2.0;
1413  bounds.x1=0.0;
1414  bounds.y1=0.0;
1415  bounds.x2=0.0;
1416  bounds.y2=0.0;
1417  if (polygon_info != (PolygonInfo *) NULL)
1418  {
1419  bounds=polygon_info->edges[0].bounds;
1420  for (i=1; i < (ssize_t) polygon_info->number_edges; i++)
1421  {
1422  if (polygon_info->edges[i].bounds.x1 < (double) bounds.x1)
1423  bounds.x1=polygon_info->edges[i].bounds.x1;
1424  if (polygon_info->edges[i].bounds.y1 < (double) bounds.y1)
1425  bounds.y1=polygon_info->edges[i].bounds.y1;
1426  if (polygon_info->edges[i].bounds.x2 > (double) bounds.x2)
1427  bounds.x2=polygon_info->edges[i].bounds.x2;
1428  if (polygon_info->edges[i].bounds.y2 > (double) bounds.y2)
1429  bounds.y2=polygon_info->edges[i].bounds.y2;
1430  }
1431  bounds.x1-=mid;
1432  bounds.x1=bounds.x1 < 0.0 ? 0.0 : bounds.x1 >= (double)
1433  image->columns ? (double) image->columns-1 : bounds.x1;
1434  bounds.y1-=mid;
1435  bounds.y1=bounds.y1 < 0.0 ? 0.0 : bounds.y1 >= (double)
1436  image->rows ? (double) image->rows-1 : bounds.y1;
1437  bounds.x2+=mid;
1438  bounds.x2=bounds.x2 < 0.0 ? 0.0 : bounds.x2 >= (double)
1439  image->columns ? (double) image->columns-1 : bounds.x2;
1440  bounds.y2+=mid;
1441  bounds.y2=bounds.y2 < 0.0 ? 0.0 : bounds.y2 >= (double)
1442  image->rows ? (double) image->rows-1 : bounds.y2;
1443  for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
1444  {
1445  if (polygon_info->edges[i].direction != 0)
1446  status=QueryColorDatabase("#f00",&clone_info->stroke,
1447  &image->exception);
1448  else
1449  status=QueryColorDatabase("#0f0",&clone_info->stroke,
1450  &image->exception);
1451  if (status == MagickFalse)
1452  break;
1453  start.x=(double) (polygon_info->edges[i].bounds.x1-mid);
1454  start.y=(double) (polygon_info->edges[i].bounds.y1-mid);
1455  end.x=(double) (polygon_info->edges[i].bounds.x2+mid);
1456  end.y=(double) (polygon_info->edges[i].bounds.y2+mid);
1457  primitive_info[0].primitive=RectanglePrimitive;
1458  status&=TraceRectangle(primitive_info,start,end);
1459  primitive_info[0].method=ReplaceMethod;
1460  coordinates=(ssize_t) primitive_info[0].coordinates;
1461  primitive_info[coordinates].primitive=UndefinedPrimitive;
1462  status=DrawPrimitive(image,clone_info,primitive_info);
1463  if (status == MagickFalse)
1464  break;
1465  }
1466  if (i < (ssize_t) polygon_info->number_edges)
1467  {
1468  clone_info=DestroyDrawInfo(clone_info);
1469  return(status == 0 ? MagickFalse : MagickTrue);
1470  }
1471  }
1472  status=QueryColorDatabase("#00f",&clone_info->stroke,&image->exception);
1473  if (status == MagickFalse)
1474  {
1475  clone_info=DestroyDrawInfo(clone_info);
1476  return(MagickFalse);
1477  }
1478  start.x=(double) (bounds.x1-mid);
1479  start.y=(double) (bounds.y1-mid);
1480  end.x=(double) (bounds.x2+mid);
1481  end.y=(double) (bounds.y2+mid);
1482  primitive_info[0].primitive=RectanglePrimitive;
1483  status&=TraceRectangle(primitive_info,start,end);
1484  primitive_info[0].method=ReplaceMethod;
1485  coordinates=(ssize_t) primitive_info[0].coordinates;
1486  primitive_info[coordinates].primitive=UndefinedPrimitive;
1487  status=DrawPrimitive(image,clone_info,primitive_info);
1488  clone_info=DestroyDrawInfo(clone_info);
1489  return(status == 0 ? MagickFalse : MagickTrue);
1490 }
1491 
1492 /*
1493 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1494 % %
1495 % %
1496 % %
1497 % D r a w C l i p P a t h %
1498 % %
1499 % %
1500 % %
1501 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1502 %
1503 % DrawClipPath() draws the clip path on the image mask.
1504 %
1505 % The format of the DrawClipPath method is:
1506 %
1507 % MagickBooleanType DrawClipPath(Image *image,const DrawInfo *draw_info,
1508 % const char *id)
1509 %
1510 % A description of each parameter follows:
1511 %
1512 % o image: the image.
1513 %
1514 % o draw_info: the draw info.
1515 %
1516 % o id: the clip path id.
1517 %
1518 */
1519 MagickExport MagickBooleanType DrawClipPath(Image *image,
1520  const DrawInfo *draw_info,const char *id)
1521 {
1522  const char
1523  *clip_path;
1524 
1525  Image
1526  *clipping_mask;
1527 
1528  MagickBooleanType
1529  status;
1530 
1531  clip_path=GetImageArtifact(image,id);
1532  if (clip_path == (const char *) NULL)
1533  return(MagickFalse);
1534  clipping_mask=DrawClippingMask(image,draw_info,draw_info->clip_mask,clip_path,
1535  &image->exception);
1536  if (clipping_mask == (Image *) NULL)
1537  return(MagickFalse);
1538  status=SetImageClipMask(image,clipping_mask);
1539  clipping_mask=DestroyImage(clipping_mask);
1540  return(status);
1541 }
1542 
1543 /*
1544 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1545 % %
1546 % %
1547 % %
1548 % D r a w C l i p p i n g M a s k %
1549 % %
1550 % %
1551 % %
1552 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1553 %
1554 % DrawClippingMask() draws the clip path and returns it as an image clipping
1555 % mask.
1556 %
1557 % The format of the DrawClippingMask method is:
1558 %
1559 % Image *DrawClippingMask(Image *image,const DrawInfo *draw_info,
1560 % const char *id,const char *clip_path,ExceptionInfo *exception)
1561 %
1562 % A description of each parameter follows:
1563 %
1564 % o image: the image.
1565 %
1566 % o draw_info: the draw info.
1567 %
1568 % o id: the clip path id.
1569 %
1570 % o clip_path: the clip path.
1571 %
1572 % o exception: return any errors or warnings in this structure.
1573 %
1574 */
1575 static Image *DrawClippingMask(Image *image,const DrawInfo *draw_info,
1576  const char *id,const char *clip_path,ExceptionInfo *exception)
1577 {
1578  DrawInfo
1579  *clone_info;
1580 
1581  Image
1582  *clip_mask;
1583 
1584  MagickStatusType
1585  status;
1586 
1587  /*
1588  Draw a clip path.
1589  */
1590  assert(image != (Image *) NULL);
1591  assert(image->signature == MagickCoreSignature);
1592  assert(draw_info != (const DrawInfo *) NULL);
1593  if (IsEventLogging() != MagickFalse)
1594  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1595  clip_mask=AcquireImage((const ImageInfo *) NULL);
1596  status=SetImageExtent(clip_mask,image->columns,image->rows);
1597  if (status == MagickFalse)
1598  return(DestroyImage(clip_mask));
1599  status=SetImageClipMask(image,(Image *) NULL);
1600  status=QueryColorCompliance("#0000",AllCompliance,
1601  &clip_mask->background_color,exception);
1602  clip_mask->background_color.opacity=(Quantum) TransparentOpacity;
1603  status=SetImageBackgroundColor(clip_mask);
1604  if (draw_info->debug != MagickFalse)
1605  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"\nbegin clip-path %s",
1606  id);
1607  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1608  (void) CloneString(&clone_info->primitive,clip_path);
1609  status=QueryColorCompliance("#ffffff",AllCompliance,&clone_info->fill,
1610  exception);
1611  if (clone_info->clip_mask != (char *) NULL)
1612  clone_info->clip_mask=DestroyString(clone_info->clip_mask);
1613  (void) QueryColorCompliance("#00000000",AllCompliance,&clone_info->stroke,
1614  exception);
1615  clone_info->stroke_width=0.0;
1616  clone_info->opacity=OpaqueOpacity;
1617  clone_info->clip_path=MagickTrue;
1618  status=RenderMVGContent(clip_mask,clone_info,0);
1619  clone_info=DestroyDrawInfo(clone_info);
1620  status&=SeparateImageChannel(clip_mask,TrueAlphaChannel);
1621  if (draw_info->compliance != SVGCompliance)
1622  status&=NegateImage(clip_mask,MagickFalse);
1623  if (status == MagickFalse)
1624  clip_mask=DestroyImage(clip_mask);
1625  if (draw_info->debug != MagickFalse)
1626  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end clip-path");
1627  return(clip_mask);
1628 }
1629 
1630 /*
1631 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1632 % %
1633 % %
1634 % %
1635 % D r a w C o m p o s i t e M a s k %
1636 % %
1637 % %
1638 % %
1639 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1640 %
1641 % DrawCompositeMask() draws the mask path and returns it as an image mask.
1642 %
1643 % The format of the DrawCompositeMask method is:
1644 %
1645 % Image *DrawCompositeMask(Image *image,const DrawInfo *draw_info,
1646 % const char *id,const char *mask_path,ExceptionInfo *exception)
1647 %
1648 % A description of each parameter follows:
1649 %
1650 % o image: the image.
1651 %
1652 % o draw_info: the draw info.
1653 %
1654 % o id: the mask path id.
1655 %
1656 % o mask_path: the mask path.
1657 %
1658 % o exception: return any errors or warnings in this structure.
1659 %
1660 */
1661 static Image *DrawCompositeMask(Image *image,const DrawInfo *draw_info,
1662  const char *id,const char *mask_path,ExceptionInfo *exception)
1663 {
1664  Image
1665  *composite_mask;
1666 
1667  DrawInfo
1668  *clone_info;
1669 
1670  MagickStatusType
1671  status;
1672 
1673  /*
1674  Draw a mask path.
1675  */
1676  assert(image != (Image *) NULL);
1677  assert(image->signature == MagickCoreSignature);
1678  assert(draw_info != (const DrawInfo *) NULL);
1679  if (IsEventLogging() != MagickFalse)
1680  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1681  composite_mask=AcquireImage((const ImageInfo *) NULL);
1682  status=SetImageExtent(composite_mask,image->columns,image->rows);
1683  if (status == MagickFalse)
1684  return(DestroyImage(composite_mask));
1685  status=SetImageMask(image,(Image *) NULL);
1686  status=QueryColorCompliance("#0000",AllCompliance,
1687  &composite_mask->background_color,exception);
1688  composite_mask->background_color.opacity=(Quantum) TransparentOpacity;
1689  (void) SetImageBackgroundColor(composite_mask);
1690  if (draw_info->debug != MagickFalse)
1691  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"\nbegin mask-path %s",
1692  id);
1693  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1694  (void) CloneString(&clone_info->primitive,mask_path);
1695  status=QueryColorCompliance("#ffffff",AllCompliance,&clone_info->fill,
1696  exception);
1697  status=QueryColorCompliance("#00000000",AllCompliance,&clone_info->stroke,
1698  exception);
1699  clone_info->stroke_width=0.0;
1700  clone_info->opacity=OpaqueOpacity;
1701  status=RenderMVGContent(composite_mask,clone_info,0);
1702  clone_info=DestroyDrawInfo(clone_info);
1703  status&=SeparateImageChannel(composite_mask,TrueAlphaChannel);
1704  status&=NegateImage(composite_mask,MagickFalse);
1705  if (status == MagickFalse)
1706  composite_mask=DestroyImage(composite_mask);
1707  if (draw_info->debug != MagickFalse)
1708  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end mask-path");
1709  return(composite_mask);
1710 }
1711 
1712 /*
1713 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1714 % %
1715 % %
1716 % %
1717 + D r a w D a s h P o l y g o n %
1718 % %
1719 % %
1720 % %
1721 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1722 %
1723 % DrawDashPolygon() draws a dashed polygon (line, rectangle, ellipse) on the
1724 % image while respecting the dash offset and dash pattern attributes.
1725 %
1726 % The format of the DrawDashPolygon method is:
1727 %
1728 % MagickBooleanType DrawDashPolygon(const DrawInfo *draw_info,
1729 % const PrimitiveInfo *primitive_info,Image *image)
1730 %
1731 % A description of each parameter follows:
1732 %
1733 % o draw_info: the draw info.
1734 %
1735 % o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
1736 %
1737 % o image: the image.
1738 %
1739 %
1740 */
1741 static MagickBooleanType DrawDashPolygon(const DrawInfo *draw_info,
1742  const PrimitiveInfo *primitive_info,Image *image)
1743 {
1744  double
1745  dx,
1746  dy,
1747  length,
1748  maximum_length,
1749  offset,
1750  scale,
1751  total_length;
1752 
1753  DrawInfo
1754  *clone_info;
1755 
1756  MagickStatusType
1757  status;
1758 
1760  *dash_polygon;
1761 
1762  ssize_t
1763  i;
1764 
1765  size_t
1766  number_vertices;
1767 
1768  ssize_t
1769  j,
1770  n;
1771 
1772  assert(draw_info != (const DrawInfo *) NULL);
1773  if (draw_info->debug != MagickFalse)
1774  (void) LogMagickEvent(DrawEvent,GetMagickModule()," begin draw-dash");
1775  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
1776  number_vertices=(size_t) i;
1777  dash_polygon=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
1778  (2UL*number_vertices+32UL),sizeof(*dash_polygon));
1779  if (dash_polygon == (PrimitiveInfo *) NULL)
1780  {
1781  (void) ThrowMagickException(&image->exception,GetMagickModule(),
1782  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
1783  return(MagickFalse);
1784  }
1785  (void) memset(dash_polygon,0,(2UL*number_vertices+32UL)*
1786  sizeof(*dash_polygon));
1787  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
1788  clone_info->miterlimit=0;
1789  dash_polygon[0]=primitive_info[0];
1790  scale=ExpandAffine(&draw_info->affine);
1791  length=scale*draw_info->dash_pattern[0];
1792  offset=fabs(draw_info->dash_offset) >= MagickEpsilon ?
1793  scale*draw_info->dash_offset : 0.0;
1794  j=1;
1795  for (n=0; offset > 0.0; j=0)
1796  {
1797  if (draw_info->dash_pattern[n] <= 0.0)
1798  break;
1799  length=scale*(draw_info->dash_pattern[n]+(n == 0 ? -0.5 : 0.5));
1800  if (offset > length)
1801  {
1802  offset-=length;
1803  n++;
1804  length=scale*draw_info->dash_pattern[n];
1805  continue;
1806  }
1807  if (offset < length)
1808  {
1809  length-=offset;
1810  offset=0.0;
1811  break;
1812  }
1813  offset=0.0;
1814  n++;
1815  }
1816  status=MagickTrue;
1817  maximum_length=0.0;
1818  total_length=0.0;
1819  for (i=1; (i < (ssize_t) number_vertices) && (length >= 0.0); i++)
1820  {
1821  dx=primitive_info[i].point.x-primitive_info[i-1].point.x;
1822  dy=primitive_info[i].point.y-primitive_info[i-1].point.y;
1823  maximum_length=hypot(dx,dy);
1824  if (maximum_length > (double) (MaxBezierCoordinates >> 2))
1825  continue;
1826  if (fabs(length) < MagickEpsilon)
1827  {
1828  if (fabs(draw_info->dash_pattern[n]) >= MagickEpsilon)
1829  n++;
1830  if (fabs(draw_info->dash_pattern[n]) < MagickEpsilon)
1831  n=0;
1832  length=scale*draw_info->dash_pattern[n];
1833  }
1834  for (total_length=0.0; (length >= 0.0) && (maximum_length >= (total_length+length)); )
1835  {
1836  total_length+=length;
1837  if ((n & 0x01) != 0)
1838  {
1839  dash_polygon[0]=primitive_info[0];
1840  dash_polygon[0].point.x=(double) (primitive_info[i-1].point.x+dx*
1841  total_length*PerceptibleReciprocal(maximum_length));
1842  dash_polygon[0].point.y=(double) (primitive_info[i-1].point.y+dy*
1843  total_length*PerceptibleReciprocal(maximum_length));
1844  j=1;
1845  }
1846  else
1847  {
1848  if ((j+1) > (ssize_t) number_vertices)
1849  break;
1850  dash_polygon[j]=primitive_info[i-1];
1851  dash_polygon[j].point.x=(double) (primitive_info[i-1].point.x+dx*
1852  total_length*PerceptibleReciprocal(maximum_length));
1853  dash_polygon[j].point.y=(double) (primitive_info[i-1].point.y+dy*
1854  total_length*PerceptibleReciprocal(maximum_length));
1855  dash_polygon[j].coordinates=1;
1856  j++;
1857  dash_polygon[0].coordinates=(size_t) j;
1858  dash_polygon[j].primitive=UndefinedPrimitive;
1859  status&=DrawStrokePolygon(image,clone_info,dash_polygon);
1860  if (status == MagickFalse)
1861  break;
1862  }
1863  if (fabs(draw_info->dash_pattern[n]) >= MagickEpsilon)
1864  n++;
1865  if (fabs(draw_info->dash_pattern[n]) < MagickEpsilon)
1866  n=0;
1867  length=scale*draw_info->dash_pattern[n];
1868  }
1869  length-=(maximum_length-total_length);
1870  if ((n & 0x01) != 0)
1871  continue;
1872  dash_polygon[j]=primitive_info[i];
1873  dash_polygon[j].coordinates=1;
1874  j++;
1875  }
1876  if ((status != MagickFalse) && (total_length < maximum_length) &&
1877  ((n & 0x01) == 0) && (j > 1))
1878  {
1879  dash_polygon[j]=primitive_info[i-1];
1880  dash_polygon[j].point.x+=MagickEpsilon;
1881  dash_polygon[j].point.y+=MagickEpsilon;
1882  dash_polygon[j].coordinates=1;
1883  j++;
1884  dash_polygon[0].coordinates=(size_t) j;
1885  dash_polygon[j].primitive=UndefinedPrimitive;
1886  status&=DrawStrokePolygon(image,clone_info,dash_polygon);
1887  }
1888  dash_polygon=(PrimitiveInfo *) RelinquishMagickMemory(dash_polygon);
1889  clone_info=DestroyDrawInfo(clone_info);
1890  if (draw_info->debug != MagickFalse)
1891  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end draw-dash");
1892  return(status != 0 ? MagickTrue : MagickFalse);
1893 }
1894 
1895 /*
1896 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1897 % %
1898 % %
1899 % %
1900 % D r a w G r a d i e n t I m a g e %
1901 % %
1902 % %
1903 % %
1904 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1905 %
1906 % DrawGradientImage() draws a linear gradient on the image.
1907 %
1908 % The format of the DrawGradientImage method is:
1909 %
1910 % MagickBooleanType DrawGradientImage(Image *image,
1911 % const DrawInfo *draw_info)
1912 %
1913 % A description of each parameter follows:
1914 %
1915 % o image: the image.
1916 %
1917 % o draw_info: the draw info.
1918 %
1919 */
1920 
1921 static inline double GetStopColorOffset(const GradientInfo *gradient,
1922  const ssize_t x,const ssize_t y)
1923 {
1924  switch (gradient->type)
1925  {
1926  case UndefinedGradient:
1927  case LinearGradient:
1928  {
1929  double
1930  gamma,
1931  length,
1932  offset,
1933  scale;
1934 
1935  PointInfo
1936  p,
1937  q;
1938 
1939  const SegmentInfo
1940  *gradient_vector;
1941 
1942  gradient_vector=(&gradient->gradient_vector);
1943  p.x=gradient_vector->x2-gradient_vector->x1;
1944  p.y=gradient_vector->y2-gradient_vector->y1;
1945  q.x=(double) x-gradient_vector->x1;
1946  q.y=(double) y-gradient_vector->y1;
1947  length=sqrt(q.x*q.x+q.y*q.y);
1948  gamma=sqrt(p.x*p.x+p.y*p.y)*length;
1949  gamma=PerceptibleReciprocal(gamma);
1950  scale=p.x*q.x+p.y*q.y;
1951  offset=gamma*scale*length;
1952  return(offset);
1953  }
1954  case RadialGradient:
1955  {
1956  PointInfo
1957  v;
1958 
1959  if (gradient->spread == RepeatSpread)
1960  {
1961  v.x=(double) x-gradient->center.x;
1962  v.y=(double) y-gradient->center.y;
1963  return(sqrt(v.x*v.x+v.y*v.y));
1964  }
1965  v.x=(double) (((x-gradient->center.x)*cos(DegreesToRadians(
1966  gradient->angle)))+((y-gradient->center.y)*sin(DegreesToRadians(
1967  gradient->angle))))*PerceptibleReciprocal(gradient->radii.x);
1968  v.y=(double) (((x-gradient->center.x)*sin(DegreesToRadians(
1969  gradient->angle)))-((y-gradient->center.y)*cos(DegreesToRadians(
1970  gradient->angle))))*PerceptibleReciprocal(gradient->radii.y);
1971  return(sqrt(v.x*v.x+v.y*v.y));
1972  }
1973  }
1974  return(0.0);
1975 }
1976 
1977 MagickExport MagickBooleanType DrawGradientImage(Image *image,
1978  const DrawInfo *draw_info)
1979 {
1980  CacheView
1981  *image_view;
1982 
1983  const GradientInfo
1984  *gradient;
1985 
1986  const SegmentInfo
1987  *gradient_vector;
1988 
1989  double
1990  length;
1991 
1993  *exception;
1994 
1995  MagickBooleanType
1996  status;
1997 
1999  zero;
2000 
2001  PointInfo
2002  point;
2003 
2005  bounding_box;
2006 
2007  ssize_t
2008  y;
2009 
2010  /*
2011  Draw linear or radial gradient on image.
2012  */
2013  assert(image != (Image *) NULL);
2014  assert(image->signature == MagickCoreSignature);
2015  assert(draw_info != (const DrawInfo *) NULL);
2016  if (IsEventLogging() != MagickFalse)
2017  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
2018  gradient=(&draw_info->gradient);
2019  gradient_vector=(&gradient->gradient_vector);
2020  point.x=gradient_vector->x2-gradient_vector->x1;
2021  point.y=gradient_vector->y2-gradient_vector->y1;
2022  length=sqrt(point.x*point.x+point.y*point.y);
2023  bounding_box=gradient->bounding_box;
2024  status=MagickTrue;
2025  exception=(&image->exception);
2026  GetMagickPixelPacket(image,&zero);
2027  image_view=AcquireAuthenticCacheView(image,exception);
2028 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2029  #pragma omp parallel for schedule(static) shared(status) \
2030  magick_number_threads(image,image,bounding_box.height-bounding_box.y,1)
2031 #endif
2032  for (y=bounding_box.y; y < (ssize_t) bounding_box.height; y++)
2033  {
2034  double
2035  alpha,
2036  offset;
2037 
2039  composite,
2040  pixel;
2041 
2042  IndexPacket
2043  *magick_restrict indexes;
2044 
2045  ssize_t
2046  i,
2047  x;
2048 
2049  PixelPacket
2050  *magick_restrict q;
2051 
2052  ssize_t
2053  j;
2054 
2055  if (status == MagickFalse)
2056  continue;
2057  q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
2058  if (q == (PixelPacket *) NULL)
2059  {
2060  status=MagickFalse;
2061  continue;
2062  }
2063  indexes=GetCacheViewAuthenticIndexQueue(image_view);
2064  pixel=zero;
2065  composite=zero;
2066  offset=GetStopColorOffset(gradient,0,y);
2067  if (gradient->type != RadialGradient)
2068  offset*=PerceptibleReciprocal(length);
2069  for (x=bounding_box.x; x < (ssize_t) bounding_box.width; x++)
2070  {
2071  SetMagickPixelPacket(image,q,indexes+x,&pixel);
2072  switch (gradient->spread)
2073  {
2074  case UndefinedSpread:
2075  case PadSpread:
2076  {
2077  if ((x != CastDoubleToLong(ceil(gradient_vector->x1-0.5))) ||
2078  (y != CastDoubleToLong(ceil(gradient_vector->y1-0.5))))
2079  {
2080  offset=GetStopColorOffset(gradient,x,y);
2081  if (gradient->type != RadialGradient)
2082  offset*=PerceptibleReciprocal(length);
2083  }
2084  for (i=0; i < (ssize_t) gradient->number_stops; i++)
2085  if (offset < gradient->stops[i].offset)
2086  break;
2087  if ((offset < 0.0) || (i == 0))
2088  composite=gradient->stops[0].color;
2089  else
2090  if ((offset > 1.0) || (i == (ssize_t) gradient->number_stops))
2091  composite=gradient->stops[gradient->number_stops-1].color;
2092  else
2093  {
2094  j=i;
2095  i--;
2096  alpha=(offset-gradient->stops[i].offset)/
2097  (gradient->stops[j].offset-gradient->stops[i].offset);
2098  MagickPixelCompositeBlend(&gradient->stops[i].color,1.0-alpha,
2099  &gradient->stops[j].color,alpha,&composite);
2100  }
2101  break;
2102  }
2103  case ReflectSpread:
2104  {
2105  if ((x != CastDoubleToLong(ceil(gradient_vector->x1-0.5))) ||
2106  (y != CastDoubleToLong(ceil(gradient_vector->y1-0.5))))
2107  {
2108  offset=GetStopColorOffset(gradient,x,y);
2109  if (gradient->type != RadialGradient)
2110  offset*=PerceptibleReciprocal(length);
2111  }
2112  if (offset < 0.0)
2113  offset=(-offset);
2114  if ((ssize_t) fmod(offset,2.0) == 0)
2115  offset=fmod(offset,1.0);
2116  else
2117  offset=1.0-fmod(offset,1.0);
2118  for (i=0; i < (ssize_t) gradient->number_stops; i++)
2119  if (offset < gradient->stops[i].offset)
2120  break;
2121  if (i == 0)
2122  composite=gradient->stops[0].color;
2123  else
2124  if (i == (ssize_t) gradient->number_stops)
2125  composite=gradient->stops[gradient->number_stops-1].color;
2126  else
2127  {
2128  j=i;
2129  i--;
2130  alpha=(offset-gradient->stops[i].offset)/
2131  (gradient->stops[j].offset-gradient->stops[i].offset);
2132  MagickPixelCompositeBlend(&gradient->stops[i].color,1.0-alpha,
2133  &gradient->stops[j].color,alpha,&composite);
2134  }
2135  break;
2136  }
2137  case RepeatSpread:
2138  {
2139  double
2140  repeat;
2141 
2142  MagickBooleanType
2143  antialias;
2144 
2145  antialias=MagickFalse;
2146  repeat=0.0;
2147  if ((x != CastDoubleToLong(ceil(gradient_vector->x1-0.5))) ||
2148  (y != CastDoubleToLong(ceil(gradient_vector->y1-0.5))))
2149  {
2150  offset=GetStopColorOffset(gradient,x,y);
2151  if (gradient->type == LinearGradient)
2152  {
2153  repeat=fmod(offset,length);
2154  if (repeat < 0.0)
2155  repeat=length-fmod(-repeat,length);
2156  else
2157  repeat=fmod(offset,length);
2158  antialias=(repeat < length) && ((repeat+1.0) > length) ?
2159  MagickTrue : MagickFalse;
2160  offset=PerceptibleReciprocal(length)*repeat;
2161  }
2162  else
2163  {
2164  repeat=fmod(offset,(double) gradient->radius);
2165  if (repeat < 0.0)
2166  repeat=gradient->radius-fmod(-repeat,
2167  (double) gradient->radius);
2168  else
2169  repeat=fmod(offset,(double) gradient->radius);
2170  antialias=repeat+1.0 > gradient->radius ? MagickTrue :
2171  MagickFalse;
2172  offset=repeat*PerceptibleReciprocal(gradient->radius);
2173  }
2174  }
2175  for (i=0; i < (ssize_t) gradient->number_stops; i++)
2176  if (offset < gradient->stops[i].offset)
2177  break;
2178  if (i == 0)
2179  composite=gradient->stops[0].color;
2180  else
2181  if (i == (ssize_t) gradient->number_stops)
2182  composite=gradient->stops[gradient->number_stops-1].color;
2183  else
2184  {
2185  j=i;
2186  i--;
2187  alpha=(offset-gradient->stops[i].offset)/
2188  (gradient->stops[j].offset-gradient->stops[i].offset);
2189  if (antialias != MagickFalse)
2190  {
2191  if (gradient->type == LinearGradient)
2192  alpha=length-repeat;
2193  else
2194  alpha=gradient->radius-repeat;
2195  i=0;
2196  j=(ssize_t) gradient->number_stops-1L;
2197  }
2198  MagickPixelCompositeBlend(&gradient->stops[i].color,1.0-alpha,
2199  &gradient->stops[j].color,alpha,&composite);
2200  }
2201  break;
2202  }
2203  }
2204  MagickPixelCompositeOver(&composite,composite.opacity,&pixel,
2205  pixel.opacity,&pixel);
2206  SetPixelPacket(image,&pixel,q,indexes+x);
2207  q++;
2208  }
2209  if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
2210  status=MagickFalse;
2211  }
2212  image_view=DestroyCacheView(image_view);
2213  return(status);
2214 }
2215 
2216 /*
2217 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2218 % %
2219 % %
2220 % %
2221 % D r a w I m a g e %
2222 % %
2223 % %
2224 % %
2225 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2226 %
2227 % DrawImage() draws a graphic primitive on your image. The primitive
2228 % may be represented as a string or filename. Precede the filename with an
2229 % "at" sign (@) and the contents of the file are drawn on the image. You
2230 % can affect how text is drawn by setting one or more members of the draw
2231 % info structure.
2232 %
2233 % The format of the DrawImage method is:
2234 %
2235 % MagickBooleanType DrawImage(Image *image,const DrawInfo *draw_info)
2236 %
2237 % A description of each parameter follows:
2238 %
2239 % o image: the image.
2240 %
2241 % o draw_info: the draw info.
2242 %
2243 */
2244 
2245 static MagickBooleanType CheckPrimitiveExtent(MVGInfo *mvg_info,
2246  const double pad)
2247 {
2248  char
2249  **text = (char **) NULL;
2250 
2251  double
2252  extent;
2253 
2254  size_t
2255  quantum;
2256 
2257  ssize_t
2258  i;
2259 
2260  /*
2261  Check if there is enough storage for drawing primitives.
2262  */
2263  quantum=sizeof(**mvg_info->primitive_info);
2264  extent=(double) mvg_info->offset+pad+(PrimitiveExtentPad+1)*quantum;
2265  if (extent <= (double) *mvg_info->extent)
2266  return(MagickTrue);
2267  if ((extent >= (double) MAGICK_SSIZE_MAX) || (IsNaN(extent) != 0))
2268  return(MagickFalse);
2269  if (mvg_info->offset > 0)
2270  {
2271  text=(char **) AcquireQuantumMemory(mvg_info->offset,sizeof(*text));
2272  if (text == (char **) NULL)
2273  return(MagickFalse);
2274  for (i=0; i < mvg_info->offset; i++)
2275  text[i]=(*mvg_info->primitive_info)[i].text;
2276  }
2277  *mvg_info->primitive_info=(PrimitiveInfo *) ResizeQuantumMemory(
2278  *mvg_info->primitive_info,(size_t) (extent+1),quantum);
2279  if (*mvg_info->primitive_info != (PrimitiveInfo *) NULL)
2280  {
2281  if (text != (char **) NULL)
2282  text=(char **) RelinquishMagickMemory(text);
2283  *mvg_info->extent=(size_t) extent;
2284  for (i=mvg_info->offset+1; i <= (ssize_t) extent; i++)
2285  {
2286  (*mvg_info->primitive_info)[i].primitive=UndefinedPrimitive;
2287  (*mvg_info->primitive_info)[i].text=(char *) NULL;
2288  }
2289  return(MagickTrue);
2290  }
2291  /*
2292  Reallocation failed, allocate a primitive to facilitate unwinding.
2293  */
2294  if (text != (char **) NULL)
2295  {
2296  for (i=0; i < mvg_info->offset; i++)
2297  if (text[i] != (char *) NULL)
2298  text[i]=DestroyString(text[i]);
2299  text=(char **) RelinquishMagickMemory(text);
2300  }
2301  (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
2302  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
2303  *mvg_info->primitive_info=(PrimitiveInfo *) AcquireCriticalMemory((size_t)
2304  (PrimitiveExtentPad+1)*quantum);
2305  (void) memset(*mvg_info->primitive_info,0,(size_t) ((PrimitiveExtentPad+1)*
2306  quantum));
2307  *mvg_info->extent=1;
2308  mvg_info->offset=0;
2309  return(MagickFalse);
2310 }
2311 
2312 static inline double GetDrawValue(const char *magick_restrict string,
2313  char **magick_restrict sentinel)
2314 {
2315  char
2316  **magick_restrict q;
2317 
2318  double
2319  value;
2320 
2321  q=sentinel;
2322  value=InterpretLocaleValue(string,q);
2323  sentinel=q;
2324  return(value);
2325 }
2326 
2327 static int MVGMacroCompare(const void *target,const void *source)
2328 {
2329  const char
2330  *p,
2331  *q;
2332 
2333  p=(const char *) target;
2334  q=(const char *) source;
2335  return(strcmp(p,q));
2336 }
2337 
2338 static SplayTreeInfo *GetMVGMacros(const char *primitive)
2339 {
2340  char
2341  *macro,
2342  *token;
2343 
2344  const char
2345  *q;
2346 
2347  size_t
2348  extent;
2349 
2351  *macros;
2352 
2353  /*
2354  Scan graphic primitives for definitions and classes.
2355  */
2356  if (primitive == (const char *) NULL)
2357  return((SplayTreeInfo *) NULL);
2358  macros=NewSplayTree(MVGMacroCompare,RelinquishMagickMemory,
2359  RelinquishMagickMemory);
2360  macro=AcquireString(primitive);
2361  token=AcquireString(primitive);
2362  extent=strlen(token)+MagickPathExtent;
2363  for (q=primitive; *q != '\0'; )
2364  {
2365  if (GetNextToken(q,&q,extent,token) < 1)
2366  break;
2367  if (*token == '\0')
2368  break;
2369  if (LocaleCompare("push",token) == 0)
2370  {
2371  const char
2372  *end,
2373  *start;
2374 
2375  (void) GetNextToken(q,&q,extent,token);
2376  if (*q == '"')
2377  {
2378  char
2379  name[MagickPathExtent];
2380 
2381  const char
2382  *p;
2383 
2384  ssize_t
2385  n;
2386 
2387  /*
2388  Named macro (e.g. push graphic-context "wheel").
2389  */
2390  (void) GetNextToken(q,&q,extent,token);
2391  start=q;
2392  end=q;
2393  (void) CopyMagickString(name,token,MagickPathExtent);
2394  n=1;
2395  for (p=q; *p != '\0'; )
2396  {
2397  if (GetNextToken(p,&p,extent,token) < 1)
2398  break;
2399  if (*token == '\0')
2400  break;
2401  if (LocaleCompare(token,"pop") == 0)
2402  {
2403  end=p-strlen(token)-1;
2404  n--;
2405  }
2406  if (LocaleCompare(token,"push") == 0)
2407  n++;
2408  if ((n == 0) && (end > start))
2409  {
2410  /*
2411  Extract macro.
2412  */
2413  (void) GetNextToken(p,&p,extent,token);
2414  (void) CopyMagickString(macro,start,(size_t) (end-start));
2415  (void) AddValueToSplayTree(macros,ConstantString(name),
2416  ConstantString(macro));
2417  break;
2418  }
2419  }
2420  }
2421  }
2422  }
2423  token=DestroyString(token);
2424  macro=DestroyString(macro);
2425  return(macros);
2426 }
2427 
2428 static inline MagickBooleanType IsPoint(const char *point)
2429 {
2430  char
2431  *p;
2432 
2433  double
2434  value;
2435 
2436  value=GetDrawValue(point,&p);
2437  return((fabs(value) < MagickEpsilon) && (p == point) ? MagickFalse :
2438  MagickTrue);
2439 }
2440 
2441 static inline MagickBooleanType TracePoint(PrimitiveInfo *primitive_info,
2442  const PointInfo point)
2443 {
2444  primitive_info->coordinates=1;
2445  primitive_info->closed_subpath=MagickFalse;
2446  primitive_info->point=point;
2447  return(MagickTrue);
2448 }
2449 
2450 static MagickBooleanType RenderMVGContent(Image *image,
2451  const DrawInfo *draw_info,const size_t depth)
2452 {
2453 #define RenderImageTag "Render/Image"
2454 
2455  AffineMatrix
2456  affine,
2457  current;
2458 
2459  char
2460  key[2*MaxTextExtent],
2461  keyword[MaxTextExtent],
2462  geometry[MaxTextExtent],
2463  name[MaxTextExtent],
2464  *next_token,
2465  pattern[MaxTextExtent],
2466  *primitive,
2467  *token;
2468 
2469  const char
2470  *q;
2471 
2472  double
2473  angle,
2474  coordinates,
2475  cursor,
2476  factor,
2477  primitive_extent;
2478 
2479  DrawInfo
2480  *clone_info,
2481  **graphic_context;
2482 
2483  MagickBooleanType
2484  proceed;
2485 
2486  MagickStatusType
2487  status;
2488 
2489  MVGInfo
2490  mvg_info;
2491 
2492  PointInfo
2493  point;
2494 
2495  PixelPacket
2496  start_color;
2497 
2499  *primitive_info;
2500 
2501  PrimitiveType
2502  primitive_type;
2503 
2504  const char
2505  *p;
2506 
2507  ssize_t
2508  i,
2509  x;
2510 
2511  SegmentInfo
2512  bounds;
2513 
2514  size_t
2515  extent,
2516  number_points;
2517 
2519  *macros;
2520 
2521  ssize_t
2522  defsDepth,
2523  j,
2524  k,
2525  n,
2526  symbolDepth;
2527 
2528  TypeMetric
2529  metrics;
2530 
2531  assert(image != (Image *) NULL);
2532  assert(image->signature == MagickCoreSignature);
2533  assert(draw_info != (DrawInfo *) NULL);
2534  assert(draw_info->signature == MagickCoreSignature);
2535  if (IsEventLogging() != MagickFalse)
2536  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
2537  if (depth > MagickMaxRecursionDepth)
2538  ThrowBinaryImageException(DrawError,"VectorGraphicsNestedTooDeeply",
2539  image->filename);
2540  if ((draw_info->primitive == (char *) NULL) ||
2541  (*draw_info->primitive == '\0'))
2542  return(MagickFalse);
2543  if (draw_info->debug != MagickFalse)
2544  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"begin draw-image");
2545  if (SetImageStorageClass(image,DirectClass) == MagickFalse)
2546  return(MagickFalse);
2547  if (image->matte == MagickFalse)
2548  {
2549  status=SetImageAlphaChannel(image,OpaqueAlphaChannel);
2550  if (status == MagickFalse)
2551  return(MagickFalse);
2552  }
2553  primitive=(char *) NULL;
2554  if ((*draw_info->primitive == '@') && (strlen(draw_info->primitive) > 1) &&
2555  (*(draw_info->primitive+1) != '-') && (depth == 0))
2556  primitive=FileToString(draw_info->primitive,~0UL,&image->exception);
2557  else
2558  primitive=AcquireString(draw_info->primitive);
2559  if (primitive == (char *) NULL)
2560  return(MagickFalse);
2561  primitive_extent=(double) strlen(primitive);
2562  (void) SetImageArtifact(image,"mvg:vector-graphics",primitive);
2563  n=0;
2564  /*
2565  Allocate primitive info memory.
2566  */
2567  graphic_context=(DrawInfo **) AcquireMagickMemory(sizeof(*graphic_context));
2568  if (graphic_context == (DrawInfo **) NULL)
2569  {
2570  primitive=DestroyString(primitive);
2571  ThrowBinaryImageException(ResourceLimitError,"MemoryAllocationFailed",
2572  image->filename);
2573  }
2574  number_points=(size_t) PrimitiveExtentPad;
2575  primitive_info=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
2576  (number_points+1),sizeof(*primitive_info));
2577  if (primitive_info == (PrimitiveInfo *) NULL)
2578  {
2579  primitive=DestroyString(primitive);
2580  for ( ; n >= 0; n--)
2581  graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
2582  graphic_context=(DrawInfo **) RelinquishMagickMemory(graphic_context);
2583  ThrowBinaryImageException(ResourceLimitError,"MemoryAllocationFailed",
2584  image->filename);
2585  }
2586  (void) memset(primitive_info,0,(size_t) (number_points+1)*
2587  sizeof(*primitive_info));
2588  (void) memset(&mvg_info,0,sizeof(mvg_info));
2589  mvg_info.primitive_info=(&primitive_info);
2590  mvg_info.extent=(&number_points);
2591  mvg_info.exception=(&image->exception);
2592  graphic_context[n]=CloneDrawInfo((ImageInfo *) NULL,draw_info);
2593  graphic_context[n]->viewbox=image->page;
2594  if ((image->page.width == 0) || (image->page.height == 0))
2595  {
2596  graphic_context[n]->viewbox.width=image->columns;
2597  graphic_context[n]->viewbox.height=image->rows;
2598  }
2599  token=AcquireString(primitive);
2600  extent=strlen(token)+MaxTextExtent;
2601  cursor=0.0;
2602  defsDepth=0;
2603  symbolDepth=0;
2604  macros=GetMVGMacros(primitive);
2605  status=QueryColorDatabase("#000000",&start_color,&image->exception);
2606  for (q=primitive; *q != '\0'; )
2607  {
2608  /*
2609  Interpret graphic primitive.
2610  */
2611  if (GetNextToken(q,&q,MaxTextExtent,keyword) < 1)
2612  break;
2613  if (*keyword == '\0')
2614  break;
2615  if (*keyword == '#')
2616  {
2617  /*
2618  Comment.
2619  */
2620  while ((*q != '\n') && (*q != '\0'))
2621  q++;
2622  continue;
2623  }
2624  p=q-strlen(keyword)-1;
2625  primitive_type=UndefinedPrimitive;
2626  current=graphic_context[n]->affine;
2627  GetAffineMatrix(&affine);
2628  *token='\0';
2629  switch (*keyword)
2630  {
2631  case ';':
2632  break;
2633  case 'a':
2634  case 'A':
2635  {
2636  if (LocaleCompare("affine",keyword) == 0)
2637  {
2638  (void) GetNextToken(q,&q,extent,token);
2639  affine.sx=GetDrawValue(token,&next_token);
2640  if (token == next_token)
2641  ThrowPointExpectedException(image,token);
2642  (void) GetNextToken(q,&q,extent,token);
2643  if (*token == ',')
2644  (void) GetNextToken(q,&q,extent,token);
2645  affine.rx=GetDrawValue(token,&next_token);
2646  if (token == next_token)
2647  ThrowPointExpectedException(image,token);
2648  (void) GetNextToken(q,&q,extent,token);
2649  if (*token == ',')
2650  (void) GetNextToken(q,&q,extent,token);
2651  affine.ry=GetDrawValue(token,&next_token);
2652  if (token == next_token)
2653  ThrowPointExpectedException(image,token);
2654  (void) GetNextToken(q,&q,extent,token);
2655  if (*token == ',')
2656  (void) GetNextToken(q,&q,extent,token);
2657  affine.sy=GetDrawValue(token,&next_token);
2658  if (token == next_token)
2659  ThrowPointExpectedException(image,token);
2660  (void) GetNextToken(q,&q,extent,token);
2661  if (*token == ',')
2662  (void) GetNextToken(q,&q,extent,token);
2663  affine.tx=GetDrawValue(token,&next_token);
2664  if (token == next_token)
2665  ThrowPointExpectedException(image,token);
2666  (void) GetNextToken(q,&q,extent,token);
2667  if (*token == ',')
2668  (void) GetNextToken(q,&q,extent,token);
2669  affine.ty=GetDrawValue(token,&next_token);
2670  if (token == next_token)
2671  ThrowPointExpectedException(image,token);
2672  break;
2673  }
2674  if (LocaleCompare("arc",keyword) == 0)
2675  {
2676  primitive_type=ArcPrimitive;
2677  break;
2678  }
2679  status=MagickFalse;
2680  break;
2681  }
2682  case 'b':
2683  case 'B':
2684  {
2685  if (LocaleCompare("bezier",keyword) == 0)
2686  {
2687  primitive_type=BezierPrimitive;
2688  break;
2689  }
2690  if (LocaleCompare("border-color",keyword) == 0)
2691  {
2692  (void) GetNextToken(q,&q,extent,token);
2693  status&=QueryColorDatabase(token,&graphic_context[n]->border_color,
2694  &image->exception);
2695  break;
2696  }
2697  status=MagickFalse;
2698  break;
2699  }
2700  case 'c':
2701  case 'C':
2702  {
2703  if (LocaleCompare("class",keyword) == 0)
2704  {
2705  const char
2706  *mvg_class;
2707 
2708  (void) GetNextToken(q,&q,extent,token);
2709  if ((*token == '\0') || (*token == ';'))
2710  {
2711  status=MagickFalse;
2712  break;
2713  }
2714  if (LocaleCompare(token,graphic_context[n]->id) == 0)
2715  break;
2716  mvg_class=(const char *) GetValueFromSplayTree(macros,token);
2717  if ((graphic_context[n]->render != MagickFalse) &&
2718  (mvg_class != (const char *) NULL) && (p > primitive))
2719  {
2720  char
2721  *elements;
2722 
2723  ssize_t
2724  offset;
2725 
2726  /*
2727  Inject class elements in stream.
2728  */
2729  offset=(ssize_t) (p-primitive);
2730  elements=AcquireString(primitive);
2731  elements[offset]='\0';
2732  (void) ConcatenateString(&elements,mvg_class);
2733  (void) ConcatenateString(&elements,"\n");
2734  (void) ConcatenateString(&elements,q);
2735  primitive=DestroyString(primitive);
2736  primitive=elements;
2737  q=primitive+offset;
2738  }
2739  break;
2740  }
2741  if (LocaleCompare("clip-path",keyword) == 0)
2742  {
2743  const char
2744  *clip_path;
2745 
2746  /*
2747  Take a node from within the MVG document, and duplicate it here.
2748  */
2749  (void) GetNextToken(q,&q,extent,token);
2750  if (*token == '\0')
2751  {
2752  status=MagickFalse;
2753  break;
2754  }
2755  (void) CloneString(&graphic_context[n]->clip_mask,token);
2756  clip_path=(const char *) GetValueFromSplayTree(macros,token);
2757  if (clip_path != (const char *) NULL)
2758  {
2759  if (graphic_context[n]->clipping_mask != (Image *) NULL)
2760  graphic_context[n]->clipping_mask=
2761  DestroyImage(graphic_context[n]->clipping_mask);
2762  graphic_context[n]->clipping_mask=DrawClippingMask(image,
2763  graphic_context[n],token,clip_path,&image->exception);
2764  if (graphic_context[n]->compliance != SVGCompliance)
2765  {
2766  const char
2767  *clip_path;
2768 
2769  clip_path=(const char *) GetValueFromSplayTree(macros,
2770  graphic_context[n]->clip_mask);
2771  if (clip_path != (const char *) NULL)
2772  (void) SetImageArtifact(image,
2773  graphic_context[n]->clip_mask,clip_path);
2774  status&=DrawClipPath(image,graphic_context[n],
2775  graphic_context[n]->clip_mask);
2776  }
2777  }
2778  break;
2779  }
2780  if (LocaleCompare("clip-rule",keyword) == 0)
2781  {
2782  ssize_t
2783  fill_rule;
2784 
2785  (void) GetNextToken(q,&q,extent,token);
2786  fill_rule=ParseCommandOption(MagickFillRuleOptions,MagickFalse,
2787  token);
2788  if (fill_rule == -1)
2789  {
2790  status=MagickFalse;
2791  break;
2792  }
2793  graphic_context[n]->fill_rule=(FillRule) fill_rule;
2794  break;
2795  }
2796  if (LocaleCompare("clip-units",keyword) == 0)
2797  {
2798  ssize_t
2799  clip_units;
2800 
2801  (void) GetNextToken(q,&q,extent,token);
2802  clip_units=ParseCommandOption(MagickClipPathOptions,MagickFalse,
2803  token);
2804  if (clip_units == -1)
2805  {
2806  status=MagickFalse;
2807  break;
2808  }
2809  graphic_context[n]->clip_units=(ClipPathUnits) clip_units;
2810  if (clip_units == ObjectBoundingBox)
2811  {
2812  GetAffineMatrix(&current);
2813  affine.sx=draw_info->bounds.x2;
2814  affine.sy=draw_info->bounds.y2;
2815  affine.tx=draw_info->bounds.x1;
2816  affine.ty=draw_info->bounds.y1;
2817  break;
2818  }
2819  break;
2820  }
2821  if (LocaleCompare("circle",keyword) == 0)
2822  {
2823  primitive_type=CirclePrimitive;
2824  break;
2825  }
2826  if (LocaleCompare("color",keyword) == 0)
2827  {
2828  primitive_type=ColorPrimitive;
2829  break;
2830  }
2831  if (LocaleCompare("compliance",keyword) == 0)
2832  {
2833  /*
2834  MVG compliance associates a clipping mask with an image; SVG
2835  compliance associates a clipping mask with a graphics context.
2836  */
2837  (void) GetNextToken(q,&q,extent,token);
2838  graphic_context[n]->compliance=(ComplianceType) ParseCommandOption(
2839  MagickComplianceOptions,MagickFalse,token);
2840  break;
2841  }
2842  if (LocaleCompare("currentColor",keyword) == 0)
2843  {
2844  (void) GetNextToken(q,&q,extent,token);
2845  break;
2846  }
2847  status=MagickFalse;
2848  break;
2849  }
2850  case 'd':
2851  case 'D':
2852  {
2853  if (LocaleCompare("decorate",keyword) == 0)
2854  {
2855  ssize_t
2856  decorate;
2857 
2858  (void) GetNextToken(q,&q,extent,token);
2859  decorate=ParseCommandOption(MagickDecorateOptions,MagickFalse,
2860  token);
2861  if (decorate == -1)
2862  {
2863  status=MagickFalse;
2864  break;
2865  }
2866  graphic_context[n]->decorate=(DecorationType) decorate;
2867  break;
2868  }
2869  if (LocaleCompare("density",keyword) == 0)
2870  {
2871  (void) GetNextToken(q,&q,extent,token);
2872  (void) CloneString(&graphic_context[n]->density,token);
2873  break;
2874  }
2875  if (LocaleCompare("direction",keyword) == 0)
2876  {
2877  ssize_t
2878  direction;
2879 
2880  (void) GetNextToken(q,&q,extent,token);
2881  direction=ParseCommandOption(MagickDirectionOptions,MagickFalse,
2882  token);
2883  if (direction == -1)
2884  status=MagickFalse;
2885  else
2886  graphic_context[n]->direction=(DirectionType) direction;
2887  break;
2888  }
2889  status=MagickFalse;
2890  break;
2891  }
2892  case 'e':
2893  case 'E':
2894  {
2895  if (LocaleCompare("ellipse",keyword) == 0)
2896  {
2897  primitive_type=EllipsePrimitive;
2898  break;
2899  }
2900  if (LocaleCompare("encoding",keyword) == 0)
2901  {
2902  (void) GetNextToken(q,&q,extent,token);
2903  (void) CloneString(&graphic_context[n]->encoding,token);
2904  break;
2905  }
2906  status=MagickFalse;
2907  break;
2908  }
2909  case 'f':
2910  case 'F':
2911  {
2912  if (LocaleCompare("fill",keyword) == 0)
2913  {
2914  const char
2915  *mvg_class;
2916 
2917  (void) GetNextToken(q,&q,extent,token);
2918  if (graphic_context[n]->clip_path != MagickFalse)
2919  break;
2920  mvg_class=(const char *) GetValueFromSplayTree(macros,token);
2921  if (mvg_class != (const char *) NULL)
2922  {
2923  (void) DrawPatternPath(image,draw_info,mvg_class,
2924  &graphic_context[n]->fill_pattern);
2925  break;
2926  }
2927  (void) FormatLocaleString(pattern,MaxTextExtent,"%s",token);
2928  if (GetImageArtifact(image,pattern) != (const char *) NULL)
2929  {
2930  (void) DrawPatternPath(image,draw_info,token,
2931  &graphic_context[n]->fill_pattern);
2932  break;
2933  }
2934  status&=QueryColorDatabase(token,&graphic_context[n]->fill,
2935  &image->exception);
2936  if (graphic_context[n]->fill_opacity != OpaqueOpacity)
2937  graphic_context[n]->fill.opacity=ClampToQuantum(
2938  graphic_context[n]->fill_opacity);
2939  break;
2940  }
2941  if (LocaleCompare("fill-opacity",keyword) == 0)
2942  {
2943  double
2944  opacity;
2945 
2946  (void) GetNextToken(q,&q,extent,token);
2947  if (graphic_context[n]->clip_path != MagickFalse)
2948  break;
2949  factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
2950  opacity=MagickMin(MagickMax(factor*
2951  GetDrawValue(token,&next_token),0.0),1.0);
2952  if (token == next_token)
2953  ThrowPointExpectedException(image,token);
2954  if (graphic_context[n]->compliance == SVGCompliance)
2955  graphic_context[n]->fill_opacity*=(1.0-opacity);
2956  else
2957  graphic_context[n]->fill_opacity=(QuantumRange-
2958  graphic_context[n]->fill_opacity)*(1.0-opacity);
2959  if (graphic_context[n]->fill.opacity != TransparentOpacity)
2960  graphic_context[n]->fill.opacity=(Quantum)
2961  graphic_context[n]->fill_opacity;
2962  else
2963  graphic_context[n]->fill.opacity=ClampToQuantum(QuantumRange*
2964  (1.0-opacity));
2965  break;
2966  }
2967  if (LocaleCompare("fill-rule",keyword) == 0)
2968  {
2969  ssize_t
2970  fill_rule;
2971 
2972  (void) GetNextToken(q,&q,extent,token);
2973  fill_rule=ParseCommandOption(MagickFillRuleOptions,MagickFalse,
2974  token);
2975  if (fill_rule == -1)
2976  {
2977  status=MagickFalse;
2978  break;
2979  }
2980  graphic_context[n]->fill_rule=(FillRule) fill_rule;
2981  break;
2982  }
2983  if (LocaleCompare("font",keyword) == 0)
2984  {
2985  (void) GetNextToken(q,&q,extent,token);
2986  (void) CloneString(&graphic_context[n]->font,token);
2987  if (LocaleCompare("none",token) == 0)
2988  graphic_context[n]->font=(char *) RelinquishMagickMemory(
2989  graphic_context[n]->font);
2990  break;
2991  }
2992  if (LocaleCompare("font-family",keyword) == 0)
2993  {
2994  (void) GetNextToken(q,&q,extent,token);
2995  (void) CloneString(&graphic_context[n]->family,token);
2996  break;
2997  }
2998  if (LocaleCompare("font-size",keyword) == 0)
2999  {
3000  (void) GetNextToken(q,&q,extent,token);
3001  graphic_context[n]->pointsize=GetDrawValue(token,&next_token);
3002  if (token == next_token)
3003  ThrowPointExpectedException(image,token);
3004  break;
3005  }
3006  if (LocaleCompare("font-stretch",keyword) == 0)
3007  {
3008  ssize_t
3009  stretch;
3010 
3011  (void) GetNextToken(q,&q,extent,token);
3012  stretch=ParseCommandOption(MagickStretchOptions,MagickFalse,token);
3013  if (stretch == -1)
3014  {
3015  status=MagickFalse;
3016  break;
3017  }
3018  graphic_context[n]->stretch=(StretchType) stretch;
3019  break;
3020  }
3021  if (LocaleCompare("font-style",keyword) == 0)
3022  {
3023  ssize_t
3024  style;
3025 
3026  (void) GetNextToken(q,&q,extent,token);
3027  style=ParseCommandOption(MagickStyleOptions,MagickFalse,token);
3028  if (style == -1)
3029  {
3030  status=MagickFalse;
3031  break;
3032  }
3033  graphic_context[n]->style=(StyleType) style;
3034  break;
3035  }
3036  if (LocaleCompare("font-weight",keyword) == 0)
3037  {
3038  ssize_t
3039  weight;
3040 
3041  (void) GetNextToken(q,&q,extent,token);
3042  weight=ParseCommandOption(MagickWeightOptions,MagickFalse,token);
3043  if (weight == -1)
3044  weight=(ssize_t) StringToUnsignedLong(token);
3045  graphic_context[n]->weight=(size_t) weight;
3046  break;
3047  }
3048  status=MagickFalse;
3049  break;
3050  }
3051  case 'g':
3052  case 'G':
3053  {
3054  if (LocaleCompare("gradient-units",keyword) == 0)
3055  {
3056  (void) GetNextToken(q,&q,extent,token);
3057  break;
3058  }
3059  if (LocaleCompare("gravity",keyword) == 0)
3060  {
3061  ssize_t
3062  gravity;
3063 
3064  (void) GetNextToken(q,&q,extent,token);
3065  gravity=ParseCommandOption(MagickGravityOptions,MagickFalse,token);
3066  if (gravity == -1)
3067  {
3068  status=MagickFalse;
3069  break;
3070  }
3071  graphic_context[n]->gravity=(GravityType) gravity;
3072  break;
3073  }
3074  status=MagickFalse;
3075  break;
3076  }
3077  case 'i':
3078  case 'I':
3079  {
3080  if (LocaleCompare("image",keyword) == 0)
3081  {
3082  ssize_t
3083  compose;
3084 
3085  primitive_type=ImagePrimitive;
3086  (void) GetNextToken(q,&q,extent,token);
3087  compose=ParseCommandOption(MagickComposeOptions,MagickFalse,token);
3088  if (compose == -1)
3089  {
3090  status=MagickFalse;
3091  break;
3092  }
3093  graphic_context[n]->compose=(CompositeOperator) compose;
3094  break;
3095  }
3096  if (LocaleCompare("interline-spacing",keyword) == 0)
3097  {
3098  (void) GetNextToken(q,&q,extent,token);
3099  graphic_context[n]->interline_spacing=GetDrawValue(token,
3100  &next_token);
3101  if (token == next_token)
3102  ThrowPointExpectedException(image,token);
3103  break;
3104  }
3105  if (LocaleCompare("interword-spacing",keyword) == 0)
3106  {
3107  (void) GetNextToken(q,&q,extent,token);
3108  graphic_context[n]->interword_spacing=GetDrawValue(token,
3109  &next_token);
3110  if (token == next_token)
3111  ThrowPointExpectedException(image,token);
3112  break;
3113  }
3114  status=MagickFalse;
3115  break;
3116  }
3117  case 'k':
3118  case 'K':
3119  {
3120  if (LocaleCompare("kerning",keyword) == 0)
3121  {
3122  (void) GetNextToken(q,&q,extent,token);
3123  graphic_context[n]->kerning=GetDrawValue(token,&next_token);
3124  if (token == next_token)
3125  ThrowPointExpectedException(image,token);
3126  break;
3127  }
3128  status=MagickFalse;
3129  break;
3130  }
3131  case 'l':
3132  case 'L':
3133  {
3134  if (LocaleCompare("letter-spacing",keyword) == 0)
3135  {
3136  (void) GetNextToken(q,&q,extent,token);
3137  if (IsPoint(token) == MagickFalse)
3138  break;
3139  clone_info=CloneDrawInfo((ImageInfo *) NULL,graphic_context[n]);
3140  clone_info->text=AcquireString(" ");
3141  status&=GetTypeMetrics(image,clone_info,&metrics);
3142  graphic_context[n]->kerning=metrics.width*
3143  GetDrawValue(token,&next_token);
3144  clone_info=DestroyDrawInfo(clone_info);
3145  if (token == next_token)
3146  ThrowPointExpectedException(image,token);
3147  break;
3148  }
3149  if (LocaleCompare("line",keyword) == 0)
3150  {
3151  primitive_type=LinePrimitive;
3152  break;
3153  }
3154  status=MagickFalse;
3155  break;
3156  }
3157  case 'm':
3158  case 'M':
3159  {
3160  if (LocaleCompare("mask",keyword) == 0)
3161  {
3162  const char
3163  *mask_path;
3164 
3165  /*
3166  Take a node from within the MVG document, and duplicate it here.
3167  */
3168  (void) GetNextToken(q,&q,extent,token);
3169  mask_path=(const char *) GetValueFromSplayTree(macros,token);
3170  if (mask_path != (const char *) NULL)
3171  {
3172  if (graphic_context[n]->composite_mask != (Image *) NULL)
3173  graphic_context[n]->composite_mask=
3174  DestroyImage(graphic_context[n]->composite_mask);
3175  graphic_context[n]->composite_mask=DrawCompositeMask(image,
3176  graphic_context[n],token,mask_path,&image->exception);
3177  if (graphic_context[n]->compliance != SVGCompliance)
3178  status=SetImageMask(image,graphic_context[n]->composite_mask);
3179  }
3180  break;
3181  }
3182  if (LocaleCompare("matte",keyword) == 0)
3183  {
3184  primitive_type=MattePrimitive;
3185  break;
3186  }
3187  status=MagickFalse;
3188  break;
3189  }
3190  case 'o':
3191  case 'O':
3192  {
3193  if (LocaleCompare("offset",keyword) == 0)
3194  {
3195  (void) GetNextToken(q,&q,extent,token);
3196  break;
3197  }
3198  if (LocaleCompare("opacity",keyword) == 0)
3199  {
3200  double
3201  opacity;
3202 
3203  (void) GetNextToken(q,&q,extent,token);
3204  if (graphic_context[n]->clip_path != MagickFalse)
3205  break;
3206  factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
3207  opacity=MagickMin(MagickMax(factor*
3208  GetDrawValue(token,&next_token),0.0),1.0);
3209  if (token == next_token)
3210  ThrowPointExpectedException(image,token);
3211  if (graphic_context[n]->compliance == SVGCompliance)
3212  {
3213  graphic_context[n]->fill_opacity*=(1.0-opacity);
3214  graphic_context[n]->stroke_opacity*=(1.0-opacity);
3215  }
3216  else
3217  {
3218  graphic_context[n]->fill_opacity=(QuantumRange-
3219  graphic_context[n]->fill_opacity)*(1.0-opacity);
3220  graphic_context[n]->stroke_opacity=(QuantumRange-
3221  graphic_context[n]->stroke_opacity)*(1.0-opacity);
3222  }
3223  break;
3224  }
3225  status=MagickFalse;
3226  break;
3227  }
3228  case 'p':
3229  case 'P':
3230  {
3231  if (LocaleCompare("path",keyword) == 0)
3232  {
3233  primitive_type=PathPrimitive;
3234  break;
3235  }
3236  if (LocaleCompare("point",keyword) == 0)
3237  {
3238  primitive_type=PointPrimitive;
3239  break;
3240  }
3241  if (LocaleCompare("polyline",keyword) == 0)
3242  {
3243  primitive_type=PolylinePrimitive;
3244  break;
3245  }
3246  if (LocaleCompare("polygon",keyword) == 0)
3247  {
3248  primitive_type=PolygonPrimitive;
3249  break;
3250  }
3251  if (LocaleCompare("pop",keyword) == 0)
3252  {
3253  if (GetNextToken(q,&q,extent,token) < 1)
3254  break;
3255  if (LocaleCompare("class",token) == 0)
3256  break;
3257  if (LocaleCompare("clip-path",token) == 0)
3258  break;
3259  if (LocaleCompare("defs",token) == 0)
3260  {
3261  defsDepth--;
3262  graphic_context[n]->render=defsDepth > 0 ? MagickFalse :
3263  MagickTrue;
3264  break;
3265  }
3266  if (LocaleCompare("gradient",token) == 0)
3267  break;
3268  if (LocaleCompare("graphic-context",token) == 0)
3269  {
3270  if (n <= 0)
3271  {
3272  (void) ThrowMagickException(&image->exception,
3273  GetMagickModule(),DrawError,
3274  "UnbalancedGraphicContextPushPop","`%s'",token);
3275  status=MagickFalse;
3276  n=0;
3277  break;
3278  }
3279  if ((graphic_context[n]->clip_mask != (char *) NULL) &&
3280  (graphic_context[n]->compliance != SVGCompliance))
3281  if (LocaleCompare(graphic_context[n]->clip_mask,
3282  graphic_context[n-1]->clip_mask) != 0)
3283  status=SetImageClipMask(image,(Image *) NULL);
3284  graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
3285  n--;
3286  break;
3287  }
3288  if (LocaleCompare("mask",token) == 0)
3289  break;
3290  if (LocaleCompare("pattern",token) == 0)
3291  break;
3292  if (LocaleCompare("symbol",token) == 0)
3293  {
3294  symbolDepth--;
3295  graphic_context[n]->render=symbolDepth > 0 ? MagickFalse :
3296  MagickTrue;
3297  break;
3298  }
3299  status=MagickFalse;
3300  break;
3301  }
3302  if (LocaleCompare("push",keyword) == 0)
3303  {
3304  if (GetNextToken(q,&q,extent,token) < 1)
3305  break;
3306  if (LocaleCompare("class",token) == 0)
3307  {
3308  /*
3309  Class context.
3310  */
3311  for (p=q; *q != '\0'; )
3312  {
3313  if (GetNextToken(q,&q,extent,token) < 1)
3314  break;
3315  if (LocaleCompare(token,"pop") != 0)
3316  continue;
3317  (void) GetNextToken(q,(const char **) NULL,extent,token);
3318  if (LocaleCompare(token,"class") != 0)
3319  continue;
3320  break;
3321  }
3322  (void) GetNextToken(q,&q,extent,token);
3323  break;
3324  }
3325  if (LocaleCompare("clip-path",token) == 0)
3326  {
3327  (void) GetNextToken(q,&q,extent,token);
3328  for (p=q; *q != '\0'; )
3329  {
3330  if (GetNextToken(q,&q,extent,token) < 1)
3331  break;
3332  if (LocaleCompare(token,"pop") != 0)
3333  continue;
3334  (void) GetNextToken(q,(const char **) NULL,extent,token);
3335  if (LocaleCompare(token,"clip-path") != 0)
3336  continue;
3337  break;
3338  }
3339  if ((q == (char *) NULL) || (p == (char *) NULL) || ((q-4) < p))
3340  {
3341  status=MagickFalse;
3342  break;
3343  }
3344  (void) GetNextToken(q,&q,extent,token);
3345  break;
3346  }
3347  if (LocaleCompare("defs",token) == 0)
3348  {
3349  defsDepth++;
3350  graphic_context[n]->render=defsDepth > 0 ? MagickFalse :
3351  MagickTrue;
3352  break;
3353  }
3354  if (LocaleCompare("gradient",token) == 0)
3355  {
3356  char
3357  key[2*MaxTextExtent],
3358  name[MaxTextExtent],
3359  type[MaxTextExtent];
3360 
3361  SegmentInfo
3362  segment;
3363 
3364  (void) GetNextToken(q,&q,extent,token);
3365  (void) CopyMagickString(name,token,MaxTextExtent);
3366  (void) GetNextToken(q,&q,extent,token);
3367  (void) CopyMagickString(type,token,MaxTextExtent);
3368  (void) GetNextToken(q,&q,extent,token);
3369  segment.x1=GetDrawValue(token,&next_token);
3370  if (token == next_token)
3371  ThrowPointExpectedException(image,token);
3372  (void) GetNextToken(q,&q,extent,token);
3373  if (*token == ',')
3374  (void) GetNextToken(q,&q,extent,token);
3375  segment.y1=GetDrawValue(token,&next_token);
3376  if (token == next_token)
3377  ThrowPointExpectedException(image,token);
3378  (void) GetNextToken(q,&q,extent,token);
3379  if (*token == ',')
3380  (void) GetNextToken(q,&q,extent,token);
3381  segment.x2=GetDrawValue(token,&next_token);
3382  if (token == next_token)
3383  ThrowPointExpectedException(image,token);
3384  (void) GetNextToken(q,&q,extent,token);
3385  if (*token == ',')
3386  (void) GetNextToken(q,&q,extent,token);
3387  segment.y2=GetDrawValue(token,&next_token);
3388  if (token == next_token)
3389  ThrowPointExpectedException(image,token);
3390  if (LocaleCompare(type,"radial") == 0)
3391  {
3392  (void) GetNextToken(q,&q,extent,token);
3393  if (*token == ',')
3394  (void) GetNextToken(q,&q,extent,token);
3395  }
3396  for (p=q; *q != '\0'; )
3397  {
3398  if (GetNextToken(q,&q,extent,token) < 1)
3399  break;
3400  if (LocaleCompare(token,"pop") != 0)
3401  continue;
3402  (void) GetNextToken(q,(const char **) NULL,extent,token);
3403  if (LocaleCompare(token,"gradient") != 0)
3404  continue;
3405  break;
3406  }
3407  if ((q == (char *) NULL) || (*q == '\0') ||
3408  (p == (char *) NULL) || ((q-4) < p))
3409  {
3410  status=MagickFalse;
3411  break;
3412  }
3413  (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
3414  bounds.x1=graphic_context[n]->affine.sx*segment.x1+
3415  graphic_context[n]->affine.ry*segment.y1+
3416  graphic_context[n]->affine.tx;
3417  bounds.y1=graphic_context[n]->affine.rx*segment.x1+
3418  graphic_context[n]->affine.sy*segment.y1+
3419  graphic_context[n]->affine.ty;
3420  bounds.x2=graphic_context[n]->affine.sx*segment.x2+
3421  graphic_context[n]->affine.ry*segment.y2+
3422  graphic_context[n]->affine.tx;
3423  bounds.y2=graphic_context[n]->affine.rx*segment.x2+
3424  graphic_context[n]->affine.sy*segment.y2+
3425  graphic_context[n]->affine.ty;
3426  (void) FormatLocaleString(key,MaxTextExtent,"%s",name);
3427  (void) SetImageArtifact(image,key,token);
3428  (void) FormatLocaleString(key,MaxTextExtent,"%s-type",name);
3429  (void) SetImageArtifact(image,key,type);
3430  (void) FormatLocaleString(key,MaxTextExtent,"%s-geometry",name);
3431  (void) FormatLocaleString(geometry,MaxTextExtent,
3432  "%gx%g%+.15g%+.15g",
3433  MagickMax(fabs(bounds.x2-bounds.x1+1.0),1.0),
3434  MagickMax(fabs(bounds.y2-bounds.y1+1.0),1.0),
3435  bounds.x1,bounds.y1);
3436  (void) SetImageArtifact(image,key,geometry);
3437  (void) GetNextToken(q,&q,extent,token);
3438  break;
3439  }
3440  if (LocaleCompare("graphic-context",token) == 0)
3441  {
3442  n++;
3443  graphic_context=(DrawInfo **) ResizeQuantumMemory(
3444  graphic_context,(size_t) (n+1),sizeof(*graphic_context));
3445  if (graphic_context == (DrawInfo **) NULL)
3446  {
3447  (void) ThrowMagickException(&image->exception,
3448  GetMagickModule(),ResourceLimitError,
3449  "MemoryAllocationFailed","`%s'",image->filename);
3450  break;
3451  }
3452  graphic_context[n]=CloneDrawInfo((ImageInfo *) NULL,
3453  graphic_context[n-1]);
3454  if (*q == '"')
3455  {
3456  (void) GetNextToken(q,&q,extent,token);
3457  (void) CloneString(&graphic_context[n]->id,token);
3458  }
3459  break;
3460  }
3461  if (LocaleCompare("mask",token) == 0)
3462  {
3463  (void) GetNextToken(q,&q,extent,token);
3464  break;
3465  }
3466  if (LocaleCompare("pattern",token) == 0)
3467  {
3469  bounds;
3470 
3471  (void) GetNextToken(q,&q,extent,token);
3472  (void) CopyMagickString(name,token,MaxTextExtent);
3473  (void) GetNextToken(q,&q,extent,token);
3474  bounds.x=CastDoubleToLong(ceil(GetDrawValue(token,
3475  &next_token)-0.5));
3476  if (token == next_token)
3477  ThrowPointExpectedException(image,token);
3478  (void) GetNextToken(q,&q,extent,token);
3479  if (*token == ',')
3480  (void) GetNextToken(q,&q,extent,token);
3481  bounds.y=CastDoubleToLong(ceil(GetDrawValue(token,
3482  &next_token)-0.5));
3483  if (token == next_token)
3484  ThrowPointExpectedException(image,token);
3485  (void) GetNextToken(q,&q,extent,token);
3486  if (*token == ',')
3487  (void) GetNextToken(q,&q,extent,token);
3488  bounds.width=CastDoubleToUnsigned(GetDrawValue(token,&next_token)+
3489  0.5);
3490  if (token == next_token)
3491  ThrowPointExpectedException(image,token);
3492  (void) GetNextToken(q,&q,extent,token);
3493  if (*token == ',')
3494  (void) GetNextToken(q,&q,extent,token);
3495  bounds.height=CastDoubleToUnsigned(GetDrawValue(token,&next_token)+
3496  0.5);
3497  if (token == next_token)
3498  ThrowPointExpectedException(image,token);
3499  for (p=q; *q != '\0'; )
3500  {
3501  if (GetNextToken(q,&q,extent,token) < 1)
3502  break;
3503  if (LocaleCompare(token,"pop") != 0)
3504  continue;
3505  (void) GetNextToken(q,(const char **) NULL,extent,token);
3506  if (LocaleCompare(token,"pattern") != 0)
3507  continue;
3508  break;
3509  }
3510  if ((q == (char *) NULL) || (p == (char *) NULL) || ((q-4) < p))
3511  {
3512  status=MagickFalse;
3513  break;
3514  }
3515  (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
3516  (void) FormatLocaleString(key,MaxTextExtent,"%s",name);
3517  (void) SetImageArtifact(image,key,token);
3518  (void) FormatLocaleString(key,MaxTextExtent,"%s-geometry",name);
3519  (void) FormatLocaleString(geometry,MaxTextExtent,
3520  "%.20gx%.20g%+.20g%+.20g",(double) bounds.width,(double)
3521  bounds.height,(double) bounds.x,(double) bounds.y);
3522  (void) SetImageArtifact(image,key,geometry);
3523  (void) GetNextToken(q,&q,extent,token);
3524  break;
3525  }
3526  if (LocaleCompare("symbol",token) == 0)
3527  {
3528  symbolDepth++;
3529  graphic_context[n]->render=symbolDepth > 0 ? MagickFalse :
3530  MagickTrue;
3531  break;
3532  }
3533  status=MagickFalse;
3534  break;
3535  }
3536  status=MagickFalse;
3537  break;
3538  }
3539  case 'r':
3540  case 'R':
3541  {
3542  if (LocaleCompare("rectangle",keyword) == 0)
3543  {
3544  primitive_type=RectanglePrimitive;
3545  break;
3546  }
3547  if (LocaleCompare("rotate",keyword) == 0)
3548  {
3549  (void) GetNextToken(q,&q,extent,token);
3550  angle=GetDrawValue(token,&next_token);
3551  if (token == next_token)
3552  ThrowPointExpectedException(image,token);
3553  affine.sx=cos(DegreesToRadians(fmod((double) angle,360.0)));
3554  affine.rx=sin(DegreesToRadians(fmod((double) angle,360.0)));
3555  affine.ry=(-sin(DegreesToRadians(fmod((double) angle,360.0))));
3556  affine.sy=cos(DegreesToRadians(fmod((double) angle,360.0)));
3557  break;
3558  }
3559  if (LocaleCompare("roundRectangle",keyword) == 0)
3560  {
3561  primitive_type=RoundRectanglePrimitive;
3562  break;
3563  }
3564  status=MagickFalse;
3565  break;
3566  }
3567  case 's':
3568  case 'S':
3569  {
3570  if (LocaleCompare("scale",keyword) == 0)
3571  {
3572  (void) GetNextToken(q,&q,extent,token);
3573  affine.sx=GetDrawValue(token,&next_token);
3574  if (token == next_token)
3575  ThrowPointExpectedException(image,token);
3576  (void) GetNextToken(q,&q,extent,token);
3577  if (*token == ',')
3578  (void) GetNextToken(q,&q,extent,token);
3579  affine.sy=GetDrawValue(token,&next_token);
3580  if (token == next_token)
3581  ThrowPointExpectedException(image,token);
3582  break;
3583  }
3584  if (LocaleCompare("skewX",keyword) == 0)
3585  {
3586  (void) GetNextToken(q,&q,extent,token);
3587  angle=GetDrawValue(token,&next_token);
3588  if (token == next_token)
3589  ThrowPointExpectedException(image,token);
3590  affine.ry=sin(DegreesToRadians(angle));
3591  break;
3592  }
3593  if (LocaleCompare("skewY",keyword) == 0)
3594  {
3595  (void) GetNextToken(q,&q,extent,token);
3596  angle=GetDrawValue(token,&next_token);
3597  if (token == next_token)
3598  ThrowPointExpectedException(image,token);
3599  affine.rx=(-tan(DegreesToRadians(angle)/2.0));
3600  break;
3601  }
3602  if (LocaleCompare("stop-color",keyword) == 0)
3603  {
3604  GradientType
3605  type;
3606 
3607  PixelPacket
3608  stop_color;
3609 
3610  (void) GetNextToken(q,&q,extent,token);
3611  status&=QueryColorDatabase(token,&stop_color,&image->exception);
3612  type=LinearGradient;
3613  if (draw_info->gradient.type == RadialGradient)
3614  type=RadialGradient;
3615  (void) GradientImage(image,type,PadSpread,&start_color,&stop_color);
3616  start_color=stop_color;
3617  (void) GetNextToken(q,&q,extent,token);
3618  break;
3619  }
3620  if (LocaleCompare("stroke",keyword) == 0)
3621  {
3622  const char
3623  *mvg_class;
3624 
3625  (void) GetNextToken(q,&q,extent,token);
3626  if (graphic_context[n]->clip_path != MagickFalse)
3627  break;
3628  mvg_class=(const char *) GetValueFromSplayTree(macros,token);
3629  if (mvg_class != (const char *) NULL)
3630  {
3631  (void) DrawPatternPath(image,draw_info,mvg_class,
3632  &graphic_context[n]->stroke_pattern);
3633  break;
3634  }
3635  (void) FormatLocaleString(pattern,MaxTextExtent,"%s",token);
3636  if (GetImageArtifact(image,pattern) != (const char *) NULL)
3637  {
3638  (void) DrawPatternPath(image,draw_info,token,
3639  &graphic_context[n]->stroke_pattern);
3640  break;
3641  }
3642  status&=QueryColorDatabase(token,&graphic_context[n]->stroke,
3643  &image->exception);
3644  if (graphic_context[n]->stroke_opacity != OpaqueOpacity)
3645  graphic_context[n]->stroke.opacity=ClampToQuantum(
3646  graphic_context[n]->stroke_opacity);
3647  break;
3648  }
3649  if (LocaleCompare("stroke-antialias",keyword) == 0)
3650  {
3651  (void) GetNextToken(q,&q,extent,token);
3652  graphic_context[n]->stroke_antialias=StringToLong(token) != 0 ?
3653  MagickTrue : MagickFalse;
3654  break;
3655  }
3656  if (LocaleCompare("stroke-dasharray",keyword) == 0)
3657  {
3658  if (graphic_context[n]->dash_pattern != (double *) NULL)
3659  graphic_context[n]->dash_pattern=(double *)
3660  RelinquishMagickMemory(graphic_context[n]->dash_pattern);
3661  if (IsPoint(q) != MagickFalse)
3662  {
3663  const char
3664  *p;
3665 
3666  p=q;
3667  (void) GetNextToken(p,&p,extent,token);
3668  if (*token == ',')
3669  (void) GetNextToken(p,&p,extent,token);
3670  for (x=0; IsPoint(token) != MagickFalse; x++)
3671  {
3672  (void) GetNextToken(p,&p,extent,token);
3673  if (*token == ',')
3674  (void) GetNextToken(p,&p,extent,token);
3675  }
3676  graphic_context[n]->dash_pattern=(double *)
3677  AcquireQuantumMemory((size_t) (2*x+2),
3678  sizeof(*graphic_context[n]->dash_pattern));
3679  if (graphic_context[n]->dash_pattern == (double *) NULL)
3680  {
3681  (void) ThrowMagickException(&image->exception,
3682  GetMagickModule(),ResourceLimitError,
3683  "MemoryAllocationFailed","`%s'",image->filename);
3684  status=MagickFalse;
3685  break;
3686  }
3687  (void) memset(graphic_context[n]->dash_pattern,0,(size_t)
3688  (2*x+2)*sizeof(*graphic_context[n]->dash_pattern));
3689  for (j=0; j < x; j++)
3690  {
3691  (void) GetNextToken(q,&q,extent,token);
3692  if (*token == ',')
3693  (void) GetNextToken(q,&q,extent,token);
3694  graphic_context[n]->dash_pattern[j]=GetDrawValue(token,
3695  &next_token);
3696  if (token == next_token)
3697  ThrowPointExpectedException(image,token);
3698  if (graphic_context[n]->dash_pattern[j] <= 0.0)
3699  status=MagickFalse;
3700  }
3701  if ((x & 0x01) != 0)
3702  for ( ; j < (2*x); j++)
3703  graphic_context[n]->dash_pattern[j]=
3704  graphic_context[n]->dash_pattern[j-x];
3705  graphic_context[n]->dash_pattern[j]=0.0;
3706  break;
3707  }
3708  (void) GetNextToken(q,&q,extent,token);
3709  break;
3710  }
3711  if (LocaleCompare("stroke-dashoffset",keyword) == 0)
3712  {
3713  (void) GetNextToken(q,&q,extent,token);
3714  graphic_context[n]->dash_offset=GetDrawValue(token,&next_token);
3715  if (token == next_token)
3716  ThrowPointExpectedException(image,token);
3717  break;
3718  }
3719  if (LocaleCompare("stroke-linecap",keyword) == 0)
3720  {
3721  ssize_t
3722  linecap;
3723 
3724  (void) GetNextToken(q,&q,extent,token);
3725  linecap=ParseCommandOption(MagickLineCapOptions,MagickFalse,token);
3726  if (linecap == -1)
3727  {
3728  status=MagickFalse;
3729  break;
3730  }
3731  graphic_context[n]->linecap=(LineCap) linecap;
3732  break;
3733  }
3734  if (LocaleCompare("stroke-linejoin",keyword) == 0)
3735  {
3736  ssize_t
3737  linejoin;
3738 
3739  (void) GetNextToken(q,&q,extent,token);
3740  linejoin=ParseCommandOption(MagickLineJoinOptions,MagickFalse,
3741  token);
3742  if (linejoin == -1)
3743  {
3744  status=MagickFalse;
3745  break;
3746  }
3747  graphic_context[n]->linejoin=(LineJoin) linejoin;
3748  break;
3749  }
3750  if (LocaleCompare("stroke-miterlimit",keyword) == 0)
3751  {
3752  (void) GetNextToken(q,&q,extent,token);
3753  graphic_context[n]->miterlimit=StringToUnsignedLong(token);
3754  break;
3755  }
3756  if (LocaleCompare("stroke-opacity",keyword) == 0)
3757  {
3758  double
3759  opacity;
3760 
3761  (void) GetNextToken(q,&q,extent,token);
3762  if (graphic_context[n]->clip_path != MagickFalse)
3763  break;
3764  factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
3765  opacity=MagickMin(MagickMax(factor*
3766  GetDrawValue(token,&next_token),0.0),1.0);
3767  if (token == next_token)
3768  ThrowPointExpectedException(image,token);
3769  if (graphic_context[n]->compliance == SVGCompliance)
3770  graphic_context[n]->stroke_opacity*=(1.0-opacity);
3771  else
3772  graphic_context[n]->stroke_opacity=(QuantumRange-
3773  graphic_context[n]->stroke_opacity)*(1.0-opacity);
3774  if (graphic_context[n]->stroke.opacity != TransparentOpacity)
3775  graphic_context[n]->stroke.opacity=(Quantum)
3776  graphic_context[n]->stroke_opacity;
3777  else
3778  graphic_context[n]->stroke.opacity=ClampToQuantum(QuantumRange*
3779  opacity);
3780  break;
3781  }
3782  if (LocaleCompare("stroke-width",keyword) == 0)
3783  {
3784  (void) GetNextToken(q,&q,extent,token);
3785  if (graphic_context[n]->clip_path != MagickFalse)
3786  break;
3787  graphic_context[n]->stroke_width=GetDrawValue(token,&next_token);
3788  if ((token == next_token) ||
3789  (graphic_context[n]->stroke_width < 0.0))
3790  ThrowPointExpectedException(image,token);
3791  break;
3792  }
3793  status=MagickFalse;
3794  break;
3795  }
3796  case 't':
3797  case 'T':
3798  {
3799  if (LocaleCompare("text",keyword) == 0)
3800  {
3801  primitive_type=TextPrimitive;
3802  cursor=0.0;
3803  break;
3804  }
3805  if (LocaleCompare("text-align",keyword) == 0)
3806  {
3807  ssize_t
3808  align;
3809 
3810  (void) GetNextToken(q,&q,extent,token);
3811  align=ParseCommandOption(MagickAlignOptions,MagickFalse,token);
3812  if (align == -1)
3813  {
3814  status=MagickFalse;
3815  break;
3816  }
3817  graphic_context[n]->align=(AlignType) align;
3818  break;
3819  }
3820  if (LocaleCompare("text-anchor",keyword) == 0)
3821  {
3822  ssize_t
3823  align;
3824 
3825  (void) GetNextToken(q,&q,extent,token);
3826  align=ParseCommandOption(MagickAlignOptions,MagickFalse,token);
3827  if (align == -1)
3828  {
3829  status=MagickFalse;
3830  break;
3831  }
3832  graphic_context[n]->align=(AlignType) align;
3833  break;
3834  }
3835  if (LocaleCompare("text-antialias",keyword) == 0)
3836  {
3837  (void) GetNextToken(q,&q,extent,token);
3838  graphic_context[n]->text_antialias=StringToLong(token) != 0 ?
3839  MagickTrue : MagickFalse;
3840  break;
3841  }
3842  if (LocaleCompare("text-undercolor",keyword) == 0)
3843  {
3844  (void) GetNextToken(q,&q,extent,token);
3845  status&=QueryColorDatabase(token,&graphic_context[n]->undercolor,
3846  &image->exception);
3847  break;
3848  }
3849  if (LocaleCompare("translate",keyword) == 0)
3850  {
3851  (void) GetNextToken(q,&q,extent,token);
3852  affine.tx=GetDrawValue(token,&next_token);
3853  if (token == next_token)
3854  ThrowPointExpectedException(image,token);
3855  (void) GetNextToken(q,&q,extent,token);
3856  if (*token == ',')
3857  (void) GetNextToken(q,&q,extent,token);
3858  affine.ty=GetDrawValue(token,&next_token);
3859  if (token == next_token)
3860  ThrowPointExpectedException(image,token);
3861  break;
3862  }
3863  status=MagickFalse;
3864  break;
3865  }
3866  case 'u':
3867  case 'U':
3868  {
3869  if (LocaleCompare("use",keyword) == 0)
3870  {
3871  const char
3872  *use;
3873 
3874  /*
3875  Get a macro from the MVG document, and "use" it here.
3876  */
3877  (void) GetNextToken(q,&q,extent,token);
3878  use=(const char *) GetValueFromSplayTree(macros,token);
3879  if (use != (const char *) NULL)
3880  {
3881  clone_info=CloneDrawInfo((ImageInfo *) NULL,graphic_context[n]);
3882  (void) CloneString(&clone_info->primitive,use);
3883  status=RenderMVGContent(image,clone_info,depth+1);
3884  clone_info=DestroyDrawInfo(clone_info);
3885  }
3886  break;
3887  }
3888  status=MagickFalse;
3889  break;
3890  }
3891  case 'v':
3892  case 'V':
3893  {
3894  if (LocaleCompare("viewbox",keyword) == 0)
3895  {
3896  (void) GetNextToken(q,&q,extent,token);
3897  graphic_context[n]->viewbox.x=CastDoubleToLong(ceil(
3898  GetDrawValue(token,&next_token)-0.5));
3899  if (token == next_token)
3900  ThrowPointExpectedException(image,token);
3901  (void) GetNextToken(q,&q,extent,token);
3902  if (*token == ',')
3903  (void) GetNextToken(q,&q,extent,token);
3904  graphic_context[n]->viewbox.y=CastDoubleToLong(ceil(
3905  GetDrawValue(token,&next_token)-0.5));
3906  if (token == next_token)
3907  ThrowPointExpectedException(image,token);
3908  (void) GetNextToken(q,&q,extent,token);
3909  if (*token == ',')
3910  (void) GetNextToken(q,&q,extent,token);
3911  graphic_context[n]->viewbox.width=CastDoubleToUnsigned(GetDrawValue(
3912  token,&next_token)+0.5);
3913  if (token == next_token)
3914  ThrowPointExpectedException(image,token);
3915  (void) GetNextToken(q,&q,extent,token);
3916  if (*token == ',')
3917  (void) GetNextToken(q,&q,extent,token);
3918  graphic_context[n]->viewbox.height=CastDoubleToUnsigned(GetDrawValue(
3919  token,&next_token)+0.5);
3920  if (token == next_token)
3921  ThrowPointExpectedException(image,token);
3922  break;
3923  }
3924  status=MagickFalse;
3925  break;
3926  }
3927  case 'w':
3928  case 'W':
3929  {
3930  if (LocaleCompare("word-spacing",keyword) == 0)
3931  {
3932  (void) GetNextToken(q,&q,extent,token);
3933  graphic_context[n]->interword_spacing=GetDrawValue(token,
3934  &next_token);
3935  if (token == next_token)
3936  ThrowPointExpectedException(image,token);
3937  break;
3938  }
3939  status=MagickFalse;
3940  break;
3941  }
3942  default:
3943  {
3944  status=MagickFalse;
3945  break;
3946  }
3947  }
3948  if (status == MagickFalse)
3949  break;
3950  if ((fabs(affine.sx-1.0) >= MagickEpsilon) ||
3951  (fabs(affine.rx) >= MagickEpsilon) || (fabs(affine.ry) >= MagickEpsilon) ||
3952  (fabs(affine.sy-1.0) >= MagickEpsilon) ||
3953  (fabs(affine.tx) >= MagickEpsilon) || (fabs(affine.ty) >= MagickEpsilon))
3954  {
3955  graphic_context[n]->affine.sx=current.sx*affine.sx+current.ry*affine.rx;
3956  graphic_context[n]->affine.rx=current.rx*affine.sx+current.sy*affine.rx;
3957  graphic_context[n]->affine.ry=current.sx*affine.ry+current.ry*affine.sy;
3958  graphic_context[n]->affine.sy=current.rx*affine.ry+current.sy*affine.sy;
3959  graphic_context[n]->affine.tx=current.sx*affine.tx+current.ry*affine.ty+
3960  current.tx;
3961  graphic_context[n]->affine.ty=current.rx*affine.tx+current.sy*affine.ty+
3962  current.ty;
3963  }
3964  if (primitive_type == UndefinedPrimitive)
3965  {
3966  if ((draw_info->debug != MagickFalse) && (q > p))
3967  (void) LogMagickEvent(DrawEvent,GetMagickModule()," %.*s",(int)
3968  (q-p-1),p);
3969  continue;
3970  }
3971  /*
3972  Parse the primitive attributes.
3973  */
3974  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
3975  if (primitive_info[i].text != (char *) NULL)
3976  primitive_info[i].text=DestroyString(primitive_info[i].text);
3977  i=0;
3978  mvg_info.offset=i;
3979  j=0;
3980  primitive_info[0].point.x=0.0;
3981  primitive_info[0].point.y=0.0;
3982  primitive_info[0].coordinates=0;
3983  primitive_info[0].method=FloodfillMethod;
3984  primitive_info[0].closed_subpath=MagickFalse;
3985  for (x=0; *q != '\0'; x++)
3986  {
3987  /*
3988  Define points.
3989  */
3990  if (IsPoint(q) == MagickFalse)
3991  break;
3992  (void) GetNextToken(q,&q,extent,token);
3993  point.x=GetDrawValue(token,&next_token);
3994  if (token == next_token)
3995  ThrowPointExpectedException(image,token);
3996  (void) GetNextToken(q,&q,extent,token);
3997  if (*token == ',')
3998  (void) GetNextToken(q,&q,extent,token);
3999  point.y=GetDrawValue(token,&next_token);
4000  if (token == next_token)
4001  ThrowPointExpectedException(image,token);
4002  (void) GetNextToken(q,(const char **) NULL,extent,token);
4003  if (*token == ',')
4004  (void) GetNextToken(q,&q,extent,token);
4005  primitive_info[i].primitive=primitive_type;
4006  primitive_info[i].point=point;
4007  primitive_info[i].coordinates=0;
4008  primitive_info[i].method=FloodfillMethod;
4009  primitive_info[i].closed_subpath=MagickFalse;
4010  i++;
4011  mvg_info.offset=i;
4012  if (i < (ssize_t) number_points)
4013  continue;
4014  status&=CheckPrimitiveExtent(&mvg_info,(double) number_points);
4015  primitive_info=(*mvg_info.primitive_info);
4016  }
4017  if (status == MagickFalse)
4018  break;
4019  if (primitive_info[j].text != (char *) NULL)
4020  primitive_info[j].text=DestroyString(primitive_info[j].text);
4021  primitive_info[j].primitive=primitive_type;
4022  primitive_info[j].coordinates=(size_t) x;
4023  primitive_info[j].method=FloodfillMethod;
4024  primitive_info[j].closed_subpath=MagickFalse;
4025  /*
4026  Circumscribe primitive within a circle.
4027  */
4028  bounds.x1=primitive_info[j].point.x;
4029  bounds.y1=primitive_info[j].point.y;
4030  bounds.x2=primitive_info[j].point.x;
4031  bounds.y2=primitive_info[j].point.y;
4032  for (k=1; k < (ssize_t) primitive_info[j].coordinates; k++)
4033  {
4034  point=primitive_info[j+k].point;
4035  if (point.x < bounds.x1)
4036  bounds.x1=point.x;
4037  if (point.y < bounds.y1)
4038  bounds.y1=point.y;
4039  if (point.x > bounds.x2)
4040  bounds.x2=point.x;
4041  if (point.y > bounds.y2)
4042  bounds.y2=point.y;
4043  }
4044  /*
4045  Speculate how many points our primitive might consume.
4046  */
4047  coordinates=(double) primitive_info[j].coordinates;
4048  switch (primitive_type)
4049  {
4050  case RectanglePrimitive:
4051  {
4052  coordinates*=5.0;
4053  break;
4054  }
4055  case RoundRectanglePrimitive:
4056  {
4057  double
4058  alpha,
4059  beta,
4060  radius;
4061 
4062  alpha=bounds.x2-bounds.x1;
4063  beta=bounds.y2-bounds.y1;
4064  radius=hypot(alpha,beta);
4065  coordinates*=5.0;
4066  coordinates+=2.0*((size_t) ceil((double) MagickPI*radius))+6.0*
4067  BezierQuantum+360.0;
4068  break;
4069  }
4070  case BezierPrimitive:
4071  {
4072  coordinates=(BezierQuantum*(double) primitive_info[j].coordinates);
4073  break;
4074  }
4075  case PathPrimitive:
4076  {
4077  char
4078  *s,
4079  *t;
4080 
4081  (void) GetNextToken(q,&q,extent,token);
4082  coordinates=1.0;
4083  t=token;
4084  for (s=token; *s != '\0'; s=t)
4085  {
4086  double
4087  value;
4088 
4089  value=GetDrawValue(s,&t);
4090  (void) value;
4091  if (s == t)
4092  {
4093  t++;
4094  continue;
4095  }
4096  coordinates++;
4097  }
4098  for (s=token; *s != '\0'; s++)
4099  if (strspn(s,"AaCcQqSsTt") != 0)
4100  coordinates+=(20.0*BezierQuantum)+360.0;
4101  break;
4102  }
4103  default:
4104  break;
4105  }
4106  if (status == MagickFalse)
4107  break;
4108  if (((size_t) (i+coordinates)) >= number_points)
4109  {
4110  /*
4111  Resize based on speculative points required by primitive.
4112  */
4113  number_points+=coordinates+1;
4114  if (number_points < (size_t) coordinates)
4115  {
4116  (void) ThrowMagickException(&image->exception,GetMagickModule(),
4117  ResourceLimitError,"MemoryAllocationFailed","`%s'",
4118  image->filename);
4119  break;
4120  }
4121  mvg_info.offset=i;
4122  status&=CheckPrimitiveExtent(&mvg_info,(double) number_points);
4123  primitive_info=(*mvg_info.primitive_info);
4124  }
4125  status&=CheckPrimitiveExtent(&mvg_info,PrimitiveExtentPad);
4126  primitive_info=(*mvg_info.primitive_info);
4127  if (status == MagickFalse)
4128  break;
4129  mvg_info.offset=j;
4130  switch (primitive_type)
4131  {
4132  case PointPrimitive:
4133  default:
4134  {
4135  if (primitive_info[j].coordinates != 1)
4136  {
4137  status=MagickFalse;
4138  break;
4139  }
4140  status&=TracePoint(primitive_info+j,primitive_info[j].point);
4141  primitive_info=(*mvg_info.primitive_info);
4142  i=(ssize_t) (j+primitive_info[j].coordinates);
4143  break;
4144  }
4145  case LinePrimitive:
4146  {
4147  if (primitive_info[j].coordinates != 2)
4148  {
4149  status=MagickFalse;
4150  break;
4151  }
4152  status&=TraceLine(primitive_info+j,primitive_info[j].point,
4153  primitive_info[j+1].point);
4154  primitive_info=(*mvg_info.primitive_info);
4155  i=(ssize_t) (j+primitive_info[j].coordinates);
4156  break;
4157  }
4158  case RectanglePrimitive:
4159  {
4160  if (primitive_info[j].coordinates != 2)
4161  {
4162  status=MagickFalse;
4163  break;
4164  }
4165  status&=TraceRectangle(primitive_info+j,primitive_info[j].point,
4166  primitive_info[j+1].point);
4167  primitive_info=(*mvg_info.primitive_info);
4168  i=(ssize_t) (j+primitive_info[j].coordinates);
4169  break;
4170  }
4171  case RoundRectanglePrimitive:
4172  {
4173  if (primitive_info[j].coordinates != 3)
4174  {
4175  status=MagickFalse;
4176  break;
4177  }
4178  if ((primitive_info[j+2].point.x < 0.0) ||
4179  (primitive_info[j+2].point.y < 0.0))
4180  {
4181  status=MagickFalse;
4182  break;
4183  }
4184  if ((primitive_info[j+1].point.x-primitive_info[j].point.x) < 0.0)
4185  {
4186  status=MagickFalse;
4187  break;
4188  }
4189  if ((primitive_info[j+1].point.y-primitive_info[j].point.y) < 0.0)
4190  {
4191  status=MagickFalse;
4192  break;
4193  }
4194  status&=TraceRoundRectangle(&mvg_info,primitive_info[j].point,
4195  primitive_info[j+1].point,primitive_info[j+2].point);
4196  primitive_info=(*mvg_info.primitive_info);
4197  i=(ssize_t) (j+primitive_info[j].coordinates);
4198  break;
4199  }
4200  case ArcPrimitive:
4201  {
4202  if (primitive_info[j].coordinates != 3)
4203  {
4204  status=MagickFalse;
4205  break;
4206  }
4207  status&=TraceArc(&mvg_info,primitive_info[j].point,
4208  primitive_info[j+1].point,primitive_info[j+2].point);
4209  primitive_info=(*mvg_info.primitive_info);
4210  i=(ssize_t) (j+primitive_info[j].coordinates);
4211  break;
4212  }
4213  case EllipsePrimitive:
4214  {
4215  if (primitive_info[j].coordinates != 3)
4216  {
4217  status=MagickFalse;
4218  break;
4219  }
4220  if ((primitive_info[j+1].point.x < 0.0) ||
4221  (primitive_info[j+1].point.y < 0.0))
4222  {
4223  status=MagickFalse;
4224  break;
4225  }
4226  status&=TraceEllipse(&mvg_info,primitive_info[j].point,
4227  primitive_info[j+1].point,primitive_info[j+2].point);
4228  primitive_info=(*mvg_info.primitive_info);
4229  i=(ssize_t) (j+primitive_info[j].coordinates);
4230  break;
4231  }
4232  case CirclePrimitive:
4233  {
4234  if (primitive_info[j].coordinates != 2)
4235  {
4236  status=MagickFalse;
4237  break;
4238  }
4239  status&=TraceCircle(&mvg_info,primitive_info[j].point,
4240  primitive_info[j+1].point);
4241  primitive_info=(*mvg_info.primitive_info);
4242  i=(ssize_t) (j+primitive_info[j].coordinates);
4243  break;
4244  }
4245  case PolylinePrimitive:
4246  {
4247  if (primitive_info[j].coordinates < 1)
4248  {
4249  status=MagickFalse;
4250  break;
4251  }
4252  break;
4253  }
4254  case PolygonPrimitive:
4255  {
4256  if (primitive_info[j].coordinates < 3)
4257  {
4258  status=MagickFalse;
4259  break;
4260  }
4261  primitive_info[i]=primitive_info[j];
4262  primitive_info[i].coordinates=0;
4263  primitive_info[j].coordinates++;
4264  primitive_info[j].closed_subpath=MagickTrue;
4265  i++;
4266  break;
4267  }
4268  case BezierPrimitive:
4269  {
4270  if (primitive_info[j].coordinates < 3)
4271  {
4272  status=MagickFalse;
4273  break;
4274  }
4275  status&=TraceBezier(&mvg_info,primitive_info[j].coordinates);
4276  primitive_info=(*mvg_info.primitive_info);
4277  i=(ssize_t) (j+primitive_info[j].coordinates);
4278  break;
4279  }
4280  case PathPrimitive:
4281  {
4282  coordinates=(double) TracePath(image,&mvg_info,token);
4283  primitive_info=(*mvg_info.primitive_info);
4284  if (coordinates < 0.0)
4285  {
4286  status=MagickFalse;
4287  break;
4288  }
4289  i=(ssize_t) (j+coordinates);
4290  break;
4291  }
4292  case ColorPrimitive:
4293  case MattePrimitive:
4294  {
4295  ssize_t
4296  method;
4297 
4298  if (primitive_info[j].coordinates != 1)
4299  {
4300  status=MagickFalse;
4301  break;
4302  }
4303  (void) GetNextToken(q,&q,extent,token);
4304  method=ParseCommandOption(MagickMethodOptions,MagickFalse,token);
4305  if (method == -1)
4306  {
4307  status=MagickFalse;
4308  break;
4309  }
4310  primitive_info[j].method=(PaintMethod) method;
4311  break;
4312  }
4313  case TextPrimitive:
4314  {
4315  char
4316  geometry[MagickPathExtent];
4317 
4318  if (primitive_info[j].coordinates != 1)
4319  {
4320  status=MagickFalse;
4321  break;
4322  }
4323  if (*token != ',')
4324  (void) GetNextToken(q,&q,extent,token);
4325  (void) CloneString(&primitive_info[j].text,token);
4326  /*
4327  Compute text cursor offset.
4328  */
4329  clone_info=CloneDrawInfo((ImageInfo *) NULL,graphic_context[n]);
4330  if ((fabs(mvg_info.point.x-primitive_info->point.x) < MagickEpsilon) &&
4331  (fabs(mvg_info.point.y-primitive_info->point.y) < MagickEpsilon))
4332  {
4333  mvg_info.point=primitive_info->point;
4334  primitive_info->point.x+=cursor;
4335  }
4336  else
4337  {
4338  mvg_info.point=primitive_info->point;
4339  cursor=0.0;
4340  }
4341  (void) FormatLocaleString(geometry,MagickPathExtent,"%+f%+f",
4342  primitive_info->point.x,primitive_info->point.y);
4343  clone_info->render=MagickFalse;
4344  clone_info->text=AcquireString(token);
4345  status&=GetTypeMetrics(image,clone_info,&metrics);
4346  clone_info=DestroyDrawInfo(clone_info);
4347  cursor+=metrics.width;
4348  if (graphic_context[n]->compliance != SVGCompliance)
4349  cursor=0.0;
4350  break;
4351  }
4352  case ImagePrimitive:
4353  {
4354  if (primitive_info[j].coordinates != 2)
4355  {
4356  status=MagickFalse;
4357  break;
4358  }
4359  (void) GetNextToken(q,&q,extent,token);
4360  (void) CloneString(&primitive_info[j].text,token);
4361  break;
4362  }
4363  }
4364  mvg_info.offset=i;
4365  if (status == 0)
4366  break;
4367  primitive_info[i].primitive=UndefinedPrimitive;
4368  if ((draw_info->debug != MagickFalse) && (q > p))
4369  (void) LogMagickEvent(DrawEvent,GetMagickModule()," %.*s",(int) (q-p),p);
4370  /*
4371  Sanity check.
4372  */
4373  status&=CheckPrimitiveExtent(&mvg_info,ExpandAffine(
4374  &graphic_context[n]->affine));
4375  primitive_info=(*mvg_info.primitive_info);
4376  if (status == 0)
4377  break;
4378  status&=CheckPrimitiveExtent(&mvg_info,(double)
4379  graphic_context[n]->stroke_width);
4380  primitive_info=(*mvg_info.primitive_info);
4381  if (status == 0)
4382  break;
4383  if (i == 0)
4384  continue;
4385  /*
4386  Transform points.
4387  */
4388  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
4389  {
4390  point=primitive_info[i].point;
4391  primitive_info[i].point.x=graphic_context[n]->affine.sx*point.x+
4392  graphic_context[n]->affine.ry*point.y+graphic_context[n]->affine.tx;
4393  primitive_info[i].point.y=graphic_context[n]->affine.rx*point.x+
4394  graphic_context[n]->affine.sy*point.y+graphic_context[n]->affine.ty;
4395  point=primitive_info[i].point;
4396  if (point.x < graphic_context[n]->bounds.x1)
4397  graphic_context[n]->bounds.x1=point.x;
4398  if (point.y < graphic_context[n]->bounds.y1)
4399  graphic_context[n]->bounds.y1=point.y;
4400  if (point.x > graphic_context[n]->bounds.x2)
4401  graphic_context[n]->bounds.x2=point.x;
4402  if (point.y > graphic_context[n]->bounds.y2)
4403  graphic_context[n]->bounds.y2=point.y;
4404  if (primitive_info[i].primitive == ImagePrimitive)
4405  break;
4406  if (i >= (ssize_t) number_points)
4407  ThrowFatalException(ResourceLimitFatalError,"MemoryAllocationFailed");
4408  }
4409  if (graphic_context[n]->render != MagickFalse)
4410  {
4411  if ((n != 0) && (graphic_context[n]->compliance != SVGCompliance) &&
4412  (graphic_context[n]->clip_mask != (char *) NULL) &&
4413  (LocaleCompare(graphic_context[n]->clip_mask,
4414  graphic_context[n-1]->clip_mask) != 0))
4415  {
4416  const char
4417  *clip_path;
4418 
4419  clip_path=(const char *) GetValueFromSplayTree(macros,
4420  graphic_context[n]->clip_mask);
4421  if (clip_path != (const char *) NULL)
4422  (void) SetImageArtifact(image,graphic_context[n]->clip_mask,
4423  clip_path);
4424  status&=DrawClipPath(image,graphic_context[n],
4425  graphic_context[n]->clip_mask);
4426  }
4427  status&=DrawPrimitive(image,graphic_context[n],primitive_info);
4428  }
4429  proceed=SetImageProgress(image,RenderImageTag,q-primitive,(MagickSizeType)
4430  primitive_extent);
4431  if (proceed == MagickFalse)
4432  break;
4433  if (status == 0)
4434  break;
4435  }
4436  if (draw_info->debug != MagickFalse)
4437  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end draw-image");
4438  /*
4439  Relinquish resources.
4440  */
4441  macros=DestroySplayTree(macros);
4442  token=DestroyString(token);
4443  if (primitive_info != (PrimitiveInfo *) NULL)
4444  {
4445  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
4446  if (primitive_info[i].text != (char *) NULL)
4447  primitive_info[i].text=DestroyString(primitive_info[i].text);
4448  primitive_info=(PrimitiveInfo *) RelinquishMagickMemory(primitive_info);
4449  }
4450  primitive=DestroyString(primitive);
4451  for ( ; n >= 0; n--)
4452  graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
4453  graphic_context=(DrawInfo **) RelinquishMagickMemory(graphic_context);
4454  if (status == MagickFalse)
4455  ThrowBinaryImageException(DrawError,
4456  "NonconformingDrawingPrimitiveDefinition",keyword);
4457  return(status != 0 ? MagickTrue : MagickFalse);
4458 }
4459 
4460 MagickExport MagickBooleanType DrawImage(Image *image,const DrawInfo *draw_info)
4461 {
4462  return(RenderMVGContent(image,draw_info,0));
4463 }
4464 
4465 /*
4466 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4467 % %
4468 % %
4469 % %
4470 % D r a w P a t t e r n P a t h %
4471 % %
4472 % %
4473 % %
4474 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4475 %
4476 % DrawPatternPath() draws a pattern.
4477 %
4478 % The format of the DrawPatternPath method is:
4479 %
4480 % MagickBooleanType DrawPatternPath(Image *image,const DrawInfo *draw_info,
4481 % const char *name,Image **pattern)
4482 %
4483 % A description of each parameter follows:
4484 %
4485 % o image: the image.
4486 %
4487 % o draw_info: the draw info.
4488 %
4489 % o name: the pattern name.
4490 %
4491 % o image: the image.
4492 %
4493 */
4494 MagickExport MagickBooleanType DrawPatternPath(Image *image,
4495  const DrawInfo *draw_info,const char *name,Image **pattern)
4496 {
4497  char
4498  property[MaxTextExtent];
4499 
4500  const char
4501  *geometry,
4502  *path,
4503  *type;
4504 
4505  DrawInfo
4506  *clone_info;
4507 
4508  ImageInfo
4509  *image_info;
4510 
4511  MagickBooleanType
4512  status;
4513 
4514  assert(image != (Image *) NULL);
4515  assert(image->signature == MagickCoreSignature);
4516  assert(draw_info != (const DrawInfo *) NULL);
4517  if (IsEventLogging() != MagickFalse)
4518  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
4519  assert(name != (const char *) NULL);
4520  (void) FormatLocaleString(property,MaxTextExtent,"%s",name);
4521  path=GetImageArtifact(image,property);
4522  if (path == (const char *) NULL)
4523  return(MagickFalse);
4524  (void) FormatLocaleString(property,MaxTextExtent,"%s-geometry",name);
4525  geometry=GetImageArtifact(image,property);
4526  if (geometry == (const char *) NULL)
4527  return(MagickFalse);
4528  if ((*pattern) != (Image *) NULL)
4529  *pattern=DestroyImage(*pattern);
4530  image_info=AcquireImageInfo();
4531  image_info->size=AcquireString(geometry);
4532  *pattern=AcquireImage(image_info);
4533  image_info=DestroyImageInfo(image_info);
4534  (void) QueryColorDatabase("#00000000",&(*pattern)->background_color,
4535  &image->exception);
4536  (void) SetImageBackgroundColor(*pattern);
4537  if (draw_info->debug != MagickFalse)
4538  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
4539  "begin pattern-path %s %s",name,geometry);
4540  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
4541  if (clone_info->fill_pattern != (Image *) NULL)
4542  clone_info->fill_pattern=DestroyImage(clone_info->fill_pattern);
4543  if (clone_info->stroke_pattern != (Image *) NULL)
4544  clone_info->stroke_pattern=DestroyImage(clone_info->stroke_pattern);
4545  (void) FormatLocaleString(property,MaxTextExtent,"%s-type",name);
4546  type=GetImageArtifact(image,property);
4547  if (type != (const char *) NULL)
4548  clone_info->gradient.type=(GradientType) ParseCommandOption(
4549  MagickGradientOptions,MagickFalse,type);
4550  (void) CloneString(&clone_info->primitive,path);
4551  status=RenderMVGContent(*pattern,clone_info,0);
4552  clone_info=DestroyDrawInfo(clone_info);
4553  if (draw_info->debug != MagickFalse)
4554  (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end pattern-path");
4555  return(status);
4556 }
4557 
4558 /*
4559 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4560 % %
4561 % %
4562 % %
4563 + D r a w P o l y g o n P r i m i t i v e %
4564 % %
4565 % %
4566 % %
4567 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4568 %
4569 % DrawPolygonPrimitive() draws a polygon on the image.
4570 %
4571 % The format of the DrawPolygonPrimitive method is:
4572 %
4573 % MagickBooleanType DrawPolygonPrimitive(Image *image,
4574 % const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
4575 %
4576 % A description of each parameter follows:
4577 %
4578 % o image: the image.
4579 %
4580 % o draw_info: the draw info.
4581 %
4582 % o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
4583 %
4584 */
4585 
4586 static PolygonInfo **DestroyPolygonTLS(PolygonInfo **polygon_info)
4587 {
4588  ssize_t
4589  i;
4590 
4591  assert(polygon_info != (PolygonInfo **) NULL);
4592  for (i=0; i < (ssize_t) GetMagickResourceLimit(ThreadResource); i++)
4593  if (polygon_info[i] != (PolygonInfo *) NULL)
4594  polygon_info[i]=DestroyPolygonInfo(polygon_info[i]);
4595  polygon_info=(PolygonInfo **) RelinquishMagickMemory(polygon_info);
4596  return(polygon_info);
4597 }
4598 
4599 static PolygonInfo **AcquirePolygonTLS(const DrawInfo *draw_info,
4600  const PrimitiveInfo *primitive_info,ExceptionInfo *exception)
4601 {
4602  PathInfo
4603  *magick_restrict path_info;
4604 
4605  PolygonInfo
4606  **polygon_info;
4607 
4608  size_t
4609  number_threads;
4610 
4611  number_threads=(size_t) GetMagickResourceLimit(ThreadResource);
4612  polygon_info=(PolygonInfo **) AcquireQuantumMemory(number_threads,
4613  sizeof(*polygon_info));
4614  if (polygon_info == (PolygonInfo **) NULL)
4615  {
4616  (void) ThrowMagickException(exception,GetMagickModule(),
4617  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
4618  return((PolygonInfo **) NULL);
4619  }
4620  (void) memset(polygon_info,0,number_threads*sizeof(*polygon_info));
4621  path_info=ConvertPrimitiveToPath(draw_info,primitive_info,exception);
4622  if (path_info == (PathInfo *) NULL)
4623  return(DestroyPolygonTLS(polygon_info));
4624  polygon_info[0]=ConvertPathToPolygon(path_info,exception);
4625  if (polygon_info[0] == (PolygonInfo *) NULL)
4626  {
4627  (void) ThrowMagickException(exception,GetMagickModule(),
4628  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
4629  return(DestroyPolygonTLS(polygon_info));
4630  }
4631  path_info=(PathInfo *) RelinquishMagickMemory(path_info);
4632  return(polygon_info);
4633 }
4634 
4635 static MagickBooleanType AcquirePolygonEdgesTLS(PolygonInfo **polygon_info,
4636  const size_t number_threads,ExceptionInfo *exception)
4637 {
4638  ssize_t
4639  i;
4640 
4641  for (i=1; i < (ssize_t) number_threads; i++)
4642  {
4643  EdgeInfo
4644  *edge_info;
4645 
4646  ssize_t
4647  j;
4648 
4649  polygon_info[i]=(PolygonInfo *) AcquireMagickMemory(
4650  sizeof(*polygon_info[i]));
4651  if (polygon_info[i] == (PolygonInfo *) NULL)
4652  {
4653  (void) ThrowMagickException(exception,GetMagickModule(),
4654  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
4655  return(MagickFalse);
4656  }
4657  polygon_info[i]->number_edges=0;
4658  edge_info=polygon_info[0]->edges;
4659  polygon_info[i]->edges=(EdgeInfo *) AcquireQuantumMemory(
4660  polygon_info[0]->number_edges,sizeof(*edge_info));
4661  if (polygon_info[i]->edges == (EdgeInfo *) NULL)
4662  {
4663  (void) ThrowMagickException(exception,GetMagickModule(),
4664  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
4665  return(MagickFalse);
4666  }
4667  (void) memcpy(polygon_info[i]->edges,edge_info,
4668  polygon_info[0]->number_edges*sizeof(*edge_info));
4669  for (j=0; j < (ssize_t) polygon_info[i]->number_edges; j++)
4670  polygon_info[i]->edges[j].points=(PointInfo *) NULL;
4671  polygon_info[i]->number_edges=polygon_info[0]->number_edges;
4672  for (j=0; j < (ssize_t) polygon_info[i]->number_edges; j++)
4673  {
4674  edge_info=polygon_info[0]->edges+j;
4675  polygon_info[i]->edges[j].points=(PointInfo *) AcquireQuantumMemory(
4676  edge_info->number_points,sizeof(*edge_info));
4677  if (polygon_info[i]->edges[j].points == (PointInfo *) NULL)
4678  {
4679  (void) ThrowMagickException(exception,GetMagickModule(),
4680  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
4681  return(MagickFalse);
4682  }
4683  (void) memcpy(polygon_info[i]->edges[j].points,edge_info->points,
4684  edge_info->number_points*sizeof(*edge_info->points));
4685  }
4686  }
4687  return(MagickTrue);
4688 }
4689 
4690 static size_t DestroyEdge(PolygonInfo *polygon_info,const ssize_t edge)
4691 {
4692  assert(edge < (ssize_t) polygon_info->number_edges);
4693  polygon_info->edges[edge].points=(PointInfo *) RelinquishMagickMemory(
4694  polygon_info->edges[edge].points);
4695  polygon_info->number_edges--;
4696  if (edge < (ssize_t) polygon_info->number_edges)
4697  (void) memmove(polygon_info->edges+edge,polygon_info->edges+edge+1,
4698  (size_t) (polygon_info->number_edges-edge)*sizeof(*polygon_info->edges));
4699  return(polygon_info->number_edges);
4700 }
4701 
4702 static double GetOpacityPixel(PolygonInfo *polygon_info,const double mid,
4703  const MagickBooleanType fill,const FillRule fill_rule,const ssize_t x,
4704  const ssize_t y,double *stroke_opacity)
4705 {
4706  double
4707  alpha,
4708  beta,
4709  distance,
4710  subpath_opacity;
4711 
4712  PointInfo
4713  delta;
4714 
4715  EdgeInfo
4716  *p;
4717 
4718  const PointInfo
4719  *q;
4720 
4721  ssize_t
4722  i;
4723 
4724  ssize_t
4725  j,
4726  winding_number;
4727 
4728  /*
4729  Compute fill & stroke opacity for this (x,y) point.
4730  */
4731  *stroke_opacity=0.0;
4732  subpath_opacity=0.0;
4733  p=polygon_info->edges;
4734  for (j=0; j < (ssize_t) polygon_info->number_edges; j++, p++)
4735  {
4736  if ((double) y <= (p->bounds.y1-mid-0.5))
4737  break;
4738  if ((double) y > (p->bounds.y2+mid+0.5))
4739  {
4740  p--;
4741  (void) DestroyEdge(polygon_info,j--);
4742  continue;
4743  }
4744  if (((double) x <= (p->bounds.x1-mid-0.5)) ||
4745  ((double) x > (p->bounds.x2+mid+0.5)))
4746  continue;
4747  i=(ssize_t) MagickMax((double) p->highwater,1.0);
4748  for ( ; i < (ssize_t) p->number_points; i++)
4749  {
4750  if ((double) y <= (p->points[i-1].y-mid-0.5))
4751  break;
4752  if ((double) y > (p->points[i].y+mid+0.5))
4753  continue;
4754  if (p->scanline != (double) y)
4755  {
4756  p->scanline=(double) y;
4757  p->highwater=(size_t) i;
4758  }
4759  /*
4760  Compute distance between a point and an edge.
4761  */
4762  q=p->points+i-1;
4763  delta.x=(q+1)->x-q->x;
4764  delta.y=(q+1)->y-q->y;
4765  beta=delta.x*(x-q->x)+delta.y*(y-q->y);
4766  if (beta <= 0.0)
4767  {
4768  delta.x=(double) x-q->x;
4769  delta.y=(double) y-q->y;
4770  distance=delta.x*delta.x+delta.y*delta.y;
4771  }
4772  else
4773  {
4774  alpha=delta.x*delta.x+delta.y*delta.y;
4775  if (beta >= alpha)
4776  {
4777  delta.x=(double) x-(q+1)->x;
4778  delta.y=(double) y-(q+1)->y;
4779  distance=delta.x*delta.x+delta.y*delta.y;
4780  }
4781  else
4782  {
4783  alpha=PerceptibleReciprocal(alpha);
4784  beta=delta.x*(y-q->y)-delta.y*(x-q->x);
4785  distance=alpha*beta*beta;
4786  }
4787  }
4788  /*
4789  Compute stroke & subpath opacity.
4790  */
4791  beta=0.0;
4792  if (p->ghostline == MagickFalse)
4793  {
4794  alpha=mid+0.5;
4795  if ((*stroke_opacity < 1.0) &&
4796  (distance <= ((alpha+0.25)*(alpha+0.25))))
4797  {
4798  alpha=mid-0.5;
4799  if (distance <= ((alpha+0.25)*(alpha+0.25)))
4800  *stroke_opacity=1.0;
4801  else
4802  {
4803  beta=1.0;
4804  if (fabs(distance-1.0) >= MagickEpsilon)
4805  beta=sqrt((double) distance);
4806  alpha=beta-mid-0.5;
4807  if (*stroke_opacity < ((alpha-0.25)*(alpha-0.25)))
4808  *stroke_opacity=(alpha-0.25)*(alpha-0.25);
4809  }
4810  }
4811  }
4812  if ((fill == MagickFalse) || (distance > 1.0) || (subpath_opacity >= 1.0))
4813  continue;
4814  if (distance <= 0.0)
4815  {
4816  subpath_opacity=1.0;
4817  continue;
4818  }
4819  if (distance > 1.0)
4820  continue;
4821  if (fabs(beta) < MagickEpsilon)
4822  {
4823  beta=1.0;
4824  if (fabs(distance-1.0) >= MagickEpsilon)
4825  beta=sqrt(distance);
4826  }
4827  alpha=beta-1.0;
4828  if (subpath_opacity < (alpha*alpha))
4829  subpath_opacity=alpha*alpha;
4830  }
4831  }
4832  /*
4833  Compute fill opacity.
4834  */
4835  if (fill == MagickFalse)
4836  return(0.0);
4837  if (subpath_opacity >= 1.0)
4838  return(1.0);
4839  /*
4840  Determine winding number.
4841  */
4842  winding_number=0;
4843  p=polygon_info->edges;
4844  for (j=0; j < (ssize_t) polygon_info->number_edges; j++, p++)
4845  {
4846  if ((double) y <= p->bounds.y1)
4847  break;
4848  if (((double) y > p->bounds.y2) || ((double) x <= p->bounds.x1))
4849  continue;
4850  if ((double) x > p->bounds.x2)
4851  {
4852  winding_number+=p->direction != 0 ? 1 : -1;
4853  continue;
4854  }
4855  i=(ssize_t) MagickMax((double) p->highwater,1.0);
4856  for ( ; i < (ssize_t) (p->number_points-1); i++)
4857  if ((double) y <= p->points[i].y)
4858  break;
4859  q=p->points+i-1;
4860  if ((((q+1)->x-q->x)*(y-q->y)) <= (((q+1)->y-q->y)*(x-q->x)))
4861  winding_number+=p->direction != 0 ? 1 : -1;
4862  }
4863  if (fill_rule != NonZeroRule)
4864  {
4865  if ((MagickAbsoluteValue(winding_number) & 0x01) != 0)
4866  return(1.0);
4867  }
4868  else
4869  if (MagickAbsoluteValue(winding_number) != 0)
4870  return(1.0);
4871  return(subpath_opacity);
4872 }
4873 
4874 static MagickBooleanType DrawPolygonPrimitive(Image *image,
4875  const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
4876 {
4877  typedef struct _ExtentInfo
4878  {
4879  ssize_t
4880  x1,
4881  y1,
4882  x2,
4883  y2;
4884  } ExtentInfo;
4885 
4886  CacheView
4887  *image_view;
4888 
4889  const char
4890  *artifact;
4891 
4892  double
4893  mid;
4894 
4896  *exception;
4897 
4898  ExtentInfo
4899  poly_extent;
4900 
4901  MagickBooleanType
4902  fill,
4903  status;
4904 
4905  PolygonInfo
4906  **magick_restrict polygon_info;
4907 
4908  EdgeInfo
4909  *p;
4910 
4911  SegmentInfo
4912  bounds;
4913 
4914  size_t
4915  number_threads = 1;
4916 
4917  ssize_t
4918  i,
4919  y;
4920 
4921  assert(image != (Image *) NULL);
4922  assert(image->signature == MagickCoreSignature);
4923  assert(draw_info != (DrawInfo *) NULL);
4924  assert(draw_info->signature == MagickCoreSignature);
4925  assert(primitive_info != (PrimitiveInfo *) NULL);
4926  if (IsEventLogging() != MagickFalse)
4927  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
4928  if (primitive_info->coordinates <= 1)
4929  return(MagickTrue);
4930  /*
4931  Compute bounding box.
4932  */
4933  polygon_info=AcquirePolygonTLS(draw_info,primitive_info,&image->exception);
4934  if (polygon_info == (PolygonInfo **) NULL)
4935  return(MagickFalse);
4936  if (draw_info->debug != MagickFalse)
4937  (void) LogMagickEvent(DrawEvent,GetMagickModule()," begin draw-polygon");
4938  fill=(primitive_info->method == FillToBorderMethod) ||
4939  (primitive_info->method == FloodfillMethod) ? MagickTrue : MagickFalse;
4940  mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
4941  bounds=polygon_info[0]->edges[0].bounds;
4942  artifact=GetImageArtifact(image,"draw:render-bounding-rectangles");
4943  if (IsStringTrue(artifact) != MagickFalse)
4944  (void) DrawBoundingRectangles(image,draw_info,polygon_info[0]);
4945  for (i=1; i < (ssize_t) polygon_info[0]->number_edges; i++)
4946  {
4947  p=polygon_info[0]->edges+i;
4948  if (p->bounds.x1 < bounds.x1)
4949  bounds.x1=p->bounds.x1;
4950  if (p->bounds.y1 < bounds.y1)
4951  bounds.y1=p->bounds.y1;
4952  if (p->bounds.x2 > bounds.x2)
4953  bounds.x2=p->bounds.x2;
4954  if (p->bounds.y2 > bounds.y2)
4955  bounds.y2=p->bounds.y2;
4956  }
4957  bounds.x1-=(mid+1.0);
4958  bounds.y1-=(mid+1.0);
4959  bounds.x2+=(mid+1.0);
4960  bounds.y2+=(mid+1.0);
4961  if ((bounds.x1 >= (double) image->columns) ||
4962  (bounds.y1 >= (double) image->rows) ||
4963  (bounds.x2 <= 0.0) || (bounds.y2 <= 0.0))
4964  {
4965  polygon_info=DestroyPolygonTLS(polygon_info);
4966  return(MagickTrue); /* virtual polygon */
4967  }
4968  bounds.x1=bounds.x1 < 0.0 ? 0.0 : bounds.x1 >= (double) image->columns-1.0 ?
4969  (double) image->columns-1.0 : bounds.x1;
4970  bounds.y1=bounds.y1 < 0.0 ? 0.0 : bounds.y1 >= (double) image->rows-1.0 ?
4971  (double) image->rows-1.0 : bounds.y1;
4972  bounds.x2=bounds.x2 < 0.0 ? 0.0 : bounds.x2 >= (double) image->columns-1.0 ?
4973  (double) image->columns-1.0 : bounds.x2;
4974  bounds.y2=bounds.y2 < 0.0 ? 0.0 : bounds.y2 >= (double) image->rows-1.0 ?
4975  (double) image->rows-1.0 : bounds.y2;
4976  poly_extent.x1=CastDoubleToLong(ceil(bounds.x1-0.5));
4977  poly_extent.y1=CastDoubleToLong(ceil(bounds.y1-0.5));
4978  poly_extent.x2=CastDoubleToLong(floor(bounds.x2+0.5));
4979  poly_extent.y2=CastDoubleToLong(floor(bounds.y2+0.5));
4980  number_threads=GetMagickNumberThreads(image,image,poly_extent.y2-
4981  poly_extent.y1+1,1);
4982  status=AcquirePolygonEdgesTLS(polygon_info,number_threads,&image->exception);
4983  if (status == MagickFalse)
4984  {
4985  polygon_info=DestroyPolygonTLS(polygon_info);
4986  return(status);
4987  }
4988  status=MagickTrue;
4989  exception=(&image->exception);
4990  image_view=AcquireAuthenticCacheView(image,exception);
4991  if ((primitive_info->coordinates == 1) ||
4992  (polygon_info[0]->number_edges == 0))
4993  {
4994  /*
4995  Draw point.
4996  */
4997 #if defined(MAGICKCORE_OPENMP_SUPPORT)
4998  #pragma omp parallel for schedule(static) shared(status) \
4999  num_threads(number_threads)
5000 #endif
5001  for (y=poly_extent.y1; y <= poly_extent.y2; y++)
5002  {
5003  MagickBooleanType
5004  sync;
5005 
5006  PixelPacket
5007  *magick_restrict q;
5008 
5009  ssize_t
5010  x;
5011 
5012  if (status == MagickFalse)
5013  continue;
5014  x=poly_extent.x1;
5015  q=GetCacheViewAuthenticPixels(image_view,x,y,(size_t) (poly_extent.x2-
5016  x+1),1,exception);
5017  if (q == (PixelPacket *) NULL)
5018  {
5019  status=MagickFalse;
5020  continue;
5021  }
5022  for ( ; x <= poly_extent.x2; x++)
5023  {
5024  if ((x == CastDoubleToLong(ceil(primitive_info->point.x-0.5))) &&
5025  (y == CastDoubleToLong(ceil(primitive_info->point.y-0.5))))
5026  (void) GetFillColor(draw_info,x-poly_extent.x1,y-poly_extent.y1,q);
5027  q++;
5028  }
5029  sync=SyncCacheViewAuthenticPixels(image_view,exception);
5030  if (sync == MagickFalse)
5031  status=MagickFalse;
5032  }
5033  image_view=DestroyCacheView(image_view);
5034  polygon_info=DestroyPolygonTLS(polygon_info);
5035  if (draw_info->debug != MagickFalse)
5036  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5037  " end draw-polygon");
5038  return(status);
5039  }
5040  /*
5041  Draw polygon or line.
5042  */
5043  poly_extent.y1=CastDoubleToLong(ceil(bounds.y1-0.5));
5044  poly_extent.y2=CastDoubleToLong(floor(bounds.y2+0.5));
5045 #if defined(MAGICKCORE_OPENMP_SUPPORT)
5046  #pragma omp parallel for schedule(static) shared(status) \
5047  num_threads(number_threads)
5048 #endif
5049  for (y=poly_extent.y1; y <= poly_extent.y2; y++)
5050  {
5051  const int
5052  id = GetOpenMPThreadId();
5053 
5054  PixelPacket
5055  fill_color,
5056  stroke_color;
5057 
5058  PixelPacket
5059  *magick_restrict q;
5060 
5061  ssize_t
5062  x;
5063 
5064  if (status == MagickFalse)
5065  continue;
5066  q=GetCacheViewAuthenticPixels(image_view,poly_extent.x1,y,(size_t)
5067  (poly_extent.x2-poly_extent.x1+1),1,exception);
5068  if (q == (PixelPacket *) NULL)
5069  {
5070  status=MagickFalse;
5071  continue;
5072  }
5073  for (x=poly_extent.x1; x <= poly_extent.x2; x++)
5074  {
5075  double
5076  fill_opacity,
5077  stroke_opacity;
5078 
5079  /*
5080  Fill and/or stroke.
5081  */
5082  fill_opacity=GetOpacityPixel(polygon_info[id],mid,fill,
5083  draw_info->fill_rule,x,y,&stroke_opacity);
5084  if (draw_info->stroke_antialias == MagickFalse)
5085  {
5086  fill_opacity=fill_opacity >= AntialiasThreshold ? 1.0 : 0.0;
5087  stroke_opacity=stroke_opacity >= AntialiasThreshold ? 1.0 : 0.0;
5088  }
5089  (void) GetFillColor(draw_info,x-poly_extent.x1,y-poly_extent.y1,
5090  &fill_color);
5091  fill_opacity=(double) (QuantumRange-fill_opacity*(QuantumRange-
5092  fill_color.opacity));
5093  MagickCompositeOver(&fill_color,(MagickRealType) fill_opacity,q,
5094  (MagickRealType) q->opacity,q);
5095  (void) GetStrokeColor(draw_info,x-poly_extent.x1,y-poly_extent.y1,
5096  &stroke_color);
5097  stroke_opacity=(double) (QuantumRange-stroke_opacity*(QuantumRange-
5098  stroke_color.opacity));
5099  MagickCompositeOver(&stroke_color,(MagickRealType) stroke_opacity,q,
5100  (MagickRealType) q->opacity,q);
5101  q++;
5102  }
5103  if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
5104  status=MagickFalse;
5105  }
5106  image_view=DestroyCacheView(image_view);
5107  polygon_info=DestroyPolygonTLS(polygon_info);
5108  if (draw_info->debug != MagickFalse)
5109  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end draw-polygon");
5110  return(status);
5111 }
5112 
5113 /*
5114 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5115 % %
5116 % %
5117 % %
5118 % D r a w P r i m i t i v e %
5119 % %
5120 % %
5121 % %
5122 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5123 %
5124 % DrawPrimitive() draws a primitive (line, rectangle, ellipse) on the image.
5125 %
5126 % The format of the DrawPrimitive method is:
5127 %
5128 % MagickBooleanType DrawPrimitive(Image *image,const DrawInfo *draw_info,
5129 % PrimitiveInfo *primitive_info)
5130 %
5131 % A description of each parameter follows:
5132 %
5133 % o image: the image.
5134 %
5135 % o draw_info: the draw info.
5136 %
5137 % o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
5138 %
5139 */
5140 static void LogPrimitiveInfo(const PrimitiveInfo *primitive_info)
5141 {
5142  const char
5143  *methods[] =
5144  {
5145  "point",
5146  "replace",
5147  "floodfill",
5148  "filltoborder",
5149  "reset",
5150  "?"
5151  };
5152 
5153  PointInfo
5154  p,
5155  q,
5156  point;
5157 
5158  ssize_t
5159  i,
5160  x;
5161 
5162  ssize_t
5163  coordinates,
5164  y;
5165 
5166  x=CastDoubleToLong(ceil(primitive_info->point.x-0.5));
5167  y=CastDoubleToLong(ceil(primitive_info->point.y-0.5));
5168  switch (primitive_info->primitive)
5169  {
5170  case PointPrimitive:
5171  {
5172  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5173  "PointPrimitive %.20g,%.20g %s",(double) x,(double) y,
5174  methods[primitive_info->method]);
5175  return;
5176  }
5177  case ColorPrimitive:
5178  {
5179  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5180  "ColorPrimitive %.20g,%.20g %s",(double) x,(double) y,
5181  methods[primitive_info->method]);
5182  return;
5183  }
5184  case MattePrimitive:
5185  {
5186  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5187  "MattePrimitive %.20g,%.20g %s",(double) x,(double) y,
5188  methods[primitive_info->method]);
5189  return;
5190  }
5191  case TextPrimitive:
5192  {
5193  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5194  "TextPrimitive %.20g,%.20g",(double) x,(double) y);
5195  return;
5196  }
5197  case ImagePrimitive:
5198  {
5199  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5200  "ImagePrimitive %.20g,%.20g",(double) x,(double) y);
5201  return;
5202  }
5203  default:
5204  break;
5205  }
5206  coordinates=0;
5207  p=primitive_info[0].point;
5208  q.x=(-1.0);
5209  q.y=(-1.0);
5210  for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
5211  {
5212  point=primitive_info[i].point;
5213  if (coordinates <= 0)
5214  {
5215  coordinates=(ssize_t) primitive_info[i].coordinates;
5216  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5217  " begin open (%.20g)",(double) coordinates);
5218  p=point;
5219  }
5220  point=primitive_info[i].point;
5221  if ((fabs(q.x-point.x) >= MagickEpsilon) ||
5222  (fabs(q.y-point.y) >= MagickEpsilon))
5223  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5224  " %.20g: %.18g,%.18g",(double) coordinates,point.x,point.y);
5225  else
5226  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5227  " %.20g: %g %g (duplicate)",(double) coordinates,point.x,point.y);
5228  q=point;
5229  coordinates--;
5230  if (coordinates > 0)
5231  continue;
5232  if ((fabs(p.x-point.x) >= MagickEpsilon) ||
5233  (fabs(p.y-point.y) >= MagickEpsilon))
5234  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end last (%.20g)",
5235  (double) coordinates);
5236  else
5237  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end open (%.20g)",
5238  (double) coordinates);
5239  }
5240 }
5241 
5242 MagickExport MagickBooleanType DrawPrimitive(Image *image,
5243  const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
5244 {
5245  CacheView
5246  *image_view;
5247 
5249  *exception;
5250 
5251  MagickStatusType
5252  status;
5253 
5254  ssize_t
5255  i,
5256  x;
5257 
5258  ssize_t
5259  y;
5260 
5261  if (draw_info->debug != MagickFalse)
5262  {
5263  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5264  " begin draw-primitive");
5265  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5266  " affine: %g,%g,%g,%g,%g,%g",draw_info->affine.sx,
5267  draw_info->affine.rx,draw_info->affine.ry,draw_info->affine.sy,
5268  draw_info->affine.tx,draw_info->affine.ty);
5269  }
5270  exception=(&image->exception);
5271  status=MagickTrue;
5272  if ((IsGrayColorspace(image->colorspace) != MagickFalse) &&
5273  ((IsPixelGray(&draw_info->fill) == MagickFalse) ||
5274  (IsPixelGray(&draw_info->stroke) == MagickFalse)))
5275  status=SetImageColorspace(image,sRGBColorspace);
5276  if (draw_info->compliance == SVGCompliance)
5277  {
5278  status&=SetImageClipMask(image,draw_info->clipping_mask);
5279  status&=SetImageMask(image,draw_info->composite_mask);
5280  }
5281  x=CastDoubleToLong(ceil(primitive_info->point.x-0.5));
5282  y=CastDoubleToLong(ceil(primitive_info->point.y-0.5));
5283  image_view=AcquireAuthenticCacheView(image,exception);
5284  switch (primitive_info->primitive)
5285  {
5286  case ColorPrimitive:
5287  {
5288  switch (primitive_info->method)
5289  {
5290  case PointMethod:
5291  default:
5292  {
5293  PixelPacket
5294  *q;
5295 
5296  q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
5297  if (q == (PixelPacket *) NULL)
5298  break;
5299  (void) GetFillColor(draw_info,x,y,q);
5300  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5301  break;
5302  }
5303  case ReplaceMethod:
5304  {
5305  PixelPacket
5306  target;
5307 
5308  status&=GetOneCacheViewVirtualPixel(image_view,x,y,&target,exception);
5309  for (y=0; y < (ssize_t) image->rows; y++)
5310  {
5311  PixelPacket
5312  *magick_restrict q;
5313 
5314  q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
5315  exception);
5316  if (q == (PixelPacket *) NULL)
5317  break;
5318  for (x=0; x < (ssize_t) image->columns; x++)
5319  {
5320  if (IsColorSimilar(image,q,&target) == MagickFalse)
5321  {
5322  q++;
5323  continue;
5324  }
5325  (void) GetFillColor(draw_info,x,y,q);
5326  q++;
5327  }
5328  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5329  if (status == MagickFalse)
5330  break;
5331  }
5332  break;
5333  }
5334  case FloodfillMethod:
5335  case FillToBorderMethod:
5336  {
5338  target;
5339 
5340  status&=GetOneVirtualMagickPixel(image,x,y,&target,exception);
5341  if (primitive_info->method == FillToBorderMethod)
5342  {
5343  target.red=(MagickRealType) draw_info->border_color.red;
5344  target.green=(MagickRealType) draw_info->border_color.green;
5345  target.blue=(MagickRealType) draw_info->border_color.blue;
5346  }
5347  status&=FloodfillPaintImage(image,DefaultChannels,draw_info,&target,x,
5348  y,primitive_info->method == FloodfillMethod ? MagickFalse :
5349  MagickTrue);
5350  break;
5351  }
5352  case ResetMethod:
5353  {
5354  for (y=0; y < (ssize_t) image->rows; y++)
5355  {
5356  PixelPacket
5357  *magick_restrict q;
5358 
5359  ssize_t
5360  x;
5361 
5362  q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
5363  exception);
5364  if (q == (PixelPacket *) NULL)
5365  break;
5366  for (x=0; x < (ssize_t) image->columns; x++)
5367  {
5368  (void) GetFillColor(draw_info,x,y,q);
5369  q++;
5370  }
5371  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5372  if (status == MagickFalse)
5373  break;
5374  }
5375  break;
5376  }
5377  }
5378  break;
5379  }
5380  case MattePrimitive:
5381  {
5382  if (image->matte == MagickFalse)
5383  status&=SetImageAlphaChannel(image,OpaqueAlphaChannel);
5384  switch (primitive_info->method)
5385  {
5386  case PointMethod:
5387  default:
5388  {
5389  PixelPacket
5390  pixel;
5391 
5392  PixelPacket
5393  *q;
5394 
5395  q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
5396  if (q == (PixelPacket *) NULL)
5397  break;
5398  (void) GetFillColor(draw_info,x,y,&pixel);
5399  SetPixelOpacity(q,pixel.opacity);
5400  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5401  break;
5402  }
5403  case ReplaceMethod:
5404  {
5405  PixelPacket
5406  pixel,
5407  target;
5408 
5409  status&=GetOneCacheViewVirtualPixel(image_view,x,y,&target,exception);
5410  for (y=0; y < (ssize_t) image->rows; y++)
5411  {
5412  PixelPacket
5413  *magick_restrict q;
5414 
5415  ssize_t
5416  x;
5417 
5418  q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
5419  exception);
5420  if (q == (PixelPacket *) NULL)
5421  break;
5422  for (x=0; x < (ssize_t) image->columns; x++)
5423  {
5424  if (IsColorSimilar(image,q,&target) == MagickFalse)
5425  {
5426  q++;
5427  continue;
5428  }
5429  (void) GetFillColor(draw_info,x,y,&pixel);
5430  SetPixelOpacity(q,pixel.opacity);
5431  q++;
5432  }
5433  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5434  if (status == MagickFalse)
5435  break;
5436  }
5437  break;
5438  }
5439  case FloodfillMethod:
5440  case FillToBorderMethod:
5441  {
5443  target;
5444 
5445  status&=GetOneVirtualMagickPixel(image,x,y,&target,exception);
5446  if (primitive_info->method == FillToBorderMethod)
5447  {
5448  target.red=(MagickRealType) draw_info->border_color.red;
5449  target.green=(MagickRealType) draw_info->border_color.green;
5450  target.blue=(MagickRealType) draw_info->border_color.blue;
5451  }
5452  status&=FloodfillPaintImage(image,OpacityChannel,draw_info,&target,x,
5453  y,primitive_info->method == FloodfillMethod ? MagickFalse :
5454  MagickTrue);
5455  break;
5456  }
5457  case ResetMethod:
5458  {
5459  PixelPacket
5460  pixel;
5461 
5462  for (y=0; y < (ssize_t) image->rows; y++)
5463  {
5464  PixelPacket
5465  *magick_restrict q;
5466 
5467  ssize_t
5468  x;
5469 
5470  q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
5471  exception);
5472  if (q == (PixelPacket *) NULL)
5473  break;
5474  for (x=0; x < (ssize_t) image->columns; x++)
5475  {
5476  (void) GetFillColor(draw_info,x,y,&pixel);
5477  SetPixelOpacity(q,pixel.opacity);
5478  q++;
5479  }
5480  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5481  if (status == MagickFalse)
5482  break;
5483  }
5484  break;
5485  }
5486  }
5487  break;
5488  }
5489  case ImagePrimitive:
5490  {
5491  AffineMatrix
5492  affine;
5493 
5494  char
5495  composite_geometry[MaxTextExtent];
5496 
5497  Image
5498  *composite_image,
5499  *composite_images;
5500 
5501  ImageInfo
5502  *clone_info;
5503 
5505  geometry;
5506 
5507  ssize_t
5508  x1,
5509  y1;
5510 
5511  if (primitive_info->text == (char *) NULL)
5512  break;
5513  clone_info=AcquireImageInfo();
5514  composite_images=(Image *) NULL;
5515  if (LocaleNCompare(primitive_info->text,"data:",5) == 0)
5516  composite_images=ReadInlineImage(clone_info,primitive_info->text,
5517  &image->exception);
5518  else
5519  if (*primitive_info->text != '\0')
5520  {
5521  const MagickInfo
5522  *magick_info;
5523 
5524  MagickStatusType
5525  path_status;
5526 
5527  struct stat
5528  attributes;
5529 
5530  /*
5531  Read composite image.
5532  */
5533  (void) CopyMagickString(clone_info->filename,primitive_info->text,
5534  MagickPathExtent);
5535  (void) SetImageInfo(clone_info,1,exception);
5536  magick_info=GetMagickInfo(clone_info->magick,exception);
5537  if ((magick_info != (const MagickInfo*) NULL) &&
5538  (LocaleCompare(magick_info->magick_module,"SVG") == 0))
5539  {
5540  (void) ThrowMagickException(exception,GetMagickModule(),
5541  CorruptImageError,"ImageTypeNotSupported","`%s'",
5542  clone_info->filename);
5543  clone_info=DestroyImageInfo(clone_info);
5544  break;
5545  }
5546  (void) CopyMagickString(clone_info->filename,primitive_info->text,
5547  MagickPathExtent);
5548  if (clone_info->size != (char *) NULL)
5549  clone_info->size=DestroyString(clone_info->size);
5550  if (clone_info->extract != (char *) NULL)
5551  clone_info->extract=DestroyString(clone_info->extract);
5552  path_status=GetPathAttributes(clone_info->filename,&attributes);
5553  if (path_status != MagickFalse)
5554  {
5555  if (S_ISCHR(attributes.st_mode) == 0)
5556  composite_images=ReadImage(clone_info,exception);
5557  else
5558  (void) ThrowMagickException(exception,GetMagickModule(),
5559  FileOpenError,"UnableToOpenFile","`%s'",
5560  clone_info->filename);
5561  }
5562  else
5563  if ((LocaleCompare(clone_info->magick,"ftp") != 0) &&
5564  (LocaleCompare(clone_info->magick,"http") != 0) &&
5565  (LocaleCompare(clone_info->magick,"https") != 0) &&
5566  (LocaleCompare(clone_info->magick,"vid") != 0))
5567  composite_images=ReadImage(clone_info,exception);
5568  else
5569  (void) ThrowMagickException(exception,GetMagickModule(),
5570  FileOpenError,"UnableToOpenFile","`%s'",clone_info->filename);
5571  }
5572  clone_info=DestroyImageInfo(clone_info);
5573  if (composite_images == (Image *) NULL)
5574  {
5575  status=0;
5576  break;
5577  }
5578  composite_image=RemoveFirstImageFromList(&composite_images);
5579  composite_images=DestroyImageList(composite_images);
5580  (void) SetImageProgressMonitor(composite_image,(MagickProgressMonitor)
5581  NULL,(void *) NULL);
5582  x1=CastDoubleToLong(ceil(primitive_info[1].point.x-0.5));
5583  y1=CastDoubleToLong(ceil(primitive_info[1].point.y-0.5));
5584  if (((x1 != 0L) && (x1 != (ssize_t) composite_image->columns)) ||
5585  ((y1 != 0L) && (y1 != (ssize_t) composite_image->rows)))
5586  {
5587  char
5588  geometry[MaxTextExtent];
5589 
5590  /*
5591  Resize image.
5592  */
5593  (void) FormatLocaleString(geometry,MaxTextExtent,"%gx%g!",
5594  primitive_info[1].point.x,primitive_info[1].point.y);
5595  composite_image->filter=image->filter;
5596  status&=TransformImage(&composite_image,(char *) NULL,geometry);
5597  }
5598  if (composite_image->matte == MagickFalse)
5599  status&=SetImageAlphaChannel(composite_image,OpaqueAlphaChannel);
5600  if (draw_info->opacity != OpaqueOpacity)
5601  status&=SetImageOpacity(composite_image,draw_info->opacity);
5602  SetGeometry(image,&geometry);
5603  image->gravity=draw_info->gravity;
5604  geometry.x=x;
5605  geometry.y=y;
5606  (void) FormatLocaleString(composite_geometry,MaxTextExtent,
5607  "%.20gx%.20g%+.20g%+.20g",(double) composite_image->columns,(double)
5608  composite_image->rows,(double) geometry.x,(double) geometry.y);
5609  (void) ParseGravityGeometry(image,composite_geometry,&geometry,
5610  &image->exception);
5611  affine=draw_info->affine;
5612  affine.tx=(double) geometry.x;
5613  affine.ty=(double) geometry.y;
5614  composite_image->interpolate=image->interpolate;
5615  if ((draw_info->compose == OverCompositeOp) ||
5616  (draw_info->compose == SrcOverCompositeOp))
5617  status&=DrawAffineImage(image,composite_image,&affine);
5618  else
5619  status&=CompositeImage(image,draw_info->compose,composite_image,
5620  geometry.x,geometry.y);
5621  composite_image=DestroyImage(composite_image);
5622  break;
5623  }
5624  case PointPrimitive:
5625  {
5626  PixelPacket
5627  fill_color;
5628 
5629  PixelPacket
5630  *q;
5631 
5632  if ((y < 0) || (y >= (ssize_t) image->rows))
5633  break;
5634  if ((x < 0) || (x >= (ssize_t) image->columns))
5635  break;
5636  q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
5637  if (q == (PixelPacket *) NULL)
5638  break;
5639  (void) GetFillColor(draw_info,x,y,&fill_color);
5640  MagickCompositeOver(&fill_color,(MagickRealType) fill_color.opacity,q,
5641  (MagickRealType) q->opacity,q);
5642  status&=SyncCacheViewAuthenticPixels(image_view,exception);
5643  break;
5644  }
5645  case TextPrimitive:
5646  {
5647  char
5648  geometry[MaxTextExtent];
5649 
5650  DrawInfo
5651  *clone_info;
5652 
5653  if (primitive_info->text == (char *) NULL)
5654  break;
5655  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
5656  (void) CloneString(&clone_info->text,primitive_info->text);
5657  (void) FormatLocaleString(geometry,MaxTextExtent,"%+f%+f",
5658  primitive_info->point.x,primitive_info->point.y);
5659  (void) CloneString(&clone_info->geometry,geometry);
5660  status&=AnnotateImage(image,clone_info);
5661  clone_info=DestroyDrawInfo(clone_info);
5662  break;
5663  }
5664  default:
5665  {
5666  double
5667  mid,
5668  scale;
5669 
5670  DrawInfo
5671  *clone_info;
5672 
5673  if (IsEventLogging() != MagickFalse)
5674  LogPrimitiveInfo(primitive_info);
5675  scale=ExpandAffine(&draw_info->affine);
5676  if ((draw_info->dash_pattern != (double *) NULL) &&
5677  (fabs(draw_info->dash_pattern[0]) >= MagickEpsilon) &&
5678  (fabs(scale*draw_info->stroke_width) >= MagickEpsilon) &&
5679  (draw_info->stroke.opacity != (Quantum) TransparentOpacity))
5680  {
5681  /*
5682  Draw dash polygon.
5683  */
5684  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
5685  clone_info->stroke_width=0.0;
5686  clone_info->stroke.opacity=(Quantum) TransparentOpacity;
5687  status&=DrawPolygonPrimitive(image,clone_info,primitive_info);
5688  clone_info=DestroyDrawInfo(clone_info);
5689  if (status != MagickFalse)
5690  status&=DrawDashPolygon(draw_info,primitive_info,image);
5691  break;
5692  }
5693  mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
5694  if ((mid > 1.0) &&
5695  ((draw_info->stroke.opacity != (Quantum) TransparentOpacity) ||
5696  (draw_info->stroke_pattern != (Image *) NULL)))
5697  {
5698  double
5699  x,
5700  y;
5701 
5702  MagickBooleanType
5703  closed_path;
5704 
5705  /*
5706  Draw strokes while respecting line cap/join attributes.
5707  */
5708  closed_path=primitive_info[0].closed_subpath;
5709  i=(ssize_t) primitive_info[0].coordinates;
5710  x=fabs(primitive_info[i-1].point.x-primitive_info[0].point.x);
5711  y=fabs(primitive_info[i-1].point.y-primitive_info[0].point.y);
5712  if ((x < MagickEpsilon) && (y < MagickEpsilon))
5713  closed_path=MagickTrue;
5714  if ((((draw_info->linecap == RoundCap) ||
5715  (closed_path != MagickFalse)) &&
5716  (draw_info->linejoin == RoundJoin)) ||
5717  (primitive_info[i].primitive != UndefinedPrimitive))
5718  {
5719  status&=DrawPolygonPrimitive(image,draw_info,primitive_info);
5720  break;
5721  }
5722  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
5723  clone_info->stroke_width=0.0;
5724  clone_info->stroke.opacity=(Quantum) TransparentOpacity;
5725  status&=DrawPolygonPrimitive(image,clone_info,primitive_info);
5726  clone_info=DestroyDrawInfo(clone_info);
5727  if (status != MagickFalse)
5728  status&=DrawStrokePolygon(image,draw_info,primitive_info);
5729  break;
5730  }
5731  status&=DrawPolygonPrimitive(image,draw_info,primitive_info);
5732  break;
5733  }
5734  }
5735  image_view=DestroyCacheView(image_view);
5736  if (draw_info->compliance == SVGCompliance)
5737  {
5738  status&=SetImageClipMask(image,(Image *) NULL);
5739  status&=SetImageMask(image,(Image *) NULL);
5740  }
5741  if (draw_info->debug != MagickFalse)
5742  (void) LogMagickEvent(DrawEvent,GetMagickModule()," end draw-primitive");
5743  return(status != 0 ? MagickTrue : MagickFalse);
5744 }
5745 
5746 /*
5747 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5748 % %
5749 % %
5750 % %
5751 + D r a w S t r o k e P o l y g o n %
5752 % %
5753 % %
5754 % %
5755 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5756 %
5757 % DrawStrokePolygon() draws a stroked polygon (line, rectangle, ellipse) on
5758 % the image while respecting the line cap and join attributes.
5759 %
5760 % The format of the DrawStrokePolygon method is:
5761 %
5762 % MagickBooleanType DrawStrokePolygon(Image *image,
5763 % const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
5764 %
5765 % A description of each parameter follows:
5766 %
5767 % o image: the image.
5768 %
5769 % o draw_info: the draw info.
5770 %
5771 % o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
5772 %
5773 %
5774 */
5775 
5776 static MagickBooleanType DrawRoundLinecap(Image *image,
5777  const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
5778 {
5780  linecap[5];
5781 
5782  ssize_t
5783  i;
5784 
5785  for (i=0; i < 4; i++)
5786  linecap[i]=(*primitive_info);
5787  linecap[0].coordinates=4;
5788  linecap[1].point.x+=2.0*MagickEpsilon;
5789  linecap[2].point.x+=2.0*MagickEpsilon;
5790  linecap[2].point.y+=2.0*MagickEpsilon;
5791  linecap[3].point.y+=2.0*MagickEpsilon;
5792  linecap[4].primitive=UndefinedPrimitive;
5793  return(DrawPolygonPrimitive(image,draw_info,linecap));
5794 }
5795 
5796 static MagickBooleanType DrawStrokePolygon(Image *image,
5797  const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
5798 {
5799  DrawInfo
5800  *clone_info;
5801 
5802  MagickBooleanType
5803  closed_path;
5804 
5805  MagickStatusType
5806  status;
5807 
5809  *stroke_polygon;
5810 
5811  const PrimitiveInfo
5812  *p,
5813  *q;
5814 
5815  /*
5816  Draw stroked polygon.
5817  */
5818  if (draw_info->debug != MagickFalse)
5819  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5820  " begin draw-stroke-polygon");
5821  clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
5822  clone_info->fill=draw_info->stroke;
5823  if (clone_info->fill_pattern != (Image *) NULL)
5824  clone_info->fill_pattern=DestroyImage(clone_info->fill_pattern);
5825  if (clone_info->stroke_pattern != (Image *) NULL)
5826  clone_info->fill_pattern=CloneImage(clone_info->stroke_pattern,0,0,
5827  MagickTrue,&clone_info->stroke_pattern->exception);
5828  clone_info->stroke.opacity=(Quantum) TransparentOpacity;
5829  clone_info->stroke_width=0.0;
5830  clone_info->fill_rule=NonZeroRule;
5831  status=MagickTrue;
5832  for (p=primitive_info; p->primitive != UndefinedPrimitive; p+=p->coordinates)
5833  {
5834  if (p->coordinates == 1)
5835  continue;
5836  stroke_polygon=TraceStrokePolygon(draw_info,p,&image->exception);
5837  if (stroke_polygon == (PrimitiveInfo *) NULL)
5838  {
5839  status=0;
5840  break;
5841  }
5842  status&=DrawPolygonPrimitive(image,clone_info,stroke_polygon);
5843  stroke_polygon=(PrimitiveInfo *) RelinquishMagickMemory(stroke_polygon);
5844  if (status == 0)
5845  break;
5846  q=p+p->coordinates-1;
5847  closed_path=p->closed_subpath;
5848  if ((draw_info->linecap == RoundCap) && (closed_path == MagickFalse))
5849  {
5850  status&=DrawRoundLinecap(image,draw_info,p);
5851  status&=DrawRoundLinecap(image,draw_info,q);
5852  }
5853  }
5854  clone_info=DestroyDrawInfo(clone_info);
5855  if (draw_info->debug != MagickFalse)
5856  (void) LogMagickEvent(DrawEvent,GetMagickModule(),
5857  " end draw-stroke-polygon");
5858  return(status != 0 ? MagickTrue : MagickFalse);
5859 }
5860 
5861 /*
5862 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5863 % %
5864 % %
5865 % %
5866 % G e t A f f i n e M a t r i x %
5867 % %
5868 % %
5869 % %
5870 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5871 %
5872 % GetAffineMatrix() returns an AffineMatrix initialized to the identity
5873 % matrix.
5874 %
5875 % The format of the GetAffineMatrix method is:
5876 %
5877 % void GetAffineMatrix(AffineMatrix *affine_matrix)
5878 %
5879 % A description of each parameter follows:
5880 %
5881 % o affine_matrix: the affine matrix.
5882 %
5883 */
5884 MagickExport void GetAffineMatrix(AffineMatrix *affine_matrix)
5885 {
5886  assert(affine_matrix != (AffineMatrix *) NULL);
5887  if (IsEventLogging() != MagickFalse)
5888  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
5889  (void) memset(affine_matrix,0,sizeof(*affine_matrix));
5890  affine_matrix->sx=1.0;
5891  affine_matrix->sy=1.0;
5892 }
5893 
5894 /*
5895 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5896 % %
5897 % %
5898 % %
5899 + G e t D r a w I n f o %
5900 % %
5901 % %
5902 % %
5903 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5904 %
5905 % GetDrawInfo() initializes draw_info to default values from image_info.
5906 %
5907 % The format of the GetDrawInfo method is:
5908 %
5909 % void GetDrawInfo(const ImageInfo *image_info,DrawInfo *draw_info)
5910 %
5911 % A description of each parameter follows:
5912 %
5913 % o image_info: the image info..
5914 %
5915 % o draw_info: the draw info.
5916 %
5917 */
5918 MagickExport void GetDrawInfo(const ImageInfo *image_info,DrawInfo *draw_info)
5919 {
5920  char
5921  *next_token;
5922 
5923  const char
5924  *option;
5925 
5927  *exception;
5928 
5929  /*
5930  Initialize draw attributes.
5931  */
5932  assert(draw_info != (DrawInfo *) NULL);
5933  if (IsEventLogging() != MagickFalse)
5934  (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
5935  (void) memset(draw_info,0,sizeof(*draw_info));
5936  draw_info->image_info=CloneImageInfo(image_info);
5937  GetAffineMatrix(&draw_info->affine);
5938  exception=AcquireExceptionInfo();
5939  (void) QueryColorDatabase("#000F",&draw_info->fill,exception);
5940  (void) QueryColorDatabase("#FFF0",&draw_info->stroke,exception);
5941  draw_info->stroke_antialias=draw_info->image_info->antialias;
5942  draw_info->stroke_width=1.0;
5943  draw_info->fill_rule=EvenOddRule;
5944  draw_info->opacity=OpaqueOpacity;
5945  draw_info->fill_opacity=OpaqueOpacity;
5946  draw_info->stroke_opacity=OpaqueOpacity;
5947  draw_info->linecap=ButtCap;
5948  draw_info->linejoin=MiterJoin;
5949  draw_info->miterlimit=10;
5950  draw_info->decorate=NoDecoration;
5951  if (draw_info->image_info->font != (char *) NULL)
5952  draw_info->font=AcquireString(draw_info->image_info->font);
5953  if (draw_info->image_info->density != (char *) NULL)
5954  draw_info->density=AcquireString(draw_info->image_info->density);
5955  draw_info->text_antialias=draw_info->image_info->antialias;
5956  draw_info->pointsize=12.0;
5957  if (fabs(draw_info->image_info->pointsize) >= MagickEpsilon)
5958  draw_info->pointsize=draw_info->image_info->pointsize;
5959  draw_info->undercolor.opacity=(Quantum) TransparentOpacity;
5960  draw_info->border_color=draw_info->image_info->border_color;
5961  draw_info->compose=OverCompositeOp;
5962  if (draw_info->image_info->server_name != (char *) NULL)
5963  draw_info->server_name=AcquireString(draw_info->image_info->server_name);
5964  draw_info->render=MagickTrue;
5965  draw_info->clip_path=MagickFalse;
5966  draw_info->debug=(GetLogEventMask() & (DrawEvent | AnnotateEvent)) != 0 ?
5967  MagickTrue : MagickFalse;
5968  option=GetImageOption(draw_info->image_info,"direction");
5969  if (option != (const char *) NULL)
5970  draw_info->direction=(DirectionType) ParseCommandOption(
5971  MagickDirectionOptions,MagickFalse,option);
5972  else
5973  draw_info->direction=UndefinedDirection;
5974  option=GetImageOption(draw_info->image_info,"encoding");
5975  if (option != (const char *) NULL)
5976  (void) CloneString(&draw_info->encoding,option);
5977  option=GetImageOption(draw_info->image_info,"family");
5978  if (option != (const char *) NULL)
5979  (void) CloneString(&draw_info->family,option);
5980  option=GetImageOption(draw_info->image_info,"fill");
5981  if (option != (const char *) NULL)
5982  (void) QueryColorDatabase(option,&draw_info->fill,exception);
5983  option=GetImageOption(draw_info->image_info,"gravity");
5984  if (option != (const char *) NULL)
5985  draw_info->gravity=(GravityType) ParseCommandOption(MagickGravityOptions,
5986  MagickFalse,option);
5987  option=GetImageOption(draw_info->image_info,"interline-spacing");
5988  if (option != (const char *) NULL)
5989  draw_info->interline_spacing=GetDrawValue(option,&next_token);
5990  option=GetImageOption(draw_info->image_info,"interword-spacing");
5991  if (option != (const char *) NULL)
5992  draw_info->interword_spacing=GetDrawValue(option,&next_token);
5993  option=GetImageOption(draw_info->image_info,"kerning");
5994  if (option != (const char *) NULL)
5995  draw_info->kerning=GetDrawValue(option,&next_token);
5996  option=GetImageOption(draw_info->image_info,"stroke");
5997  if (option != (const char *) NULL)
5998  (void) QueryColorDatabase(option,&draw_info->stroke,exception);
5999  option=GetImageOption(draw_info->image_info,"strokewidth");
6000  if (option != (const char *) NULL)
6001  draw_info->stroke_width=GetDrawValue(option,&next_token);
6002  option=GetImageOption(draw_info->image_info,"style");
6003  if (option != (const char *) NULL)
6004  draw_info->style=(StyleType) ParseCommandOption(MagickStyleOptions,
6005  MagickFalse,option);
6006  option=GetImageOption(draw_info->image_info,"undercolor");
6007  if (option != (const char *) NULL)
6008  (void) QueryColorDatabase(option,&draw_info->undercolor,exception);
6009  option=GetImageOption(draw_info->image_info,"weight");
6010  if (option != (const char *) NULL)
6011  {
6012  ssize_t
6013  weight;
6014 
6015  weight=ParseCommandOption(MagickWeightOptions,MagickFalse,option);
6016  if (weight == -1)
6017  weight=(ssize_t) StringToUnsignedLong(option);
6018  draw_info->weight=(size_t) weight;
6019  }
6020  exception=DestroyExceptionInfo(exception);
6021  draw_info->signature=MagickCoreSignature;
6022 }
6023 
6024 /*
6025 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6026 % %
6027 % %
6028 % %
6029 + P e r m u t a t e %
6030 % %
6031 % %
6032 % %
6033 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6034 %
6035 % Permutate() returns the permutation of the (n,k).
6036 %
6037 % The format of the Permutate method is:
6038 %
6039 % void Permutate(ssize_t n,ssize_t k)
6040 %
6041 % A description of each parameter follows:
6042 %
6043 % o n:
6044 %
6045 % o k:
6046 %
6047 %
6048 */
6049 static inline double Permutate(const ssize_t n,const ssize_t k)
6050 {
6051  double
6052  r;
6053 
6054  ssize_t
6055  i;
6056 
6057  r=1.0;
6058  for (i=k+1; i <= n; i++)
6059  r*=i;
6060  for (i=1; i <= (n-k); i++)
6061  r/=i;
6062  return(r);
6063 }
6064 
6065 /*
6066 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6067 % %
6068 % %
6069 % %
6070 + T r a c e P r i m i t i v e %
6071 % %
6072 % %
6073 % %
6074 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6075 %
6076 % TracePrimitive is a collection of methods for generating graphic
6077 % primitives such as arcs, ellipses, paths, etc.
6078 %
6079 */
6080 
6081 static MagickBooleanType TraceArc(MVGInfo *mvg_info,const PointInfo start,
6082  const PointInfo end,const PointInfo degrees)
6083 {
6084  PointInfo
6085  center,
6086  radius;
6087 
6088  center.x=0.5*(end.x+start.x);
6089  center.y=0.5*(end.y+start.y);
6090  radius.x=fabs(center.x-start.x);
6091  radius.y=fabs(center.y-start.y);
6092  return(TraceEllipse(mvg_info,center,radius,degrees));
6093 }
6094 
6095 static MagickBooleanType TraceArcPath(MVGInfo *mvg_info,const PointInfo start,
6096  const PointInfo end,const PointInfo arc,const double angle,
6097  const MagickBooleanType large_arc,const MagickBooleanType sweep)
6098 {
6099  double
6100  alpha,
6101  beta,
6102  delta,
6103  factor,
6104  gamma,
6105  theta;
6106 
6107  MagickStatusType
6108  status;
6109 
6110  PointInfo
6111  center,
6112  points[3],
6113  radii;
6114 
6115  double
6116  cosine,
6117  sine;
6118 
6120  *primitive_info;
6121 
6123  *p;
6124 
6125  ssize_t
6126  i;
6127 
6128  size_t
6129  arc_segments;
6130 
6131  ssize_t
6132  offset;
6133 
6134  offset=mvg_info->offset;
6135  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6136  primitive_info->coordinates=0;
6137  if ((fabs(start.x-end.x) < MagickEpsilon) &&
6138  (fabs(start.y-end.y) < MagickEpsilon))
6139  return(TracePoint(primitive_info,end));
6140  radii.x=fabs(arc.x);
6141  radii.y=fabs(arc.y);
6142  if ((radii.x < MagickEpsilon) || (radii.y < MagickEpsilon))
6143  return(TraceLine(primitive_info,start,end));
6144  cosine=cos(DegreesToRadians(fmod((double) angle,360.0)));
6145  sine=sin(DegreesToRadians(fmod((double) angle,360.0)));
6146  center.x=(double) (cosine*(end.x-start.x)/2+sine*(end.y-start.y)/2);
6147  center.y=(double) (cosine*(end.y-start.y)/2-sine*(end.x-start.x)/2);
6148  delta=(center.x*center.x)/(radii.x*radii.x)+(center.y*center.y)/
6149  (radii.y*radii.y);
6150  if (delta < MagickEpsilon)
6151  return(TraceLine(primitive_info,start,end));
6152  if (delta > 1.0)
6153  {
6154  radii.x*=sqrt((double) delta);
6155  radii.y*=sqrt((double) delta);
6156  }
6157  points[0].x=(double) (cosine*start.x/radii.x+sine*start.y/radii.x);
6158  points[0].y=(double) (cosine*start.y/radii.y-sine*start.x/radii.y);
6159  points[1].x=(double) (cosine*end.x/radii.x+sine*end.y/radii.x);
6160  points[1].y=(double) (cosine*end.y/radii.y-sine*end.x/radii.y);
6161  alpha=points[1].x-points[0].x;
6162  beta=points[1].y-points[0].y;
6163  if (fabs(alpha*alpha+beta*beta) < MagickEpsilon)
6164  return(TraceLine(primitive_info,start,end));
6165  factor=PerceptibleReciprocal(alpha*alpha+beta*beta)-0.25;
6166  if (factor <= 0.0)
6167  factor=0.0;
6168  else
6169  {
6170  factor=sqrt((double) factor);
6171  if (sweep == large_arc)
6172  factor=(-factor);
6173  }
6174  center.x=(double) ((points[0].x+points[1].x)/2-factor*beta);
6175  center.y=(double) ((points[0].y+points[1].y)/2+factor*alpha);
6176  alpha=atan2(points[0].y-center.y,points[0].x-center.x);
6177  theta=atan2(points[1].y-center.y,points[1].x-center.x)-alpha;
6178  if ((theta < 0.0) && (sweep != MagickFalse))
6179  theta+=2.0*MagickPI;
6180  else
6181  if ((theta > 0.0) && (sweep == MagickFalse))
6182  theta-=2.0*MagickPI;
6183  arc_segments=(size_t) CastDoubleToLong(ceil(fabs((double) (theta/(0.5*
6184  MagickPI+MagickEpsilon)))));
6185  p=primitive_info;
6186  status=MagickTrue;
6187  for (i=0; i < (ssize_t) arc_segments; i++)
6188  {
6189  beta=0.5*((alpha+(i+1)*theta/arc_segments)-(alpha+i*theta/arc_segments));
6190  gamma=(8.0/3.0)*sin(fmod((double) (0.5*beta),DegreesToRadians(360.0)))*
6191  sin(fmod((double) (0.5*beta),DegreesToRadians(360.0)))/
6192  sin(fmod((double) beta,DegreesToRadians(360.0)));
6193  points[0].x=(double) (center.x+cos(fmod((double) (alpha+(double) i*theta/
6194  arc_segments),DegreesToRadians(360.0)))-gamma*sin(fmod((double) (alpha+
6195  (double) i*theta/arc_segments),DegreesToRadians(360.0))));
6196  points[0].y=(double) (center.y+sin(fmod((double) (alpha+(double) i*theta/
6197  arc_segments),DegreesToRadians(360.0)))+gamma*cos(fmod((double) (alpha+
6198  (double) i*theta/arc_segments),DegreesToRadians(360.0))));
6199  points[2].x=(double) (center.x+cos(fmod((double) (alpha+(double) (i+1)*
6200  theta/arc_segments),DegreesToRadians(360.0))));
6201  points[2].y=(double) (center.y+sin(fmod((double) (alpha+(double) (i+1)*
6202  theta/arc_segments),DegreesToRadians(360.0))));
6203  points[1].x=(double) (points[2].x+gamma*sin(fmod((double) (alpha+(double)
6204  (i+1)*theta/arc_segments),DegreesToRadians(360.0))));
6205  points[1].y=(double) (points[2].y-gamma*cos(fmod((double) (alpha+(double)
6206  (i+1)*theta/arc_segments),DegreesToRadians(360.0))));
6207  p->point.x=(p == primitive_info) ? start.x : (p-1)->point.x;
6208  p->point.y=(p == primitive_info) ? start.y : (p-1)->point.y;
6209  (p+1)->point.x=(double) (cosine*radii.x*points[0].x-sine*radii.y*
6210  points[0].y);
6211  (p+1)->point.y=(double) (sine*radii.x*points[0].x+cosine*radii.y*
6212  points[0].y);
6213  (p+2)->point.x=(double) (cosine*radii.x*points[1].x-sine*radii.y*
6214  points[1].y);
6215  (p+2)->point.y=(double) (sine*radii.x*points[1].x+cosine*radii.y*
6216  points[1].y);
6217  (p+3)->point.x=(double) (cosine*radii.x*points[2].x-sine*radii.y*
6218  points[2].y);
6219  (p+3)->point.y=(double) (sine*radii.x*points[2].x+cosine*radii.y*
6220  points[2].y);
6221  if (i == (ssize_t) (arc_segments-1))
6222  (p+3)->point=end;
6223  status&=TraceBezier(mvg_info,4);
6224  if (status == 0)
6225  break;
6226  p=(*mvg_info->primitive_info)+mvg_info->offset;
6227  mvg_info->offset+=p->coordinates;
6228  p+=p->coordinates;
6229  }
6230  if (status == 0)
6231  return(MagickFalse);
6232  mvg_info->offset=offset;
6233  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6234  primitive_info->coordinates=(size_t) (p-primitive_info);
6235  primitive_info->closed_subpath=MagickFalse;
6236  for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
6237  {
6238  p->primitive=primitive_info->primitive;
6239  p--;
6240  }
6241  return(MagickTrue);
6242 }
6243 
6244 static MagickBooleanType TraceBezier(MVGInfo *mvg_info,
6245  const size_t number_coordinates)
6246 {
6247  double
6248  alpha,
6249  *coefficients,
6250  weight;
6251 
6252  PointInfo
6253  end,
6254  point,
6255  *points;
6256 
6258  *primitive_info;
6259 
6261  *p;
6262 
6263  ssize_t
6264  i,
6265  j;
6266 
6267  size_t
6268  control_points,
6269  quantum;
6270 
6271  /*
6272  Allocate coefficients.
6273  */
6274  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6275  quantum=number_coordinates;
6276  for (i=0; i < (ssize_t) number_coordinates; i++)
6277  {
6278  for (j=i+1; j < (ssize_t) number_coordinates; j++)
6279  {
6280  alpha=fabs(primitive_info[j].point.x-primitive_info[i].point.x);
6281  if (alpha > (double) MAGICK_SSIZE_MAX)
6282  {
6283  (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
6284  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
6285  return(MagickFalse);
6286  }
6287  if (alpha > (double) quantum)
6288  quantum=(size_t) alpha;
6289  alpha=fabs(primitive_info[j].point.y-primitive_info[i].point.y);
6290  if (alpha > (double) MAGICK_SSIZE_MAX)
6291  {
6292  (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
6293  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
6294  return(MagickFalse);
6295  }
6296  if (alpha > (double) quantum)
6297  quantum=(size_t) alpha;
6298  }
6299  }
6300  coefficients=(double *) AcquireQuantumMemory(number_coordinates,
6301  sizeof(*coefficients));
6302  quantum=MagickMin(quantum/number_coordinates,BezierQuantum);
6303  points=(PointInfo *) AcquireQuantumMemory(quantum,number_coordinates*
6304  sizeof(*points));
6305  if ((coefficients == (double *) NULL) || (points == (PointInfo *) NULL))
6306  {
6307  if (points != (PointInfo *) NULL)
6308  points=(PointInfo *) RelinquishMagickMemory(points);
6309  if (coefficients != (double *) NULL)
6310  coefficients=(double *) RelinquishMagickMemory(coefficients);
6311  (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
6312  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
6313  return(MagickFalse);
6314  }
6315  control_points=quantum*number_coordinates;
6316  if (CheckPrimitiveExtent(mvg_info,(double) control_points+1) == MagickFalse)
6317  {
6318  points=(PointInfo *) RelinquishMagickMemory(points);
6319  coefficients=(double *) RelinquishMagickMemory(coefficients);
6320  return(MagickFalse);
6321  }
6322  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6323  /*
6324  Compute bezier points.
6325  */
6326  end=primitive_info[number_coordinates-1].point;
6327  for (i=0; i < (ssize_t) number_coordinates; i++)
6328  coefficients[i]=Permutate((ssize_t) number_coordinates-1,i);
6329  weight=0.0;
6330  for (i=0; i < (ssize_t) control_points; i++)
6331  {
6332  p=primitive_info;
6333  point.x=0.0;
6334  point.y=0.0;
6335  alpha=pow((double) (1.0-weight),(double) number_coordinates-1.0);
6336  for (j=0; j < (ssize_t) number_coordinates; j++)
6337  {
6338  point.x+=alpha*coefficients[j]*p->point.x;
6339  point.y+=alpha*coefficients[j]*p->point.y;
6340  alpha*=weight/(1.0-weight);
6341  p++;
6342  }
6343  points[i]=point;
6344  weight+=1.0/control_points;
6345  }
6346  /*
6347  Bezier curves are just short segmented polys.
6348  */
6349  p=primitive_info;
6350  for (i=0; i < (ssize_t) control_points; i++)
6351  {
6352  if (TracePoint(p,points[i]) == MagickFalse)
6353  {
6354  points=(PointInfo *) RelinquishMagickMemory(points);
6355  coefficients=(double *) RelinquishMagickMemory(coefficients);
6356  return(MagickFalse);
6357  }
6358  p+=p->coordinates;
6359  }
6360  if (TracePoint(p,end) == MagickFalse)
6361  {
6362  points=(PointInfo *) RelinquishMagickMemory(points);
6363  coefficients=(double *) RelinquishMagickMemory(coefficients);
6364  return(MagickFalse);
6365  }
6366  p+=p->coordinates;
6367  primitive_info->coordinates=(size_t) (p-primitive_info);
6368  primitive_info->closed_subpath=MagickFalse;
6369  for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
6370  {
6371  p->primitive=primitive_info->primitive;
6372  p--;
6373  }
6374  points=(PointInfo *) RelinquishMagickMemory(points);
6375  coefficients=(double *) RelinquishMagickMemory(coefficients);
6376  return(MagickTrue);
6377 }
6378 
6379 static MagickBooleanType TraceCircle(MVGInfo *mvg_info,const PointInfo start,
6380  const PointInfo end)
6381 {
6382  double
6383  alpha,
6384  beta,
6385  radius;
6386 
6387  PointInfo
6388  offset,
6389  degrees;
6390 
6391  alpha=end.x-start.x;
6392  beta=end.y-start.y;
6393  radius=hypot((double) alpha,(double) beta);
6394  offset.x=(double) radius;
6395  offset.y=(double) radius;
6396  degrees.x=0.0;
6397  degrees.y=360.0;
6398  return(TraceEllipse(mvg_info,start,offset,degrees));
6399 }
6400 
6401 static MagickBooleanType TraceEllipse(MVGInfo *mvg_info,const PointInfo center,
6402  const PointInfo radii,const PointInfo arc)
6403 {
6404  double
6405  coordinates,
6406  delta,
6407  step,
6408  x,
6409  y;
6410 
6411  PointInfo
6412  angle,
6413  point;
6414 
6416  *primitive_info;
6417 
6419  *p;
6420 
6421  ssize_t
6422  i;
6423 
6424  /*
6425  Ellipses are just short segmented polys.
6426  */
6427  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6428  primitive_info->coordinates=0;
6429  if ((fabs(radii.x) < MagickEpsilon) || (fabs(radii.y) < MagickEpsilon))
6430  return(MagickTrue);
6431  delta=2.0*PerceptibleReciprocal(MagickMax(radii.x,radii.y));
6432  step=MagickPI/8.0;
6433  if ((delta >= 0.0) && (delta < (MagickPI/8.0)))
6434  step=MagickPI/4.0/(MagickPI*PerceptibleReciprocal(delta)/2.0);
6435  angle.x=DegreesToRadians(arc.x);
6436  y=arc.y;
6437  while (y < arc.x)
6438  y+=360.0;
6439  angle.y=DegreesToRadians(y);
6440  coordinates=ceil((angle.y-angle.x)/step+1.0);
6441  if (CheckPrimitiveExtent(mvg_info,coordinates) == MagickFalse)
6442  return(MagickFalse);
6443  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6444  for (p=primitive_info; angle.x < angle.y; angle.x+=step)
6445  {
6446  point.x=cos(fmod(angle.x,DegreesToRadians(360.0)))*radii.x+center.x;
6447  point.y=sin(fmod(angle.x,DegreesToRadians(360.0)))*radii.y+center.y;
6448  if (TracePoint(p,point) == MagickFalse)
6449  return(MagickFalse);
6450  p+=p->coordinates;
6451  }
6452  point.x=cos(fmod(angle.y,DegreesToRadians(360.0)))*radii.x+center.x;
6453  point.y=sin(fmod(angle.y,DegreesToRadians(360.0)))*radii.y+center.y;
6454  if (TracePoint(p,point) == MagickFalse)
6455  return(MagickFalse);
6456  p+=p->coordinates;
6457  primitive_info->coordinates=(size_t) (p-primitive_info);
6458  primitive_info->closed_subpath=MagickFalse;
6459  x=fabs(primitive_info[0].point.x-
6460  primitive_info[primitive_info->coordinates-1].point.x);
6461  y=fabs(primitive_info[0].point.y-
6462  primitive_info[primitive_info->coordinates-1].point.y);
6463  if ((x < MagickEpsilon) && (y < MagickEpsilon))
6464  primitive_info->closed_subpath=MagickTrue;
6465  for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
6466  {
6467  p->primitive=primitive_info->primitive;
6468  p--;
6469  }
6470  return(MagickTrue);
6471 }
6472 
6473 static MagickBooleanType TraceLine(PrimitiveInfo *primitive_info,
6474  const PointInfo start,const PointInfo end)
6475 {
6476  if (TracePoint(primitive_info,start) == MagickFalse)
6477  return(MagickFalse);
6478  if (TracePoint(primitive_info+1,end) == MagickFalse)
6479  return(MagickFalse);
6480  (primitive_info+1)->primitive=primitive_info->primitive;
6481  primitive_info->coordinates=2;
6482  primitive_info->closed_subpath=MagickFalse;
6483  return(MagickTrue);
6484 }
6485 
6486 static ssize_t TracePath(Image *image,MVGInfo *mvg_info,const char *path)
6487 {
6488  char
6489  *next_token,
6490  token[MaxTextExtent] = "";
6491 
6492  const char
6493  *p;
6494 
6495  double
6496  x,
6497  y;
6498 
6499  int
6500  attribute,
6501  last_attribute;
6502 
6503  MagickStatusType
6504  status;
6505 
6506  PointInfo
6507  end = {0.0, 0.0},
6508  points[4] = { {0.0, 0.0}, {0.0, 0.0}, {0.0, 0.0}, {0.0, 0.0} },
6509  point = {0.0, 0.0},
6510  start = {0.0, 0.0};
6511 
6513  *primitive_info;
6514 
6515  PrimitiveType
6516  primitive_type;
6517 
6519  *q;
6520 
6521  ssize_t
6522  i;
6523 
6524  size_t
6525  number_coordinates,
6526  z_count;
6527 
6528  ssize_t
6529  subpath_offset;
6530 
6531  subpath_offset=mvg_info->offset;
6532  primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
6533  status=MagickTrue;
6534  attribute=0;
6535  number_coordinates=0;
6536  z_count=0;
6537  *token='\0';
6538  primitive_type=primitive_info->primitive;
6539  q=primitive_info;
6540  for (p=path; *p != '\0'; )
6541  {
6542  if (status == MagickFalse)
6543  break;
6544  while (isspace((int) ((unsigned char) *p)) != 0)
6545  p++;
6546  if (*p == '\0')
6547  break;
6548  last_attribute=attribute;
6549  attribute=(int) (*p++);
6550  switch (attribute)
6551  {
6552  case 'a':
6553  case 'A':
6554  {
6555  double
6556  angle = 0.0;
6557 
6558  MagickBooleanType
6559  large_arc = MagickFalse,
6560  sweep = MagickFalse;
6561 
6562  PointInfo
6563  arc = {0.0, 0.0};
6564 
6565  /*
6566  Elliptical arc.
6567  */
6568  do
6569  {
6570  (void) GetNextToken(p,&p,MaxTextExtent,token);
6571  if (*token == ',')
6572  (void) GetNextToken(p,&p,MaxTextExtent,token);
6573  arc.x=GetDrawValue(token,&next_token);
6574  if (token == next_token)
6575  ThrowPointExpectedException(image,token);
6576  (void) GetNextToken(p,&p,MaxTextExtent,token);
6577  if (*token == ',')
6578  (void) GetNextToken(p,&p,MaxTextExtent,token);
6579  arc.y=GetDrawValue(token,&next_token);
6580  if (token == next_token)
6581  ThrowPointExpectedException(image,token);
6582  (void) GetNextToken(p,&p,MaxTextExtent,token);
6583  if (*token == ',')
6584  (void) GetNextToken(p,&p,MaxTextExtent,token);
6585  angle=GetDrawValue(token,&next_token);
6586  if (token == next_token)
6587  ThrowPointExpectedException(image,token);
6588  (void) GetNextToken(p,&p,MaxTextExtent,token);
6589  if (*token == ',')
6590  (void) GetNextToken(p,&p,MaxTextExtent,token);
6591  large_arc=StringToLong(token) != 0 ? MagickTrue : MagickFalse;
6592  (void) GetNextToken(p,&p,MaxTextExtent,token);
6593  if (*token == ',')
6594  (void) GetNextToken(p,&p,MaxTextExtent,token);
6595  sweep=StringToLong(token) != 0 ? MagickTrue : MagickFalse;
6596  if (*token == ',')
6597  (void) GetNextToken(p,&p,MaxTextExtent,token);
6598  (void) GetNextToken(p,&p,MaxTextExtent,token);
6599  if (*token == ',')
6600  (void) GetNextToken(p,&p,MaxTextExtent,token);
6601  x=GetDrawValue(token,&next_token);
6602  if (token == next_token)
6603  ThrowPointExpectedException(image,token);
6604  (void) GetNextToken(p,&p,MaxTextExtent,token);
6605  if (*token == ',')
6606  (void) GetNextToken(p,&p,MaxTextExtent,token);
6607  y=GetDrawValue(token,&next_token);
6608  if (token == next_token)
6609  ThrowPointExpectedException(image,token);
6610  end.x=(double) (attribute == (int) 'A' ? x : point.x+x);
6611  end.y=(double) (attribute == (int) 'A' ? y : point.y+y);
6612  status&=TraceArcPath(mvg_info,point,end,arc,angle,large_arc,sweep);
6613  q=(*mvg_info->primitive_info)+mvg_info->offset;
6614  mvg_info->offset+=q->coordinates;
6615  q+=q->coordinates;
6616  point=end;
6617  while (isspace((int) ((unsigned char) *p)) != 0)
6618  p++;
6619  if (*p == ',')
6620  p++;
6621  } while (IsPoint(p) != MagickFalse);
6622  break;
6623  }
6624  case 'c':
6625  case 'C':
6626  {
6627  /*
6628  Cubic Bézier curve.
6629  */
6630  do
6631  {
6632  points[0]=point;
6633  for (i=1; i < 4; i++)
6634  {
6635  (void) GetNextToken(p,&p,MaxTextExtent,token);
6636  if (*token == ',')
6637  (void) GetNextToken(p,&p,MaxTextExtent,token);
6638  x=GetDrawValue(token,&next_token);
6639  if (token == next_token)
6640  ThrowPointExpectedException(image,token);
6641  (void) GetNextToken(p,&p,MaxTextExtent,token);
6642  if (*token == ',')
6643  (void) GetNextToken(p,&p,MaxTextExtent,token);
6644  y=GetDrawValue(token,&next_token);
6645  if (token == next_token)
6646  ThrowPointExpectedException(image,token);
6647  end.x=(double) (attribute == (int) 'C' ? x : point.x+x);
6648  end.y=(double) (attribute == (int) 'C' ? y : point.y+y);
6649  points[i]=end;
6650  }
6651  for (i=0; i < 4; i++)
6652  (q+i)->point=points[i];
6653  if (TraceBezier(mvg_info,4) == MagickFalse)
6654  return(-1);
6655  q=(*mvg_info->primitive_info)+mvg_info->offset;
6656  mvg_info->offset+=q->coordinates;
6657  q+=q->coordinates;
6658  point=end;
6659  while (isspace((int) ((unsigned char) *p)) != 0)
6660  p++;
6661  if (*p == ',')
6662  p++;
6663  } while (IsPoint(p) != MagickFalse);
6664  break;
6665  }
6666  case 'H':
6667  case 'h':
6668  {
6669  do
6670  {
6671  (void) GetNextToken(p,&p,MaxTextExtent,token);
6672  if (*token == ',')
6673  (void) GetNextToken(p,&p,MaxTextExtent,token);
6674  x=GetDrawValue(token,&next_token);
6675  if (token == next_token)
6676  ThrowPointExpectedException(image,token);
6677  point.x=(double) (attribute == (int) 'H' ? x: point.x+x);
6678  if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
6679  return(-1);
6680  q=(*mvg_info->primitive_info)+mvg_info->offset;
6681  if (TracePoint(q,point) == MagickFalse)
6682  return(-1);
6683  mvg_info->offset+=q->coordinates;
6684  q+=q->coordinates;
6685  while (isspace((int) ((unsigned char) *p)) != 0)
6686  p++;
6687  if (*p == ',')
6688  p++;
6689  } while (IsPoint(p) != MagickFalse);
6690  break;
6691  }
6692  case 'l':
6693  case 'L':
6694  {
6695  /*
6696  Line to.
6697  */
6698  do
6699  {
6700  (void) GetNextToken(p,&p,MaxTextExtent,token);
6701  if (*token == ',')
6702  (void) GetNextToken(p,&p,MaxTextExtent,token);
6703  x=GetDrawValue(token,&next_token);
6704  if (token == next_token)
6705  ThrowPointExpectedException(image,token);
6706  (void) GetNextToken(p,&p,MaxTextExtent,token);
6707  if (*token == ',')
6708  (void) GetNextToken(p,&p,MaxTextExtent,token);
6709  y=GetDrawValue(token,&next_token);
6710  if (token == next_token)
6711  ThrowPointExpectedException(image,token);
6712  point.x=(double) (attribute == (int) 'L' ? x : point.x+x);
6713  point.y=(double) (attribute == (int) 'L' ? y : point.y+y);
6714  if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
6715  return(-1);
6716  q=(*mvg_info->primitive_info)+mvg_info->offset;
6717  if (TracePoint(q,point) == MagickFalse)
6718  return(-1);
6719  mvg_info->offset+=q->coordinates;
6720  q+=q->coordinates;
6721  while (isspace((int) ((unsigned char) *p)) != 0)
6722  p++;
6723  if (*p == ',')
6724  p++;
6725  } while (IsPoint(p) != MagickFalse);
6726  break;
6727  }
6728  case 'M':
6729  case 'm':
6730  {
6731  /*
6732  Move to.
6733  */
6734  if (mvg_info->offset != subpath_offset)
6735  {
6736  primitive_info=(*mvg_info->primitive_info)+subpath_offset;
6737  primitive_info->coordinates=(size_t) (q-primitive_info);
6738  number_coordinates+=primitive_info->coordinates;
6739  primitive_info=q;
6740  subpath_offset=mvg_info->offset;
6741  }
6742  i=0;
6743  do
6744  {
6745  (void) GetNextToken(p,&p,MaxTextExtent,token);
6746  if (*token == ',')
6747  (void) GetNextToken(p,&p,MaxTextExtent,token);
6748  x=GetDrawValue(token,&next_token);
6749  if (token == next_token)
6750  ThrowPointExpectedException(image,token);
6751  (void) GetNextToken(p,&p,MaxTextExtent,token);
6752  if (*token == ',')
6753  (void) GetNextToken(p,&p,MaxTextExtent,token);
6754  y=GetDrawValue(token,&next_token);
6755  if (token == next_token)
6756  ThrowPointExpectedException(image,token);
6757  point.x=(double) (attribute == (int) 'M' ? x : point.x+x);
6758  point.y=(double) (attribute == (int) 'M' ? y : point.y+y);
6759  if (i == 0)
6760  start=point;
6761  i++;
6762  if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
6763  return(-1);
6764  q=(*mvg_info->primitive_info)+mvg_info->offset;
6765  if (TracePoint(q,point) == MagickFalse)
6766  return(-1);
6767  mvg_info->offset+=q->coordinates;
6768  q+=q->coordinates;
6769  while (isspace((int) ((unsigned char) *p)) != 0)
6770  p++;
6771  if (*p == ',')
6772  p++;
6773  } while (IsPoint(p) != MagickFalse);
6774  break;
6775  }
6776  case 'q':
6777  case 'Q':
6778  {
6779  /*
6780  Quadratic Bézier curve.
6781  */
6782  do
6783  {
6784  points[0]=point;
6785  for (i=1; i < 3; i++)
6786  {
6787  (void) GetNextToken(p,&p,MaxTextExtent,token);
6788  if (*token == ',')
6789  (void) GetNextToken(p,&p,MaxTextExtent,token);
6790  x=GetDrawValue(token,&next_token);
6791  if (token == next_token)
6792  ThrowPointExpectedException(image,token);
6793  (void) GetNextToken(p,&p,MaxTextExtent,token);
6794  if (*token == ',')
6795  (void) GetNextToken(p,&p,MaxTextExtent,token);
6796  y=GetDrawValue(token,&next_token);
6797  if (token == next_token)
6798  ThrowPointExpectedException(image,token);
6799  if (*p == ',')
6800  p++;
6801  end.x=(double) (attribute == (int) 'Q' ? x : point.x+x);
6802  end.y=(double) (attribute == (int) 'Q' ? y : point.y+y);
6803  points[i]=end;
6804  }
6805  for (i=0; i < 3; i++)
6806  (q+i)->point=points[i];
6807  if (TraceBezier(mvg_info,3) == MagickFalse)
6808  return(-1);
6809  q=(*mvg_info->primitive_info)+mvg_info->offset;
6810  mvg_info->offset+=q->coordinates;
6811  q+=q->coordinates;
6812  point=end;
6813  while (isspace((int) ((unsigned char) *p)) != 0)
6814  p++;
6815  if (*p == ',')
6816  p++;
6817  } while (IsPoint(p) != MagickFalse);
6818  break;
6819  }
6820  case 's':
6821  case 'S':
6822  {
6823  /*
6824  Cubic Bézier curve.
6825  */
6826  do
6827  {
6828  points[0]=points[3];
6829  points[1].x=2.0*points[3].x-points[2].x;
6830  points[1].y=2.0*points[3].y-points[2].y;
6831  for (i=2; i < 4; i++)
6832  {
6833  (void) GetNextToken(p,&p,MaxTextExtent,token);
6834  if (*token == ',')
6835  (void) GetNextToken(p,&p,MaxTextExtent,token);
6836  x=GetDrawValue(token,&next_token);
6837  if (token == next_token)
6838  ThrowPointExpectedException(image,token);
6839  (void) GetNextToken(p,&p,MaxTextExtent,token);
6840  if (*token == ',')
6841  (void) GetNextToken(p,&p,MaxTextExtent,token);
6842  y=GetDrawValue(token,&next_token);
6843  if (token == next_token)
6844  ThrowPointExpectedException(image,token);
6845  if (*p == ',')
6846  p++;
6847  end.x=(double) (attribute == (int) 'S' ? x : point.x+x);
6848  end.y=(double) (attribute == (int) 'S' ? y : point.y+y);
6849  points[i]=end;
6850  }
6851  if (strchr("CcSs",last_attribute) == (char *) NULL)
6852  {
6853  points[0]=point;
6854  points[1]=point;
6855  }
6856  for (i=0; i < 4; i++)
6857  (q+i)->point=points[i];
6858  if (TraceBezier(mvg_info,4) == MagickFalse)
6859  return(-1);
6860  q=(*mvg_info->primitive_info)+mvg_info->offset;
6861  mvg_info->offset+=q->coordinates;
6862  q+=q->coordinates;
6863  point=end;
6864  last_attribute=attribute;
6865  while (isspace((int) ((unsigned char) *p)) != 0)
6866  p++;
6867  if (*p == ',')
6868  p++;
6869  } while (IsPoint(p) != MagickFalse);
6870  break;
6871  }
6872  case 't':
6873  case 'T':
6874  {
6875  /*
6876  Quadratic Bézier curve.
6877  */
6878  do
6879  {
6880  points[0]=points[2];
6881  points[1].x=2.0*points[2].x-points[1].x;
6882  points[1].y=2.0*points[2].y-points[1].y;
6883  for (i=2; i < 3; i++)
6884  {
6885  (void) GetNextToken(p,&p,MaxTextExtent,token);
6886  if (*token == ',')
6887  (void) GetNextToken(p,&p,MaxTextExtent,token);
6888  x=GetDrawValue(token,&next_token);
6889  if (token == next_token)
6890  ThrowPointExpectedException(image,token);
6891  (void) GetNextToken(p,&p,MaxTextExtent,token);
6892  if (*token == ',')
6893  (void) GetNextToken(p,&p,MaxTextExtent,token);
6894  y=GetDrawValue(token,&next_token);
6895  if (token == next_token)
6896  ThrowPointExpectedException(image,token);
6897  end.x=(double) (attribute == (int) 'T' ? x : point.x+x);
6898  end.y=(double) (attribute == (int) 'T' ? y : point.y+y);
6899  points[i]=end;
6900  }
6901  if (status == MagickFalse)
6902  break;
6903  if (strchr("QqTt",last_attribute) == (char *) NULL)
6904  {
6905  points[0]=point;
6906  points[1]=point;
6907  }
6908  for (i=0; i < 3; i++)
6909  (q+i)->point=points[i];
6910  if (TraceBezier(mvg_info,3) == MagickFalse)
6911  return(-1);
6912  q=(*mvg_info->primitive_info)+mvg_info->offset;
6913  mvg_info->offset+=q->coordinates;
6914  q+=q->coordinates;
6915  point=end;
6916  last_attribute=attribute;
6917  while (isspace((int) ((unsigned char) *p)) != 0)
6918  p++;
6919  if (*p == ',')
6920  p++;
6921  } while (IsPoint(p) != MagickFalse);
6922  break;
6923  }
6924  case 'v':
6925  case 'V':
6926  {
6927  /*
6928  Line to.
6929  */
6930  do
6931  {
6932  (void) GetNextToken(p,&p,MaxTextExtent,token);
6933  if (*token == ',')
6934  (void) GetNextToken(p,&p,MaxTextExtent,token);
6935  y=GetDrawValue(token,&next_token);
6936  if (token == next_token)
6937  ThrowPointExpectedException(image,token);
6938  point.y=(double) (attribute == (int) 'V' ? y : point.y+y);
6939  if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
6940  return(-1);
6941  q=(*mvg_info->primitive_info)+mvg_info->offset;
6942  if (TracePoint(q,point) == MagickFalse)
6943  return(-1);
6944  mvg_info->offset+=q->coordinates;
6945  q+=q->coordinates;
6946  while (isspace((int) ((unsigned char) *p)) != 0)
6947  p++;
6948  if (*p == ',')
6949  p++;
6950  } while (IsPoint(p) != MagickFalse);
6951  break;
6952  }
6953  case 'z':
6954  case 'Z':
6955  {
6956  /*
6957  Close path.
6958  */
6959  point=start;
6960  if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
6961  return(-1);
6962  q=(*mvg_info->primitive_info)+mvg_info->offset;
6963  if (TracePoint(q,point) == MagickFalse)
6964  return(-1);
6965  mvg_info->offset+=q->coordinates;
6966  q+=q->coordinates;
6967  primitive_info=(*mvg_info->primitive_info)+subpath_offset;
6968  primitive_info->coordinates=(size_t) (q-primitive_info);
6969  primitive_info->closed_subpath=MagickTrue;
6970  number_coordinates+=primitive_info->coordinates;
6971  primitive_info=q;
6972  subpath_offset=mvg_info->offset;
6973  z_count++;
6974  break;
6975  }
6976  default:
6977  {
6978  ThrowPointExpectedException(image,token);
6979  break;
6980  }
6981  }
6982  }
6983  if (status == MagickFalse)
6984  return(-1);
6985  primitive_info=(*mvg_info->primitive_info)+subpath_offset;
6986  primitive_info->coordinates=(size_t) (q-primitive_info);
6987  number_coordinates+=primitive_info->coordinates;
6988  for (i=0; i < (ssize_t) number_coordinates; i++)
6989  {
6990  q--;
6991  q->primitive=primitive_type;
6992  if (z_count > 1)
6993  q->method=FillToBorderMethod;
6994  }
6995  q=primitive_info;
6996  return((ssize_t) number_coordinates);
6997 }
6998 
6999 static MagickBooleanType TraceRectangle(PrimitiveInfo *primitive_info,
7000  const PointInfo start,const PointInfo end)
7001 {
7002  PointInfo
7003  point;
7004 
7006  *p;
7007 
7008  ssize_t
7009  i;
7010 
7011  p=primitive_info;
7012  if (TracePoint(p,start) == MagickFalse)
7013  return(MagickFalse);
7014  p+=p->coordinates;
7015  point.x=start.x;
7016  point.y=end.y;
7017  if (TracePoint(p,point) == MagickFalse)
7018  return(MagickFalse);
7019  p+=p->coordinates;
7020  if (TracePoint(p,end) == MagickFalse)
7021  return(MagickFalse);
7022  p+=p->coordinates;
7023  point.x=end.x;
7024  point.y=start.y;
7025  if (TracePoint(p,point) == MagickFalse)
7026  return(MagickFalse);
7027  p+=p->coordinates;
7028  if (TracePoint(p,start) == MagickFalse)
7029  return(MagickFalse);
7030  p+=p->coordinates;
7031  primitive_info->coordinates=(size_t) (p-primitive_info);
7032  primitive_info->closed_subpath=MagickTrue;
7033  for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
7034  {
7035  p->primitive=primitive_info->primitive;
7036  p--;
7037  }
7038  return(MagickTrue);
7039 }
7040 
7041 static MagickBooleanType TraceRoundRectangle(MVGInfo *mvg_info,
7042  const PointInfo start,const PointInfo end,PointInfo arc)
7043 {
7044  PointInfo
7045  degrees,
7046  point,
7047  segment;
7048 
7050  *primitive_info;
7051 
7053  *p;
7054 
7055  ssize_t
7056  i;
7057 
7058  ssize_t
7059  offset;
7060 
7061  offset=mvg_info->offset;
7062  segment.x=fabs(end.x-start.x);
7063  segment.y=fabs(end.y-start.y);
7064  if ((segment.x < MagickEpsilon) || (segment.y < MagickEpsilon))
7065  {
7066  (*mvg_info->primitive_info+mvg_info->offset)->coordinates=0;
7067  return(MagickTrue);
7068  }
7069  if (arc.x > (0.5*segment.x))
7070  arc.x=0.5*segment.x;
7071  if (arc.y > (0.5*segment.y))
7072  arc.y=0.5*segment.y;
7073  point.x=start.x+segment.x-arc.x;
7074  point.y=start.y+arc.y;
7075  degrees.x=270.0;
7076  degrees.y=360.0;
7077  if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
7078  return(MagickFalse);
7079  p=(*mvg_info->primitive_info)+mvg_info->offset;
7080  mvg_info->offset+=p->coordinates;
7081  point.x=start.x+segment.x-arc.x;
7082  point.y=start.y+segment.y-arc.y;
7083  degrees.x=0.0;
7084  degrees.y=90.0;
7085  if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
7086  return(MagickFalse);
7087  p=(*mvg_info->primitive_info)+mvg_info->offset;
7088  mvg_info->offset+=p->coordinates;
7089  point.x=start.x+arc.x;
7090  point.y=start.y+segment.y-arc.y;
7091  degrees.x=90.0;
7092  degrees.y=180.0;
7093  if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
7094  return(MagickFalse);
7095  p=(*mvg_info->primitive_info)+mvg_info->offset;
7096  mvg_info->offset+=p->coordinates;
7097  point.x=start.x+arc.x;
7098  point.y=start.y+arc.y;
7099  degrees.x=180.0;
7100  degrees.y=270.0;
7101  if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
7102  return(MagickFalse);
7103  p=(*mvg_info->primitive_info)+mvg_info->offset;
7104  mvg_info->offset+=p->coordinates;
7105  if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
7106  return(MagickFalse);
7107  p=(*mvg_info->primitive_info)+mvg_info->offset;
7108  if (TracePoint(p,(*mvg_info->primitive_info+offset)->point) == MagickFalse)
7109  return(MagickFalse);
7110  p+=p->coordinates;
7111  mvg_info->offset=offset;
7112  primitive_info=(*mvg_info->primitive_info)+offset;
7113  primitive_info->coordinates=(size_t) (p-primitive_info);
7114  primitive_info->closed_subpath=MagickTrue;
7115  for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
7116  {
7117  p->primitive=primitive_info->primitive;
7118  p--;
7119  }
7120  return(MagickTrue);
7121 }
7122 
7123 static MagickBooleanType TraceSquareLinecap(PrimitiveInfo *primitive_info,
7124  const size_t number_vertices,const double offset)
7125 {
7126  double
7127  distance;
7128 
7129  double
7130  dx,
7131  dy;
7132 
7133  ssize_t
7134  i;
7135 
7136  ssize_t
7137  j;
7138 
7139  dx=0.0;
7140  dy=0.0;
7141  for (i=1; i < (ssize_t) number_vertices; i++)
7142  {
7143  dx=primitive_info[0].point.x-primitive_info[i].point.x;
7144  dy=primitive_info[0].point.y-primitive_info[i].point.y;
7145  if ((fabs((double) dx) >= MagickEpsilon) ||
7146  (fabs((double) dy) >= MagickEpsilon))
7147  break;
7148  }
7149  if (i == (ssize_t) number_vertices)
7150  i=(ssize_t) number_vertices-1L;
7151  distance=hypot((double) dx,(double) dy);
7152  primitive_info[0].point.x=(double) (primitive_info[i].point.x+
7153  dx*(distance+offset)/distance);
7154  primitive_info[0].point.y=(double) (primitive_info[i].point.y+
7155  dy*(distance+offset)/distance);
7156  for (j=(ssize_t) number_vertices-2; j >= 0; j--)
7157  {
7158  dx=primitive_info[number_vertices-1].point.x-primitive_info[j].point.x;
7159  dy=primitive_info[number_vertices-1].point.y-primitive_info[j].point.y;
7160  if ((fabs((double) dx) >= MagickEpsilon) ||
7161  (fabs((double) dy) >= MagickEpsilon))
7162  break;
7163  }
7164  distance=hypot((double) dx,(double) dy);
7165  primitive_info[number_vertices-1].point.x=(double) (primitive_info[j].point.x+
7166  dx*(distance+offset)/distance);
7167  primitive_info[number_vertices-1].point.y=(double) (primitive_info[j].point.y+
7168  dy*(distance+offset)/distance);
7169  return(MagickTrue);
7170 }
7171 
7172 static PrimitiveInfo *TraceStrokePolygon(const DrawInfo *draw_info,
7173  const PrimitiveInfo *primitive_info,ExceptionInfo *exception)
7174 {
7175 #define MaxStrokePad (6*BezierQuantum+360)
7176 #define CheckPathExtent(pad_p,pad_q) \
7177 { \
7178  if ((pad_p) > MaxBezierCoordinates) \
7179  stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p); \
7180  else \
7181  if ((ssize_t) (p+(pad_p)) >= (ssize_t) extent_p) \
7182  { \
7183  if (~extent_p < (pad_p)) \
7184  stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p); \
7185  else \
7186  { \
7187  extent_p+=(pad_p); \
7188  stroke_p=(PointInfo *) ResizeQuantumMemory(stroke_p,extent_p+ \
7189  MaxStrokePad,sizeof(*stroke_p)); \
7190  } \
7191  } \
7192  if ((pad_q) > MaxBezierCoordinates) \
7193  stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q); \
7194  else \
7195  if ((ssize_t) (q+(pad_q)) >= (ssize_t) extent_q) \
7196  { \
7197  if (~extent_q < (pad_q)) \
7198  stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q); \
7199  else \
7200  { \
7201  extent_q+=(pad_q); \
7202  stroke_q=(PointInfo *) ResizeQuantumMemory(stroke_q,extent_q+ \
7203  MaxStrokePad,sizeof(*stroke_q)); \
7204  } \
7205  } \
7206  if ((stroke_p == (PointInfo *) NULL) || (stroke_q == (PointInfo *) NULL)) \
7207  { \
7208  if (stroke_p != (PointInfo *) NULL) \
7209  stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p); \
7210  if (stroke_q != (PointInfo *) NULL) \
7211  stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q); \
7212  polygon_primitive=(PrimitiveInfo *) \
7213  RelinquishMagickMemory(polygon_primitive); \
7214  (void) ThrowMagickException(exception,GetMagickModule(), \
7215  ResourceLimitError,"MemoryAllocationFailed","`%s'",""); \
7216  return((PrimitiveInfo *) NULL); \
7217  } \
7218 }
7219 
7220  typedef struct _StrokeSegment
7221  {
7222  double
7223  p,
7224  q;
7225  } StrokeSegment;
7226 
7227  double
7228  delta_theta,
7229  dot_product,
7230  mid,
7231  miterlimit;
7232 
7233  MagickBooleanType
7234  closed_path;
7235 
7236  PointInfo
7237  box_p[5],
7238  box_q[5],
7239  center,
7240  offset,
7241  *stroke_p,
7242  *stroke_q;
7243 
7245  *polygon_primitive,
7246  *stroke_polygon;
7247 
7248  ssize_t
7249  i;
7250 
7251  size_t
7252  arc_segments,
7253  extent_p,
7254  extent_q,
7255  number_vertices;
7256 
7257  ssize_t
7258  j,
7259  n,
7260  p,
7261  q;
7262 
7263  StrokeSegment
7264  dx = {0.0, 0.0},
7265  dy = {0.0, 0.0},
7266  inverse_slope = {0.0, 0.0},
7267  slope = {0.0, 0.0},
7268  theta = {0.0, 0.0};
7269 
7270  /*
7271  Allocate paths.
7272  */
7273  number_vertices=primitive_info->coordinates;
7274  polygon_primitive=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
7275  number_vertices+2UL,sizeof(*polygon_primitive));
7276  if (polygon_primitive == (PrimitiveInfo *) NULL)
7277  {
7278  (void) ThrowMagickException(exception,GetMagickModule(),
7279  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
7280  return((PrimitiveInfo *) NULL);
7281  }
7282  (void) memcpy(polygon_primitive,primitive_info,(size_t) number_vertices*
7283  sizeof(*polygon_primitive));
7284  offset.x=primitive_info[number_vertices-1].point.x-primitive_info[0].point.x;
7285  offset.y=primitive_info[number_vertices-1].point.y-primitive_info[0].point.y;
7286  closed_path=(fabs(offset.x) < MagickEpsilon) &&
7287  (fabs(offset.y) < MagickEpsilon) ? MagickTrue : MagickFalse;
7288  if (((draw_info->linejoin == RoundJoin) ||
7289  (draw_info->linejoin == MiterJoin)) && (closed_path != MagickFalse))
7290  {
7291  polygon_primitive[number_vertices]=primitive_info[1];
7292  number_vertices++;
7293  }
7294  polygon_primitive[number_vertices].primitive=UndefinedPrimitive;
7295  /*
7296  Compute the slope for the first line segment, p.
7297  */
7298  dx.p=0.0;
7299  dy.p=0.0;
7300  for (n=1; n < (ssize_t) number_vertices; n++)
7301  {
7302  dx.p=polygon_primitive[n].point.x-polygon_primitive[0].point.x;
7303  dy.p=polygon_primitive[n].point.y-polygon_primitive[0].point.y;
7304  if ((fabs(dx.p) >= MagickEpsilon) || (fabs(dy.p) >= MagickEpsilon))
7305  break;
7306  }
7307  if (n == (ssize_t) number_vertices)
7308  {
7309  if ((draw_info->linecap != RoundCap) || (closed_path != MagickFalse))
7310  {
7311  /*
7312  Zero length subpath.
7313  */
7314  stroke_polygon=(PrimitiveInfo *) AcquireCriticalMemory(
7315  sizeof(*stroke_polygon));
7316  stroke_polygon[0]=polygon_primitive[0];
7317  stroke_polygon[0].coordinates=0;
7318  polygon_primitive=(PrimitiveInfo *) RelinquishMagickMemory(
7319  polygon_primitive);
7320  return(stroke_polygon);
7321  }
7322  n=(ssize_t) number_vertices-1L;
7323  }
7324  extent_p=2*number_vertices;
7325  extent_q=2*number_vertices;
7326  stroke_p=(PointInfo *) AcquireQuantumMemory((size_t) extent_p+MaxStrokePad,
7327  sizeof(*stroke_p));
7328  stroke_q=(PointInfo *) AcquireQuantumMemory((size_t) extent_q+MaxStrokePad,
7329  sizeof(*stroke_q));
7330  if ((stroke_p == (PointInfo *) NULL) || (stroke_q == (PointInfo *) NULL))
7331  {
7332  if (stroke_p != (PointInfo *) NULL)
7333  stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p);
7334  if (stroke_q != (PointInfo *) NULL)
7335  stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q);
7336  polygon_primitive=(PrimitiveInfo *)
7337  RelinquishMagickMemory(polygon_primitive);
7338  (void) ThrowMagickException(exception,GetMagickModule(),
7339  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
7340  return((PrimitiveInfo *) NULL);
7341  }
7342  slope.p=0.0;
7343  inverse_slope.p=0.0;
7344  if (fabs(dx.p) < MagickEpsilon)
7345  {
7346  if (dx.p >= 0.0)
7347  slope.p=dy.p < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
7348  else
7349  slope.p=dy.p < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
7350  }
7351  else
7352  if (fabs(dy.p) < MagickEpsilon)
7353  {
7354  if (dy.p >= 0.0)
7355  inverse_slope.p=dx.p < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
7356  else
7357  inverse_slope.p=dx.p < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
7358  }
7359  else
7360  {
7361  slope.p=dy.p/dx.p;
7362  inverse_slope.p=(-1.0*PerceptibleReciprocal(slope.p));
7363  }
7364  mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
7365  miterlimit=(double) (draw_info->miterlimit*draw_info->miterlimit*mid*mid);
7366  if ((draw_info->linecap == SquareCap) && (closed_path == MagickFalse))
7367  (void) TraceSquareLinecap(polygon_primitive,number_vertices,mid);
7368  offset.x=sqrt((double) (mid*mid/(inverse_slope.p*inverse_slope.p+1.0)));
7369  offset.y=(double) (offset.x*inverse_slope.p);
7370  if ((dy.p*offset.x-dx.p*offset.y) > 0.0)
7371  {
7372  box_p[0].x=polygon_primitive[0].point.x-offset.x;
7373  box_p[0].y=polygon_primitive[0].point.y-offset.x*inverse_slope.p;
7374  box_p[1].x=polygon_primitive[n].point.x-offset.x;
7375  box_p[1].y=polygon_primitive[n].point.y-offset.x*inverse_slope.p;
7376  box_q[0].x=polygon_primitive[0].point.x+offset.x;
7377  box_q[0].y=polygon_primitive[0].point.y+offset.x*inverse_slope.p;
7378  box_q[1].x=polygon_primitive[n].point.x+offset.x;
7379  box_q[1].y=polygon_primitive[n].point.y+offset.x*inverse_slope.p;
7380  }
7381  else
7382  {
7383  box_p[0].x=polygon_primitive[0].point.x+offset.x;
7384  box_p[0].y=polygon_primitive[0].point.y+offset.y;
7385  box_p[1].x=polygon_primitive[n].point.x+offset.x;
7386  box_p[1].y=polygon_primitive[n].point.y+offset.y;
7387  box_q[0].x=polygon_primitive[0].point.x-offset.x;
7388  box_q[0].y=polygon_primitive[0].point.y-offset.y;
7389  box_q[1].x=polygon_primitive[n].point.x-offset.x;
7390  box_q[1].y=polygon_primitive[n].point.y-offset.y;
7391  }
7392  /*
7393  Create strokes for the line join attribute: bevel, miter, round.
7394  */
7395  p=0;
7396  q=0;
7397  stroke_q[p++]=box_q[0];
7398  stroke_p[q++]=box_p[0];
7399  for (i=(ssize_t) n+1; i < (ssize_t) number_vertices; i++)
7400  {
7401  /*
7402  Compute the slope for this line segment, q.
7403  */
7404  dx.q=polygon_primitive[i].point.x-polygon_primitive[n].point.x;
7405  dy.q=polygon_primitive[i].point.y-polygon_primitive[n].point.y;
7406  dot_product=dx.q*dx.q+dy.q*dy.q;
7407  if (dot_product < 0.25)
7408  continue;
7409  slope.q=0.0;
7410  inverse_slope.q=0.0;
7411  if (fabs(dx.q) < MagickEpsilon)
7412  {
7413  if (dx.q >= 0.0)
7414  slope.q=dy.q < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
7415  else
7416  slope.q=dy.q < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
7417  }
7418  else
7419  if (fabs(dy.q) < MagickEpsilon)
7420  {
7421  if (dy.q >= 0.0)
7422  inverse_slope.q=dx.q < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
7423  else
7424  inverse_slope.q=dx.q < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
7425  }
7426  else
7427  {
7428  slope.q=dy.q/dx.q;
7429  inverse_slope.q=(-1.0*PerceptibleReciprocal(slope.q));
7430  }
7431  offset.x=sqrt((double) (mid*mid/(inverse_slope.q*inverse_slope.q+1.0)));
7432  offset.y=(double) (offset.x*inverse_slope.q);
7433  dot_product=dy.q*offset.x-dx.q*offset.y;
7434  if (dot_product > 0.0)
7435  {
7436  box_p[2].x=polygon_primitive[n].point.x-offset.x;
7437  box_p[2].y=polygon_primitive[n].point.y-offset.y;
7438  box_p[3].x=polygon_primitive[i].point.x-offset.x;
7439  box_p[3].y=polygon_primitive[i].point.y-offset.y;
7440  box_q[2].x=polygon_primitive[n].point.x+offset.x;
7441  box_q[2].y=polygon_primitive[n].point.y+offset.y;
7442  box_q[3].x=polygon_primitive[i].point.x+offset.x;
7443  box_q[3].y=polygon_primitive[i].point.y+offset.y;
7444  }
7445  else
7446  {
7447  box_p[2].x=polygon_primitive[n].point.x+offset.x;
7448  box_p[2].y=polygon_primitive[n].point.y+offset.y;
7449  box_p[3].x=polygon_primitive[i].point.x+offset.x;
7450  box_p[3].y=polygon_primitive[i].point.y+offset.y;
7451  box_q[2].x=polygon_primitive[n].point.x-offset.x;
7452  box_q[2].y=polygon_primitive[n].point.y-offset.y;
7453  box_q[3].x=polygon_primitive[i].point.x-offset.x;
7454  box_q[3].y=polygon_primitive[i].point.y-offset.y;
7455  }
7456  if (fabs((double) (slope.p-slope.q)) < MagickEpsilon)
7457  {
7458  box_p[4]=box_p[1];
7459  box_q[4]=box_q[1];
7460  }
7461  else
7462  {
7463  box_p[4].x=(double) ((slope.p*box_p[0].x-box_p[0].y-slope.q*box_p[3].x+
7464  box_p[3].y)/(slope.p-slope.q));
7465  box_p[4].y=(double) (slope.p*(box_p[4].x-box_p[0].x)+box_p[0].y);
7466  box_q[4].x=(double) ((slope.p*box_q[0].x-box_q[0].y-slope.q*box_q[3].x+
7467  box_q[3].y)/(slope.p-slope.q));
7468  box_q[4].y=(double) (slope.p*(box_q[4].x-box_q[0].x)+box_q[0].y);
7469  }
7470  CheckPathExtent(MaxStrokePad,MaxStrokePad);
7471  dot_product=dx.q*dy.p-dx.p*dy.q;
7472  if (dot_product <= 0.0)
7473  switch (draw_info->linejoin)
7474  {
7475  case BevelJoin:
7476  {
7477  stroke_q[q++]=box_q[1];
7478  stroke_q[q++]=box_q[2];
7479  dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
7480  (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
7481  if (dot_product <= miterlimit)
7482  stroke_p[p++]=box_p[4];
7483  else
7484  {
7485  stroke_p[p++]=box_p[1];
7486  stroke_p[p++]=box_p[2];
7487  }
7488  break;
7489  }
7490  case MiterJoin:
7491  {
7492  dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
7493  (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
7494  if (dot_product <= miterlimit)
7495  {
7496  stroke_q[q++]=box_q[4];
7497  stroke_p[p++]=box_p[4];
7498  }
7499  else
7500  {
7501  stroke_q[q++]=box_q[1];
7502  stroke_q[q++]=box_q[2];
7503  stroke_p[p++]=box_p[1];
7504  stroke_p[p++]=box_p[2];
7505  }
7506  break;
7507  }
7508  case RoundJoin:
7509  {
7510  dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
7511  (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
7512  if (dot_product <= miterlimit)
7513  stroke_p[p++]=box_p[4];
7514  else
7515  {
7516  stroke_p[p++]=box_p[1];
7517  stroke_p[p++]=box_p[2];
7518  }
7519  center=polygon_primitive[n].point;
7520  theta.p=atan2(box_q[1].y-center.y,box_q[1].x-center.x);
7521  theta.q=atan2(box_q[2].y-center.y,box_q[2].x-center.x);
7522  if (theta.q < theta.p)
7523  theta.q+=2.0*MagickPI;
7524  arc_segments=(size_t) CastDoubleToLong(ceil((double) ((theta.q-
7525  theta.p)/(2.0*sqrt(PerceptibleReciprocal(mid))))));
7526  CheckPathExtent(MaxStrokePad,arc_segments+MaxStrokePad);
7527  stroke_q[q].x=box_q[1].x;
7528  stroke_q[q].y=box_q[1].y;
7529  q++;
7530  for (j=1; j < (ssize_t) arc_segments; j++)
7531  {
7532  delta_theta=(double) (j*(theta.q-theta.p)/arc_segments);
7533  stroke_q[q].x=(double) (center.x+mid*cos(fmod((double)
7534  (theta.p+delta_theta),DegreesToRadians(360.0))));
7535  stroke_q[q].y=(double) (center.y+mid*sin(fmod((double)
7536  (theta.p+delta_theta),DegreesToRadians(360.0))));
7537  q++;
7538  }
7539  stroke_q[q++]=box_q[2];
7540  break;
7541  }
7542  default:
7543  break;
7544  }
7545  else
7546  switch (draw_info->linejoin)
7547  {
7548  case BevelJoin:
7549  {
7550  stroke_p[p++]=box_p[1];
7551  stroke_p[p++]=box_p[2];
7552  dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
7553  (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
7554  if (dot_product <= miterlimit)
7555  stroke_q[q++]=box_q[4];
7556  else
7557  {
7558  stroke_q[q++]=box_q[1];
7559  stroke_q[q++]=box_q[2];
7560  }
7561  break;
7562  }
7563  case MiterJoin:
7564  {
7565  dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
7566  (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
7567  if (dot_product <= miterlimit)
7568  {
7569  stroke_q[q++]=box_q[4];
7570  stroke_p[p++]=box_p[4];
7571  }
7572  else
7573  {
7574  stroke_q[q++]=box_q[1];
7575  stroke_q[q++]=box_q[2];
7576  stroke_p[p++]=box_p[1];
7577  stroke_p[p++]=box_p[2];
7578  }
7579  break;
7580  }
7581  case RoundJoin:
7582  {
7583  dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
7584  (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
7585  if (dot_product <= miterlimit)
7586  stroke_q[q++]=box_q[4];
7587  else
7588  {
7589  stroke_q[q++]=box_q[1];
7590  stroke_q[q++]=box_q[2];
7591  }
7592  center=polygon_primitive[n].point;
7593  theta.p=atan2(box_p[1].y-center.y,box_p[1].x-center.x);
7594  theta.q=atan2(box_p[2].y-center.y,box_p[2].x-center.x);
7595  if (theta.p < theta.q)
7596  theta.p+=2.0*MagickPI;
7597  arc_segments=(size_t) CastDoubleToLong(ceil((double) ((theta.p-
7598  theta.q)/(2.0*sqrt((double) (PerceptibleReciprocal(mid)))))));
7599  CheckPathExtent(arc_segments+MaxStrokePad,MaxStrokePad);
7600  stroke_p[p++]=box_p[1];
7601  for (j=1; j < (ssize_t) arc_segments; j++)
7602  {
7603  delta_theta=(double) (j*(theta.q-theta.p)/arc_segments);
7604  stroke_p[p].x=(double) (center.x+mid*cos(fmod((double)
7605  (theta.p+delta_theta),DegreesToRadians(360.0))));
7606  stroke_p[p].y=(double) (center.y+mid*sin(fmod((double)
7607  (theta.p+delta_theta),DegreesToRadians(360.0))));
7608  p++;
7609  }
7610  stroke_p[p++]=box_p[2];
7611  break;
7612  }
7613  default:
7614  break;
7615  }
7616  slope.p=slope.q;
7617  inverse_slope.p=inverse_slope.q;
7618  box_p[0]=box_p[2];
7619  box_p[1]=box_p[3];
7620  box_q[0]=box_q[2];
7621  box_q[1]=box_q[3];
7622  dx.p=dx.q;
7623  dy.p=dy.q;
7624  n=i;
7625  }
7626  stroke_p[p++]=box_p[1];
7627  stroke_q[q++]=box_q[1];
7628  /*
7629  Trace stroked polygon.
7630  */
7631  stroke_polygon=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
7632  (p+q+2UL*closed_path+2UL),sizeof(*stroke_polygon));
7633  if (stroke_polygon == (PrimitiveInfo *) NULL)
7634  {
7635  (void) ThrowMagickException(exception,GetMagickModule(),
7636  ResourceLimitError,"MemoryAllocationFailed","`%s'","");
7637  stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p);
7638  stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q);
7639  polygon_primitive=(PrimitiveInfo *) RelinquishMagickMemory(
7640  polygon_primitive);
7641  return(stroke_polygon);
7642  }
7643  for (i=0; i < (ssize_t) p; i++)
7644  {
7645  stroke_polygon[i]=polygon_primitive[0];
7646  stroke_polygon[i].point=stroke_p[i];
7647  }
7648  if (closed_path != MagickFalse)
7649  {
7650  stroke_polygon[i]=polygon_primitive[0];
7651  stroke_polygon[i].point=stroke_polygon[0].point;
7652  i++;
7653  }
7654  for ( ; i < (ssize_t) (p+q+closed_path); i++)
7655  {
7656  stroke_polygon[i]=polygon_primitive[0];
7657  stroke_polygon[i].point=stroke_q[p+q+closed_path-(i+1)];
7658  }
7659  if (closed_path != MagickFalse)
7660  {
7661  stroke_polygon[i]=polygon_primitive[0];
7662  stroke_polygon[i].point=stroke_polygon[p+closed_path].point;
7663  i++;
7664  }
7665  stroke_polygon[i]=polygon_primitive[0];
7666  stroke_polygon[i].point=stroke_polygon[0].point;
7667  i++;
7668  stroke_polygon[i].primitive=UndefinedPrimitive;
7669  stroke_polygon[0].coordinates=(size_t) (p+q+2*closed_path+1);
7670  stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p);
7671  stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q);
7672  polygon_primitive=(PrimitiveInfo *) RelinquishMagickMemory(polygon_primitive);
7673  return(stroke_polygon);
7674 }
Definition: draw.c:145
Definition: image.h:152