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

Comparing CV/CV.xs (file contents):
Revision 1.25 by root, Fri Aug 19 00:34:54 2005 UTC vs.
Revision 1.61 by root, Sun Jul 25 11:24:48 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
9#include <glib.h> 14#include <glib.h>
10#include <gtk/gtk.h> 15#include <gtk/gtk.h>
16#include <gdk/gdkx.h>
11#include <gdk-pixbuf/gdk-pixbuf.h> 17#include <gdk-pixbuf/gdk-pixbuf.h>
12 18
13#include <gperl.h> 19#include <gperl.h>
14#include <gtk2perl.h> 20#include <gtk2perl.h>
15 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
16#define IW 80 /* MUST match Schnauer.pm! */ 36#define IW 80 /* MUST match Schnauzer.pm! */
17#define IH 60 /* MUST match Schnauer.pm! */ 37#define IH 60 /* MUST match Schnauzer.pm! */
18
19#define RAND (seed = (seed + 7141) * 54773 % 134456)
20 38
21#define LINELENGTH 240 39#define LINELENGTH 240
22 40
23#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 */
24 46
25struct jpg_err_mgr 47struct jpg_err_mgr
26{ 48{
27 struct jpeg_error_mgr err; 49 struct jpeg_error_mgr err;
28 jmp_buf setjmp_buffer; 50 jmp_buf setjmp_buffer;
29}; 51};
30 52
31static void 53static void
32cv_error_exit (j_common_ptr cinfo) 54cv_error_exit (j_common_ptr cinfo)
33{ 55{
56 cinfo->err->output_message (cinfo);
34 longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99); 57 longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
35} 58}
36 59
37static void 60static void
38cv_error_output (j_common_ptr cinfo) 61cv_error_output (j_common_ptr cinfo)
39{ 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);
40 return; 68 return;
41} 69}
42 70
43static void 71static void
44rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, 72rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
65 93
66 *h &= 255; 94 *h &= 255;
67 } 95 }
68} 96}
69 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
70static guint32 a85_val; 153static guint32 a85_val;
71static guint a85_cnt; 154static guint a85_cnt;
72static guchar a85_buf[LINELENGTH], *a85_ptr; 155static guchar a85_buf[LINELENGTH], *a85_ptr;
73 156
74static void 157static void
86 if (!--a85_cnt) 169 if (!--a85_cnt)
87 { 170 {
88 a85_cnt = 4; 171 a85_cnt = 4;
89 if (a85_val) 172 if (a85_val)
90 { 173 {
91 a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; 174 a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
92 a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; 175 a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
93 a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; 176 a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
94 a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; 177 a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
95 a85_ptr[0] = (a85_val ) + 33; 178 a85_ptr[0] = (a85_val ) + 33;
96 179
97 a85_ptr += 5; 180 a85_ptr += 5;
98 } 181 }
121 204
122 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); 205 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
123} 206}
124 207
125///////////////////////////////////////////////////////////////////////////// 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/////////////////////////////////////////////////////////////////////////////
126 303
127MODULE = Gtk2::CV PACKAGE = Gtk2::CV 304MODULE = Gtk2::CV PACKAGE = Gtk2::CV
128 305
129PROTOTYPES: ENABLE 306PROTOTYPES: ENABLE
130 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
131# missing in Gtk2 perl module 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);
132 363
133gboolean 364gboolean
134gdk_net_wm_supports (GdkAtom property) 365gdk_net_wm_supports (GdkAtom property)
135 CODE: 366 CODE:
136#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) 367#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
140#endif 371#endif
141 OUTPUT: 372 OUTPUT:
142 RETVAL 373 RETVAL
143 374
144GdkPixbuf_noinc * 375GdkPixbuf_noinc *
145transpose (GdkPixbuf *pb) 376dealpha_expose (GdkPixbuf *pb)
146 CODE: 377 CODE:
378 perlinterp_release ();
147{ 379{
148 int w = gdk_pixbuf_get_width (pb); 380 int w = gdk_pixbuf_get_width (pb);
149 int h = gdk_pixbuf_get_height (pb); 381 int h = gdk_pixbuf_get_height (pb);
150 int bpp = gdk_pixbuf_get_n_channels (pb); 382 int bpp = gdk_pixbuf_get_n_channels (pb);
151 int x, y, i; 383 int x, y, i;
152 guchar *src = gdk_pixbuf_get_pixels (pb), *dst; 384 guchar *src = gdk_pixbuf_get_pixels (pb), *dst;
153 int sstr = gdk_pixbuf_get_rowstride (pb), dstr; 385 int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
154 386
155 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w); 387 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
156 388
157 dst = gdk_pixbuf_get_pixels (RETVAL); 389 dst = gdk_pixbuf_get_pixels (RETVAL);
158 dstr = gdk_pixbuf_get_rowstride (RETVAL); 390 dstr = gdk_pixbuf_get_rowstride (RETVAL);
159 391
160 for (y = 0; y < h; y++)
161 for (x = 0; x < w; x++) 392 for (x = 0; x < w; x++)
393 for (y = 0; y < h; y++)
162 for (i = 0; i < bpp; i++) 394 for (i = 0; i < 3; i++)
163 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];
164} 396}
397 perlinterp_acquire ();
165 OUTPUT: 398 OUTPUT:
166 RETVAL 399 RETVAL
167 400
168GdkPixbuf_noinc * 401GdkPixbuf_noinc *
169flop (GdkPixbuf *pb) 402rotate (GdkPixbuf *pb, int angle)
170 CODE: 403 CODE:
171{ 404 perlinterp_release ();
172 int w = gdk_pixbuf_get_width (pb); 405 if (angle < 0)
173 int h = gdk_pixbuf_get_height (pb); 406 angle += 360;
174 int bpp = gdk_pixbuf_get_n_channels (pb); 407 RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
175 int x, y, i; 408 : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
176 guchar *src = gdk_pixbuf_get_pixels (pb), *dst; 409 : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
177 int sstr = gdk_pixbuf_get_rowstride (pb), dstr; 410 : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
178 411 : angle);
179 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h); 412 perlinterp_acquire ();
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[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
188}
189 OUTPUT: 413 OUTPUT:
190 RETVAL 414 RETVAL
191 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 if (data_len >= 13
457 && data[0] == 'G'
458 && data[1] == 'I'
459 && data[2] == 'F'
460 && data[3] == '8'
461 //&& (data[4] == '7' || data[4] == '9')
462 && data[5] == 'a')
463 {
464 RETVAL = "image/gif";
465
466 // now see if its animated - we require the netscape application header for this
467 int ofs = 13;
468
469 if (data[10] & 0x80)
470 ofs += (1 << ((data[10] & 7) + 1)) * 3;
471
472 if (data_len >= ofs + 2 + 1 + 11)
473 {
474 // skip a rare first graphic control extension block
475 // we assume there is exactly one block
476 if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9)
477 ofs += 3 + data[ofs + 2] + 1;
478
479 if (data_len >= ofs + 2 + 1 + 11)
480 if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1))
481 RETVAL = "video/gif";
482 }
483 }
484
485 else
486 XSRETURN_UNDEF;
487}
488 OUTPUT:
489 RETVAL
490
192GdkPixbuf_noinc * 491GdkPixbuf_noinc *
193load_jpeg (SV *path, int thumbnail=0) 492decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
493 CODE:
494{
495#if WEBP
496 STRLEN data_size;
497 int alpha;
498 WebPData data;
499 WebPDemuxer *demux;
500 WebPIterator iter;
501 WebPDecoderConfig config;
502 int inw, inh;
503
504 data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
505 data.size = data_size;
506
507 perlinterp_release ();
508
509 RETVAL = 0;
510
511 if (!(demux = WebPDemux (&data)))
512 goto err_demux;
513
514 if (!WebPDemuxGetFrame (demux, 1, &iter))
515 goto err_iter;
516
517 if (!WebPInitDecoderConfig (&config))
518 goto err_iter;
519
520 config.options.use_threads = 1;
521
522 if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
523 goto err_iter;
524
525 inw = config.input.width;
526 inh = config.input.height;
527
528 if (thumbnail)
529 {
530 if (inw * ih > inh * iw)
531 ih = (iw * inh + inw - 1) / inw;
532 else
533 iw = (ih * inw + inh - 1) / inh;
534
535 config.options.bypass_filtering = 1;
536 config.options.no_fancy_upsampling = 1;
537
538 config.options.use_scaling = 1;
539 config.options.scaled_width = iw;
540 config.options.scaled_height = ih;
541 }
542 else
543 {
544 iw = inw;
545 ih = inh;
546 }
547
548 alpha = !!config.input.has_alpha;
549
550 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
551 if (!RETVAL)
552 goto err_iter;
553
554 config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
555 config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
556 config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
557 config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
558 config.output.is_external_memory = 1;
559
560 if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
561 {
562 g_object_unref (RETVAL);
563 RETVAL = 0;
564 goto err_iter;
565 }
566
567 err_iter:
568 WebPDemuxReleaseIterator (&iter);
569 err_demux:
570 WebPDemuxDelete (demux);
571
572 perlinterp_acquire ();
573#else
574 croak ("load_webp: webp not enabled at compile time");
575#endif
576}
577 OUTPUT:
578 RETVAL
579
580GdkPixbuf_noinc *
581decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
582 CODE:
583{
584#if JXL
585 JxlDecoder *dec = JxlDecoderCreate (0);
586 JxlBasicInfo info;
587 const uint8_t *next_in = (uint8_t *)SvPVbyte_nolen (image_data);
588 size_t avail_in = SvCUR (image_data);
589 const char *error = 0;
590 void *runner = 0;
591 struct bmff_box box;
592
593 perlinterp_release ();
594
595 RETVAL = 0;
596
597 error = "JxlDecoderCreate failed";
598 if (!dec)
599 goto done;
600
601 runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
602
603 error = "JxlDecoderSetParallelRunner failed";
604 if (JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner) != JXL_DEC_SUCCESS)
605 goto done;
606
607 error = "JxlDecoderSubscribeEvents failed";
608 if (JxlDecoderSubscribeEvents (dec, JXL_DEC_FULL_IMAGE | JXL_DEC_BASIC_INFO) != JXL_DEC_SUCCESS)
609 goto done;
610
611 error = "JxlDecoderSetInput failed";
612 if (JxlDecoderSetInput (dec, next_in, avail_in) != JXL_DEC_SUCCESS)
613 goto done;
614
615 for (;;)
616 {
617 JxlDecoderStatus status = JxlDecoderProcessInput (dec);
618
619 printf ("status %d\n",status);
620
621 switch (status)
622 {
623 case JXL_DEC_ERROR:
624 error = "JxlDecoderProcessInput failed";
625 goto done;
626
627 case JXL_DEC_NEED_MORE_INPUT:
628 error = "incomplete file";
629 goto done;
630
631 case JXL_DEC_SUCCESS:
632 goto done;
633
634 case JXL_DEC_NEED_IMAGE_OUT_BUFFER:
635 {
636 error = "JxlDecoderGetBasicInfo failed";
637 if (JxlDecoderGetBasicInfo (dec, &info) != JXL_DEC_SUCCESS)
638 goto done;
639
640 RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
641 error = "unable to allocate pixbuf";
642 if (!RETVAL)
643 goto done;
644
645 JxlPixelFormat format = {
646 info.alpha_bits ? 4 : 3,
647 JXL_TYPE_UINT8,
648 JXL_NATIVE_ENDIAN,
649 gdk_pixbuf_get_rowstride (RETVAL)
650 };
651
652 error = "JxlDecoderSetImageOutBuffer failed";
653 if (JxlDecoderSetImageOutBuffer (
654 dec,
655 &format,
656 gdk_pixbuf_get_pixels (RETVAL),
657 gdk_pixbuf_get_byte_length (RETVAL)
658 ) != JXL_DEC_SUCCESS)
659 goto done;
660 }
661 break;
662
663 default:
664 error = "unexpected event";
665 goto done;
666 }
667 }
668
669 done:
670 if (dec)
671 JxlDecoderDestroy (dec);
672
673 if (runner)
674 JxlThreadParallelRunnerDestroy (runner);
675
676 perlinterp_acquire ();
677
678 if (error)
679 {
680 if (RETVAL)
681 g_object_unref (RETVAL);
682
683 croak ("load_jxl: %s", error);
684 }
685#else
686 croak ("load_jxl: jpeg-xl not enabled at compile time");
687#endif
688}
689 OUTPUT:
690 RETVAL
691
692GdkPixbuf_noinc *
693decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
194 CODE: 694 CODE:
195{ 695{
196 struct jpeg_decompress_struct cinfo; 696 struct jpeg_decompress_struct cinfo;
197 struct jpg_err_mgr jerr; 697 struct jpg_err_mgr jerr;
198 guchar *data;
199 int rs; 698 int rs;
200 FILE *fp;
201 volatile GdkPixbuf *pb = 0; 699 volatile GdkPixbuf *pb = 0;
202 gchar *filename; 700 STRLEN data_len;
701 guchar *data = SvPVbyte (image_data, data_len);
203 702
204 RETVAL = 0; 703 RETVAL = 0;
205 704
206 filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0); 705 perlinterp_release ();
207 fp = fopen (filename, "rb");
208 g_free (filename);
209
210 if (!fp)
211 XSRETURN_UNDEF;
212 706
213 cinfo.err = jpeg_std_error (&jerr.err); 707 cinfo.err = jpeg_std_error (&jerr.err);
214 708
215 jerr.err.error_exit = cv_error_exit; 709 jerr.err.error_exit = cv_error_exit;
216 jerr.err.output_message = cv_error_output; 710 jerr.err.output_message = cv_error_output;
217 711
218 if ((rs = setjmp (jerr.setjmp_buffer))) 712 if ((rs = setjmp (jerr.setjmp_buffer)))
219 { 713 {
220 fclose (fp);
221 jpeg_destroy_decompress (&cinfo); 714 jpeg_destroy_decompress (&cinfo);
222 715
223 if (pb) 716 if (pb)
224 g_object_unref ((gpointer)pb); 717 g_object_unref ((gpointer)pb);
225 718
719 perlinterp_acquire ();
226 XSRETURN_UNDEF; 720 XSRETURN_UNDEF;
227 } 721 }
228 722
723 if (!data_len)
724 longjmp (jerr.setjmp_buffer, 4);
725
229 jpeg_create_decompress (&cinfo); 726 jpeg_create_decompress (&cinfo);
727 cv_jpeg_mem_src (&cinfo, data, data_len);
230 728
231 jpeg_stdio_src (&cinfo, fp);
232 jpeg_read_header (&cinfo, TRUE); 729 jpeg_read_header (&cinfo, TRUE);
233 730
234 cinfo.dct_method = JDCT_DEFAULT; 731 cinfo.dct_method = JDCT_DEFAULT;
235 cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ 732 cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
236 cinfo.do_block_smoothing = FALSE; 733 cinfo.do_block_smoothing = FALSE;
246 { 743 {
247 cinfo.dct_method = JDCT_FASTEST; 744 cinfo.dct_method = JDCT_FASTEST;
248 cinfo.do_fancy_upsampling = FALSE; 745 cinfo.do_fancy_upsampling = FALSE;
249 746
250 while (cinfo.scale_denom < 8 747 while (cinfo.scale_denom < 8
251 && cinfo.output_width >= IW*4 748 && cinfo.output_width >= iw*4
252 && cinfo.output_height >= IH*4) 749 && cinfo.output_height >= ih*4)
253 { 750 {
254 cinfo.scale_denom <<= 1; 751 cinfo.scale_denom <<= 1;
255 jpeg_calc_output_dimensions (&cinfo); 752 jpeg_calc_output_dimensions (&cinfo);
256 } 753 }
257 } 754 }
258 755
756 if (cinfo.output_components != 3)
757 longjmp (jerr.setjmp_buffer, 3);
758
759 if (cinfo.jpeg_color_space == JCS_YCCK || cinfo.jpeg_color_space == JCS_CMYK)
760 {
761 cinfo.out_color_space = JCS_CMYK;
762 cinfo.output_components = 4;
763 }
764
259 pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height); 765 pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
260 if (!RETVAL) 766 if (!RETVAL)
261 longjmp (jerr.setjmp_buffer, 2); 767 longjmp (jerr.setjmp_buffer, 2);
262 768
263 data = gdk_pixbuf_get_pixels (RETVAL); 769 data = gdk_pixbuf_get_pixels (RETVAL);
264 rs = gdk_pixbuf_get_rowstride (RETVAL); 770 rs = gdk_pixbuf_get_rowstride (RETVAL);
265
266 if (cinfo.output_components != 3)
267 longjmp (jerr.setjmp_buffer, 3);
268 771
269 jpeg_start_decompress (&cinfo); 772 jpeg_start_decompress (&cinfo);
270 773
271 while (cinfo.output_scanline < cinfo.output_height) 774 while (cinfo.output_scanline < cinfo.output_height)
272 { 775 {
279 rp [3] = (guchar *)rp [2] + rs; 782 rp [3] = (guchar *)rp [2] + rs;
280 783
281 jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); 784 jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
282 } 785 }
283 786
787 if (cinfo.out_color_space == JCS_CMYK)
788 {
789 guchar *end = data + cinfo.output_height * rs;
790
791 while (data < end)
792 {
793 U32 c = data [0];
794 U32 m = data [1];
795 U32 y = data [2];
796 U32 k = data [3];
797
798 if (0)
799 if (cinfo.Adobe_transform == 2)
800 {
801 c ^= 0xff;
802 m ^= 0xff;
803 y ^= 0xff;
804 k ^= 0xff;
805 }
806
807 data [0] = (c * k + 0x80) / 0xff;
808 data [1] = (m * k + 0x80) / 0xff;
809 data [2] = (y * k + 0x80) / 0xff;
810 data [3] = 0xff;
811
812 data += 4;
813 }
814 }
815
284 jpeg_finish_decompress (&cinfo); 816 jpeg_finish_decompress (&cinfo);
285 fclose (fp);
286 jpeg_destroy_decompress (&cinfo); 817 jpeg_destroy_decompress (&cinfo);
818 perlinterp_acquire ();
287} 819}
288 OUTPUT: 820 OUTPUT:
289 RETVAL 821 RETVAL
290 822
823void
824compare (GdkPixbuf *a, GdkPixbuf *b)
825 PPCODE:
826 perlinterp_release ();
827{
828 int w = gdk_pixbuf_get_width (a);
829 int h = gdk_pixbuf_get_height (a);
830
831 int sa = gdk_pixbuf_get_rowstride (a);
832 int sb = gdk_pixbuf_get_rowstride (b);
833
834 guchar *pa = gdk_pixbuf_get_pixels (a);
835 guchar *pb = gdk_pixbuf_get_pixels (b);
836
837 int x, y;
838
839 assert (w == gdk_pixbuf_get_width (b));
840 assert (h == gdk_pixbuf_get_height (b));
841
842 assert (gdk_pixbuf_get_n_channels (a) == 3);
843 assert (gdk_pixbuf_get_n_channels (b) == 3);
844
845 double diff = 0.;
846 int peak = 0;
847
848 if (w && h)
849 for (y = 0; y < h; y++)
850 {
851 guchar *pa_ = pa + y * sa;
852 guchar *pb_ = pb + y * sb;
853
854 for (x = 0; x < w; x++)
855 {
856 int d;
857
858 d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d));
859 d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d));
860 d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d));
861 }
862 }
863
864 perlinterp_acquire ();
865
866 EXTEND (SP, 2);
867 PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
868 PUSHs (sv_2mortal (newSVnv (peak / 255.)));
869}
870
291############################################################################# 871#############################################################################
292 872
293MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer 873MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
874
875# currently only works for filenames (octet strings)
294 876
295SV * 877SV *
296foldcase (SV *pathsv) 878foldcase (SV *pathsv)
297 PROTOTYPE: $ 879 PROTOTYPE: $
298 CODE: 880 CODE:
299{ 881{
300 STRLEN plen; 882 STRLEN plen;
301 U8 *path = (U8 *)SvPVutf8 (pathsv, plen); 883 U8 *path = (U8 *)SvPV (pathsv, plen);
302 U8 *pend = path + plen; 884 U8 *pend = path + plen;
303 U8 dst [plen * 6 * 3], *dstp = dst; 885 U8 dst [plen * 8 * 3], *dstp = dst;
304 886
305 while (path < pend) 887 while (path < pend)
306 { 888 {
307 U8 ch = *path; 889 U8 ch = *path;
308 890
309 if (ch >= 'a' && ch <= 'z') 891 if (ch >= 'a' && ch <= 'z')
310 *dstp++ = *path++; 892 *dstp++ = *path++;
893 else if (ch >= 'A' && ch <= 'Z')
894 *dstp++ = *path++ + ('a' - 'A');
311 else if (ch >= '0' && ch <= '9') 895 else if (ch >= '0' && ch <= '9')
312 { 896 {
897 /* version sort, up to 8 digits */
313 STRLEN el, nl = 0; 898 STRLEN el, nl = 0;
314 while (*path >= '0' && *path <= '9' && path < pend) 899 while (*path >= '0' && *path <= '9' && path < pend)
315 path++, nl++; 900 path++, nl++;
316 901
317 for (el = nl; el < 6; el++) 902 for (el = nl; el < 8; el++)
318 *dstp++ = '0'; 903 *dstp++ = '0';
319 904
320 memcpy (dstp, path - nl, nl); 905 memcpy (dstp, path - nl, nl);
321 dstp += nl; 906 dstp += nl;
322 } 907 }
908 else
909 *dstp++ = *path++;
910#if 0
323 else 911 else
324 { 912 {
325 STRLEN cl; 913 STRLEN cl;
326 to_utf8_fold (path, dstp, &cl); 914 to_utf8_fold (path, dstp, &cl);
327 dstp += cl; 915 dstp += cl;
328 path += is_utf8_char (path); 916 path += is_utf8_char (path);
329 } 917 }
918#endif
330 } 919 }
331 920
332 RETVAL = newSVpvn ((const char *)dst, dstp - dst); 921 RETVAL = newSVpvn ((const char *)dst, dstp - dst);
333} 922}
334 OUTPUT: 923 OUTPUT:
359 } 948 }
360} 949}
361 OUTPUT: 950 OUTPUT:
362 RETVAL 951 RETVAL
363 952
953#############################################################################
954
955MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
956
957void
958dump_ascii85 (PerlIO *fp, GdkPixbuf *pb)
959 CODE:
960{
961 int w = gdk_pixbuf_get_width (pb);
962 int h = gdk_pixbuf_get_height (pb);
963 int x, y, i;
964 guchar *dst;
965 int bpp = gdk_pixbuf_get_n_channels (pb);
966 guchar *src = gdk_pixbuf_get_pixels (pb);
967 int sstr = gdk_pixbuf_get_rowstride (pb);
968
969 a85_init ();
970
971 for (y = 0; y < h; y++)
972 for (x = 0; x < w; x++)
973 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
974 a85_push (fp, src [x * bpp + y * sstr + i]);
975
976 a85_finish (fp);
977}
978
979void
980dump_binary (PerlIO *fp, GdkPixbuf *pb)
981 CODE:
982{
983 int w = gdk_pixbuf_get_width (pb);
984 int h = gdk_pixbuf_get_height (pb);
985 int x, y, i;
986 guchar *dst;
987 int bpp = gdk_pixbuf_get_n_channels (pb);
988 guchar *src = gdk_pixbuf_get_pixels (pb);
989 int sstr = gdk_pixbuf_get_rowstride (pb);
990
991 for (y = 0; y < h; y++)
992 for (x = 0; x < w; x++)
993 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
994 PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
995}
996
997#############################################################################
998
999MODULE = Gtk2::CV PACKAGE = Gtk2::CV
1000
364SV * 1001SV *
365pb_to_p7 (GdkPixbuf *pb) 1002pb_to_hv84 (GdkPixbuf *pb)
366 CODE: 1003 CODE:
367{ 1004{
368 int w = gdk_pixbuf_get_width (pb); 1005 int w = gdk_pixbuf_get_width (pb);
369 int h = gdk_pixbuf_get_height (pb); 1006 int h = gdk_pixbuf_get_height (pb);
370 int x, y; 1007 int x, y;
371 guchar *dst; 1008 guchar *dst;
372 int bpp = gdk_pixbuf_get_n_channels (pb); 1009 int bpp = gdk_pixbuf_get_n_channels (pb);
373 guchar *src = gdk_pixbuf_get_pixels (pb); 1010 guchar *src = gdk_pixbuf_get_pixels (pb);
374 int sstr = gdk_pixbuf_get_rowstride (pb); 1011 int sstr = gdk_pixbuf_get_rowstride (pb);
375 int Er[IW], Eg[IW], Eb[IW];
376 int seed = 77;
377
378 RETVAL = newSV (w * h);
379 SvPOK_only (RETVAL);
380 SvCUR_set (RETVAL, w * h);
381
382 dst = (guchar *)SvPVX (RETVAL);
383
384 memset (Er, 0, sizeof (int) * IW);
385 memset (Eg, 0, sizeof (int) * IW);
386 memset (Eb, 0, sizeof (int) * IW);
387
388 /* some primitive error distribution + random dithering */
389
390 for (y = 0; y < h; y++)
391 {
392 int er = 0, eg = 0, eb = 0;
393
394 for (x = 0; x < w; x++)
395 {
396 int r, g, b;
397 guchar *p = src + x * bpp + y * sstr;
398
399 r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255;
400 g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255;
401 b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255;
402
403 r = r > 7 ? 7 : r < 0 ? 0 : r;
404 g = g > 7 ? 7 : g < 0 ? 0 : g;
405 b = b > 3 ? 3 : b < 0 ? 0 : b;
406
407 er += p[0] - (r * 255 + 4) / 7;
408 eg += p[1] - (g * 255 + 4) / 7;
409 eb += p[2] - (b * 255 + 2) / 3;
410
411 Er[x] = er >> 1; er -= (er + 1) >> 1;
412 Eg[x] = eg >> 1; eg -= (eg + 1) >> 1;
413 Eb[x] = eb >> 1; eb -= (eb + 1) >> 1;
414
415 *dst++ = r << 5 | g << 2 | b;
416 }
417 }
418}
419 OUTPUT:
420 RETVAL
421
422#############################################################################
423
424MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
425
426void
427dump_ascii85 (PerlIO *fp, GdkPixbuf *pb)
428 CODE:
429{
430 int w = gdk_pixbuf_get_width (pb);
431 int h = gdk_pixbuf_get_height (pb);
432 int x, y, i;
433 guchar *dst;
434 int bpp = gdk_pixbuf_get_n_channels (pb);
435 guchar *src = gdk_pixbuf_get_pixels (pb);
436 int sstr = gdk_pixbuf_get_rowstride (pb);
437
438 a85_init ();
439
440 for (y = 0; y < h; y++)
441 for (x = 0; x < w; x++)
442 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
443 a85_push (fp, src [x * bpp + y * sstr + i]);
444
445 a85_finish (fp);
446}
447
448void
449dump_binary (PerlIO *fp, GdkPixbuf *pb)
450 CODE:
451{
452 int w = gdk_pixbuf_get_width (pb);
453 int h = gdk_pixbuf_get_height (pb);
454 int x, y, i;
455 guchar *dst;
456 int bpp = gdk_pixbuf_get_n_channels (pb);
457 guchar *src = gdk_pixbuf_get_pixels (pb);
458 int sstr = gdk_pixbuf_get_rowstride (pb);
459
460 for (y = 0; y < h; y++)
461 for (x = 0; x < w; x++)
462 for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
463 PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
464}
465
466#############################################################################
467
468MODULE = Gtk2::CV PACKAGE = Gtk2::CV
469
470SV *
471pb_to_hv84 (GdkPixbuf *pb)
472 CODE:
473{
474 int w = gdk_pixbuf_get_width (pb);
475 int h = gdk_pixbuf_get_height (pb);
476 int x, y;
477 guchar *dst;
478 int bpp = gdk_pixbuf_get_n_channels (pb);
479 guchar *src = gdk_pixbuf_get_pixels (pb);
480 int sstr = gdk_pixbuf_get_rowstride (pb);
481 1012
482 RETVAL = newSV (6 * 8 * 12 / 8); 1013 RETVAL = newSV (6 * 8 * 12 / 8);
483 SvPOK_only (RETVAL); 1014 SvPOK_only (RETVAL);
484 SvCUR_set (RETVAL, 6 * 8 * 12 / 8); 1015 SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
485 1016
554############################################################################# 1085#############################################################################
555 1086
556MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster 1087MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
557 1088
558SV * 1089SV *
559make_histograms (SV *ar) 1090extract_features (SV *ar)
560 CODE: 1091 CODE:
561{ 1092{
562 int i; 1093 int i;
563 AV *av, *result; 1094 AV *av, *result;
564 1095
565 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) 1096 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
566 croak ("Not an array ref as first argument to make_histogram"); 1097 croak ("Not an array ref as first argument to extract_features");
567 1098
568 av = (AV *) SvRV (ar); 1099 av = (AV *) SvRV (ar);
569 result = newAV (); 1100 result = newAV ();
570 1101
571 for (i = 0; i <= av_len (av); ++i) 1102 for (i = 0; i <= av_len (av); ++i)
572 { 1103 {
573 const int HISTSIZE = 64;
574
575 int j;
576 SV *sv = *av_fetch (av, i, 1); 1104 SV *sv = *av_fetch (av, i, 1);
577 STRLEN len;
578 char *buf = SvPVbyte (sv, len);
579
580 int tmphist[HISTSIZE];
581 float *hist;
582
583 SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); 1105 SV *histsv = newSV (9 * sizeof (float) + 1);
1106
584 SvPOK_on (histsv); 1107 SvPOK_on (histsv);
585 SvCUR_set (histsv, HISTSIZE * sizeof (float)); 1108 SvCUR_set (histsv, 9 * sizeof (float));
586 hist = (float *)SvPVX (histsv); 1109 float *hist = (float *)SvPVX (histsv);
587 1110
588 Zero (tmphist, sizeof (tmphist), char); 1111 struct feature f_h, f_s, f_v;
1112 feature_init (&f_h);
1113 feature_init (&f_s);
1114 feature_init (&f_v);
589 1115
590 for (j = len; j--; )
591 {
592 unsigned int idx
593 = ((*buf & 0xc0) >> 2)
594 | ((*buf & 0x18) >> 1)
595 | (*buf & 0x03);
596
597 ++tmphist[idx];
598 ++buf;
599 }
600 1116 {
601 for (j = 0; j < HISTSIZE; ++j) 1117 STRLEN len;
602 hist[j] = (float)tmphist[j] / (len + 1e-30); 1118 unsigned char *buf = (unsigned char *)SvPVbyte (sv, len);
1119 while (len >= 3)
1120 {
1121 unsigned int r, g, b, h, s, v;
1122 r = *buf++; g = *buf++; b = *buf++;
1123 rgb_to_hsv (r, g, b, &h, &s, &v);
1124
1125 feature_update_pass_1 (&f_h, h);
1126 feature_update_pass_1 (&f_s, s);
1127 feature_update_pass_1 (&f_v, v);
1128
1129 len -= 3;
1130 }
1131
1132 feature_finish_pass_1 (&f_h);
1133 feature_finish_pass_1 (&f_s);
1134 feature_finish_pass_1 (&f_v);
1135 }
1136
1137 {
1138 STRLEN len;
1139 unsigned char *buf = (unsigned char *)SvPVbyte (sv, len);
1140 while (len >= 3)
1141 {
1142 unsigned int r, g, b, h, s, v;
1143 r = *buf++; g = *buf++; b = *buf++;
1144 rgb_to_hsv (r, g, b, &h, &s, &v);
1145
1146 feature_update_pass_2 (&f_h, h);
1147 feature_update_pass_2 (&f_s, s);
1148 feature_update_pass_2 (&f_v, v);
1149
1150 len -= 3;
1151 }
1152
1153 feature_finish_pass_2 (&f_h);
1154 feature_finish_pass_2 (&f_s);
1155 feature_finish_pass_2 (&f_v);
1156 }
1157
1158 hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
1159 hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
1160 hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
603 1161
604 av_push (result, histsv); 1162 av_push (result, histsv);
605 } 1163 }
606 1164
607 RETVAL = newRV_noinc ((SV *)result); 1165 RETVAL = newRV_noinc ((SV *)result);
608} 1166}
609 OUTPUT: 1167 OUTPUT:
610 RETVAL 1168 RETVAL
611 1169
612

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines