2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6 % EEEEE N N H H AAA N N CCCC EEEEE %
7 % E NN N H H A A NN N C E %
8 % EEE N N N HHHHH AAAAA N N N C EEE %
9 % E N NN H H A A N NN C E %
10 % EEEEE N N H H A A N N CCCC EEEEE %
13 % MagickCore Image Enhancement Methods %
20 % Copyright 1999-2012 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 % http://www.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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
43 #include "MagickCore/studio.h"
44 #include "MagickCore/artifact.h"
45 #include "MagickCore/cache.h"
46 #include "MagickCore/cache-view.h"
47 #include "MagickCore/color.h"
48 #include "MagickCore/color-private.h"
49 #include "MagickCore/colorspace.h"
50 #include "MagickCore/composite-private.h"
51 #include "MagickCore/enhance.h"
52 #include "MagickCore/exception.h"
53 #include "MagickCore/exception-private.h"
54 #include "MagickCore/fx.h"
55 #include "MagickCore/gem.h"
56 #include "MagickCore/gem-private.h"
57 #include "MagickCore/geometry.h"
58 #include "MagickCore/histogram.h"
59 #include "MagickCore/image.h"
60 #include "MagickCore/image-private.h"
61 #include "MagickCore/memory_.h"
62 #include "MagickCore/monitor.h"
63 #include "MagickCore/monitor-private.h"
64 #include "MagickCore/option.h"
65 #include "MagickCore/pixel-accessor.h"
66 #include "MagickCore/quantum.h"
67 #include "MagickCore/quantum-private.h"
68 #include "MagickCore/resample.h"
69 #include "MagickCore/resample-private.h"
70 #include "MagickCore/statistic.h"
71 #include "MagickCore/string_.h"
72 #include "MagickCore/string-private.h"
73 #include "MagickCore/thread-private.h"
74 #include "MagickCore/token.h"
75 #include "MagickCore/xml-tree.h"
76 #include "MagickCore/xml-tree-private.h"
79 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
83 % A u t o G a m m a I m a g e %
87 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
89 % AutoGammaImage() extract the 'mean' from the image and adjust the image
90 % to try make set its gamma appropriatally.
92 % The format of the AutoGammaImage method is:
94 % MagickBooleanType AutoGammaImage(Image *image,ExceptionInfo *exception)
96 % A description of each parameter follows:
98 % o image: The image to auto-level
100 % o exception: return any errors or warnings in this structure.
103 MagickExport MagickBooleanType AutoGammaImage(Image *image,
104 ExceptionInfo *exception)
119 if (image->channel_mask == DefaultChannels)
122 Apply gamma correction equally across all given channels.
124 (void) GetImageMean(image,&mean,&sans,exception);
125 gamma=log(mean*QuantumScale)/log_mean;
126 return(LevelImage(image,0.0,(double) QuantumRange,gamma,exception));
129 Auto-gamma each channel separately.
132 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
143 channel=GetPixelChannelMapChannel(image,i);
144 traits=GetPixelChannelMapTraits(image,channel);
145 if ((traits & UpdatePixelTrait) == 0)
147 channel_mask=SetPixelChannelMask(image,(ChannelType) (1 << i));
148 status=GetImageMean(image,&mean,&sans,exception);
149 gamma=log(mean*QuantumScale)/log_mean;
150 status&=LevelImage(image,0.0,(double) QuantumRange,gamma,exception);
151 (void) SetPixelChannelMask(image,channel_mask);
152 if (status == MagickFalse)
155 return(status != 0 ? MagickTrue : MagickFalse);
159 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
163 % A u t o L e v e l I m a g e %
167 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
169 % AutoLevelImage() adjusts the levels of a particular image channel by
170 % scaling the minimum and maximum values to the full quantum range.
172 % The format of the LevelImage method is:
174 % MagickBooleanType AutoLevelImage(Image *image,ExceptionInfo *exception)
176 % A description of each parameter follows:
178 % o image: The image to auto-level
180 % o exception: return any errors or warnings in this structure.
183 MagickExport MagickBooleanType AutoLevelImage(Image *image,
184 ExceptionInfo *exception)
186 return(MinMaxStretchImage(image,0.0,0.0,1.0,exception));
190 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
194 % B r i g h t n e s s C o n t r a s t I m a g e %
198 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
200 % BrightnessContrastImage() changes the brightness and/or contrast of an
201 % image. It converts the brightness and contrast parameters into slope and
202 % intercept and calls a polynomical function to apply to the image.
204 % The format of the BrightnessContrastImage method is:
206 % MagickBooleanType BrightnessContrastImage(Image *image,
207 % const double brightness,const double contrast,ExceptionInfo *exception)
209 % A description of each parameter follows:
211 % o image: the image.
213 % o brightness: the brightness percent (-100 .. 100).
215 % o contrast: the contrast percent (-100 .. 100).
217 % o exception: return any errors or warnings in this structure.
220 MagickExport MagickBooleanType BrightnessContrastImage(Image *image,
221 const double brightness,const double contrast,ExceptionInfo *exception)
223 #define BrightnessContastImageTag "BrightnessContast/Image"
235 Compute slope and intercept.
237 assert(image != (Image *) NULL);
238 assert(image->signature == MagickSignature);
239 if (image->debug != MagickFalse)
240 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
242 slope=tan((double) (MagickPI*(alpha/100.0+1.0)/4.0));
245 intercept=brightness/100.0+((100-brightness)/200.0)*(1.0-slope);
246 coefficients[0]=slope;
247 coefficients[1]=intercept;
248 status=FunctionImage(image,PolynomialFunction,2,coefficients,exception);
253 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
257 % C l u t I m a g e %
261 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
263 % ClutImage() replaces each color value in the given image, by using it as an
264 % index to lookup a replacement color value in a Color Look UP Table in the
265 % form of an image. The values are extracted along a diagonal of the CLUT
266 % image so either a horizontal or vertial gradient image can be used.
268 % Typically this is used to either re-color a gray-scale image according to a
269 % color gradient in the CLUT image, or to perform a freeform histogram
270 % (level) adjustment according to the (typically gray-scale) gradient in the
273 % When the 'channel' mask includes the matte/alpha transparency channel but
274 % one image has no such channel it is assumed that that image is a simple
275 % gray-scale image that will effect the alpha channel values, either for
276 % gray-scale coloring (with transparent or semi-transparent colors), or
277 % a histogram adjustment of existing alpha channel values. If both images
278 % have matte channels, direct and normal indexing is applied, which is rarely
281 % The format of the ClutImage method is:
283 % MagickBooleanType ClutImage(Image *image,Image *clut_image,
284 % const PixelInterpolateMethod method,ExceptionInfo *exception)
286 % A description of each parameter follows:
288 % o image: the image, which is replaced by indexed CLUT values
290 % o clut_image: the color lookup table image for replacement color values.
292 % o method: the pixel interpolation method.
294 % o exception: return any errors or warnings in this structure.
297 MagickExport MagickBooleanType ClutImage(Image *image,const Image *clut_image,
298 const PixelInterpolateMethod method,ExceptionInfo *exception)
300 #define ClutImageTag "Clut/Image"
322 assert(image != (Image *) NULL);
323 assert(image->signature == MagickSignature);
324 if (image->debug != MagickFalse)
325 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
326 assert(clut_image != (Image *) NULL);
327 assert(clut_image->signature == MagickSignature);
328 if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
330 clut_map=(double *) AcquireQuantumMemory(MaxMap+1UL,GetPixelChannels(image)*
332 if (clut_map == (double *) NULL)
333 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
340 adjust=(ssize_t) (clut_image->interpolate == IntegerInterpolatePixel ? 0 : 1);
341 clut_view=AcquireCacheView(clut_image);
342 #if defined(MAGICKCORE_OPENMP_SUPPORT)
343 #pragma omp parallel for schedule(static,4)
345 for (x=0; x <= (ssize_t) MaxMap; x++)
350 for (i=0; i < (ssize_t) GetPixelChannels(clut_image); i++)
351 (void) InterpolatePixelChannel(clut_image,clut_view,(PixelChannel) i,
352 method,QuantumScale*x*(clut_image->columns-adjust),QuantumScale*x*
353 (clut_image->rows-adjust),clut_map+x*GetPixelChannels(clut_image)+i,
356 clut_view=DestroyCacheView(clut_view);
357 image_view=AcquireCacheView(image);
358 #if defined(MAGICKCORE_OPENMP_SUPPORT)
359 #pragma omp parallel for schedule(static,4) shared(progress,status)
361 for (y=0; y < (ssize_t) image->rows; y++)
369 if (status == MagickFalse)
371 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
372 if (q == (Quantum *) NULL)
377 for (x=0; x < (ssize_t) image->columns; x++)
382 for (i=0; i < (ssize_t) GetPixelChannels(clut_image); i++)
391 channel=GetPixelChannelMapChannel(clut_image,i);
392 clut_traits=GetPixelChannelMapTraits(clut_image,channel);
393 traits=GetPixelChannelMapTraits(clut_image,channel);
394 if ((traits == UndefinedPixelTrait) ||
395 (clut_traits == UndefinedPixelTrait) ||
396 ((traits & UpdatePixelTrait) == 0))
398 SetPixelChannel(clut_image,channel,ClampToQuantum(clut_map[
399 ScaleQuantumToMap(GetPixelChannel(clut_image,channel,q))*
400 GetPixelChannels(clut_image)+channel]),q);
402 q+=GetPixelChannels(image);
404 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
406 if (image->progress_monitor != (MagickProgressMonitor) NULL)
411 #if defined(MAGICKCORE_OPENMP_SUPPORT)
412 #pragma omp critical (MagickCore_ClutImage)
414 proceed=SetImageProgress(image,ClutImageTag,progress++,image->rows);
415 if (proceed == MagickFalse)
419 image_view=DestroyCacheView(image_view);
420 clut_map=(double *) RelinquishMagickMemory(clut_map);
421 if ((clut_image->matte != MagickFalse) &&
422 ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0))
423 (void) SetImageAlphaChannel(image,ActivateAlphaChannel,exception);
428 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
432 % C o l o r D e c i s i o n L i s t I m a g e %
436 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
438 % ColorDecisionListImage() accepts a lightweight Color Correction Collection
439 % (CCC) file which solely contains one or more color corrections and applies
440 % the correction to the image. Here is a sample CCC file:
442 % <ColorCorrectionCollection xmlns="urn:ASC:CDL:v1.2">
443 % <ColorCorrection id="cc03345">
445 % <Slope> 0.9 1.2 0.5 </Slope>
446 % <Offset> 0.4 -0.5 0.6 </Offset>
447 % <Power> 1.0 0.8 1.5 </Power>
450 % <Saturation> 0.85 </Saturation>
453 % </ColorCorrectionCollection>
455 % which includes the slop, offset, and power for each of the RGB channels
456 % as well as the saturation.
458 % The format of the ColorDecisionListImage method is:
460 % MagickBooleanType ColorDecisionListImage(Image *image,
461 % const char *color_correction_collection,ExceptionInfo *exception)
463 % A description of each parameter follows:
465 % o image: the image.
467 % o color_correction_collection: the color correction collection in XML.
469 % o exception: return any errors or warnings in this structure.
472 MagickExport MagickBooleanType ColorDecisionListImage(Image *image,
473 const char *color_correction_collection,ExceptionInfo *exception)
475 #define ColorDecisionListCorrectImageTag "ColorDecisionList/Image"
477 typedef struct _Correction
485 typedef struct _ColorCorrection
500 token[MaxTextExtent];
531 Allocate and initialize cdl maps.
533 assert(image != (Image *) NULL);
534 assert(image->signature == MagickSignature);
535 if (image->debug != MagickFalse)
536 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
537 if (color_correction_collection == (const char *) NULL)
539 ccc=NewXMLTree((const char *) color_correction_collection,exception);
540 if (ccc == (XMLTreeInfo *) NULL)
542 cc=GetXMLTreeChild(ccc,"ColorCorrection");
543 if (cc == (XMLTreeInfo *) NULL)
545 ccc=DestroyXMLTree(ccc);
548 color_correction.red.slope=1.0;
549 color_correction.red.offset=0.0;
550 color_correction.red.power=1.0;
551 color_correction.green.slope=1.0;
552 color_correction.green.offset=0.0;
553 color_correction.green.power=1.0;
554 color_correction.blue.slope=1.0;
555 color_correction.blue.offset=0.0;
556 color_correction.blue.power=1.0;
557 color_correction.saturation=0.0;
558 sop=GetXMLTreeChild(cc,"SOPNode");
559 if (sop != (XMLTreeInfo *) NULL)
566 slope=GetXMLTreeChild(sop,"Slope");
567 if (slope != (XMLTreeInfo *) NULL)
569 content=GetXMLTreeContent(slope);
570 p=(const char *) content;
571 for (i=0; (*p != '\0') && (i < 3); i++)
573 GetMagickToken(p,&p,token);
575 GetMagickToken(p,&p,token);
580 color_correction.red.slope=StringToDouble(token,(char **) NULL);
585 color_correction.green.slope=StringToDouble(token,
591 color_correction.blue.slope=StringToDouble(token,
598 offset=GetXMLTreeChild(sop,"Offset");
599 if (offset != (XMLTreeInfo *) NULL)
601 content=GetXMLTreeContent(offset);
602 p=(const char *) content;
603 for (i=0; (*p != '\0') && (i < 3); i++)
605 GetMagickToken(p,&p,token);
607 GetMagickToken(p,&p,token);
612 color_correction.red.offset=StringToDouble(token,
618 color_correction.green.offset=StringToDouble(token,
624 color_correction.blue.offset=StringToDouble(token,
631 power=GetXMLTreeChild(sop,"Power");
632 if (power != (XMLTreeInfo *) NULL)
634 content=GetXMLTreeContent(power);
635 p=(const char *) content;
636 for (i=0; (*p != '\0') && (i < 3); i++)
638 GetMagickToken(p,&p,token);
640 GetMagickToken(p,&p,token);
645 color_correction.red.power=StringToDouble(token,(char **) NULL);
650 color_correction.green.power=StringToDouble(token,
656 color_correction.blue.power=StringToDouble(token,
664 sat=GetXMLTreeChild(cc,"SATNode");
665 if (sat != (XMLTreeInfo *) NULL)
670 saturation=GetXMLTreeChild(sat,"Saturation");
671 if (saturation != (XMLTreeInfo *) NULL)
673 content=GetXMLTreeContent(saturation);
674 p=(const char *) content;
675 GetMagickToken(p,&p,token);
676 color_correction.saturation=StringToDouble(token,
680 ccc=DestroyXMLTree(ccc);
681 if (image->debug != MagickFalse)
683 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
684 " Color Correction Collection:");
685 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
686 " color_correction.red.slope: %g",color_correction.red.slope);
687 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
688 " color_correction.red.offset: %g",color_correction.red.offset);
689 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
690 " color_correction.red.power: %g",color_correction.red.power);
691 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
692 " color_correction.green.slope: %g",color_correction.green.slope);
693 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
694 " color_correction.green.offset: %g",color_correction.green.offset);
695 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
696 " color_correction.green.power: %g",color_correction.green.power);
697 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
698 " color_correction.blue.slope: %g",color_correction.blue.slope);
699 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
700 " color_correction.blue.offset: %g",color_correction.blue.offset);
701 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
702 " color_correction.blue.power: %g",color_correction.blue.power);
703 (void) LogMagickEvent(TransformEvent,GetMagickModule(),
704 " color_correction.saturation: %g",color_correction.saturation);
706 cdl_map=(PixelInfo *) AcquireQuantumMemory(MaxMap+1UL,sizeof(*cdl_map));
707 if (cdl_map == (PixelInfo *) NULL)
708 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
710 #if defined(MAGICKCORE_OPENMP_SUPPORT)
711 #pragma omp parallel for schedule(static,4)
713 for (i=0; i <= (ssize_t) MaxMap; i++)
715 cdl_map[i].red=(MagickRealType) ScaleMapToQuantum((MagickRealType)
716 (MaxMap*(pow(color_correction.red.slope*i/MaxMap+
717 color_correction.red.offset,color_correction.red.power))));
718 cdl_map[i].green=(MagickRealType) ScaleMapToQuantum((MagickRealType)
719 (MaxMap*(pow(color_correction.green.slope*i/MaxMap+
720 color_correction.green.offset,color_correction.green.power))));
721 cdl_map[i].blue=(MagickRealType) ScaleMapToQuantum((MagickRealType)
722 (MaxMap*(pow(color_correction.blue.slope*i/MaxMap+
723 color_correction.blue.offset,color_correction.blue.power))));
725 if (image->storage_class == PseudoClass)
728 Apply transfer function to colormap.
730 #if defined(MAGICKCORE_OPENMP_SUPPORT)
731 #pragma omp parallel for schedule(static,4) shared(progress,status)
733 for (i=0; i < (ssize_t) image->colors; i++)
738 luma=0.2126*image->colormap[i].red+0.7152*image->colormap[i].green+
739 0.0722*image->colormap[i].blue;
740 image->colormap[i].red=luma+color_correction.saturation*
741 cdl_map[ScaleQuantumToMap(ClampToQuantum(
742 image->colormap[i].red))].red-luma;
743 image->colormap[i].green=luma+color_correction.saturation*
744 cdl_map[ScaleQuantumToMap(ClampToQuantum(
745 image->colormap[i].green))].green-luma;
746 image->colormap[i].blue=luma+color_correction.saturation*
747 cdl_map[ScaleQuantumToMap(ClampToQuantum(
748 image->colormap[i].blue))].blue-luma;
752 Apply transfer function to image.
756 image_view=AcquireCacheView(image);
757 #if defined(MAGICKCORE_OPENMP_SUPPORT)
758 #pragma omp parallel for schedule(static,4) shared(progress,status)
760 for (y=0; y < (ssize_t) image->rows; y++)
771 if (status == MagickFalse)
773 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
774 if (q == (Quantum *) NULL)
779 for (x=0; x < (ssize_t) image->columns; x++)
781 luma=0.2126*GetPixelRed(image,q)+0.7152*GetPixelGreen(image,q)+0.0722*
782 GetPixelBlue(image,q);
783 SetPixelRed(image,ClampToQuantum(luma+color_correction.saturation*
784 (cdl_map[ScaleQuantumToMap(GetPixelRed(image,q))].red-luma)),q);
785 SetPixelGreen(image,ClampToQuantum(luma+color_correction.saturation*
786 (cdl_map[ScaleQuantumToMap(GetPixelGreen(image,q))].green-luma)),q);
787 SetPixelBlue(image,ClampToQuantum(luma+color_correction.saturation*
788 (cdl_map[ScaleQuantumToMap(GetPixelBlue(image,q))].blue-luma)),q);
789 q+=GetPixelChannels(image);
791 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
793 if (image->progress_monitor != (MagickProgressMonitor) NULL)
798 #if defined(MAGICKCORE_OPENMP_SUPPORT)
799 #pragma omp critical (MagickCore_ColorDecisionListImageChannel)
801 proceed=SetImageProgress(image,ColorDecisionListCorrectImageTag,
802 progress++,image->rows);
803 if (proceed == MagickFalse)
807 image_view=DestroyCacheView(image_view);
808 cdl_map=(PixelInfo *) RelinquishMagickMemory(cdl_map);
813 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
817 % C o n t r a s t I m a g e %
821 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
823 % ContrastImage() enhances the intensity differences between the lighter and
824 % darker elements of the image. Set sharpen to a MagickTrue to increase the
825 % image contrast otherwise the contrast is reduced.
827 % The format of the ContrastImage method is:
829 % MagickBooleanType ContrastImage(Image *image,
830 % const MagickBooleanType sharpen,ExceptionInfo *exception)
832 % A description of each parameter follows:
834 % o image: the image.
836 % o sharpen: Increase or decrease image contrast.
838 % o exception: return any errors or warnings in this structure.
842 static void Contrast(const int sign,double *red,double *green,double *blue)
850 Enhance contrast: dark color become darker, light color become lighter.
852 assert(red != (double *) NULL);
853 assert(green != (double *) NULL);
854 assert(blue != (double *) NULL);
858 ConvertRGBToHSB(*red,*green,*blue,&hue,&saturation,&brightness);
859 brightness+=0.5*sign*(0.5*(sin((double) (MagickPI*(brightness-0.5)))+1.0)-
861 if (brightness > 1.0)
864 if (brightness < 0.0)
866 ConvertHSBToRGB(hue,saturation,brightness,red,green,blue);
869 MagickExport MagickBooleanType ContrastImage(Image *image,
870 const MagickBooleanType sharpen,ExceptionInfo *exception)
872 #define ContrastImageTag "Contrast/Image"
892 assert(image != (Image *) NULL);
893 assert(image->signature == MagickSignature);
894 if (image->debug != MagickFalse)
895 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
896 sign=sharpen != MagickFalse ? 1 : -1;
897 if (image->storage_class == PseudoClass)
900 Contrast enhance colormap.
902 for (i=0; i < (ssize_t) image->colors; i++)
903 Contrast(sign,&image->colormap[i].red,&image->colormap[i].green,
904 &image->colormap[i].blue);
907 Contrast enhance image.
911 image_view=AcquireCacheView(image);
912 #if defined(MAGICKCORE_OPENMP_SUPPORT)
913 #pragma omp parallel for schedule(static,4) shared(progress,status)
915 for (y=0; y < (ssize_t) image->rows; y++)
928 if (status == MagickFalse)
930 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
931 if (q == (Quantum *) NULL)
936 for (x=0; x < (ssize_t) image->columns; x++)
938 red=(double) GetPixelRed(image,q);
939 green=(double) GetPixelGreen(image,q);
940 blue=(double) GetPixelBlue(image,q);
941 Contrast(sign,&red,&green,&blue);
942 SetPixelRed(image,ClampToQuantum(red),q);
943 SetPixelGreen(image,ClampToQuantum(green),q);
944 SetPixelBlue(image,ClampToQuantum(blue),q);
945 q+=GetPixelChannels(image);
947 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
949 if (image->progress_monitor != (MagickProgressMonitor) NULL)
954 #if defined(MAGICKCORE_OPENMP_SUPPORT)
955 #pragma omp critical (MagickCore_ContrastImage)
957 proceed=SetImageProgress(image,ContrastImageTag,progress++,image->rows);
958 if (proceed == MagickFalse)
962 image_view=DestroyCacheView(image_view);
967 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
971 % C o n t r a s t S t r e t c h I m a g e %
975 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
977 % ContrastStretchImage() is a simple image enhancement technique that attempts
978 % to improve the contrast in an image by `stretching' the range of intensity
979 % values it contains to span a desired range of values. It differs from the
980 % more sophisticated histogram equalization in that it can only apply a
981 % linear scaling function to the image pixel values. As a result the
982 % `enhancement' is less harsh.
984 % The format of the ContrastStretchImage method is:
986 % MagickBooleanType ContrastStretchImage(Image *image,
987 % const char *levels,ExceptionInfo *exception)
989 % A description of each parameter follows:
991 % o image: the image.
993 % o black_point: the black point.
995 % o white_point: the white point.
997 % o levels: Specify the levels where the black and white points have the
998 % range of 0 to number-of-pixels (e.g. 1%, 10x90%, etc.).
1000 % o exception: return any errors or warnings in this structure.
1003 MagickExport MagickBooleanType ContrastStretchImage(Image *image,
1004 const double black_point,const double white_point,ExceptionInfo *exception)
1006 #define MaxRange(color) ((MagickRealType) ScaleQuantumToMap((Quantum) (color)))
1007 #define ContrastStretchImageTag "ContrastStretch/Image"
1031 Allocate histogram and stretch map.
1033 assert(image != (Image *) NULL);
1034 assert(image->signature == MagickSignature);
1035 if (image->debug != MagickFalse)
1036 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1037 black=(double *) AcquireQuantumMemory(GetPixelChannels(image),sizeof(*black));
1038 white=(double *) AcquireQuantumMemory(GetPixelChannels(image),sizeof(*white));
1039 histogram=(double *) AcquireQuantumMemory(MaxMap+1UL,
1040 GetPixelChannels(image)*sizeof(*histogram));
1041 stretch_map=(double *) AcquireQuantumMemory(MaxMap+1UL,
1042 GetPixelChannels(image)*sizeof(*stretch_map));
1043 if ((black == (double *) NULL) || (white == (double *) NULL) ||
1044 (histogram == (double *) NULL) || (stretch_map == (double *) NULL))
1046 if (stretch_map != (double *) NULL)
1047 stretch_map=(double *) RelinquishMagickMemory(stretch_map);
1048 if (histogram != (double *) NULL)
1049 histogram=(double *) RelinquishMagickMemory(histogram);
1050 if (white != (double *) NULL)
1051 white=(double *) RelinquishMagickMemory(white);
1052 if (black != (double *) NULL)
1053 black=(double *) RelinquishMagickMemory(black);
1054 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
1061 (void) ResetMagickMemory(histogram,0,(MaxMap+1)*GetPixelChannels(image)*
1062 sizeof(*histogram));
1063 image_view=AcquireCacheView(image);
1064 for (y=0; y < (ssize_t) image->rows; y++)
1066 register const Quantum
1072 if (status == MagickFalse)
1074 p=GetCacheViewVirtualPixels(image_view,0,y,image->columns,1,exception);
1075 if (p == (const Quantum *) NULL)
1080 for (x=0; x < (ssize_t) image->columns; x++)
1085 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1086 histogram[GetPixelChannels(image)*ScaleQuantumToMap(p[i])+i]++;
1087 p+=GetPixelChannels(image);
1091 Find the histogram boundaries by locating the black/white levels.
1093 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1094 #pragma omp parallel for schedule(static,4) shared(progress,status)
1096 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1105 white[i]=MaxRange(QuantumRange);
1107 for (j=0; j <= (ssize_t) MaxMap; j++)
1109 intensity+=histogram[GetPixelChannels(image)*j+i];
1110 if (intensity > black_point)
1113 black[i]=(MagickRealType) j;
1115 for (j=(ssize_t) MaxMap; j != 0; j--)
1117 intensity+=histogram[GetPixelChannels(image)*j+i];
1118 if (intensity > ((double) image->columns*image->rows-white_point))
1121 white[i]=(MagickRealType) j;
1123 histogram=(double *) RelinquishMagickMemory(histogram);
1125 Stretch the histogram to create the stretched image mapping.
1127 (void) ResetMagickMemory(stretch_map,0,(MaxMap+1)*GetPixelChannels(image)*
1128 sizeof(*stretch_map));
1129 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1130 #pragma omp parallel for schedule(static,4) shared(progress,status)
1132 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1137 for (j=0; j <= (ssize_t) MaxMap; j++)
1139 if (j < (ssize_t) black[i])
1140 stretch_map[GetPixelChannels(image)*j+i]=0.0;
1142 if (j > (ssize_t) white[i])
1143 stretch_map[GetPixelChannels(image)*j+i]=(MagickRealType)
1146 if (black[i] != white[i])
1147 stretch_map[GetPixelChannels(image)*j+i]=(MagickRealType)
1148 ScaleMapToQuantum((MagickRealType) (MaxMap*(j-black[i])/
1149 (white[i]-black[i])));
1152 if (image->storage_class == PseudoClass)
1158 Stretch-contrast colormap.
1160 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1161 #pragma omp parallel for schedule(static,4) shared(progress,status)
1163 for (j=0; j < (ssize_t) image->colors; j++)
1165 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
1167 i=GetPixelChannelMapChannel(image,RedPixelChannel);
1168 if (black[i] != white[i])
1169 image->colormap[j].red=stretch_map[GetPixelChannels(image)*
1170 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].red))]+i;
1172 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
1174 i=GetPixelChannelMapChannel(image,GreenPixelChannel);
1175 if (black[i] != white[i])
1176 image->colormap[j].green=stretch_map[GetPixelChannels(image)*
1177 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].green))]+i;
1179 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
1181 i=GetPixelChannelMapChannel(image,BluePixelChannel);
1182 if (black[i] != white[i])
1183 image->colormap[j].blue=stretch_map[GetPixelChannels(image)*
1184 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].blue))]+i;
1186 if ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0)
1188 i=GetPixelChannelMapChannel(image,AlphaPixelChannel);
1189 if (black[i] != white[i])
1190 image->colormap[j].alpha=stretch_map[GetPixelChannels(image)*
1191 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].alpha))]+i;
1196 Stretch-contrast image.
1200 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1201 #pragma omp parallel for schedule(static,4) shared(progress,status)
1203 for (y=0; y < (ssize_t) image->rows; y++)
1211 if (status == MagickFalse)
1213 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
1214 if (q == (Quantum *) NULL)
1219 for (x=0; x < (ssize_t) image->columns; x++)
1224 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1232 channel=GetPixelChannelMapChannel(image,i);
1233 traits=GetPixelChannelMapTraits(image,channel);
1234 if (((traits & UpdatePixelTrait) != 0) && (black[i] != white[i]))
1235 q[i]=ClampToQuantum(stretch_map[GetPixelChannels(image)*
1236 ScaleQuantumToMap(q[i])+i]);
1238 q+=GetPixelChannels(image);
1240 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1242 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1247 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1248 #pragma omp critical (MagickCore_ContrastStretchImage)
1250 proceed=SetImageProgress(image,ContrastStretchImageTag,progress++,
1252 if (proceed == MagickFalse)
1256 image_view=DestroyCacheView(image_view);
1257 stretch_map=(double *) RelinquishMagickMemory(stretch_map);
1258 white=(double *) RelinquishMagickMemory(white);
1259 black=(double *) RelinquishMagickMemory(black);
1264 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1268 % E n h a n c e I m a g e %
1272 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1274 % EnhanceImage() applies a digital filter that improves the quality of a
1277 % The format of the EnhanceImage method is:
1279 % Image *EnhanceImage(const Image *image,ExceptionInfo *exception)
1281 % A description of each parameter follows:
1283 % o image: the image.
1285 % o exception: return any errors or warnings in this structure.
1288 MagickExport Image *EnhanceImage(const Image *image,ExceptionInfo *exception)
1290 #define EnhancePixel(weight) \
1291 mean=((MagickRealType) r[i]+GetPixelChannel(enhance_image,channel,q))/2.0; \
1292 distance=(MagickRealType) r[i]-(MagickRealType) GetPixelChannel( \
1293 enhance_image,channel,q); \
1294 distance_squared=QuantumScale*(2.0*((MagickRealType) QuantumRange+1.0)+ \
1295 mean)*distance*distance; \
1296 if (distance_squared < ((MagickRealType) QuantumRange*(MagickRealType) \
1297 QuantumRange/25.0f)) \
1299 aggregate+=(weight)*r[i]; \
1300 total_weight+=(weight); \
1302 r+=GetPixelChannels(image);
1303 #define EnhanceImageTag "Enhance/Image"
1322 Initialize enhanced image attributes.
1324 assert(image != (const Image *) NULL);
1325 assert(image->signature == MagickSignature);
1326 if (image->debug != MagickFalse)
1327 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1328 assert(exception != (ExceptionInfo *) NULL);
1329 assert(exception->signature == MagickSignature);
1330 enhance_image=CloneImage(image,image->columns,image->rows,MagickTrue,
1332 if (enhance_image == (Image *) NULL)
1333 return((Image *) NULL);
1334 if (SetImageStorageClass(enhance_image,DirectClass,exception) == MagickFalse)
1336 enhance_image=DestroyImage(enhance_image);
1337 return((Image *) NULL);
1344 image_view=AcquireCacheView(image);
1345 enhance_view=AcquireCacheView(enhance_image);
1346 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1347 #pragma omp parallel for schedule(static,4) shared(progress,status)
1349 for (y=0; y < (ssize_t) image->rows; y++)
1351 register const Quantum
1363 if (status == MagickFalse)
1365 p=GetCacheViewVirtualPixels(image_view,-2,y-2,image->columns+4,5,exception);
1366 q=QueueCacheViewAuthenticPixels(enhance_view,0,y,enhance_image->columns,1,
1368 if ((p == (const Quantum *) NULL) || (q == (Quantum *) NULL))
1373 center=(ssize_t) GetPixelChannels(image)*(2*(image->columns+4)+2);
1374 for (x=0; x < (ssize_t) image->columns; x++)
1379 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1395 register const Quantum
1398 channel=GetPixelChannelMapChannel(image,i);
1399 traits=GetPixelChannelMapTraits(image,channel);
1400 enhance_traits=GetPixelChannelMapTraits(enhance_image,channel);
1401 if ((traits == UndefinedPixelTrait) ||
1402 (enhance_traits == UndefinedPixelTrait))
1404 SetPixelChannel(enhance_image,channel,p[center+i],q);
1405 if ((enhance_traits & CopyPixelTrait) != 0)
1408 Compute weighted average of target pixel color components.
1413 EnhancePixel(5.0); EnhancePixel(8.0); EnhancePixel(10.0);
1414 EnhancePixel(8.0); EnhancePixel(5.0);
1415 r=p+1*GetPixelChannels(image)*(image->columns+4);
1416 EnhancePixel(8.0); EnhancePixel(20.0); EnhancePixel(40.0);
1417 EnhancePixel(20.0); EnhancePixel(8.0);
1418 r=p+2*GetPixelChannels(image)*(image->columns+4);
1419 EnhancePixel(10.0); EnhancePixel(40.0); EnhancePixel(80.0);
1420 EnhancePixel(40.0); EnhancePixel(10.0);
1421 r=p+3*GetPixelChannels(image)*(image->columns+4);
1422 EnhancePixel(8.0); EnhancePixel(20.0); EnhancePixel(40.0);
1423 EnhancePixel(20.0); EnhancePixel(8.0);
1424 r=p+4*GetPixelChannels(image)*(image->columns+4);
1425 EnhancePixel(5.0); EnhancePixel(8.0); EnhancePixel(10.0);
1426 EnhancePixel(8.0); EnhancePixel(5.0);
1427 SetPixelChannel(enhance_image,channel,ClampToQuantum(aggregate/
1430 p+=GetPixelChannels(image);
1431 q+=GetPixelChannels(enhance_image);
1433 if (SyncCacheViewAuthenticPixels(enhance_view,exception) == MagickFalse)
1435 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1440 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1441 #pragma omp critical (MagickCore_EnhanceImage)
1443 proceed=SetImageProgress(image,EnhanceImageTag,progress++,image->rows);
1444 if (proceed == MagickFalse)
1448 enhance_view=DestroyCacheView(enhance_view);
1449 image_view=DestroyCacheView(image_view);
1450 return(enhance_image);
1454 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1458 % E q u a l i z e I m a g e %
1462 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1464 % EqualizeImage() applies a histogram equalization to the image.
1466 % The format of the EqualizeImage method is:
1468 % MagickBooleanType EqualizeImage(Image *image,ExceptionInfo *exception)
1470 % A description of each parameter follows:
1472 % o image: the image.
1474 % o exception: return any errors or warnings in this structure.
1477 MagickExport MagickBooleanType EqualizeImage(Image *image,
1478 ExceptionInfo *exception)
1480 #define EqualizeImageTag "Equalize/Image"
1492 black[CompositePixelChannel],
1496 white[CompositePixelChannel];
1505 Allocate and initialize histogram arrays.
1507 assert(image != (Image *) NULL);
1508 assert(image->signature == MagickSignature);
1509 if (image->debug != MagickFalse)
1510 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1511 equalize_map=(MagickRealType *) AcquireQuantumMemory(MaxMap+1UL,
1512 GetPixelChannels(image)*sizeof(*equalize_map));
1513 histogram=(MagickRealType *) AcquireQuantumMemory(MaxMap+1UL,
1514 GetPixelChannels(image)*sizeof(*histogram));
1515 map=(MagickRealType *) AcquireQuantumMemory(MaxMap+1UL,
1516 GetPixelChannels(image)*sizeof(*map));
1517 if ((equalize_map == (MagickRealType *) NULL) ||
1518 (histogram == (MagickRealType *) NULL) ||
1519 (map == (MagickRealType *) NULL))
1521 if (map != (MagickRealType *) NULL)
1522 map=(MagickRealType *) RelinquishMagickMemory(map);
1523 if (histogram != (MagickRealType *) NULL)
1524 histogram=(MagickRealType *) RelinquishMagickMemory(histogram);
1525 if (equalize_map != (MagickRealType *) NULL)
1526 equalize_map=(MagickRealType *) RelinquishMagickMemory(equalize_map);
1527 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
1534 (void) ResetMagickMemory(histogram,0,(MaxMap+1)*GetPixelChannels(image)*
1535 sizeof(*histogram));
1536 image_view=AcquireCacheView(image);
1537 for (y=0; y < (ssize_t) image->rows; y++)
1539 register const Quantum
1545 if (status == MagickFalse)
1547 p=GetCacheViewVirtualPixels(image_view,0,y,image->columns,1,exception);
1548 if (p == (const Quantum *) NULL)
1553 for (x=0; x < (ssize_t) image->columns; x++)
1558 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1559 histogram[GetPixelChannels(image)*ScaleQuantumToMap(p[i])+i]++;
1560 p+=GetPixelChannels(image);
1564 Integrate the histogram to get the equalization map.
1566 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1567 #pragma omp parallel for schedule(static,4) shared(progress,status)
1569 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1578 for (j=0; j <= (ssize_t) MaxMap; j++)
1580 intensity+=histogram[GetPixelChannels(image)*j+i];
1581 map[GetPixelChannels(image)*j+i]=intensity;
1584 (void) ResetMagickMemory(equalize_map,0,(MaxMap+1)*GetPixelChannels(image)*
1585 sizeof(*equalize_map));
1586 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1587 #pragma omp parallel for schedule(static,4) shared(progress,status)
1589 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1595 white[i]=map[GetPixelChannels(image)*MaxMap+i];
1596 if (black[i] != white[i])
1597 for (j=0; j <= (ssize_t) MaxMap; j++)
1598 equalize_map[GetPixelChannels(image)*j+i]=(MagickRealType)
1599 ScaleMapToQuantum((MagickRealType) ((MaxMap*(map[
1600 GetPixelChannels(image)*j+i]-black[i]))/(white[i]-black[i])));
1602 histogram=(MagickRealType *) RelinquishMagickMemory(histogram);
1603 map=(MagickRealType *) RelinquishMagickMemory(map);
1604 if (image->storage_class == PseudoClass)
1615 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1616 #pragma omp parallel for schedule(static,4) shared(progress,status)
1618 for (j=0; j < (ssize_t) image->colors; j++)
1620 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
1622 channel=GetPixelChannelMapChannel(image,RedPixelChannel);
1623 if (black[channel] != white[channel])
1624 image->colormap[j].red=equalize_map[GetPixelChannels(image)*
1625 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].red))]+
1628 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
1630 channel=GetPixelChannelMapChannel(image,GreenPixelChannel);
1631 if (black[channel] != white[channel])
1632 image->colormap[j].green=equalize_map[GetPixelChannels(image)*
1633 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].green))]+
1636 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
1638 channel=GetPixelChannelMapChannel(image,BluePixelChannel);
1639 if (black[channel] != white[channel])
1640 image->colormap[j].blue=equalize_map[GetPixelChannels(image)*
1641 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].blue))]+
1644 if ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0)
1646 channel=GetPixelChannelMapChannel(image,AlphaPixelChannel);
1647 if (black[channel] != white[channel])
1648 image->colormap[j].alpha=equalize_map[GetPixelChannels(image)*
1649 ScaleQuantumToMap(ClampToQuantum(image->colormap[j].alpha))]+
1658 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1659 #pragma omp parallel for schedule(static,4) shared(progress,status)
1661 for (y=0; y < (ssize_t) image->rows; y++)
1669 if (status == MagickFalse)
1671 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
1672 if (q == (Quantum *) NULL)
1677 for (x=0; x < (ssize_t) image->columns; x++)
1682 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1690 channel=GetPixelChannelMapChannel(image,i);
1691 traits=GetPixelChannelMapTraits(image,channel);
1692 if (((traits & UpdatePixelTrait) != 0) && (black[i] != white[i]))
1693 q[i]=ClampToQuantum(equalize_map[GetPixelChannels(image)*
1694 ScaleQuantumToMap(q[i])+i]);
1696 q+=GetPixelChannels(image);
1698 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1700 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1705 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1706 #pragma omp critical (MagickCore_EqualizeImage)
1708 proceed=SetImageProgress(image,EqualizeImageTag,progress++,image->rows);
1709 if (proceed == MagickFalse)
1713 image_view=DestroyCacheView(image_view);
1714 equalize_map=(MagickRealType *) RelinquishMagickMemory(equalize_map);
1719 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1723 % G a m m a I m a g e %
1727 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1729 % GammaImage() gamma-corrects a particular image channel. The same
1730 % image viewed on different devices will have perceptual differences in the
1731 % way the image's intensities are represented on the screen. Specify
1732 % individual gamma levels for the red, green, and blue channels, or adjust
1733 % all three with the gamma parameter. Values typically range from 0.8 to 2.3.
1735 % You can also reduce the influence of a particular channel with a gamma
1738 % The format of the GammaImage method is:
1740 % MagickBooleanType GammaImage(Image *image,const double gamma,
1741 % ExceptionInfo *exception)
1743 % A description of each parameter follows:
1745 % o image: the image.
1747 % o level: the image gamma as a string (e.g. 1.6,1.2,1.0).
1749 % o gamma: the image gamma.
1752 MagickExport MagickBooleanType GammaImage(Image *image,const double gamma,
1753 ExceptionInfo *exception)
1755 #define GammaCorrectImageTag "GammaCorrect/Image"
1776 Allocate and initialize gamma maps.
1778 assert(image != (Image *) NULL);
1779 assert(image->signature == MagickSignature);
1780 if (image->debug != MagickFalse)
1781 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1784 gamma_map=(double *) AcquireQuantumMemory(MaxMap+1UL,sizeof(*gamma_map));
1785 if (gamma_map == (double *) NULL)
1786 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
1788 (void) ResetMagickMemory(gamma_map,0,(MaxMap+1)*sizeof(*gamma_map));
1790 #if defined(MAGICKCORE_OPENMP_SUPPORT) && (MaxMap > 256)
1791 #pragma omp parallel for
1793 for (i=0; i <= (ssize_t) MaxMap; i++)
1794 gamma_map[i]=(MagickRealType) ScaleMapToQuantum((
1795 MagickRealType) (MaxMap*pow((double) i/MaxMap,1.0/gamma)));
1796 if (image->storage_class == PseudoClass)
1799 Gamma-correct colormap.
1801 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1802 #pragma omp parallel for schedule(static,4) shared(progress,status)
1804 for (i=0; i < (ssize_t) image->colors; i++)
1806 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
1807 image->colormap[i].red=gamma_map[
1808 ScaleQuantumToMap(ClampToQuantum(image->colormap[i].red))];
1809 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
1810 image->colormap[i].green=gamma_map[
1811 ScaleQuantumToMap(ClampToQuantum(image->colormap[i].green))];
1812 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
1813 image->colormap[i].blue=gamma_map[
1814 ScaleQuantumToMap(ClampToQuantum(image->colormap[i].blue))];
1815 if ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0)
1816 image->colormap[i].alpha=gamma_map[
1817 ScaleQuantumToMap(ClampToQuantum(image->colormap[i].alpha))];
1821 Gamma-correct image.
1825 image_view=AcquireCacheView(image);
1826 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1827 #pragma omp parallel for schedule(static,4) shared(progress,status)
1829 for (y=0; y < (ssize_t) image->rows; y++)
1837 if (status == MagickFalse)
1839 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
1840 if (q == (Quantum *) NULL)
1845 for (x=0; x < (ssize_t) image->columns; x++)
1850 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
1858 channel=GetPixelChannelMapChannel(image,i);
1859 traits=GetPixelChannelMapTraits(image,channel);
1860 if ((traits & UpdatePixelTrait) != 0)
1861 q[i]=ClampToQuantum(gamma_map[ScaleQuantumToMap(q[i])]);
1863 q+=GetPixelChannels(image);
1865 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
1867 if (image->progress_monitor != (MagickProgressMonitor) NULL)
1872 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1873 #pragma omp critical (MagickCore_GammaImage)
1875 proceed=SetImageProgress(image,GammaCorrectImageTag,progress++,
1877 if (proceed == MagickFalse)
1881 image_view=DestroyCacheView(image_view);
1882 gamma_map=(double *) RelinquishMagickMemory(gamma_map);
1883 if (image->gamma != 0.0)
1884 image->gamma*=gamma;
1889 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1893 % H a l d C l u t I m a g e %
1897 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1899 % HaldClutImage() applies a Hald color lookup table to the image. A Hald
1900 % color lookup table is a 3-dimensional color cube mapped to 2 dimensions.
1901 % Create it with the HALD coder. You can apply any color transformation to
1902 % the Hald image and then use this method to apply the transform to the
1905 % The format of the HaldClutImage method is:
1907 % MagickBooleanType HaldClutImage(Image *image,Image *hald_image,
1908 % ExceptionInfo *exception)
1910 % A description of each parameter follows:
1912 % o image: the image, which is replaced by indexed CLUT values
1914 % o hald_image: the color lookup table image for replacement color values.
1916 % o exception: return any errors or warnings in this structure.
1920 static inline size_t MagickMin(const size_t x,const size_t y)
1927 MagickExport MagickBooleanType HaldClutImage(Image *image,
1928 const Image *hald_image,ExceptionInfo *exception)
1930 #define HaldClutImageTag "Clut/Image"
1932 typedef struct _HaldInfo
1964 assert(image != (Image *) NULL);
1965 assert(image->signature == MagickSignature);
1966 if (image->debug != MagickFalse)
1967 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
1968 assert(hald_image != (Image *) NULL);
1969 assert(hald_image->signature == MagickSignature);
1970 if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
1971 return(MagickFalse);
1972 if (image->matte == MagickFalse)
1973 (void) SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
1979 length=MagickMin(hald_image->columns,hald_image->rows);
1980 for (level=2; (level*level*level) < length; level++) ;
1982 cube_size=level*level;
1983 width=(double) hald_image->columns;
1984 GetPixelInfo(hald_image,&zero);
1985 image_view=AcquireCacheView(image);
1986 hald_view=AcquireCacheView(hald_image);
1987 #if defined(MAGICKCORE_OPENMP_SUPPORT)
1988 #pragma omp parallel for schedule(static,4) shared(progress,status)
1990 for (y=0; y < (ssize_t) image->rows; y++)
1998 if (status == MagickFalse)
2000 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
2001 if (q == (Quantum *) NULL)
2006 for (x=0; x < (ssize_t) image->columns; x++)
2021 point.x=QuantumScale*(level-1.0)*GetPixelRed(image,q);
2022 point.y=QuantumScale*(level-1.0)*GetPixelGreen(image,q);
2023 point.z=QuantumScale*(level-1.0)*GetPixelBlue(image,q);
2024 offset=point.x+level*floor(point.y)+cube_size*floor(point.z);
2025 point.x-=floor(point.x);
2026 point.y-=floor(point.y);
2027 point.z-=floor(point.z);
2029 (void) InterpolatePixelInfo(image,hald_view,image->interpolate,
2030 fmod(offset,width),floor(offset/width),&pixel1,exception);
2032 (void) InterpolatePixelInfo(image,hald_view,image->interpolate,
2033 fmod(offset+level,width),floor((offset+level)/width),&pixel2,exception);
2035 CompositePixelInfoAreaBlend(&pixel1,pixel1.alpha,&pixel2,pixel2.alpha,
2038 (void) InterpolatePixelInfo(image,hald_view,image->interpolate,
2039 fmod(offset,width),floor(offset/width),&pixel1,exception);
2040 (void) InterpolatePixelInfo(image,hald_view,image->interpolate,
2041 fmod(offset+level,width),floor((offset+level)/width),&pixel2,exception);
2043 CompositePixelInfoAreaBlend(&pixel1,pixel1.alpha,&pixel2,pixel2.alpha,
2046 CompositePixelInfoAreaBlend(&pixel3,pixel3.alpha,&pixel4,pixel4.alpha,
2048 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
2049 SetPixelRed(image,ClampToQuantum(pixel.red),q);
2050 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
2051 SetPixelGreen(image,ClampToQuantum(pixel.green),q);
2052 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
2053 SetPixelBlue(image,ClampToQuantum(pixel.blue),q);
2054 if (((GetPixelBlackTraits(image) & UpdatePixelTrait) != 0) &&
2055 (image->colorspace == CMYKColorspace))
2056 SetPixelBlack(image,ClampToQuantum(pixel.black),q);
2057 if (((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0) &&
2058 (image->matte != MagickFalse))
2059 SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
2060 q+=GetPixelChannels(image);
2062 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
2064 if (image->progress_monitor != (MagickProgressMonitor) NULL)
2069 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2070 #pragma omp critical (MagickCore_HaldClutImage)
2072 proceed=SetImageProgress(image,HaldClutImageTag,progress++,image->rows);
2073 if (proceed == MagickFalse)
2077 hald_view=DestroyCacheView(hald_view);
2078 image_view=DestroyCacheView(image_view);
2083 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2087 % L e v e l I m a g e %
2091 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2093 % LevelImage() adjusts the levels of a particular image channel by
2094 % scaling the colors falling between specified white and black points to
2095 % the full available quantum range.
2097 % The parameters provided represent the black, and white points. The black
2098 % point specifies the darkest color in the image. Colors darker than the
2099 % black point are set to zero. White point specifies the lightest color in
2100 % the image. Colors brighter than the white point are set to the maximum
2103 % If a '!' flag is given, map black and white colors to the given levels
2104 % rather than mapping those levels to black and white. See
2105 % LevelizeImage() below.
2107 % Gamma specifies a gamma correction to apply to the image.
2109 % The format of the LevelImage method is:
2111 % MagickBooleanType LevelImage(Image *image,const double black_point,
2112 % const double white_point,const double gamma,ExceptionInfo *exception)
2114 % A description of each parameter follows:
2116 % o image: the image.
2118 % o black_point: The level to map zero (black) to.
2120 % o white_point: The level to map QuantumRange (white) to.
2122 % o exception: return any errors or warnings in this structure.
2126 static inline MagickRealType LevelPixel(const double black_point,
2127 const double white_point,const double gamma,const MagickRealType pixel)
2133 if (pixel < black_point)
2135 if (pixel > white_point)
2136 return((MagickRealType) QuantumRange);
2137 scale=(white_point != black_point) ? 1.0/(white_point-black_point) : 1.0;
2138 level_pixel=(MagickRealType) QuantumRange*pow(scale*((double) pixel-
2139 black_point),1.0/gamma);
2140 return(level_pixel);
2143 MagickExport MagickBooleanType LevelImage(Image *image,const double black_point,
2144 const double white_point,const double gamma,ExceptionInfo *exception)
2146 #define LevelImageTag "Level/Image"
2164 Allocate and initialize levels map.
2166 assert(image != (Image *) NULL);
2167 assert(image->signature == MagickSignature);
2168 if (image->debug != MagickFalse)
2169 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
2170 if (image->storage_class == PseudoClass)
2171 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2172 #pragma omp parallel for schedule(static,4) shared(progress,status)
2174 for (i=0; i < (ssize_t) image->colors; i++)
2179 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
2180 image->colormap[i].red=(double) ClampToQuantum(LevelPixel(black_point,
2181 white_point,gamma,image->colormap[i].red));
2182 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
2183 image->colormap[i].green=(double) ClampToQuantum(LevelPixel(black_point,
2184 white_point,gamma,image->colormap[i].green));
2185 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
2186 image->colormap[i].blue=(double) ClampToQuantum(LevelPixel(black_point,
2187 white_point,gamma,image->colormap[i].blue));
2188 if ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0)
2189 image->colormap[i].alpha=(double) ClampToQuantum(LevelPixel(black_point,
2190 white_point,gamma,image->colormap[i].alpha));
2197 image_view=AcquireCacheView(image);
2198 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2199 #pragma omp parallel for schedule(static,4) shared(progress,status)
2201 for (y=0; y < (ssize_t) image->rows; y++)
2209 if (status == MagickFalse)
2211 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
2212 if (q == (Quantum *) NULL)
2217 for (x=0; x < (ssize_t) image->columns; x++)
2222 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
2230 channel=GetPixelChannelMapChannel(image,i);
2231 traits=GetPixelChannelMapTraits(image,channel);
2232 if (traits == UndefinedPixelTrait)
2234 q[i]=ClampToQuantum(LevelPixel(black_point,white_point,gamma,
2235 (MagickRealType) q[i]));
2237 q+=GetPixelChannels(image);
2239 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
2241 if (image->progress_monitor != (MagickProgressMonitor) NULL)
2246 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2247 #pragma omp critical (MagickCore_LevelImage)
2249 proceed=SetImageProgress(image,LevelImageTag,progress++,image->rows);
2250 if (proceed == MagickFalse)
2254 image_view=DestroyCacheView(image_view);
2259 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2263 % L e v e l i z e I m a g e %
2267 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2269 % LevelizeImage() applies the reversed LevelImage() operation to just
2270 % the specific channels specified. It compresses the full range of color
2271 % values, so that they lie between the given black and white points. Gamma is
2272 % applied before the values are mapped.
2274 % LevelizeImage() can be called with by using a +level command line
2275 % API option, or using a '!' on a -level or LevelImage() geometry string.
2277 % It can be used for example de-contrast a greyscale image to the exact
2278 % levels specified. Or by using specific levels for each channel of an image
2279 % you can convert a gray-scale image to any linear color gradient, according
2282 % The format of the LevelizeImage method is:
2284 % MagickBooleanType LevelizeImage(Image *image,const double black_point,
2285 % const double white_point,const double gamma,ExceptionInfo *exception)
2287 % A description of each parameter follows:
2289 % o image: the image.
2291 % o black_point: The level to map zero (black) to.
2293 % o white_point: The level to map QuantumRange (white) to.
2295 % o gamma: adjust gamma by this factor before mapping values.
2297 % o exception: return any errors or warnings in this structure.
2300 MagickExport MagickBooleanType LevelizeImage(Image *image,
2301 const double black_point,const double white_point,const double gamma,
2302 ExceptionInfo *exception)
2304 #define LevelizeImageTag "Levelize/Image"
2305 #define LevelizeValue(x) (ClampToQuantum(((MagickRealType) \
2306 pow((double) (QuantumScale*(x)),1.0/gamma))*(white_point-black_point)+ \
2325 Allocate and initialize levels map.
2327 assert(image != (Image *) NULL);
2328 assert(image->signature == MagickSignature);
2329 if (image->debug != MagickFalse)
2330 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
2331 if (image->storage_class == PseudoClass)
2332 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2333 #pragma omp parallel for schedule(static,4) shared(progress,status)
2335 for (i=0; i < (ssize_t) image->colors; i++)
2340 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
2341 image->colormap[i].red=(double) LevelizeValue(
2342 image->colormap[i].red);
2343 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
2344 image->colormap[i].green=(double) LevelizeValue(
2345 image->colormap[i].green);
2346 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
2347 image->colormap[i].blue=(double) LevelizeValue(
2348 image->colormap[i].blue);
2349 if ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0)
2350 image->colormap[i].alpha=(double) LevelizeValue(
2351 image->colormap[i].alpha);
2358 image_view=AcquireCacheView(image);
2359 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2360 #pragma omp parallel for schedule(static,4) shared(progress,status)
2362 for (y=0; y < (ssize_t) image->rows; y++)
2370 if (status == MagickFalse)
2372 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
2373 if (q == (Quantum *) NULL)
2378 for (x=0; x < (ssize_t) image->columns; x++)
2383 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
2391 channel=GetPixelChannelMapChannel(image,i);
2392 traits=GetPixelChannelMapTraits(image,channel);
2393 if (traits == UndefinedPixelTrait)
2395 q[i]=LevelizeValue(q[i]);
2397 q+=GetPixelChannels(image);
2399 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
2401 if (image->progress_monitor != (MagickProgressMonitor) NULL)
2406 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2407 #pragma omp critical (MagickCore_LevelizeImage)
2409 proceed=SetImageProgress(image,LevelizeImageTag,progress++,image->rows);
2410 if (proceed == MagickFalse)
2414 image_view=DestroyCacheView(image_view);
2419 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2423 % L e v e l I m a g e C o l o r s %
2427 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2429 % LevelImageColors() maps the given color to "black" and "white" values,
2430 % linearly spreading out the colors, and level values on a channel by channel
2431 % bases, as per LevelImage(). The given colors allows you to specify
2432 % different level ranges for each of the color channels separately.
2434 % If the boolean 'invert' is set true the image values will modifyed in the
2435 % reverse direction. That is any existing "black" and "white" colors in the
2436 % image will become the color values given, with all other values compressed
2437 % appropriatally. This effectivally maps a greyscale gradient into the given
2440 % The format of the LevelImageColors method is:
2442 % MagickBooleanType LevelImageColors(Image *image,
2443 % const PixelInfo *black_color,const PixelInfo *white_color,
2444 % const MagickBooleanType invert,ExceptionInfo *exception)
2446 % A description of each parameter follows:
2448 % o image: the image.
2450 % o black_color: The color to map black to/from
2452 % o white_point: The color to map white to/from
2454 % o invert: if true map the colors (levelize), rather than from (level)
2456 % o exception: return any errors or warnings in this structure.
2459 MagickExport MagickBooleanType LevelImageColors(Image *image,
2460 const PixelInfo *black_color,const PixelInfo *white_color,
2461 const MagickBooleanType invert,ExceptionInfo *exception)
2470 Allocate and initialize levels map.
2472 assert(image != (Image *) NULL);
2473 assert(image->signature == MagickSignature);
2474 if (image->debug != MagickFalse)
2475 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
2477 if (invert == MagickFalse)
2479 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
2481 channel_mask=SetPixelChannelMask(image,RedChannel);
2482 status|=LevelImage(image,black_color->red,white_color->red,1.0,
2484 (void) SetPixelChannelMask(image,channel_mask);
2486 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
2488 channel_mask=SetPixelChannelMask(image,GreenChannel);
2489 status|=LevelImage(image,black_color->green,white_color->green,1.0,
2491 (void) SetPixelChannelMask(image,channel_mask);
2493 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
2495 channel_mask=SetPixelChannelMask(image,BlueChannel);
2496 status|=LevelImage(image,black_color->blue,white_color->blue,1.0,
2498 (void) SetPixelChannelMask(image,channel_mask);
2500 if (((GetPixelBlackTraits(image) & UpdatePixelTrait) != 0) &&
2501 (image->colorspace == CMYKColorspace))
2503 channel_mask=SetPixelChannelMask(image,BlackChannel);
2504 status|=LevelImage(image,black_color->black,white_color->black,1.0,
2506 (void) SetPixelChannelMask(image,channel_mask);
2508 if (((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0) &&
2509 (image->matte == MagickTrue))
2511 channel_mask=SetPixelChannelMask(image,AlphaChannel);
2512 status|=LevelImage(image,black_color->alpha,white_color->alpha,1.0,
2514 (void) SetPixelChannelMask(image,channel_mask);
2519 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
2521 channel_mask=SetPixelChannelMask(image,RedChannel);
2522 status|=LevelizeImage(image,black_color->red,white_color->red,1.0,
2524 (void) SetPixelChannelMask(image,channel_mask);
2526 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
2528 channel_mask=SetPixelChannelMask(image,GreenChannel);
2529 status|=LevelizeImage(image,black_color->green,white_color->green,1.0,
2531 (void) SetPixelChannelMask(image,channel_mask);
2533 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
2535 channel_mask=SetPixelChannelMask(image,BlueChannel);
2536 status|=LevelizeImage(image,black_color->blue,white_color->blue,1.0,
2538 (void) SetPixelChannelMask(image,channel_mask);
2540 if (((GetPixelBlackTraits(image) & UpdatePixelTrait) != 0) &&
2541 (image->colorspace == CMYKColorspace))
2543 channel_mask=SetPixelChannelMask(image,BlackChannel);
2544 status|=LevelizeImage(image,black_color->black,white_color->black,1.0,
2546 (void) SetPixelChannelMask(image,channel_mask);
2548 if (((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0) &&
2549 (image->matte == MagickTrue))
2551 channel_mask=SetPixelChannelMask(image,AlphaChannel);
2552 status|=LevelizeImage(image,black_color->alpha,white_color->alpha,1.0,
2554 (void) SetPixelChannelMask(image,channel_mask);
2557 return(status == 0 ? MagickFalse : MagickTrue);
2561 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2565 % L i n e a r S t r e t c h I m a g e %
2569 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2571 % LinearStretchImage() discards any pixels below the black point and above
2572 % the white point and levels the remaining pixels.
2574 % The format of the LinearStretchImage method is:
2576 % MagickBooleanType LinearStretchImage(Image *image,
2577 % const double black_point,const double white_point,
2578 % ExceptionInfo *exception)
2580 % A description of each parameter follows:
2582 % o image: the image.
2584 % o black_point: the black point.
2586 % o white_point: the white point.
2588 % o exception: return any errors or warnings in this structure.
2591 MagickExport MagickBooleanType LinearStretchImage(Image *image,
2592 const double black_point,const double white_point,ExceptionInfo *exception)
2594 #define LinearStretchImageTag "LinearStretch/Image"
2612 Allocate histogram and linear map.
2614 assert(image != (Image *) NULL);
2615 assert(image->signature == MagickSignature);
2616 histogram=(MagickRealType *) AcquireQuantumMemory(MaxMap+1UL,
2617 sizeof(*histogram));
2618 if (histogram == (MagickRealType *) NULL)
2619 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
2624 (void) ResetMagickMemory(histogram,0,(MaxMap+1)*sizeof(*histogram));
2625 image_view=AcquireCacheView(image);
2626 for (y=0; y < (ssize_t) image->rows; y++)
2628 register const Quantum
2634 p=GetCacheViewVirtualPixels(image_view,0,y,image->columns,1,exception);
2635 if (p == (const Quantum *) NULL)
2637 for (x=0; x < (ssize_t) image->columns; x++)
2639 histogram[ScaleQuantumToMap(GetPixelIntensity(image,p))]++;
2640 p+=GetPixelChannels(image);
2643 image_view=DestroyCacheView(image_view);
2645 Find the histogram boundaries by locating the black and white point levels.
2648 for (black=0; black < (ssize_t) MaxMap; black++)
2650 intensity+=histogram[black];
2651 if (intensity >= black_point)
2655 for (white=(ssize_t) MaxMap; white != 0; white--)
2657 intensity+=histogram[white];
2658 if (intensity >= white_point)
2661 histogram=(MagickRealType *) RelinquishMagickMemory(histogram);
2662 status=LevelImage(image,(double) black,(double) white,1.0,exception);
2667 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2671 % M o d u l a t e I m a g e %
2675 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2677 % ModulateImage() lets you control the brightness, saturation, and hue
2678 % of an image. Modulate represents the brightness, saturation, and hue
2679 % as one parameter (e.g. 90,150,100). If the image colorspace is HSL, the
2680 % modulation is lightness, saturation, and hue. And if the colorspace is
2681 % HWB, use blackness, whiteness, and hue.
2683 % The format of the ModulateImage method is:
2685 % MagickBooleanType ModulateImage(Image *image,const char *modulate,
2686 % ExceptionInfo *exception)
2688 % A description of each parameter follows:
2690 % o image: the image.
2692 % o modulate: Define the percent change in brightness, saturation, and hue.
2694 % o exception: return any errors or warnings in this structure.
2698 static void ModulateHSB(const double percent_hue,
2699 const double percent_saturation,const double percent_brightness,double *red,
2700 double *green,double *blue)
2708 Increase or decrease color brightness, saturation, or hue.
2710 assert(red != (double *) NULL);
2711 assert(green != (double *) NULL);
2712 assert(blue != (double *) NULL);
2713 ConvertRGBToHSB(*red,*green,*blue,&hue,&saturation,&brightness);
2714 hue+=0.5*(0.01*percent_hue-1.0);
2719 saturation*=0.01*percent_saturation;
2720 brightness*=0.01*percent_brightness;
2721 ConvertHSBToRGB(hue,saturation,brightness,red,green,blue);
2724 static void ModulateHSL(const double percent_hue,
2725 const double percent_saturation,const double percent_lightness,double *red,
2726 double *green,double *blue)
2734 Increase or decrease color lightness, saturation, or hue.
2736 assert(red != (double *) NULL);
2737 assert(green != (double *) NULL);
2738 assert(blue != (double *) NULL);
2739 ConvertRGBToHSL(*red,*green,*blue,&hue,&saturation,&lightness);
2740 hue+=0.5*(0.01*percent_hue-1.0);
2745 saturation*=0.01*percent_saturation;
2746 lightness*=0.01*percent_lightness;
2747 ConvertHSLToRGB(hue,saturation,lightness,red,green,blue);
2750 static void ModulateHWB(const double percent_hue,const double percent_whiteness, const double percent_blackness,double *red,double *green,double *blue)
2758 Increase or decrease color blackness, whiteness, or hue.
2760 assert(red != (double *) NULL);
2761 assert(green != (double *) NULL);
2762 assert(blue != (double *) NULL);
2763 ConvertRGBToHWB(*red,*green,*blue,&hue,&whiteness,&blackness);
2764 hue+=0.5*(0.01*percent_hue-1.0);
2769 blackness*=0.01*percent_blackness;
2770 whiteness*=0.01*percent_whiteness;
2771 ConvertHWBToRGB(hue,whiteness,blackness,red,green,blue);
2774 MagickExport MagickBooleanType ModulateImage(Image *image,const char *modulate,
2775 ExceptionInfo *exception)
2777 #define ModulateImageTag "Modulate/Image"
2812 Initialize modulate table.
2814 assert(image != (Image *) NULL);
2815 assert(image->signature == MagickSignature);
2816 if (image->debug != MagickFalse)
2817 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
2818 if (modulate == (char *) NULL)
2819 return(MagickFalse);
2820 flags=ParseGeometry(modulate,&geometry_info);
2821 percent_brightness=geometry_info.rho;
2822 percent_saturation=geometry_info.sigma;
2823 if ((flags & SigmaValue) == 0)
2824 percent_saturation=100.0;
2825 percent_hue=geometry_info.xi;
2826 if ((flags & XiValue) == 0)
2828 colorspace=UndefinedColorspace;
2829 artifact=GetImageArtifact(image,"modulate:colorspace");
2830 if (artifact != (const char *) NULL)
2831 colorspace=(ColorspaceType) ParseCommandOption(MagickColorspaceOptions,
2832 MagickFalse,artifact);
2833 if (image->storage_class == PseudoClass)
2838 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2839 #pragma omp parallel for schedule(static,4) shared(progress,status)
2841 for (i=0; i < (ssize_t) image->colors; i++)
2846 ModulateHSB(percent_hue,percent_saturation,percent_brightness,
2847 &image->colormap[i].red,&image->colormap[i].green,
2848 &image->colormap[i].blue);
2854 ModulateHSL(percent_hue,percent_saturation,percent_brightness,
2855 &image->colormap[i].red,&image->colormap[i].green,
2856 &image->colormap[i].blue);
2861 ModulateHWB(percent_hue,percent_saturation,percent_brightness,
2862 &image->colormap[i].red,&image->colormap[i].green,
2863 &image->colormap[i].blue);
2873 image_view=AcquireCacheView(image);
2874 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2875 #pragma omp parallel for schedule(static,4) shared(progress,status)
2877 for (y=0; y < (ssize_t) image->rows; y++)
2890 if (status == MagickFalse)
2892 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
2893 if (q == (Quantum *) NULL)
2898 for (x=0; x < (ssize_t) image->columns; x++)
2900 red=(double) GetPixelRed(image,q);
2901 green=(double) GetPixelGreen(image,q);
2902 blue=(double) GetPixelBlue(image,q);
2907 ModulateHSB(percent_hue,percent_saturation,percent_brightness,
2914 ModulateHSL(percent_hue,percent_saturation,percent_brightness,
2920 ModulateHWB(percent_hue,percent_saturation,percent_brightness,
2925 SetPixelRed(image,ClampToQuantum(red),q);
2926 SetPixelGreen(image,ClampToQuantum(green),q);
2927 SetPixelBlue(image,ClampToQuantum(blue),q);
2928 q+=GetPixelChannels(image);
2930 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
2932 if (image->progress_monitor != (MagickProgressMonitor) NULL)
2937 #if defined(MAGICKCORE_OPENMP_SUPPORT)
2938 #pragma omp critical (MagickCore_ModulateImage)
2940 proceed=SetImageProgress(image,ModulateImageTag,progress++,image->rows);
2941 if (proceed == MagickFalse)
2945 image_view=DestroyCacheView(image_view);
2950 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2954 % N e g a t e I m a g e %
2958 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2960 % NegateImage() negates the colors in the reference image. The grayscale
2961 % option means that only grayscale values within the image are negated.
2963 % The format of the NegateImage method is:
2965 % MagickBooleanType NegateImage(Image *image,
2966 % const MagickBooleanType grayscale,ExceptionInfo *exception)
2968 % A description of each parameter follows:
2970 % o image: the image.
2972 % o grayscale: If MagickTrue, only negate grayscale pixels within the image.
2974 % o exception: return any errors or warnings in this structure.
2977 MagickExport MagickBooleanType NegateImage(Image *image,
2978 const MagickBooleanType grayscale,ExceptionInfo *exception)
2980 #define NegateImageTag "Negate/Image"
2997 assert(image != (Image *) NULL);
2998 assert(image->signature == MagickSignature);
2999 if (image->debug != MagickFalse)
3000 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
3001 if (image->storage_class == PseudoClass)
3006 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3007 #pragma omp parallel for schedule(static,4) shared(progress,status)
3009 for (i=0; i < (ssize_t) image->colors; i++)
3011 if (grayscale != MagickFalse)
3012 if ((image->colormap[i].red != image->colormap[i].green) ||
3013 (image->colormap[i].green != image->colormap[i].blue))
3015 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
3016 image->colormap[i].red=(Quantum) QuantumRange-
3017 image->colormap[i].red;
3018 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
3019 image->colormap[i].green=(Quantum) QuantumRange-
3020 image->colormap[i].green;
3021 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
3022 image->colormap[i].blue=(Quantum) QuantumRange-
3023 image->colormap[i].blue;
3031 image_view=AcquireCacheView(image);
3032 if (grayscale != MagickFalse)
3034 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3035 #pragma omp parallel for schedule(static,4) shared(progress,status)
3037 for (y=0; y < (ssize_t) image->rows; y++)
3048 if (status == MagickFalse)
3050 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
3052 if (q == (Quantum *) NULL)
3057 for (x=0; x < (ssize_t) image->columns; x++)
3062 if (IsPixelGray(image,q) != MagickFalse)
3064 q+=GetPixelChannels(image);
3067 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
3075 channel=GetPixelChannelMapChannel(image,i);
3076 traits=GetPixelChannelMapTraits(image,channel);
3077 if ((traits & UpdatePixelTrait) != 0)
3078 q[i]=QuantumRange-q[i];
3080 q+=GetPixelChannels(image);
3082 sync=SyncCacheViewAuthenticPixels(image_view,exception);
3083 if (sync == MagickFalse)
3085 if (image->progress_monitor != (MagickProgressMonitor) NULL)
3090 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3091 #pragma omp critical (MagickCore_NegateImage)
3093 proceed=SetImageProgress(image,NegateImageTag,progress++,
3095 if (proceed == MagickFalse)
3099 image_view=DestroyCacheView(image_view);
3105 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3106 #pragma omp parallel for schedule(static,4) shared(progress,status)
3108 for (y=0; y < (ssize_t) image->rows; y++)
3116 if (status == MagickFalse)
3118 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
3119 if (q == (Quantum *) NULL)
3124 for (x=0; x < (ssize_t) image->columns; x++)
3129 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
3137 channel=GetPixelChannelMapChannel(image,i);
3138 traits=GetPixelChannelMapTraits(image,channel);
3139 if ((traits & UpdatePixelTrait) != 0)
3140 q[i]=QuantumRange-q[i];
3142 q+=GetPixelChannels(image);
3144 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
3146 if (image->progress_monitor != (MagickProgressMonitor) NULL)
3151 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3152 #pragma omp critical (MagickCore_NegateImage)
3154 proceed=SetImageProgress(image,NegateImageTag,progress++,image->rows);
3155 if (proceed == MagickFalse)
3159 image_view=DestroyCacheView(image_view);
3164 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3168 % N o r m a l i z e I m a g e %
3172 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3174 % NormalizeImage() enhances the contrast of a color image by mapping the
3175 % darkest 2 percent of all pixel to black and the brightest 1 percent to white.
3177 % The format of the NormalizeImage method is:
3179 % MagickBooleanType NormalizeImage(Image *image,ExceptionInfo *exception)
3181 % A description of each parameter follows:
3183 % o image: the image.
3185 % o exception: return any errors or warnings in this structure.
3188 MagickExport MagickBooleanType NormalizeImage(Image *image,
3189 ExceptionInfo *exception)
3195 black_point=(double) image->columns*image->rows*0.0015;
3196 white_point=(double) image->columns*image->rows*0.9995;
3197 return(ContrastStretchImage(image,black_point,white_point,exception));
3201 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3205 % S i g m o i d a l C o n t r a s t I m a g e %
3209 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3211 % SigmoidalContrastImage() adjusts the contrast of an image with a non-linear
3212 % sigmoidal contrast algorithm. Increase the contrast of the image using a
3213 % sigmoidal transfer function without saturating highlights or shadows.
3214 % Contrast indicates how much to increase the contrast (0 is none; 3 is
3215 % typical; 20 is pushing it); mid-point indicates where midtones fall in the
3216 % resultant image (0 is white; 50% is middle-gray; 100% is black). Set
3217 % sharpen to MagickTrue to increase the image contrast otherwise the contrast
3220 % The format of the SigmoidalContrastImage method is:
3222 % MagickBooleanType SigmoidalContrastImage(Image *image,
3223 % const MagickBooleanType sharpen,const char *levels,
3224 % ExceptionInfo *exception)
3226 % A description of each parameter follows:
3228 % o image: the image.
3230 % o sharpen: Increase or decrease image contrast.
3232 % o alpha: strength of the contrast, the larger the number the more
3233 % 'threshold-like' it becomes.
3235 % o beta: midpoint of the function as a color value 0 to QuantumRange.
3237 % o exception: return any errors or warnings in this structure.
3240 MagickExport MagickBooleanType SigmoidalContrastImage(Image *image,
3241 const MagickBooleanType sharpen,const double contrast,const double midpoint,
3242 ExceptionInfo *exception)
3244 #define SigmoidalContrastImageTag "SigmoidalContrast/Image"
3265 Allocate and initialize sigmoidal maps.
3267 assert(image != (Image *) NULL);
3268 assert(image->signature == MagickSignature);
3269 if (image->debug != MagickFalse)
3270 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
3271 sigmoidal_map=(MagickRealType *) AcquireQuantumMemory(MaxMap+1UL,
3272 sizeof(*sigmoidal_map));
3273 if (sigmoidal_map == (MagickRealType *) NULL)
3274 ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
3276 (void) ResetMagickMemory(sigmoidal_map,0,(MaxMap+1)*sizeof(*sigmoidal_map));
3277 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3278 #pragma omp parallel for schedule(static,4) shared(progress,status)
3280 for (i=0; i <= (ssize_t) MaxMap; i++)
3282 if (sharpen != MagickFalse)
3284 sigmoidal_map[i]=(MagickRealType) ScaleMapToQuantum((MagickRealType)
3285 (MaxMap*((1.0/(1.0+exp(contrast*(midpoint/(double) QuantumRange-
3286 (double) i/MaxMap))))-(1.0/(1.0+exp(contrast*(midpoint/(double)
3287 QuantumRange)))))/((1.0/(1.0+exp(contrast*(midpoint/(double)
3288 QuantumRange-1.0))))-(1.0/(1.0+exp(contrast*(midpoint/(double)
3289 QuantumRange)))))+0.5));
3292 sigmoidal_map[i]=(MagickRealType) ScaleMapToQuantum((MagickRealType)
3293 (MaxMap*(QuantumScale*midpoint-log((1.0-(1.0/(1.0+exp(midpoint/(double)
3294 QuantumRange*contrast))+((double) i/MaxMap)*((1.0/(1.0+exp(contrast*(
3295 midpoint/(double) QuantumRange-1.0))))-(1.0/(1.0+exp(midpoint/(double)
3296 QuantumRange*contrast))))))/(1.0/(1.0+exp(midpoint/(double) QuantumRange*
3297 contrast))+((double) i/MaxMap)*((1.0/(1.0+exp(contrast*(midpoint/(double)
3298 QuantumRange-1.0))))-(1.0/(1.0+exp(midpoint/(double) QuantumRange*
3299 contrast))))))/contrast)));
3301 if (image->storage_class == PseudoClass)
3304 Sigmoidal-contrast enhance colormap.
3306 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3307 #pragma omp parallel for schedule(static,4) shared(progress,status)
3309 for (i=0; i < (ssize_t) image->colors; i++)
3311 if ((GetPixelRedTraits(image) & UpdatePixelTrait) != 0)
3312 image->colormap[i].red=sigmoidal_map[ScaleQuantumToMap(
3313 ClampToQuantum(image->colormap[i].red))];
3314 if ((GetPixelGreenTraits(image) & UpdatePixelTrait) != 0)
3315 image->colormap[i].green=sigmoidal_map[ScaleQuantumToMap(
3316 ClampToQuantum(image->colormap[i].green))];
3317 if ((GetPixelBlueTraits(image) & UpdatePixelTrait) != 0)
3318 image->colormap[i].blue=sigmoidal_map[ScaleQuantumToMap(
3319 ClampToQuantum(image->colormap[i].blue))];
3320 if ((GetPixelAlphaTraits(image) & UpdatePixelTrait) != 0)
3321 image->colormap[i].alpha=sigmoidal_map[ScaleQuantumToMap(
3322 ClampToQuantum(image->colormap[i].alpha))];
3326 Sigmoidal-contrast enhance image.
3330 image_view=AcquireCacheView(image);
3331 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3332 #pragma omp parallel for schedule(static,4) shared(progress,status)
3334 for (y=0; y < (ssize_t) image->rows; y++)
3342 if (status == MagickFalse)
3344 q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
3345 if (q == (Quantum *) NULL)
3350 for (x=0; x < (ssize_t) image->columns; x++)
3355 for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
3363 channel=GetPixelChannelMapChannel(image,i);
3364 traits=GetPixelChannelMapTraits(image,channel);
3365 if ((traits & UpdatePixelTrait) != 0)
3366 q[i]=ClampToQuantum(sigmoidal_map[ScaleQuantumToMap(q[i])]);
3368 q+=GetPixelChannels(image);
3370 if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
3372 if (image->progress_monitor != (MagickProgressMonitor) NULL)
3377 #if defined(MAGICKCORE_OPENMP_SUPPORT)
3378 #pragma omp critical (MagickCore_SigmoidalContrastImage)
3380 proceed=SetImageProgress(image,SigmoidalContrastImageTag,progress++,
3382 if (proceed == MagickFalse)
3386 image_view=DestroyCacheView(image_view);
3387 sigmoidal_map=(MagickRealType *) RelinquishMagickMemory(sigmoidal_map);