app: add gimp-gegl-loops.[ch] which for now contains a convolve function
This commit is contained in:
Michael Natterer
@ -22,6 +22,8 @@ libappgegl_a_sources = \
|
||||
gimp-gegl.h \
|
||||
gimp-gegl-config-proxy.c \
|
||||
gimp-gegl-config-proxy.h \
|
||||
gimp-gegl-loops.c \
|
||||
gimp-gegl-loops.h \
|
||||
gimp-gegl-nodes.c \
|
||||
gimp-gegl-nodes.h \
|
||||
gimp-gegl-utils.c \
|
||||
|
||||
191
app/gegl/gimp-gegl-loops.c
Normal file
191
app/gegl/gimp-gegl-loops.c
Normal file
@ -0,0 +1,191 @@
|
||||
/* GIMP - The GNU Image Manipulation Program
|
||||
* Copyright (C) 1995 Spencer Kimball and Peter Mattis
|
||||
*
|
||||
* gimp-gegl-loops.c
|
||||
* Copyright (C) 2012 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 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/>.
|
||||
*/
|
||||
|
||||
#include "config.h"
|
||||
|
||||
#include <gegl.h>
|
||||
|
||||
#include "libgimpmath/gimpmath.h"
|
||||
|
||||
#include "gimp-gegl-types.h"
|
||||
|
||||
#include "gimp-gegl-loops.h"
|
||||
|
||||
|
||||
void
|
||||
gimp_gegl_convolve (GeglBuffer *src_buffer,
|
||||
const GeglRectangle *src_rect,
|
||||
GeglBuffer *dest_buffer,
|
||||
const GeglRectangle *dest_rect,
|
||||
const gfloat *matrix,
|
||||
gint size,
|
||||
gdouble divisor,
|
||||
GimpConvolutionType mode,
|
||||
gboolean alpha_weighting)
|
||||
{
|
||||
GeglBufferIterator *iter;
|
||||
GeglRectangle *src_roi;
|
||||
GeglRectangle *dest_roi;
|
||||
const Babl *src_format;
|
||||
const Babl *dest_format;
|
||||
gint src_bpp;
|
||||
gint dest_bpp;
|
||||
|
||||
src_format = gegl_buffer_get_format (src_buffer);
|
||||
dest_format = gegl_buffer_get_format (dest_buffer);
|
||||
|
||||
src_bpp = babl_format_get_bytes_per_pixel (src_format);
|
||||
dest_bpp = babl_format_get_bytes_per_pixel (dest_format);
|
||||
|
||||
iter = gegl_buffer_iterator_new (src_buffer, src_rect, 0, NULL,
|
||||
GEGL_BUFFER_READ, GEGL_ABYSS_NONE);
|
||||
src_roi = &iter->roi[0];
|
||||
|
||||
gegl_buffer_iterator_add (iter, dest_buffer, dest_rect, 0, NULL,
|
||||
GEGL_BUFFER_WRITE, GEGL_ABYSS_NONE);
|
||||
dest_roi = &iter->roi[1];
|
||||
|
||||
while (gegl_buffer_iterator_next (iter))
|
||||
{
|
||||
/* Convolve the src image using the convolution matrix, writing
|
||||
* to dest Convolve is not tile-enabled--use accordingly
|
||||
*/
|
||||
const guchar *src = iter->data[0];
|
||||
guchar *dest = iter->data[1];
|
||||
const gint bytes = src_bpp;
|
||||
const gint a_byte = bytes - 1;
|
||||
const gint rowstride = src_bpp * src_roi->width;
|
||||
const gint margin = size / 2;
|
||||
const gint x1 = src_roi->x;
|
||||
const gint y1 = src_roi->y;
|
||||
const gint x2 = src_roi->x + src_roi->width - 1;
|
||||
const gint y2 = src_roi->y + src_roi->height - 1;
|
||||
gint x, y;
|
||||
gint offset;
|
||||
|
||||
/* If the mode is NEGATIVE_CONVOL, the offset should be 128 */
|
||||
if (mode == GIMP_NEGATIVE_CONVOL)
|
||||
{
|
||||
offset = 128;
|
||||
mode = GIMP_NORMAL_CONVOL;
|
||||
}
|
||||
else
|
||||
{
|
||||
offset = 0;
|
||||
}
|
||||
|
||||
for (y = 0; y < dest_roi->height; y++)
|
||||
{
|
||||
guchar *d = dest;
|
||||
|
||||
if (alpha_weighting)
|
||||
{
|
||||
for (x = 0; x < dest_roi->width; x++)
|
||||
{
|
||||
const gfloat *m = matrix;
|
||||
gdouble total[4] = { 0.0, 0.0, 0.0, 0.0 };
|
||||
gdouble weighted_divisor = 0.0;
|
||||
gint i, j, b;
|
||||
|
||||
for (j = y - margin; j <= y + margin; j++)
|
||||
{
|
||||
for (i = x - margin; i <= x + margin; i++, m++)
|
||||
{
|
||||
gint xx = CLAMP (i, x1, x2);
|
||||
gint yy = CLAMP (j, y1, y2);
|
||||
const guchar *s = src + yy * rowstride + xx * bytes;
|
||||
const guchar a = s[a_byte];
|
||||
|
||||
if (a)
|
||||
{
|
||||
gdouble mult_alpha = *m * a;
|
||||
|
||||
weighted_divisor += mult_alpha;
|
||||
|
||||
for (b = 0; b < a_byte; b++)
|
||||
total[b] += mult_alpha * s[b];
|
||||
|
||||
total[a_byte] += mult_alpha;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
if (weighted_divisor == 0.0)
|
||||
weighted_divisor = divisor;
|
||||
|
||||
for (b = 0; b < a_byte; b++)
|
||||
total[b] /= weighted_divisor;
|
||||
|
||||
total[a_byte] /= divisor;
|
||||
|
||||
for (b = 0; b < bytes; b++)
|
||||
{
|
||||
total[b] += offset;
|
||||
|
||||
if (mode != GIMP_NORMAL_CONVOL && total[b] < 0.0)
|
||||
total[b] = - total[b];
|
||||
|
||||
if (total[b] < 0.0)
|
||||
*d++ = 0;
|
||||
else
|
||||
*d++ = (total[b] > 255.0) ? 255 : (guchar) ROUND (total[b]);
|
||||
}
|
||||
}
|
||||
}
|
||||
else
|
||||
{
|
||||
for (x = 0; x < dest_roi->width; x++)
|
||||
{
|
||||
const gfloat *m = matrix;
|
||||
gdouble total[4] = { 0.0, 0.0, 0.0, 0.0 };
|
||||
gint i, j, b;
|
||||
|
||||
for (j = y - margin; j <= y + margin; j++)
|
||||
{
|
||||
for (i = x - margin; i <= x + margin; i++, m++)
|
||||
{
|
||||
gint xx = CLAMP (i, x1, x2);
|
||||
gint yy = CLAMP (j, y1, y2);
|
||||
const guchar *s = src + yy * rowstride + xx * bytes;
|
||||
|
||||
for (b = 0; b < bytes; b++)
|
||||
total[b] += *m * s[b];
|
||||
}
|
||||
}
|
||||
|
||||
for (b = 0; b < bytes; b++)
|
||||
{
|
||||
total[b] = total[b] / divisor + offset;
|
||||
|
||||
if (mode != GIMP_NORMAL_CONVOL && total[b] < 0.0)
|
||||
total[b] = - total[b];
|
||||
|
||||
if (total[b] < 0.0)
|
||||
*d++ = 0.0;
|
||||
else
|
||||
*d++ = (total[b] > 255.0) ? 255 : (guchar) ROUND (total[b]);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
dest += dest_roi->width * dest_bpp;
|
||||
}
|
||||
}
|
||||
}
|
||||
39
app/gegl/gimp-gegl-loops.h
Normal file
39
app/gegl/gimp-gegl-loops.h
Normal file
@ -0,0 +1,39 @@
|
||||
/* GIMP - The GNU Image Manipulation Program
|
||||
* Copyright (C) 1995 Spencer Kimball and Peter Mattis
|
||||
*
|
||||
* gimp-gegl-loops.h
|
||||
* Copyright (C) 2012 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 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/>.
|
||||
*/
|
||||
|
||||
#ifndef __GIMP_GEGL_LOOPS_H__
|
||||
#define __GIMP_GEGL_LOOPS_H__
|
||||
|
||||
|
||||
/* this is a pretty stupid port of concolve_region() that only works
|
||||
* on a linear source buffer
|
||||
*/
|
||||
void gimp_gegl_convolve (GeglBuffer *src_buffer,
|
||||
const GeglRectangle *src_rect,
|
||||
GeglBuffer *dest_buffer,
|
||||
const GeglRectangle *dest_rect,
|
||||
const gfloat *matrix,
|
||||
gint size,
|
||||
gdouble divisor,
|
||||
GimpConvolutionType mode,
|
||||
gboolean alpha_weighting);
|
||||
|
||||
|
||||
#endif /* __GIMP_GEGL_LOOPS_H__ */
|
||||
Reference in New Issue
Block a user