2023-09-26 00:35:21 +02:00
|
|
|
/* LIC 0.14 -- image filter plug-in for PIKA
|
|
|
|
* Copyright (C) 1996 Tom Bech
|
|
|
|
*
|
|
|
|
* E-mail: tomb@gimp.org
|
|
|
|
* You can contact the original PIKA authors at pika@xcf.berkeley.edu
|
|
|
|
*
|
|
|
|
* 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/>.
|
|
|
|
*
|
|
|
|
* In other words, you can't sue me for whatever happens while using this ;)
|
|
|
|
*
|
|
|
|
* Changes (post 0.10):
|
|
|
|
* -> 0.11: Fixed a bug in the convolution kernels (Tom).
|
|
|
|
* -> 0.12: Added Quartic's bilinear interpolation stuff (Tom).
|
|
|
|
* -> 0.13 Changed some UI stuff causing trouble with the 0.60 release, added
|
|
|
|
* the (PIKA) tags and changed random() calls to rand() (Tom)
|
|
|
|
* -> 0.14 Ported to 0.99.11 (Tom)
|
|
|
|
*
|
|
|
|
* This plug-in implements the Line Integral Convolution (LIC) as described in
|
|
|
|
* Cabral et al. "Imaging vector fields using line integral convolution" in the
|
|
|
|
* Proceedings of ACM SIGGRAPH 93. Publ. by ACM, New York, NY, USA. p. 263-270.
|
|
|
|
* (See http://www8.cs.umu.se/kurser/TDBD13/VT00/extra/p263-cabral.pdf)
|
|
|
|
*
|
|
|
|
* Some of the code is based on code by Steinar Haugen (thanks!), the Perlin
|
|
|
|
* noise function is practically ripped as is :)
|
|
|
|
*/
|
|
|
|
|
|
|
|
#include "config.h"
|
|
|
|
|
|
|
|
#include <libpika/pika.h>
|
|
|
|
#include <libpika/pikaui.h>
|
|
|
|
|
|
|
|
#include "libpika/stdplugins-intl.h"
|
|
|
|
|
|
|
|
|
|
|
|
/************/
|
|
|
|
/* Typedefs */
|
|
|
|
/************/
|
|
|
|
|
|
|
|
#define numx 40 /* Pseudo-random vector grid size */
|
|
|
|
#define numy 40
|
|
|
|
|
|
|
|
#define PLUG_IN_PROC "plug-in-lic"
|
|
|
|
#define PLUG_IN_BINARY "van-gogh-lic"
|
|
|
|
#define PLUG_IN_ROLE "pika-van-gogh-lic"
|
|
|
|
|
|
|
|
typedef enum
|
|
|
|
{
|
|
|
|
LIC_HUE,
|
|
|
|
LIC_SATURATION,
|
|
|
|
LIC_BRIGHTNESS
|
|
|
|
} LICEffectChannel;
|
|
|
|
|
|
|
|
|
|
|
|
typedef struct _Lic Lic;
|
|
|
|
typedef struct _LicClass LicClass;
|
|
|
|
|
|
|
|
struct _Lic
|
|
|
|
{
|
|
|
|
PikaPlugIn parent_instance;
|
|
|
|
};
|
|
|
|
|
|
|
|
struct _LicClass
|
|
|
|
{
|
|
|
|
PikaPlugInClass parent_class;
|
|
|
|
};
|
|
|
|
|
|
|
|
|
|
|
|
#define LIC_TYPE (lic_get_type ())
|
|
|
|
#define LIC (obj) (G_TYPE_CHECK_INSTANCE_CAST ((obj), LIC_TYPE, Lic))
|
|
|
|
|
|
|
|
GType lic_get_type (void) G_GNUC_CONST;
|
|
|
|
|
|
|
|
static GList * lic_query_procedures (PikaPlugIn *plug_in);
|
|
|
|
static PikaProcedure * lic_create_procedure (PikaPlugIn *plug_in,
|
|
|
|
const gchar *name);
|
|
|
|
|
|
|
|
static PikaValueArray * lic_run (PikaProcedure *procedure,
|
|
|
|
PikaRunMode run_mode,
|
|
|
|
PikaImage *image,
|
|
|
|
gint n_drawables,
|
|
|
|
PikaDrawable **drawables,
|
2023-10-30 23:55:30 +01:00
|
|
|
PikaProcedureConfig *config,
|
2023-09-26 00:35:21 +02:00
|
|
|
gpointer run_data);
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
static gboolean create_main_dialog (PikaProcedure *procedure,
|
|
|
|
PikaProcedureConfig *config);
|
|
|
|
static void compute_image (PikaProcedureConfig *config,
|
|
|
|
PikaDrawable *drawable);
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
G_DEFINE_TYPE (Lic, lic, PIKA_TYPE_PLUG_IN)
|
|
|
|
|
|
|
|
PIKA_MAIN (LIC_TYPE)
|
|
|
|
DEFINE_STD_SET_I18N
|
|
|
|
|
|
|
|
|
|
|
|
static void
|
|
|
|
lic_class_init (LicClass *klass)
|
|
|
|
{
|
|
|
|
PikaPlugInClass *plug_in_class = PIKA_PLUG_IN_CLASS (klass);
|
|
|
|
|
|
|
|
plug_in_class->query_procedures = lic_query_procedures;
|
|
|
|
plug_in_class->create_procedure = lic_create_procedure;
|
|
|
|
plug_in_class->set_i18n = STD_SET_I18N;
|
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
lic_init (Lic *lic)
|
|
|
|
{
|
|
|
|
}
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
static GList *
|
|
|
|
lic_query_procedures (PikaPlugIn *plug_in)
|
|
|
|
{
|
|
|
|
return g_list_append (NULL, g_strdup (PLUG_IN_PROC));
|
|
|
|
}
|
|
|
|
|
|
|
|
static PikaProcedure *
|
|
|
|
lic_create_procedure (PikaPlugIn *plug_in,
|
|
|
|
const gchar *name)
|
|
|
|
{
|
|
|
|
PikaProcedure *procedure = NULL;
|
|
|
|
|
|
|
|
if (! strcmp (name, PLUG_IN_PROC))
|
|
|
|
{
|
|
|
|
procedure = pika_image_procedure_new (plug_in, name,
|
|
|
|
PIKA_PDB_PROC_TYPE_PLUGIN,
|
|
|
|
lic_run, NULL, NULL);
|
|
|
|
|
|
|
|
pika_procedure_set_image_types (procedure, "RGB*");
|
|
|
|
pika_procedure_set_sensitivity_mask (procedure,
|
|
|
|
PIKA_PROCEDURE_SENSITIVE_DRAWABLE);
|
|
|
|
|
|
|
|
pika_procedure_set_menu_label (procedure, _("_Van Gogh (LIC)..."));
|
|
|
|
pika_procedure_add_menu_path (procedure, "<Image>/Filters/Artistic");
|
|
|
|
|
|
|
|
pika_procedure_set_documentation (procedure,
|
|
|
|
_("Special effects that nobody "
|
|
|
|
"understands"),
|
|
|
|
"No help yet",
|
|
|
|
name);
|
|
|
|
pika_procedure_set_attribution (procedure,
|
|
|
|
"Tom Bech & Federico Mena Quintero",
|
|
|
|
"Tom Bech & Federico Mena Quintero",
|
|
|
|
"Version 0.14, September 24 1997");
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_CHOICE (procedure, "effect-channel",
|
|
|
|
_("E_ffect Channel"),
|
|
|
|
_("Effect Channel"),
|
|
|
|
pika_choice_new_with_values ("hue", LIC_HUE, _("Hue"), NULL,
|
|
|
|
"saturation", LIC_SATURATION, _("Saturation"), NULL,
|
|
|
|
"brightness", LIC_BRIGHTNESS, _("Brightness"), NULL,
|
|
|
|
NULL),
|
|
|
|
"brightness",
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_CHOICE (procedure, "effect-operator",
|
|
|
|
_("Effect O_perator"),
|
|
|
|
_("Effect Operator"),
|
|
|
|
pika_choice_new_with_values ("derivative", 0, _("Derivative"), NULL,
|
|
|
|
"gradient", 1, _("Gradient"), NULL,
|
|
|
|
NULL),
|
|
|
|
"gradient",
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_CHOICE (procedure, "effect-convolve",
|
|
|
|
_("Con_volve"),
|
|
|
|
_("Convolve"),
|
|
|
|
pika_choice_new_with_values ("with-white-noise", 0, _("With white noise"), NULL,
|
|
|
|
"with-source-image", 1, _("With source image"), NULL,
|
|
|
|
NULL),
|
|
|
|
"with-source-image",
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_DRAWABLE (procedure, "effect-image",
|
|
|
|
_("Effect i_mage"),
|
|
|
|
_("Effect image"),
|
|
|
|
TRUE,
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_DOUBLE (procedure, "filter-length",
|
|
|
|
_("Fil_ter length"),
|
|
|
|
_("Filter length"),
|
|
|
|
0.1, 64.0, 5.0,
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_DOUBLE (procedure, "noise-magnitude",
|
|
|
|
_("_Noise Magnitude"),
|
|
|
|
_("Noise Magnitude"),
|
|
|
|
1.0, 5.0, 2.0,
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_DOUBLE (procedure, "integration-steps",
|
|
|
|
_("Inte_gration steps"),
|
|
|
|
_("Integration steps"),
|
|
|
|
1.0, 40.0, 25.0,
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_DOUBLE (procedure, "min-value",
|
|
|
|
_("Minimum v_alue"),
|
|
|
|
_("Minimum value"),
|
|
|
|
-100.0, 0, -25.0,
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
|
|
|
|
PIKA_PROC_ARG_DOUBLE (procedure, "max-value",
|
|
|
|
_("Ma_ximum value"),
|
|
|
|
_("Maximum value"),
|
|
|
|
0.0, 100, 25.0,
|
|
|
|
G_PARAM_READWRITE);
|
|
|
|
}
|
|
|
|
|
|
|
|
return procedure;
|
|
|
|
}
|
|
|
|
|
|
|
|
static PikaValueArray *
|
|
|
|
lic_run (PikaProcedure *procedure,
|
|
|
|
PikaRunMode run_mode,
|
|
|
|
PikaImage *image,
|
|
|
|
gint n_drawables,
|
|
|
|
PikaDrawable **drawables,
|
|
|
|
PikaProcedureConfig *config,
|
|
|
|
gpointer run_data)
|
|
|
|
{
|
|
|
|
PikaDrawable *drawable;
|
|
|
|
|
|
|
|
gegl_init (NULL, NULL);
|
|
|
|
|
|
|
|
if (n_drawables != 1)
|
|
|
|
{
|
|
|
|
GError *error = NULL;
|
|
|
|
|
|
|
|
g_set_error (&error, PIKA_PLUG_IN_ERROR, 0,
|
|
|
|
_("Procedure '%s' only works with one drawable."),
|
|
|
|
PLUG_IN_PROC);
|
|
|
|
|
|
|
|
return pika_procedure_new_return_values (procedure,
|
|
|
|
PIKA_PDB_CALLING_ERROR,
|
|
|
|
error);
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
drawable = drawables[0];
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Make sure that the drawable is RGBA or RGB color */
|
|
|
|
/* ================================================ */
|
|
|
|
|
|
|
|
if (pika_drawable_is_rgb (drawable))
|
|
|
|
{
|
|
|
|
switch (run_mode)
|
|
|
|
{
|
|
|
|
case PIKA_RUN_INTERACTIVE:
|
|
|
|
if (! create_main_dialog (procedure, config))
|
|
|
|
{
|
|
|
|
return pika_procedure_new_return_values (procedure,
|
|
|
|
PIKA_PDB_CANCEL,
|
|
|
|
NULL);
|
|
|
|
}
|
|
|
|
/* fallthrough */
|
|
|
|
|
|
|
|
case PIKA_RUN_WITH_LAST_VALS:
|
|
|
|
compute_image (config, drawable);
|
|
|
|
break;
|
|
|
|
|
|
|
|
default:
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
return pika_procedure_new_return_values (procedure,
|
|
|
|
PIKA_PDB_EXECUTION_ERROR,
|
|
|
|
NULL);
|
|
|
|
}
|
|
|
|
|
|
|
|
return pika_procedure_new_return_values (procedure, PIKA_PDB_SUCCESS, NULL);
|
|
|
|
}
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
/*****************************/
|
|
|
|
/* Global variables and such */
|
|
|
|
/*****************************/
|
|
|
|
|
|
|
|
static gdouble G[numx][numy][2];
|
|
|
|
|
|
|
|
static gdouble l = 10.0;
|
|
|
|
static gdouble dx = 2.0;
|
|
|
|
static gdouble dy = 2.0;
|
|
|
|
static gdouble minv = -2.5;
|
|
|
|
static gdouble maxv = 2.5;
|
|
|
|
static gdouble isteps = 20.0;
|
|
|
|
|
|
|
|
static gboolean source_drw_has_alpha = FALSE;
|
|
|
|
|
|
|
|
static gint effect_width, effect_height;
|
|
|
|
static gint border_x, border_y, border_w, border_h;
|
|
|
|
|
|
|
|
/************************/
|
|
|
|
/* Convenience routines */
|
|
|
|
/************************/
|
|
|
|
|
|
|
|
static void
|
|
|
|
peek (GeglBuffer *buffer,
|
|
|
|
gint x,
|
|
|
|
gint y,
|
|
|
|
PikaRGB *color)
|
|
|
|
{
|
|
|
|
gegl_buffer_sample (buffer, x, y, NULL,
|
|
|
|
color, babl_format ("R'G'B'A double"),
|
|
|
|
GEGL_SAMPLER_NEAREST, GEGL_ABYSS_NONE);
|
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
poke (GeglBuffer *buffer,
|
|
|
|
gint x,
|
|
|
|
gint y,
|
|
|
|
PikaRGB *color)
|
|
|
|
{
|
|
|
|
gegl_buffer_set (buffer, GEGL_RECTANGLE (x, y, 1, 1), 0,
|
|
|
|
babl_format ("R'G'B'A double"), color,
|
|
|
|
GEGL_AUTO_ROWSTRIDE);
|
|
|
|
}
|
|
|
|
|
|
|
|
static gint
|
|
|
|
peekmap (const guchar *image,
|
|
|
|
gint x,
|
|
|
|
gint y)
|
|
|
|
{
|
|
|
|
while (x < 0)
|
|
|
|
x += effect_width;
|
|
|
|
x %= effect_width;
|
|
|
|
|
|
|
|
while (y < 0)
|
|
|
|
y += effect_height;
|
|
|
|
y %= effect_height;
|
|
|
|
|
|
|
|
return (gint) image[x + effect_width * y];
|
|
|
|
}
|
|
|
|
|
|
|
|
/*************/
|
|
|
|
/* Main part */
|
|
|
|
/*************/
|
|
|
|
|
|
|
|
/***************************************************/
|
|
|
|
/* Compute the derivative in the x and y direction */
|
|
|
|
/* We use these convolution kernels: */
|
|
|
|
/* |1 0 -1| | 1 2 1| */
|
|
|
|
/* DX: |2 0 -2| DY: | 0 0 0| */
|
|
|
|
/* |1 0 -1| | -1 -2 -1| */
|
|
|
|
/* (It's a variation of the Sobel kernels, really) */
|
|
|
|
/***************************************************/
|
|
|
|
|
|
|
|
static gint
|
|
|
|
gradx (const guchar *image,
|
|
|
|
gint x,
|
|
|
|
gint y)
|
|
|
|
{
|
|
|
|
gint val = 0;
|
|
|
|
|
|
|
|
val = val + peekmap (image, x-1, y-1);
|
|
|
|
val = val - peekmap (image, x+1, y-1);
|
|
|
|
|
|
|
|
val = val + 2 * peekmap (image, x-1, y);
|
|
|
|
val = val - 2 * peekmap (image, x+1, y);
|
|
|
|
|
|
|
|
val = val + peekmap (image, x-1, y+1);
|
|
|
|
val = val - peekmap (image, x+1, y+1);
|
|
|
|
|
|
|
|
return val;
|
|
|
|
}
|
|
|
|
|
|
|
|
static gint
|
|
|
|
grady (const guchar *image,
|
|
|
|
gint x,
|
|
|
|
gint y)
|
|
|
|
{
|
|
|
|
gint val = 0;
|
|
|
|
|
|
|
|
val = val + peekmap (image, x-1, y-1);
|
|
|
|
val = val + 2 * peekmap (image, x, y-1);
|
|
|
|
val = val + peekmap (image, x+1, y-1);
|
|
|
|
|
|
|
|
val = val - peekmap (image, x-1, y+1);
|
|
|
|
val = val - 2 * peekmap (image, x, y+1);
|
|
|
|
val = val - peekmap (image, x+1, y+1);
|
|
|
|
|
|
|
|
return val;
|
|
|
|
}
|
|
|
|
|
|
|
|
/************************************/
|
|
|
|
/* A nice 2nd order cubic spline :) */
|
|
|
|
/************************************/
|
|
|
|
|
|
|
|
static gdouble
|
|
|
|
cubic (gdouble t)
|
|
|
|
{
|
|
|
|
gdouble at = fabs (t);
|
|
|
|
|
|
|
|
return (at < 1.0) ? at * at * (2.0 * at - 3.0) + 1.0 : 0.0;
|
|
|
|
}
|
|
|
|
|
|
|
|
static gdouble
|
|
|
|
omega (gdouble u,
|
|
|
|
gdouble v,
|
|
|
|
gint i,
|
|
|
|
gint j)
|
|
|
|
{
|
|
|
|
while (i < 0)
|
|
|
|
i += numx;
|
|
|
|
|
|
|
|
while (j < 0)
|
|
|
|
j += numy;
|
|
|
|
|
|
|
|
i %= numx;
|
|
|
|
j %= numy;
|
|
|
|
|
|
|
|
return cubic (u) * cubic (v) * (G[i][j][0]*u + G[i][j][1]*v);
|
|
|
|
}
|
|
|
|
|
|
|
|
/*************************************************************/
|
|
|
|
/* The noise function (2D variant of Perlins noise function) */
|
|
|
|
/*************************************************************/
|
|
|
|
|
|
|
|
static gdouble
|
|
|
|
noise (gdouble x,
|
|
|
|
gdouble y)
|
|
|
|
{
|
|
|
|
gint i, sti = (gint) floor (x / dx);
|
|
|
|
gint j, stj = (gint) floor (y / dy);
|
|
|
|
|
|
|
|
gdouble sum = 0.0;
|
|
|
|
|
|
|
|
/* Calculate the gdouble sum */
|
|
|
|
/* ======================== */
|
|
|
|
|
|
|
|
for (i = sti; i <= sti + 1; i++)
|
|
|
|
for (j = stj; j <= stj + 1; j++)
|
|
|
|
sum += omega ((x - (gdouble) i * dx) / dx,
|
|
|
|
(y - (gdouble) j * dy) / dy,
|
|
|
|
i, j);
|
|
|
|
|
|
|
|
return sum;
|
|
|
|
}
|
|
|
|
|
|
|
|
/*************************************************/
|
|
|
|
/* Generates pseudo-random vectors with length 1 */
|
|
|
|
/*************************************************/
|
|
|
|
|
|
|
|
static void
|
|
|
|
generatevectors (void)
|
|
|
|
{
|
|
|
|
gdouble alpha;
|
|
|
|
gint i, j;
|
|
|
|
GRand *gr;
|
|
|
|
|
|
|
|
gr = g_rand_new();
|
|
|
|
|
|
|
|
for (i = 0; i < numx; i++)
|
|
|
|
{
|
|
|
|
for (j = 0; j < numy; j++)
|
|
|
|
{
|
|
|
|
alpha = g_rand_double_range (gr, 0, 2) * G_PI;
|
|
|
|
G[i][j][0] = cos (alpha);
|
|
|
|
G[i][j][1] = sin (alpha);
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
g_rand_free (gr);
|
|
|
|
}
|
|
|
|
|
|
|
|
/* A simple triangle filter */
|
|
|
|
/* ======================== */
|
|
|
|
|
|
|
|
static gdouble
|
|
|
|
filter (gdouble u)
|
|
|
|
{
|
|
|
|
gdouble f = 1.0 - fabs (u) / l;
|
|
|
|
|
|
|
|
return (f < 0.0) ? 0.0 : f;
|
|
|
|
}
|
|
|
|
|
|
|
|
/******************************************************/
|
|
|
|
/* Compute the Line Integral Convolution (LIC) at x,y */
|
|
|
|
/******************************************************/
|
|
|
|
|
|
|
|
static gdouble
|
|
|
|
lic_noise (gint x,
|
|
|
|
gint y,
|
|
|
|
gdouble vx,
|
|
|
|
gdouble vy)
|
|
|
|
{
|
|
|
|
gdouble i = 0.0;
|
|
|
|
gdouble f1 = 0.0, f2 = 0.0;
|
|
|
|
gdouble u, step = 2.0 * l / isteps;
|
|
|
|
gdouble xx = (gdouble) x, yy = (gdouble) y;
|
|
|
|
gdouble c, s;
|
|
|
|
|
|
|
|
/* Get vector at x,y */
|
|
|
|
/* ================= */
|
|
|
|
|
|
|
|
c = vx;
|
|
|
|
s = vy;
|
|
|
|
|
|
|
|
/* Calculate integral numerically */
|
|
|
|
/* ============================== */
|
|
|
|
|
|
|
|
f1 = filter (-l) * noise (xx + l * c , yy + l * s);
|
|
|
|
|
|
|
|
for (u = -l + step; u <= l; u += step)
|
|
|
|
{
|
|
|
|
f2 = filter (u) * noise ( xx - u * c , yy - u * s);
|
|
|
|
i += (f1 + f2) * 0.5 * step;
|
|
|
|
f1 = f2;
|
|
|
|
}
|
|
|
|
|
|
|
|
i = (i - minv) / (maxv - minv);
|
|
|
|
|
|
|
|
i = CLAMP (i, 0.0, 1.0);
|
|
|
|
|
|
|
|
i = (i / 2.0) + 0.5;
|
|
|
|
|
|
|
|
return i;
|
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
getpixel (GeglBuffer *buffer,
|
|
|
|
PikaRGB *p,
|
|
|
|
gdouble u,
|
|
|
|
gdouble v)
|
|
|
|
{
|
|
|
|
register gint x1, y1, x2, y2;
|
|
|
|
gint width, height;
|
|
|
|
static PikaRGB pp[4];
|
|
|
|
|
|
|
|
width = border_w;
|
|
|
|
height = border_h;
|
|
|
|
|
|
|
|
x1 = (gint)u;
|
|
|
|
y1 = (gint)v;
|
|
|
|
|
|
|
|
if (x1 < 0)
|
|
|
|
x1 = width - (-x1 % width);
|
|
|
|
else
|
|
|
|
x1 = x1 % width;
|
|
|
|
|
|
|
|
if (y1 < 0)
|
|
|
|
y1 = height - (-y1 % height);
|
|
|
|
else
|
|
|
|
y1 = y1 % height;
|
|
|
|
|
|
|
|
x2 = (x1 + 1) % width;
|
|
|
|
y2 = (y1 + 1) % height;
|
|
|
|
|
|
|
|
peek (buffer, x1, y1, &pp[0]);
|
|
|
|
peek (buffer, x2, y1, &pp[1]);
|
|
|
|
peek (buffer, x1, y2, &pp[2]);
|
|
|
|
peek (buffer, x2, y2, &pp[3]);
|
|
|
|
|
|
|
|
if (source_drw_has_alpha)
|
|
|
|
*p = pika_bilinear_rgba (u, v, pp);
|
|
|
|
else
|
|
|
|
*p = pika_bilinear_rgb (u, v, pp);
|
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
lic_image (GeglBuffer *buffer,
|
|
|
|
gint x,
|
|
|
|
gint y,
|
|
|
|
gdouble vx,
|
|
|
|
gdouble vy,
|
|
|
|
PikaRGB *color)
|
|
|
|
{
|
|
|
|
gdouble u, step = 2.0 * l / isteps;
|
|
|
|
gdouble xx = (gdouble) x, yy = (gdouble) y;
|
|
|
|
gdouble c, s;
|
|
|
|
PikaRGB col = { 0, 0, 0, 0 };
|
|
|
|
PikaRGB col1, col2, col3;
|
|
|
|
|
|
|
|
/* Get vector at x,y */
|
|
|
|
/* ================= */
|
|
|
|
|
|
|
|
c = vx;
|
|
|
|
s = vy;
|
|
|
|
|
|
|
|
/* Calculate integral numerically */
|
|
|
|
/* ============================== */
|
|
|
|
|
|
|
|
getpixel (buffer, &col1, xx + l * c, yy + l * s);
|
|
|
|
|
|
|
|
if (source_drw_has_alpha)
|
|
|
|
pika_rgba_multiply (&col1, filter (-l));
|
|
|
|
else
|
|
|
|
pika_rgb_multiply (&col1, filter (-l));
|
|
|
|
|
|
|
|
for (u = -l + step; u <= l; u += step)
|
|
|
|
{
|
|
|
|
getpixel (buffer, &col2, xx - u * c, yy - u * s);
|
|
|
|
|
|
|
|
if (source_drw_has_alpha)
|
|
|
|
{
|
|
|
|
pika_rgba_multiply (&col2, filter (u));
|
|
|
|
|
|
|
|
col3 = col1;
|
|
|
|
pika_rgba_add (&col3, &col2);
|
|
|
|
pika_rgba_multiply (&col3, 0.5 * step);
|
|
|
|
pika_rgba_add (&col, &col3);
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
pika_rgb_multiply (&col2, filter (u));
|
|
|
|
|
|
|
|
col3 = col1;
|
|
|
|
pika_rgb_add (&col3, &col2);
|
|
|
|
pika_rgb_multiply (&col3, 0.5 * step);
|
|
|
|
pika_rgb_add (&col, &col3);
|
|
|
|
}
|
|
|
|
col1 = col2;
|
|
|
|
}
|
|
|
|
if (source_drw_has_alpha)
|
|
|
|
pika_rgba_multiply (&col, 1.0 / l);
|
|
|
|
else
|
|
|
|
pika_rgb_multiply (&col, 1.0 / l);
|
|
|
|
pika_rgb_clamp (&col);
|
|
|
|
|
|
|
|
*color = col;
|
|
|
|
}
|
|
|
|
|
|
|
|
static guchar *
|
|
|
|
rgb_to_hsl (PikaDrawable *drawable,
|
|
|
|
LICEffectChannel effect_channel)
|
|
|
|
{
|
|
|
|
GeglBuffer *buffer;
|
|
|
|
guchar *themap, data[4];
|
|
|
|
gint x, y;
|
|
|
|
PikaRGB color;
|
|
|
|
PikaHSL color_hsl;
|
|
|
|
gdouble val = 0.0;
|
|
|
|
glong maxc, index = 0;
|
|
|
|
GRand *gr;
|
|
|
|
|
|
|
|
gr = g_rand_new ();
|
|
|
|
|
|
|
|
maxc = pika_drawable_get_width (drawable) * pika_drawable_get_height (drawable);
|
|
|
|
|
|
|
|
buffer = pika_drawable_get_buffer (drawable);
|
|
|
|
|
|
|
|
themap = g_new (guchar, maxc);
|
|
|
|
|
|
|
|
for (y = 0; y < border_h; y++)
|
|
|
|
{
|
|
|
|
for (x = 0; x < border_w; x++)
|
|
|
|
{
|
|
|
|
data[3] = 255;
|
|
|
|
|
|
|
|
gegl_buffer_sample (buffer, x, y, NULL,
|
|
|
|
data, babl_format ("R'G'B'A u8"),
|
|
|
|
GEGL_SAMPLER_NEAREST, GEGL_ABYSS_NONE);
|
|
|
|
|
|
|
|
pika_rgba_set_uchar (&color, data[0], data[1], data[2], data[3]);
|
|
|
|
pika_rgb_to_hsl (&color, &color_hsl);
|
|
|
|
|
|
|
|
switch (effect_channel)
|
|
|
|
{
|
|
|
|
case LIC_HUE:
|
|
|
|
val = color_hsl.h * 255;
|
|
|
|
break;
|
|
|
|
case LIC_SATURATION:
|
|
|
|
val = color_hsl.s * 255;
|
|
|
|
break;
|
|
|
|
case LIC_BRIGHTNESS:
|
|
|
|
val = color_hsl.l * 255;
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* add some random to avoid unstructured areas. */
|
|
|
|
val += g_rand_double_range (gr, -1.0, 1.0);
|
|
|
|
|
|
|
|
themap[index++] = (guchar) CLAMP0255 (RINT (val));
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
g_object_unref (buffer);
|
|
|
|
|
|
|
|
g_rand_free (gr);
|
|
|
|
|
|
|
|
return themap;
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
static void
|
2023-10-30 23:55:30 +01:00
|
|
|
compute_lic (PikaProcedureConfig *config,
|
|
|
|
PikaDrawable *drawable,
|
|
|
|
const guchar *scalarfield,
|
|
|
|
gboolean rotate)
|
2023-09-26 00:35:21 +02:00
|
|
|
{
|
|
|
|
GeglBuffer *src_buffer;
|
|
|
|
GeglBuffer *dest_buffer;
|
|
|
|
gint xcount, ycount;
|
|
|
|
PikaRGB color;
|
|
|
|
gdouble vx, vy, tmp;
|
2023-10-30 23:55:30 +01:00
|
|
|
gint effect_convolve;
|
|
|
|
|
|
|
|
effect_convolve = pika_procedure_config_get_choice_id (config,
|
|
|
|
"effect-convolve");
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
src_buffer = pika_drawable_get_buffer (drawable);
|
|
|
|
dest_buffer = pika_drawable_get_shadow_buffer (drawable);
|
|
|
|
|
|
|
|
for (ycount = 0; ycount < border_h; ycount++)
|
|
|
|
{
|
|
|
|
for (xcount = 0; xcount < border_w; xcount++)
|
|
|
|
{
|
|
|
|
/* Get derivative at (x,y) and normalize it */
|
|
|
|
/* ============================================================== */
|
|
|
|
|
|
|
|
vx = gradx (scalarfield, xcount, ycount);
|
|
|
|
vy = grady (scalarfield, xcount, ycount);
|
|
|
|
|
|
|
|
/* Rotate if needed */
|
|
|
|
if (rotate)
|
|
|
|
{
|
|
|
|
tmp = vy;
|
|
|
|
vy = -vx;
|
|
|
|
vx = tmp;
|
|
|
|
}
|
|
|
|
|
|
|
|
tmp = sqrt (vx * vx + vy * vy);
|
|
|
|
if (tmp >= 0.000001)
|
|
|
|
{
|
|
|
|
tmp = 1.0 / tmp;
|
|
|
|
vx *= tmp;
|
|
|
|
vy *= tmp;
|
|
|
|
}
|
|
|
|
|
|
|
|
/* Convolve with the LIC at (x,y) */
|
|
|
|
/* ============================== */
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
if (effect_convolve == 0)
|
2023-09-26 00:35:21 +02:00
|
|
|
{
|
|
|
|
peek (src_buffer, xcount, ycount, &color);
|
|
|
|
|
|
|
|
tmp = lic_noise (xcount, ycount, vx, vy);
|
|
|
|
|
|
|
|
if (source_drw_has_alpha)
|
|
|
|
pika_rgba_multiply (&color, tmp);
|
|
|
|
else
|
|
|
|
pika_rgb_multiply (&color, tmp);
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
lic_image (src_buffer, xcount, ycount, vx, vy, &color);
|
|
|
|
}
|
|
|
|
|
|
|
|
poke (dest_buffer, xcount, ycount, &color);
|
|
|
|
}
|
|
|
|
|
|
|
|
pika_progress_update ((gfloat) ycount / (gfloat) border_h);
|
|
|
|
}
|
|
|
|
|
|
|
|
g_object_unref (src_buffer);
|
|
|
|
g_object_unref (dest_buffer);
|
|
|
|
|
|
|
|
pika_progress_update (1.0);
|
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
2023-10-30 23:55:30 +01:00
|
|
|
compute_image (PikaProcedureConfig *config,
|
|
|
|
PikaDrawable *drawable)
|
2023-09-26 00:35:21 +02:00
|
|
|
{
|
|
|
|
PikaDrawable *effect_image;
|
|
|
|
guchar *scalarfield = NULL;
|
2023-10-30 23:55:30 +01:00
|
|
|
gdouble filtlen;
|
|
|
|
gdouble noisemag;
|
|
|
|
gdouble intsteps;
|
|
|
|
gdouble minv;
|
|
|
|
gdouble maxv;
|
|
|
|
gint effect_channel;
|
|
|
|
gint effect_operator;
|
|
|
|
gint effect_convolve;
|
|
|
|
|
|
|
|
g_object_get (config,
|
|
|
|
"filter-length", &filtlen,
|
|
|
|
"noise-magnitude", &noisemag,
|
|
|
|
"effect-image", &effect_image,
|
|
|
|
"integration-steps", &intsteps,
|
|
|
|
"min-value", &minv,
|
|
|
|
"max-value", &maxv,
|
|
|
|
NULL);
|
|
|
|
effect_channel = pika_procedure_config_get_choice_id (config,
|
|
|
|
"effect-channel");
|
|
|
|
effect_operator = pika_procedure_config_get_choice_id (config,
|
|
|
|
"effect-operator");
|
|
|
|
effect_convolve = pika_procedure_config_get_choice_id (config,
|
|
|
|
"effect-convolve");
|
|
|
|
|
|
|
|
if (! effect_image)
|
|
|
|
return;
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
/* Get some useful info on the input drawable */
|
|
|
|
/* ========================================== */
|
|
|
|
if (! pika_drawable_mask_intersect (drawable,
|
|
|
|
&border_x, &border_y,
|
|
|
|
&border_w, &border_h))
|
|
|
|
return;
|
|
|
|
|
|
|
|
pika_progress_init (_("Van Gogh (LIC)"));
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
if (effect_convolve == 0)
|
2023-09-26 00:35:21 +02:00
|
|
|
generatevectors ();
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
if (filtlen < 0.1)
|
|
|
|
filtlen = 0.1;
|
2023-09-26 00:35:21 +02:00
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
l = filtlen;
|
|
|
|
dx = dy = noisemag;
|
|
|
|
minv = minv / 10.0;
|
|
|
|
maxv = maxv / 10.0;
|
|
|
|
isteps = intsteps;
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
source_drw_has_alpha = pika_drawable_has_alpha (drawable);
|
|
|
|
|
|
|
|
effect_width = pika_drawable_get_width (effect_image);
|
|
|
|
effect_height = pika_drawable_get_height (effect_image);
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
switch (effect_channel)
|
2023-09-26 00:35:21 +02:00
|
|
|
{
|
2023-10-30 23:55:30 +01:00
|
|
|
case LIC_HUE:
|
2023-09-26 00:35:21 +02:00
|
|
|
scalarfield = rgb_to_hsl (effect_image, LIC_HUE);
|
|
|
|
break;
|
2023-10-30 23:55:30 +01:00
|
|
|
case LIC_SATURATION:
|
2023-09-26 00:35:21 +02:00
|
|
|
scalarfield = rgb_to_hsl (effect_image, LIC_SATURATION);
|
|
|
|
break;
|
2023-10-30 23:55:30 +01:00
|
|
|
case LIC_BRIGHTNESS:
|
2023-09-26 00:35:21 +02:00
|
|
|
scalarfield = rgb_to_hsl (effect_image, LIC_BRIGHTNESS);
|
|
|
|
break;
|
|
|
|
}
|
2023-10-30 23:55:30 +01:00
|
|
|
g_clear_object (&effect_image);
|
2023-09-26 00:35:21 +02:00
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
compute_lic (config, drawable, scalarfield, effect_operator);
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
g_free (scalarfield);
|
|
|
|
|
|
|
|
/* Update image */
|
|
|
|
/* ============ */
|
|
|
|
|
|
|
|
pika_drawable_merge_shadow (drawable, TRUE);
|
|
|
|
pika_drawable_update (drawable, border_x, border_y, border_w, border_h);
|
|
|
|
|
|
|
|
pika_displays_flush ();
|
|
|
|
}
|
|
|
|
|
|
|
|
/**************************/
|
|
|
|
/* Below is only UI stuff */
|
|
|
|
/**************************/
|
|
|
|
|
|
|
|
static gboolean
|
2023-10-30 23:55:30 +01:00
|
|
|
create_main_dialog (PikaProcedure *procedure,
|
|
|
|
PikaProcedureConfig *config)
|
2023-09-26 00:35:21 +02:00
|
|
|
{
|
2023-10-30 23:55:30 +01:00
|
|
|
GtkWidget *dialog;
|
2023-09-26 00:35:21 +02:00
|
|
|
gboolean run;
|
|
|
|
|
|
|
|
pika_ui_init (PLUG_IN_BINARY);
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
dialog = pika_procedure_dialog_new (procedure,
|
|
|
|
PIKA_PROCEDURE_CONFIG (config),
|
|
|
|
_("Van Gogh (LIC)"));
|
|
|
|
|
|
|
|
pika_procedure_dialog_get_scale_entry (PIKA_PROCEDURE_DIALOG (dialog),
|
|
|
|
"filter-length", 1.0);
|
|
|
|
pika_procedure_dialog_get_scale_entry (PIKA_PROCEDURE_DIALOG (dialog),
|
|
|
|
"noise-magnitude", 1.0);
|
|
|
|
pika_procedure_dialog_get_scale_entry (PIKA_PROCEDURE_DIALOG (dialog),
|
|
|
|
"integration-steps", 1.0);
|
|
|
|
pika_procedure_dialog_get_scale_entry (PIKA_PROCEDURE_DIALOG (dialog),
|
|
|
|
"min-value", 1.0);
|
|
|
|
pika_procedure_dialog_get_scale_entry (PIKA_PROCEDURE_DIALOG (dialog),
|
|
|
|
"max-value", 1.0);
|
|
|
|
|
|
|
|
pika_procedure_dialog_fill (PIKA_PROCEDURE_DIALOG (dialog), NULL);
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
gtk_widget_show (dialog);
|
|
|
|
|
2023-10-30 23:55:30 +01:00
|
|
|
run = pika_procedure_dialog_run (PIKA_PROCEDURE_DIALOG (dialog));
|
2023-09-26 00:35:21 +02:00
|
|
|
|
|
|
|
gtk_widget_destroy (dialog);
|
|
|
|
|
|
|
|
return run;
|
|
|
|
}
|