2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6 % PPPP AAA IIIII N N TTTTT %
8 % PPPP AAAAA I N N N T %
13 % Methods to Paint on an Image %
20 % Copyright 1999-2019 ImageMagick Studio LLC, a non-profit organization %
21 % dedicated to making software imaging solutions freely available. %
23 % You may not use this file except in compliance with the License. You may %
24 % obtain a copy of the License at %
26 % https://imagemagick.org/script/license.php %
28 % Unless required by applicable law or agreed to in writing, software %
29 % distributed under the License is distributed on an "AS IS" BASIS, %
30 % WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. %
31 % See the License for the specific language governing permissions and %
32 % limitations under the License. %
34 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
42 #include "MagickCore/studio.h"
43 #include "MagickCore/artifact.h"
44 #include "MagickCore/channel.h"
45 #include "MagickCore/color.h"
46 #include "MagickCore/color-private.h"
47 #include "MagickCore/colorspace-private.h"
48 #include "MagickCore/composite.h"
49 #include "MagickCore/composite-private.h"
50 #include "MagickCore/draw.h"
51 #include "MagickCore/draw-private.h"
52 #include "MagickCore/exception.h"
53 #include "MagickCore/exception-private.h"
54 #include "MagickCore/gem.h"
55 #include "MagickCore/gem-private.h"
56 #include "MagickCore/monitor.h"
57 #include "MagickCore/monitor-private.h"
58 #include "MagickCore/option.h"
59 #include "MagickCore/paint.h"
60 #include "MagickCore/pixel-accessor.h"
61 #include "MagickCore/resource_.h"
62 #include "MagickCore/statistic.h"
63 #include "MagickCore/string_.h"
64 #include "MagickCore/string-private.h"
65 #include "MagickCore/thread-private.h"
68 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
72 % F l o o d f i l l P a i n t I m a g e %
76 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
78 % FloodfillPaintImage() changes the color value of any pixel that matches
79 % target and is an immediate neighbor. If the method FillToBorderMethod is
80 % specified, the color value is changed for any neighbor pixel that does not
81 % match the bordercolor member of image.
83 % By default target must match a particular pixel color exactly. However,
84 % in many cases two colors may differ by a small amount. The fuzz member of
85 % image defines how much tolerance is acceptable to consider two colors as
86 % the same. For example, set fuzz to 10 and the color red at intensities of
87 % 100 and 102 respectively are now interpreted as the same color for the
88 % purposes of the floodfill.
90 % The format of the FloodfillPaintImage method is:
92 % MagickBooleanType FloodfillPaintImage(Image *image,
93 % const DrawInfo *draw_info,const PixelInfo target,
94 % const ssize_t x_offset,const ssize_t y_offset,
95 % const MagickBooleanType invert,ExceptionInfo *exception)
97 % A description of each parameter follows:
101 % o draw_info: the draw info.
103 % o target: the RGB value of the target color.
105 % o x_offset,y_offset: the starting location of the operation.
107 % o invert: paint any pixel that does not match the target color.
109 % o exception: return any errors or warnings in this structure.
112 MagickExport MagickBooleanType FloodfillPaintImage(Image *image,
113 const DrawInfo *draw_info,const PixelInfo *target,const ssize_t x_offset,
114 const ssize_t y_offset,const MagickBooleanType invert,
115 ExceptionInfo *exception)
117 #define MaxStacksize 524288UL
118 #define PushSegmentStack(up,left,right,delta) \
120 if (s >= (segment_stack+MaxStacksize)) \
121 ThrowBinaryException(DrawError,"SegmentStackOverflow",image->filename) \
124 if ((((up)+(delta)) >= 0) && (((up)+(delta)) < (ssize_t) image->rows)) \
126 s->x1=(double) (left); \
127 s->y1=(double) (up); \
128 s->x2=(double) (right); \
129 s->y2=(double) (delta); \
167 Check boundary conditions.
169 assert(image != (Image *) NULL);
170 assert(image->signature == MagickCoreSignature);
171 if (image->debug != MagickFalse)
172 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
173 assert(draw_info != (DrawInfo *) NULL);
174 assert(draw_info->signature == MagickCoreSignature);
175 if ((x_offset < 0) || (x_offset >= (ssize_t) image->columns))
177 if ((y_offset < 0) || (y_offset >= (ssize_t) image->rows))
179 if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
181 if (IsGrayColorspace(image->colorspace) != MagickFalse)
182 (void) SetImageColorspace(image,sRGBColorspace,exception);
183 if ((image->alpha_trait == UndefinedPixelTrait) &&
184 (draw_info->fill.alpha_trait != UndefinedPixelTrait))
185 (void) SetImageAlpha(image,OpaqueAlpha,exception);
189 floodplane_image=CloneImage(image,0,0,MagickTrue,
191 if (floodplane_image == (Image *) NULL)
193 floodplane_image->alpha_trait=UndefinedPixelTrait;
194 floodplane_image->colorspace=GRAYColorspace;
195 (void) QueryColorCompliance("#000",AllCompliance,
196 &floodplane_image->background_color,exception);
197 (void) SetImageBackgroundColor(floodplane_image,exception);
198 segment_info=AcquireVirtualMemory(MaxStacksize,sizeof(*segment_stack));
199 if (segment_info == (MemoryInfo *) NULL)
201 floodplane_image=DestroyImage(floodplane_image);
202 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
205 segment_stack=(SegmentInfo *) GetVirtualMemoryBlob(segment_info);
207 Push initial segment on stack.
212 PushSegmentStack(y_offset,x_offset,x_offset,1);
213 PushSegmentStack(y_offset+1,x_offset,x_offset,-1);
214 GetPixelInfo(image,&pixel);
215 image_view=AcquireVirtualCacheView(image,exception);
216 floodplane_view=AcquireAuthenticCacheView(floodplane_image,exception);
217 while (s > segment_stack)
219 register const Quantum
229 Pop segment off stack.
234 offset=(ssize_t) s->y2;
235 y=(ssize_t) s->y1+offset;
237 Recolor neighboring pixels.
239 p=GetCacheViewVirtualPixels(image_view,0,y,(size_t) (x1+1),1,exception);
240 q=GetCacheViewAuthenticPixels(floodplane_view,0,y,(size_t) (x1+1),1,
242 if ((p == (const Quantum *) NULL) || (q == (Quantum *) NULL))
244 p+=x1*GetPixelChannels(image);
245 q+=x1*GetPixelChannels(floodplane_image);
246 for (x=x1; x >= 0; x--)
248 if (GetPixelGray(floodplane_image,q) != 0)
250 GetPixelInfoPixel(image,p,&pixel);
251 if (IsFuzzyEquivalencePixelInfo(&pixel,target) == invert)
253 SetPixelGray(floodplane_image,QuantumRange,q);
254 p-=GetPixelChannels(image);
255 q-=GetPixelChannels(floodplane_image);
257 if (SyncCacheViewAuthenticPixels(floodplane_view,exception) == MagickFalse)
259 skip=x >= x1 ? MagickTrue : MagickFalse;
260 if (skip == MagickFalse)
264 PushSegmentStack(y,start,x1-1,-offset);
269 if (skip == MagickFalse)
271 if (x < (ssize_t) image->columns)
273 p=GetCacheViewVirtualPixels(image_view,x,y,image->columns-x,1,
275 q=GetCacheViewAuthenticPixels(floodplane_view,x,y,image->columns-
277 if ((p == (const Quantum *) NULL) || (q == (Quantum *) NULL))
279 for ( ; x < (ssize_t) image->columns; x++)
281 if (GetPixelGray(floodplane_image,q) != 0)
283 GetPixelInfoPixel(image,p,&pixel);
284 if (IsFuzzyEquivalencePixelInfo(&pixel,target) == invert)
286 SetPixelGray(floodplane_image,QuantumRange,q);
287 p+=GetPixelChannels(image);
288 q+=GetPixelChannels(floodplane_image);
290 status=SyncCacheViewAuthenticPixels(floodplane_view,exception);
291 if (status == MagickFalse)
294 PushSegmentStack(y,start,x-1,offset);
296 PushSegmentStack(y,x2+1,x-1,-offset);
302 p=GetCacheViewVirtualPixels(image_view,x,y,(size_t) (x2-x+1),1,
304 q=GetCacheViewAuthenticPixels(floodplane_view,x,y,(size_t) (x2-x+1),1,
306 if ((p == (const Quantum *) NULL) || (q == (Quantum *) NULL))
308 for ( ; x <= x2; x++)
310 if (GetPixelGray(floodplane_image,q) != 0)
312 GetPixelInfoPixel(image,p,&pixel);
313 if (IsFuzzyEquivalencePixelInfo(&pixel,target) != invert)
315 p+=GetPixelChannels(image);
316 q+=GetPixelChannels(floodplane_image);
323 for (y=0; y < (ssize_t) image->rows; y++)
325 register const Quantum
335 Tile fill color onto floodplane.
337 if (status == MagickFalse)
339 p=GetCacheViewVirtualPixels(floodplane_view,0,y,image->columns,1,exception);
340 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
341 if ((p == (const Quantum *) NULL) || (q == (Quantum *) NULL))
346 for (x=0; x < (ssize_t) image->columns; x++)
348 if (GetPixelGray(floodplane_image,p) != 0)
350 GetFillColor(draw_info,x,y,&fill_color,exception);
351 SetPixelViaPixelInfo(image,&fill_color,q);
353 p+=GetPixelChannels(floodplane_image);
354 q+=GetPixelChannels(image);
356 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
359 floodplane_view=DestroyCacheView(floodplane_view);
360 image_view=DestroyCacheView(image_view);
361 segment_info=RelinquishVirtualMemory(segment_info);
362 floodplane_image=DestroyImage(floodplane_image);
367 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
371 + G r a d i e n t I m a g e %
375 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
377 % GradientImage() applies a continuously smooth color transitions along a
378 % vector from one color to another.
380 % Note, the interface of this method will change in the future to support
381 % more than one transistion.
383 % The format of the GradientImage method is:
385 % MagickBooleanType GradientImage(Image *image,const GradientType type,
386 % const SpreadMethod method,const PixelInfo *start_color,
387 % const PixelInfo *stop_color,ExceptionInfo *exception)
389 % A description of each parameter follows:
391 % o image: the image.
393 % o type: the gradient type: linear or radial.
395 % o spread: the gradient spread meathod: pad, reflect, or repeat.
397 % o start_color: the start color.
399 % o stop_color: the stop color.
401 % o exception: return any errors or warnings in this structure.
404 MagickExport MagickBooleanType GradientImage(Image *image,
405 const GradientType type,const SpreadMethod method,const StopInfo *stops,
406 const size_t number_stops,ExceptionInfo *exception)
421 Set gradient start-stop end points.
423 assert(image != (const Image *) NULL);
424 assert(image->signature == MagickCoreSignature);
425 if (image->debug != MagickFalse)
426 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
427 assert(stops != (const StopInfo *) NULL);
428 assert(number_stops > 0);
429 draw_info=AcquireDrawInfo();
430 gradient=(&draw_info->gradient);
432 gradient->bounding_box.width=image->columns;
433 gradient->bounding_box.height=image->rows;
434 artifact=GetImageArtifact(image,"gradient:bounding-box");
435 if (artifact != (const char *) NULL)
436 (void) ParseAbsoluteGeometry(artifact,&gradient->bounding_box);
437 gradient->gradient_vector.x2=(double) image->columns-1;
438 gradient->gradient_vector.y2=(double) image->rows-1;
439 artifact=GetImageArtifact(image,"gradient:direction");
440 if (artifact != (const char *) NULL)
445 direction=(GravityType) ParseCommandOption(MagickGravityOptions,
446 MagickFalse,artifact);
449 case NorthWestGravity:
451 gradient->gradient_vector.x1=(double) image->columns-1;
452 gradient->gradient_vector.y1=(double) image->rows-1;
453 gradient->gradient_vector.x2=0.0;
454 gradient->gradient_vector.y2=0.0;
459 gradient->gradient_vector.x1=0.0;
460 gradient->gradient_vector.y1=(double) image->rows-1;
461 gradient->gradient_vector.x2=0.0;
462 gradient->gradient_vector.y2=0.0;
465 case NorthEastGravity:
467 gradient->gradient_vector.x1=0.0;
468 gradient->gradient_vector.y1=(double) image->rows-1;
469 gradient->gradient_vector.x2=(double) image->columns-1;
470 gradient->gradient_vector.y2=0.0;
475 gradient->gradient_vector.x1=(double) image->columns-1;
476 gradient->gradient_vector.y1=0.0;
477 gradient->gradient_vector.x2=0.0;
478 gradient->gradient_vector.y2=0.0;
483 gradient->gradient_vector.x1=0.0;
484 gradient->gradient_vector.y1=0.0;
485 gradient->gradient_vector.x2=(double) image->columns-1;
486 gradient->gradient_vector.y2=0.0;
489 case SouthWestGravity:
491 gradient->gradient_vector.x1=(double) image->columns-1;
492 gradient->gradient_vector.y1=0.0;
493 gradient->gradient_vector.x2=0.0;
494 gradient->gradient_vector.y2=(double) image->rows-1;
499 gradient->gradient_vector.x1=0.0;
500 gradient->gradient_vector.y1=0.0;
501 gradient->gradient_vector.x2=0.0;
502 gradient->gradient_vector.y2=(double) image->columns-1;
505 case SouthEastGravity:
507 gradient->gradient_vector.x1=0.0;
508 gradient->gradient_vector.y1=0.0;
509 gradient->gradient_vector.x2=(double) image->columns-1;
510 gradient->gradient_vector.y2=(double) image->rows-1;
517 artifact=GetImageArtifact(image,"gradient:angle");
518 if (artifact != (const char *) NULL)
519 gradient->angle=StringToDouble(artifact,(char **) NULL);
520 artifact=GetImageArtifact(image,"gradient:vector");
521 if (artifact != (const char *) NULL)
522 (void) sscanf(artifact,"%lf%*[ ,]%lf%*[ ,]%lf%*[ ,]%lf",
523 &gradient->gradient_vector.x1,&gradient->gradient_vector.y1,
524 &gradient->gradient_vector.x2,&gradient->gradient_vector.y2);
525 if ((GetImageArtifact(image,"gradient:angle") == (const char *) NULL) &&
526 (GetImageArtifact(image,"gradient:direction") == (const char *) NULL) &&
527 (GetImageArtifact(image,"gradient:extent") == (const char *) NULL) &&
528 (GetImageArtifact(image,"gradient:vector") == (const char *) NULL))
529 if ((type == LinearGradient) && (gradient->gradient_vector.y2 != 0.0))
530 gradient->gradient_vector.x2=0.0;
531 gradient->center.x=(double) gradient->gradient_vector.x2/2.0;
532 gradient->center.y=(double) gradient->gradient_vector.y2/2.0;
533 artifact=GetImageArtifact(image,"gradient:center");
534 if (artifact != (const char *) NULL)
535 (void) sscanf(artifact,"%lf%*[ ,]%lf",&gradient->center.x,
536 &gradient->center.y);
537 artifact=GetImageArtifact(image,"gradient:angle");
538 if ((type == LinearGradient) && (artifact != (const char *) NULL))
546 Reference https://drafts.csswg.org/css-images-3/#linear-gradients.
548 sine=sin((double) DegreesToRadians(gradient->angle-90.0));
549 cosine=cos((double) DegreesToRadians(gradient->angle-90.0));
550 distance=fabs((double) (image->columns-1.0)*cosine)+
551 fabs((double) (image->rows-1.0)*sine);
552 gradient->gradient_vector.x1=0.5*((image->columns-1.0)-distance*cosine);
553 gradient->gradient_vector.y1=0.5*((image->rows-1.0)-distance*sine);
554 gradient->gradient_vector.x2=0.5*((image->columns-1.0)+distance*cosine);
555 gradient->gradient_vector.y2=0.5*((image->rows-1.0)+distance*sine);
557 gradient->radii.x=(double) MagickMax((image->columns-1.0),(image->rows-1.0))/
559 gradient->radii.y=gradient->radii.x;
560 artifact=GetImageArtifact(image,"gradient:extent");
561 if (artifact != (const char *) NULL)
563 if (LocaleCompare(artifact,"Circle") == 0)
565 gradient->radii.x=(double) MagickMax((image->columns-1.0),
566 (image->rows-1.0))/2.0;
567 gradient->radii.y=gradient->radii.x;
569 if (LocaleCompare(artifact,"Diagonal") == 0)
571 gradient->radii.x=(double) (sqrt((double) (image->columns-1.0)*
572 (image->columns-1.0)+(image->rows-1.0)*(image->rows-1.0)))/2.0;
573 gradient->radii.y=gradient->radii.x;
575 if (LocaleCompare(artifact,"Ellipse") == 0)
577 gradient->radii.x=(double) (image->columns-1.0)/2.0;
578 gradient->radii.y=(double) (image->rows-1.0)/2.0;
580 if (LocaleCompare(artifact,"Maximum") == 0)
582 gradient->radii.x=(double) MagickMax((image->columns-1.0),
583 (image->rows-1.0))/2.0;
584 gradient->radii.y=gradient->radii.x;
586 if (LocaleCompare(artifact,"Minimum") == 0)
588 gradient->radii.x=(double) (MagickMin((image->columns-1.0),
589 (image->rows-1.0)))/2.0;
590 gradient->radii.y=gradient->radii.x;
593 artifact=GetImageArtifact(image,"gradient:radii");
594 if (artifact != (const char *) NULL)
595 (void) sscanf(artifact,"%lf%*[ ,]%lf",&gradient->radii.x,
597 gradient->radius=MagickMax(gradient->radii.x,gradient->radii.y);
598 gradient->spread=method;
600 Define the gradient to fill between the stops.
602 gradient->number_stops=number_stops;
603 gradient->stops=(StopInfo *) AcquireQuantumMemory(gradient->number_stops,
604 sizeof(*gradient->stops));
605 if (gradient->stops == (StopInfo *) NULL)
606 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
608 (void) memcpy(gradient->stops,stops,(size_t) number_stops*
611 Draw a gradient on the image.
613 status=DrawGradientImage(image,draw_info,exception);
614 draw_info=DestroyDrawInfo(draw_info);
619 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
623 % O i l P a i n t I m a g e %
627 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
629 % OilPaintImage() applies a special effect filter that simulates an oil
630 % painting. Each pixel is replaced by the most frequent color occurring
631 % in a circular region defined by radius.
633 % The format of the OilPaintImage method is:
635 % Image *OilPaintImage(const Image *image,const double radius,
636 % const double sigma,ExceptionInfo *exception)
638 % A description of each parameter follows:
640 % o image: the image.
642 % o radius: the radius of the circular neighborhood.
644 % o sigma: the standard deviation of the Gaussian, in pixels.
646 % o exception: return any errors or warnings in this structure.
650 static size_t **DestroyHistogramThreadSet(size_t **histogram)
655 assert(histogram != (size_t **) NULL);
656 for (i=0; i < (ssize_t) GetMagickResourceLimit(ThreadResource); i++)
657 if (histogram[i] != (size_t *) NULL)
658 histogram[i]=(size_t *) RelinquishMagickMemory(histogram[i]);
659 histogram=(size_t **) RelinquishMagickMemory(histogram);
663 static size_t **AcquireHistogramThreadSet(const size_t count)
672 number_threads=(size_t) GetMagickResourceLimit(ThreadResource);
673 histogram=(size_t **) AcquireQuantumMemory(number_threads,sizeof(*histogram));
674 if (histogram == (size_t **) NULL)
675 return((size_t **) NULL);
676 (void) memset(histogram,0,number_threads*sizeof(*histogram));
677 for (i=0; i < (ssize_t) number_threads; i++)
679 histogram[i]=(size_t *) AcquireQuantumMemory(count,sizeof(**histogram));
680 if (histogram[i] == (size_t *) NULL)
681 return(DestroyHistogramThreadSet(histogram));
686 MagickExport Image *OilPaintImage(const Image *image,const double radius,
687 const double sigma,ExceptionInfo *exception)
689 #define NumberPaintBins 256
690 #define OilPaintImageTag "OilPaint/Image"
715 Initialize painted image attributes.
717 assert(image != (const Image *) NULL);
718 assert(image->signature == MagickCoreSignature);
719 if (image->debug != MagickFalse)
720 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
721 assert(exception != (ExceptionInfo *) NULL);
722 assert(exception->signature == MagickCoreSignature);
723 width=GetOptimalKernelWidth2D(radius,sigma);
724 linear_image=CloneImage(image,0,0,MagickTrue,exception);
725 paint_image=CloneImage(image,0,0,MagickTrue,exception);
726 if ((linear_image == (Image *) NULL) || (paint_image == (Image *) NULL))
728 if (linear_image != (Image *) NULL)
729 linear_image=DestroyImage(linear_image);
730 if (paint_image != (Image *) NULL)
731 linear_image=DestroyImage(paint_image);
732 return((Image *) NULL);
734 if (SetImageStorageClass(paint_image,DirectClass,exception) == MagickFalse)
736 linear_image=DestroyImage(linear_image);
737 paint_image=DestroyImage(paint_image);
738 return((Image *) NULL);
740 histograms=AcquireHistogramThreadSet(NumberPaintBins);
741 if (histograms == (size_t **) NULL)
743 linear_image=DestroyImage(linear_image);
744 paint_image=DestroyImage(paint_image);
745 ThrowImageException(ResourceLimitError,"MemoryAllocationFailed");
752 center=(ssize_t) GetPixelChannels(linear_image)*(linear_image->columns+width)*
753 (width/2L)+GetPixelChannels(linear_image)*(width/2L);
754 image_view=AcquireVirtualCacheView(linear_image,exception);
755 paint_view=AcquireAuthenticCacheView(paint_image,exception);
756 #if defined(MAGICKCORE_OPENMP_SUPPORT)
757 #pragma omp parallel for schedule(static) shared(progress,status) \
758 magick_number_threads(linear_image,paint_image,linear_image->rows,1)
760 for (y=0; y < (ssize_t) linear_image->rows; y++)
762 register const Quantum
774 if (status == MagickFalse)
776 p=GetCacheViewVirtualPixels(image_view,-((ssize_t) width/2L),y-(ssize_t)
777 (width/2L),linear_image->columns+width,width,exception);
778 q=QueueCacheViewAuthenticPixels(paint_view,0,y,paint_image->columns,1,
780 if ((p == (const Quantum *) NULL) || (q == (Quantum *) NULL))
785 histogram=histograms[GetOpenMPThreadId()];
786 for (x=0; x < (ssize_t) linear_image->columns; x++)
802 Assign most frequent color.
807 (void) memset(histogram,0,NumberPaintBins* sizeof(*histogram));
808 for (v=0; v < (ssize_t) width; v++)
810 for (u=0; u < (ssize_t) width; u++)
812 n=(ssize_t) ScaleQuantumToChar(ClampToQuantum(GetPixelIntensity(
813 linear_image,p+GetPixelChannels(linear_image)*(u+k))));
815 if (histogram[n] > count)
821 k+=(ssize_t) (linear_image->columns+width);
823 for (i=0; i < (ssize_t) GetPixelChannels(linear_image); i++)
825 PixelChannel channel = GetPixelChannelChannel(linear_image,i);
826 PixelTrait traits = GetPixelChannelTraits(linear_image,channel);
827 PixelTrait paint_traits=GetPixelChannelTraits(paint_image,channel);
828 if ((traits == UndefinedPixelTrait) ||
829 (paint_traits == UndefinedPixelTrait))
831 if ((paint_traits & CopyPixelTrait) != 0)
833 SetPixelChannel(paint_image,channel,p[center+i],q);
836 SetPixelChannel(paint_image,channel,p[j*GetPixelChannels(linear_image)+
839 p+=GetPixelChannels(linear_image);
840 q+=GetPixelChannels(paint_image);
842 if (SyncCacheViewAuthenticPixels(paint_view,exception) == MagickFalse)
844 if (linear_image->progress_monitor != (MagickProgressMonitor) NULL)
849 #if defined(MAGICKCORE_OPENMP_SUPPORT)
853 proceed=SetImageProgress(linear_image,OilPaintImageTag,progress,
855 if (proceed == MagickFalse)
859 paint_view=DestroyCacheView(paint_view);
860 image_view=DestroyCacheView(image_view);
861 histograms=DestroyHistogramThreadSet(histograms);
862 linear_image=DestroyImage(linear_image);
863 if (status == MagickFalse)
864 paint_image=DestroyImage(paint_image);
869 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
873 % O p a q u e P a i n t I m a g e %
877 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
879 % OpaquePaintImage() changes any pixel that matches color with the color
880 % defined by fill argument.
882 % By default color must match a particular pixel color exactly. However, in
883 % many cases two colors may differ by a small amount. Fuzz defines how much
884 % tolerance is acceptable to consider two colors as the same. For example,
885 % set fuzz to 10 and the color red at intensities of 100 and 102 respectively
886 % are now interpreted as the same color.
888 % The format of the OpaquePaintImage method is:
890 % MagickBooleanType OpaquePaintImage(Image *image,const PixelInfo *target,
891 % const PixelInfo *fill,const MagickBooleanType invert,
892 % ExceptionInfo *exception)
894 % A description of each parameter follows:
896 % o image: the image.
898 % o target: the RGB value of the target color.
900 % o fill: the replacement color.
902 % o invert: paint any pixel that does not match the target color.
904 % o exception: return any errors or warnings in this structure.
907 MagickExport MagickBooleanType OpaquePaintImage(Image *image,
908 const PixelInfo *target,const PixelInfo *fill,const MagickBooleanType invert,
909 ExceptionInfo *exception)
911 #define OpaquePaintImageTag "Opaque/Image"
930 assert(image != (Image *) NULL);
931 assert(image->signature == MagickCoreSignature);
932 assert(target != (PixelInfo *) NULL);
933 assert(fill != (PixelInfo *) NULL);
934 if (image->debug != MagickFalse)
935 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
936 if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
938 ConformPixelInfo(image,fill,&conform_fill,exception);
939 ConformPixelInfo(image,target,&conform_target,exception);
941 Make image color opaque.
945 GetPixelInfo(image,&zero);
946 image_view=AcquireAuthenticCacheView(image,exception);
947 #if defined(MAGICKCORE_OPENMP_SUPPORT)
948 #pragma omp parallel for schedule(static) shared(progress,status) \
949 magick_number_threads(image,image,image->rows,1)
951 for (y=0; y < (ssize_t) image->rows; y++)
962 if (status == MagickFalse)
964 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
965 if (q == (Quantum *) NULL)
971 for (x=0; x < (ssize_t) image->columns; x++)
973 GetPixelInfoPixel(image,q,&pixel);
974 if (IsFuzzyEquivalencePixelInfo(&pixel,&conform_target) != invert)
979 traits=GetPixelChannelTraits(image,RedPixelChannel);
980 if ((traits & UpdatePixelTrait) != 0)
981 SetPixelRed(image,(Quantum) conform_fill.red,q);
982 traits=GetPixelChannelTraits(image,GreenPixelChannel);
983 if ((traits & UpdatePixelTrait) != 0)
984 SetPixelGreen(image,(Quantum) conform_fill.green,q);
985 traits=GetPixelChannelTraits(image,BluePixelChannel);
986 if ((traits & UpdatePixelTrait) != 0)
987 SetPixelBlue(image,(Quantum) conform_fill.blue,q);
988 traits=GetPixelChannelTraits(image,BlackPixelChannel);
989 if ((traits & UpdatePixelTrait) != 0)
990 SetPixelBlack(image,(Quantum) conform_fill.black,q);
991 traits=GetPixelChannelTraits(image,AlphaPixelChannel);
992 if ((traits & UpdatePixelTrait) != 0)
993 SetPixelAlpha(image,(Quantum) conform_fill.alpha,q);
995 q+=GetPixelChannels(image);
997 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
999 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1004 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1008 proceed=SetImageProgress(image,OpaquePaintImageTag,progress,
1010 if (proceed == MagickFalse)
1014 image_view=DestroyCacheView(image_view);
1019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1023 % T r a n s p a r e n t P a i n t I m a g e %
1027 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1029 % TransparentPaintImage() changes the opacity value associated with any pixel
1030 % that matches color to the value defined by opacity.
1032 % By default color must match a particular pixel color exactly. However, in
1033 % many cases two colors may differ by a small amount. Fuzz defines how much
1034 % tolerance is acceptable to consider two colors as the same. For example,
1035 % set fuzz to 10 and the color red at intensities of 100 and 102 respectively
1036 % are now interpreted as the same color.
1038 % The format of the TransparentPaintImage method is:
1040 % MagickBooleanType TransparentPaintImage(Image *image,
1041 % const PixelInfo *target,const Quantum opacity,
1042 % const MagickBooleanType invert,ExceptionInfo *exception)
1044 % A description of each parameter follows:
1046 % o image: the image.
1048 % o target: the target color.
1050 % o opacity: the replacement opacity value.
1052 % o invert: paint any pixel that does not match the target color.
1054 % o exception: return any errors or warnings in this structure.
1057 MagickExport MagickBooleanType TransparentPaintImage(Image *image,
1058 const PixelInfo *target,const Quantum opacity,const MagickBooleanType invert,
1059 ExceptionInfo *exception)
1061 #define TransparentPaintImageTag "Transparent/Image"
1078 assert(image != (Image *) NULL);
1079 assert(image->signature == MagickCoreSignature);
1080 assert(target != (PixelInfo *) NULL);
1081 if (image->debug != MagickFalse)
1082 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1083 if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
1084 return(MagickFalse);
1085 if (image->alpha_trait == UndefinedPixelTrait)
1086 (void) SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
1088 Make image color transparent.
1092 GetPixelInfo(image,&zero);
1093 image_view=AcquireAuthenticCacheView(image,exception);
1094 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1095 #pragma omp parallel for schedule(static) shared(progress,status) \
1096 magick_number_threads(image,image,image->rows,1)
1098 for (y=0; y < (ssize_t) image->rows; y++)
1109 if (status == MagickFalse)
1111 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
1112 if (q == (Quantum *) NULL)
1118 for (x=0; x < (ssize_t) image->columns; x++)
1120 GetPixelInfoPixel(image,q,&pixel);
1121 if (IsFuzzyEquivalencePixelInfo(&pixel,target) != invert)
1122 SetPixelAlpha(image,opacity,q);
1123 q+=GetPixelChannels(image);
1125 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1127 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1132 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1136 proceed=SetImageProgress(image,TransparentPaintImageTag,progress,
1138 if (proceed == MagickFalse)
1142 image_view=DestroyCacheView(image_view);
1147 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1151 % T r a n s p a r e n t P a i n t I m a g e C h r o m a %
1155 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1157 % TransparentPaintImageChroma() changes the opacity value associated with any
1158 % pixel that matches color to the value defined by opacity.
1160 % As there is one fuzz value for the all the channels, TransparentPaintImage()
1161 % is not suitable for the operations like chroma, where the tolerance for
1162 % similarity of two color component (RGB) can be different. Thus we define
1163 % this method to take two target pixels (one low and one high) and all the
1164 % pixels of an image which are lying between these two pixels are made
1167 % The format of the TransparentPaintImageChroma method is:
1169 % MagickBooleanType TransparentPaintImageChroma(Image *image,
1170 % const PixelInfo *low,const PixelInfo *high,const Quantum opacity,
1171 % const MagickBooleanType invert,ExceptionInfo *exception)
1173 % A description of each parameter follows:
1175 % o image: the image.
1177 % o low: the low target color.
1179 % o high: the high target color.
1181 % o opacity: the replacement opacity value.
1183 % o invert: paint any pixel that does not match the target color.
1185 % o exception: return any errors or warnings in this structure.
1188 MagickExport MagickBooleanType TransparentPaintImageChroma(Image *image,
1189 const PixelInfo *low,const PixelInfo *high,const Quantum opacity,
1190 const MagickBooleanType invert,ExceptionInfo *exception)
1192 #define TransparentPaintImageTag "Transparent/Image"
1206 assert(image != (Image *) NULL);
1207 assert(image->signature == MagickCoreSignature);
1208 assert(high != (PixelInfo *) NULL);
1209 assert(low != (PixelInfo *) NULL);
1210 if (image->debug != MagickFalse)
1211 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1212 if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
1213 return(MagickFalse);
1214 if (image->alpha_trait == UndefinedPixelTrait)
1215 (void) SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
1217 Make image color transparent.
1221 image_view=AcquireAuthenticCacheView(image,exception);
1222 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1223 #pragma omp parallel for schedule(static) shared(progress,status) \
1224 magick_number_threads(image,image,image->rows,1)
1226 for (y=0; y < (ssize_t) image->rows; y++)
1240 if (status == MagickFalse)
1242 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
1243 if (q == (Quantum *) NULL)
1248 GetPixelInfo(image,&pixel);
1249 for (x=0; x < (ssize_t) image->columns; x++)
1251 GetPixelInfoPixel(image,q,&pixel);
1252 match=((pixel.red >= low->red) && (pixel.red <= high->red) &&
1253 (pixel.green >= low->green) && (pixel.green <= high->green) &&
1254 (pixel.blue >= low->blue) && (pixel.blue <= high->blue)) ? MagickTrue :
1256 if (match != invert)
1257 SetPixelAlpha(image,opacity,q);
1258 q+=GetPixelChannels(image);
1260 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1262 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1267 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1271 proceed=SetImageProgress(image,TransparentPaintImageTag,progress,
1273 if (proceed == MagickFalse)
1277 image_view=DestroyCacheView(image_view);