granicus.if.org Git - imagemagick/blob - filters/analyze.c

   1 /*
   2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   3 %                                                                             %
   4 %                                                                             %
   5 %                AAA   N   N   AAA   L      Y   Y  ZZZZZ  EEEEE               %
   6 %               A   A  NN  N  A   A  L       Y Y      ZZ  E                   %
   7 %               AAAAA  N N N  AAAAA  L        Y     ZZZ   EEE                 %
   8 %               A   A  N  NN  A   A  L        Y    ZZ     E                   %
   9 %               A   A  N   N  A   A  LLLLL    Y    ZZZZZ  EEEEE               %
  10 %                                                                             %
  11 %                             Analyze An Image                                %
  12 %                                                                             %
  13 %                             Software Design                                 %
  14 %                               Bill Corbis                                   %
  15 %                              December 1998                                  %
  16 %                                                                             %
  17 %                                                                             %
  18 %  Copyright 1999-2012 ImageMagick Studio LLC, a non-profit organization      %
  19 %  dedicated to making software imaging solutions freely available.           %
  20 %                                                                             %
  21 %  You may not use this file except in compliance with the License.  You may  %
  22 %  obtain a copy of the License at                                            %
  23 %                                                                             %
  24 %    http://www.imagemagick.org/script/license.php                            %
  25 %                                                                             %
  26 %  Unless required by applicable law or agreed to in writing, software        %
  27 %  distributed under the License is distributed on an "AS IS" BASIS,          %
  28 %  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.   %
  29 %  See the License for the specific language governing permissions and        %
  30 %  limitations under the License.                                             %
  31 %                                                                             %
  32 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  33 %
  34 */
  35 \f
  36 /*
  37   Include declarations.
  38 */
  39 #include <stdio.h>
  40 #include <stdlib.h>
  41 #include <string.h>
  42 #include <time.h>
  43 #include <assert.h>
  44 #include <math.h>
  45 #include "MagickCore/MagickCore.h"
  46 \f
  47 /*
  48 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  49 %                                                                             %
  50 %                                                                             %
  51 %                                                                             %
  52 %   a n a l y z e I m a g e                                                   %
  53 %                                                                             %
  54 %                                                                             %
  55 %                                                                             %
  56 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  57 %
  58 %  analyzeImage() computes the brightness and saturation mean,  standard
  59 %  deviation, kurtosis and skewness and stores these values as attributes
  60 %  of the image.
  61 %
  62 %  The format of the analyzeImage method is:
  63 %
  64 %      size_t analyzeImage(Image *images,const int argc,
  65 %        char **argv,ExceptionInfo *exception)
  66 %
  67 %  A description of each parameter follows:
  68 %
  69 %    o image: the address of a structure of type Image.
  70 %
  71 %    o argc: Specifies a pointer to an integer describing the number of
  72 %      elements in the argument vector.
  73 %
  74 %    o argv: Specifies a pointer to a text array containing the command line
  75 %      arguments.
  76 %
  77 %    o exception: return any errors or warnings in this structure.
  78 %
  79 */
  80 ModuleExport size_t analyzeImage(Image **images,const int argc,
  81   const char **argv,ExceptionInfo *exception)
  82 {
  83   char
  84     text[MaxTextExtent];
  85
  86   double
  87     area,
  88     brightness,
  89     brightness_mean,
  90     brightness_standard_deviation,
  91     brightness_kurtosis,
  92     brightness_skewness,
  93     brightness_sum_x,
  94     brightness_sum_x2,
  95     brightness_sum_x3,
  96     brightness_sum_x4,
  97     hue,
  98     saturation,
  99     saturation_mean,
 100     saturation_standard_deviation,
 101     saturation_kurtosis,
 102     saturation_skewness,
 103     saturation_sum_x,
 104     saturation_sum_x2,
 105     saturation_sum_x3,
 106     saturation_sum_x4;
 107
 108   Image
 109     *image;
 110
 111   assert(images != (Image **) NULL);
 112   assert(*images != (Image *) NULL);
 113   assert((*images)->signature == MagickSignature);
 114   (void) argc;
 115   (void) argv;
 116   image=(*images);
 117   for ( ; image != (Image *) NULL; image=GetNextImageInList(image))
 118   {
 119     CacheView
 120       *image_view;
 121
 122     ssize_t
 123       y;
 124
 125     MagickBooleanType
 126       status;
 127
 128     brightness_sum_x=0.0;
 129     brightness_sum_x2=0.0;
 130     brightness_sum_x3=0.0;
 131     brightness_sum_x4=0.0;
 132     brightness_mean=0.0;
 133     brightness_standard_deviation=0.0;
 134     brightness_kurtosis=0.0;
 135     brightness_skewness=0.0;
 136     saturation_sum_x=0.0;
 137     saturation_sum_x2=0.0;
 138     saturation_sum_x3=0.0;
 139     saturation_sum_x4=0.0;
 140     saturation_mean=0.0;
 141     saturation_standard_deviation=0.0;
 142     saturation_kurtosis=0.0;
 143     saturation_skewness=0.0;
 144     area=0.0;
 145     status=MagickTrue;
 146     image_view=AcquireCacheView(image);
 147 #if defined(MAGICKCORE_OPENMP_SUPPORT)
 148     #pragma omp parallel for schedule(static,4) shared(status)
 149 #endif
 150     for (y=0; y < (ssize_t) image->rows; y++)
 151     {
 152       register const Quantum
 153         *p;
 154
 155       register ssize_t
 156         x;
 157
 158       if (status == MagickFalse)
 159         continue;
 160       p=GetCacheViewVirtualPixels(image_view,0,y,image->columns,1,exception);
 161       if (p == (const Quantum *) NULL)
 162         {
 163           status=MagickFalse;
 164           continue;
 165         }
 166       for (x=0; x < (ssize_t) image->columns; x++)
 167       {
 168         ConvertRGBToHSL(GetPixelRed(image,p),GetPixelGreen(image,p),
 169           GetPixelBlue(image,p),&hue,&saturation,&brightness);
 170         brightness*=QuantumRange;
 171         brightness_sum_x+=brightness;
 172         brightness_sum_x2+=brightness*brightness;
 173         brightness_sum_x3+=brightness*brightness*brightness;
 174         brightness_sum_x4+=brightness*brightness*brightness*brightness;
 175         saturation*=QuantumRange;
 176         saturation_sum_x+=saturation;
 177         saturation_sum_x2+=saturation*saturation;
 178         saturation_sum_x3+=saturation*saturation*saturation;
 179         saturation_sum_x4+=saturation*saturation*saturation*saturation;
 180         area++;
 181         p+=GetPixelChannels(image);
 182       }
 183     }
 184     image_view=DestroyCacheView(image_view);
 185     if (area <= 0.0)
 186       break;
 187     brightness_mean=brightness_sum_x/area;
 188     (void) FormatLocaleString(text,MaxTextExtent,"%g",brightness_mean);
 189     (void) SetImageProperty(image,"filter:brightness:mean",text,
 190       exception);
 191     brightness_standard_deviation=sqrt(brightness_sum_x2/area-(brightness_sum_x/
 192       area*brightness_sum_x/area));
 193     (void) FormatLocaleString(text,MaxTextExtent,"%g",
 194       brightness_standard_deviation);
 195     (void) SetImageProperty(image,"filter:brightness:standard-deviation",text,
 196       exception);
 197     if (fabs(brightness_standard_deviation) >= MagickEpsilon)
 198       brightness_kurtosis=(brightness_sum_x4/area-4.0*brightness_mean*
 199         brightness_sum_x3/area+6.0*brightness_mean*brightness_mean*
 200         brightness_sum_x2/area-3.0*brightness_mean*brightness_mean*
 201         brightness_mean*brightness_mean)/(brightness_standard_deviation*
 202         brightness_standard_deviation*brightness_standard_deviation*
 203         brightness_standard_deviation)-3.0;
 204     (void) FormatLocaleString(text,MaxTextExtent,"%g",brightness_kurtosis);
 205     (void) SetImageProperty(image,"filter:brightness:kurtosis",text,
 206       exception);
 207     if (brightness_standard_deviation != 0)
 208       brightness_skewness=(brightness_sum_x3/area-3.0*brightness_mean*
 209         brightness_sum_x2/area+2.0*brightness_mean*brightness_mean*
 210         brightness_mean)/(brightness_standard_deviation*
 211         brightness_standard_deviation*brightness_standard_deviation);
 212     (void) FormatLocaleString(text,MaxTextExtent,"%g",brightness_skewness);
 213     (void) SetImageProperty(image,"filter:brightness:skewness",text,
 214       exception);
 215     saturation_mean=saturation_sum_x/area;
 216     (void) FormatLocaleString(text,MaxTextExtent,"%g",saturation_mean);
 217     (void) SetImageProperty(image,"filter:saturation:mean",text,
 218       exception);
 219     saturation_standard_deviation=sqrt(saturation_sum_x2/area-(saturation_sum_x/
 220       area*saturation_sum_x/area));
 221     (void) FormatLocaleString(text,MaxTextExtent,"%g",
 222       saturation_standard_deviation);
 223     (void) SetImageProperty(image,"filter:saturation:standard-deviation",text,
 224       exception);
 225     if (fabs(saturation_standard_deviation) >= MagickEpsilon)
 226       saturation_kurtosis=(saturation_sum_x4/area-4.0*saturation_mean*
 227         saturation_sum_x3/area+6.0*saturation_mean*saturation_mean*
 228         saturation_sum_x2/area-3.0*saturation_mean*saturation_mean*
 229         saturation_mean*saturation_mean)/(saturation_standard_deviation*
 230         saturation_standard_deviation*saturation_standard_deviation*
 231         saturation_standard_deviation)-3.0;
 232     (void) FormatLocaleString(text,MaxTextExtent,"%g",saturation_kurtosis);
 233     (void) SetImageProperty(image,"filter:saturation:kurtosis",text,
 234       exception);
 235     if (fabs(saturation_standard_deviation) >= MagickEpsilon)
 236       saturation_skewness=(saturation_sum_x3/area-3.0*saturation_mean*
 237         saturation_sum_x2/area+2.0*saturation_mean*saturation_mean*
 238         saturation_mean)/(saturation_standard_deviation*
 239         saturation_standard_deviation*saturation_standard_deviation);
 240     (void) FormatLocaleString(text,MaxTextExtent,"%g",saturation_skewness);
 241     (void) SetImageProperty(image,"filter:saturation:skewness",text,
 242       exception);
 243   }
 244   return(MagickImageFilterSignature);
 245 }