… | |
… | |
7 | #include <math.h> |
7 | #include <math.h> |
8 | |
8 | |
9 | #include <magic.h> |
9 | #include <magic.h> |
10 | |
10 | |
11 | #include <jpeglib.h> |
11 | #include <jpeglib.h> |
|
|
12 | #include <jerror.h> |
|
|
13 | |
12 | #include <glib.h> |
14 | #include <glib.h> |
13 | #include <gtk/gtk.h> |
15 | #include <gtk/gtk.h> |
|
|
16 | #include <gdk/gdkx.h> |
14 | #include <gdk-pixbuf/gdk-pixbuf.h> |
17 | #include <gdk-pixbuf/gdk-pixbuf.h> |
15 | |
18 | |
16 | #include <gperl.h> |
19 | #include <gperl.h> |
17 | #include <gtk2perl.h> |
20 | #include <gtk2perl.h> |
18 | |
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 | |
19 | #define IW 80 /* MUST match Schnauzer.pm! */ |
36 | #define IW 80 /* MUST match Schnauzer.pm! */ |
20 | #define IH 60 /* MUST match Schnauzer.pm! */ |
37 | #define IH 60 /* MUST match Schnauzer.pm! */ |
21 | |
38 | |
22 | #define RAND (seed = (seed + 7141) * 54773 % 134456) |
|
|
23 | |
|
|
24 | #define LINELENGTH 240 |
39 | #define LINELENGTH 240 |
25 | |
40 | |
26 | #define ELLIPSIS "\xe2\x80\xa6" |
41 | #define ELLIPSIS "\xe2\x80\xa6" |
|
|
42 | |
|
|
43 | typedef char *octet_string; |
|
|
44 | |
|
|
45 | static magic_t magic_cookie[2]; /* !mime, mime */ |
27 | |
46 | |
28 | struct jpg_err_mgr |
47 | struct jpg_err_mgr |
29 | { |
48 | { |
30 | struct jpeg_error_mgr err; |
49 | struct jpeg_error_mgr err; |
31 | jmp_buf setjmp_buffer; |
50 | jmp_buf setjmp_buffer; |
32 | }; |
51 | }; |
33 | |
52 | |
34 | static void |
53 | static void |
35 | cv_error_exit (j_common_ptr cinfo) |
54 | cv_error_exit (j_common_ptr cinfo) |
36 | { |
55 | { |
|
|
56 | cinfo->err->output_message (cinfo); |
37 | longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99); |
57 | longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99); |
38 | } |
58 | } |
39 | |
59 | |
40 | static void |
60 | static void |
41 | cv_error_output (j_common_ptr cinfo) |
61 | cv_error_output (j_common_ptr cinfo) |
42 | { |
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); |
43 | return; |
68 | return; |
44 | } |
69 | } |
45 | |
70 | |
46 | static void |
71 | static void |
47 | rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, |
72 | rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b, |
… | |
… | |
144 | if (!--a85_cnt) |
169 | if (!--a85_cnt) |
145 | { |
170 | { |
146 | a85_cnt = 4; |
171 | a85_cnt = 4; |
147 | if (a85_val) |
172 | if (a85_val) |
148 | { |
173 | { |
149 | a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
174 | a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
150 | a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
175 | a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
151 | a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
176 | a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
152 | a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
177 | a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
153 | a85_ptr[0] = (a85_val ) + 33; |
178 | a85_ptr[0] = (a85_val ) + 33; |
154 | |
179 | |
155 | a85_ptr += 5; |
180 | a85_ptr += 5; |
156 | } |
181 | } |
… | |
… | |
179 | |
204 | |
180 | PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
205 | PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
181 | } |
206 | } |
182 | |
207 | |
183 | ///////////////////////////////////////////////////////////////////////////// |
208 | ///////////////////////////////////////////////////////////////////////////// |
|
|
209 | // memory source for libjpeg |
|
|
210 | |
|
|
211 | static void cv_ms_init (j_decompress_ptr cinfo) |
|
|
212 | { |
|
|
213 | } |
|
|
214 | |
|
|
215 | static void cv_ms_term (j_decompress_ptr cinfo) |
|
|
216 | { |
|
|
217 | } |
|
|
218 | |
|
|
219 | static 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 | |
|
|
235 | static 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 | |
|
|
243 | static 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 | |
|
|
267 | struct bmff_box |
|
|
268 | { |
|
|
269 | char type[4]; |
|
|
270 | const uint8_t *ptr; |
|
|
271 | size_t size; |
|
|
272 | }; |
|
|
273 | |
|
|
274 | static int |
|
|
275 | bmff_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 | ///////////////////////////////////////////////////////////////////////////// |
184 | |
303 | |
185 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
304 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
186 | |
305 | |
187 | PROTOTYPES: ENABLE |
306 | PROTOTYPES: ENABLE |
188 | |
307 | |
|
|
308 | void |
|
|
309 | _exit (int code) |
|
|
310 | |
|
|
311 | # calculate the common prefix length of two strings |
189 | # missing function in perl. really :) |
312 | # missing function in perl. really :) |
190 | int |
313 | int |
191 | common_prefix_length (a, b) |
314 | common_prefix_length (a, b) |
192 | unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
315 | unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
193 | unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
316 | unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
… | |
… | |
201 | } |
324 | } |
202 | |
325 | |
203 | OUTPUT: |
326 | OUTPUT: |
204 | RETVAL |
327 | RETVAL |
205 | |
328 | |
|
|
329 | int |
|
|
330 | common_prefix_length_byte (a, b) |
|
|
331 | unsigned char *a = (unsigned char *)SvPVbyte_nolen ($arg); |
|
|
332 | unsigned char *b = (unsigned char *)SvPVbyte_nolen ($arg); |
|
|
333 | CODE: |
|
|
334 | RETVAL = 0; |
|
|
335 | while (*a == *b && *a) |
|
|
336 | a++, b++, RETVAL++; |
|
|
337 | OUTPUT: |
|
|
338 | RETVAL |
|
|
339 | |
206 | const char * |
340 | const char * |
207 | magic (const char *path) |
341 | magic (SV *path_or_data) |
|
|
342 | ALIAS: |
|
|
343 | magic = 0 |
|
|
344 | magic_mime = 1 |
|
|
345 | magic_buffer = 2 |
|
|
346 | magic_buffer_mime = 3 |
208 | CODE: |
347 | CODE: |
209 | { |
348 | { |
210 | static magic_t cookie; |
349 | STRLEN len; |
|
|
350 | char *data = SvPVbyte (path_or_data, len); |
211 | |
351 | |
212 | if (!cookie) |
352 | if (!magic_cookie[0]) |
|
|
353 | { |
|
|
354 | magic_cookie[0] = magic_open (MAGIC_SYMLINK); |
|
|
355 | magic_cookie[1] = magic_open (MAGIC_SYMLINK | MAGIC_MIME_TYPE); |
|
|
356 | magic_load (magic_cookie[0], 0); |
|
|
357 | magic_load (magic_cookie[1], 0); |
213 | { |
358 | } |
214 | cookie = magic_open (MAGIC_NONE); |
|
|
215 | |
359 | |
216 | if (cookie) |
360 | perlinterp_release (); |
217 | magic_load (cookie, 0); |
|
|
218 | else |
|
|
219 | XSRETURN_UNDEF; |
|
|
220 | } |
|
|
221 | |
361 | |
222 | RETVAL = magic_file (cookie, path); |
362 | RETVAL = ix & 2 |
|
|
363 | ? magic_buffer (magic_cookie[ix & 1], data, len) |
|
|
364 | : magic_file (magic_cookie[ix & 1], data); |
|
|
365 | |
|
|
366 | perlinterp_acquire (); |
223 | } |
367 | } |
224 | OUTPUT: |
368 | OUTPUT: |
225 | RETVAL |
369 | RETVAL |
226 | |
370 | |
227 | const char * |
|
|
228 | magic_mime (const char *path) |
|
|
229 | CODE: |
|
|
230 | { |
|
|
231 | static magic_t cookie; |
|
|
232 | |
|
|
233 | if (!cookie) |
|
|
234 | { |
|
|
235 | cookie = magic_open (MAGIC_MIME); |
|
|
236 | |
|
|
237 | if (cookie) |
|
|
238 | magic_load (cookie, 0); |
|
|
239 | else |
|
|
240 | XSRETURN_UNDEF; |
|
|
241 | } |
|
|
242 | |
|
|
243 | RETVAL = magic_file (cookie, path); |
|
|
244 | } |
|
|
245 | OUTPUT: |
|
|
246 | RETVAL |
|
|
247 | |
|
|
248 | # missing in Gtk2 perl module |
371 | # missing/broken in Gtk2 perl module |
|
|
372 | |
|
|
373 | void |
|
|
374 | gdk_window_clear_hints (GdkWindow *window) |
|
|
375 | CODE: |
|
|
376 | gdk_window_set_geometry_hints (window, 0, 0); |
249 | |
377 | |
250 | gboolean |
378 | gboolean |
251 | gdk_net_wm_supports (GdkAtom property) |
379 | gdk_net_wm_supports (GdkAtom property) |
252 | CODE: |
380 | CODE: |
253 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
381 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
… | |
… | |
259 | RETVAL |
387 | RETVAL |
260 | |
388 | |
261 | GdkPixbuf_noinc * |
389 | GdkPixbuf_noinc * |
262 | dealpha_expose (GdkPixbuf *pb) |
390 | dealpha_expose (GdkPixbuf *pb) |
263 | CODE: |
391 | CODE: |
|
|
392 | perlinterp_release (); |
264 | { |
393 | { |
265 | int w = gdk_pixbuf_get_width (pb); |
394 | int w = gdk_pixbuf_get_width (pb); |
266 | int h = gdk_pixbuf_get_height (pb); |
395 | int h = gdk_pixbuf_get_height (pb); |
267 | int bpp = gdk_pixbuf_get_n_channels (pb); |
396 | int bpp = gdk_pixbuf_get_n_channels (pb); |
268 | int x, y, i; |
397 | int x, y, i; |
… | |
… | |
277 | for (x = 0; x < w; x++) |
406 | for (x = 0; x < w; x++) |
278 | for (y = 0; y < h; y++) |
407 | for (y = 0; y < h; y++) |
279 | for (i = 0; i < 3; i++) |
408 | for (i = 0; i < 3; i++) |
280 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
409 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
281 | } |
410 | } |
|
|
411 | perlinterp_acquire (); |
282 | OUTPUT: |
412 | OUTPUT: |
283 | RETVAL |
413 | RETVAL |
284 | |
414 | |
285 | GdkPixbuf_noinc * |
415 | GdkPixbuf_noinc * |
286 | rotate (GdkPixbuf *pb, int angle) |
416 | rotate (GdkPixbuf *pb, int angle) |
287 | CODE: |
417 | CODE: |
|
|
418 | perlinterp_release (); |
|
|
419 | if (angle < 0) |
|
|
420 | angle += 360; |
288 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
421 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
289 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
422 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
290 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
423 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
291 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
424 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
292 | : angle); |
425 | : angle); |
|
|
426 | perlinterp_acquire (); |
293 | OUTPUT: |
427 | OUTPUT: |
294 | RETVAL |
428 | RETVAL |
295 | |
429 | |
|
|
430 | const char * |
|
|
431 | filetype (SV *image_data) |
|
|
432 | CODE: |
|
|
433 | { |
|
|
434 | STRLEN data_len; |
|
|
435 | U8 *data = SvPVbyte (image_data, data_len); |
|
|
436 | static const unsigned char jxl_header[] = { |
|
|
437 | 0, 0, 0, 0x0c, 0x4a, 0x58, 0x4c, 0x20, 0x0d, 0xa, 0x87, 0x0a |
|
|
438 | }; |
|
|
439 | |
|
|
440 | if (data_len >= 20 |
|
|
441 | && data[0] == 0xff |
|
|
442 | && data[1] == 0xd8 |
|
|
443 | && data[2] == 0xff) |
|
|
444 | RETVAL = "image/jpeg"; |
|
|
445 | else if (data_len >= 12 |
|
|
446 | && data[ 0] == (U8)'R' |
|
|
447 | && data[ 1] == (U8)'I' |
|
|
448 | && data[ 2] == (U8)'F' |
|
|
449 | && data[ 3] == (U8)'F' |
|
|
450 | && data[ 8] == (U8)'W' |
|
|
451 | && data[ 9] == (U8)'E' |
|
|
452 | && data[10] == (U8)'B' |
|
|
453 | && data[11] == (U8)'P') |
|
|
454 | RETVAL = "image/webp"; |
|
|
455 | else if (data_len >= 16 |
|
|
456 | && data[ 0] == 0x89 |
|
|
457 | && data[ 1] == (U8)'P' |
|
|
458 | && data[ 2] == (U8)'N' |
|
|
459 | && data[ 3] == (U8)'G' |
|
|
460 | && data[ 4] == 0x0d |
|
|
461 | && data[ 5] == 0x0a |
|
|
462 | && data[ 6] == 0x1a |
|
|
463 | && data[ 7] == 0x0a) |
|
|
464 | RETVAL = "image/png"; |
|
|
465 | else if (data_len >= sizeof (jxl_header) && memcmp (data, jxl_header, sizeof (jxl_header)) == 0) |
|
|
466 | RETVAL = "image/jxl"; // todo: might want to use JxlSignatureCheck |
|
|
467 | else if (data_len >= 2 |
|
|
468 | && data[0] == 0xff |
|
|
469 | && data[1] == 0x0a) |
|
|
470 | RETVAL = "image/jxl"; |
|
|
471 | else if (data_len >= 13 |
|
|
472 | && data[0] == 'G' |
|
|
473 | && data[1] == 'I' |
|
|
474 | && data[2] == 'F' |
|
|
475 | && data[3] == '8' |
|
|
476 | //&& (data[4] == '7' || data[4] == '9') |
|
|
477 | && data[5] == 'a') |
|
|
478 | { |
|
|
479 | RETVAL = "image/gif"; |
|
|
480 | |
|
|
481 | // now see if its animated - we require the netscape application header for this |
|
|
482 | int ofs = 13; |
|
|
483 | |
|
|
484 | if (data[10] & 0x80) |
|
|
485 | ofs += (1 << ((data[10] & 7) + 1)) * 3; |
|
|
486 | |
|
|
487 | if (data_len >= ofs + 2 + 1 + 11) |
|
|
488 | { |
|
|
489 | |
|
|
490 | // skip a graphic control extension block. we assume |
|
|
491 | // there is at most one such block - while the NAB |
|
|
492 | // has to come firstz, some files do not obey this |
|
|
493 | if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9) |
|
|
494 | ofs += 3 + data[ofs + 2] + 1; |
|
|
495 | |
|
|
496 | if (data_len >= ofs + 2 + 1 + 11) |
|
|
497 | if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1)) |
|
|
498 | RETVAL = "video/gif"; |
|
|
499 | } |
|
|
500 | } |
|
|
501 | else if (data_len >= 0x8000 + 6 |
|
|
502 | && data[0x8000+1] == (U8)'B' |
|
|
503 | && data[0x8000+2] == (U8)'E' |
|
|
504 | && data[0x8000+3] == (U8)'A' |
|
|
505 | && data[0x8000+4] == (U8)'0' |
|
|
506 | && data[0x8000+5] == (U8)'1') |
|
|
507 | RETVAL = "video/iso-bluray"; |
|
|
508 | else if (data_len >= 0x8000 + 6 |
|
|
509 | && data[0x8000+1] == (U8)'C' |
|
|
510 | && data[0x8000+2] == (U8)'D' |
|
|
511 | && data[0x8000+3] == (U8)'0' |
|
|
512 | && data[0x8000+4] == (U8)'0' |
|
|
513 | && data[0x8000+5] == (U8)'1') |
|
|
514 | RETVAL = "video/iso9660"; |
|
|
515 | |
|
|
516 | else |
|
|
517 | XSRETURN_UNDEF; |
|
|
518 | } |
|
|
519 | OUTPUT: |
|
|
520 | RETVAL |
|
|
521 | |
296 | GdkPixbuf_noinc * |
522 | GdkPixbuf_noinc * |
297 | load_jpeg (SV *path, int thumbnail=0) |
523 | decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
|
|
524 | CODE: |
|
|
525 | { |
|
|
526 | #if WEBP |
|
|
527 | STRLEN data_size; |
|
|
528 | int alpha; |
|
|
529 | WebPData data; |
|
|
530 | WebPDemuxer *demux; |
|
|
531 | WebPIterator iter; |
|
|
532 | WebPDecoderConfig config; |
|
|
533 | int inw, inh; |
|
|
534 | |
|
|
535 | data.bytes = (uint8_t *)SvPVbyte (image_data, data_size); |
|
|
536 | data.size = data_size; |
|
|
537 | |
|
|
538 | perlinterp_release (); |
|
|
539 | |
|
|
540 | RETVAL = 0; |
|
|
541 | |
|
|
542 | if (!(demux = WebPDemux (&data))) |
|
|
543 | goto err_demux; |
|
|
544 | |
|
|
545 | if (!WebPDemuxGetFrame (demux, 1, &iter)) |
|
|
546 | goto err_iter; |
|
|
547 | |
|
|
548 | if (!WebPInitDecoderConfig (&config)) |
|
|
549 | goto err_iter; |
|
|
550 | |
|
|
551 | config.options.use_threads = 1; |
|
|
552 | |
|
|
553 | if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK) |
|
|
554 | goto err_iter; |
|
|
555 | |
|
|
556 | inw = config.input.width; |
|
|
557 | inh = config.input.height; |
|
|
558 | |
|
|
559 | if (thumbnail) |
|
|
560 | { |
|
|
561 | if (inw * ih > inh * iw) |
|
|
562 | ih = (iw * inh + inw - 1) / inw; |
|
|
563 | else |
|
|
564 | iw = (ih * inw + inh - 1) / inh; |
|
|
565 | |
|
|
566 | config.options.bypass_filtering = 1; |
|
|
567 | config.options.no_fancy_upsampling = 1; |
|
|
568 | |
|
|
569 | config.options.use_scaling = 1; |
|
|
570 | config.options.scaled_width = iw; |
|
|
571 | config.options.scaled_height = ih; |
|
|
572 | } |
|
|
573 | else |
|
|
574 | { |
|
|
575 | iw = inw; |
|
|
576 | ih = inh; |
|
|
577 | } |
|
|
578 | |
|
|
579 | alpha = !!config.input.has_alpha; |
|
|
580 | |
|
|
581 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih); |
|
|
582 | if (!RETVAL) |
|
|
583 | goto err_iter; |
|
|
584 | |
|
|
585 | config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB; |
|
|
586 | config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL); |
|
|
587 | config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL); |
|
|
588 | config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL); |
|
|
589 | config.output.is_external_memory = 1; |
|
|
590 | |
|
|
591 | if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK) |
|
|
592 | { |
|
|
593 | g_object_unref (RETVAL); |
|
|
594 | RETVAL = 0; |
|
|
595 | goto err_iter; |
|
|
596 | } |
|
|
597 | |
|
|
598 | err_iter: |
|
|
599 | WebPDemuxReleaseIterator (&iter); |
|
|
600 | err_demux: |
|
|
601 | WebPDemuxDelete (demux); |
|
|
602 | |
|
|
603 | perlinterp_acquire (); |
|
|
604 | #else |
|
|
605 | croak ("load_webp: webp not enabled at compile time"); |
|
|
606 | #endif |
|
|
607 | } |
|
|
608 | OUTPUT: |
|
|
609 | RETVAL |
|
|
610 | |
|
|
611 | GdkPixbuf_noinc * |
|
|
612 | decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
|
|
613 | CODE: |
|
|
614 | { |
|
|
615 | #if JXL |
|
|
616 | JxlDecoder *dec; |
|
|
617 | JxlBasicInfo info; |
|
|
618 | const uint8_t *next_in = (uint8_t *)SvPVbyte_nolen (image_data); |
|
|
619 | size_t avail_in = SvCUR (image_data); |
|
|
620 | const char *error = 0; |
|
|
621 | JxlDecoderStatus status; |
|
|
622 | static void *runner_cache; |
|
|
623 | void *runner = 0; |
|
|
624 | |
|
|
625 | RETVAL = 0; |
|
|
626 | |
|
|
627 | if (runner_cache) |
|
|
628 | runner = runner_cache; |
|
|
629 | else |
|
|
630 | runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ()); |
|
|
631 | |
|
|
632 | runner_cache = 0; |
|
|
633 | |
|
|
634 | perlinterp_release (); |
|
|
635 | |
|
|
636 | dec = JxlDecoderCreate (0); |
|
|
637 | |
|
|
638 | error = "JxlDecoderCreate failed"; |
|
|
639 | if (!dec) |
|
|
640 | goto done; |
|
|
641 | |
|
|
642 | status = JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner); |
|
|
643 | error = "JxlDecoderSetParallelRunner failed"; |
|
|
644 | if (status != JXL_DEC_SUCCESS) |
|
|
645 | goto done; |
|
|
646 | |
|
|
647 | error = "JxlDecoderSubscribeEvents failed"; |
|
|
648 | status = JxlDecoderSubscribeEvents (dec, JXL_DEC_BASIC_INFO | JXL_DEC_FULL_IMAGE); |
|
|
649 | if (status != JXL_DEC_SUCCESS) |
|
|
650 | goto done; |
|
|
651 | |
|
|
652 | status = JxlDecoderSetInput (dec, next_in, avail_in); |
|
|
653 | error = "JxlDecoderSetInput failed"; |
|
|
654 | if (status != JXL_DEC_SUCCESS) |
|
|
655 | goto done; |
|
|
656 | |
|
|
657 | for (;;) |
|
|
658 | { |
|
|
659 | status = JxlDecoderProcessInput (dec); |
|
|
660 | |
|
|
661 | switch (status) |
|
|
662 | { |
|
|
663 | case JXL_DEC_FULL_IMAGE: |
|
|
664 | error = 0; |
|
|
665 | goto done; |
|
|
666 | |
|
|
667 | case JXL_DEC_ERROR: |
|
|
668 | error = "JxlDecoderProcessInput failed"; |
|
|
669 | goto done; |
|
|
670 | |
|
|
671 | case JXL_DEC_NEED_MORE_INPUT: |
|
|
672 | error = "incomplete file"; |
|
|
673 | goto done; |
|
|
674 | |
|
|
675 | case JXL_DEC_SUCCESS: |
|
|
676 | error = "incomplete decode"; |
|
|
677 | goto done; |
|
|
678 | |
|
|
679 | case JXL_DEC_BASIC_INFO: |
|
|
680 | { |
|
|
681 | status = JxlDecoderGetBasicInfo (dec, &info); |
|
|
682 | error = "JxlDecoderGetBasicInfo failed"; |
|
|
683 | if (status != JXL_DEC_SUCCESS) |
|
|
684 | goto done; |
|
|
685 | |
|
|
686 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize); |
|
|
687 | error = "unable to allocate pixbuf"; |
|
|
688 | if (!RETVAL) |
|
|
689 | goto done; |
|
|
690 | |
|
|
691 | JxlPixelFormat format = { |
|
|
692 | info.alpha_bits ? 4 : 3, |
|
|
693 | JXL_TYPE_UINT8, |
|
|
694 | JXL_NATIVE_ENDIAN, |
|
|
695 | gdk_pixbuf_get_rowstride (RETVAL) |
|
|
696 | }; |
|
|
697 | |
|
|
698 | // cannot use gdk_pixbuf_get_byte_length because that does |
|
|
699 | // not return the size of the buffer, but the size of the buffer without |
|
|
700 | // the last padding bytes. the internal buffer is rowstride * ysize, |
|
|
701 | // and this is what the jxl decoder needs. none of this is documented |
|
|
702 | // in either library, of course. |
|
|
703 | |
|
|
704 | status = JxlDecoderSetImageOutBuffer ( |
|
|
705 | dec, |
|
|
706 | &format, |
|
|
707 | gdk_pixbuf_get_pixels (RETVAL), |
|
|
708 | gdk_pixbuf_get_rowstride (RETVAL) * info.ysize |
|
|
709 | ); |
|
|
710 | error = "JxlDecoderSetImageOutBuffer failed"; |
|
|
711 | if (status != JXL_DEC_SUCCESS) |
|
|
712 | goto done; |
|
|
713 | } |
|
|
714 | break; |
|
|
715 | |
|
|
716 | default: |
|
|
717 | error = "unexpected event"; |
|
|
718 | goto done; |
|
|
719 | } |
|
|
720 | } |
|
|
721 | |
|
|
722 | done: |
|
|
723 | if (dec) |
|
|
724 | JxlDecoderDestroy (dec); |
|
|
725 | |
|
|
726 | perlinterp_acquire (); |
|
|
727 | |
|
|
728 | if (runner_cache) |
|
|
729 | JxlThreadParallelRunnerDestroy (runner); |
|
|
730 | |
|
|
731 | runner_cache = runner; |
|
|
732 | |
|
|
733 | if (error) |
|
|
734 | { |
|
|
735 | if (RETVAL) |
|
|
736 | g_object_unref (RETVAL); |
|
|
737 | |
|
|
738 | croak ("load_jxl: %s (status %d)", error, status); |
|
|
739 | } |
|
|
740 | #else |
|
|
741 | croak ("load_jxl: jpeg-xl not enabled at compile time"); |
|
|
742 | #endif |
|
|
743 | } |
|
|
744 | OUTPUT: |
|
|
745 | RETVAL |
|
|
746 | |
|
|
747 | GdkPixbuf_noinc * |
|
|
748 | decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
298 | CODE: |
749 | CODE: |
299 | { |
750 | { |
300 | struct jpeg_decompress_struct cinfo; |
751 | struct jpeg_decompress_struct cinfo; |
301 | struct jpg_err_mgr jerr; |
752 | struct jpg_err_mgr jerr; |
302 | guchar *data; |
|
|
303 | int rs; |
753 | int rs; |
304 | FILE *fp; |
|
|
305 | volatile GdkPixbuf *pb = 0; |
754 | volatile GdkPixbuf *pb = 0; |
|
|
755 | STRLEN data_len; |
|
|
756 | guchar *data = SvPVbyte (image_data, data_len); |
306 | |
757 | |
307 | RETVAL = 0; |
758 | RETVAL = 0; |
308 | |
759 | |
309 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
760 | perlinterp_release (); |
310 | |
|
|
311 | if (!fp) |
|
|
312 | XSRETURN_UNDEF; |
|
|
313 | |
761 | |
314 | cinfo.err = jpeg_std_error (&jerr.err); |
762 | cinfo.err = jpeg_std_error (&jerr.err); |
315 | |
763 | |
316 | jerr.err.error_exit = cv_error_exit; |
764 | jerr.err.error_exit = cv_error_exit; |
317 | jerr.err.output_message = cv_error_output; |
765 | jerr.err.output_message = cv_error_output; |
318 | |
766 | |
319 | if ((rs = setjmp (jerr.setjmp_buffer))) |
767 | if ((rs = setjmp (jerr.setjmp_buffer))) |
320 | { |
768 | { |
321 | fclose (fp); |
|
|
322 | jpeg_destroy_decompress (&cinfo); |
769 | jpeg_destroy_decompress (&cinfo); |
323 | |
770 | |
324 | if (pb) |
771 | if (pb) |
325 | g_object_unref ((gpointer)pb); |
772 | g_object_unref ((gpointer)pb); |
326 | |
773 | |
|
|
774 | perlinterp_acquire (); |
327 | XSRETURN_UNDEF; |
775 | XSRETURN_UNDEF; |
328 | } |
776 | } |
329 | |
777 | |
|
|
778 | if (!data_len) |
|
|
779 | longjmp (jerr.setjmp_buffer, 4); |
|
|
780 | |
330 | jpeg_create_decompress (&cinfo); |
781 | jpeg_create_decompress (&cinfo); |
|
|
782 | cv_jpeg_mem_src (&cinfo, data, data_len); |
331 | |
783 | |
332 | jpeg_stdio_src (&cinfo, fp); |
|
|
333 | jpeg_read_header (&cinfo, TRUE); |
784 | jpeg_read_header (&cinfo, TRUE); |
334 | |
785 | |
335 | cinfo.dct_method = JDCT_DEFAULT; |
786 | cinfo.dct_method = JDCT_DEFAULT; |
336 | cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ |
787 | cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ |
337 | cinfo.do_block_smoothing = FALSE; |
788 | cinfo.do_block_smoothing = FALSE; |
… | |
… | |
347 | { |
798 | { |
348 | cinfo.dct_method = JDCT_FASTEST; |
799 | cinfo.dct_method = JDCT_FASTEST; |
349 | cinfo.do_fancy_upsampling = FALSE; |
800 | cinfo.do_fancy_upsampling = FALSE; |
350 | |
801 | |
351 | while (cinfo.scale_denom < 8 |
802 | while (cinfo.scale_denom < 8 |
352 | && cinfo.output_width >= IW*4 |
803 | && cinfo.output_width >= iw*4 |
353 | && cinfo.output_height >= IH*4) |
804 | && cinfo.output_height >= ih*4) |
354 | { |
805 | { |
355 | cinfo.scale_denom <<= 1; |
806 | cinfo.scale_denom <<= 1; |
356 | jpeg_calc_output_dimensions (&cinfo); |
807 | jpeg_calc_output_dimensions (&cinfo); |
357 | } |
808 | } |
358 | } |
809 | } |
359 | |
810 | |
|
|
811 | if (cinfo.output_components != 3) |
|
|
812 | longjmp (jerr.setjmp_buffer, 3); |
|
|
813 | |
|
|
814 | if (cinfo.jpeg_color_space == JCS_YCCK || cinfo.jpeg_color_space == JCS_CMYK) |
|
|
815 | { |
|
|
816 | cinfo.out_color_space = JCS_CMYK; |
|
|
817 | cinfo.output_components = 4; |
|
|
818 | } |
|
|
819 | |
360 | pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height); |
820 | pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height); |
361 | if (!RETVAL) |
821 | if (!RETVAL) |
362 | longjmp (jerr.setjmp_buffer, 2); |
822 | longjmp (jerr.setjmp_buffer, 2); |
363 | |
823 | |
364 | data = gdk_pixbuf_get_pixels (RETVAL); |
824 | data = gdk_pixbuf_get_pixels (RETVAL); |
365 | rs = gdk_pixbuf_get_rowstride (RETVAL); |
825 | rs = gdk_pixbuf_get_rowstride (RETVAL); |
366 | |
|
|
367 | if (cinfo.output_components != 3) |
|
|
368 | longjmp (jerr.setjmp_buffer, 3); |
|
|
369 | |
826 | |
370 | jpeg_start_decompress (&cinfo); |
827 | jpeg_start_decompress (&cinfo); |
371 | |
828 | |
372 | while (cinfo.output_scanline < cinfo.output_height) |
829 | while (cinfo.output_scanline < cinfo.output_height) |
373 | { |
830 | { |
… | |
… | |
380 | rp [3] = (guchar *)rp [2] + rs; |
837 | rp [3] = (guchar *)rp [2] + rs; |
381 | |
838 | |
382 | jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); |
839 | jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); |
383 | } |
840 | } |
384 | |
841 | |
|
|
842 | if (cinfo.out_color_space == JCS_CMYK) |
|
|
843 | { |
|
|
844 | guchar *end = data + cinfo.output_height * rs; |
|
|
845 | |
|
|
846 | while (data < end) |
|
|
847 | { |
|
|
848 | U32 c = data [0]; |
|
|
849 | U32 m = data [1]; |
|
|
850 | U32 y = data [2]; |
|
|
851 | U32 k = data [3]; |
|
|
852 | |
|
|
853 | if (0) |
|
|
854 | if (cinfo.Adobe_transform == 2) |
|
|
855 | { |
|
|
856 | c ^= 0xff; |
|
|
857 | m ^= 0xff; |
|
|
858 | y ^= 0xff; |
|
|
859 | k ^= 0xff; |
|
|
860 | } |
|
|
861 | |
|
|
862 | data [0] = (c * k + 0x80) / 0xff; |
|
|
863 | data [1] = (m * k + 0x80) / 0xff; |
|
|
864 | data [2] = (y * k + 0x80) / 0xff; |
|
|
865 | data [3] = 0xff; |
|
|
866 | |
|
|
867 | data += 4; |
|
|
868 | } |
|
|
869 | } |
|
|
870 | |
385 | jpeg_finish_decompress (&cinfo); |
871 | jpeg_finish_decompress (&cinfo); |
386 | fclose (fp); |
|
|
387 | jpeg_destroy_decompress (&cinfo); |
872 | jpeg_destroy_decompress (&cinfo); |
|
|
873 | perlinterp_acquire (); |
388 | } |
874 | } |
389 | OUTPUT: |
875 | OUTPUT: |
390 | RETVAL |
876 | RETVAL |
391 | |
877 | |
392 | void |
878 | void |
393 | compare (GdkPixbuf *a, GdkPixbuf *b) |
879 | compare (GdkPixbuf *a, GdkPixbuf *b) |
394 | PPCODE: |
880 | PPCODE: |
|
|
881 | perlinterp_release (); |
395 | { |
882 | { |
396 | int w = gdk_pixbuf_get_width (a); |
883 | int w = gdk_pixbuf_get_width (a); |
397 | int h = gdk_pixbuf_get_height (a); |
884 | int h = gdk_pixbuf_get_height (a); |
|
|
885 | |
398 | int sa = gdk_pixbuf_get_rowstride (a); |
886 | int sa = gdk_pixbuf_get_rowstride (a); |
399 | int sb = gdk_pixbuf_get_rowstride (b); |
887 | int sb = gdk_pixbuf_get_rowstride (b); |
400 | |
888 | |
401 | guchar *pa = gdk_pixbuf_get_pixels (a); |
889 | guchar *pa = gdk_pixbuf_get_pixels (a); |
402 | guchar *pb = gdk_pixbuf_get_pixels (b); |
890 | guchar *pb = gdk_pixbuf_get_pixels (b); |
… | |
… | |
426 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
914 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
427 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
915 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
428 | } |
916 | } |
429 | } |
917 | } |
430 | |
918 | |
|
|
919 | perlinterp_acquire (); |
|
|
920 | |
431 | EXTEND (SP, 2); |
921 | EXTEND (SP, 2); |
432 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
922 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
433 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
923 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
434 | } |
924 | } |
435 | |
925 | |
436 | ############################################################################# |
926 | ############################################################################# |
437 | |
927 | |
438 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
928 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
439 | |
929 | |
|
|
930 | # currently only works for filenames (octet strings) |
|
|
931 | |
440 | SV * |
932 | SV * |
441 | foldcase (SV *pathsv) |
933 | foldcase (SV *pathsv) |
442 | PROTOTYPE: $ |
934 | PROTOTYPE: $ |
443 | CODE: |
935 | CODE: |
444 | { |
936 | { |
445 | STRLEN plen; |
937 | STRLEN plen; |
446 | U8 *path = (U8 *)SvPVutf8 (pathsv, plen); |
938 | U8 *path = (U8 *)SvPV (pathsv, plen); |
447 | U8 *pend = path + plen; |
939 | U8 *pend = path + plen; |
448 | U8 dst [plen * 6 * 3], *dstp = dst; |
940 | U8 dst [plen * 8 * 3], *dstp = dst; |
449 | |
941 | |
450 | while (path < pend) |
942 | while (path < pend) |
451 | { |
943 | { |
452 | U8 ch = *path; |
944 | U8 ch = *path; |
453 | |
945 | |
454 | if (ch >= 'a' && ch <= 'z') |
946 | if (ch >= 'a' && ch <= 'z') |
455 | *dstp++ = *path++; |
947 | *dstp++ = *path++; |
|
|
948 | else if (ch >= 'A' && ch <= 'Z') |
|
|
949 | *dstp++ = *path++ + ('a' - 'A'); |
456 | else if (ch >= '0' && ch <= '9') |
950 | else if (ch >= '0' && ch <= '9') |
457 | { |
951 | { |
|
|
952 | /* version sort, up to 8 digits */ |
458 | STRLEN el, nl = 0; |
953 | STRLEN el, nl = 0; |
459 | while (*path >= '0' && *path <= '9' && path < pend) |
954 | while (*path >= '0' && *path <= '9' && path < pend) |
460 | path++, nl++; |
955 | path++, nl++; |
461 | |
956 | |
462 | for (el = nl; el < 6; el++) |
957 | for (el = nl; el < 8; el++) |
463 | *dstp++ = '0'; |
958 | *dstp++ = '0'; |
464 | |
959 | |
465 | memcpy (dstp, path - nl, nl); |
960 | memcpy (dstp, path - nl, nl); |
466 | dstp += nl; |
961 | dstp += nl; |
467 | } |
962 | } |
|
|
963 | else |
|
|
964 | *dstp++ = *path++; |
|
|
965 | #if 0 |
468 | else |
966 | else |
469 | { |
967 | { |
470 | STRLEN cl; |
968 | STRLEN cl; |
471 | to_utf8_fold (path, dstp, &cl); |
969 | to_utf8_fold (path, dstp, &cl); |
472 | dstp += cl; |
970 | dstp += cl; |
473 | path += is_utf8_char (path); |
971 | path += is_utf8_char (path); |
474 | } |
972 | } |
|
|
973 | #endif |
475 | } |
974 | } |
476 | |
975 | |
477 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
976 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
478 | } |
977 | } |
479 | OUTPUT: |
978 | OUTPUT: |
… | |
… | |
648 | { |
1147 | { |
649 | int i; |
1148 | int i; |
650 | AV *av, *result; |
1149 | AV *av, *result; |
651 | |
1150 | |
652 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
1151 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
653 | croak ("Not an array ref as first argument to make_histogram"); |
1152 | croak ("Not an array ref as first argument to extract_features"); |
654 | |
1153 | |
655 | av = (AV *) SvRV (ar); |
1154 | av = (AV *) SvRV (ar); |
656 | result = newAV (); |
1155 | result = newAV (); |
657 | |
1156 | |
658 | for (i = 0; i <= av_len (av); ++i) |
1157 | for (i = 0; i <= av_len (av); ++i) |
… | |
… | |
711 | feature_finish_pass_2 (&f_v); |
1210 | feature_finish_pass_2 (&f_v); |
712 | } |
1211 | } |
713 | |
1212 | |
714 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
1213 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
715 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
1214 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
716 | hist [6] = f_v.v1 ; hist [7] = f_v.v2 ; hist [8] = f_v.v3 ; |
1215 | hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
717 | |
1216 | |
718 | av_push (result, histsv); |
1217 | av_push (result, histsv); |
719 | } |
1218 | } |
720 | |
1219 | |
721 | RETVAL = newRV_noinc ((SV *)result); |
1220 | RETVAL = newRV_noinc ((SV *)result); |
722 | } |
1221 | } |
723 | OUTPUT: |
1222 | OUTPUT: |
724 | RETVAL |
1223 | RETVAL |
725 | |
1224 | |
726 | SV * |
|
|
727 | make_histograms (SV *ar) |
|
|
728 | CODE: |
|
|
729 | { |
|
|
730 | int i; |
|
|
731 | AV *av, *result; |
|
|
732 | |
|
|
733 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
|
|
734 | croak ("Not an array ref as first argument to make_histogram"); |
|
|
735 | |
|
|
736 | av = (AV *) SvRV (ar); |
|
|
737 | result = newAV (); |
|
|
738 | |
|
|
739 | for (i = 0; i <= av_len (av); ++i) |
|
|
740 | { |
|
|
741 | const int HISTSIZE = 64; |
|
|
742 | |
|
|
743 | int j; |
|
|
744 | SV *sv = *av_fetch (av, i, 1); |
|
|
745 | STRLEN len; |
|
|
746 | char *buf = SvPVbyte (sv, len); |
|
|
747 | |
|
|
748 | int tmphist[HISTSIZE]; |
|
|
749 | float *hist; |
|
|
750 | |
|
|
751 | SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); |
|
|
752 | SvPOK_on (histsv); |
|
|
753 | SvCUR_set (histsv, HISTSIZE * sizeof (float)); |
|
|
754 | hist = (float *)SvPVX (histsv); |
|
|
755 | |
|
|
756 | Zero (tmphist, sizeof (tmphist), char); |
|
|
757 | |
|
|
758 | for (j = len; j--; ) |
|
|
759 | { |
|
|
760 | unsigned int idx |
|
|
761 | = ((*buf & 0xc0) >> 2) |
|
|
762 | | ((*buf & 0x18) >> 1) |
|
|
763 | | (*buf & 0x03); |
|
|
764 | |
|
|
765 | ++tmphist[idx]; |
|
|
766 | ++buf; |
|
|
767 | } |
|
|
768 | |
|
|
769 | for (j = 0; j < HISTSIZE; ++j) |
|
|
770 | hist[j] = (float)tmphist[j] / (len + 1e-30); |
|
|
771 | |
|
|
772 | av_push (result, histsv); |
|
|
773 | } |
|
|
774 | |
|
|
775 | RETVAL = newRV_noinc ((SV *)result); |
|
|
776 | } |
|
|
777 | OUTPUT: |
|
|
778 | RETVAL |
|
|
779 | |
|
|
780 | |
|
|