3 * Copyright (C) 2008 The Geeqie Team
5 * Author: Vladimir Nadvornik
6 * based on a patch by Uwe Ohse
8 * This software is released under the GNU General Public License (GNU GPL).
9 * Please read the included file COPYING for more information.
10 * This software comes with no warranty of any kind, use at your own risk!
14 #include "histogram.h"
16 #include "pixbuf_util.h"
21 *----------------------------------------------------------------------------
23 *----------------------------------------------------------------------------
26 #define HISTOGRAM_SIZE 256
29 gulong histmap[HISTOGRAM_SIZE*4];
31 gint histogram_logmode;
35 Histogram *histogram_new()
39 histogram = g_new0(Histogram, 1);
40 histogram->histogram_chan = HCHAN_RGB;
41 histogram->histogram_logmode = 1;
46 void histogram_free(Histogram *histogram)
52 gint histogram_set_channel(Histogram *histogram, gint chan)
54 if (!histogram) return 0;
55 histogram->histogram_chan = chan;
59 gint histogram_get_channel(Histogram *histogram)
61 if (!histogram) return 0;
62 return histogram->histogram_chan;
65 gint histogram_set_mode(Histogram *histogram, gint mode)
67 if (!histogram) return 0;
68 histogram->histogram_logmode = mode;
72 gint histogram_get_mode(Histogram *histogram)
74 if (!histogram) return 0;
75 return histogram->histogram_logmode;
78 const gchar *histogram_label(Histogram *histogram)
81 if (!histogram) return NULL;
83 if (histogram->histogram_logmode)
84 switch (histogram->histogram_chan)
86 case HCHAN_R: t1 = _("logarithmical histogram on red"); break;
87 case HCHAN_G: t1 = _("logarithmical histogram on green"); break;
88 case HCHAN_B: t1 = _("logarithmical histogram on blue"); break;
89 case HCHAN_VAL: t1 = _("logarithmical histogram on value"); break;
90 case HCHAN_RGB: t1 = _("logarithmical histogram on RGB"); break;
91 case HCHAN_MAX: t1 = _("logarithmical histogram on max value"); break;
94 switch (histogram->histogram_chan)
96 case HCHAN_R: t1 = _("linear histogram on red"); break;
97 case HCHAN_G: t1 = _("linear histogram on green"); break;
98 case HCHAN_B: t1 = _("linear histogram on blue"); break;
99 case HCHAN_VAL: t1 = _("linear histogram on value"); break;
100 case HCHAN_RGB: t1 = _("linear histogram on RGB"); break;
101 case HCHAN_MAX: t1 = _("linear histogram on max value"); break;
106 gulong histogram_read(Histogram *histogram, GdkPixbuf *imgpixbuf)
108 gint w, h, i, j, srs, has_alpha, step;
111 if (!histogram) return 0;
113 w = gdk_pixbuf_get_width(imgpixbuf);
114 h = gdk_pixbuf_get_height(imgpixbuf);
115 srs = gdk_pixbuf_get_rowstride(imgpixbuf);
116 s_pix = gdk_pixbuf_get_pixels(imgpixbuf);
117 has_alpha = gdk_pixbuf_get_has_alpha(imgpixbuf);
119 memset(histogram->histmap, 0, sizeof(histogram->histmap));
121 /* code duplication is here to speed up the calculation */
122 step = 3 + !!(has_alpha);
123 if (histogram->histogram_chan == HCHAN_MAX)
125 for (i = 0; i < h; i++)
127 guchar *sp = s_pix + (i * srs); /* 8bit */
128 for (j = 0; j < w; j++)
131 if (sp[1]>t) t = sp[1];
132 if (sp[2]>t) t = sp[2];
134 histogram->histmap[sp[0] + 0 * HISTOGRAM_SIZE]++;
135 histogram->histmap[sp[1] + 1 * HISTOGRAM_SIZE]++;
136 histogram->histmap[sp[2] + 2 * HISTOGRAM_SIZE]++;
137 histogram->histmap[t + 3 * HISTOGRAM_SIZE]++;
144 for (i = 0; i < h; i++)
146 guchar *sp = s_pix + (i * srs); /* 8bit */
147 for (j = 0; j < w; j++)
149 histogram->histmap[sp[0] + 0 * HISTOGRAM_SIZE]++;
150 histogram->histmap[sp[1] + 1 * HISTOGRAM_SIZE]++;
151 histogram->histmap[sp[2] + 2 * HISTOGRAM_SIZE]++;
152 histogram->histmap[3 * HISTOGRAM_SIZE + (sp[0]+sp[1]+sp[2])/3]++;
161 gint histogram_draw(Histogram *histogram, GdkPixbuf *pixbuf, gint x, gint y, gint width, gint height)
163 /* FIXME: use the coordinates correctly */
168 if (!histogram) return 0;
170 for (i = 0; i < 1024; i++) {
173 switch (histogram->histogram_chan)
175 case HCHAN_RGB: if ((i%4) != 3) flag = 1; break;
176 case HCHAN_R: if ((i%4) == 0) flag = 1; break;
177 case HCHAN_G: if ((i%4) == 1) flag = 1; break;
178 case HCHAN_B: if ((i%4) == 2) flag = 1; break;
179 case HCHAN_VAL: if ((i%4) == 3) flag = 1; break;
180 case HCHAN_MAX: if ((i%4) == 3) flag = 1; break;
182 if (flag && histogram->histmap[i] > max) max = histogram->histmap[i];
186 for (i = 0; i < width; i++)
189 glong v[4] = {0, 0, 0, 0};
193 gint ii = i * HISTOGRAM_SIZE / width;
194 gint combine = (HISTOGRAM_SIZE - 1) / width + 1;
196 for (j = 0; j < combine; j++)
198 v[0] += histogram->histmap[ii + j + 0 * HISTOGRAM_SIZE]; // r
199 v[1] += histogram->histmap[ii + j + 1 * HISTOGRAM_SIZE]; // g
200 v[2] += histogram->histmap[ii + j + 2 * HISTOGRAM_SIZE]; // b
201 v[3] += histogram->histmap[ii + j + 3 * HISTOGRAM_SIZE]; // value, max
204 for (j = 0; j < 4; j++)
209 for (k = 1; k < 4; k++)
210 if (v[k] > v[max2]) max2 = k;
212 if (histogram->histogram_chan >= HCHAN_RGB
213 || max2 == histogram->histogram_chan)
222 case HCHAN_R: rplus = r = 255; break;
223 case HCHAN_G: gplus = g = 255; break;
224 case HCHAN_B: bplus = b = 255; break;
227 switch (histogram->histogram_chan)
230 if (r == 255 && g == 255 && b == 255)
235 case HCHAN_R: b = 0; g = 0; break;
236 case HCHAN_G: r = 0; b = 0; break;
237 case HCHAN_B: r = 0; g = 0; break;
239 case HCHAN_VAL: r = 0; b = 0; g = 0; break;
244 else if (histogram->histogram_logmode)
245 pt = ((float)log(v[max2])) / logmax * (height - 1);
247 pt = ((float)v[max2])/ max * (height - 1);
249 pixbuf_draw_line(pixbuf,
251 x + i, y + height, x + i, y + height - pt,