… | |
… | |
2 | #include "perl.h" |
2 | #include "perl.h" |
3 | #include "XSUB.h" |
3 | #include "XSUB.h" |
4 | |
4 | |
5 | #include <string.h> |
5 | #include <string.h> |
6 | #include <setjmp.h> |
6 | #include <setjmp.h> |
|
|
7 | #include <math.h> |
|
|
8 | |
|
|
9 | #include <magic.h> |
7 | |
10 | |
8 | #include <jpeglib.h> |
11 | #include <jpeglib.h> |
9 | #include <glib.h> |
12 | #include <glib.h> |
10 | #include <gtk/gtk.h> |
13 | #include <gtk/gtk.h> |
11 | #include <gdk-pixbuf/gdk-pixbuf.h> |
14 | #include <gdk-pixbuf/gdk-pixbuf.h> |
12 | |
15 | |
13 | #include <gperl.h> |
16 | #include <gperl.h> |
14 | #include <gtk2perl.h> |
17 | #include <gtk2perl.h> |
15 | |
18 | |
|
|
19 | #include <assert.h> |
|
|
20 | |
|
|
21 | #include "perlmulticore.h" |
|
|
22 | |
16 | #define IW 80 /* MUST match Schnauer.pm! */ |
23 | #define IW 80 /* MUST match Schnauzer.pm! */ |
17 | #define IH 60 /* MUST match Schnauer.pm! */ |
24 | #define IH 60 /* MUST match Schnauzer.pm! */ |
18 | |
25 | |
19 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
26 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
20 | |
27 | |
21 | #define LINELENGTH 240 |
28 | #define LINELENGTH 240 |
22 | |
29 | |
23 | #define ELLIPSIS "\xe2\x80\xa6" |
30 | #define ELLIPSIS "\xe2\x80\xa6" |
|
|
31 | |
|
|
32 | typedef char *octet_string; |
24 | |
33 | |
25 | struct jpg_err_mgr |
34 | struct jpg_err_mgr |
26 | { |
35 | { |
27 | struct jpeg_error_mgr err; |
36 | struct jpeg_error_mgr err; |
28 | jmp_buf setjmp_buffer; |
37 | jmp_buf setjmp_buffer; |
… | |
… | |
65 | |
74 | |
66 | *h &= 255; |
75 | *h &= 255; |
67 | } |
76 | } |
68 | } |
77 | } |
69 | |
78 | |
|
|
79 | struct feature { |
|
|
80 | float v1, v2, v3; // mean, square, cube |
|
|
81 | int n; |
|
|
82 | }; |
|
|
83 | |
|
|
84 | static void |
|
|
85 | feature_init (struct feature *f) |
|
|
86 | { |
|
|
87 | f->v1 = 0.; |
|
|
88 | f->v2 = 0.; |
|
|
89 | f->v3 = 0.; |
|
|
90 | f->n = 0; |
|
|
91 | } |
|
|
92 | |
|
|
93 | // didn't find an algorithm to neatly do mean, variance and skew in one pass. |
|
|
94 | // elmex ist schuld. |
|
|
95 | static void |
|
|
96 | feature_update_pass_1 (struct feature *f, unsigned int v) |
|
|
97 | { |
|
|
98 | f->v1 += v; |
|
|
99 | f->n += 1; |
|
|
100 | } |
|
|
101 | |
|
|
102 | static void |
|
|
103 | feature_finish_pass_1 (struct feature *f) |
|
|
104 | { |
|
|
105 | if (f->n < 1) |
|
|
106 | return; |
|
|
107 | |
|
|
108 | f->v1 /= f->n; |
|
|
109 | } |
|
|
110 | |
|
|
111 | static void |
|
|
112 | feature_update_pass_2 (struct feature *f, unsigned int v) |
|
|
113 | { |
|
|
114 | float d = v - f->v1; |
|
|
115 | |
|
|
116 | f->v2 += d * d; |
|
|
117 | f->v3 += d * d * d; |
|
|
118 | } |
|
|
119 | |
|
|
120 | static void |
|
|
121 | feature_finish_pass_2 (struct feature *f) |
|
|
122 | { |
|
|
123 | if (f->n < 1) |
|
|
124 | return; |
|
|
125 | |
|
|
126 | f->v2 /= f->n; |
|
|
127 | f->v3 /= f->n; |
|
|
128 | |
|
|
129 | f->v1 /= 255.; |
|
|
130 | f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2); |
|
|
131 | f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.); |
|
|
132 | } |
|
|
133 | |
70 | static guint32 a85_val; |
134 | static guint32 a85_val; |
71 | static guint a85_cnt; |
135 | static guint a85_cnt; |
72 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
136 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
73 | |
137 | |
74 | static void |
138 | static void |
… | |
… | |
143 | } |
207 | } |
144 | |
208 | |
145 | OUTPUT: |
209 | OUTPUT: |
146 | RETVAL |
210 | RETVAL |
147 | |
211 | |
|
|
212 | const char * |
|
|
213 | magic (octet_string path) |
|
|
214 | CODE: |
|
|
215 | { |
|
|
216 | static magic_t cookie; |
|
|
217 | |
|
|
218 | if (!cookie) |
|
|
219 | { |
|
|
220 | cookie = magic_open (MAGIC_SYMLINK); |
|
|
221 | |
|
|
222 | if (cookie) |
|
|
223 | magic_load (cookie, 0); |
|
|
224 | else |
|
|
225 | XSRETURN_UNDEF; |
|
|
226 | } |
|
|
227 | |
|
|
228 | RETVAL = magic_file (cookie, path); |
|
|
229 | } |
|
|
230 | OUTPUT: |
|
|
231 | RETVAL |
|
|
232 | |
|
|
233 | const char * |
|
|
234 | magic_mime (octet_string path) |
|
|
235 | CODE: |
|
|
236 | { |
|
|
237 | static magic_t cookie; |
|
|
238 | |
|
|
239 | if (!cookie) |
|
|
240 | { |
|
|
241 | cookie = magic_open (MAGIC_MIME | MAGIC_SYMLINK); |
|
|
242 | |
|
|
243 | if (cookie) |
|
|
244 | magic_load (cookie, 0); |
|
|
245 | else |
|
|
246 | XSRETURN_UNDEF; |
|
|
247 | } |
|
|
248 | |
|
|
249 | perlinterp_release (); |
|
|
250 | RETVAL = magic_file (cookie, path); |
|
|
251 | perlinterp_acquire (); |
|
|
252 | } |
|
|
253 | OUTPUT: |
|
|
254 | RETVAL |
|
|
255 | |
148 | # missing in Gtk2 perl module |
256 | # missing/broken in Gtk2 perl module |
|
|
257 | |
|
|
258 | void |
|
|
259 | gdk_window_clear_hints (GdkWindow *window) |
|
|
260 | CODE: |
|
|
261 | gdk_window_set_geometry_hints (window, 0, 0); |
149 | |
262 | |
150 | gboolean |
263 | gboolean |
151 | gdk_net_wm_supports (GdkAtom property) |
264 | gdk_net_wm_supports (GdkAtom property) |
152 | CODE: |
265 | CODE: |
153 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
266 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
… | |
… | |
159 | RETVAL |
272 | RETVAL |
160 | |
273 | |
161 | GdkPixbuf_noinc * |
274 | GdkPixbuf_noinc * |
162 | dealpha_expose (GdkPixbuf *pb) |
275 | dealpha_expose (GdkPixbuf *pb) |
163 | CODE: |
276 | CODE: |
|
|
277 | perlinterp_release (); |
164 | { |
278 | { |
165 | int w = gdk_pixbuf_get_width (pb); |
279 | int w = gdk_pixbuf_get_width (pb); |
166 | int h = gdk_pixbuf_get_height (pb); |
280 | int h = gdk_pixbuf_get_height (pb); |
167 | int bpp = gdk_pixbuf_get_n_channels (pb); |
281 | int bpp = gdk_pixbuf_get_n_channels (pb); |
168 | int x, y, i; |
282 | int x, y, i; |
… | |
… | |
177 | for (x = 0; x < w; x++) |
291 | for (x = 0; x < w; x++) |
178 | for (y = 0; y < h; y++) |
292 | for (y = 0; y < h; y++) |
179 | for (i = 0; i < 3; i++) |
293 | for (i = 0; i < 3; i++) |
180 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
294 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
181 | } |
295 | } |
|
|
296 | perlinterp_acquire (); |
182 | OUTPUT: |
297 | OUTPUT: |
183 | RETVAL |
298 | RETVAL |
184 | |
299 | |
185 | GdkPixbuf_noinc * |
300 | GdkPixbuf_noinc * |
186 | rotate (GdkPixbuf *pb, int angle) |
301 | rotate (GdkPixbuf *pb, int angle) |
187 | CODE: |
302 | CODE: |
|
|
303 | perlinterp_release (); |
|
|
304 | if (angle < 0) |
|
|
305 | angle += 360; |
188 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
306 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
189 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
307 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
190 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
308 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
191 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
309 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
192 | : angle); |
310 | : angle); |
|
|
311 | perlinterp_acquire (); |
193 | OUTPUT: |
312 | OUTPUT: |
194 | RETVAL |
313 | RETVAL |
195 | |
314 | |
196 | GdkPixbuf_noinc * |
315 | GdkPixbuf_noinc * |
197 | load_jpeg (SV *path, int thumbnail=0) |
316 | load_jpeg (SV *path, int thumbnail=0) |
… | |
… | |
209 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
328 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
210 | |
329 | |
211 | if (!fp) |
330 | if (!fp) |
212 | XSRETURN_UNDEF; |
331 | XSRETURN_UNDEF; |
213 | |
332 | |
|
|
333 | perlinterp_release (); |
|
|
334 | |
214 | cinfo.err = jpeg_std_error (&jerr.err); |
335 | cinfo.err = jpeg_std_error (&jerr.err); |
215 | |
336 | |
216 | jerr.err.error_exit = cv_error_exit; |
337 | jerr.err.error_exit = cv_error_exit; |
217 | jerr.err.output_message = cv_error_output; |
338 | jerr.err.output_message = cv_error_output; |
218 | |
339 | |
… | |
… | |
222 | jpeg_destroy_decompress (&cinfo); |
343 | jpeg_destroy_decompress (&cinfo); |
223 | |
344 | |
224 | if (pb) |
345 | if (pb) |
225 | g_object_unref ((gpointer)pb); |
346 | g_object_unref ((gpointer)pb); |
226 | |
347 | |
|
|
348 | perlinterp_acquire (); |
227 | XSRETURN_UNDEF; |
349 | XSRETURN_UNDEF; |
228 | } |
350 | } |
229 | |
351 | |
230 | jpeg_create_decompress (&cinfo); |
352 | jpeg_create_decompress (&cinfo); |
231 | |
353 | |
… | |
… | |
283 | } |
405 | } |
284 | |
406 | |
285 | jpeg_finish_decompress (&cinfo); |
407 | jpeg_finish_decompress (&cinfo); |
286 | fclose (fp); |
408 | fclose (fp); |
287 | jpeg_destroy_decompress (&cinfo); |
409 | jpeg_destroy_decompress (&cinfo); |
|
|
410 | perlinterp_acquire (); |
288 | } |
411 | } |
289 | OUTPUT: |
412 | OUTPUT: |
290 | RETVAL |
413 | RETVAL |
|
|
414 | |
|
|
415 | void |
|
|
416 | compare (GdkPixbuf *a, GdkPixbuf *b) |
|
|
417 | PPCODE: |
|
|
418 | perlinterp_release (); |
|
|
419 | { |
|
|
420 | int w = gdk_pixbuf_get_width (a); |
|
|
421 | int h = gdk_pixbuf_get_height (a); |
|
|
422 | |
|
|
423 | int sa = gdk_pixbuf_get_rowstride (a); |
|
|
424 | int sb = gdk_pixbuf_get_rowstride (b); |
|
|
425 | |
|
|
426 | guchar *pa = gdk_pixbuf_get_pixels (a); |
|
|
427 | guchar *pb = gdk_pixbuf_get_pixels (b); |
|
|
428 | |
|
|
429 | int x, y; |
|
|
430 | |
|
|
431 | assert (w == gdk_pixbuf_get_width (b)); |
|
|
432 | assert (h == gdk_pixbuf_get_height (b)); |
|
|
433 | |
|
|
434 | assert (gdk_pixbuf_get_n_channels (a) == 3); |
|
|
435 | assert (gdk_pixbuf_get_n_channels (b) == 3); |
|
|
436 | |
|
|
437 | double diff = 0.; |
|
|
438 | int peak = 0; |
|
|
439 | |
|
|
440 | if (w && h) |
|
|
441 | for (y = 0; y < h; y++) |
|
|
442 | { |
|
|
443 | guchar *pa_ = pa + y * sa; |
|
|
444 | guchar *pb_ = pb + y * sb; |
|
|
445 | |
|
|
446 | for (x = 0; x < w; x++) |
|
|
447 | { |
|
|
448 | int d; |
|
|
449 | |
|
|
450 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
451 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
452 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
453 | } |
|
|
454 | } |
|
|
455 | |
|
|
456 | perlinterp_acquire (); |
|
|
457 | |
|
|
458 | EXTEND (SP, 2); |
|
|
459 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
|
|
460 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
|
|
461 | } |
291 | |
462 | |
292 | ############################################################################# |
463 | ############################################################################# |
293 | |
464 | |
294 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
465 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
|
|
466 | |
|
|
467 | # currently only works for filenames (octet strings) |
295 | |
468 | |
296 | SV * |
469 | SV * |
297 | foldcase (SV *pathsv) |
470 | foldcase (SV *pathsv) |
298 | PROTOTYPE: $ |
471 | PROTOTYPE: $ |
299 | CODE: |
472 | CODE: |
300 | { |
473 | { |
301 | STRLEN plen; |
474 | STRLEN plen; |
302 | U8 *path = (U8 *)SvPVutf8 (pathsv, plen); |
475 | U8 *path = (U8 *)SvPV (pathsv, plen); |
303 | U8 *pend = path + plen; |
476 | U8 *pend = path + plen; |
304 | U8 dst [plen * 6 * 3], *dstp = dst; |
477 | U8 dst [plen * 6 * 3], *dstp = dst; |
305 | |
478 | |
306 | while (path < pend) |
479 | while (path < pend) |
307 | { |
480 | { |
308 | U8 ch = *path; |
481 | U8 ch = *path; |
309 | |
482 | |
310 | if (ch >= 'a' && ch <= 'z') |
483 | if (ch >= 'a' && ch <= 'z') |
311 | *dstp++ = *path++; |
484 | *dstp++ = *path++; |
|
|
485 | else if (ch >= 'A' && ch <= 'Z') |
|
|
486 | *dstp++ = *path++ + ('a' - 'A'); |
312 | else if (ch >= '0' && ch <= '9') |
487 | else if (ch >= '0' && ch <= '9') |
313 | { |
488 | { |
314 | STRLEN el, nl = 0; |
489 | STRLEN el, nl = 0; |
315 | while (*path >= '0' && *path <= '9' && path < pend) |
490 | while (*path >= '0' && *path <= '9' && path < pend) |
316 | path++, nl++; |
491 | path++, nl++; |
… | |
… | |
319 | *dstp++ = '0'; |
494 | *dstp++ = '0'; |
320 | |
495 | |
321 | memcpy (dstp, path - nl, nl); |
496 | memcpy (dstp, path - nl, nl); |
322 | dstp += nl; |
497 | dstp += nl; |
323 | } |
498 | } |
|
|
499 | else |
|
|
500 | *dstp++ = *path++; |
|
|
501 | #if 0 |
324 | else |
502 | else |
325 | { |
503 | { |
326 | STRLEN cl; |
504 | STRLEN cl; |
327 | to_utf8_fold (path, dstp, &cl); |
505 | to_utf8_fold (path, dstp, &cl); |
328 | dstp += cl; |
506 | dstp += cl; |
329 | path += is_utf8_char (path); |
507 | path += is_utf8_char (path); |
330 | } |
508 | } |
|
|
509 | #endif |
331 | } |
510 | } |
332 | |
511 | |
333 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
512 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
334 | } |
513 | } |
335 | OUTPUT: |
514 | OUTPUT: |
… | |
… | |
497 | ############################################################################# |
676 | ############################################################################# |
498 | |
677 | |
499 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
678 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
500 | |
679 | |
501 | SV * |
680 | SV * |
502 | make_histograms (SV *ar) |
681 | extract_features (SV *ar) |
503 | CODE: |
682 | CODE: |
504 | { |
683 | { |
505 | int i; |
684 | int i; |
506 | AV *av, *result; |
685 | AV *av, *result; |
507 | |
686 | |
508 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
687 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
509 | croak ("Not an array ref as first argument to make_histogram"); |
688 | croak ("Not an array ref as first argument to extract_features"); |
510 | |
689 | |
511 | av = (AV *) SvRV (ar); |
690 | av = (AV *) SvRV (ar); |
512 | result = newAV (); |
691 | result = newAV (); |
513 | |
692 | |
514 | for (i = 0; i <= av_len (av); ++i) |
693 | for (i = 0; i <= av_len (av); ++i) |
515 | { |
694 | { |
516 | const int HISTSIZE = 64; |
|
|
517 | |
|
|
518 | int j; |
|
|
519 | SV *sv = *av_fetch (av, i, 1); |
695 | SV *sv = *av_fetch (av, i, 1); |
520 | STRLEN len; |
|
|
521 | char *buf = SvPVbyte (sv, len); |
|
|
522 | |
|
|
523 | int tmphist[HISTSIZE]; |
|
|
524 | float *hist; |
|
|
525 | |
|
|
526 | SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); |
696 | SV *histsv = newSV (9 * sizeof (float) + 1); |
|
|
697 | |
527 | SvPOK_on (histsv); |
698 | SvPOK_on (histsv); |
528 | SvCUR_set (histsv, HISTSIZE * sizeof (float)); |
699 | SvCUR_set (histsv, 9 * sizeof (float)); |
529 | hist = (float *)SvPVX (histsv); |
700 | float *hist = (float *)SvPVX (histsv); |
530 | |
701 | |
531 | Zero (tmphist, sizeof (tmphist), char); |
702 | struct feature f_h, f_s, f_v; |
|
|
703 | feature_init (&f_h); |
|
|
704 | feature_init (&f_s); |
|
|
705 | feature_init (&f_v); |
532 | |
706 | |
533 | for (j = len; j--; ) |
|
|
534 | { |
|
|
535 | unsigned int idx |
|
|
536 | = ((*buf & 0xc0) >> 2) |
|
|
537 | | ((*buf & 0x18) >> 1) |
|
|
538 | | (*buf & 0x03); |
|
|
539 | |
|
|
540 | ++tmphist[idx]; |
|
|
541 | ++buf; |
|
|
542 | } |
|
|
543 | |
707 | { |
544 | for (j = 0; j < HISTSIZE; ++j) |
708 | STRLEN len; |
545 | hist[j] = (float)tmphist[j] / (len + 1e-30); |
709 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
710 | while (len >= 3) |
|
|
711 | { |
|
|
712 | unsigned int r, g, b, h, s, v; |
|
|
713 | r = *buf++; g = *buf++; b = *buf++; |
|
|
714 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
715 | |
|
|
716 | feature_update_pass_1 (&f_h, h); |
|
|
717 | feature_update_pass_1 (&f_s, s); |
|
|
718 | feature_update_pass_1 (&f_v, v); |
|
|
719 | |
|
|
720 | len -= 3; |
|
|
721 | } |
|
|
722 | |
|
|
723 | feature_finish_pass_1 (&f_h); |
|
|
724 | feature_finish_pass_1 (&f_s); |
|
|
725 | feature_finish_pass_1 (&f_v); |
|
|
726 | } |
|
|
727 | |
|
|
728 | { |
|
|
729 | STRLEN len; |
|
|
730 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
731 | while (len >= 3) |
|
|
732 | { |
|
|
733 | unsigned int r, g, b, h, s, v; |
|
|
734 | r = *buf++; g = *buf++; b = *buf++; |
|
|
735 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
736 | |
|
|
737 | feature_update_pass_2 (&f_h, h); |
|
|
738 | feature_update_pass_2 (&f_s, s); |
|
|
739 | feature_update_pass_2 (&f_v, v); |
|
|
740 | |
|
|
741 | len -= 3; |
|
|
742 | } |
|
|
743 | |
|
|
744 | feature_finish_pass_2 (&f_h); |
|
|
745 | feature_finish_pass_2 (&f_s); |
|
|
746 | feature_finish_pass_2 (&f_v); |
|
|
747 | } |
|
|
748 | |
|
|
749 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
750 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
751 | hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
546 | |
752 | |
547 | av_push (result, histsv); |
753 | av_push (result, histsv); |
548 | } |
754 | } |
549 | |
755 | |
550 | RETVAL = newRV_noinc ((SV *)result); |
756 | RETVAL = newRV_noinc ((SV *)result); |
551 | } |
757 | } |
552 | OUTPUT: |
758 | OUTPUT: |
553 | RETVAL |
759 | RETVAL |
554 | |
760 | |
555 | |
|
|