ViewVC Help
View File | Revision Log | Show Annotations | Download File
/cvs/CV/CV.xs
(Generate patch)

Comparing CV/CV.xs (file contents):
Revision 1.14 by root, Tue Feb 8 12:51:12 2005 UTC vs.
Revision 1.49 by root, Wed May 24 19:22:32 2017 UTC

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)
13 27
14#define LINELENGTH 240 28#define LINELENGTH 240
15 29
16#define ELLIPSIS "\xe2\x80\xa6" 30#define ELLIPSIS "\xe2\x80\xa6"
17 31
18static guint32 a85_val; 32typedef char *octet_string;
19static guint a85_cnt;
20static guchar a85_buf[LINELENGTH], *a85_ptr;
21 33
22static void 34struct jpg_err_mgr
23a85_init (void)
24{ 35{
25 a85_cnt = 4; 36 struct jpeg_error_mgr err;
26 a85_ptr = a85_buf; 37 jmp_buf setjmp_buffer;
27} 38};
28 39
29static void 40static void
30a85_push (PerlIO *fp, guchar c) 41cv_error_exit (j_common_ptr cinfo)
31{ 42{
32 a85_val = a85_val << 8 | c; 43 cinfo->err->output_message (cinfo);
33 44 longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
34 if (!--a85_cnt)
35 {
36 a85_cnt = 4;
37 if (a85_val)
38 {
39 a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
40 a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
41 a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
42 a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
43 a85_ptr[0] = (a85_val ) + 33;
44
45 a85_ptr += 5;
46 }
47 else
48 *a85_ptr++ = 'z';
49
50 if (a85_ptr >= a85_buf + sizeof (a85_buf) - 7)
51 {
52 *a85_ptr++ = '\n';
53 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
54 a85_ptr = a85_buf;
55 }
56 }
57
58} 45}
59 46
60static void 47static void
61a85_finish (PerlIO *fp) 48cv_error_output (j_common_ptr cinfo)
62{ 49{
63 while (a85_cnt != 4) 50 char msg[JMSG_LENGTH_MAX];
64 a85_push (fp, 0);
65 51
66 *a85_ptr++ = '~'; // probably buggy end-marker 52 cinfo->err->format_message (cinfo, msg);
67 *a85_ptr++ = '>'; // probably buggy end-marker
68 *a85_ptr++ = '\n';
69 53
70 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); 54 fprintf (stderr, "JPEG decoding error: %s\n", msg);
55 return;
71} 56}
72 57
73static void 58static void
74rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, 59rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
75 unsigned int *h, unsigned int *s, unsigned int *v) 60 unsigned int *h, unsigned int *s, unsigned int *v)
95 80
96 *h &= 255; 81 *h &= 255;
97 } 82 }
98} 83}
99 84
85struct feature {
86 float v1, v2, v3; // mean, square, cube
87 int n;
88};
89
90static void
91feature_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.
101static void
102feature_update_pass_1 (struct feature *f, unsigned int v)
103{
104 f->v1 += v;
105 f->n += 1;
106}
107
108static void
109feature_finish_pass_1 (struct feature *f)
110{
111 if (f->n < 1)
112 return;
113
114 f->v1 /= f->n;
115}
116
117static void
118feature_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
126static void
127feature_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}
139
140static guint32 a85_val;
141static guint a85_cnt;
142static guchar a85_buf[LINELENGTH], *a85_ptr;
143
144static void
145a85_init (void)
146{
147 a85_cnt = 4;
148 a85_ptr = a85_buf;
149}
150
151static void
152a85_push (PerlIO *fp, guchar c)
153{
154 a85_val = a85_val << 8 | c;
155
156 if (!--a85_cnt)
157 {
158 a85_cnt = 4;
159 if (a85_val)
160 {
161 a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
162 a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
163 a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
164 a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
165 a85_ptr[0] = (a85_val ) + 33;
166
167 a85_ptr += 5;
168 }
169 else
170 *a85_ptr++ = 'z';
171
172 if (a85_ptr >= a85_buf + sizeof (a85_buf) - 7)
173 {
174 *a85_ptr++ = '\n';
175 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
176 a85_ptr = a85_buf;
177 }
178 }
179
180}
181
182static void
183a85_finish (PerlIO *fp)
184{
185 while (a85_cnt != 4)
186 a85_push (fp, 0);
187
188 *a85_ptr++ = '~'; // probably buggy end-marker
189 *a85_ptr++ = '>'; // probably buggy end-marker
190 *a85_ptr++ = '\n';
191
192 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
193}
194
100///////////////////////////////////////////////////////////////////////////// 195/////////////////////////////////////////////////////////////////////////////
101 196
102MODULE = Gtk2::CV PACKAGE = Gtk2::CV 197MODULE = Gtk2::CV PACKAGE = Gtk2::CV
103 198
104PROTOTYPES: ENABLE 199PROTOTYPES: ENABLE
105 200
201# missing function in perl. really :)
202int
203common_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
218const char *
219magic (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
239const char *
240magic_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
264void
265gdk_window_clear_hints (GdkWindow *window)
266 CODE:
267 gdk_window_set_geometry_hints (window, 0, 0);
268
269gboolean
270gdk_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
106GdkPixbuf_noinc * 280GdkPixbuf_noinc *
107transpose (GdkPixbuf *pb) 281dealpha_expose (GdkPixbuf *pb)
108 CODE: 282 CODE:
283 perlinterp_release ();
109{ 284{
110 int w = gdk_pixbuf_get_width (pb); 285 int w = gdk_pixbuf_get_width (pb);
111 int h = gdk_pixbuf_get_height (pb); 286 int h = gdk_pixbuf_get_height (pb);
112 int bpp = gdk_pixbuf_get_n_channels (pb); 287 int bpp = gdk_pixbuf_get_n_channels (pb);
113 int x, y, i; 288 int x, y, i;
114 guchar *src = gdk_pixbuf_get_pixels (pb), *dst; 289 guchar *src = gdk_pixbuf_get_pixels (pb), *dst;
115 int sstr = gdk_pixbuf_get_rowstride (pb), dstr; 290 int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
116 291
117 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w); 292 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
118 293
119 dst = gdk_pixbuf_get_pixels (RETVAL); 294 dst = gdk_pixbuf_get_pixels (RETVAL);
120 dstr = gdk_pixbuf_get_rowstride (RETVAL); 295 dstr = gdk_pixbuf_get_rowstride (RETVAL);
121 296
122 for (y = 0; y < h; y++)
123 for (x = 0; x < w; x++) 297 for (x = 0; x < w; x++)
298 for (y = 0; y < h; y++)
124 for (i = 0; i < bpp; i++) 299 for (i = 0; i < 3; i++)
125 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];
126} 301}
302 perlinterp_acquire ();
127 OUTPUT: 303 OUTPUT:
128 RETVAL 304 RETVAL
129 305
130GdkPixbuf_noinc * 306GdkPixbuf_noinc *
131flop (GdkPixbuf *pb) 307rotate (GdkPixbuf *pb, int angle)
132 CODE: 308 CODE:
133{ 309 perlinterp_release ();
134 int w = gdk_pixbuf_get_width (pb); 310 if (angle < 0)
135 int h = gdk_pixbuf_get_height (pb); 311 angle += 360;
136 int bpp = gdk_pixbuf_get_n_channels (pb); 312 RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
137 int x, y, i; 313 : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
138 guchar *src = gdk_pixbuf_get_pixels (pb), *dst; 314 : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
139 int sstr = gdk_pixbuf_get_rowstride (pb), dstr; 315 : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
140 316 : angle);
141 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h); 317 perlinterp_acquire ();
142
143 dst = gdk_pixbuf_get_pixels (RETVAL);
144 dstr = gdk_pixbuf_get_rowstride (RETVAL);
145
146 for (y = 0; y < h; y++)
147 for (x = 0; x < w; x++)
148 for (i = 0; i < bpp; i++)
149 dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
150}
151 OUTPUT: 318 OUTPUT:
152 RETVAL 319 RETVAL
153 320
321GdkPixbuf_noinc *
322load_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
456void
457compare (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
154############################################################################# 504#############################################################################
155 505
156MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer 506MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
157 507
508# currently only works for filenames (octet strings)
509
510SV *
511foldcase (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
158GdkPixbuf_noinc * 558GdkPixbuf_noinc *
159p7_to_pb (int w, int h, guchar *src) 559p7_to_pb (int w, int h, SV *src_sv)
560 PROTOTYPE: @
160 CODE: 561 CODE:
161{ 562{
162 int x, y; 563 int x, y;
163 guchar *dst, *d; 564 guchar *dst, *d;
164 int dstr; 565 int dstr;
566 guchar *src = (guchar *)SvPVbyte_nolen (src_sv);
165 567
166 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); 568 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
167 dst = gdk_pixbuf_get_pixels (RETVAL); 569 dst = gdk_pixbuf_get_pixels (RETVAL);
168 dstr = gdk_pixbuf_get_rowstride (RETVAL); 570 dstr = gdk_pixbuf_get_rowstride (RETVAL);
169 571
178 } 580 }
179} 581}
180 OUTPUT: 582 OUTPUT:
181 RETVAL 583 RETVAL
182 584
585#############################################################################
586
587MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
588
589void
590dump_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
611void
612dump_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
631MODULE = Gtk2::CV PACKAGE = Gtk2::CV
632
183SV * 633SV *
184pb_to_p7 (GdkPixbuf *pb) 634pb_to_hv84 (GdkPixbuf *pb)
185 CODE: 635 CODE:
186{ 636{
187 int w = gdk_pixbuf_get_width (pb); 637 int w = gdk_pixbuf_get_width (pb);
188 int h = gdk_pixbuf_get_height (pb); 638 int h = gdk_pixbuf_get_height (pb);
189 int x, y; 639 int x, y;
190 guchar *dst; 640 guchar *dst;
191 int bpp = gdk_pixbuf_get_n_channels (pb); 641 int bpp = gdk_pixbuf_get_n_channels (pb);
192 guchar *src = gdk_pixbuf_get_pixels (pb); 642 guchar *src = gdk_pixbuf_get_pixels (pb);
193 int sstr = gdk_pixbuf_get_rowstride (pb); 643 int sstr = gdk_pixbuf_get_rowstride (pb);
194 int Er[IW], Eg[IW], Eb[IW];
195 int seed = 77;
196
197 RETVAL = newSV (w * h);
198 SvPOK_only (RETVAL);
199 SvCUR_set (RETVAL, w * h);
200
201 dst = SvPVX (RETVAL);
202
203 memset (Er, 0, sizeof (int) * IW);
204 memset (Eg, 0, sizeof (int) * IW);
205 memset (Eb, 0, sizeof (int) * IW);
206
207 /* some primitive error distribution + random dithering */
208
209 for (y = 0; y < h; y++)
210 {
211 int er = 0, eg = 0, eb = 0;
212
213 for (x = 0; x < w; x++)
214 {
215 int r, g, b;
216 guchar *p = src + x * bpp + y * sstr;
217
218 r = ((p[0] + er + Er[x]) * 7 + 128) / 255;
219 g = ((p[1] + eg + Eg[x]) * 7 + 128) / 255;
220 b = ((p[2] + eb + Eb[x]) * 3 + 128) / 255;
221
222 r = r > 7 ? 7 : r < 0 ? 0 : r;
223 g = g > 7 ? 7 : g < 0 ? 0 : g;
224 b = b > 3 ? 3 : b < 0 ? 0 : b;
225
226 er += p[0] - (r * 255 + 4) / 7;
227 eg += p[1] - (g * 255 + 4) / 7;
228 eb += p[2] - (b * 255 + 2) / 3;
229
230 Er[x] = er / 2; er -= er / 2 + RAND % 7 - 3;
231 Eg[x] = eg / 2; eg -= eg / 2 + RAND % 7 - 3;
232 Eb[x] = eb / 2; eb -= eb / 2 + RAND % 7 - 3;
233
234 *dst++ = r << 5 | g << 2 | b;
235 }
236 }
237}
238 OUTPUT:
239 RETVAL
240
241SV *
242make_histogram (SV *ar)
243 CODE:
244{
245 int i;
246 AV *av, *result;
247
248 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
249 croak ("Not an array ref as first argument to make_histogram");
250
251 av = (AV *) SvRV (ar);
252 result = newAV ();
253
254 for (i = 0; i <= av_len (av); ++i)
255 {
256 int bigst, j;
257 SV *sv = *av_fetch (av, i, 1);
258 STRLEN len;
259 char *buf = SvPVbyte (sv, len);
260
261 int tmphist[256];
262 char *hist;
263
264 SV *histsv = newSV (257);
265 SvPOK_on (histsv);
266 SvCUR_set (histsv, 256);
267 hist = SvPVX (histsv);
268
269 Zero (tmphist, 256, int);
270
271 j = len;
272
273 while (len--)
274 ++tmphist[(unsigned char)*buf++ & 0xda];
275
276 len = j;
277
278 bigst = 1;
279 for (j = 0; j < 256; ++j)
280 if (tmphist[j] > bigst)
281 bigst = tmphist[j];
282
283 for (j = 0; j < 256; ++j)
284 hist[j] = tmphist[j] * 255 / (len == 0 ? 1 : len);//bigst;
285
286 av_push (result, histsv);
287 }
288
289 RETVAL = newRV_noinc ((SV *)result);
290}
291 OUTPUT:
292 RETVAL
293
294
295
296#############################################################################
297
298MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
299
300void
301dump_ascii85 (PerlIO *fp, GdkPixbuf *pb)
302 CODE:
303{
304 int w = gdk_pixbuf_get_width (pb);
305 int h = gdk_pixbuf_get_height (pb);
306 int x, y, i;
307 guchar *dst;
308 int bpp = gdk_pixbuf_get_n_channels (pb);
309 guchar *src = gdk_pixbuf_get_pixels (pb);
310 int sstr = gdk_pixbuf_get_rowstride (pb);
311
312 a85_init ();
313
314 for (y = 0; y < h; y++)
315 for (x = 0; x < w; x++)
316 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
317 a85_push (fp, src [x * bpp + y * sstr + i]);
318
319 a85_finish (fp);
320}
321
322void
323dump_binary (PerlIO *fp, GdkPixbuf *pb)
324 CODE:
325{
326 int w = gdk_pixbuf_get_width (pb);
327 int h = gdk_pixbuf_get_height (pb);
328 int x, y, i;
329 guchar *dst;
330 int bpp = gdk_pixbuf_get_n_channels (pb);
331 guchar *src = gdk_pixbuf_get_pixels (pb);
332 int sstr = gdk_pixbuf_get_rowstride (pb);
333
334 for (y = 0; y < h; y++)
335 for (x = 0; x < w; x++)
336 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
337 PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
338}
339
340#############################################################################
341
342MODULE = Gtk2::CV PACKAGE = Gtk2::CV
343
344SV *
345pb_to_hv84 (GdkPixbuf *pb)
346 CODE:
347{
348 int w = gdk_pixbuf_get_width (pb);
349 int h = gdk_pixbuf_get_height (pb);
350 int x, y;
351 guchar *dst;
352 int bpp = gdk_pixbuf_get_n_channels (pb);
353 guchar *src = gdk_pixbuf_get_pixels (pb);
354 int sstr = gdk_pixbuf_get_rowstride (pb);
355 644
356 RETVAL = newSV (6 * 8 * 12 / 8); 645 RETVAL = newSV (6 * 8 * 12 / 8);
357 SvPOK_only (RETVAL); 646 SvPOK_only (RETVAL);
358 SvCUR_set (RETVAL, 6 * 8 * 12 / 8); 647 SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
359 648
360 dst = SvPVX (RETVAL); 649 dst = (guchar *)SvPVX (RETVAL);
361 650
362 /* some primitive error distribution + random dithering */ 651 /* some primitive error distribution + random dithering */
363 652
364 for (y = 0; y < h; y++) 653 for (y = 0; y < h; y++)
365 { 654 {
423 } 712 }
424} 713}
425 OUTPUT: 714 OUTPUT:
426 RETVAL 715 RETVAL
427 716
717#############################################################################
428 718
719MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
429 720
721SV *
722extract_features (SV *ar)
723 CODE:
724{
725 int i;
726 AV *av, *result;
727
728 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
729 croak ("Not an array ref as first argument to extract_features");
730
731 av = (AV *) SvRV (ar);
732 result = newAV ();
733
734 for (i = 0; i <= av_len (av); ++i)
735 {
736 SV *sv = *av_fetch (av, i, 1);
737 SV *histsv = newSV (9 * sizeof (float) + 1);
738
739 SvPOK_on (histsv);
740 SvCUR_set (histsv, 9 * sizeof (float));
741 float *hist = (float *)SvPVX (histsv);
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);
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);
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

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines