… | |
… | |
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> |
7 | |
8 | |
8 | #include <magic.h> |
9 | #include <magic.h> |
9 | |
10 | |
10 | #include <jpeglib.h> |
11 | #include <jpeglib.h> |
11 | #include <glib.h> |
12 | #include <glib.h> |
… | |
… | |
13 | #include <gdk-pixbuf/gdk-pixbuf.h> |
14 | #include <gdk-pixbuf/gdk-pixbuf.h> |
14 | |
15 | |
15 | #include <gperl.h> |
16 | #include <gperl.h> |
16 | #include <gtk2perl.h> |
17 | #include <gtk2perl.h> |
17 | |
18 | |
|
|
19 | #include <assert.h> |
|
|
20 | |
18 | #define IW 80 /* MUST match Schnauer.pm! */ |
21 | #define IW 80 /* MUST match Schnauzer.pm! */ |
19 | #define IH 60 /* MUST match Schnauer.pm! */ |
22 | #define IH 60 /* MUST match Schnauzer.pm! */ |
20 | |
23 | |
21 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
24 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
22 | |
25 | |
23 | #define LINELENGTH 240 |
26 | #define LINELENGTH 240 |
24 | |
27 | |
25 | #define ELLIPSIS "\xe2\x80\xa6" |
28 | #define ELLIPSIS "\xe2\x80\xa6" |
|
|
29 | |
|
|
30 | typedef char *octet_string; |
26 | |
31 | |
27 | struct jpg_err_mgr |
32 | struct jpg_err_mgr |
28 | { |
33 | { |
29 | struct jpeg_error_mgr err; |
34 | struct jpeg_error_mgr err; |
30 | jmp_buf setjmp_buffer; |
35 | jmp_buf setjmp_buffer; |
… | |
… | |
67 | |
72 | |
68 | *h &= 255; |
73 | *h &= 255; |
69 | } |
74 | } |
70 | } |
75 | } |
71 | |
76 | |
|
|
77 | struct feature { |
|
|
78 | float v1, v2, v3; // mean, square, cube |
|
|
79 | int n; |
|
|
80 | }; |
|
|
81 | |
|
|
82 | static void |
|
|
83 | feature_init (struct feature *f) |
|
|
84 | { |
|
|
85 | f->v1 = 0.; |
|
|
86 | f->v2 = 0.; |
|
|
87 | f->v3 = 0.; |
|
|
88 | f->n = 0; |
|
|
89 | } |
|
|
90 | |
|
|
91 | // didn't find an algorithm to neatly do mean, variance and skew in one pass. |
|
|
92 | // elmex ist schuld. |
|
|
93 | static void |
|
|
94 | feature_update_pass_1 (struct feature *f, unsigned int v) |
|
|
95 | { |
|
|
96 | f->v1 += v; |
|
|
97 | f->n += 1; |
|
|
98 | } |
|
|
99 | |
|
|
100 | static void |
|
|
101 | feature_finish_pass_1 (struct feature *f) |
|
|
102 | { |
|
|
103 | if (f->n < 1) |
|
|
104 | return; |
|
|
105 | |
|
|
106 | f->v1 /= f->n; |
|
|
107 | } |
|
|
108 | |
|
|
109 | static void |
|
|
110 | feature_update_pass_2 (struct feature *f, unsigned int v) |
|
|
111 | { |
|
|
112 | float d = v - f->v1; |
|
|
113 | |
|
|
114 | f->v2 += d * d; |
|
|
115 | f->v3 += d * d * d; |
|
|
116 | } |
|
|
117 | |
|
|
118 | static void |
|
|
119 | feature_finish_pass_2 (struct feature *f) |
|
|
120 | { |
|
|
121 | if (f->n < 1) |
|
|
122 | return; |
|
|
123 | |
|
|
124 | f->v2 /= f->n; |
|
|
125 | f->v3 /= f->n; |
|
|
126 | |
|
|
127 | f->v1 /= 255.; |
|
|
128 | f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2); |
|
|
129 | f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.); |
|
|
130 | } |
|
|
131 | |
72 | static guint32 a85_val; |
132 | static guint32 a85_val; |
73 | static guint a85_cnt; |
133 | static guint a85_cnt; |
74 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
134 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
75 | |
135 | |
76 | static void |
136 | static void |
… | |
… | |
146 | |
206 | |
147 | OUTPUT: |
207 | OUTPUT: |
148 | RETVAL |
208 | RETVAL |
149 | |
209 | |
150 | const char * |
210 | const char * |
151 | magic (const char *path) |
211 | magic (octet_string path) |
152 | CODE: |
212 | CODE: |
153 | { |
213 | { |
154 | static magic_t cookie; |
214 | static magic_t cookie; |
155 | |
215 | |
156 | if (!cookie) |
216 | if (!cookie) |
… | |
… | |
167 | } |
227 | } |
168 | OUTPUT: |
228 | OUTPUT: |
169 | RETVAL |
229 | RETVAL |
170 | |
230 | |
171 | const char * |
231 | const char * |
172 | magic_mime (const char *path) |
232 | magic_mime (octet_string path) |
173 | CODE: |
233 | CODE: |
174 | { |
234 | { |
175 | static magic_t cookie; |
235 | static magic_t cookie; |
176 | |
236 | |
177 | if (!cookie) |
237 | if (!cookie) |
… | |
… | |
187 | RETVAL = magic_file (cookie, path); |
247 | RETVAL = magic_file (cookie, path); |
188 | } |
248 | } |
189 | OUTPUT: |
249 | OUTPUT: |
190 | RETVAL |
250 | RETVAL |
191 | |
251 | |
192 | # missing in Gtk2 perl module |
252 | # missing/broken in Gtk2 perl module |
|
|
253 | |
|
|
254 | void |
|
|
255 | gdk_window_clear_hints (GdkWindow *window) |
|
|
256 | CODE: |
|
|
257 | gdk_window_set_geometry_hints (window, 0, 0); |
193 | |
258 | |
194 | gboolean |
259 | gboolean |
195 | gdk_net_wm_supports (GdkAtom property) |
260 | gdk_net_wm_supports (GdkAtom property) |
196 | CODE: |
261 | CODE: |
197 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
262 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
… | |
… | |
227 | RETVAL |
292 | RETVAL |
228 | |
293 | |
229 | GdkPixbuf_noinc * |
294 | GdkPixbuf_noinc * |
230 | rotate (GdkPixbuf *pb, int angle) |
295 | rotate (GdkPixbuf *pb, int angle) |
231 | CODE: |
296 | CODE: |
|
|
297 | if (angle < 0) |
|
|
298 | angle += 360; |
232 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
299 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
233 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
300 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
234 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
301 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
235 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
302 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
236 | : angle); |
303 | : angle); |
… | |
… | |
331 | jpeg_destroy_decompress (&cinfo); |
398 | jpeg_destroy_decompress (&cinfo); |
332 | } |
399 | } |
333 | OUTPUT: |
400 | OUTPUT: |
334 | RETVAL |
401 | RETVAL |
335 | |
402 | |
|
|
403 | void |
|
|
404 | compare (GdkPixbuf *a, GdkPixbuf *b) |
|
|
405 | PPCODE: |
|
|
406 | { |
|
|
407 | int w = gdk_pixbuf_get_width (a); |
|
|
408 | int h = gdk_pixbuf_get_height (a); |
|
|
409 | |
|
|
410 | int sa = gdk_pixbuf_get_rowstride (a); |
|
|
411 | int sb = gdk_pixbuf_get_rowstride (b); |
|
|
412 | |
|
|
413 | guchar *pa = gdk_pixbuf_get_pixels (a); |
|
|
414 | guchar *pb = gdk_pixbuf_get_pixels (b); |
|
|
415 | |
|
|
416 | int x, y; |
|
|
417 | |
|
|
418 | assert (w == gdk_pixbuf_get_width (b)); |
|
|
419 | assert (h == gdk_pixbuf_get_height (b)); |
|
|
420 | |
|
|
421 | assert (gdk_pixbuf_get_n_channels (a) == 3); |
|
|
422 | assert (gdk_pixbuf_get_n_channels (b) == 3); |
|
|
423 | |
|
|
424 | double diff = 0.; |
|
|
425 | int peak = 0; |
|
|
426 | |
|
|
427 | if (w && h) |
|
|
428 | for (y = 0; y < h; y++) |
|
|
429 | { |
|
|
430 | guchar *pa_ = pa + y * sa; |
|
|
431 | guchar *pb_ = pb + y * sb; |
|
|
432 | |
|
|
433 | for (x = 0; x < w; x++) |
|
|
434 | { |
|
|
435 | int d; |
|
|
436 | |
|
|
437 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
438 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
439 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
440 | } |
|
|
441 | } |
|
|
442 | |
|
|
443 | EXTEND (SP, 2); |
|
|
444 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
|
|
445 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
|
|
446 | } |
|
|
447 | |
336 | ############################################################################# |
448 | ############################################################################# |
337 | |
449 | |
338 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
450 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
|
|
451 | |
|
|
452 | # currently only works for filenames (octet strings) |
339 | |
453 | |
340 | SV * |
454 | SV * |
341 | foldcase (SV *pathsv) |
455 | foldcase (SV *pathsv) |
342 | PROTOTYPE: $ |
456 | PROTOTYPE: $ |
343 | CODE: |
457 | CODE: |
344 | { |
458 | { |
345 | STRLEN plen; |
459 | STRLEN plen; |
346 | U8 *path = (U8 *)SvPVutf8 (pathsv, plen); |
460 | U8 *path = (U8 *)SvPV (pathsv, plen); |
347 | U8 *pend = path + plen; |
461 | U8 *pend = path + plen; |
348 | U8 dst [plen * 6 * 3], *dstp = dst; |
462 | U8 dst [plen * 6 * 3], *dstp = dst; |
349 | |
463 | |
350 | while (path < pend) |
464 | while (path < pend) |
351 | { |
465 | { |
352 | U8 ch = *path; |
466 | U8 ch = *path; |
353 | |
467 | |
354 | if (ch >= 'a' && ch <= 'z') |
468 | if (ch >= 'a' && ch <= 'z') |
355 | *dstp++ = *path++; |
469 | *dstp++ = *path++; |
|
|
470 | else if (ch >= 'A' && ch <= 'Z') |
|
|
471 | *dstp++ = *path++ + ('a' - 'A'); |
356 | else if (ch >= '0' && ch <= '9') |
472 | else if (ch >= '0' && ch <= '9') |
357 | { |
473 | { |
358 | STRLEN el, nl = 0; |
474 | STRLEN el, nl = 0; |
359 | while (*path >= '0' && *path <= '9' && path < pend) |
475 | while (*path >= '0' && *path <= '9' && path < pend) |
360 | path++, nl++; |
476 | path++, nl++; |
… | |
… | |
363 | *dstp++ = '0'; |
479 | *dstp++ = '0'; |
364 | |
480 | |
365 | memcpy (dstp, path - nl, nl); |
481 | memcpy (dstp, path - nl, nl); |
366 | dstp += nl; |
482 | dstp += nl; |
367 | } |
483 | } |
|
|
484 | else |
|
|
485 | *dstp++ = *path++; |
|
|
486 | #if 0 |
368 | else |
487 | else |
369 | { |
488 | { |
370 | STRLEN cl; |
489 | STRLEN cl; |
371 | to_utf8_fold (path, dstp, &cl); |
490 | to_utf8_fold (path, dstp, &cl); |
372 | dstp += cl; |
491 | dstp += cl; |
373 | path += is_utf8_char (path); |
492 | path += is_utf8_char (path); |
374 | } |
493 | } |
|
|
494 | #endif |
375 | } |
495 | } |
376 | |
496 | |
377 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
497 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
378 | } |
498 | } |
379 | OUTPUT: |
499 | OUTPUT: |
… | |
… | |
541 | ############################################################################# |
661 | ############################################################################# |
542 | |
662 | |
543 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
663 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
544 | |
664 | |
545 | SV * |
665 | SV * |
546 | make_histograms (SV *ar) |
666 | extract_features (SV *ar) |
547 | CODE: |
667 | CODE: |
548 | { |
668 | { |
549 | int i; |
669 | int i; |
550 | AV *av, *result; |
670 | AV *av, *result; |
551 | |
671 | |
552 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
672 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
553 | croak ("Not an array ref as first argument to make_histogram"); |
673 | croak ("Not an array ref as first argument to extract_features"); |
554 | |
674 | |
555 | av = (AV *) SvRV (ar); |
675 | av = (AV *) SvRV (ar); |
556 | result = newAV (); |
676 | result = newAV (); |
557 | |
677 | |
558 | for (i = 0; i <= av_len (av); ++i) |
678 | for (i = 0; i <= av_len (av); ++i) |
559 | { |
679 | { |
560 | const int HISTSIZE = 64; |
|
|
561 | |
|
|
562 | int j; |
|
|
563 | SV *sv = *av_fetch (av, i, 1); |
680 | SV *sv = *av_fetch (av, i, 1); |
564 | STRLEN len; |
|
|
565 | char *buf = SvPVbyte (sv, len); |
|
|
566 | |
|
|
567 | int tmphist[HISTSIZE]; |
|
|
568 | float *hist; |
|
|
569 | |
|
|
570 | SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); |
681 | SV *histsv = newSV (9 * sizeof (float) + 1); |
|
|
682 | |
571 | SvPOK_on (histsv); |
683 | SvPOK_on (histsv); |
572 | SvCUR_set (histsv, HISTSIZE * sizeof (float)); |
684 | SvCUR_set (histsv, 9 * sizeof (float)); |
573 | hist = (float *)SvPVX (histsv); |
685 | float *hist = (float *)SvPVX (histsv); |
574 | |
686 | |
575 | Zero (tmphist, sizeof (tmphist), char); |
687 | struct feature f_h, f_s, f_v; |
|
|
688 | feature_init (&f_h); |
|
|
689 | feature_init (&f_s); |
|
|
690 | feature_init (&f_v); |
576 | |
691 | |
577 | for (j = len; j--; ) |
|
|
578 | { |
|
|
579 | unsigned int idx |
|
|
580 | = ((*buf & 0xc0) >> 2) |
|
|
581 | | ((*buf & 0x18) >> 1) |
|
|
582 | | (*buf & 0x03); |
|
|
583 | |
|
|
584 | ++tmphist[idx]; |
|
|
585 | ++buf; |
|
|
586 | } |
|
|
587 | |
692 | { |
588 | for (j = 0; j < HISTSIZE; ++j) |
693 | STRLEN len; |
589 | hist[j] = (float)tmphist[j] / (len + 1e-30); |
694 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
695 | while (len >= 3) |
|
|
696 | { |
|
|
697 | unsigned int r, g, b, h, s, v; |
|
|
698 | r = *buf++; g = *buf++; b = *buf++; |
|
|
699 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
700 | |
|
|
701 | feature_update_pass_1 (&f_h, h); |
|
|
702 | feature_update_pass_1 (&f_s, s); |
|
|
703 | feature_update_pass_1 (&f_v, v); |
|
|
704 | |
|
|
705 | len -= 3; |
|
|
706 | } |
|
|
707 | |
|
|
708 | feature_finish_pass_1 (&f_h); |
|
|
709 | feature_finish_pass_1 (&f_s); |
|
|
710 | feature_finish_pass_1 (&f_v); |
|
|
711 | } |
|
|
712 | |
|
|
713 | { |
|
|
714 | STRLEN len; |
|
|
715 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
716 | while (len >= 3) |
|
|
717 | { |
|
|
718 | unsigned int r, g, b, h, s, v; |
|
|
719 | r = *buf++; g = *buf++; b = *buf++; |
|
|
720 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
721 | |
|
|
722 | feature_update_pass_2 (&f_h, h); |
|
|
723 | feature_update_pass_2 (&f_s, s); |
|
|
724 | feature_update_pass_2 (&f_v, v); |
|
|
725 | |
|
|
726 | len -= 3; |
|
|
727 | } |
|
|
728 | |
|
|
729 | feature_finish_pass_2 (&f_h); |
|
|
730 | feature_finish_pass_2 (&f_s); |
|
|
731 | feature_finish_pass_2 (&f_v); |
|
|
732 | } |
|
|
733 | |
|
|
734 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
735 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
736 | hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
590 | |
737 | |
591 | av_push (result, histsv); |
738 | av_push (result, histsv); |
592 | } |
739 | } |
593 | |
740 | |
594 | RETVAL = newRV_noinc ((SV *)result); |
741 | RETVAL = newRV_noinc ((SV *)result); |
595 | } |
742 | } |
596 | OUTPUT: |
743 | OUTPUT: |
597 | RETVAL |
744 | RETVAL |
598 | |
745 | |
599 | |
|
|