GNOME.org
 

[gegl] operations: add gegl:diffraction-patterns

commit 7589e6107db6b94f92213e0faedeb9bca146a5db
Author: Thomas Manni <thomas manni free fr>
Date:   Wed Jan 14 16:03:10 2015 +0100

    operations: add gegl:diffraction-patterns

 operations/common/Makefile.am                      |    1 +
 operations/common/diffraction-patterns.c           |  259 ++++++++++++++++++++
 po/POTFILES.in                                     |    1 +
 tests/compositions/Makefile.am                     |    1 +
 tests/compositions/diffraction-patterns.xml        |   17 ++
 .../reference/diffraction-patterns.png             |  Bin 0 -> 212438 bytes
 6 files changed, 279 insertions(+), 0 deletions(-)
---
diff --git a/operations/common/Makefile.am b/operations/common/Makefile.am
index 2833161..1663c89 100644
--- a/operations/common/Makefile.am
+++ b/operations/common/Makefile.am
@@ -37,6 +37,7 @@ op_LTLIBRARIES = \
        cubism.la \
        deinterlace.la \
        difference-of-gaussians.la \
+       diffraction-patterns.la \
        display.la \
        distance-transform.la \
        dropshadow.la \
diff --git a/operations/common/diffraction-patterns.c b/operations/common/diffraction-patterns.c
new file mode 100644
index 0000000..5b2d935
--- /dev/null
+++ b/operations/common/diffraction-patterns.c
@@ -0,0 +1,259 @@
+/* This file is an image processing operation for GEGL
+ *
+ * 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 <http://www.gnu.org/licenses/>.
+ *
+ * Copyright (C) 1997 Federico Mena Quintero and David Bleecker
+ * federico nuclecu unam mx
+ * bleecker math hawaii edu
+ *
+ * GEGL port: Thomas Manni <thomas manni free fr>
+ *
+ */
+
+#include "config.h"
+#include <glib/gi18n-lib.h>
+#include <math.h>
+
+#ifdef GEGL_PROPERTIES
+
+property_double (red_frequency, _("Red frenquency"), 0.815)
+    description (_("Light frequency (red)"))
+    value_range (0.0, 20.0)
+
+property_double (green_frequency, _("Green frenquency"), 1.221)
+    description (_("Light frequency (green)"))
+    value_range (0.0, 20.0)
+
+property_double (blue_frequency, _("Blue frenquency"), 1.123)
+    description (_("Light frequency (blue)"))
+    value_range (0.0, 20.0)
+
+property_double (red_contours, _("Red contours"), 0.821)
+    description (_("Number of contours (red)"))
+    value_range (0.0, 10.0)
+
+property_double (green_contours, _("Green contours"), 0.821)
+    description (_("Number of contours (green)"))
+    value_range (0.0, 10.0)
+
+property_double (blue_contours, _("Blue contours"), 0.974)
+    description (_("Number of contours (blue)"))
+    value_range (0.0, 10.0)
+
+property_double (red_sedges, _("Red sharp edges"), 0.610)
+    description (_("Number of sharp edges (red)"))
+    value_range (0.0, 1.0)
+
+property_double (green_sedges, _("Green sharp edges"), 0.677)
+    description (_("Number of sharp edges (green)"))
+    value_range (0.0, 1.0)
+
+property_double (blue_sedges, _("Blue sharp edges"), 0.636)
+    description (_("Number of sharp edges (blue)"))
+    value_range (0.0, 1.0)
+
+property_double (brightness, _("Brightness"), 0.066)
+    description (_("Brightness and shifting/fattening of contours"))
+    value_range (0.0, 1.0)
+
+property_double (scattering, _("Scattering"), 37.126)
+    description (_("Scattering (Speed vs. quality)"))
+    value_range (0.0, 100.0)
+
+property_double (polarization, _("Polarization"), -0.473)
+    description (_("Polarization"))
+    value_range (-1.0, 1.0)
+
+property_int    (width, _("Width"), 200)
+    description (_("Width of the generated buffer"))
+    value_range (0, G_MAXINT)
+    ui_range    (0, 4096)
+    ui_meta     ("unit", "pixel-distance")
+    ui_meta     ("axis", "x")
+    ui_meta     ("role", "output-extent")
+
+property_int    (height, _("Height"), 200)
+    description (_("Height of the generated buffer"))
+    value_range (0, G_MAXINT)
+    ui_range    (0, 4096)
+    ui_meta     ("unit", "pixel-distance")
+    ui_meta     ("axis", "y")
+    ui_meta     ("role", "output-extent")
+
+#else
+
+#define GEGL_OP_POINT_RENDER
+#define GEGL_OP_C_FILE "diffraction-patterns.c"
+
+#include "gegl-op.h"
+
+#define ITERATIONS      100
+#define WEIRD_FACTOR    0.04
+
+static gdouble cos_lut[ITERATIONS + 1];
+static gdouble lut1[ITERATIONS + 1];
+static gdouble lut2[ITERATIONS + 1];
+
+static void
+init_luts (void)
+{
+  gint    i;
+  gdouble a;
+  gdouble sina;
+
+  a = -G_PI;
+
+  for (i = 0; i <= ITERATIONS; i++)
+    {
+      sina = sin (a);
+      cos_lut[i] = cos (a);
+
+      lut1[i] = 0.75 * sina;
+      lut2[i] = 0.5 * (4.0 * cos_lut[i] * cos_lut[i] + sina * sina);
+
+      a += (2.0 * G_PI / ITERATIONS);
+    }
+}
+
+static inline gdouble
+diff_intensity (gdouble         x,
+                gdouble         y,
+                gdouble         lam,
+                GeglProperties *o)
+{
+  gint    i;
+  gdouble cxy, sxy;
+  gdouble p;
+  gdouble polpi2;
+  gdouble cospolpi2, sinpolpi2;
+
+  cxy = 0.0;
+  sxy = 0.0;
+
+  lam *= 4.0;
+
+  for (i = 0; i <= ITERATIONS; i++)
+    {
+      p = lam * (cos_lut[i] * x + lut1[i] * y - lut2[i]);
+
+      cxy += cos (p);
+      sxy += sin (p);
+    }
+
+  cxy *= WEIRD_FACTOR;
+  sxy *= WEIRD_FACTOR;
+
+  polpi2 = o->polarization * (G_PI / 2.0);
+
+  cospolpi2 = cos (polpi2);
+  sinpolpi2 = sin (polpi2);
+
+  return o->scattering * ((cospolpi2 + sinpolpi2) * cxy * cxy +
+                          (cospolpi2 - sinpolpi2) * sxy * sxy);
+}
+
+static void
+diffract (gdouble         x,
+          gdouble         y,
+          gdouble        *r,
+          gdouble        *g,
+          gdouble        *b,
+          GeglProperties *o)
+{
+  *r = fabs (o->red_sedges * sin (o->red_contours * atan (o->brightness *
+                                             diff_intensity (x, y, o->red_frequency, o))));
+  *g = fabs (o->green_sedges * sin (o->green_contours * atan (o->brightness *
+                                             diff_intensity (x, y, o->green_frequency, o))));
+  *b = fabs (o->blue_sedges * sin (o->blue_contours * atan (o->brightness *
+                                             diff_intensity (x, y, o->blue_frequency, o))));
+}
+
+static void
+prepare (GeglOperation *operation)
+{
+  gegl_operation_set_format (operation, "output", babl_format ("R'G'B' float"));
+}
+
+static GeglRectangle
+get_bounding_box (GeglOperation *operation)
+{
+  return gegl_rectangle_infinite_plane ();
+}
+
+static gboolean
+process (GeglOperation       *operation,
+         void                *out_buf,
+         glong                n_pixels,
+         const GeglRectangle *roi,
+         gint                 level)
+{
+  GeglProperties *o = GEGL_PROPERTIES (operation);
+  gfloat         *out;
+  gdouble         r, g, b;
+  gint            x, y;
+  gdouble         px, py;
+  gdouble         dh, dv;
+
+  dh = 10.0 / (o->width - 1);
+  dv = -10.0 / (o->height - 1);
+  out = out_buf;
+
+  for (y = roi->y; y < roi->y + roi->height; y++)
+    {
+      for (x = roi->x; x < roi->x + roi->width; x++)
+        {
+          px = dh * x - 5.0;
+          py = dv * y + 5.0;
+
+          diffract (px, py, &r, &g, &b, o);
+
+          out[0] = r;
+          out[1] = g;
+          out[2] = b;
+
+          out += 3;
+        }
+    }
+
+  return TRUE;
+}
+
+static void
+gegl_op_class_init (GeglOpClass *klass)
+{
+  GeglOperationClass            *operation_class;
+  GeglOperationPointRenderClass *point_render_class;
+
+  init_luts();
+
+  operation_class    = GEGL_OPERATION_CLASS (klass);
+  point_render_class = GEGL_OPERATION_POINT_RENDER_CLASS (klass);
+
+  point_render_class->process = process;
+
+  operation_class->get_bounding_box = get_bounding_box;
+  operation_class->prepare = prepare;
+  operation_class->opencl_support = FALSE;
+
+  gegl_operation_class_set_keys (operation_class,
+    "name",               "gegl:diffraction-patterns",
+    "title",              _("Diffraction Patterns"),
+    "categories",         "render",
+    "position-dependent", "true",
+    "license",            "GPL3+",
+    "description",        _("Generate diffraction patterns"),
+    NULL);
+}
+
+#endif
diff --git a/po/POTFILES.in b/po/POTFILES.in
index 3aa29ae..d187d4e 100644
--- a/po/POTFILES.in
+++ b/po/POTFILES.in
@@ -31,6 +31,7 @@ operations/common/copy-buffer.c
 operations/common/cubism.c
 operations/common/deinterlace.c
 operations/common/difference-of-gaussians.c
+operations/common/diffraction-patterns.c
 operations/common/display.c
 operations/common/distance-transform.c
 operations/common/dropshadow.c
diff --git a/tests/compositions/Makefile.am b/tests/compositions/Makefile.am
index d7d65b9..cfd84f5 100644
--- a/tests/compositions/Makefile.am
+++ b/tests/compositions/Makefile.am
@@ -26,6 +26,7 @@ TESTS = \
   color-to-alpha.xml              \
   composite-transform.xml         \
   contrast-curve.xml              \
+  diffraction-patterns.xml        \
   edge-laplace.xml                \
   edge-sobel.xml                  \
   engrave.xml                     \
diff --git a/tests/compositions/diffraction-patterns.xml b/tests/compositions/diffraction-patterns.xml
new file mode 100644
index 0000000..c99ccd2
--- /dev/null
+++ b/tests/compositions/diffraction-patterns.xml
@@ -0,0 +1,17 @@
+<?xml version='1.0' encoding='UTF-8'?>
+<gegl>
+  <node operation='gegl:crop'>
+      <params>
+        <param name='x'>0</param>
+        <param name='y'>0</param>
+        <param name='width'>200</param>
+        <param name='height'>200</param>
+      </params>
+  </node>
+  <node operation='gegl:diffraction-patterns'>
+    <params>
+        <param name="width">200</param>
+        <param name="height">200</param>
+      </params>
+  </node>
+</gegl>
diff --git a/tests/compositions/reference/diffraction-patterns.png 
b/tests/compositions/reference/diffraction-patterns.png
new file mode 100644
index 0000000..0e499ae
Binary files /dev/null and b/tests/compositions/reference/diffraction-patterns.png differ
[Date Prev][Date Next]   [Thread Prev][Thread Next]   [Thread Index] [Date Index] [Author Index]