… | |
… | |
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 | |
16 | #define IW 80 /* MUST match Schnauer.pm! */ |
19 | #define IW 80 /* MUST match Schnauzer.pm! */ |
17 | #define IH 60 /* MUST match Schnauer.pm! */ |
20 | #define IH 60 /* MUST match Schnauzer.pm! */ |
18 | |
21 | |
19 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
22 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
20 | |
23 | |
21 | #define LINELENGTH 240 |
24 | #define LINELENGTH 240 |
22 | |
25 | |
… | |
… | |
65 | |
68 | |
66 | *h &= 255; |
69 | *h &= 255; |
67 | } |
70 | } |
68 | } |
71 | } |
69 | |
72 | |
|
|
73 | struct feature { |
|
|
74 | float v1, v2, v3; // mean, square, cube |
|
|
75 | int n; |
|
|
76 | }; |
|
|
77 | |
|
|
78 | static void |
|
|
79 | feature_init (struct feature *f) |
|
|
80 | { |
|
|
81 | f->v1 = 0.; |
|
|
82 | f->v2 = 0.; |
|
|
83 | f->v3 = 0.; |
|
|
84 | f->n = 0; |
|
|
85 | } |
|
|
86 | |
|
|
87 | // didn't find an algorithm to neatly do mean, variance and skew in one pass. |
|
|
88 | // elmex ist schuld. |
|
|
89 | static void |
|
|
90 | feature_update_pass_1 (struct feature *f, unsigned int v) |
|
|
91 | { |
|
|
92 | f->v1 += v; |
|
|
93 | f->n += 1; |
|
|
94 | } |
|
|
95 | |
|
|
96 | static void |
|
|
97 | feature_finish_pass_1 (struct feature *f) |
|
|
98 | { |
|
|
99 | if (f->n < 1) |
|
|
100 | return; |
|
|
101 | |
|
|
102 | f->v1 /= f->n; |
|
|
103 | } |
|
|
104 | |
|
|
105 | static void |
|
|
106 | feature_update_pass_2 (struct feature *f, unsigned int v) |
|
|
107 | { |
|
|
108 | float d = v - f->v1; |
|
|
109 | |
|
|
110 | f->v2 += d * d; |
|
|
111 | f->v3 += d * d * d; |
|
|
112 | } |
|
|
113 | |
|
|
114 | static void |
|
|
115 | feature_finish_pass_2 (struct feature *f) |
|
|
116 | { |
|
|
117 | if (f->n < 1) |
|
|
118 | return; |
|
|
119 | |
|
|
120 | f->v2 /= f->n; |
|
|
121 | f->v3 /= f->n; |
|
|
122 | |
|
|
123 | // this is a peculiar weightung, but it works fine |
|
|
124 | f->v1 /= 255.; |
|
|
125 | f->v2 /= 128. * 128.; |
|
|
126 | f->v3 /= 128. * 128. * 64.; |
|
|
127 | } |
|
|
128 | |
70 | static guint32 a85_val; |
129 | static guint32 a85_val; |
71 | static guint a85_cnt; |
130 | static guint a85_cnt; |
72 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
131 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
73 | |
132 | |
74 | static void |
133 | static void |
… | |
… | |
143 | } |
202 | } |
144 | |
203 | |
145 | OUTPUT: |
204 | OUTPUT: |
146 | RETVAL |
205 | RETVAL |
147 | |
206 | |
|
|
207 | const char * |
|
|
208 | magic (const char *path) |
|
|
209 | CODE: |
|
|
210 | { |
|
|
211 | static magic_t cookie; |
|
|
212 | |
|
|
213 | if (!cookie) |
|
|
214 | { |
|
|
215 | cookie = magic_open (MAGIC_NONE); |
|
|
216 | |
|
|
217 | if (cookie) |
|
|
218 | magic_load (cookie, 0); |
|
|
219 | else |
|
|
220 | XSRETURN_UNDEF; |
|
|
221 | } |
|
|
222 | |
|
|
223 | RETVAL = magic_file (cookie, path); |
|
|
224 | } |
|
|
225 | OUTPUT: |
|
|
226 | RETVAL |
|
|
227 | |
|
|
228 | const char * |
|
|
229 | magic_mime (const char *path) |
|
|
230 | CODE: |
|
|
231 | { |
|
|
232 | static magic_t cookie; |
|
|
233 | |
|
|
234 | if (!cookie) |
|
|
235 | { |
|
|
236 | cookie = magic_open (MAGIC_MIME); |
|
|
237 | |
|
|
238 | if (cookie) |
|
|
239 | magic_load (cookie, 0); |
|
|
240 | else |
|
|
241 | XSRETURN_UNDEF; |
|
|
242 | } |
|
|
243 | |
|
|
244 | RETVAL = magic_file (cookie, path); |
|
|
245 | } |
|
|
246 | OUTPUT: |
|
|
247 | RETVAL |
|
|
248 | |
148 | # missing in Gtk2 perl module |
249 | # missing in Gtk2 perl module |
149 | |
250 | |
150 | gboolean |
251 | gboolean |
151 | gdk_net_wm_supports (GdkAtom property) |
252 | gdk_net_wm_supports (GdkAtom property) |
152 | CODE: |
253 | CODE: |
… | |
… | |
201 | struct jpg_err_mgr jerr; |
302 | struct jpg_err_mgr jerr; |
202 | guchar *data; |
303 | guchar *data; |
203 | int rs; |
304 | int rs; |
204 | FILE *fp; |
305 | FILE *fp; |
205 | volatile GdkPixbuf *pb = 0; |
306 | volatile GdkPixbuf *pb = 0; |
206 | gchar *filename; |
|
|
207 | |
307 | |
208 | RETVAL = 0; |
308 | RETVAL = 0; |
209 | |
309 | |
210 | filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0); |
|
|
211 | fp = fopen (filename, "rb"); |
310 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
212 | g_free (filename); |
|
|
213 | |
311 | |
214 | if (!fp) |
312 | if (!fp) |
215 | XSRETURN_UNDEF; |
313 | XSRETURN_UNDEF; |
216 | |
314 | |
217 | cinfo.err = jpeg_std_error (&jerr.err); |
315 | cinfo.err = jpeg_std_error (&jerr.err); |
… | |
… | |
289 | fclose (fp); |
387 | fclose (fp); |
290 | jpeg_destroy_decompress (&cinfo); |
388 | jpeg_destroy_decompress (&cinfo); |
291 | } |
389 | } |
292 | OUTPUT: |
390 | OUTPUT: |
293 | RETVAL |
391 | RETVAL |
|
|
392 | |
|
|
393 | void |
|
|
394 | compare (GdkPixbuf *a, GdkPixbuf *b) |
|
|
395 | PPCODE: |
|
|
396 | { |
|
|
397 | int w = gdk_pixbuf_get_width (a); |
|
|
398 | int h = gdk_pixbuf_get_height (a); |
|
|
399 | int sa = gdk_pixbuf_get_rowstride (a); |
|
|
400 | int sb = gdk_pixbuf_get_rowstride (b); |
|
|
401 | |
|
|
402 | guchar *pa = gdk_pixbuf_get_pixels (a); |
|
|
403 | guchar *pb = gdk_pixbuf_get_pixels (b); |
|
|
404 | |
|
|
405 | int x, y; |
|
|
406 | |
|
|
407 | assert (w == gdk_pixbuf_get_width (b)); |
|
|
408 | assert (h == gdk_pixbuf_get_height (b)); |
|
|
409 | |
|
|
410 | assert (gdk_pixbuf_get_n_channels (a) == 3); |
|
|
411 | assert (gdk_pixbuf_get_n_channels (b) == 3); |
|
|
412 | |
|
|
413 | double diff = 0.; |
|
|
414 | int peak = 0; |
|
|
415 | |
|
|
416 | if (w && h) |
|
|
417 | for (y = 0; y < h; y++) |
|
|
418 | { |
|
|
419 | guchar *pa_ = pa + y * sa; |
|
|
420 | guchar *pb_ = pb + y * sb; |
|
|
421 | |
|
|
422 | for (x = 0; x < w; x++) |
|
|
423 | { |
|
|
424 | int d; |
|
|
425 | |
|
|
426 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
427 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
428 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
429 | } |
|
|
430 | } |
|
|
431 | |
|
|
432 | EXTEND (SP, 2); |
|
|
433 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
|
|
434 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
|
|
435 | } |
294 | |
436 | |
295 | ############################################################################# |
437 | ############################################################################# |
296 | |
438 | |
297 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
439 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
298 | |
440 | |
… | |
… | |
500 | ############################################################################# |
642 | ############################################################################# |
501 | |
643 | |
502 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
644 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
503 | |
645 | |
504 | SV * |
646 | SV * |
|
|
647 | extract_features (SV *ar) |
|
|
648 | CODE: |
|
|
649 | { |
|
|
650 | int i; |
|
|
651 | AV *av, *result; |
|
|
652 | |
|
|
653 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
|
|
654 | croak ("Not an array ref as first argument to make_histogram"); |
|
|
655 | |
|
|
656 | av = (AV *) SvRV (ar); |
|
|
657 | result = newAV (); |
|
|
658 | |
|
|
659 | for (i = 0; i <= av_len (av); ++i) |
|
|
660 | { |
|
|
661 | SV *sv = *av_fetch (av, i, 1); |
|
|
662 | SV *histsv = newSV (9 * sizeof (float) + 1); |
|
|
663 | |
|
|
664 | SvPOK_on (histsv); |
|
|
665 | SvCUR_set (histsv, 9 * sizeof (float)); |
|
|
666 | float *hist = (float *)SvPVX (histsv); |
|
|
667 | |
|
|
668 | struct feature f_h, f_s, f_v; |
|
|
669 | feature_init (&f_h); |
|
|
670 | feature_init (&f_s); |
|
|
671 | feature_init (&f_v); |
|
|
672 | |
|
|
673 | { |
|
|
674 | STRLEN len; |
|
|
675 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
676 | while (len >= 3) |
|
|
677 | { |
|
|
678 | unsigned int r, g, b, h, s, v; |
|
|
679 | r = *buf++; g = *buf++; b = *buf++; |
|
|
680 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
681 | |
|
|
682 | feature_update_pass_1 (&f_h, h); |
|
|
683 | feature_update_pass_1 (&f_s, s); |
|
|
684 | feature_update_pass_1 (&f_v, v); |
|
|
685 | |
|
|
686 | len -= 3; |
|
|
687 | } |
|
|
688 | |
|
|
689 | feature_finish_pass_1 (&f_h); |
|
|
690 | feature_finish_pass_1 (&f_s); |
|
|
691 | feature_finish_pass_1 (&f_v); |
|
|
692 | } |
|
|
693 | |
|
|
694 | { |
|
|
695 | STRLEN len; |
|
|
696 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
697 | while (len >= 3) |
|
|
698 | { |
|
|
699 | unsigned int r, g, b, h, s, v; |
|
|
700 | r = *buf++; g = *buf++; b = *buf++; |
|
|
701 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
702 | |
|
|
703 | feature_update_pass_2 (&f_h, h); |
|
|
704 | feature_update_pass_2 (&f_s, s); |
|
|
705 | feature_update_pass_2 (&f_v, v); |
|
|
706 | |
|
|
707 | len -= 3; |
|
|
708 | } |
|
|
709 | |
|
|
710 | feature_finish_pass_2 (&f_h); |
|
|
711 | feature_finish_pass_2 (&f_s); |
|
|
712 | feature_finish_pass_2 (&f_v); |
|
|
713 | } |
|
|
714 | |
|
|
715 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
716 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
717 | hist [6] = f_v.v1 ; hist [7] = f_v.v2 ; hist [8] = f_v.v3 ; |
|
|
718 | |
|
|
719 | av_push (result, histsv); |
|
|
720 | } |
|
|
721 | |
|
|
722 | RETVAL = newRV_noinc ((SV *)result); |
|
|
723 | } |
|
|
724 | OUTPUT: |
|
|
725 | RETVAL |
|
|
726 | |
|
|
727 | SV * |
505 | make_histograms (SV *ar) |
728 | make_histograms (SV *ar) |
506 | CODE: |
729 | CODE: |
507 | { |
730 | { |
508 | int i; |
731 | int i; |
509 | AV *av, *result; |
732 | AV *av, *result; |