AlderconeStudio/PIKApp
AlderconeStudio
/
PIKApp
0
0
Fork 0
Code Issues 1 Pull Requests Packages Projects Releases 1 Wiki Activity
PIKApp/modules/display-filter-color-blind.c

500 lines
17 KiB
C
Raw Permalink Blame History

/* PIKA - Photo and Image Kooker Application
* a rebranding of The GNU Image Manipulation Program (created with heckimp)
* A derived work which may be trivial. However, any changes may be (C)2023 by Aldercone Studio
*
* Original copyright, applying to most contents (license remains unchanged):
* Copyright (C) 1995-1997 Spencer Kimball and Peter Mattis
*
* cdisplay_colorblind.c
* Copyright (C) 2002-2003 Michael Natterer <mitch@gimp.org>,
* Sven Neumann <sven@gimp.org>,
* Robert Dougherty <bob@vischeck.com> and
* Alex Wade <alex@vischeck.com>
*
* This code is an implementation of an algorithm described by Hans Brettel,
* Francoise Vienot and John Mollon in the Journal of the Optical Society of
* America V14(10), pg 2647. (See http://vischeck.com/ for more info.)
*
* 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 <string.h>
#include <gegl.h>
#include <gtk/gtk.h>
#include "libpikacolor/pikacolor.h"
#include "libpikaconfig/pikaconfig.h"
#include "libpikamath/pikamath.h"
#include "libpikamodule/pikamodule.h"
#include "libpikawidgets/pikawidgets.h"
#include "libpika/libpika-intl.h"
typedef enum
{
COLORBLIND_DEFICIENCY_PROTANOPIA,
COLORBLIND_DEFICIENCY_DEUTERANOPIA,
COLORBLIND_DEFICIENCY_TRITANOPIA
} ColorblindDeficiencyType;
#define CDISPLAY_TYPE_COLORBLIND_DEFICIENCY_TYPE (cdisplay_colorblind_deficiency_type_type)
static GType cdisplay_colorblind_deficiency_type_register_type (GTypeModule *module);
static const GEnumValue enum_values[] =
{
{ COLORBLIND_DEFICIENCY_PROTANOPIA,
"COLORBLIND_DEFICIENCY_PROTANOPIA", "protanopia" },
{ COLORBLIND_DEFICIENCY_DEUTERANOPIA,
"COLORBLIND_DEFICIENCY_DEUTERANOPIA", "deuteranopia" },
{ COLORBLIND_DEFICIENCY_TRITANOPIA,
"COLORBLIND_DEFICIENCY_TRITANOPIA", "tritanopia" },
{ 0, NULL, NULL }
};
static const PikaEnumDesc enum_descs[] =
{
{ COLORBLIND_DEFICIENCY_PROTANOPIA,
N_("Protanopia (insensitivity to red)"), NULL },
{ COLORBLIND_DEFICIENCY_DEUTERANOPIA,
N_("Deuteranopia (insensitivity to green)"), NULL },
{ COLORBLIND_DEFICIENCY_TRITANOPIA,
N_("Tritanopia (insensitivity to blue)"), NULL },
{ 0, NULL, NULL }
};
#define DEFAULT_DEFICIENCY_TYPE COLORBLIND_DEFICIENCY_DEUTERANOPIA
#define COLOR_CACHE_SIZE 1021
#define CDISPLAY_TYPE_COLORBLIND (cdisplay_colorblind_get_type ())
#define CDISPLAY_COLORBLIND(obj) (G_TYPE_CHECK_INSTANCE_CAST ((obj), CDISPLAY_TYPE_COLORBLIND, CdisplayColorblind))
#define CDISPLAY_COLORBLIND_CLASS(klass) (G_TYPE_CHECK_CLASS_CAST ((klass), CDISPLAY_TYPE_COLORBLIND, CdisplayColorblindClass))
#define CDISPLAY_IS_COLORBLIND(obj) (G_TYPE_CHECK_INSTANCE_TYPE ((obj), CDISPLAY_TYPE_COLORBLIND))
#define CDISPLAY_IS_COLORBLIND_CLASS(klass) (G_TYPE_CHECK_CLASS_TYPE ((klass), CDISPLAY_TYPE_COLORBLIND))
typedef struct _CdisplayColorblind CdisplayColorblind;
typedef struct _CdisplayColorblindClass CdisplayColorblindClass;
struct _CdisplayColorblind
{
PikaColorDisplay parent_instance;
ColorblindDeficiencyType type;
gfloat a1, b1, c1;
gfloat a2, b2, c2;
gfloat inflection;
};
struct _CdisplayColorblindClass
{
PikaColorDisplayClass parent_instance;
};
enum
{
PROP_0,
PROP_TYPE
};
GType cdisplay_colorblind_get_type (void);
static void cdisplay_colorblind_set_property (GObject *object,
guint property_id,
const GValue *value,
GParamSpec *pspec);
static void cdisplay_colorblind_get_property (GObject *object,
guint property_id,
GValue *value,
GParamSpec *pspec);
static void cdisplay_colorblind_convert_buffer (PikaColorDisplay *display,
GeglBuffer *buffer,
GeglRectangle *area);
static void cdisplay_colorblind_changed (PikaColorDisplay *display);
static void cdisplay_colorblind_set_type (CdisplayColorblind *colorblind,
ColorblindDeficiencyType value);
/* The RGB<->LMS transforms above are computed from the human cone
* photo-pigment absorption spectra and the monitor phosphor
* emission spectra. These parameters are fairly constant for most
* humans and most monitors (at least for modern CRTs). However,
* gamma will vary quite a bit, as it is a property of the monitor
* (eg. amplifier gain), the video card, and even the
* software. Further, users can adjust their gammas (either via
* adjusting the monitor amp gains or in software). That said, the
* following are the gamma estimates that we have used in the
* Vischeck code. Many colorblind users have viewed our simulations
* and told us that they "work" (simulated and original images are
* indistinguishable).
*/
#if 0
/* Gamma conversion is now handled by simply asking for a linear buffer */
static const gfloat gammaRGB = 2.1;
#endif
/* For most modern Cathode-Ray Tube monitors (CRTs), the following
* are good estimates of the RGB->LMS and LMS->RGB transform
* matrices. They are based on spectra measured on a typical CRT
* with a PhotoResearch PR650 spectral photometer and the Stockman
* human cone fundamentals. NOTE: these estimates will NOT work well
* for LCDs!
*/
static const gfloat rgb2lms[9] =
{
0.05059983,
0.08585369,
0.00952420,
0.01893033,
0.08925308,
0.01370054,
0.00292202,
0.00975732,
0.07145979
};
static const gfloat lms2rgb[9] =
{
30.830854,
-29.832659,
1.610474,
-6.481468,
17.715578,
-2.532642,
-0.375690,
-1.199062,
14.273846
};
static const PikaModuleInfo cdisplay_colorblind_info =
{
PIKA_MODULE_ABI_VERSION,
N_("Color deficit simulation filter (Brettel-Vienot-Mollon algorithm)"),
"Michael Natterer <mitch@gimp.org>, Bob Dougherty <bob@vischeck.com>, "
"Alex Wade <alex@vischeck.com>",
"v0.2",
"(c) 2002-2004, released under the GPL",
"January 22, 2003"
};
G_DEFINE_DYNAMIC_TYPE (CdisplayColorblind, cdisplay_colorblind,
PIKA_TYPE_COLOR_DISPLAY)
static GType cdisplay_colorblind_deficiency_type_type = 0;
G_MODULE_EXPORT const PikaModuleInfo *
pika_module_query (GTypeModule *module)
{
return &cdisplay_colorblind_info;
}
G_MODULE_EXPORT gboolean
pika_module_register (GTypeModule *module)
{
cdisplay_colorblind_register_type (module);
cdisplay_colorblind_deficiency_type_register_type (module);
return TRUE;
}
static GType
cdisplay_colorblind_deficiency_type_register_type (GTypeModule *module)
{
if (! cdisplay_colorblind_deficiency_type_type)
{
cdisplay_colorblind_deficiency_type_type =
g_type_module_register_enum (module, "CDisplayColorblindDeficiencyType",
enum_values);
pika_type_set_translation_domain (cdisplay_colorblind_deficiency_type_type,
GETTEXT_PACKAGE "-libpika");
pika_enum_set_value_descriptions (cdisplay_colorblind_deficiency_type_type,
enum_descs);
}
return cdisplay_colorblind_deficiency_type_type;
}
static void
cdisplay_colorblind_class_init (CdisplayColorblindClass *klass)
{
GObjectClass *object_class = G_OBJECT_CLASS (klass);
PikaColorDisplayClass *display_class = PIKA_COLOR_DISPLAY_CLASS (klass);
object_class->get_property = cdisplay_colorblind_get_property;
object_class->set_property = cdisplay_colorblind_set_property;
PIKA_CONFIG_PROP_ENUM (object_class, PROP_TYPE,
"type",
_("Type"),
_("Color vision deficiency type"),
CDISPLAY_TYPE_COLORBLIND_DEFICIENCY_TYPE,
DEFAULT_DEFICIENCY_TYPE,
0);
display_class->name = _("Color Deficient Vision");
display_class->help_id = "pika-colordisplay-colorblind";
display_class->icon_name = PIKA_ICON_DISPLAY_FILTER_COLORBLIND;
display_class->convert_buffer = cdisplay_colorblind_convert_buffer;
display_class->changed = cdisplay_colorblind_changed;
}
static void
cdisplay_colorblind_class_finalize (CdisplayColorblindClass *klass)
{
}
static void
cdisplay_colorblind_init (CdisplayColorblind *colorblind)
{
}
static void
cdisplay_colorblind_get_property (GObject *object,
guint property_id,
GValue *value,
GParamSpec *pspec)
{
CdisplayColorblind *colorblind = CDISPLAY_COLORBLIND (object);
switch (property_id)
{
case PROP_TYPE:
g_value_set_enum (value, colorblind->type);
break;
default:
G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
break;
}
}
static void
cdisplay_colorblind_set_property (GObject *object,
guint property_id,
const GValue *value,
GParamSpec *pspec)
{
CdisplayColorblind *colorblind = CDISPLAY_COLORBLIND (object);
switch (property_id)
{
case PROP_TYPE:
cdisplay_colorblind_set_type (colorblind, g_value_get_enum (value));
break;
default:
G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
break;
}
}
static void
cdisplay_colorblind_convert_buffer (PikaColorDisplay *display,
GeglBuffer *buffer,
GeglRectangle *area)
{
CdisplayColorblind *colorblind = CDISPLAY_COLORBLIND (display);
GeglBufferIterator *iter;
const gfloat a1 = colorblind->a1;
const gfloat b1 = colorblind->b1;
const gfloat c1 = colorblind->c1;
const gfloat a2 = colorblind->a2;
const gfloat b2 = colorblind->b2;
const gfloat c2 = colorblind->c2;
iter = gegl_buffer_iterator_new (buffer, area, 0,
babl_format ("RGBA float") /* linear! */,
GEGL_ACCESS_READWRITE, GEGL_ABYSS_NONE, 1);
while (gegl_buffer_iterator_next (iter))
{
gfloat *data = iter->items[0].data;
gint count = iter->length;
while (count--)
{
gfloat tmp;
gfloat red, green, blue;
gfloat redOld, greenOld;
red = data[0];
green = data[1];
blue = data[2];
/* Convert to LMS (dot product with transform matrix) */
redOld = red;
greenOld = green;
red = redOld * rgb2lms[0] + greenOld * rgb2lms[1] + blue * rgb2lms[2];
green = redOld * rgb2lms[3] + greenOld * rgb2lms[4] + blue * rgb2lms[5];
blue = redOld * rgb2lms[6] + greenOld * rgb2lms[7] + blue * rgb2lms[8];
switch (colorblind->type)
{
case COLORBLIND_DEFICIENCY_DEUTERANOPIA:
tmp = blue / red;
/* See which side of the inflection line we fall... */
if (tmp < colorblind->inflection)
green = -(a1 * red + c1 * blue) / b1;
else
green = -(a2 * red + c2 * blue) / b2;
break;
case COLORBLIND_DEFICIENCY_PROTANOPIA:
tmp = blue / green;
/* See which side of the inflection line we fall... */
if (tmp < colorblind->inflection)
red = -(b1 * green + c1 * blue) / a1;
else
red = -(b2 * green + c2 * blue) / a2;
break;
case COLORBLIND_DEFICIENCY_TRITANOPIA:
tmp = green / red;
/* See which side of the inflection line we fall... */
if (tmp < colorblind->inflection)
blue = -(a1 * red + b1 * green) / c1;
else
blue = -(a2 * red + b2 * green) / c2;
break;
default:
break;
}
/* Convert back to RGB (cross product with transform matrix) */
redOld = red;
greenOld = green;
red = redOld * lms2rgb[0] + greenOld * lms2rgb[1] + blue * lms2rgb[2];
green = redOld * lms2rgb[3] + greenOld * lms2rgb[4] + blue * lms2rgb[5];
blue = redOld * lms2rgb[6] + greenOld * lms2rgb[7] + blue * lms2rgb[8];
data[0] = red;
data[1] = green;
data[2] = blue;
data += 4;
}
}
}
static void
cdisplay_colorblind_changed (PikaColorDisplay *display)
{
CdisplayColorblind *colorblind = CDISPLAY_COLORBLIND (display);
gfloat anchor_e[3];
gfloat anchor[12];
/* This function performs initialisations that are dependent
* on the type of color deficiency.
*/
/* Performs protan, deutan or tritan color image simulation based on
* Brettel, Vienot and Mollon JOSA 14/10 1997
* L,M,S for lambda=475,485,575,660
*
* Load the LMS anchor-point values for lambda = 475 & 485 nm (for
* protans & deutans) and the LMS values for lambda = 575 & 660 nm
* (for tritans)
*/
anchor[0] = 0.08008; anchor[1] = 0.1579; anchor[2] = 0.5897;
anchor[3] = 0.1284; anchor[4] = 0.2237; anchor[5] = 0.3636;
anchor[6] = 0.9856; anchor[7] = 0.7325; anchor[8] = 0.001079;
anchor[9] = 0.0914; anchor[10] = 0.007009; anchor[11] = 0.0;
/* We also need LMS for RGB=(1,1,1)- the equal-energy point (one of
* our anchors) (we can just peel this out of the rgb2lms transform
* matrix)
*/
anchor_e[0] = rgb2lms[0] + rgb2lms[1] + rgb2lms[2];
anchor_e[1] = rgb2lms[3] + rgb2lms[4] + rgb2lms[5];
anchor_e[2] = rgb2lms[6] + rgb2lms[7] + rgb2lms[8];
switch (colorblind->type)
{
case COLORBLIND_DEFICIENCY_DEUTERANOPIA:
/* find a,b,c for lam=575nm and lam=475 */
colorblind->a1 = anchor_e[1] * anchor[8] - anchor_e[2] * anchor[7];
colorblind->b1 = anchor_e[2] * anchor[6] - anchor_e[0] * anchor[8];
colorblind->c1 = anchor_e[0] * anchor[7] - anchor_e[1] * anchor[6];
colorblind->a2 = anchor_e[1] * anchor[2] - anchor_e[2] * anchor[1];
colorblind->b2 = anchor_e[2] * anchor[0] - anchor_e[0] * anchor[2];
colorblind->c2 = anchor_e[0] * anchor[1] - anchor_e[1] * anchor[0];
colorblind->inflection = (anchor_e[2] / anchor_e[0]);
break;
case COLORBLIND_DEFICIENCY_PROTANOPIA:
/* find a,b,c for lam=575nm and lam=475 */
colorblind->a1 = anchor_e[1] * anchor[8] - anchor_e[2] * anchor[7];
colorblind->b1 = anchor_e[2] * anchor[6] - anchor_e[0] * anchor[8];
colorblind->c1 = anchor_e[0] * anchor[7] - anchor_e[1] * anchor[6];
colorblind->a2 = anchor_e[1] * anchor[2] - anchor_e[2] * anchor[1];
colorblind->b2 = anchor_e[2] * anchor[0] - anchor_e[0] * anchor[2];
colorblind->c2 = anchor_e[0] * anchor[1] - anchor_e[1] * anchor[0];
colorblind->inflection = (anchor_e[2] / anchor_e[1]);
break;
case COLORBLIND_DEFICIENCY_TRITANOPIA:
/* Set 1: regions where lambda_a=575, set 2: lambda_a=475 */
colorblind->a1 = anchor_e[1] * anchor[11] - anchor_e[2] * anchor[10];
colorblind->b1 = anchor_e[2] * anchor[9] - anchor_e[0] * anchor[11];
colorblind->c1 = anchor_e[0] * anchor[10] - anchor_e[1] * anchor[9];
colorblind->a2 = anchor_e[1] * anchor[5] - anchor_e[2] * anchor[4];
colorblind->b2 = anchor_e[2] * anchor[3] - anchor_e[0] * anchor[5];
colorblind->c2 = anchor_e[0] * anchor[4] - anchor_e[1] * anchor[3];
colorblind->inflection = (anchor_e[1] / anchor_e[0]);
break;
}
}
static void
cdisplay_colorblind_set_type (CdisplayColorblind *colorblind,
ColorblindDeficiencyType value)
{
if (value != colorblind->type)
{
GEnumClass *enum_class;
enum_class = g_type_class_peek (CDISPLAY_TYPE_COLORBLIND_DEFICIENCY_TYPE);
if (! g_enum_get_value (enum_class, value))
return;
colorblind->type = value;
g_object_notify (G_OBJECT (colorblind), "type");
pika_color_display_changed (PIKA_COLOR_DISPLAY (colorblind));
}
}
Reference in New Issue View Git Blame Copy Permalink