summaryrefslogtreecommitdiffstats
path: root/app/operations/gimpoperationcagecoefcalc.c
diff options
context:
space:
mode:
Diffstat (limited to 'app/operations/gimpoperationcagecoefcalc.c')
-rw-r--r--app/operations/gimpoperationcagecoefcalc.c294
1 files changed, 294 insertions, 0 deletions
diff --git a/app/operations/gimpoperationcagecoefcalc.c b/app/operations/gimpoperationcagecoefcalc.c
new file mode 100644
index 0000000..ba8cd3f
--- /dev/null
+++ b/app/operations/gimpoperationcagecoefcalc.c
@@ -0,0 +1,294 @@
+/* GIMP - The GNU Image Manipulation Program
+ *
+ * gimpoperationcagecoefcalc.c
+ * Copyright (C) 2010 Michael Muré <batolettre@gmail.com>
+ *
+ * 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 3 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, see <https://www.gnu.org/licenses/>.
+ */
+
+#include "config.h"
+
+#include <gdk-pixbuf/gdk-pixbuf.h>
+#include <gegl.h>
+
+#include "libgimpmath/gimpmath.h"
+
+#include "operations-types.h"
+
+#include "gimpoperationcagecoefcalc.h"
+#include "gimpcageconfig.h"
+
+#include "gimp-intl.h"
+
+
+static void gimp_operation_cage_coef_calc_finalize (GObject *object);
+static void gimp_operation_cage_coef_calc_get_property (GObject *object,
+ guint property_id,
+ GValue *value,
+ GParamSpec *pspec);
+static void gimp_operation_cage_coef_calc_set_property (GObject *object,
+ guint property_id,
+ const GValue *value,
+ GParamSpec *pspec);
+
+static void gimp_operation_cage_coef_calc_prepare (GeglOperation *operation);
+static GeglRectangle gimp_operation_cage_coef_calc_get_bounding_box (GeglOperation *operation);
+static gboolean gimp_operation_cage_coef_calc_process (GeglOperation *operation,
+ GeglBuffer *output,
+ const GeglRectangle *roi,
+ gint level);
+
+
+G_DEFINE_TYPE (GimpOperationCageCoefCalc, gimp_operation_cage_coef_calc,
+ GEGL_TYPE_OPERATION_SOURCE)
+
+#define parent_class gimp_operation_cage_coef_calc_parent_class
+
+
+static void
+gimp_operation_cage_coef_calc_class_init (GimpOperationCageCoefCalcClass *klass)
+{
+ GObjectClass *object_class = G_OBJECT_CLASS (klass);
+ GeglOperationSourceClass *source_class = GEGL_OPERATION_SOURCE_CLASS (klass);
+ GeglOperationClass *operation_class = GEGL_OPERATION_CLASS (klass);
+
+ gegl_operation_class_set_keys (operation_class,
+ "name", "gimp:cage-coef-calc",
+ "categories", "transform",
+ "description", _("Compute a set of coefficient buffer for the GIMP cage tool"),
+ NULL);
+
+ operation_class->prepare = gimp_operation_cage_coef_calc_prepare;
+ operation_class->get_bounding_box = gimp_operation_cage_coef_calc_get_bounding_box;
+ operation_class->cache_policy = GEGL_CACHE_POLICY_ALWAYS;
+ operation_class->get_cached_region = NULL;
+
+ source_class->process = gimp_operation_cage_coef_calc_process;
+
+ object_class->get_property = gimp_operation_cage_coef_calc_get_property;
+ object_class->set_property = gimp_operation_cage_coef_calc_set_property;
+ object_class->finalize = gimp_operation_cage_coef_calc_finalize;
+
+ g_object_class_install_property (object_class,
+ GIMP_OPERATION_CAGE_COEF_CALC_PROP_CONFIG,
+ g_param_spec_object ("config",
+ "Config",
+ "A GimpCageConfig object, that define the transformation",
+ GIMP_TYPE_CAGE_CONFIG,
+ G_PARAM_READWRITE |
+ G_PARAM_CONSTRUCT));
+}
+
+static void
+gimp_operation_cage_coef_calc_init (GimpOperationCageCoefCalc *self)
+{
+}
+
+static void
+gimp_operation_cage_coef_calc_finalize (GObject *object)
+{
+ GimpOperationCageCoefCalc *self = GIMP_OPERATION_CAGE_COEF_CALC (object);
+
+ g_clear_object (&self->config);
+
+ G_OBJECT_CLASS (parent_class)->finalize (object);
+}
+
+static void
+gimp_operation_cage_coef_calc_get_property (GObject *object,
+ guint property_id,
+ GValue *value,
+ GParamSpec *pspec)
+{
+ GimpOperationCageCoefCalc *self = GIMP_OPERATION_CAGE_COEF_CALC (object);
+
+ switch (property_id)
+ {
+ case GIMP_OPERATION_CAGE_COEF_CALC_PROP_CONFIG:
+ g_value_set_object (value, self->config);
+ break;
+
+ default:
+ G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
+ break;
+ }
+}
+
+static void
+gimp_operation_cage_coef_calc_set_property (GObject *object,
+ guint property_id,
+ const GValue *value,
+ GParamSpec *pspec)
+{
+ GimpOperationCageCoefCalc *self = GIMP_OPERATION_CAGE_COEF_CALC (object);
+
+ switch (property_id)
+ {
+ case GIMP_OPERATION_CAGE_COEF_CALC_PROP_CONFIG:
+ if (self->config)
+ g_object_unref (self->config);
+ self->config = g_value_dup_object (value);
+ break;
+
+ default:
+ G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
+ break;
+ }
+}
+
+static gboolean
+gimp_operation_cage_coef_calc_is_on_straight (GimpVector2 *d1,
+ GimpVector2 *d2,
+ GimpVector2 *p)
+{
+ GimpVector2 v1, v2;
+ gfloat deter;
+
+ v1.x = p->x - d1->x;
+ v1.y = p->y - d1->y;
+ v2.x = d2->x - d1->x;
+ v2.y = d2->y - d1->y;
+
+ gimp_vector2_normalize (&v1);
+ gimp_vector2_normalize (&v2);
+
+ deter = v1.x * v2.y - v2.x * v1.y;
+
+ return (deter < 0.000000001) && (deter > -0.000000001);
+}
+
+static void
+gimp_operation_cage_coef_calc_prepare (GeglOperation *operation)
+{
+ GimpOperationCageCoefCalc *occc = GIMP_OPERATION_CAGE_COEF_CALC (operation);
+ GimpCageConfig *config = GIMP_CAGE_CONFIG (occc->config);
+
+ gegl_operation_set_format (operation,
+ "output",
+ babl_format_n (babl_type ("float"),
+ 2 * gimp_cage_config_get_n_points (config)));
+}
+
+static GeglRectangle
+gimp_operation_cage_coef_calc_get_bounding_box (GeglOperation *operation)
+{
+ GimpOperationCageCoefCalc *occc = GIMP_OPERATION_CAGE_COEF_CALC (operation);
+ GimpCageConfig *config = GIMP_CAGE_CONFIG (occc->config);
+
+ return gimp_cage_config_get_bounding_box (config);
+}
+
+static gboolean
+gimp_operation_cage_coef_calc_process (GeglOperation *operation,
+ GeglBuffer *output,
+ const GeglRectangle *roi,
+ gint level)
+{
+ GimpOperationCageCoefCalc *occc = GIMP_OPERATION_CAGE_COEF_CALC (operation);
+ GimpCageConfig *config = GIMP_CAGE_CONFIG (occc->config);
+
+ const Babl *format;
+
+ GeglBufferIterator *it;
+ guint n_cage_vertices;
+ GimpCagePoint *current, *last;
+
+ if (! config)
+ return FALSE;
+
+ format = babl_format_n (babl_type ("float"), 2 * gimp_cage_config_get_n_points (config));
+
+ n_cage_vertices = gimp_cage_config_get_n_points (config);
+
+ it = gegl_buffer_iterator_new (output, roi, 0, format,
+ GEGL_ACCESS_READWRITE, GEGL_ABYSS_NONE, 1);
+
+ while (gegl_buffer_iterator_next (it))
+ {
+ /* iterate inside the roi */
+ gint n_pixels = it->length;
+ gint x = it->items[0].roi.x; /* initial x */
+ gint y = it->items[0].roi.y; /* and y coordinates */
+ gint j;
+
+ gfloat *coef = it->items[0].data;
+
+ while(n_pixels--)
+ {
+ if (gimp_cage_config_point_inside(config, x, y))
+ {
+ last = &(g_array_index (config->cage_points, GimpCagePoint, 0));
+
+ for( j = 0; j < n_cage_vertices; j++)
+ {
+ GimpVector2 v1,v2,a,b,p;
+ gdouble BA,SRT,L0,L1,A0,A1,A10,L10, Q,S,R, absa;
+
+ current = &(g_array_index (config->cage_points, GimpCagePoint, (j+1) % n_cage_vertices));
+ v1 = last->src_point;
+ v2 = current->src_point;
+ p.x = x;
+ p.y = y;
+ a.x = v2.x - v1.x;
+ a.y = v2.y - v1.y;
+ absa = gimp_vector2_length (&a);
+
+ b.x = v1.x - x;
+ b.y = v1.y - y;
+ Q = a.x * a.x + a.y * a.y;
+ S = b.x * b.x + b.y * b.y;
+ R = 2.0 * (a.x * b.x + a.y * b.y);
+ BA = b.x * a.y - b.y * a.x;
+ SRT = sqrt(4.0 * S * Q - R * R);
+
+ L0 = log(S);
+ L1 = log(S + Q + R);
+ A0 = atan2(R, SRT) / SRT;
+ A1 = atan2(2.0 * Q + R, SRT) / SRT;
+ A10 = A1 - A0;
+ L10 = L1 - L0;
+
+ /* edge coef */
+ coef[j + n_cage_vertices] = (-absa / (4.0 * G_PI)) * ((4.0*S-(R*R)/Q) * A10 + (R / (2.0 * Q)) * L10 + L1 - 2.0);
+
+ if (isnan(coef[j + n_cage_vertices]))
+ {
+ coef[j + n_cage_vertices] = 0.0;
+ }
+
+ /* vertice coef */
+ if (!gimp_operation_cage_coef_calc_is_on_straight (&v1, &v2, &p))
+ {
+ coef[j] += (BA / (2.0 * G_PI)) * (L10 /(2.0*Q) - A10 * (2.0 + R / Q));
+ coef[(j+1)%n_cage_vertices] -= (BA / (2.0 * G_PI)) * (L10 / (2.0 * Q) - A10 * (R / Q));
+ }
+
+ last = current;
+ }
+ }
+
+ coef += 2 * n_cage_vertices;
+
+ /* update x and y coordinates */
+ x++;
+ if (x >= (it->items[0].roi.x + it->items[0].roi.width))
+ {
+ x = it->items[0].roi.x;
+ y++;
+ }
+ }
+ }
+
+ return TRUE;
+}