gimpoperationlevels.c 6.47 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
/* GIMP - The GNU Image Manipulation Program
 * Copyright (C) 1995 Spencer Kimball and Peter Mattis
 *
 * gimpoperationlevels.c
 * Copyright (C) 2007 Michael Natterer <mitch@gimp.org>
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
 */

#include "config.h"

24
#include <gegl.h>
25 26 27 28 29 30

#include "libgimpcolor/gimpcolor.h"
#include "libgimpmath/gimpmath.h"

#include "gegl-types.h"

31
#include "gimplevelsconfig.h"
32 33 34
#include "gimpoperationlevels.h"


35 36 37 38 39
static gboolean gimp_operation_levels_process (GeglOperation       *operation,
                                               void                *in_buf,
                                               void                *out_buf,
                                               glong                samples,
                                               const GeglRectangle *roi);
40 41 42


G_DEFINE_TYPE (GimpOperationLevels, gimp_operation_levels,
43
               GIMP_TYPE_OPERATION_POINT_FILTER)
44 45 46 47 48

#define parent_class gimp_operation_levels_parent_class


static void
49
gimp_operation_levels_class_init (GimpOperationLevelsClass *klass)
50 51 52 53 54
{
  GObjectClass                  *object_class    = G_OBJECT_CLASS (klass);
  GeglOperationClass            *operation_class = GEGL_OPERATION_CLASS (klass);
  GeglOperationPointFilterClass *point_class     = GEGL_OPERATION_POINT_FILTER_CLASS (klass);

55 56
  object_class->set_property   = gimp_operation_point_filter_set_property;
  object_class->get_property   = gimp_operation_point_filter_get_property;
57

58
  operation_class->name        = "gimp:levels";
59 60
  operation_class->categories  = "color";
  operation_class->description = "GIMP Levels operation";
61

62
  point_class->process         = gimp_operation_levels_process;
63

64 65
  g_object_class_install_property (object_class,
                                   GIMP_OPERATION_POINT_FILTER_PROP_CONFIG,
66 67 68 69
                                   g_param_spec_object ("config",
                                                        "Config",
                                                        "The config object",
                                                        GIMP_TYPE_LEVELS_CONFIG,
70 71
                                                        G_PARAM_READWRITE |
                                                        G_PARAM_CONSTRUCT));
72 73 74 75 76
}

static void
gimp_operation_levels_init (GimpOperationLevels *self)
{
77
}
78

79 80
static inline gdouble
gimp_operation_levels_map (gdouble value,
81
                           gdouble inv_gamma,
82 83 84 85
                           gdouble low_input,
                           gdouble high_input,
                           gdouble low_output,
                           gdouble high_output)
86 87 88
{
  /*  determine input intensity  */
  if (high_input != low_input)
89
    value = (value - low_input) / (high_input - low_input);
90
  else
91
    value = (value - low_input);
92

93
  if (inv_gamma != 1.0)
94 95
    {
      if (value >= 0.0)
96
        value =  pow ( value, inv_gamma);
97
      else
98
        value = -pow (-value, inv_gamma);
99 100 101 102 103 104 105 106 107 108 109 110
    }

  /*  determine the output intensity  */
  if (high_output >= low_output)
    value = value * (high_output - low_output) + low_output;
  else if (high_output < low_output)
    value = low_output - value * (low_output - high_output);

  return value;
}

static gboolean
111 112 113 114 115
gimp_operation_levels_process (GeglOperation       *operation,
                               void                *in_buf,
                               void                *out_buf,
                               glong                samples,
                               const GeglRectangle *roi)
116
{
117 118 119 120
  GimpOperationPointFilter *point  = GIMP_OPERATION_POINT_FILTER (operation);
  GimpLevelsConfig         *config = GIMP_LEVELS_CONFIG (point->config);
  gfloat                   *src    = in_buf;
  gfloat                   *dest   = out_buf;
121 122
  gfloat                    inv_gamma[5];
  gint                      channel;
123

124 125 126
  if (! config)
    return FALSE;

127
  for (channel = 0; channel < 5; channel++)
128
    {
129 130 131 132
      g_return_val_if_fail (config->gamma[channel] != 0.0, FALSE);

      inv_gamma[channel] = 1.0 / config->gamma[channel];
    }
133

134 135
  while (samples--)
    {
136 137
      for (channel = 0; channel < 4; channel++)
        {
138
          gdouble value;
139 140

          value = gimp_operation_levels_map (src[channel],
141
                                             inv_gamma[channel + 1],
142 143 144 145
                                             config->low_input[channel + 1],
                                             config->high_input[channel + 1],
                                             config->low_output[channel + 1],
                                             config->high_output[channel + 1]);
146 147

          /* don't apply the overall curve to the alpha channel */
148
          if (channel != ALPHA)
149
            value = gimp_operation_levels_map (value,
150
                                               inv_gamma[0],
151 152 153 154
                                               config->low_input[0],
                                               config->high_input[0],
                                               config->low_output[0],
                                               config->high_output[0]);
155 156 157 158 159 160 161 162 163 164

          dest[channel] = value;
        }

      src  += 4;
      dest += 4;
    }

  return TRUE;
}
165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191


/*  public functions  */

gdouble
gimp_operation_levels_map_input (GimpLevelsConfig     *config,
                                 GimpHistogramChannel  channel,
                                 gdouble               value)
{
  g_return_val_if_fail (GIMP_IS_LEVELS_CONFIG (config), 0.0);

  /*  determine input intensity  */
  if (config->high_input[channel] != config->low_input[channel])
    value = ((value - config->low_input[channel]) /
             (config->high_input[channel] - config->low_input[channel]));
  else
    value = (value - config->low_input[channel]);

  value = CLAMP (value, 0.0, 1.0);

  if (config->gamma[channel] != 0.0)
    {
      value = pow (value, 1.0 / config->gamma[channel]);
    }

  return value;
}