]> granicus.if.org Git - imagemagick/blobdiff - MagickCore/composite.c
sigmoidal-contrast: prevent argument out of range and remove unnecessary ClampToQuantum
[imagemagick] / MagickCore / composite.c
index 3cb275b9f7f096890ffcf84c0d758e4feaa8095a..788bae850f048999970b5d918d906bc6e725938f 100644 (file)
@@ -72,6 +72,7 @@
 #include "MagickCore/resource_.h"
 #include "MagickCore/string_.h"
 #include "MagickCore/thread-private.h"
+#include "MagickCore/threshold.h"
 #include "MagickCore/token.h"
 #include "MagickCore/utility.h"
 #include "MagickCore/utility-private.h"
        performed as 'pure' mathematical operations, rather than as image
        operations.
 */
-static void CompositeHSB(const double red,const double green,
-  const double blue,double *hue,double *saturation,double *brightness)
+
+static inline double MagickMin(const double x,const double y)
 {
-  double
-    delta,
-    max,
-    min;
+  if (x < y)
+    return(x);
+  return(y);
+}
 
-  /*
-    Convert RGB to HSB colorspace.
-  */
-  assert(hue != (double *) NULL);
-  assert(saturation != (double *) NULL);
-  assert(brightness != (double *) NULL);
-  max=(red > green ? red : green);
-  if (blue > max)
-    max=blue;
-  min=(red < green ? red : green);
-  if (blue < min)
-    min=blue;
-  *hue=0.0;
-  *saturation=0.0;
-  *brightness=(double) (QuantumScale*max);
-  if (fabs((double) max) < MagickEpsilon)
-    return;
-  *saturation=(double) (1.0-min/max);
-  delta=(MagickRealType) max-min;
-  if (fabs(delta) < MagickEpsilon)
-    return;
-  if (fabs((double) red-max) < MagickEpsilon)
-    *hue=(double) ((green-blue)/delta);
-  else
-    if (fabs((double) green-max) < MagickEpsilon)
-      *hue=(double) (2.0+(blue-red)/delta);
-    else
-      if (fabs((double) blue-max) < MagickEpsilon)
-        *hue=(double) (4.0+(red-green)/delta);
-  *hue/=6.0;
-  if (*hue < 0.0)
-    *hue+=1.0;
+static inline double MagickMax(const double x,const double y)
+{
+  if (x > y)
+    return(x);
+  return(y);
+}
+
+static inline double ConvertHueToRGB(double m1,
+  double m2,double hue)
+{
+  if (hue < 0.0)
+    hue+=1.0;
+  if (hue > 1.0)
+    hue-=1.0;
+  if ((6.0*hue) < 1.0)
+    return(m1+6.0*(m2-m1)*hue);
+  if ((2.0*hue) < 1.0)
+    return(m2);
+  if ((3.0*hue) < 2.0)
+    return(m1+6.0*(m2-m1)*(2.0/3.0-hue));
+  return(m1);
 }
 
-static void HSBComposite(const double hue,const double saturation,
-  const double brightness,double *red,double *green,double *blue)
+static void HCLComposite(const double hue,const double chroma,const double luma,
+  double *red,double *green,double *blue)
 {
   double
-    f,
+    b,
+    c,
+    g,
     h,
-    p,
-    q,
-    t;
+    m,
+    r,
+    x,
+    z;
 
   /*
-    Convert HSB to RGB colorspace.
+    Convert HCL to RGB colorspace.
   */
   assert(red != (double *) NULL);
   assert(green != (double *) NULL);
   assert(blue != (double *) NULL);
-  if (saturation == 0.0)
-    {
-      *red=(double) QuantumRange*brightness;
-      *green=(*red);
-      *blue=(*red);
-      return;
-    }
-  h=6.0*(hue-floor(hue));
-  f=h-floor((double) h);
-  p=brightness*(1.0-saturation);
-  q=brightness*(1.0-saturation*f);
-  t=brightness*(1.0-saturation*(1.0-f));
-  switch ((int) h)
-  {
-    case 0:
-    default:
-    {
-      *red=(double) QuantumRange*brightness;
-      *green=(double) QuantumRange*t;
-      *blue=(double) QuantumRange*p;
-      break;
-    }
-    case 1:
-    {
-      *red=(double) QuantumRange*q;
-      *green=(double) QuantumRange*brightness;
-      *blue=(double) QuantumRange*p;
-      break;
-    }
-    case 2:
-    {
-      *red=(double) QuantumRange*p;
-      *green=(double) QuantumRange*brightness;
-      *blue=(double) QuantumRange*t;
-      break;
-    }
-    case 3:
-    {
-      *red=(double) QuantumRange*p;
-      *green=(double) QuantumRange*q;
-      *blue=(double) QuantumRange*brightness;
-      break;
-    }
-    case 4:
+  h=6.0*hue;
+  c=chroma;
+  x=c*(1.0-fabs(fmod(h,2.0)-1.0));
+  r=0.0;
+  g=0.0;
+  b=0.0;
+  if ((0.0 <= h) && (h < 1.0))
     {
-      *red=(double) QuantumRange*t;
-      *green=(double) QuantumRange*p;
-      *blue=(double) QuantumRange*brightness;
-      break;
+      r=c;
+      g=x;
     }
-    case 5:
+  else
+    if ((1.0 <= h) && (h < 2.0))
+      {
+        r=x;
+        g=c;
+      }
+    else
+      if ((2.0 <= h) && (h < 3.0))
+        {
+          g=c;
+          b=x;
+        }
+      else
+        if ((3.0 <= h) && (h < 4.0))
+          {
+            g=x;
+            b=c;
+          }
+        else
+          if ((4.0 <= h) && (h < 5.0))
+            {
+              r=x;
+              b=c;
+            }
+          else
+            if ((5.0 <= h) && (h < 6.0))
+              {
+                r=c;
+                b=x;
+              }
+  m=luma-(0.298839*r+0.586811*g+0.114350*b);
+  /*
+    Choose saturation strategy to clip it into the RGB cube; hue and luma are
+    preserved and chroma may be changed.
+  */
+  z=1.0;
+  if (m < 0.0)
     {
-      *red=(double) QuantumRange*brightness;
-      *green=(double) QuantumRange*p;
-      *blue=(double) QuantumRange*q;
-      break;
+      z=luma/(luma-m);
+      m=0.0;
     }
-  }
+  else
+    if (m+c > 1.0)
+      {
+        z=(1.0-luma)/(m+c-luma);
+        m=1.0-z*c;
+      }
+  *red=QuantumRange*(z*r+m);
+  *green=QuantumRange*(z*g+m);
+  *blue=QuantumRange*(z*b+m);
 }
 
-static inline double MagickMin(const double x,const double y)
+static void CompositeHCL(const double red,const double green,const double blue,
+  double *hue,double *chroma,double *luma)
 {
-  if (x < y)
-    return(x);
-  return(y);
-}
+  double
+    b,
+    c,
+    g,
+    h,
+    max,
+    r;
 
-static inline double MagickMax(const double x,const double y)
-{
-  if (x > y)
-    return(x);
-  return(y);
+  /*
+    Convert RGB to HCL colorspace.
+  */
+  assert(hue != (double *) NULL);
+  assert(chroma != (double *) NULL);
+  assert(luma != (double *) NULL);
+  r=red;
+  g=green;
+  b=blue;
+  max=MagickMax(r,MagickMax(g,b));
+  c=max-(double) MagickMin(r,MagickMin(g,b));
+  h=0.0;
+  if (c == 0)
+    h=0.0;
+  else
+    if (red == max)
+      h=fmod(6.0+(g-b)/c,6.0);
+    else
+      if (green == max)
+        h=((b-r)/c)+2.0;
+      else
+        if (blue == max)
+          h=((r-g)/c)+4.0;
+  *hue=(h/6.0);
+  *chroma=QuantumScale*c;
+  *luma=QuantumScale*(0.298839*r+0.586811*g+0.114350*b);
 }
 
 static MagickBooleanType CompositeOverImage(Image *image,
@@ -389,7 +410,7 @@ static MagickBooleanType CompositeOverImage(Image *image,
       }
     for (x=0; x < (ssize_t) image->columns; x++)
     {
-      MagickRealType
+      double
         alpha,
         Da,
         Dc,
@@ -437,9 +458,9 @@ static MagickBooleanType CompositeOverImage(Image *image,
               composite_traits,
               traits;
 
-            channel=GetPixelChannelMapChannel(image,i);
-            traits=GetPixelChannelMapTraits(image,channel);
-            composite_traits=GetPixelChannelMapTraits(composite_image,channel);
+            channel=GetPixelChannelChannel(image,i);
+            traits=GetPixelChannelTraits(image,channel);
+            composite_traits=GetPixelChannelTraits(composite_image,channel);
             if ((traits == UndefinedPixelTrait) ||
                 (composite_traits == UndefinedPixelTrait))
               continue;
@@ -474,9 +495,9 @@ static MagickBooleanType CompositeOverImage(Image *image,
           composite_traits,
           traits;
 
-        channel=GetPixelChannelMapChannel(image,i);
-        traits=GetPixelChannelMapTraits(image,channel);
-        composite_traits=GetPixelChannelMapTraits(composite_image,channel);
+        channel=GetPixelChannelChannel(image,i);
+        traits=GetPixelChannelTraits(image,channel);
+        composite_traits=GetPixelChannelTraits(composite_image,channel);
         if ((traits == UndefinedPixelTrait) ||
             (composite_traits == UndefinedPixelTrait))
           continue;
@@ -500,8 +521,8 @@ static MagickBooleanType CompositeOverImage(Image *image,
           Sc: source color.
           Dc: destination color.
         */
-        Sc=(MagickRealType) GetPixelChannel(composite_image,channel,p);
-        Dc=(MagickRealType) q[i];
+        Sc=(double) GetPixelChannel(composite_image,channel,p);
+        Dc=(double) q[i];
         gamma=MagickEpsilonReciprocal(alpha);
         q[i]=ClampToQuantum(gamma*(Sa*Sc-Sa*Da*Dc+Da*Dc));
       }
@@ -556,12 +577,12 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   MagickOffsetType
     progress;
 
-  MagickRealType
+  double
     amount,
     destination_dissolve,
     midpoint,
-    percent_brightness,
-    percent_saturation,
+    percent_luma,
+    percent_chroma,
     source_dissolve,
     threshold;
 
@@ -582,7 +603,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   composite_image=CloneImage(composite,0,0,MagickTrue,exception);
   if (composite_image == (const Image *) NULL)
     return(MagickFalse);
-  (void) TransformImageColorspace(composite_image,image->colorspace,exception);
+  (void) SetImageColorspace(composite_image,image->colorspace,exception);
   if ((compose == OverCompositeOp) || (compose == SrcOverCompositeOp))
     {
       status=CompositeOverImage(image,composite_image,clip_to_self,x_offset,
@@ -593,8 +614,8 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   destination_image=(Image *) NULL;
   amount=0.5;
   destination_dissolve=1.0;
-  percent_brightness=100.0;
-  percent_saturation=100.0;
+  percent_luma=100.0;
+  percent_chroma=100.0;
   source_dissolve=1.0;
   threshold=0.05f;
   switch (compose)
@@ -659,9 +680,9 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               composite_traits,
               traits;
 
-            channel=GetPixelChannelMapChannel(composite_image,i);
-            composite_traits=GetPixelChannelMapTraits(composite_image,channel);
-            traits=GetPixelChannelMapTraits(image,channel);
+            channel=GetPixelChannelChannel(composite_image,i);
+            composite_traits=GetPixelChannelTraits(composite_image,channel);
+            traits=GetPixelChannelTraits(image,channel);
             if ((traits == UndefinedPixelTrait) ||
                 (composite_traits == UndefinedPixelTrait))
               continue;
@@ -700,7 +721,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
         Modify destination outside the overlaid region and require an alpha
         channel to exist, to add transparency.
       */
-      if (image->matte == MagickFalse)
+      if (image->alpha_trait != BlendPixelTrait)
         (void) SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
       break;
     }
@@ -716,7 +737,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       PixelInfo
         pixel;
 
-      MagickRealType
+      double
         angle_range,
         angle_start,
         height,
@@ -774,7 +795,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       /* rotate vectors if a rotation angle is given */
       if ((flags & XValue) != 0 )
         {
-          MagickRealType
+          double
             angle;
 
           angle=DegreesToRadians(geometry_info.xi);
@@ -840,7 +861,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
             }
           if (fabs(angle_range) > MagickEpsilon)
             {
-              MagickRealType
+              double
                 angle;
 
               angle=angle_start+angle_range*QuantumScale*
@@ -893,7 +914,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       PixelInfo
         pixel;
 
-      MagickRealType
+      double
         horizontal_scale,
         vertical_scale;
 
@@ -922,14 +943,14 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
         {
           if ((flags & AspectValue) == 0)
             {
-              horizontal_scale=(MagickRealType) (composite_image->columns-1.0)/
+              horizontal_scale=(double) (composite_image->columns-1.0)/
                 2.0;
-              vertical_scale=(MagickRealType) (composite_image->rows-1.0)/2.0;
+              vertical_scale=(double) (composite_image->rows-1.0)/2.0;
             }
           else
             {
-              horizontal_scale=(MagickRealType) (image->columns-1.0)/2.0;
-              vertical_scale=(MagickRealType) (image->rows-1.0)/2.0;
+              horizontal_scale=(double) (image->columns-1.0)/2.0;
+              vertical_scale=(double) (image->rows-1.0)/2.0;
             }
         }
       else
@@ -960,29 +981,29 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
            default = center of overlay image
            arg flag '!' = locations/percentage relative to background image
       */
-      center.x=(MagickRealType) x_offset;
-      center.y=(MagickRealType) y_offset;
+      center.x=(double) x_offset;
+      center.y=(double) y_offset;
       if (compose == DistortCompositeOp)
         {
           if ((flags & XValue) == 0)
             if ((flags & AspectValue) == 0)
-              center.x=(MagickRealType) x_offset+(composite_image->columns-1)/
+              center.x=(double) x_offset+(composite_image->columns-1)/
                 2.0;
             else
-              center.x=((MagickRealType) image->columns-1)/2.0;
+              center.x=((double) image->columns-1)/2.0;
           else
             if ((flags & AspectValue) == 0)
-              center.x=(MagickRealType) x_offset+geometry_info.xi;
+              center.x=(double) x_offset+geometry_info.xi;
             else
               center.x=geometry_info.xi;
           if ((flags & YValue) == 0)
             if ((flags & AspectValue) == 0)
-              center.y=(MagickRealType) y_offset+(composite_image->rows-1)/2.0;
+              center.y=(double) y_offset+(composite_image->rows-1)/2.0;
             else
-              center.y=((MagickRealType) image->rows-1)/2.0;
+              center.y=((double) image->rows-1)/2.0;
           else
             if ((flags & AspectValue) == 0)
-              center.y=(MagickRealType) y_offset+geometry_info.psi;
+              center.y=(double) y_offset+geometry_info.psi;
             else
               center.y=geometry_info.psi;
         }
@@ -1027,11 +1048,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
             Displace the offset.
           */
           offset.x=(horizontal_scale*(GetPixelRed(composite_image,p)-
-            (((MagickRealType) QuantumRange+1.0)/2.0)))/(((MagickRealType)
+            (((double) QuantumRange+1.0)/2.0)))/(((double)
             QuantumRange+1.0)/2.0)+center.x+((compose == DisplaceCompositeOp) ?
             x : 0);
           offset.y=(vertical_scale*(GetPixelGreen(composite_image,p)-
-            (((MagickRealType) QuantumRange+1.0)/2.0)))/(((MagickRealType)
+            (((double) QuantumRange+1.0)/2.0)))/(((double)
             QuantumRange+1.0)/2.0)+center.y+((compose == DisplaceCompositeOp) ?
             y : 0);
           (void) InterpolatePixelInfo(image,image_view,
@@ -1040,7 +1061,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           /*
             Mask with the 'invalid pixel mask' in alpha channel.
           */
-          pixel.alpha=(MagickRealType) QuantumRange*(1.0-(1.0-QuantumScale*
+          pixel.alpha=(double) QuantumRange*(1.0-(1.0-QuantumScale*
             pixel.alpha)*(1.0-QuantumScale*GetPixelAlpha(composite_image,p)));
           SetPixelInfoPixel(destination_image,&pixel,q);
           p+=GetPixelChannels(composite_image);
@@ -1134,15 +1155,15 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
         *value;
 
       /*
-        Determine the brightness and saturation scale.
+        Determine the luma and chroma scale.
       */
       value=GetImageArtifact(composite_image,"compose:args");
       if (value != (char *) NULL)
         {
           flags=ParseGeometry(value,&geometry_info);
-          percent_brightness=geometry_info.rho;
+          percent_luma=geometry_info.rho;
           if ((flags & SigmaValue) != 0)
-            percent_saturation=geometry_info.sigma;
+            percent_chroma=geometry_info.sigma;
         }
       break;
     }
@@ -1174,7 +1195,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   */
   status=MagickTrue;
   progress=0;
-  midpoint=((MagickRealType) QuantumRange+1.0)/2;
+  midpoint=((double) QuantumRange+1.0)/2;
   composite_view=AcquireVirtualCacheView(composite_image,exception);
   image_view=AcquireAuthenticCacheView(image,exception);
 #if defined(MAGICKCORE_OPENMP_SUPPORT)
@@ -1188,11 +1209,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
 
     double
       blue,
-      brightness,
+      luma,
       green,
       hue,
       red,
-      saturation;
+      chroma;
 
     PixelInfo
       destination_pixel,
@@ -1241,13 +1262,13 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
         continue;
       }
     hue=0.0;
-    saturation=0.0;
-    brightness=0.0;
+    chroma=0.0;
+    luma=0.0;
     GetPixelInfo(image,&destination_pixel);
     GetPixelInfo(composite_image,&source_pixel);
     for (x=0; x < (ssize_t) image->columns; x++)
     {
-      MagickRealType
+      double
         alpha,
         Da,
         Dc,
@@ -1293,7 +1314,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
             }
           for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
           {
-            MagickRealType
+            double
               pixel;
 
             PixelChannel
@@ -1303,9 +1324,9 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               composite_traits,
               traits;
 
-            channel=GetPixelChannelMapChannel(image,i);
-            traits=GetPixelChannelMapTraits(image,channel);
-            composite_traits=GetPixelChannelMapTraits(composite_image,channel);
+            channel=GetPixelChannelChannel(image,i);
+            traits=GetPixelChannelTraits(image,channel);
+            composite_traits=GetPixelChannelTraits(composite_image,channel);
             if ((traits == UndefinedPixelTrait) ||
                 (composite_traits == UndefinedPixelTrait))
               continue;
@@ -1322,9 +1343,9 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               case SrcInCompositeOp:
               case SrcOutCompositeOp:
               {
-                pixel=(MagickRealType) q[i];
+                pixel=(double) q[i];
                 if (channel == AlphaPixelChannel)
-                  pixel=(MagickRealType) TransparentAlpha;
+                  pixel=(double) TransparentAlpha;
                 break;
               }
               case ClearCompositeOp:
@@ -1334,7 +1355,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               {
                 if (channel == AlphaPixelChannel)
                   {
-                    pixel=(MagickRealType) TransparentAlpha;
+                    pixel=(double) TransparentAlpha;
                     break;
                   }
                 pixel=0.0;
@@ -1349,12 +1370,12 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                       source);
                     break;
                   }
-                pixel=(MagickRealType) source[channel];
+                pixel=(double) source[channel];
                 break;
               }
               default:
               {
-                pixel=(MagickRealType) source[channel];
+                pixel=(double) source[channel];
                 break;
               }
             }
@@ -1484,7 +1505,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
         double
           sans;
 
-        MagickRealType
+        double
           pixel;
 
         PixelChannel
@@ -1494,9 +1515,9 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           composite_traits,
           traits;
 
-        channel=GetPixelChannelMapChannel(image,i);
-        traits=GetPixelChannelMapTraits(image,channel);
-        composite_traits=GetPixelChannelMapTraits(composite_image,channel);
+        channel=GetPixelChannelChannel(image,i);
+        traits=GetPixelChannelTraits(image,channel);
+        composite_traits=GetPixelChannelTraits(composite_image,channel);
         if (traits == UndefinedPixelTrait)
           continue;
         if ((compose != IntensityCompositeOp) &&
@@ -1506,8 +1527,8 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           Sc: source color.
           Dc: destination color.
         */
-        Sc=(MagickRealType) GetPixelChannel(composite_image,channel,p);
-        Dc=(MagickRealType) q[i];
+        Sc=(double) GetPixelChannel(composite_image,channel,p);
+        Dc=(double) q[i];
         if ((traits & CopyPixelTrait) != 0)
           {
             if (channel != AlphaPixelChannel)
@@ -1548,23 +1569,23 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 MagickBooleanType
                   equivalent;
 
-                if (Da > ((MagickRealType) QuantumRange/2.0))
+                if (Da > ((double) QuantumRange/2.0))
                   {
-                    pixel=(MagickRealType) TransparentAlpha;
+                    pixel=(double) TransparentAlpha;
                     break;
                   }
                 equivalent=IsFuzzyEquivalencePixel(composite_image,p,image,q);
                 if (equivalent != MagickFalse)
                   {
-                    pixel=(MagickRealType) TransparentAlpha;
+                    pixel=(double) TransparentAlpha;
                     break;
                   }
-                pixel=(MagickRealType) OpaqueAlpha;
+                pixel=(double) OpaqueAlpha;
                 break;
               }
               case ClearCompositeOp:
               {
-                pixel=(MagickRealType) TransparentAlpha;
+                pixel=(double) TransparentAlpha;
                 break;
               }
               case ColorizeCompositeOp:
@@ -1593,7 +1614,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               case CopyAlphaCompositeOp:
               {
                 pixel=QuantumRange*Sa;
-                if (composite_image->matte == MagickFalse)
+                if (composite_image->alpha_trait != BlendPixelTrait)
                   pixel=GetPixelIntensity(composite_image,p);
                 break;
               }
@@ -1615,7 +1636,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               }
               case IntensityCompositeOp:
               {
-                pixel=(MagickRealType) GetPixelIntensity(composite_image,p);
+                pixel=GetPixelIntensity(composite_image,p);
                 break;
               }
               case LightenIntensityCompositeOp:
@@ -1759,11 +1780,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 pixel=Sc;
                 break;
               }
-            CompositeHSB(destination_pixel.red,destination_pixel.green,
-              destination_pixel.blue,&sans,&sans,&brightness);
-            CompositeHSB(source_pixel.red,source_pixel.green,source_pixel.blue,
-              &hue,&saturation,&sans);
-            HSBComposite(hue,saturation,brightness,&red,&green,&blue);
+            CompositeHCL(destination_pixel.red,destination_pixel.green,
+              destination_pixel.blue,&sans,&sans,&luma);
+            CompositeHCL(source_pixel.red,source_pixel.green,source_pixel.blue,
+              &hue,&chroma,&sans);
+            HCLComposite(hue,chroma,luma,&red,&green,&blue);
             switch (channel)
             {
               case RedPixelChannel: pixel=red; break;
@@ -1782,28 +1803,28 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           case CopyBlackCompositeOp:
           {
             if (channel == BlackPixelChannel)
-              pixel=(MagickRealType) GetPixelBlack(composite_image,p);
+              pixel=(double) GetPixelBlack(composite_image,p);
             break;
           }
           case CopyBlueCompositeOp:
           case CopyYellowCompositeOp:
           {
             if (channel == BluePixelChannel)
-              pixel=(MagickRealType) GetPixelBlue(composite_image,p);
+              pixel=(double) GetPixelBlue(composite_image,p);
             break;
           }
           case CopyGreenCompositeOp:
           case CopyMagentaCompositeOp:
           {
             if (channel == GreenPixelChannel)
-              pixel=(MagickRealType) GetPixelGreen(composite_image,p);
+              pixel=(double) GetPixelGreen(composite_image,p);
             break;
           }
           case CopyRedCompositeOp:
           case CopyCyanCompositeOp:
           {
             if (channel == RedPixelChannel)
-              pixel=(MagickRealType) GetPixelRed(composite_image,p);
+              pixel=(double) GetPixelRed(composite_image,p);
             break;
           }
           case DarkenCompositeOp:
@@ -1924,11 +1945,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 pixel=Sc;
                 break;
               }
-            CompositeHSB(destination_pixel.red,destination_pixel.green,
-              destination_pixel.blue,&hue,&saturation,&brightness);
-            CompositeHSB(source_pixel.red,source_pixel.green,source_pixel.blue,
+            CompositeHCL(destination_pixel.red,destination_pixel.green,
+              destination_pixel.blue,&hue,&chroma,&luma);
+            CompositeHCL(source_pixel.red,source_pixel.green,source_pixel.blue,
               &hue,&sans,&sans);
-            HSBComposite(hue,saturation,brightness,&red,&green,&blue);
+            HCLComposite(hue,chroma,luma,&red,&green,&blue);
             switch (channel)
             {
               case RedPixelChannel: pixel=red; break;
@@ -2004,11 +2025,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 pixel=Sc;
                 break;
               }
-            CompositeHSB(destination_pixel.red,destination_pixel.green,
-              destination_pixel.blue,&hue,&saturation,&brightness);
-            CompositeHSB(source_pixel.red,source_pixel.green,source_pixel.blue,
-              &sans,&sans,&brightness);
-            HSBComposite(hue,saturation,brightness,&red,&green,&blue);
+            CompositeHCL(destination_pixel.red,destination_pixel.green,
+              destination_pixel.blue,&hue,&chroma,&luma);
+            CompositeHCL(source_pixel.red,source_pixel.green,source_pixel.blue,
+              &sans,&sans,&luma);
+            HCLComposite(hue,chroma,luma,&red,&green,&blue);
             switch (channel)
             {
               case RedPixelChannel: pixel=red; break;
@@ -2075,11 +2096,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 pixel=Dc;
                 break;
               }
-            CompositeHSB(destination_pixel.red,destination_pixel.green,
-              destination_pixel.blue,&hue,&saturation,&brightness);
-            brightness+=(0.01*percent_brightness*offset)/midpoint;
-            saturation*=0.01*percent_saturation;
-            HSBComposite(hue,saturation,brightness,&red,&green,&blue);
+            CompositeHCL(destination_pixel.red,destination_pixel.green,
+              destination_pixel.blue,&hue,&chroma,&luma);
+            luma+=(0.01*percent_luma*offset)/midpoint;
+            chroma*=0.01*percent_chroma;
+            HCLComposite(hue,chroma,luma,&red,&green,&blue);
             switch (channel)
             {
               case RedPixelChannel: pixel=red; break;
@@ -2191,11 +2212,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 pixel=Sc;
                 break;
               }
-            CompositeHSB(destination_pixel.red,destination_pixel.green,
-              destination_pixel.blue,&hue,&saturation,&brightness);
-            CompositeHSB(source_pixel.red,source_pixel.green,source_pixel.blue,
-              &sans,&saturation,&sans);
-            HSBComposite(hue,saturation,brightness,&red,&green,&blue);
+            CompositeHCL(destination_pixel.red,destination_pixel.green,
+              destination_pixel.blue,&hue,&chroma,&luma);
+            CompositeHCL(source_pixel.red,source_pixel.green,source_pixel.blue,
+              &sans,&chroma,&sans);
+            HCLComposite(hue,chroma,luma,&red,&green,&blue);
             switch (channel)
             {
               case RedPixelChannel: pixel=red; break;
@@ -2239,11 +2260,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           }
           case ThresholdCompositeOp:
           {
-            MagickRealType
+            double
               delta;
 
             delta=Sc-Dc;
-            if ((MagickRealType) fabs((double) (2.0*delta)) < threshold)
+            if ((double) fabs((double) (2.0*delta)) < threshold)
               {
                 pixel=gamma*Dc;
                 break;
@@ -2315,6 +2336,8 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
     destination_image=DestroyImage(destination_image);
   else
     composite_image=DestroyImage(composite_image);
+  if (status != MagickFalse)
+    (void) ClampImage(image,exception);
   return(status);
 }
 \f
@@ -2378,8 +2401,8 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
     exception);
   status=MagickTrue;
   if ((image->compose != CopyCompositeOp) &&
-      ((image->compose != OverCompositeOp) || (image->matte != MagickFalse) ||
-       (texture_image->matte != MagickFalse)))
+      ((image->compose != OverCompositeOp) || (image->alpha_trait == BlendPixelTrait) ||
+       (texture_image->alpha_trait == BlendPixelTrait)))
     {
       /*
         Tile texture onto the image background.
@@ -2496,9 +2519,9 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
             texture_traits,
             traits;
 
-          channel=GetPixelChannelMapChannel(texture_image,i);
-          texture_traits=GetPixelChannelMapTraits(texture_image,channel);
-          traits=GetPixelChannelMapTraits(image,channel);
+          channel=GetPixelChannelChannel(texture_image,i);
+          texture_traits=GetPixelChannelTraits(texture_image,channel);
+          traits=GetPixelChannelTraits(image,channel);
           if ((traits == UndefinedPixelTrait) ||
               (texture_traits == UndefinedPixelTrait))
             continue;