| 1 | #include "EXTERN.h" |
1 | #include "EXTERN.h" |
| 2 | #include "perl.h" |
2 | #include "perl.h" |
| 3 | #include "XSUB.h" |
3 | #include "XSUB.h" |
| 4 | |
4 | |
|
|
5 | #include <string.h> |
|
|
6 | #include <setjmp.h> |
|
|
7 | #include <math.h> |
|
|
8 | |
|
|
9 | #include <magic.h> |
|
|
10 | |
|
|
11 | #include <jpeglib.h> |
|
|
12 | #include <glib.h> |
|
|
13 | #include <gtk/gtk.h> |
| 5 | #include <gdk-pixbuf/gdk-pixbuf.h> |
14 | #include <gdk-pixbuf/gdk-pixbuf.h> |
| 6 | |
15 | |
| 7 | #include <gperl.h> |
16 | #include <gperl.h> |
| 8 | #include <gtk2perl.h> |
17 | #include <gtk2perl.h> |
| 9 | |
18 | |
| 10 | #define IW 80 |
19 | #include <assert.h> |
|
|
20 | |
|
|
21 | #include "perlmulticore.h" |
|
|
22 | |
|
|
23 | #define IW 80 /* MUST match Schnauzer.pm! */ |
|
|
24 | #define IH 60 /* MUST match Schnauzer.pm! */ |
| 11 | |
25 | |
| 12 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
26 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
|
|
27 | |
|
|
28 | #define LINELENGTH 240 |
|
|
29 | |
|
|
30 | #define ELLIPSIS "\xe2\x80\xa6" |
|
|
31 | |
|
|
32 | typedef char *octet_string; |
|
|
33 | |
|
|
34 | struct jpg_err_mgr |
|
|
35 | { |
|
|
36 | struct jpeg_error_mgr err; |
|
|
37 | jmp_buf setjmp_buffer; |
|
|
38 | }; |
|
|
39 | |
|
|
40 | static void |
|
|
41 | cv_error_exit (j_common_ptr cinfo) |
|
|
42 | { |
|
|
43 | cinfo->err->output_message (cinfo); |
|
|
44 | longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99); |
|
|
45 | } |
|
|
46 | |
|
|
47 | static void |
|
|
48 | cv_error_output (j_common_ptr cinfo) |
|
|
49 | { |
|
|
50 | char msg[JMSG_LENGTH_MAX]; |
|
|
51 | |
|
|
52 | cinfo->err->format_message (cinfo, msg); |
|
|
53 | |
|
|
54 | fprintf (stderr, "JPEG decoding error: %s\n", msg); |
|
|
55 | return; |
|
|
56 | } |
|
|
57 | |
|
|
58 | static void |
|
|
59 | rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, |
|
|
60 | unsigned int *h, unsigned int *s, unsigned int *v) |
|
|
61 | { |
|
|
62 | unsigned int mx = r; if (g > mx) mx = g; if (b > mx) mx = b; |
|
|
63 | unsigned int mn = r; if (g < mn) mn = g; if (b < mn) mn = b; |
|
|
64 | unsigned int delta = mx - mn; |
|
|
65 | |
|
|
66 | *v = mx; |
|
|
67 | |
|
|
68 | *s = mx ? delta * 255 / mx : 0; |
|
|
69 | |
|
|
70 | if (delta == 0) |
|
|
71 | *h = 0; |
|
|
72 | else |
|
|
73 | { |
|
|
74 | if (r == mx) |
|
|
75 | *h = ((int)g - (int)b) * 255 / (int)(delta * 3); |
|
|
76 | else if (g == mx) |
|
|
77 | *h = ((int)b - (int)r) * 255 / (int)(delta * 3) + 52; |
|
|
78 | else if (b == mx) |
|
|
79 | *h = ((int)r - (int)g) * 255 / (int)(delta * 3) + 103; |
|
|
80 | |
|
|
81 | *h &= 255; |
|
|
82 | } |
|
|
83 | } |
|
|
84 | |
|
|
85 | struct feature { |
|
|
86 | float v1, v2, v3; // mean, square, cube |
|
|
87 | int n; |
|
|
88 | }; |
|
|
89 | |
|
|
90 | static void |
|
|
91 | feature_init (struct feature *f) |
|
|
92 | { |
|
|
93 | f->v1 = 0.; |
|
|
94 | f->v2 = 0.; |
|
|
95 | f->v3 = 0.; |
|
|
96 | f->n = 0; |
|
|
97 | } |
|
|
98 | |
|
|
99 | // didn't find an algorithm to neatly do mean, variance and skew in one pass. |
|
|
100 | // elmex ist schuld. |
|
|
101 | static void |
|
|
102 | feature_update_pass_1 (struct feature *f, unsigned int v) |
|
|
103 | { |
|
|
104 | f->v1 += v; |
|
|
105 | f->n += 1; |
|
|
106 | } |
|
|
107 | |
|
|
108 | static void |
|
|
109 | feature_finish_pass_1 (struct feature *f) |
|
|
110 | { |
|
|
111 | if (f->n < 1) |
|
|
112 | return; |
|
|
113 | |
|
|
114 | f->v1 /= f->n; |
|
|
115 | } |
|
|
116 | |
|
|
117 | static void |
|
|
118 | feature_update_pass_2 (struct feature *f, unsigned int v) |
|
|
119 | { |
|
|
120 | float d = v - f->v1; |
|
|
121 | |
|
|
122 | f->v2 += d * d; |
|
|
123 | f->v3 += d * d * d; |
|
|
124 | } |
|
|
125 | |
|
|
126 | static void |
|
|
127 | feature_finish_pass_2 (struct feature *f) |
|
|
128 | { |
|
|
129 | if (f->n < 1) |
|
|
130 | return; |
|
|
131 | |
|
|
132 | f->v2 /= f->n; |
|
|
133 | f->v3 /= f->n; |
|
|
134 | |
|
|
135 | f->v1 /= 255.; |
|
|
136 | f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2); |
|
|
137 | f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.); |
|
|
138 | } |
| 13 | |
139 | |
| 14 | static guint32 a85_val; |
140 | static guint32 a85_val; |
| 15 | static guint a85_cnt; |
141 | static guint a85_cnt; |
| 16 | static guchar a85_buf[80], *a85_ptr; |
142 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
| 17 | |
143 | |
| 18 | static void |
144 | static void |
| 19 | a85_init (void) |
145 | a85_init (void) |
| 20 | { |
146 | { |
| 21 | a85_cnt = 4; |
147 | a85_cnt = 4; |
| … | |
… | |
| 30 | if (!--a85_cnt) |
156 | if (!--a85_cnt) |
| 31 | { |
157 | { |
| 32 | a85_cnt = 4; |
158 | a85_cnt = 4; |
| 33 | if (a85_val) |
159 | if (a85_val) |
| 34 | { |
160 | { |
| 35 | a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
161 | a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
| 36 | a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
162 | a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
| 37 | a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
163 | a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
| 38 | a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
164 | a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
| 39 | a85_ptr[0] = (a85_val ) + 33; |
165 | a85_ptr[0] = (a85_val ) + 33; |
| 40 | |
166 | |
| 41 | a85_ptr += 5; |
167 | a85_ptr += 5; |
| 42 | } |
168 | } |
| … | |
… | |
| 51 | } |
177 | } |
| 52 | } |
178 | } |
| 53 | |
179 | |
| 54 | } |
180 | } |
| 55 | |
181 | |
|
|
182 | static void |
| 56 | a85_finish (PerlIO *fp) |
183 | a85_finish (PerlIO *fp) |
| 57 | { |
184 | { |
| 58 | while (a85_cnt != 4) |
185 | while (a85_cnt != 4) |
| 59 | a85_push (fp, 0); |
186 | a85_push (fp, 0); |
| 60 | |
187 | |
| … | |
… | |
| 63 | *a85_ptr++ = '\n'; |
190 | *a85_ptr++ = '\n'; |
| 64 | |
191 | |
| 65 | PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
192 | PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
| 66 | } |
193 | } |
| 67 | |
194 | |
|
|
195 | ///////////////////////////////////////////////////////////////////////////// |
|
|
196 | |
| 68 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::ImageWindow |
197 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
| 69 | |
198 | |
| 70 | PROTOTYPES: ENABLE |
199 | PROTOTYPES: ENABLE |
| 71 | |
200 | |
|
|
201 | # missing function in perl. really :) |
|
|
202 | int |
|
|
203 | common_prefix_length (a, b) |
|
|
204 | unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
|
|
205 | unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
|
|
206 | CODE: |
|
|
207 | RETVAL = 0; |
|
|
208 | |
|
|
209 | while (*a == *b && *a) |
|
|
210 | { |
|
|
211 | RETVAL += (*a & 0xc0) != 0x80; |
|
|
212 | a++, b++; |
|
|
213 | } |
|
|
214 | |
|
|
215 | OUTPUT: |
|
|
216 | RETVAL |
|
|
217 | |
|
|
218 | const char * |
|
|
219 | magic (octet_string path) |
|
|
220 | CODE: |
|
|
221 | { |
|
|
222 | static magic_t cookie; |
|
|
223 | |
|
|
224 | if (!cookie) |
|
|
225 | { |
|
|
226 | cookie = magic_open (MAGIC_SYMLINK); |
|
|
227 | |
|
|
228 | if (cookie) |
|
|
229 | magic_load (cookie, 0); |
|
|
230 | else |
|
|
231 | XSRETURN_UNDEF; |
|
|
232 | } |
|
|
233 | |
|
|
234 | RETVAL = magic_file (cookie, path); |
|
|
235 | } |
|
|
236 | OUTPUT: |
|
|
237 | RETVAL |
|
|
238 | |
|
|
239 | const char * |
|
|
240 | magic_mime (octet_string path) |
|
|
241 | CODE: |
|
|
242 | { |
|
|
243 | static magic_t cookie; |
|
|
244 | |
|
|
245 | if (!cookie) |
|
|
246 | { |
|
|
247 | cookie = magic_open (MAGIC_MIME | MAGIC_SYMLINK); |
|
|
248 | |
|
|
249 | if (cookie) |
|
|
250 | magic_load (cookie, 0); |
|
|
251 | else |
|
|
252 | XSRETURN_UNDEF; |
|
|
253 | } |
|
|
254 | |
|
|
255 | perlinterp_release (); |
|
|
256 | RETVAL = magic_file (cookie, path); |
|
|
257 | perlinterp_acquire (); |
|
|
258 | } |
|
|
259 | OUTPUT: |
|
|
260 | RETVAL |
|
|
261 | |
|
|
262 | # missing/broken in Gtk2 perl module |
|
|
263 | |
|
|
264 | void |
|
|
265 | gdk_window_clear_hints (GdkWindow *window) |
|
|
266 | CODE: |
|
|
267 | gdk_window_set_geometry_hints (window, 0, 0); |
|
|
268 | |
|
|
269 | gboolean |
|
|
270 | gdk_net_wm_supports (GdkAtom property) |
|
|
271 | CODE: |
|
|
272 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
|
|
273 | RETVAL = gdk_net_wm_supports (property); |
|
|
274 | #else |
|
|
275 | RETVAL = 0; |
|
|
276 | #endif |
|
|
277 | OUTPUT: |
|
|
278 | RETVAL |
|
|
279 | |
| 72 | GdkPixbuf_noinc * |
280 | GdkPixbuf_noinc * |
| 73 | transpose (GdkPixbuf *pb) |
281 | dealpha_expose (GdkPixbuf *pb) |
| 74 | CODE: |
282 | CODE: |
|
|
283 | perlinterp_release (); |
| 75 | { |
284 | { |
| 76 | int w = gdk_pixbuf_get_width (pb); |
285 | int w = gdk_pixbuf_get_width (pb); |
| 77 | int h = gdk_pixbuf_get_height (pb); |
286 | int h = gdk_pixbuf_get_height (pb); |
| 78 | int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3; |
287 | int bpp = gdk_pixbuf_get_n_channels (pb); |
| 79 | int x, y, i; |
288 | int x, y, i; |
| 80 | guchar *src = gdk_pixbuf_get_pixels (pb), *dst; |
289 | guchar *src = gdk_pixbuf_get_pixels (pb), *dst; |
| 81 | int sstr = gdk_pixbuf_get_rowstride (pb), dstr; |
290 | int sstr = gdk_pixbuf_get_rowstride (pb), dstr; |
| 82 | |
291 | |
| 83 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w); |
292 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); |
| 84 | |
293 | |
| 85 | dst = gdk_pixbuf_get_pixels (RETVAL); |
294 | dst = gdk_pixbuf_get_pixels (RETVAL); |
| 86 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
295 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
| 87 | |
296 | |
| 88 | for (y = 0; y < h; y++) |
|
|
| 89 | for (x = 0; x < w; x++) |
297 | for (x = 0; x < w; x++) |
|
|
298 | for (y = 0; y < h; y++) |
| 90 | for (i = 0; i < bpp; i++) |
299 | for (i = 0; i < 3; i++) |
| 91 | dst[y * bpp + x * dstr + i] = src[x * bpp + y * sstr + i]; |
300 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
| 92 | } |
301 | } |
|
|
302 | perlinterp_acquire (); |
| 93 | OUTPUT: |
303 | OUTPUT: |
| 94 | RETVAL |
304 | RETVAL |
| 95 | |
305 | |
| 96 | GdkPixbuf_noinc * |
306 | GdkPixbuf_noinc * |
| 97 | flop (GdkPixbuf *pb) |
307 | rotate (GdkPixbuf *pb, int angle) |
| 98 | CODE: |
308 | CODE: |
| 99 | { |
309 | perlinterp_release (); |
| 100 | int w = gdk_pixbuf_get_width (pb); |
310 | if (angle < 0) |
| 101 | int h = gdk_pixbuf_get_height (pb); |
311 | angle += 360; |
| 102 | int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3; |
312 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
| 103 | int x, y, i; |
313 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
| 104 | guchar *src = gdk_pixbuf_get_pixels (pb), *dst; |
314 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
| 105 | int sstr = gdk_pixbuf_get_rowstride (pb), dstr; |
315 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
| 106 | |
316 | : angle); |
| 107 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h); |
317 | perlinterp_acquire (); |
| 108 | |
|
|
| 109 | dst = gdk_pixbuf_get_pixels (RETVAL); |
|
|
| 110 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
|
|
| 111 | |
|
|
| 112 | for (y = 0; y < h; y++) |
|
|
| 113 | for (x = 0; x < w; x++) |
|
|
| 114 | for (i = 0; i < bpp; i++) |
|
|
| 115 | dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i]; |
|
|
| 116 | } |
|
|
| 117 | OUTPUT: |
318 | OUTPUT: |
| 118 | RETVAL |
319 | RETVAL |
| 119 | |
320 | |
|
|
321 | GdkPixbuf_noinc * |
|
|
322 | load_jpeg (SV *path, int thumbnail = 0, int iw = 0, int ih = 0) |
|
|
323 | CODE: |
|
|
324 | { |
|
|
325 | struct jpeg_decompress_struct cinfo; |
|
|
326 | struct jpg_err_mgr jerr; |
|
|
327 | guchar *data; |
|
|
328 | int rs; |
|
|
329 | FILE *fp; |
|
|
330 | volatile GdkPixbuf *pb = 0; |
|
|
331 | |
|
|
332 | RETVAL = 0; |
|
|
333 | |
|
|
334 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
|
|
335 | |
|
|
336 | if (!fp) |
|
|
337 | XSRETURN_UNDEF; |
|
|
338 | |
|
|
339 | perlinterp_release (); |
|
|
340 | |
|
|
341 | cinfo.err = jpeg_std_error (&jerr.err); |
|
|
342 | |
|
|
343 | jerr.err.error_exit = cv_error_exit; |
|
|
344 | jerr.err.output_message = cv_error_output; |
|
|
345 | |
|
|
346 | if ((rs = setjmp (jerr.setjmp_buffer))) |
|
|
347 | { |
|
|
348 | fclose (fp); |
|
|
349 | jpeg_destroy_decompress (&cinfo); |
|
|
350 | |
|
|
351 | if (pb) |
|
|
352 | g_object_unref ((gpointer)pb); |
|
|
353 | |
|
|
354 | perlinterp_acquire (); |
|
|
355 | XSRETURN_UNDEF; |
|
|
356 | } |
|
|
357 | |
|
|
358 | jpeg_create_decompress (&cinfo); |
|
|
359 | |
|
|
360 | jpeg_stdio_src (&cinfo, fp); |
|
|
361 | jpeg_read_header (&cinfo, TRUE); |
|
|
362 | |
|
|
363 | cinfo.dct_method = JDCT_DEFAULT; |
|
|
364 | cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ |
|
|
365 | cinfo.do_block_smoothing = FALSE; |
|
|
366 | cinfo.out_color_space = JCS_RGB; |
|
|
367 | cinfo.quantize_colors = FALSE; |
|
|
368 | |
|
|
369 | cinfo.scale_num = 1; |
|
|
370 | cinfo.scale_denom = 1; |
|
|
371 | |
|
|
372 | jpeg_calc_output_dimensions (&cinfo); |
|
|
373 | |
|
|
374 | if (thumbnail) |
|
|
375 | { |
|
|
376 | cinfo.dct_method = JDCT_FASTEST; |
|
|
377 | cinfo.do_fancy_upsampling = FALSE; |
|
|
378 | |
|
|
379 | while (cinfo.scale_denom < 8 |
|
|
380 | && cinfo.output_width >= iw*4 |
|
|
381 | && cinfo.output_height >= ih*4) |
|
|
382 | { |
|
|
383 | cinfo.scale_denom <<= 1; |
|
|
384 | jpeg_calc_output_dimensions (&cinfo); |
|
|
385 | } |
|
|
386 | } |
|
|
387 | |
|
|
388 | if (cinfo.output_components != 3) |
|
|
389 | longjmp (jerr.setjmp_buffer, 3); |
|
|
390 | |
|
|
391 | if (cinfo.jpeg_color_space == JCS_YCCK || cinfo.jpeg_color_space == JCS_CMYK) |
|
|
392 | { |
|
|
393 | cinfo.out_color_space = JCS_CMYK; |
|
|
394 | cinfo.output_components = 4; |
|
|
395 | } |
|
|
396 | |
|
|
397 | pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height); |
|
|
398 | if (!RETVAL) |
|
|
399 | longjmp (jerr.setjmp_buffer, 2); |
|
|
400 | |
|
|
401 | data = gdk_pixbuf_get_pixels (RETVAL); |
|
|
402 | rs = gdk_pixbuf_get_rowstride (RETVAL); |
|
|
403 | |
|
|
404 | jpeg_start_decompress (&cinfo); |
|
|
405 | |
|
|
406 | while (cinfo.output_scanline < cinfo.output_height) |
|
|
407 | { |
|
|
408 | int remaining = cinfo.output_height - cinfo.output_scanline; |
|
|
409 | JSAMPROW rp[4]; |
|
|
410 | |
|
|
411 | rp [0] = data + cinfo.output_scanline * rs; |
|
|
412 | rp [1] = (guchar *)rp [0] + rs; |
|
|
413 | rp [2] = (guchar *)rp [1] + rs; |
|
|
414 | rp [3] = (guchar *)rp [2] + rs; |
|
|
415 | |
|
|
416 | jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); |
|
|
417 | } |
|
|
418 | |
|
|
419 | if (cinfo.out_color_space == JCS_CMYK) |
|
|
420 | { |
|
|
421 | guchar *end = data + cinfo.output_height * rs; |
|
|
422 | |
|
|
423 | while (data < end) |
|
|
424 | { |
|
|
425 | U32 c = data [0]; |
|
|
426 | U32 m = data [1]; |
|
|
427 | U32 y = data [2]; |
|
|
428 | U32 k = data [3]; |
|
|
429 | |
|
|
430 | if (0) |
|
|
431 | if (cinfo.Adobe_transform == 2) |
|
|
432 | { |
|
|
433 | c ^= 0xff; |
|
|
434 | m ^= 0xff; |
|
|
435 | y ^= 0xff; |
|
|
436 | k ^= 0xff; |
|
|
437 | } |
|
|
438 | |
|
|
439 | data [0] = (c * k + 0x80) / 0xff; |
|
|
440 | data [1] = (m * k + 0x80) / 0xff; |
|
|
441 | data [2] = (y * k + 0x80) / 0xff; |
|
|
442 | data [3] = 0xff; |
|
|
443 | |
|
|
444 | data += 4; |
|
|
445 | } |
|
|
446 | } |
|
|
447 | |
|
|
448 | jpeg_finish_decompress (&cinfo); |
|
|
449 | fclose (fp); |
|
|
450 | jpeg_destroy_decompress (&cinfo); |
|
|
451 | perlinterp_acquire (); |
|
|
452 | } |
|
|
453 | OUTPUT: |
|
|
454 | RETVAL |
|
|
455 | |
|
|
456 | void |
|
|
457 | compare (GdkPixbuf *a, GdkPixbuf *b) |
|
|
458 | PPCODE: |
|
|
459 | perlinterp_release (); |
|
|
460 | { |
|
|
461 | int w = gdk_pixbuf_get_width (a); |
|
|
462 | int h = gdk_pixbuf_get_height (a); |
|
|
463 | |
|
|
464 | int sa = gdk_pixbuf_get_rowstride (a); |
|
|
465 | int sb = gdk_pixbuf_get_rowstride (b); |
|
|
466 | |
|
|
467 | guchar *pa = gdk_pixbuf_get_pixels (a); |
|
|
468 | guchar *pb = gdk_pixbuf_get_pixels (b); |
|
|
469 | |
|
|
470 | int x, y; |
|
|
471 | |
|
|
472 | assert (w == gdk_pixbuf_get_width (b)); |
|
|
473 | assert (h == gdk_pixbuf_get_height (b)); |
|
|
474 | |
|
|
475 | assert (gdk_pixbuf_get_n_channels (a) == 3); |
|
|
476 | assert (gdk_pixbuf_get_n_channels (b) == 3); |
|
|
477 | |
|
|
478 | double diff = 0.; |
|
|
479 | int peak = 0; |
|
|
480 | |
|
|
481 | if (w && h) |
|
|
482 | for (y = 0; y < h; y++) |
|
|
483 | { |
|
|
484 | guchar *pa_ = pa + y * sa; |
|
|
485 | guchar *pb_ = pb + y * sb; |
|
|
486 | |
|
|
487 | for (x = 0; x < w; x++) |
|
|
488 | { |
|
|
489 | int d; |
|
|
490 | |
|
|
491 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
492 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
493 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
494 | } |
|
|
495 | } |
|
|
496 | |
|
|
497 | perlinterp_acquire (); |
|
|
498 | |
|
|
499 | EXTEND (SP, 2); |
|
|
500 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
|
|
501 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
|
|
502 | } |
|
|
503 | |
|
|
504 | ############################################################################# |
|
|
505 | |
| 120 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
506 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
| 121 | |
507 | |
|
|
508 | # currently only works for filenames (octet strings) |
|
|
509 | |
|
|
510 | SV * |
|
|
511 | foldcase (SV *pathsv) |
|
|
512 | PROTOTYPE: $ |
|
|
513 | CODE: |
|
|
514 | { |
|
|
515 | STRLEN plen; |
|
|
516 | U8 *path = (U8 *)SvPV (pathsv, plen); |
|
|
517 | U8 *pend = path + plen; |
|
|
518 | U8 dst [plen * 6 * 3], *dstp = dst; |
|
|
519 | |
|
|
520 | while (path < pend) |
|
|
521 | { |
|
|
522 | U8 ch = *path; |
|
|
523 | |
|
|
524 | if (ch >= 'a' && ch <= 'z') |
|
|
525 | *dstp++ = *path++; |
|
|
526 | else if (ch >= 'A' && ch <= 'Z') |
|
|
527 | *dstp++ = *path++ + ('a' - 'A'); |
|
|
528 | else if (ch >= '0' && ch <= '9') |
|
|
529 | { |
|
|
530 | STRLEN el, nl = 0; |
|
|
531 | while (*path >= '0' && *path <= '9' && path < pend) |
|
|
532 | path++, nl++; |
|
|
533 | |
|
|
534 | for (el = nl; el < 6; el++) |
|
|
535 | *dstp++ = '0'; |
|
|
536 | |
|
|
537 | memcpy (dstp, path - nl, nl); |
|
|
538 | dstp += nl; |
|
|
539 | } |
|
|
540 | else |
|
|
541 | *dstp++ = *path++; |
|
|
542 | #if 0 |
|
|
543 | else |
|
|
544 | { |
|
|
545 | STRLEN cl; |
|
|
546 | to_utf8_fold (path, dstp, &cl); |
|
|
547 | dstp += cl; |
|
|
548 | path += is_utf8_char (path); |
|
|
549 | } |
|
|
550 | #endif |
|
|
551 | } |
|
|
552 | |
|
|
553 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
|
|
554 | } |
|
|
555 | OUTPUT: |
|
|
556 | RETVAL |
|
|
557 | |
| 122 | GdkPixbuf_noinc * |
558 | GdkPixbuf_noinc * |
| 123 | p7_to_pb (int w, int h, guchar *src) |
559 | p7_to_pb (int w, int h, SV *src_sv) |
|
|
560 | PROTOTYPE: @ |
| 124 | CODE: |
561 | CODE: |
| 125 | { |
562 | { |
| 126 | int x, y; |
563 | int x, y; |
| 127 | guchar *dst, *d; |
564 | guchar *dst, *d; |
| 128 | int dstr; |
565 | int dstr; |
|
|
566 | guchar *src = (guchar *)SvPVbyte_nolen (src_sv); |
| 129 | |
567 | |
| 130 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); |
568 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); |
| 131 | dst = gdk_pixbuf_get_pixels (RETVAL); |
569 | dst = gdk_pixbuf_get_pixels (RETVAL); |
| 132 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
570 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
| 133 | |
571 | |
| … | |
… | |
| 142 | } |
580 | } |
| 143 | } |
581 | } |
| 144 | OUTPUT: |
582 | OUTPUT: |
| 145 | RETVAL |
583 | RETVAL |
| 146 | |
584 | |
|
|
585 | ############################################################################# |
|
|
586 | |
|
|
587 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript |
|
|
588 | |
|
|
589 | void |
|
|
590 | dump_ascii85 (PerlIO *fp, GdkPixbuf *pb) |
|
|
591 | CODE: |
|
|
592 | { |
|
|
593 | int w = gdk_pixbuf_get_width (pb); |
|
|
594 | int h = gdk_pixbuf_get_height (pb); |
|
|
595 | int x, y, i; |
|
|
596 | guchar *dst; |
|
|
597 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
598 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
599 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
600 | |
|
|
601 | a85_init (); |
|
|
602 | |
|
|
603 | for (y = 0; y < h; y++) |
|
|
604 | for (x = 0; x < w; x++) |
|
|
605 | for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
|
|
606 | a85_push (fp, src [x * bpp + y * sstr + i]); |
|
|
607 | |
|
|
608 | a85_finish (fp); |
|
|
609 | } |
|
|
610 | |
|
|
611 | void |
|
|
612 | dump_binary (PerlIO *fp, GdkPixbuf *pb) |
|
|
613 | CODE: |
|
|
614 | { |
|
|
615 | int w = gdk_pixbuf_get_width (pb); |
|
|
616 | int h = gdk_pixbuf_get_height (pb); |
|
|
617 | int x, y, i; |
|
|
618 | guchar *dst; |
|
|
619 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
620 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
621 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
622 | |
|
|
623 | for (y = 0; y < h; y++) |
|
|
624 | for (x = 0; x < w; x++) |
|
|
625 | for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
|
|
626 | PerlIO_putc (fp, src [x * bpp + y * sstr + i]); |
|
|
627 | } |
|
|
628 | |
|
|
629 | ############################################################################# |
|
|
630 | |
|
|
631 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
|
|
632 | |
| 147 | SV * |
633 | SV * |
| 148 | pb_to_p7 (GdkPixbuf *pb) |
634 | pb_to_hv84 (GdkPixbuf *pb) |
| 149 | CODE: |
635 | CODE: |
| 150 | { |
636 | { |
| 151 | int w = gdk_pixbuf_get_width (pb); |
637 | int w = gdk_pixbuf_get_width (pb); |
| 152 | int h = gdk_pixbuf_get_height (pb); |
638 | int h = gdk_pixbuf_get_height (pb); |
| 153 | int x, y; |
639 | int x, y; |
| 154 | guchar *dst; |
640 | guchar *dst; |
| 155 | int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3; |
641 | int bpp = gdk_pixbuf_get_n_channels (pb); |
| 156 | guchar *src = gdk_pixbuf_get_pixels (pb); |
642 | guchar *src = gdk_pixbuf_get_pixels (pb); |
| 157 | int sstr = gdk_pixbuf_get_rowstride (pb); |
643 | int sstr = gdk_pixbuf_get_rowstride (pb); |
| 158 | int Er[IW], Eg[IW], Eb[IW]; |
|
|
| 159 | int seed = 77; |
|
|
| 160 | |
644 | |
| 161 | RETVAL = newSV (w * h); |
645 | RETVAL = newSV (6 * 8 * 12 / 8); |
| 162 | SvPOK_only (RETVAL); |
646 | SvPOK_only (RETVAL); |
| 163 | SvCUR_set (RETVAL, w * h); |
647 | SvCUR_set (RETVAL, 6 * 8 * 12 / 8); |
| 164 | |
648 | |
| 165 | dst = SvPVX (RETVAL); |
649 | dst = (guchar *)SvPVX (RETVAL); |
| 166 | |
|
|
| 167 | memset (Er, 0, sizeof (int) * IW); |
|
|
| 168 | memset (Eg, 0, sizeof (int) * IW); |
|
|
| 169 | memset (Eb, 0, sizeof (int) * IW); |
|
|
| 170 | |
650 | |
| 171 | /* some primitive error distribution + random dithering */ |
651 | /* some primitive error distribution + random dithering */ |
| 172 | |
652 | |
| 173 | for (y = 0; y < h; y++) |
653 | for (y = 0; y < h; y++) |
| 174 | { |
654 | { |
| 175 | int er = 0, eg = 0, eb = 0; |
655 | guchar *p = src + y * sstr; |
| 176 | |
656 | |
| 177 | for (x = 0; x < w; x++) |
657 | for (x = 0; x < w; x += 2) |
| 178 | { |
658 | { |
| 179 | int r, g, b; |
659 | unsigned int r, g, b, h, s, v, H, V1, V2; |
| 180 | guchar *p = src + x * bpp + y * sstr; |
|
|
| 181 | |
660 | |
| 182 | r = ((p[0] + er + Er[x]) * 7 + 128) / 255; |
661 | if (bpp == 3) |
| 183 | g = ((p[1] + eg + Eg[x]) * 7 + 128) / 255; |
662 | r = *p++, g = *p++, b = *p++; |
| 184 | b = ((p[2] + eb + Eb[x]) * 3 + 128) / 255; |
663 | else if (bpp == 1) |
|
|
664 | r = g = b = *p++; |
|
|
665 | else |
|
|
666 | abort (); |
| 185 | |
667 | |
| 186 | r = r > 7 ? 7 : r < 0 ? 0 : r; |
668 | rgb_to_hsv (r, g, b, &h, &s, &v); |
| 187 | g = g > 7 ? 7 : g < 0 ? 0 : g; |
|
|
| 188 | b = b > 3 ? 3 : b < 0 ? 0 : b; |
|
|
| 189 | |
669 | |
| 190 | er += p[0] - (r * 255 + 4) / 7; |
670 | H = (h * 15 / 255) << 4; |
| 191 | eg += p[1] - (g * 255 + 4) / 7; |
671 | V1 = v; |
| 192 | eb += p[2] - (b * 255 + 2) / 3; |
|
|
| 193 | |
672 | |
| 194 | Er[x] = er / 2; er -= er / 2 + RAND % 7 - 3; |
673 | if (bpp == 3) |
| 195 | Eg[x] = eg / 2; eg -= eg / 2 + RAND % 7 - 3; |
674 | r = *p++, g = *p++, b = *p++; |
| 196 | Eb[x] = eb / 2; eb -= eb / 2 + RAND % 7 - 3; |
675 | else if (bpp == 1) |
|
|
676 | r = g = b = *p++; |
|
|
677 | else |
|
|
678 | abort (); |
| 197 | |
679 | |
| 198 | *dst++ = r << 5 | g << 2 | b; |
680 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
681 | |
|
|
682 | H |= h * 15 / 255; |
|
|
683 | V2 = v; |
|
|
684 | |
|
|
685 | *dst++ = H; |
|
|
686 | *dst++ = V1; |
|
|
687 | *dst++ = V2; |
| 199 | } |
688 | } |
| 200 | } |
689 | } |
| 201 | } |
690 | } |
| 202 | OUTPUT: |
691 | OUTPUT: |
| 203 | RETVAL |
692 | RETVAL |
| 204 | |
693 | |
|
|
694 | SV * |
|
|
695 | hv84_to_av (unsigned char *hv84) |
|
|
696 | CODE: |
|
|
697 | { |
|
|
698 | int i = 72 / 3; |
|
|
699 | AV *av = newAV (); |
|
|
700 | |
|
|
701 | RETVAL = (SV *)newRV_noinc ((SV *)av); |
|
|
702 | while (i--) |
|
|
703 | { |
|
|
704 | int h = *hv84++; |
|
|
705 | int v1 = *hv84++; |
|
|
706 | int v2 = *hv84++; |
|
|
707 | |
|
|
708 | av_push (av, newSViv (v1)); |
|
|
709 | av_push (av, newSViv ((h >> 4) * 255 / 15)); |
|
|
710 | av_push (av, newSViv (v2)); |
|
|
711 | av_push (av, newSViv ((h & 15) * 255 / 15)); |
|
|
712 | } |
|
|
713 | } |
|
|
714 | OUTPUT: |
|
|
715 | RETVAL |
|
|
716 | |
|
|
717 | ############################################################################# |
|
|
718 | |
| 205 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript |
719 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
| 206 | |
720 | |
| 207 | void |
721 | SV * |
| 208 | dump_pb (PerlIO *fp, GdkPixbuf *pb) |
722 | extract_features (SV *ar) |
| 209 | CODE: |
723 | CODE: |
| 210 | { |
724 | { |
| 211 | int w = gdk_pixbuf_get_width (pb); |
725 | int i; |
| 212 | int h = gdk_pixbuf_get_height (pb); |
726 | AV *av, *result; |
| 213 | int x, y, i; |
|
|
| 214 | guchar *dst; |
|
|
| 215 | int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3; |
|
|
| 216 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
| 217 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
| 218 | |
727 | |
| 219 | a85_init (); |
728 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
|
|
729 | croak ("Not an array ref as first argument to extract_features"); |
| 220 | |
730 | |
| 221 | for (y = 0; y < h; y++) |
731 | av = (AV *) SvRV (ar); |
| 222 | for (x = 0; x < w; x++) |
732 | result = newAV (); |
| 223 | for (i = 0; i < 3; i++) |
|
|
| 224 | a85_push (fp, src [x * bpp + y * sstr + i]); |
|
|
| 225 | |
733 | |
| 226 | a85_finish (fp); |
734 | for (i = 0; i <= av_len (av); ++i) |
| 227 | } |
735 | { |
|
|
736 | SV *sv = *av_fetch (av, i, 1); |
|
|
737 | SV *histsv = newSV (9 * sizeof (float) + 1); |
| 228 | |
738 | |
|
|
739 | SvPOK_on (histsv); |
|
|
740 | SvCUR_set (histsv, 9 * sizeof (float)); |
|
|
741 | float *hist = (float *)SvPVX (histsv); |
| 229 | |
742 | |
|
|
743 | struct feature f_h, f_s, f_v; |
|
|
744 | feature_init (&f_h); |
|
|
745 | feature_init (&f_s); |
|
|
746 | feature_init (&f_v); |
| 230 | |
747 | |
|
|
748 | { |
|
|
749 | STRLEN len; |
|
|
750 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
751 | while (len >= 3) |
|
|
752 | { |
|
|
753 | unsigned int r, g, b, h, s, v; |
|
|
754 | r = *buf++; g = *buf++; b = *buf++; |
|
|
755 | rgb_to_hsv (r, g, b, &h, &s, &v); |
| 231 | |
756 | |
|
|
757 | feature_update_pass_1 (&f_h, h); |
|
|
758 | feature_update_pass_1 (&f_s, s); |
|
|
759 | feature_update_pass_1 (&f_v, v); |
|
|
760 | |
|
|
761 | len -= 3; |
|
|
762 | } |
|
|
763 | |
|
|
764 | feature_finish_pass_1 (&f_h); |
|
|
765 | feature_finish_pass_1 (&f_s); |
|
|
766 | feature_finish_pass_1 (&f_v); |
|
|
767 | } |
|
|
768 | |
|
|
769 | { |
|
|
770 | STRLEN len; |
|
|
771 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
772 | while (len >= 3) |
|
|
773 | { |
|
|
774 | unsigned int r, g, b, h, s, v; |
|
|
775 | r = *buf++; g = *buf++; b = *buf++; |
|
|
776 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
777 | |
|
|
778 | feature_update_pass_2 (&f_h, h); |
|
|
779 | feature_update_pass_2 (&f_s, s); |
|
|
780 | feature_update_pass_2 (&f_v, v); |
|
|
781 | |
|
|
782 | len -= 3; |
|
|
783 | } |
|
|
784 | |
|
|
785 | feature_finish_pass_2 (&f_h); |
|
|
786 | feature_finish_pass_2 (&f_s); |
|
|
787 | feature_finish_pass_2 (&f_v); |
|
|
788 | } |
|
|
789 | |
|
|
790 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
791 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
792 | hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
|
|
793 | |
|
|
794 | av_push (result, histsv); |
|
|
795 | } |
|
|
796 | |
|
|
797 | RETVAL = newRV_noinc ((SV *)result); |
|
|
798 | } |
|
|
799 | OUTPUT: |
|
|
800 | RETVAL |
|
|
801 | |
