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

Comparing CV/CV.xs (file contents):
Revision 1.18 by root, Wed Jul 20 09:06:37 2005 UTC vs.
Revision 1.59 by root, Sat Apr 10 03:32:09 2021 UTC

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>
12#include <jerror.h>
13
14#include <glib.h>
15#include <gtk/gtk.h>
16#include <gdk/gdkx.h>
9#include <gdk-pixbuf/gdk-pixbuf.h> 17#include <gdk-pixbuf/gdk-pixbuf.h>
10 18
11#include <gperl.h> 19#include <gperl.h>
12#include <gtk2perl.h> 20#include <gtk2perl.h>
13 21
22#include <assert.h>
23
24#if WEBP
25#include <webp/demux.h>
26#include <webp/decode.h>
27#endif
28
29#if JXL
30#include <jxl/decode.h>
31#include "jxl/thread_parallel_runner.h"
32#endif
33
34#include "perlmulticore.h"
35
14#define IW 80 /* MUST match Schnauer.pm! */ 36#define IW 80 /* MUST match Schnauzer.pm! */
15#define IH 60 /* MUST match Schnauer.pm! */ 37#define IH 60 /* MUST match Schnauzer.pm! */
16
17#define RAND (seed = (seed + 7141) * 54773 % 134456)
18 38
19#define LINELENGTH 240 39#define LINELENGTH 240
20 40
21#define ELLIPSIS "\xe2\x80\xa6" 41#define ELLIPSIS "\xe2\x80\xa6"
42
43typedef char *octet_string;
44
45static magic_t magic_cookie[2]; /* !mime, mime */
22 46
23struct jpg_err_mgr 47struct jpg_err_mgr
24{ 48{
25 struct jpeg_error_mgr err; 49 struct jpeg_error_mgr err;
26 jmp_buf setjmp_buffer; 50 jmp_buf setjmp_buffer;
27}; 51};
28 52
29static void 53static void
30cv_error_exit (j_common_ptr cinfo) 54cv_error_exit (j_common_ptr cinfo)
31{ 55{
56 cinfo->err->output_message (cinfo);
32 longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99); 57 longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
33} 58}
34 59
35static void 60static void
36cv_error_output (j_common_ptr cinfo) 61cv_error_output (j_common_ptr cinfo)
37{ 62{
63 char msg[JMSG_LENGTH_MAX];
64
65 cinfo->err->format_message (cinfo, msg);
66
67 fprintf (stderr, "JPEG decoding error: %s\n", msg);
38 return; 68 return;
39} 69}
40 70
41static void 71static void
42rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, 72rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
63 93
64 *h &= 255; 94 *h &= 255;
65 } 95 }
66} 96}
67 97
98struct feature {
99 float v1, v2, v3; // mean, square, cube
100 int n;
101};
102
103static void
104feature_init (struct feature *f)
105{
106 f->v1 = 0.;
107 f->v2 = 0.;
108 f->v3 = 0.;
109 f->n = 0;
110}
111
112// didn't find an algorithm to neatly do mean, variance and skew in one pass.
113// elmex ist schuld.
114static void
115feature_update_pass_1 (struct feature *f, unsigned int v)
116{
117 f->v1 += v;
118 f->n += 1;
119}
120
121static void
122feature_finish_pass_1 (struct feature *f)
123{
124 if (f->n < 1)
125 return;
126
127 f->v1 /= f->n;
128}
129
130static void
131feature_update_pass_2 (struct feature *f, unsigned int v)
132{
133 float d = v - f->v1;
134
135 f->v2 += d * d;
136 f->v3 += d * d * d;
137}
138
139static void
140feature_finish_pass_2 (struct feature *f)
141{
142 if (f->n < 1)
143 return;
144
145 f->v2 /= f->n;
146 f->v3 /= f->n;
147
148 f->v1 /= 255.;
149 f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2);
150 f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.);
151}
152
68static guint32 a85_val; 153static guint32 a85_val;
69static guint a85_cnt; 154static guint a85_cnt;
70static guchar a85_buf[LINELENGTH], *a85_ptr; 155static guchar a85_buf[LINELENGTH], *a85_ptr;
71 156
72static void 157static void
84 if (!--a85_cnt) 169 if (!--a85_cnt)
85 { 170 {
86 a85_cnt = 4; 171 a85_cnt = 4;
87 if (a85_val) 172 if (a85_val)
88 { 173 {
89 a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; 174 a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
90 a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; 175 a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
91 a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; 176 a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
92 a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; 177 a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
93 a85_ptr[0] = (a85_val ) + 33; 178 a85_ptr[0] = (a85_val ) + 33;
94 179
95 a85_ptr += 5; 180 a85_ptr += 5;
96 } 181 }
119 204
120 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); 205 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
121} 206}
122 207
123///////////////////////////////////////////////////////////////////////////// 208/////////////////////////////////////////////////////////////////////////////
209// memory source for libjpeg
210
211static void cv_ms_init (j_decompress_ptr cinfo)
212{
213}
214
215static void cv_ms_term (j_decompress_ptr cinfo)
216{
217}
218
219static boolean cv_ms_fill (j_decompress_ptr cinfo)
220{
221 // unexpected EOF, warn and generate fake EOI marker
222
223 WARNMS (cinfo, JWRN_JPEG_EOF);
224
225 struct jpeg_source_mgr *src = (struct jpeg_source_mgr *)cinfo->src;
226
227 static const JOCTET eoi[] = { 0xFF, JPEG_EOI };
228
229 src->next_input_byte = eoi;
230 src->bytes_in_buffer = sizeof (eoi);
231
232 return TRUE;
233}
234
235static void cv_ms_skip (j_decompress_ptr cinfo, long num_bytes)
236{
237 struct jpeg_source_mgr *src = (struct jpeg_source_mgr *)cinfo->src;
238
239 src->next_input_byte += num_bytes;
240 src->bytes_in_buffer -= num_bytes;
241}
242
243static void cv_jpeg_mem_src (j_decompress_ptr cinfo, void *buf, size_t buflen)
244{
245 struct jpeg_source_mgr *src;
246
247 if (!cinfo->src)
248 cinfo->src = (struct jpeg_source_mgr *)
249 (*cinfo->mem->alloc_small) (
250 (j_common_ptr) cinfo, JPOOL_PERMANENT, sizeof (struct jpeg_source_mgr)
251 );
252
253 src = (struct jpeg_source_mgr *)cinfo->src;
254 src->init_source = cv_ms_init;
255 src->fill_input_buffer = cv_ms_fill;
256 src->skip_input_data = cv_ms_skip;
257 src->resync_to_restart = jpeg_resync_to_restart;
258 src->term_source = cv_ms_term;
259 src->next_input_byte = (JOCTET *)buf;
260 src->bytes_in_buffer = buflen;
261}
262
263/////////////////////////////////////////////////////////////////////////////
264
265/* great, the jpeg-xl reference implementaton requires us to parse bmff files */
266
267struct bmff_box
268{
269 char type[4];
270 const uint8_t *ptr;
271 size_t size;
272};
273
274static int
275bmff_parse_box (struct bmff_box *box, const uint8_t **next_in, size_t *avail_in)
276{
277 if (*avail_in < 8)
278 return 0;
279
280 box->size = ((*next_in)[0] << 24)
281 | ((*next_in)[1] << 16)
282 | ((*next_in)[2] << 8)
283 | ((*next_in)[3] );
284
285 if (box->size < 8)
286 return 0;
287
288 if (*avail_in < box->size)
289 return 0;
290
291 memcpy (box->type, *next_in + 4, 4);
292 box->ptr = *next_in + 8;
293
294 *next_in += box->size;
295 *avail_in -= box->size;
296
297 box->size -= 8;
298
299 return 1;
300}
301
302/////////////////////////////////////////////////////////////////////////////
124 303
125MODULE = Gtk2::CV PACKAGE = Gtk2::CV 304MODULE = Gtk2::CV PACKAGE = Gtk2::CV
126 305
127PROTOTYPES: ENABLE 306PROTOTYPES: ENABLE
128 307
308# calculate the common prefix length of two strings
309# missing function in perl. really :)
310int
311common_prefix_length (a, b)
312 unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg);
313 unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg);
314 CODE:
315 RETVAL = 0;
316
317 while (*a == *b && *a)
318 {
319 RETVAL += (*a & 0xc0) != 0x80;
320 a++, b++;
321 }
322
323 OUTPUT:
324 RETVAL
325
326const char *
327magic (SV *path_or_data)
328 ALIAS:
329 magic = 0
330 magic_mime = 1
331 magic_buffer = 2
332 magic_buffer_mime = 3
333 CODE:
334{
335 STRLEN len;
336 char *data = SvPVbyte (path_or_data, len);
337
338 if (!magic_cookie[0])
339 {
340 magic_cookie[0] = magic_open (MAGIC_SYMLINK);
341 magic_cookie[1] = magic_open (MAGIC_SYMLINK | MAGIC_MIME_TYPE);
342 magic_load (magic_cookie[0], 0);
343 magic_load (magic_cookie[1], 0);
344 }
345
346 perlinterp_release ();
347
348 RETVAL = ix & 2
349 ? magic_buffer (magic_cookie[ix & 1], data, len)
350 : magic_file (magic_cookie[ix & 1], data);
351
352 perlinterp_acquire ();
353}
354 OUTPUT:
355 RETVAL
356
357# missing/broken in Gtk2 perl module
358
359void
360gdk_window_clear_hints (GdkWindow *window)
361 CODE:
362 gdk_window_set_geometry_hints (window, 0, 0);
363
364gboolean
365gdk_net_wm_supports (GdkAtom property)
366 CODE:
367#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
368 RETVAL = gdk_net_wm_supports (property);
369#else
370 RETVAL = 0;
371#endif
372 OUTPUT:
373 RETVAL
374
129GdkPixbuf_noinc * 375GdkPixbuf_noinc *
130transpose (GdkPixbuf *pb) 376dealpha_expose (GdkPixbuf *pb)
131 CODE: 377 CODE:
378 perlinterp_release ();
132{ 379{
133 int w = gdk_pixbuf_get_width (pb); 380 int w = gdk_pixbuf_get_width (pb);
134 int h = gdk_pixbuf_get_height (pb); 381 int h = gdk_pixbuf_get_height (pb);
135 int bpp = gdk_pixbuf_get_n_channels (pb); 382 int bpp = gdk_pixbuf_get_n_channels (pb);
136 int x, y, i; 383 int x, y, i;
137 guchar *src = gdk_pixbuf_get_pixels (pb), *dst; 384 guchar *src = gdk_pixbuf_get_pixels (pb), *dst;
138 int sstr = gdk_pixbuf_get_rowstride (pb), dstr; 385 int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
139 386
140 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w); 387 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
141 388
142 dst = gdk_pixbuf_get_pixels (RETVAL); 389 dst = gdk_pixbuf_get_pixels (RETVAL);
143 dstr = gdk_pixbuf_get_rowstride (RETVAL); 390 dstr = gdk_pixbuf_get_rowstride (RETVAL);
144 391
145 for (y = 0; y < h; y++)
146 for (x = 0; x < w; x++) 392 for (x = 0; x < w; x++)
393 for (y = 0; y < h; y++)
147 for (i = 0; i < bpp; i++) 394 for (i = 0; i < 3; i++)
148 dst[y * bpp + x * dstr + i] = src[x * bpp + y * sstr + i]; 395 dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i];
149} 396}
397 perlinterp_acquire ();
150 OUTPUT: 398 OUTPUT:
151 RETVAL 399 RETVAL
152 400
153GdkPixbuf_noinc * 401GdkPixbuf_noinc *
154flop (GdkPixbuf *pb) 402rotate (GdkPixbuf *pb, int angle)
155 CODE: 403 CODE:
156{ 404 perlinterp_release ();
157 int w = gdk_pixbuf_get_width (pb); 405 if (angle < 0)
158 int h = gdk_pixbuf_get_height (pb); 406 angle += 360;
159 int bpp = gdk_pixbuf_get_n_channels (pb); 407 RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
160 int x, y, i; 408 : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
161 guchar *src = gdk_pixbuf_get_pixels (pb), *dst; 409 : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
162 int sstr = gdk_pixbuf_get_rowstride (pb), dstr; 410 : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
163 411 : angle);
164 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h); 412 perlinterp_acquire ();
165
166 dst = gdk_pixbuf_get_pixels (RETVAL);
167 dstr = gdk_pixbuf_get_rowstride (RETVAL);
168
169 for (y = 0; y < h; y++)
170 for (x = 0; x < w; x++)
171 for (i = 0; i < bpp; i++)
172 dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
173}
174 OUTPUT: 413 OUTPUT:
175 RETVAL 414 RETVAL
176 415
416const char *
417filetype (SV *image_data)
418 CODE:
419{
420 STRLEN data_len;
421 U8 *data = SvPVbyte (image_data, data_len);
422 static const unsigned char jxl_header[] = {
423 0, 0, 0, 0x0c, 0x4a, 0x58, 0x4c, 0x20,
424 0x0d, 0xa, 0x87, 0x0a, 0, 0, 0, 0x14,
425 0x66, 0x74, 0x79, 0x70, 0x6a, 0x78, 0x6c, 0x20,
426 0, 0, 0, 0, 0x6a, 0x78, 0x6c, 0x20,
427 };
428
429 if (data_len >= 20
430 && data[0] == 0xff
431 && data[1] == 0xd8
432 && data[2] == 0xff)
433 RETVAL = "image/jpeg";
434 else if (data_len >= 12
435 && data[ 0] == (U8)'R'
436 && data[ 1] == (U8)'I'
437 && data[ 2] == (U8)'F'
438 && data[ 3] == (U8)'F'
439 && data[ 8] == (U8)'W'
440 && data[ 9] == (U8)'E'
441 && data[10] == (U8)'B'
442 && data[11] == (U8)'P')
443 RETVAL = "image/webp";
444 else if (data_len >= 16
445 && data[ 0] == 0x89
446 && data[ 1] == (U8)'P'
447 && data[ 2] == (U8)'N'
448 && data[ 3] == (U8)'G'
449 && data[ 4] == 0x0d
450 && data[ 5] == 0x0a
451 && data[ 6] == 0x1a
452 && data[ 7] == 0x0a)
453 RETVAL = "image/png";
454 else if (data_len >= sizeof (jxl_header) && memcmp (data, jxl_header, sizeof (jxl_header)) == 0)
455 RETVAL = "image/jxl";
456 else
457 XSRETURN_UNDEF;
458}
459 OUTPUT:
460 RETVAL
461
177GdkPixbuf_noinc * 462GdkPixbuf_noinc *
178load_jpeg (char *path, int thumbnail=0) 463decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
464 CODE:
465{
466#if WEBP
467 STRLEN data_size;
468 int alpha;
469 WebPData data;
470 WebPDemuxer *demux;
471 WebPIterator iter;
472 WebPDecoderConfig config;
473 int inw, inh;
474
475 data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
476 data.size = data_size;
477
478 perlinterp_release ();
479
480 RETVAL = 0;
481
482 if (!(demux = WebPDemux (&data)))
483 goto err_demux;
484
485 if (!WebPDemuxGetFrame (demux, 1, &iter))
486 goto err_iter;
487
488 if (!WebPInitDecoderConfig (&config))
489 goto err_iter;
490
491 config.options.use_threads = 1;
492
493 if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
494 goto err_iter;
495
496 inw = config.input.width;
497 inh = config.input.height;
498
499 if (thumbnail)
500 {
501 if (inw * ih > inh * iw)
502 ih = (iw * inh + inw - 1) / inw;
503 else
504 iw = (ih * inw + inh - 1) / inh;
505
506 config.options.bypass_filtering = 1;
507 config.options.no_fancy_upsampling = 1;
508
509 config.options.use_scaling = 1;
510 config.options.scaled_width = iw;
511 config.options.scaled_height = ih;
512 }
513 else
514 {
515 iw = inw;
516 ih = inh;
517 }
518
519 alpha = !!config.input.has_alpha;
520
521 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
522 if (!RETVAL)
523 goto err_iter;
524
525 config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
526 config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
527 config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
528 config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
529 config.output.is_external_memory = 1;
530
531 if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
532 {
533 g_object_unref (RETVAL);
534 RETVAL = 0;
535 goto err_iter;
536 }
537
538 err_iter:
539 WebPDemuxReleaseIterator (&iter);
540 err_demux:
541 WebPDemuxDelete (demux);
542
543 perlinterp_acquire ();
544#else
545 croak ("load_webp: webp not enabled at compile time");
546#endif
547}
548 OUTPUT:
549 RETVAL
550
551GdkPixbuf_noinc *
552decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
553 CODE:
554{
555#if JXL
556 JxlDecoder *dec = JxlDecoderCreate (0);
557 JxlBasicInfo info;
558 const uint8_t *next_in = (uint8_t *)SvPVbyte_nolen (image_data);
559 size_t avail_in = SvCUR (image_data);
560 const char *error = 0;
561 void *runner = 0;
562 struct bmff_box box;
563
564 perlinterp_release ();
565
566 RETVAL = 0;
567
568 error = "JxlDecoderCreate failed";
569 if (!dec)
570 goto done;
571
572 runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
573
574 error = "JxlDecoderSetParallelRunner failed";
575 if (JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner) != JXL_DEC_SUCCESS)
576 goto done;
577
578 error = "JxlDecoderSubscribeEvents failed";
579 if (JxlDecoderSubscribeEvents (dec, JXL_DEC_FULL_IMAGE | JXL_DEC_BASIC_INFO) != JXL_DEC_SUCCESS)
580 goto done;
581
582 error = "JxlDecoderSetInput failed";
583 if (JxlDecoderSetInput (dec, next_in, avail_in) != JXL_DEC_SUCCESS)
584 goto done;
585
586 for (;;)
587 {
588 JxlDecoderStatus status = JxlDecoderProcessInput (dec);
589
590 printf ("status %d\n",status);
591
592 switch (status)
593 {
594 case JXL_DEC_ERROR:
595 error = "JxlDecoderProcessInput failed";
596 goto done;
597
598 case JXL_DEC_NEED_MORE_INPUT:
599 error = "incomplete file";
600 goto done;
601
602 case JXL_DEC_SUCCESS:
603 goto done;
604
605 case JXL_DEC_NEED_IMAGE_OUT_BUFFER:
606 {
607 error = "JxlDecoderGetBasicInfo failed";
608 if (JxlDecoderGetBasicInfo (dec, &info) != JXL_DEC_SUCCESS)
609 goto done;
610
611 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
612 error = "unable to allocate pixbuf";
613 if (!RETVAL)
614 goto done;
615
616 JxlPixelFormat format = {
617 info.alpha_bits ? 4 : 3,
618 JXL_TYPE_UINT8,
619 JXL_NATIVE_ENDIAN,
620 gdk_pixbuf_get_rowstride (RETVAL)
621 };
622
623 error = "JxlDecoderSetImageOutBuffer failed";
624 if (JxlDecoderSetImageOutBuffer (
625 dec,
626 &format,
627 gdk_pixbuf_get_pixels (RETVAL),
628 gdk_pixbuf_get_byte_length (RETVAL)
629 ) != JXL_DEC_SUCCESS)
630 goto done;
631 }
632 break;
633
634 default:
635 error = "unexpected event";
636 goto done;
637 }
638 }
639
640 done:
641 if (dec)
642 JxlDecoderDestroy (dec);
643
644 if (runner)
645 JxlThreadParallelRunnerDestroy (runner);
646
647 perlinterp_acquire ();
648
649 if (error)
650 {
651 if (RETVAL)
652 g_object_unref (RETVAL);
653
654 croak ("load_jxl: %s", error);
655 }
656#else
657 croak ("load_jxl: jpeg-xl not enabled at compile time");
658#endif
659}
660 OUTPUT:
661 RETVAL
662
663GdkPixbuf_noinc *
664decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
179 CODE: 665 CODE:
180{ 666{
181 struct jpeg_decompress_struct cinfo; 667 struct jpeg_decompress_struct cinfo;
182 struct jpg_err_mgr jerr; 668 struct jpg_err_mgr jerr;
183 guchar *data;
184 int rs; 669 int rs;
185 FILE *fp;
186 volatile GdkPixbuf *pb = 0; 670 volatile GdkPixbuf *pb = 0;
671 STRLEN data_len;
672 guchar *data = SvPVbyte (image_data, data_len);
673
187 RETVAL = 0; 674 RETVAL = 0;
188 675
189 if (!(fp = fopen (path, "rb"))) 676 perlinterp_release ();
190 XSRETURN_UNDEF;
191 677
192 cinfo.err = jpeg_std_error (&jerr.err); 678 cinfo.err = jpeg_std_error (&jerr.err);
193 679
194 jerr.err.error_exit = cv_error_exit; 680 jerr.err.error_exit = cv_error_exit;
195 jerr.err.output_message = cv_error_output; 681 jerr.err.output_message = cv_error_output;
196 682
197 if ((rs = setjmp (jerr.setjmp_buffer))) 683 if ((rs = setjmp (jerr.setjmp_buffer)))
198 { 684 {
199 fclose (fp);
200 jpeg_destroy_decompress (&cinfo); 685 jpeg_destroy_decompress (&cinfo);
201 686
202 if (pb) 687 if (pb)
203 g_object_unref ((gpointer)pb); 688 g_object_unref ((gpointer)pb);
204 689
690 perlinterp_acquire ();
205 XSRETURN_UNDEF; 691 XSRETURN_UNDEF;
206 } 692 }
207 693
694 if (!data_len)
695 longjmp (jerr.setjmp_buffer, 4);
696
208 jpeg_create_decompress (&cinfo); 697 jpeg_create_decompress (&cinfo);
698 cv_jpeg_mem_src (&cinfo, data, data_len);
209 699
210 jpeg_stdio_src (&cinfo, fp);
211 jpeg_read_header (&cinfo, TRUE); 700 jpeg_read_header (&cinfo, TRUE);
212 701
213 cinfo.dct_method = JDCT_DEFAULT; 702 cinfo.dct_method = JDCT_DEFAULT;
214 cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ 703 cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
215 cinfo.do_block_smoothing = FALSE; 704 cinfo.do_block_smoothing = FALSE;
225 { 714 {
226 cinfo.dct_method = JDCT_FASTEST; 715 cinfo.dct_method = JDCT_FASTEST;
227 cinfo.do_fancy_upsampling = FALSE; 716 cinfo.do_fancy_upsampling = FALSE;
228 717
229 while (cinfo.scale_denom < 8 718 while (cinfo.scale_denom < 8
230 && (cinfo.output_width >> 1) >= IW 719 && cinfo.output_width >= iw*4
231 && (cinfo.output_height >> 1) >= IH) 720 && cinfo.output_height >= ih*4)
232 { 721 {
233 cinfo.scale_denom <<= 1; 722 cinfo.scale_denom <<= 1;
234 jpeg_calc_output_dimensions (&cinfo); 723 jpeg_calc_output_dimensions (&cinfo);
235 } 724 }
236 } 725 }
237 726
727 if (cinfo.output_components != 3)
728 longjmp (jerr.setjmp_buffer, 3);
729
730 if (cinfo.jpeg_color_space == JCS_YCCK || cinfo.jpeg_color_space == JCS_CMYK)
731 {
732 cinfo.out_color_space = JCS_CMYK;
733 cinfo.output_components = 4;
734 }
735
238 pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height); 736 pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
239 if (!RETVAL) 737 if (!RETVAL)
240 longjmp (jerr.setjmp_buffer, 2); 738 longjmp (jerr.setjmp_buffer, 2);
241 739
242 data = gdk_pixbuf_get_pixels (RETVAL); 740 data = gdk_pixbuf_get_pixels (RETVAL);
243 rs = gdk_pixbuf_get_rowstride (RETVAL); 741 rs = gdk_pixbuf_get_rowstride (RETVAL);
244
245 if (cinfo.output_components != 3)
246 longjmp (jerr.setjmp_buffer, 3);
247 742
248 jpeg_start_decompress (&cinfo); 743 jpeg_start_decompress (&cinfo);
249 744
250 while (cinfo.output_scanline < cinfo.output_height) 745 while (cinfo.output_scanline < cinfo.output_height)
251 { 746 {
258 rp [3] = (guchar *)rp [2] + rs; 753 rp [3] = (guchar *)rp [2] + rs;
259 754
260 jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); 755 jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
261 } 756 }
262 757
758 if (cinfo.out_color_space == JCS_CMYK)
759 {
760 guchar *end = data + cinfo.output_height * rs;
761
762 while (data < end)
763 {
764 U32 c = data [0];
765 U32 m = data [1];
766 U32 y = data [2];
767 U32 k = data [3];
768
769 if (0)
770 if (cinfo.Adobe_transform == 2)
771 {
772 c ^= 0xff;
773 m ^= 0xff;
774 y ^= 0xff;
775 k ^= 0xff;
776 }
777
778 data [0] = (c * k + 0x80) / 0xff;
779 data [1] = (m * k + 0x80) / 0xff;
780 data [2] = (y * k + 0x80) / 0xff;
781 data [3] = 0xff;
782
783 data += 4;
784 }
785 }
786
263 jpeg_finish_decompress (&cinfo); 787 jpeg_finish_decompress (&cinfo);
264 fclose (fp);
265 jpeg_destroy_decompress (&cinfo); 788 jpeg_destroy_decompress (&cinfo);
789 perlinterp_acquire ();
266} 790}
267 OUTPUT: 791 OUTPUT:
268 RETVAL 792 RETVAL
269 793
794void
795compare (GdkPixbuf *a, GdkPixbuf *b)
796 PPCODE:
797 perlinterp_release ();
798{
799 int w = gdk_pixbuf_get_width (a);
800 int h = gdk_pixbuf_get_height (a);
801
802 int sa = gdk_pixbuf_get_rowstride (a);
803 int sb = gdk_pixbuf_get_rowstride (b);
804
805 guchar *pa = gdk_pixbuf_get_pixels (a);
806 guchar *pb = gdk_pixbuf_get_pixels (b);
807
808 int x, y;
809
810 assert (w == gdk_pixbuf_get_width (b));
811 assert (h == gdk_pixbuf_get_height (b));
812
813 assert (gdk_pixbuf_get_n_channels (a) == 3);
814 assert (gdk_pixbuf_get_n_channels (b) == 3);
815
816 double diff = 0.;
817 int peak = 0;
818
819 if (w && h)
820 for (y = 0; y < h; y++)
821 {
822 guchar *pa_ = pa + y * sa;
823 guchar *pb_ = pb + y * sb;
824
825 for (x = 0; x < w; x++)
826 {
827 int d;
828
829 d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d));
830 d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d));
831 d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d));
832 }
833 }
834
835 perlinterp_acquire ();
836
837 EXTEND (SP, 2);
838 PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
839 PUSHs (sv_2mortal (newSVnv (peak / 255.)));
840}
841
270############################################################################# 842#############################################################################
271 843
272MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer 844MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
845
846# currently only works for filenames (octet strings)
273 847
274SV * 848SV *
275foldcase (SV *pathsv) 849foldcase (SV *pathsv)
276 PROTOTYPE: $ 850 PROTOTYPE: $
277 CODE: 851 CODE:
278{ 852{
279 STRLEN plen; 853 STRLEN plen;
280 U8 *path = SvPVutf8 (pathsv, plen); 854 U8 *path = (U8 *)SvPV (pathsv, plen);
281 U8 *pend = path + plen; 855 U8 *pend = path + plen;
282 U8 dst [plen * 6 * 3], *dstp = dst; 856 U8 dst [plen * 6 * 3], *dstp = dst;
283 857
284 while (path < pend) 858 while (path < pend)
285 { 859 {
860 U8 ch = *path;
861
286 if (*path >= '0' && *path <= '9') 862 if (ch >= 'a' && ch <= 'z')
863 *dstp++ = *path++;
864 else if (ch >= 'A' && ch <= 'Z')
865 *dstp++ = *path++ + ('a' - 'A');
866 else if (ch >= '0' && ch <= '9')
287 { 867 {
288 STRLEN el, nl = 0; 868 STRLEN el, nl = 0;
289 while (*path >= '0' && *path <= '9') 869 while (*path >= '0' && *path <= '9' && path < pend)
290 path++, nl++; 870 path++, nl++;
291 871
292 for (el = nl; el < 6; el++) 872 for (el = nl; el < 6; el++)
293 *dstp++ = '0'; 873 *dstp++ = '0';
294 874
295 memcpy (dstp, path - nl, nl); 875 memcpy (dstp, path - nl, nl);
296 dstp += nl; 876 dstp += nl;
297 } 877 }
878 else
879 *dstp++ = *path++;
880#if 0
298 else 881 else
299 { 882 {
300 STRLEN cl; 883 STRLEN cl;
301 to_utf8_fold (path, dstp, &cl); 884 to_utf8_fold (path, dstp, &cl);
302 dstp += cl; 885 dstp += cl;
303 path += is_utf8_char (path); 886 path += is_utf8_char (path);
304 } 887 }
888#endif
305 } 889 }
306 890
307 RETVAL = newSVpvn (dst, dstp - dst); 891 RETVAL = newSVpvn ((const char *)dst, dstp - dst);
308} 892}
309 OUTPUT: 893 OUTPUT:
310 RETVAL 894 RETVAL
311 895
312GdkPixbuf_noinc * 896GdkPixbuf_noinc *
313p7_to_pb (int w, int h, guchar *src) 897p7_to_pb (int w, int h, SV *src_sv)
898 PROTOTYPE: @
314 CODE: 899 CODE:
315{ 900{
316 int x, y; 901 int x, y;
317 guchar *dst, *d; 902 guchar *dst, *d;
318 int dstr; 903 int dstr;
904 guchar *src = (guchar *)SvPVbyte_nolen (src_sv);
319 905
320 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h); 906 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
321 dst = gdk_pixbuf_get_pixels (RETVAL); 907 dst = gdk_pixbuf_get_pixels (RETVAL);
322 dstr = gdk_pixbuf_get_rowstride (RETVAL); 908 dstr = gdk_pixbuf_get_rowstride (RETVAL);
323 909
332 } 918 }
333} 919}
334 OUTPUT: 920 OUTPUT:
335 RETVAL 921 RETVAL
336 922
923#############################################################################
924
925MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
926
927void
928dump_ascii85 (PerlIO *fp, GdkPixbuf *pb)
929 CODE:
930{
931 int w = gdk_pixbuf_get_width (pb);
932 int h = gdk_pixbuf_get_height (pb);
933 int x, y, i;
934 guchar *dst;
935 int bpp = gdk_pixbuf_get_n_channels (pb);
936 guchar *src = gdk_pixbuf_get_pixels (pb);
937 int sstr = gdk_pixbuf_get_rowstride (pb);
938
939 a85_init ();
940
941 for (y = 0; y < h; y++)
942 for (x = 0; x < w; x++)
943 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
944 a85_push (fp, src [x * bpp + y * sstr + i]);
945
946 a85_finish (fp);
947}
948
949void
950dump_binary (PerlIO *fp, GdkPixbuf *pb)
951 CODE:
952{
953 int w = gdk_pixbuf_get_width (pb);
954 int h = gdk_pixbuf_get_height (pb);
955 int x, y, i;
956 guchar *dst;
957 int bpp = gdk_pixbuf_get_n_channels (pb);
958 guchar *src = gdk_pixbuf_get_pixels (pb);
959 int sstr = gdk_pixbuf_get_rowstride (pb);
960
961 for (y = 0; y < h; y++)
962 for (x = 0; x < w; x++)
963 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
964 PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
965}
966
967#############################################################################
968
969MODULE = Gtk2::CV PACKAGE = Gtk2::CV
970
337SV * 971SV *
338pb_to_p7 (GdkPixbuf *pb) 972pb_to_hv84 (GdkPixbuf *pb)
339 CODE: 973 CODE:
340{ 974{
341 int w = gdk_pixbuf_get_width (pb); 975 int w = gdk_pixbuf_get_width (pb);
342 int h = gdk_pixbuf_get_height (pb); 976 int h = gdk_pixbuf_get_height (pb);
343 int x, y; 977 int x, y;
344 guchar *dst; 978 guchar *dst;
345 int bpp = gdk_pixbuf_get_n_channels (pb); 979 int bpp = gdk_pixbuf_get_n_channels (pb);
346 guchar *src = gdk_pixbuf_get_pixels (pb); 980 guchar *src = gdk_pixbuf_get_pixels (pb);
347 int sstr = gdk_pixbuf_get_rowstride (pb); 981 int sstr = gdk_pixbuf_get_rowstride (pb);
348 int Er[IW], Eg[IW], Eb[IW];
349 int seed = 77;
350
351 RETVAL = newSV (w * h);
352 SvPOK_only (RETVAL);
353 SvCUR_set (RETVAL, w * h);
354
355 dst = SvPVX (RETVAL);
356
357 memset (Er, 0, sizeof (int) * IW);
358 memset (Eg, 0, sizeof (int) * IW);
359 memset (Eb, 0, sizeof (int) * IW);
360
361 /* some primitive error distribution + random dithering */
362
363 for (y = 0; y < h; y++)
364 {
365 int er = 0, eg = 0, eb = 0;
366
367 for (x = 0; x < w; x++)
368 {
369 int r, g, b;
370 guchar *p = src + x * bpp + y * sstr;
371
372 r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255;
373 g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255;
374 b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255;
375
376 r = r > 7 ? 7 : r < 0 ? 0 : r;
377 g = g > 7 ? 7 : g < 0 ? 0 : g;
378 b = b > 3 ? 3 : b < 0 ? 0 : b;
379
380 er += p[0] - (r * 255 + 4) / 7;
381 eg += p[1] - (g * 255 + 4) / 7;
382 eb += p[2] - (b * 255 + 2) / 3;
383
384 Er[x] = er >> 1; er -= (er + 1) >> 1;
385 Eg[x] = eg >> 1; eg -= (eg + 1) >> 1;
386 Eb[x] = eb >> 1; eb -= (eb + 1) >> 1;
387
388 *dst++ = r << 5 | g << 2 | b;
389 }
390 }
391}
392 OUTPUT:
393 RETVAL
394
395SV *
396make_histogram (SV *ar)
397 CODE:
398{
399 int i;
400 AV *av, *result;
401
402 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
403 croak ("Not an array ref as first argument to make_histogram");
404
405 av = (AV *) SvRV (ar);
406 result = newAV ();
407
408 for (i = 0; i <= av_len (av); ++i)
409 {
410 const int HISTSIZE = 64;
411
412 int j;
413 SV *sv = *av_fetch (av, i, 1);
414 STRLEN len;
415 char *buf = SvPVbyte (sv, len);
416
417 int tmphist[HISTSIZE];
418 float *hist;
419
420 SV *histsv = newSV (HISTSIZE * sizeof (float) + 1);
421 SvPOK_on (histsv);
422 SvCUR_set (histsv, HISTSIZE * sizeof (float));
423 hist = (float *)SvPVX (histsv);
424
425 Zero (tmphist, sizeof (tmphist), char);
426
427 for (j = len; j--; )
428 {
429 unsigned int idx
430 = ((*buf & 0xc0) >> 2)
431 | ((*buf & 0x18) >> 1)
432 | (*buf & 0x03);
433
434 ++tmphist[idx];
435 ++buf;
436 }
437
438 for (j = 0; j < HISTSIZE; ++j)
439 hist[j] = (float)tmphist[j] / (len + 1e-30);
440
441 av_push (result, histsv);
442 }
443
444 RETVAL = newRV_noinc ((SV *)result);
445}
446 OUTPUT:
447 RETVAL
448
449#############################################################################
450
451MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
452
453void
454dump_ascii85 (PerlIO *fp, GdkPixbuf *pb)
455 CODE:
456{
457 int w = gdk_pixbuf_get_width (pb);
458 int h = gdk_pixbuf_get_height (pb);
459 int x, y, i;
460 guchar *dst;
461 int bpp = gdk_pixbuf_get_n_channels (pb);
462 guchar *src = gdk_pixbuf_get_pixels (pb);
463 int sstr = gdk_pixbuf_get_rowstride (pb);
464
465 a85_init ();
466
467 for (y = 0; y < h; y++)
468 for (x = 0; x < w; x++)
469 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
470 a85_push (fp, src [x * bpp + y * sstr + i]);
471
472 a85_finish (fp);
473}
474
475void
476dump_binary (PerlIO *fp, GdkPixbuf *pb)
477 CODE:
478{
479 int w = gdk_pixbuf_get_width (pb);
480 int h = gdk_pixbuf_get_height (pb);
481 int x, y, i;
482 guchar *dst;
483 int bpp = gdk_pixbuf_get_n_channels (pb);
484 guchar *src = gdk_pixbuf_get_pixels (pb);
485 int sstr = gdk_pixbuf_get_rowstride (pb);
486
487 for (y = 0; y < h; y++)
488 for (x = 0; x < w; x++)
489 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
490 PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
491}
492
493#############################################################################
494
495MODULE = Gtk2::CV PACKAGE = Gtk2::CV
496
497SV *
498pb_to_hv84 (GdkPixbuf *pb)
499 CODE:
500{
501 int w = gdk_pixbuf_get_width (pb);
502 int h = gdk_pixbuf_get_height (pb);
503 int x, y;
504 guchar *dst;
505 int bpp = gdk_pixbuf_get_n_channels (pb);
506 guchar *src = gdk_pixbuf_get_pixels (pb);
507 int sstr = gdk_pixbuf_get_rowstride (pb);
508 982
509 RETVAL = newSV (6 * 8 * 12 / 8); 983 RETVAL = newSV (6 * 8 * 12 / 8);
510 SvPOK_only (RETVAL); 984 SvPOK_only (RETVAL);
511 SvCUR_set (RETVAL, 6 * 8 * 12 / 8); 985 SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
512 986
513 dst = SvPVX (RETVAL); 987 dst = (guchar *)SvPVX (RETVAL);
514 988
515 /* some primitive error distribution + random dithering */ 989 /* some primitive error distribution + random dithering */
516 990
517 for (y = 0; y < h; y++) 991 for (y = 0; y < h; y++)
518 { 992 {
576 } 1050 }
577} 1051}
578 OUTPUT: 1052 OUTPUT:
579 RETVAL 1053 RETVAL
580 1054
1055#############################################################################
581 1056
1057MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
582 1058
1059SV *
1060extract_features (SV *ar)
1061 CODE:
1062{
1063 int i;
1064 AV *av, *result;
1065
1066 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
1067 croak ("Not an array ref as first argument to extract_features");
1068
1069 av = (AV *) SvRV (ar);
1070 result = newAV ();
1071
1072 for (i = 0; i <= av_len (av); ++i)
1073 {
1074 SV *sv = *av_fetch (av, i, 1);
1075 SV *histsv = newSV (9 * sizeof (float) + 1);
1076
1077 SvPOK_on (histsv);
1078 SvCUR_set (histsv, 9 * sizeof (float));
1079 float *hist = (float *)SvPVX (histsv);
1080
1081 struct feature f_h, f_s, f_v;
1082 feature_init (&f_h);
1083 feature_init (&f_s);
1084 feature_init (&f_v);
1085
1086 {
1087 STRLEN len;
1088 unsigned char *buf = (unsigned char *)SvPVbyte (sv, len);
1089 while (len >= 3)
1090 {
1091 unsigned int r, g, b, h, s, v;
1092 r = *buf++; g = *buf++; b = *buf++;
1093 rgb_to_hsv (r, g, b, &h, &s, &v);
1094
1095 feature_update_pass_1 (&f_h, h);
1096 feature_update_pass_1 (&f_s, s);
1097 feature_update_pass_1 (&f_v, v);
1098
1099 len -= 3;
1100 }
1101
1102 feature_finish_pass_1 (&f_h);
1103 feature_finish_pass_1 (&f_s);
1104 feature_finish_pass_1 (&f_v);
1105 }
1106
1107 {
1108 STRLEN len;
1109 unsigned char *buf = (unsigned char *)SvPVbyte (sv, len);
1110 while (len >= 3)
1111 {
1112 unsigned int r, g, b, h, s, v;
1113 r = *buf++; g = *buf++; b = *buf++;
1114 rgb_to_hsv (r, g, b, &h, &s, &v);
1115
1116 feature_update_pass_2 (&f_h, h);
1117 feature_update_pass_2 (&f_s, s);
1118 feature_update_pass_2 (&f_v, v);
1119
1120 len -= 3;
1121 }
1122
1123 feature_finish_pass_2 (&f_h);
1124 feature_finish_pass_2 (&f_s);
1125 feature_finish_pass_2 (&f_v);
1126 }
1127
1128 hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
1129 hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
1130 hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
1131
1132 av_push (result, histsv);
1133 }
1134
1135 RETVAL = newRV_noinc ((SV *)result);
1136}
1137 OUTPUT:
1138 RETVAL
1139

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines