… | |
… | |
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 | f->v1 /= 255.; |
|
|
124 | f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2); |
|
|
125 | f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.); |
|
|
126 | } |
|
|
127 | |
70 | static guint32 a85_val; |
128 | static guint32 a85_val; |
71 | static guint a85_cnt; |
129 | static guint a85_cnt; |
72 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
130 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
73 | |
131 | |
74 | static void |
132 | static void |
… | |
… | |
126 | |
184 | |
127 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
185 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
128 | |
186 | |
129 | PROTOTYPES: ENABLE |
187 | PROTOTYPES: ENABLE |
130 | |
188 | |
|
|
189 | # missing function in perl. really :) |
|
|
190 | int |
|
|
191 | common_prefix_length (a, b) |
|
|
192 | unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
|
|
193 | unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
|
|
194 | CODE: |
|
|
195 | RETVAL = 0; |
|
|
196 | |
|
|
197 | while (*a == *b && *a) |
|
|
198 | { |
|
|
199 | RETVAL += (*a & 0xc0) != 0x80; |
|
|
200 | a++, b++; |
|
|
201 | } |
|
|
202 | |
|
|
203 | OUTPUT: |
|
|
204 | RETVAL |
|
|
205 | |
|
|
206 | const char * |
|
|
207 | magic (const char *path) |
|
|
208 | CODE: |
|
|
209 | { |
|
|
210 | static magic_t cookie; |
|
|
211 | |
|
|
212 | if (!cookie) |
|
|
213 | { |
|
|
214 | cookie = magic_open (MAGIC_NONE); |
|
|
215 | |
|
|
216 | if (cookie) |
|
|
217 | magic_load (cookie, 0); |
|
|
218 | else |
|
|
219 | XSRETURN_UNDEF; |
|
|
220 | } |
|
|
221 | |
|
|
222 | RETVAL = magic_file (cookie, path); |
|
|
223 | } |
|
|
224 | OUTPUT: |
|
|
225 | RETVAL |
|
|
226 | |
|
|
227 | const char * |
|
|
228 | magic_mime (const char *path) |
|
|
229 | CODE: |
|
|
230 | { |
|
|
231 | static magic_t cookie; |
|
|
232 | |
|
|
233 | if (!cookie) |
|
|
234 | { |
|
|
235 | cookie = magic_open (MAGIC_MIME); |
|
|
236 | |
|
|
237 | if (cookie) |
|
|
238 | magic_load (cookie, 0); |
|
|
239 | else |
|
|
240 | XSRETURN_UNDEF; |
|
|
241 | } |
|
|
242 | |
|
|
243 | RETVAL = magic_file (cookie, path); |
|
|
244 | } |
|
|
245 | OUTPUT: |
|
|
246 | RETVAL |
|
|
247 | |
131 | # missing in Gtk2 perl module |
248 | # missing in Gtk2 perl module |
132 | |
249 | |
133 | gboolean |
250 | gboolean |
134 | gdk_net_wm_supports (GdkAtom property) |
251 | gdk_net_wm_supports (GdkAtom property) |
135 | CODE: |
252 | CODE: |
… | |
… | |
164 | } |
281 | } |
165 | OUTPUT: |
282 | OUTPUT: |
166 | RETVAL |
283 | RETVAL |
167 | |
284 | |
168 | GdkPixbuf_noinc * |
285 | GdkPixbuf_noinc * |
169 | transpose (GdkPixbuf *pb) |
286 | rotate (GdkPixbuf *pb, int angle) |
170 | CODE: |
287 | CODE: |
171 | { |
288 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
172 | int w = gdk_pixbuf_get_width (pb); |
289 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
173 | int h = gdk_pixbuf_get_height (pb); |
290 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
174 | int bpp = gdk_pixbuf_get_n_channels (pb); |
291 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
175 | int x, y, i; |
292 | : angle); |
176 | guchar *src = gdk_pixbuf_get_pixels (pb), *dst; |
|
|
177 | int sstr = gdk_pixbuf_get_rowstride (pb), dstr; |
|
|
178 | |
|
|
179 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w); |
|
|
180 | |
|
|
181 | dst = gdk_pixbuf_get_pixels (RETVAL); |
|
|
182 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
|
|
183 | |
|
|
184 | for (y = 0; y < h; y++) |
|
|
185 | for (x = 0; x < w; x++) |
|
|
186 | for (i = 0; i < bpp; i++) |
|
|
187 | dst[y * bpp + x * dstr + i] = src[x * bpp + y * sstr + i]; |
|
|
188 | } |
|
|
189 | OUTPUT: |
|
|
190 | RETVAL |
|
|
191 | |
|
|
192 | GdkPixbuf_noinc * |
|
|
193 | flop (GdkPixbuf *pb) |
|
|
194 | CODE: |
|
|
195 | { |
|
|
196 | int w = gdk_pixbuf_get_width (pb); |
|
|
197 | int h = gdk_pixbuf_get_height (pb); |
|
|
198 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
199 | int x, y, i; |
|
|
200 | guchar *src = gdk_pixbuf_get_pixels (pb), *dst; |
|
|
201 | int sstr = gdk_pixbuf_get_rowstride (pb), dstr; |
|
|
202 | |
|
|
203 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h); |
|
|
204 | |
|
|
205 | dst = gdk_pixbuf_get_pixels (RETVAL); |
|
|
206 | dstr = gdk_pixbuf_get_rowstride (RETVAL); |
|
|
207 | |
|
|
208 | for (y = 0; y < h; y++) |
|
|
209 | for (x = 0; x < w; x++) |
|
|
210 | for (i = 0; i < bpp; i++) |
|
|
211 | dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i]; |
|
|
212 | } |
|
|
213 | OUTPUT: |
293 | OUTPUT: |
214 | RETVAL |
294 | RETVAL |
215 | |
295 | |
216 | GdkPixbuf_noinc * |
296 | GdkPixbuf_noinc * |
217 | load_jpeg (SV *path, int thumbnail=0) |
297 | load_jpeg (SV *path, int thumbnail=0) |
… | |
… | |
221 | struct jpg_err_mgr jerr; |
301 | struct jpg_err_mgr jerr; |
222 | guchar *data; |
302 | guchar *data; |
223 | int rs; |
303 | int rs; |
224 | FILE *fp; |
304 | FILE *fp; |
225 | volatile GdkPixbuf *pb = 0; |
305 | volatile GdkPixbuf *pb = 0; |
226 | gchar *filename; |
|
|
227 | |
306 | |
228 | RETVAL = 0; |
307 | RETVAL = 0; |
229 | |
308 | |
230 | filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0); |
|
|
231 | fp = fopen (filename, "rb"); |
309 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
232 | g_free (filename); |
|
|
233 | |
310 | |
234 | if (!fp) |
311 | if (!fp) |
235 | XSRETURN_UNDEF; |
312 | XSRETURN_UNDEF; |
236 | |
313 | |
237 | cinfo.err = jpeg_std_error (&jerr.err); |
314 | cinfo.err = jpeg_std_error (&jerr.err); |
… | |
… | |
309 | fclose (fp); |
386 | fclose (fp); |
310 | jpeg_destroy_decompress (&cinfo); |
387 | jpeg_destroy_decompress (&cinfo); |
311 | } |
388 | } |
312 | OUTPUT: |
389 | OUTPUT: |
313 | RETVAL |
390 | RETVAL |
|
|
391 | |
|
|
392 | void |
|
|
393 | compare (GdkPixbuf *a, GdkPixbuf *b) |
|
|
394 | PPCODE: |
|
|
395 | { |
|
|
396 | int w = gdk_pixbuf_get_width (a); |
|
|
397 | int h = gdk_pixbuf_get_height (a); |
|
|
398 | int sa = gdk_pixbuf_get_rowstride (a); |
|
|
399 | int sb = gdk_pixbuf_get_rowstride (b); |
|
|
400 | |
|
|
401 | guchar *pa = gdk_pixbuf_get_pixels (a); |
|
|
402 | guchar *pb = gdk_pixbuf_get_pixels (b); |
|
|
403 | |
|
|
404 | int x, y; |
|
|
405 | |
|
|
406 | assert (w == gdk_pixbuf_get_width (b)); |
|
|
407 | assert (h == gdk_pixbuf_get_height (b)); |
|
|
408 | |
|
|
409 | assert (gdk_pixbuf_get_n_channels (a) == 3); |
|
|
410 | assert (gdk_pixbuf_get_n_channels (b) == 3); |
|
|
411 | |
|
|
412 | double diff = 0.; |
|
|
413 | int peak = 0; |
|
|
414 | |
|
|
415 | if (w && h) |
|
|
416 | for (y = 0; y < h; y++) |
|
|
417 | { |
|
|
418 | guchar *pa_ = pa + y * sa; |
|
|
419 | guchar *pb_ = pb + y * sb; |
|
|
420 | |
|
|
421 | for (x = 0; x < w; x++) |
|
|
422 | { |
|
|
423 | int d; |
|
|
424 | |
|
|
425 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
|
|
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 | } |
|
|
429 | } |
|
|
430 | |
|
|
431 | EXTEND (SP, 2); |
|
|
432 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
|
|
433 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
|
|
434 | } |
314 | |
435 | |
315 | ############################################################################# |
436 | ############################################################################# |
316 | |
437 | |
317 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
438 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
318 | |
439 | |
… | |
… | |
383 | } |
504 | } |
384 | } |
505 | } |
385 | OUTPUT: |
506 | OUTPUT: |
386 | RETVAL |
507 | RETVAL |
387 | |
508 | |
|
|
509 | ############################################################################# |
|
|
510 | |
|
|
511 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript |
|
|
512 | |
|
|
513 | void |
|
|
514 | dump_ascii85 (PerlIO *fp, GdkPixbuf *pb) |
|
|
515 | CODE: |
|
|
516 | { |
|
|
517 | int w = gdk_pixbuf_get_width (pb); |
|
|
518 | int h = gdk_pixbuf_get_height (pb); |
|
|
519 | int x, y, i; |
|
|
520 | guchar *dst; |
|
|
521 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
522 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
523 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
524 | |
|
|
525 | a85_init (); |
|
|
526 | |
|
|
527 | for (y = 0; y < h; y++) |
|
|
528 | for (x = 0; x < w; x++) |
|
|
529 | for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
|
|
530 | a85_push (fp, src [x * bpp + y * sstr + i]); |
|
|
531 | |
|
|
532 | a85_finish (fp); |
|
|
533 | } |
|
|
534 | |
|
|
535 | void |
|
|
536 | dump_binary (PerlIO *fp, GdkPixbuf *pb) |
|
|
537 | CODE: |
|
|
538 | { |
|
|
539 | int w = gdk_pixbuf_get_width (pb); |
|
|
540 | int h = gdk_pixbuf_get_height (pb); |
|
|
541 | int x, y, i; |
|
|
542 | guchar *dst; |
|
|
543 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
544 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
545 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
546 | |
|
|
547 | for (y = 0; y < h; y++) |
|
|
548 | for (x = 0; x < w; x++) |
|
|
549 | for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
|
|
550 | PerlIO_putc (fp, src [x * bpp + y * sstr + i]); |
|
|
551 | } |
|
|
552 | |
|
|
553 | ############################################################################# |
|
|
554 | |
|
|
555 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
|
|
556 | |
388 | SV * |
557 | SV * |
389 | pb_to_p7 (GdkPixbuf *pb) |
558 | pb_to_hv84 (GdkPixbuf *pb) |
390 | CODE: |
559 | CODE: |
391 | { |
560 | { |
392 | int w = gdk_pixbuf_get_width (pb); |
561 | int w = gdk_pixbuf_get_width (pb); |
393 | int h = gdk_pixbuf_get_height (pb); |
562 | int h = gdk_pixbuf_get_height (pb); |
394 | int x, y; |
563 | int x, y; |
395 | guchar *dst; |
564 | guchar *dst; |
396 | int bpp = gdk_pixbuf_get_n_channels (pb); |
565 | int bpp = gdk_pixbuf_get_n_channels (pb); |
397 | guchar *src = gdk_pixbuf_get_pixels (pb); |
566 | guchar *src = gdk_pixbuf_get_pixels (pb); |
398 | int sstr = gdk_pixbuf_get_rowstride (pb); |
567 | int sstr = gdk_pixbuf_get_rowstride (pb); |
399 | int Er[IW], Eg[IW], Eb[IW]; |
|
|
400 | int seed = 77; |
|
|
401 | |
|
|
402 | RETVAL = newSV (w * h); |
|
|
403 | SvPOK_only (RETVAL); |
|
|
404 | SvCUR_set (RETVAL, w * h); |
|
|
405 | |
|
|
406 | dst = (guchar *)SvPVX (RETVAL); |
|
|
407 | |
|
|
408 | memset (Er, 0, sizeof (int) * IW); |
|
|
409 | memset (Eg, 0, sizeof (int) * IW); |
|
|
410 | memset (Eb, 0, sizeof (int) * IW); |
|
|
411 | |
|
|
412 | /* some primitive error distribution + random dithering */ |
|
|
413 | |
|
|
414 | for (y = 0; y < h; y++) |
|
|
415 | { |
|
|
416 | int er = 0, eg = 0, eb = 0; |
|
|
417 | |
|
|
418 | for (x = 0; x < w; x++) |
|
|
419 | { |
|
|
420 | int r, g, b; |
|
|
421 | guchar *p = src + x * bpp + y * sstr; |
|
|
422 | |
|
|
423 | r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255; |
|
|
424 | g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255; |
|
|
425 | b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255; |
|
|
426 | |
|
|
427 | r = r > 7 ? 7 : r < 0 ? 0 : r; |
|
|
428 | g = g > 7 ? 7 : g < 0 ? 0 : g; |
|
|
429 | b = b > 3 ? 3 : b < 0 ? 0 : b; |
|
|
430 | |
|
|
431 | er += p[0] - (r * 255 + 4) / 7; |
|
|
432 | eg += p[1] - (g * 255 + 4) / 7; |
|
|
433 | eb += p[2] - (b * 255 + 2) / 3; |
|
|
434 | |
|
|
435 | Er[x] = er >> 1; er -= (er + 1) >> 1; |
|
|
436 | Eg[x] = eg >> 1; eg -= (eg + 1) >> 1; |
|
|
437 | Eb[x] = eb >> 1; eb -= (eb + 1) >> 1; |
|
|
438 | |
|
|
439 | *dst++ = r << 5 | g << 2 | b; |
|
|
440 | } |
|
|
441 | } |
|
|
442 | } |
|
|
443 | OUTPUT: |
|
|
444 | RETVAL |
|
|
445 | |
|
|
446 | ############################################################################# |
|
|
447 | |
|
|
448 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript |
|
|
449 | |
|
|
450 | void |
|
|
451 | dump_ascii85 (PerlIO *fp, GdkPixbuf *pb) |
|
|
452 | CODE: |
|
|
453 | { |
|
|
454 | int w = gdk_pixbuf_get_width (pb); |
|
|
455 | int h = gdk_pixbuf_get_height (pb); |
|
|
456 | int x, y, i; |
|
|
457 | guchar *dst; |
|
|
458 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
459 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
460 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
461 | |
|
|
462 | a85_init (); |
|
|
463 | |
|
|
464 | for (y = 0; y < h; y++) |
|
|
465 | for (x = 0; x < w; x++) |
|
|
466 | for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
|
|
467 | a85_push (fp, src [x * bpp + y * sstr + i]); |
|
|
468 | |
|
|
469 | a85_finish (fp); |
|
|
470 | } |
|
|
471 | |
|
|
472 | void |
|
|
473 | dump_binary (PerlIO *fp, GdkPixbuf *pb) |
|
|
474 | CODE: |
|
|
475 | { |
|
|
476 | int w = gdk_pixbuf_get_width (pb); |
|
|
477 | int h = gdk_pixbuf_get_height (pb); |
|
|
478 | int x, y, i; |
|
|
479 | guchar *dst; |
|
|
480 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
481 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
482 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
483 | |
|
|
484 | for (y = 0; y < h; y++) |
|
|
485 | for (x = 0; x < w; x++) |
|
|
486 | for (i = 0; i < (bpp < 3 ? 1 : 3); i++) |
|
|
487 | PerlIO_putc (fp, src [x * bpp + y * sstr + i]); |
|
|
488 | } |
|
|
489 | |
|
|
490 | ############################################################################# |
|
|
491 | |
|
|
492 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
|
|
493 | |
|
|
494 | SV * |
|
|
495 | pb_to_hv84 (GdkPixbuf *pb) |
|
|
496 | CODE: |
|
|
497 | { |
|
|
498 | int w = gdk_pixbuf_get_width (pb); |
|
|
499 | int h = gdk_pixbuf_get_height (pb); |
|
|
500 | int x, y; |
|
|
501 | guchar *dst; |
|
|
502 | int bpp = gdk_pixbuf_get_n_channels (pb); |
|
|
503 | guchar *src = gdk_pixbuf_get_pixels (pb); |
|
|
504 | int sstr = gdk_pixbuf_get_rowstride (pb); |
|
|
505 | |
568 | |
506 | RETVAL = newSV (6 * 8 * 12 / 8); |
569 | RETVAL = newSV (6 * 8 * 12 / 8); |
507 | SvPOK_only (RETVAL); |
570 | SvPOK_only (RETVAL); |
508 | SvCUR_set (RETVAL, 6 * 8 * 12 / 8); |
571 | SvCUR_set (RETVAL, 6 * 8 * 12 / 8); |
509 | |
572 | |
… | |
… | |
578 | ############################################################################# |
641 | ############################################################################# |
579 | |
642 | |
580 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
643 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
581 | |
644 | |
582 | SV * |
645 | SV * |
583 | make_histograms (SV *ar) |
646 | extract_features (SV *ar) |
584 | CODE: |
647 | CODE: |
585 | { |
648 | { |
586 | int i; |
649 | int i; |
587 | AV *av, *result; |
650 | AV *av, *result; |
588 | |
651 | |
589 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
652 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
590 | croak ("Not an array ref as first argument to make_histogram"); |
653 | croak ("Not an array ref as first argument to extract_features"); |
591 | |
654 | |
592 | av = (AV *) SvRV (ar); |
655 | av = (AV *) SvRV (ar); |
593 | result = newAV (); |
656 | result = newAV (); |
594 | |
657 | |
595 | for (i = 0; i <= av_len (av); ++i) |
658 | for (i = 0; i <= av_len (av); ++i) |
596 | { |
659 | { |
597 | const int HISTSIZE = 64; |
|
|
598 | |
|
|
599 | int j; |
|
|
600 | SV *sv = *av_fetch (av, i, 1); |
660 | SV *sv = *av_fetch (av, i, 1); |
601 | STRLEN len; |
|
|
602 | char *buf = SvPVbyte (sv, len); |
|
|
603 | |
|
|
604 | int tmphist[HISTSIZE]; |
|
|
605 | float *hist; |
|
|
606 | |
|
|
607 | SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); |
661 | SV *histsv = newSV (9 * sizeof (float) + 1); |
|
|
662 | |
608 | SvPOK_on (histsv); |
663 | SvPOK_on (histsv); |
609 | SvCUR_set (histsv, HISTSIZE * sizeof (float)); |
664 | SvCUR_set (histsv, 9 * sizeof (float)); |
610 | hist = (float *)SvPVX (histsv); |
665 | float *hist = (float *)SvPVX (histsv); |
611 | |
666 | |
612 | Zero (tmphist, sizeof (tmphist), char); |
667 | struct feature f_h, f_s, f_v; |
|
|
668 | feature_init (&f_h); |
|
|
669 | feature_init (&f_s); |
|
|
670 | feature_init (&f_v); |
613 | |
671 | |
614 | for (j = len; j--; ) |
|
|
615 | { |
|
|
616 | unsigned int idx |
|
|
617 | = ((*buf & 0xc0) >> 2) |
|
|
618 | | ((*buf & 0x18) >> 1) |
|
|
619 | | (*buf & 0x03); |
|
|
620 | |
|
|
621 | ++tmphist[idx]; |
|
|
622 | ++buf; |
|
|
623 | } |
|
|
624 | |
672 | { |
625 | for (j = 0; j < HISTSIZE; ++j) |
673 | STRLEN len; |
626 | hist[j] = (float)tmphist[j] / (len + 1e-30); |
674 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
675 | while (len >= 3) |
|
|
676 | { |
|
|
677 | unsigned int r, g, b, h, s, v; |
|
|
678 | r = *buf++; g = *buf++; b = *buf++; |
|
|
679 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
680 | |
|
|
681 | feature_update_pass_1 (&f_h, h); |
|
|
682 | feature_update_pass_1 (&f_s, s); |
|
|
683 | feature_update_pass_1 (&f_v, v); |
|
|
684 | |
|
|
685 | len -= 3; |
|
|
686 | } |
|
|
687 | |
|
|
688 | feature_finish_pass_1 (&f_h); |
|
|
689 | feature_finish_pass_1 (&f_s); |
|
|
690 | feature_finish_pass_1 (&f_v); |
|
|
691 | } |
|
|
692 | |
|
|
693 | { |
|
|
694 | STRLEN len; |
|
|
695 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
696 | while (len >= 3) |
|
|
697 | { |
|
|
698 | unsigned int r, g, b, h, s, v; |
|
|
699 | r = *buf++; g = *buf++; b = *buf++; |
|
|
700 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
701 | |
|
|
702 | feature_update_pass_2 (&f_h, h); |
|
|
703 | feature_update_pass_2 (&f_s, s); |
|
|
704 | feature_update_pass_2 (&f_v, v); |
|
|
705 | |
|
|
706 | len -= 3; |
|
|
707 | } |
|
|
708 | |
|
|
709 | feature_finish_pass_2 (&f_h); |
|
|
710 | feature_finish_pass_2 (&f_s); |
|
|
711 | feature_finish_pass_2 (&f_v); |
|
|
712 | } |
|
|
713 | |
|
|
714 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
715 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
716 | hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
627 | |
717 | |
628 | av_push (result, histsv); |
718 | av_push (result, histsv); |
629 | } |
719 | } |
630 | |
720 | |
631 | RETVAL = newRV_noinc ((SV *)result); |
721 | RETVAL = newRV_noinc ((SV *)result); |
632 | } |
722 | } |
633 | OUTPUT: |
723 | OUTPUT: |
634 | RETVAL |
724 | RETVAL |
635 | |
725 | |
636 | |
|
|