summaryrefslogtreecommitdiffstats
path: root/modules/display-filter-color-blind.c
diff options
context:
space:
mode:
authorDaniel Baumann <daniel.baumann@progress-linux.org>2024-04-07 18:30:19 +0000
committerDaniel Baumann <daniel.baumann@progress-linux.org>2024-04-07 18:30:19 +0000
commit5c1676dfe6d2f3c837a5e074117b45613fd29a72 (patch)
treecbffb45144febf451e54061db2b21395faf94bfe /modules/display-filter-color-blind.c
parentInitial commit. (diff)
downloadgimp-5c1676dfe6d2f3c837a5e074117b45613fd29a72.tar.xz
gimp-5c1676dfe6d2f3c837a5e074117b45613fd29a72.zip
Adding upstream version 2.10.34.upstream/2.10.34upstream
Signed-off-by: Daniel Baumann <daniel.baumann@progress-linux.org>
Diffstat (limited to 'modules/display-filter-color-blind.c')
-rw-r--r--modules/display-filter-color-blind.c495
1 files changed, 495 insertions, 0 deletions
diff --git a/modules/display-filter-color-blind.c b/modules/display-filter-color-blind.c
new file mode 100644
index 0000000..2d5785c
--- /dev/null
+++ b/modules/display-filter-color-blind.c
@@ -0,0 +1,495 @@
+/* GIMP - The GNU Image Manipulation Program
+ * 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 "libgimpcolor/gimpcolor.h"
+#include "libgimpconfig/gimpconfig.h"
+#include "libgimpmath/gimpmath.h"
+#include "libgimpmodule/gimpmodule.h"
+#include "libgimpwidgets/gimpwidgets.h"
+
+#include "libgimp/libgimp-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 GimpEnumDesc 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
+{
+ GimpColorDisplay parent_instance;
+
+ ColorblindDeficiencyType type;
+
+ gfloat a1, b1, c1;
+ gfloat a2, b2, c2;
+ gfloat inflection;
+};
+
+struct _CdisplayColorblindClass
+{
+ GimpColorDisplayClass parent_instance;
+};
+
+
+enum
+{
+ PROP_0,
+ PROP_TYPE
+};
+
+
+static 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 (GimpColorDisplay *display,
+ GeglBuffer *buffer,
+ GeglRectangle *area);
+static void cdisplay_colorblind_changed (GimpColorDisplay *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 GimpModuleInfo cdisplay_colorblind_info =
+{
+ GIMP_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,
+ GIMP_TYPE_COLOR_DISPLAY)
+
+static GType cdisplay_colorblind_deficiency_type_type = 0;
+
+
+G_MODULE_EXPORT const GimpModuleInfo *
+gimp_module_query (GTypeModule *module)
+{
+ return &cdisplay_colorblind_info;
+}
+
+G_MODULE_EXPORT gboolean
+gimp_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);
+
+ gimp_type_set_translation_domain (cdisplay_colorblind_deficiency_type_type,
+ GETTEXT_PACKAGE "-libgimp");
+ gimp_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);
+ GimpColorDisplayClass *display_class = GIMP_COLOR_DISPLAY_CLASS (klass);
+
+ object_class->get_property = cdisplay_colorblind_get_property;
+ object_class->set_property = cdisplay_colorblind_set_property;
+
+ GIMP_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 = "gimp-colordisplay-colorblind";
+ display_class->icon_name = GIMP_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 (GimpColorDisplay *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 (GimpColorDisplay *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");
+ gimp_color_display_changed (GIMP_COLOR_DISPLAY (colorblind));
+ }
+}