… | |
… | |
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, |
… | |
… | |
68 | |
93 | |
69 | *h &= 255; |
94 | *h &= 255; |
70 | } |
95 | } |
71 | } |
96 | } |
72 | |
97 | |
|
|
98 | struct feature { |
|
|
99 | float v1, v2, v3; // mean, square, cube |
|
|
100 | int n; |
|
|
101 | }; |
|
|
102 | |
|
|
103 | static void |
|
|
104 | feature_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. |
|
|
114 | static void |
|
|
115 | feature_update_pass_1 (struct feature *f, unsigned int v) |
|
|
116 | { |
|
|
117 | f->v1 += v; |
|
|
118 | f->n += 1; |
|
|
119 | } |
|
|
120 | |
|
|
121 | static void |
|
|
122 | feature_finish_pass_1 (struct feature *f) |
|
|
123 | { |
|
|
124 | if (f->n < 1) |
|
|
125 | return; |
|
|
126 | |
|
|
127 | f->v1 /= f->n; |
|
|
128 | } |
|
|
129 | |
|
|
130 | static void |
|
|
131 | feature_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 | |
|
|
139 | static void |
|
|
140 | feature_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 | |
73 | static guint32 a85_val; |
153 | static guint32 a85_val; |
74 | static guint a85_cnt; |
154 | static guint a85_cnt; |
75 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
155 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
76 | |
156 | |
77 | static void |
157 | static void |
… | |
… | |
89 | if (!--a85_cnt) |
169 | if (!--a85_cnt) |
90 | { |
170 | { |
91 | a85_cnt = 4; |
171 | a85_cnt = 4; |
92 | if (a85_val) |
172 | if (a85_val) |
93 | { |
173 | { |
94 | a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
174 | a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85; |
95 | a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
175 | a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85; |
96 | a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
176 | a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85; |
97 | a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
177 | a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85; |
98 | a85_ptr[0] = (a85_val ) + 33; |
178 | a85_ptr[0] = (a85_val ) + 33; |
99 | |
179 | |
100 | a85_ptr += 5; |
180 | a85_ptr += 5; |
101 | } |
181 | } |
… | |
… | |
124 | |
204 | |
125 | PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
205 | PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); |
126 | } |
206 | } |
127 | |
207 | |
128 | ///////////////////////////////////////////////////////////////////////////// |
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 | ///////////////////////////////////////////////////////////////////////////// |
129 | |
303 | |
130 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
304 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV |
131 | |
305 | |
132 | PROTOTYPES: ENABLE |
306 | PROTOTYPES: ENABLE |
133 | |
307 | |
|
|
308 | # calculate the common prefix length of two strings |
134 | # missing function in perl. really :) |
309 | # missing function in perl. really :) |
135 | int |
310 | int |
136 | common_prefix_length (a, b) |
311 | common_prefix_length (a, b) |
137 | unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
312 | unsigned char *a = (unsigned char *)SvPVutf8_nolen ($arg); |
138 | unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
313 | unsigned char *b = (unsigned char *)SvPVutf8_nolen ($arg); |
… | |
… | |
147 | |
322 | |
148 | OUTPUT: |
323 | OUTPUT: |
149 | RETVAL |
324 | RETVAL |
150 | |
325 | |
151 | const char * |
326 | const char * |
152 | magic (const char *path) |
327 | magic (SV *path_or_data) |
|
|
328 | ALIAS: |
|
|
329 | magic = 0 |
|
|
330 | magic_mime = 1 |
|
|
331 | magic_buffer = 2 |
|
|
332 | magic_buffer_mime = 3 |
153 | CODE: |
333 | CODE: |
154 | { |
334 | { |
155 | static magic_t cookie; |
335 | STRLEN len; |
|
|
336 | char *data = SvPVbyte (path_or_data, len); |
156 | |
337 | |
157 | if (!cookie) |
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); |
158 | { |
344 | } |
159 | cookie = magic_open (MAGIC_NONE); |
|
|
160 | |
345 | |
161 | if (cookie) |
346 | perlinterp_release (); |
162 | magic_load (cookie, 0); |
|
|
163 | else |
|
|
164 | XSRETURN_UNDEF; |
|
|
165 | } |
|
|
166 | |
347 | |
167 | RETVAL = magic_file (cookie, path); |
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 (); |
168 | } |
353 | } |
169 | OUTPUT: |
354 | OUTPUT: |
170 | RETVAL |
355 | RETVAL |
171 | |
356 | |
172 | const char * |
|
|
173 | magic_mime (const char *path) |
|
|
174 | CODE: |
|
|
175 | { |
|
|
176 | static magic_t cookie; |
|
|
177 | |
|
|
178 | if (!cookie) |
|
|
179 | { |
|
|
180 | cookie = magic_open (MAGIC_MIME); |
|
|
181 | |
|
|
182 | if (cookie) |
|
|
183 | magic_load (cookie, 0); |
|
|
184 | else |
|
|
185 | XSRETURN_UNDEF; |
|
|
186 | } |
|
|
187 | |
|
|
188 | RETVAL = magic_file (cookie, path); |
|
|
189 | } |
|
|
190 | OUTPUT: |
|
|
191 | RETVAL |
|
|
192 | |
|
|
193 | # missing in Gtk2 perl module |
357 | # missing/broken in Gtk2 perl module |
|
|
358 | |
|
|
359 | void |
|
|
360 | gdk_window_clear_hints (GdkWindow *window) |
|
|
361 | CODE: |
|
|
362 | gdk_window_set_geometry_hints (window, 0, 0); |
194 | |
363 | |
195 | gboolean |
364 | gboolean |
196 | gdk_net_wm_supports (GdkAtom property) |
365 | gdk_net_wm_supports (GdkAtom property) |
197 | CODE: |
366 | CODE: |
198 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
367 | #if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE) |
… | |
… | |
204 | RETVAL |
373 | RETVAL |
205 | |
374 | |
206 | GdkPixbuf_noinc * |
375 | GdkPixbuf_noinc * |
207 | dealpha_expose (GdkPixbuf *pb) |
376 | dealpha_expose (GdkPixbuf *pb) |
208 | CODE: |
377 | CODE: |
|
|
378 | perlinterp_release (); |
209 | { |
379 | { |
210 | int w = gdk_pixbuf_get_width (pb); |
380 | int w = gdk_pixbuf_get_width (pb); |
211 | int h = gdk_pixbuf_get_height (pb); |
381 | int h = gdk_pixbuf_get_height (pb); |
212 | int bpp = gdk_pixbuf_get_n_channels (pb); |
382 | int bpp = gdk_pixbuf_get_n_channels (pb); |
213 | int x, y, i; |
383 | int x, y, i; |
… | |
… | |
222 | for (x = 0; x < w; x++) |
392 | for (x = 0; x < w; x++) |
223 | for (y = 0; y < h; y++) |
393 | for (y = 0; y < h; y++) |
224 | for (i = 0; i < 3; i++) |
394 | for (i = 0; i < 3; i++) |
225 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
395 | dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i]; |
226 | } |
396 | } |
|
|
397 | perlinterp_acquire (); |
227 | OUTPUT: |
398 | OUTPUT: |
228 | RETVAL |
399 | RETVAL |
229 | |
400 | |
230 | GdkPixbuf_noinc * |
401 | GdkPixbuf_noinc * |
231 | rotate (GdkPixbuf *pb, int angle) |
402 | rotate (GdkPixbuf *pb, int angle) |
232 | CODE: |
403 | CODE: |
|
|
404 | perlinterp_release (); |
|
|
405 | if (angle < 0) |
|
|
406 | angle += 360; |
233 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
407 | RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE |
234 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
408 | : angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE |
235 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
409 | : angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN |
236 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
410 | : angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE |
237 | : angle); |
411 | : angle); |
|
|
412 | perlinterp_acquire (); |
238 | OUTPUT: |
413 | OUTPUT: |
239 | RETVAL |
414 | RETVAL |
240 | |
415 | |
|
|
416 | const char * |
|
|
417 | filetype (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, 0x0d, 0xa, 0x87, 0x0a |
|
|
424 | }; |
|
|
425 | |
|
|
426 | if (data_len >= 20 |
|
|
427 | && data[0] == 0xff |
|
|
428 | && data[1] == 0xd8 |
|
|
429 | && data[2] == 0xff) |
|
|
430 | RETVAL = "image/jpeg"; |
|
|
431 | else if (data_len >= 12 |
|
|
432 | && data[ 0] == (U8)'R' |
|
|
433 | && data[ 1] == (U8)'I' |
|
|
434 | && data[ 2] == (U8)'F' |
|
|
435 | && data[ 3] == (U8)'F' |
|
|
436 | && data[ 8] == (U8)'W' |
|
|
437 | && data[ 9] == (U8)'E' |
|
|
438 | && data[10] == (U8)'B' |
|
|
439 | && data[11] == (U8)'P') |
|
|
440 | RETVAL = "image/webp"; |
|
|
441 | else if (data_len >= 16 |
|
|
442 | && data[ 0] == 0x89 |
|
|
443 | && data[ 1] == (U8)'P' |
|
|
444 | && data[ 2] == (U8)'N' |
|
|
445 | && data[ 3] == (U8)'G' |
|
|
446 | && data[ 4] == 0x0d |
|
|
447 | && data[ 5] == 0x0a |
|
|
448 | && data[ 6] == 0x1a |
|
|
449 | && data[ 7] == 0x0a) |
|
|
450 | RETVAL = "image/png"; |
|
|
451 | else if (data_len >= sizeof (jxl_header) && memcmp (data, jxl_header, sizeof (jxl_header)) == 0) |
|
|
452 | RETVAL = "image/jxl"; // todo: might want to use JxlSignatureCheck |
|
|
453 | else if (data_len >= 2 |
|
|
454 | && data[0] == 0xff |
|
|
455 | && data[1] == 0x0a) |
|
|
456 | RETVAL = "image/jxl"; |
|
|
457 | else if (data_len >= 13 |
|
|
458 | && data[0] == 'G' |
|
|
459 | && data[1] == 'I' |
|
|
460 | && data[2] == 'F' |
|
|
461 | && data[3] == '8' |
|
|
462 | //&& (data[4] == '7' || data[4] == '9') |
|
|
463 | && data[5] == 'a') |
|
|
464 | { |
|
|
465 | RETVAL = "image/gif"; |
|
|
466 | |
|
|
467 | // now see if its animated - we require the netscape application header for this |
|
|
468 | int ofs = 13; |
|
|
469 | |
|
|
470 | if (data[10] & 0x80) |
|
|
471 | ofs += (1 << ((data[10] & 7) + 1)) * 3; |
|
|
472 | |
|
|
473 | if (data_len >= ofs + 2 + 1 + 11) |
|
|
474 | { |
|
|
475 | |
|
|
476 | // skip a graphic control extension block. we assume |
|
|
477 | // there is at most one such block - while the NAB |
|
|
478 | // has to come firstz, some files do not obey this |
|
|
479 | if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9) |
|
|
480 | ofs += 3 + data[ofs + 2] + 1; |
|
|
481 | |
|
|
482 | if (data_len >= ofs + 2 + 1 + 11) |
|
|
483 | if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1)) |
|
|
484 | RETVAL = "video/gif"; |
|
|
485 | } |
|
|
486 | } |
|
|
487 | |
|
|
488 | else |
|
|
489 | XSRETURN_UNDEF; |
|
|
490 | } |
|
|
491 | OUTPUT: |
|
|
492 | RETVAL |
|
|
493 | |
241 | GdkPixbuf_noinc * |
494 | GdkPixbuf_noinc * |
242 | load_jpeg (SV *path, int thumbnail=0) |
495 | decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
|
|
496 | CODE: |
|
|
497 | { |
|
|
498 | #if WEBP |
|
|
499 | STRLEN data_size; |
|
|
500 | int alpha; |
|
|
501 | WebPData data; |
|
|
502 | WebPDemuxer *demux; |
|
|
503 | WebPIterator iter; |
|
|
504 | WebPDecoderConfig config; |
|
|
505 | int inw, inh; |
|
|
506 | |
|
|
507 | data.bytes = (uint8_t *)SvPVbyte (image_data, data_size); |
|
|
508 | data.size = data_size; |
|
|
509 | |
|
|
510 | perlinterp_release (); |
|
|
511 | |
|
|
512 | RETVAL = 0; |
|
|
513 | |
|
|
514 | if (!(demux = WebPDemux (&data))) |
|
|
515 | goto err_demux; |
|
|
516 | |
|
|
517 | if (!WebPDemuxGetFrame (demux, 1, &iter)) |
|
|
518 | goto err_iter; |
|
|
519 | |
|
|
520 | if (!WebPInitDecoderConfig (&config)) |
|
|
521 | goto err_iter; |
|
|
522 | |
|
|
523 | config.options.use_threads = 1; |
|
|
524 | |
|
|
525 | if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK) |
|
|
526 | goto err_iter; |
|
|
527 | |
|
|
528 | inw = config.input.width; |
|
|
529 | inh = config.input.height; |
|
|
530 | |
|
|
531 | if (thumbnail) |
|
|
532 | { |
|
|
533 | if (inw * ih > inh * iw) |
|
|
534 | ih = (iw * inh + inw - 1) / inw; |
|
|
535 | else |
|
|
536 | iw = (ih * inw + inh - 1) / inh; |
|
|
537 | |
|
|
538 | config.options.bypass_filtering = 1; |
|
|
539 | config.options.no_fancy_upsampling = 1; |
|
|
540 | |
|
|
541 | config.options.use_scaling = 1; |
|
|
542 | config.options.scaled_width = iw; |
|
|
543 | config.options.scaled_height = ih; |
|
|
544 | } |
|
|
545 | else |
|
|
546 | { |
|
|
547 | iw = inw; |
|
|
548 | ih = inh; |
|
|
549 | } |
|
|
550 | |
|
|
551 | alpha = !!config.input.has_alpha; |
|
|
552 | |
|
|
553 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih); |
|
|
554 | if (!RETVAL) |
|
|
555 | goto err_iter; |
|
|
556 | |
|
|
557 | config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB; |
|
|
558 | config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL); |
|
|
559 | config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL); |
|
|
560 | config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL); |
|
|
561 | config.output.is_external_memory = 1; |
|
|
562 | |
|
|
563 | if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK) |
|
|
564 | { |
|
|
565 | g_object_unref (RETVAL); |
|
|
566 | RETVAL = 0; |
|
|
567 | goto err_iter; |
|
|
568 | } |
|
|
569 | |
|
|
570 | err_iter: |
|
|
571 | WebPDemuxReleaseIterator (&iter); |
|
|
572 | err_demux: |
|
|
573 | WebPDemuxDelete (demux); |
|
|
574 | |
|
|
575 | perlinterp_acquire (); |
|
|
576 | #else |
|
|
577 | croak ("load_webp: webp not enabled at compile time"); |
|
|
578 | #endif |
|
|
579 | } |
|
|
580 | OUTPUT: |
|
|
581 | RETVAL |
|
|
582 | |
|
|
583 | GdkPixbuf_noinc * |
|
|
584 | decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
|
|
585 | CODE: |
|
|
586 | { |
|
|
587 | #if JXL |
|
|
588 | JxlDecoder *dec; |
|
|
589 | JxlBasicInfo info; |
|
|
590 | const uint8_t *next_in = (uint8_t *)SvPVbyte_nolen (image_data); |
|
|
591 | size_t avail_in = SvCUR (image_data); |
|
|
592 | const char *error = 0; |
|
|
593 | JxlDecoderStatus status; |
|
|
594 | void *runner = 0; |
|
|
595 | struct bmff_box box; |
|
|
596 | |
|
|
597 | RETVAL = 0; |
|
|
598 | |
|
|
599 | perlinterp_release (); |
|
|
600 | |
|
|
601 | dec = JxlDecoderCreate (0); |
|
|
602 | |
|
|
603 | error = "JxlDecoderCreate failed"; |
|
|
604 | if (!dec) |
|
|
605 | goto done; |
|
|
606 | |
|
|
607 | runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ()); |
|
|
608 | |
|
|
609 | status = JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner); |
|
|
610 | error = "JxlDecoderSetParallelRunner failed"; |
|
|
611 | if (status != JXL_DEC_SUCCESS) |
|
|
612 | goto done; |
|
|
613 | |
|
|
614 | error = "JxlDecoderSubscribeEvents failed"; |
|
|
615 | status = JxlDecoderSubscribeEvents (dec, JXL_DEC_BASIC_INFO | JXL_DEC_FULL_IMAGE); |
|
|
616 | if (status != JXL_DEC_SUCCESS) |
|
|
617 | goto done; |
|
|
618 | |
|
|
619 | status = JxlDecoderSetInput (dec, next_in, avail_in); |
|
|
620 | error = "JxlDecoderSetInput failed"; |
|
|
621 | if (status != JXL_DEC_SUCCESS) |
|
|
622 | goto done; |
|
|
623 | |
|
|
624 | for (;;) |
|
|
625 | { |
|
|
626 | status = JxlDecoderProcessInput (dec); |
|
|
627 | |
|
|
628 | switch (status) |
|
|
629 | { |
|
|
630 | case JXL_DEC_FULL_IMAGE: |
|
|
631 | error = 0; |
|
|
632 | goto done; |
|
|
633 | |
|
|
634 | case JXL_DEC_ERROR: |
|
|
635 | error = "JxlDecoderProcessInput failed"; |
|
|
636 | goto done; |
|
|
637 | |
|
|
638 | case JXL_DEC_NEED_MORE_INPUT: |
|
|
639 | error = "incomplete file"; |
|
|
640 | goto done; |
|
|
641 | |
|
|
642 | case JXL_DEC_SUCCESS: |
|
|
643 | error = "incomplete decode"; |
|
|
644 | goto done; |
|
|
645 | |
|
|
646 | case JXL_DEC_BASIC_INFO: |
|
|
647 | { |
|
|
648 | status = JxlDecoderGetBasicInfo (dec, &info); |
|
|
649 | error = "JxlDecoderGetBasicInfo failed"; |
|
|
650 | if (status != JXL_DEC_SUCCESS) |
|
|
651 | goto done; |
|
|
652 | |
|
|
653 | RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize); |
|
|
654 | error = "unable to allocate pixbuf"; |
|
|
655 | if (!RETVAL) |
|
|
656 | goto done; |
|
|
657 | |
|
|
658 | JxlPixelFormat format = { |
|
|
659 | info.alpha_bits ? 4 : 3, |
|
|
660 | JXL_TYPE_UINT8, |
|
|
661 | JXL_NATIVE_ENDIAN, |
|
|
662 | gdk_pixbuf_get_rowstride (RETVAL) |
|
|
663 | }; |
|
|
664 | |
|
|
665 | // cannot use gdk_pixbuf_get_byte_length because that does |
|
|
666 | // not return the size of the buffer, but the size of the buffer without |
|
|
667 | // the last padding bytes. the internal buffer is rowstride * ysize, |
|
|
668 | // and this is what the jxl decoder needs. none of this is documented |
|
|
669 | // in either library, of course. |
|
|
670 | |
|
|
671 | status = JxlDecoderSetImageOutBuffer ( |
|
|
672 | dec, |
|
|
673 | &format, |
|
|
674 | gdk_pixbuf_get_pixels (RETVAL), |
|
|
675 | gdk_pixbuf_get_rowstride (RETVAL) * info.ysize |
|
|
676 | ); |
|
|
677 | error = "JxlDecoderSetImageOutBuffer failed"; |
|
|
678 | if (status != JXL_DEC_SUCCESS) |
|
|
679 | goto done; |
|
|
680 | } |
|
|
681 | break; |
|
|
682 | |
|
|
683 | default: |
|
|
684 | error = "unexpected event"; |
|
|
685 | goto done; |
|
|
686 | } |
|
|
687 | } |
|
|
688 | |
|
|
689 | done: |
|
|
690 | if (dec) |
|
|
691 | JxlDecoderDestroy (dec); |
|
|
692 | |
|
|
693 | if (runner) |
|
|
694 | JxlThreadParallelRunnerDestroy (runner); |
|
|
695 | |
|
|
696 | perlinterp_acquire (); |
|
|
697 | |
|
|
698 | if (error) |
|
|
699 | { |
|
|
700 | if (RETVAL) |
|
|
701 | g_object_unref (RETVAL); |
|
|
702 | |
|
|
703 | croak ("load_jxl: %s (status %d)", error, status); |
|
|
704 | } |
|
|
705 | #else |
|
|
706 | croak ("load_jxl: jpeg-xl not enabled at compile time"); |
|
|
707 | #endif |
|
|
708 | } |
|
|
709 | OUTPUT: |
|
|
710 | RETVAL |
|
|
711 | |
|
|
712 | GdkPixbuf_noinc * |
|
|
713 | decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0) |
243 | CODE: |
714 | CODE: |
244 | { |
715 | { |
245 | struct jpeg_decompress_struct cinfo; |
716 | struct jpeg_decompress_struct cinfo; |
246 | struct jpg_err_mgr jerr; |
717 | struct jpg_err_mgr jerr; |
247 | guchar *data; |
|
|
248 | int rs; |
718 | int rs; |
249 | FILE *fp; |
|
|
250 | volatile GdkPixbuf *pb = 0; |
719 | volatile GdkPixbuf *pb = 0; |
|
|
720 | STRLEN data_len; |
|
|
721 | guchar *data = SvPVbyte (image_data, data_len); |
251 | |
722 | |
252 | RETVAL = 0; |
723 | RETVAL = 0; |
253 | |
724 | |
254 | fp = fopen (SvPVbyte_nolen (path), "rb"); |
725 | perlinterp_release (); |
255 | |
|
|
256 | if (!fp) |
|
|
257 | XSRETURN_UNDEF; |
|
|
258 | |
726 | |
259 | cinfo.err = jpeg_std_error (&jerr.err); |
727 | cinfo.err = jpeg_std_error (&jerr.err); |
260 | |
728 | |
261 | jerr.err.error_exit = cv_error_exit; |
729 | jerr.err.error_exit = cv_error_exit; |
262 | jerr.err.output_message = cv_error_output; |
730 | jerr.err.output_message = cv_error_output; |
263 | |
731 | |
264 | if ((rs = setjmp (jerr.setjmp_buffer))) |
732 | if ((rs = setjmp (jerr.setjmp_buffer))) |
265 | { |
733 | { |
266 | fclose (fp); |
|
|
267 | jpeg_destroy_decompress (&cinfo); |
734 | jpeg_destroy_decompress (&cinfo); |
268 | |
735 | |
269 | if (pb) |
736 | if (pb) |
270 | g_object_unref ((gpointer)pb); |
737 | g_object_unref ((gpointer)pb); |
271 | |
738 | |
|
|
739 | perlinterp_acquire (); |
272 | XSRETURN_UNDEF; |
740 | XSRETURN_UNDEF; |
273 | } |
741 | } |
274 | |
742 | |
|
|
743 | if (!data_len) |
|
|
744 | longjmp (jerr.setjmp_buffer, 4); |
|
|
745 | |
275 | jpeg_create_decompress (&cinfo); |
746 | jpeg_create_decompress (&cinfo); |
|
|
747 | cv_jpeg_mem_src (&cinfo, data, data_len); |
276 | |
748 | |
277 | jpeg_stdio_src (&cinfo, fp); |
|
|
278 | jpeg_read_header (&cinfo, TRUE); |
749 | jpeg_read_header (&cinfo, TRUE); |
279 | |
750 | |
280 | cinfo.dct_method = JDCT_DEFAULT; |
751 | cinfo.dct_method = JDCT_DEFAULT; |
281 | cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ |
752 | cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */ |
282 | cinfo.do_block_smoothing = FALSE; |
753 | cinfo.do_block_smoothing = FALSE; |
… | |
… | |
292 | { |
763 | { |
293 | cinfo.dct_method = JDCT_FASTEST; |
764 | cinfo.dct_method = JDCT_FASTEST; |
294 | cinfo.do_fancy_upsampling = FALSE; |
765 | cinfo.do_fancy_upsampling = FALSE; |
295 | |
766 | |
296 | while (cinfo.scale_denom < 8 |
767 | while (cinfo.scale_denom < 8 |
297 | && cinfo.output_width >= IW*4 |
768 | && cinfo.output_width >= iw*4 |
298 | && cinfo.output_height >= IH*4) |
769 | && cinfo.output_height >= ih*4) |
299 | { |
770 | { |
300 | cinfo.scale_denom <<= 1; |
771 | cinfo.scale_denom <<= 1; |
301 | jpeg_calc_output_dimensions (&cinfo); |
772 | jpeg_calc_output_dimensions (&cinfo); |
302 | } |
773 | } |
303 | } |
774 | } |
304 | |
775 | |
|
|
776 | if (cinfo.output_components != 3) |
|
|
777 | longjmp (jerr.setjmp_buffer, 3); |
|
|
778 | |
|
|
779 | if (cinfo.jpeg_color_space == JCS_YCCK || cinfo.jpeg_color_space == JCS_CMYK) |
|
|
780 | { |
|
|
781 | cinfo.out_color_space = JCS_CMYK; |
|
|
782 | cinfo.output_components = 4; |
|
|
783 | } |
|
|
784 | |
305 | pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height); |
785 | pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height); |
306 | if (!RETVAL) |
786 | if (!RETVAL) |
307 | longjmp (jerr.setjmp_buffer, 2); |
787 | longjmp (jerr.setjmp_buffer, 2); |
308 | |
788 | |
309 | data = gdk_pixbuf_get_pixels (RETVAL); |
789 | data = gdk_pixbuf_get_pixels (RETVAL); |
310 | rs = gdk_pixbuf_get_rowstride (RETVAL); |
790 | rs = gdk_pixbuf_get_rowstride (RETVAL); |
311 | |
|
|
312 | if (cinfo.output_components != 3) |
|
|
313 | longjmp (jerr.setjmp_buffer, 3); |
|
|
314 | |
791 | |
315 | jpeg_start_decompress (&cinfo); |
792 | jpeg_start_decompress (&cinfo); |
316 | |
793 | |
317 | while (cinfo.output_scanline < cinfo.output_height) |
794 | while (cinfo.output_scanline < cinfo.output_height) |
318 | { |
795 | { |
… | |
… | |
325 | rp [3] = (guchar *)rp [2] + rs; |
802 | rp [3] = (guchar *)rp [2] + rs; |
326 | |
803 | |
327 | jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); |
804 | jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4); |
328 | } |
805 | } |
329 | |
806 | |
|
|
807 | if (cinfo.out_color_space == JCS_CMYK) |
|
|
808 | { |
|
|
809 | guchar *end = data + cinfo.output_height * rs; |
|
|
810 | |
|
|
811 | while (data < end) |
|
|
812 | { |
|
|
813 | U32 c = data [0]; |
|
|
814 | U32 m = data [1]; |
|
|
815 | U32 y = data [2]; |
|
|
816 | U32 k = data [3]; |
|
|
817 | |
|
|
818 | if (0) |
|
|
819 | if (cinfo.Adobe_transform == 2) |
|
|
820 | { |
|
|
821 | c ^= 0xff; |
|
|
822 | m ^= 0xff; |
|
|
823 | y ^= 0xff; |
|
|
824 | k ^= 0xff; |
|
|
825 | } |
|
|
826 | |
|
|
827 | data [0] = (c * k + 0x80) / 0xff; |
|
|
828 | data [1] = (m * k + 0x80) / 0xff; |
|
|
829 | data [2] = (y * k + 0x80) / 0xff; |
|
|
830 | data [3] = 0xff; |
|
|
831 | |
|
|
832 | data += 4; |
|
|
833 | } |
|
|
834 | } |
|
|
835 | |
330 | jpeg_finish_decompress (&cinfo); |
836 | jpeg_finish_decompress (&cinfo); |
331 | fclose (fp); |
|
|
332 | jpeg_destroy_decompress (&cinfo); |
837 | jpeg_destroy_decompress (&cinfo); |
|
|
838 | perlinterp_acquire (); |
333 | } |
839 | } |
334 | OUTPUT: |
840 | OUTPUT: |
335 | RETVAL |
841 | RETVAL |
336 | |
842 | |
337 | void |
843 | void |
338 | compare (GdkPixbuf *a, GdkPixbuf *b) |
844 | compare (GdkPixbuf *a, GdkPixbuf *b) |
339 | PPCODE: |
845 | PPCODE: |
|
|
846 | perlinterp_release (); |
340 | { |
847 | { |
341 | int w = gdk_pixbuf_get_width (a); |
848 | int w = gdk_pixbuf_get_width (a); |
342 | int h = gdk_pixbuf_get_height (a); |
849 | int h = gdk_pixbuf_get_height (a); |
|
|
850 | |
343 | int sa = gdk_pixbuf_get_rowstride (a); |
851 | int sa = gdk_pixbuf_get_rowstride (a); |
344 | int sb = gdk_pixbuf_get_rowstride (b); |
852 | int sb = gdk_pixbuf_get_rowstride (b); |
345 | |
853 | |
346 | guchar *pa = gdk_pixbuf_get_pixels (a); |
854 | guchar *pa = gdk_pixbuf_get_pixels (a); |
347 | guchar *pb = gdk_pixbuf_get_pixels (b); |
855 | guchar *pb = gdk_pixbuf_get_pixels (b); |
… | |
… | |
371 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
879 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
372 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
880 | d = ((int)*pa_++) - ((int)*pb_++); diff += d*d; peak = MAX (peak, abs (d)); |
373 | } |
881 | } |
374 | } |
882 | } |
375 | |
883 | |
|
|
884 | perlinterp_acquire (); |
|
|
885 | |
376 | EXTEND (SP, 2); |
886 | EXTEND (SP, 2); |
377 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
887 | PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.))))); |
378 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
888 | PUSHs (sv_2mortal (newSVnv (peak / 255.))); |
379 | } |
889 | } |
380 | |
890 | |
381 | ############################################################################# |
891 | ############################################################################# |
382 | |
892 | |
383 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
893 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer |
384 | |
894 | |
|
|
895 | # currently only works for filenames (octet strings) |
|
|
896 | |
385 | SV * |
897 | SV * |
386 | foldcase (SV *pathsv) |
898 | foldcase (SV *pathsv) |
387 | PROTOTYPE: $ |
899 | PROTOTYPE: $ |
388 | CODE: |
900 | CODE: |
389 | { |
901 | { |
390 | STRLEN plen; |
902 | STRLEN plen; |
391 | U8 *path = (U8 *)SvPVutf8 (pathsv, plen); |
903 | U8 *path = (U8 *)SvPV (pathsv, plen); |
392 | U8 *pend = path + plen; |
904 | U8 *pend = path + plen; |
393 | U8 dst [plen * 6 * 3], *dstp = dst; |
905 | U8 dst [plen * 8 * 3], *dstp = dst; |
394 | |
906 | |
395 | while (path < pend) |
907 | while (path < pend) |
396 | { |
908 | { |
397 | U8 ch = *path; |
909 | U8 ch = *path; |
398 | |
910 | |
399 | if (ch >= 'a' && ch <= 'z') |
911 | if (ch >= 'a' && ch <= 'z') |
400 | *dstp++ = *path++; |
912 | *dstp++ = *path++; |
|
|
913 | else if (ch >= 'A' && ch <= 'Z') |
|
|
914 | *dstp++ = *path++ + ('a' - 'A'); |
401 | else if (ch >= '0' && ch <= '9') |
915 | else if (ch >= '0' && ch <= '9') |
402 | { |
916 | { |
|
|
917 | /* version sort, up to 8 digits */ |
403 | STRLEN el, nl = 0; |
918 | STRLEN el, nl = 0; |
404 | while (*path >= '0' && *path <= '9' && path < pend) |
919 | while (*path >= '0' && *path <= '9' && path < pend) |
405 | path++, nl++; |
920 | path++, nl++; |
406 | |
921 | |
407 | for (el = nl; el < 6; el++) |
922 | for (el = nl; el < 8; el++) |
408 | *dstp++ = '0'; |
923 | *dstp++ = '0'; |
409 | |
924 | |
410 | memcpy (dstp, path - nl, nl); |
925 | memcpy (dstp, path - nl, nl); |
411 | dstp += nl; |
926 | dstp += nl; |
412 | } |
927 | } |
|
|
928 | else |
|
|
929 | *dstp++ = *path++; |
|
|
930 | #if 0 |
413 | else |
931 | else |
414 | { |
932 | { |
415 | STRLEN cl; |
933 | STRLEN cl; |
416 | to_utf8_fold (path, dstp, &cl); |
934 | to_utf8_fold (path, dstp, &cl); |
417 | dstp += cl; |
935 | dstp += cl; |
418 | path += is_utf8_char (path); |
936 | path += is_utf8_char (path); |
419 | } |
937 | } |
|
|
938 | #endif |
420 | } |
939 | } |
421 | |
940 | |
422 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
941 | RETVAL = newSVpvn ((const char *)dst, dstp - dst); |
423 | } |
942 | } |
424 | OUTPUT: |
943 | OUTPUT: |
… | |
… | |
586 | ############################################################################# |
1105 | ############################################################################# |
587 | |
1106 | |
588 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
1107 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
589 | |
1108 | |
590 | SV * |
1109 | SV * |
591 | make_histograms (SV *ar) |
1110 | extract_features (SV *ar) |
592 | CODE: |
1111 | CODE: |
593 | { |
1112 | { |
594 | int i; |
1113 | int i; |
595 | AV *av, *result; |
1114 | AV *av, *result; |
596 | |
1115 | |
597 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
1116 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
598 | croak ("Not an array ref as first argument to make_histogram"); |
1117 | croak ("Not an array ref as first argument to extract_features"); |
599 | |
1118 | |
600 | av = (AV *) SvRV (ar); |
1119 | av = (AV *) SvRV (ar); |
601 | result = newAV (); |
1120 | result = newAV (); |
602 | |
1121 | |
603 | for (i = 0; i <= av_len (av); ++i) |
1122 | for (i = 0; i <= av_len (av); ++i) |
604 | { |
1123 | { |
605 | const int HISTSIZE = 64; |
|
|
606 | |
|
|
607 | int j; |
|
|
608 | SV *sv = *av_fetch (av, i, 1); |
1124 | SV *sv = *av_fetch (av, i, 1); |
609 | STRLEN len; |
|
|
610 | char *buf = SvPVbyte (sv, len); |
|
|
611 | |
|
|
612 | int tmphist[HISTSIZE]; |
|
|
613 | float *hist; |
|
|
614 | |
|
|
615 | SV *histsv = newSV (HISTSIZE * sizeof (float) + 1); |
1125 | SV *histsv = newSV (9 * sizeof (float) + 1); |
|
|
1126 | |
616 | SvPOK_on (histsv); |
1127 | SvPOK_on (histsv); |
617 | SvCUR_set (histsv, HISTSIZE * sizeof (float)); |
1128 | SvCUR_set (histsv, 9 * sizeof (float)); |
618 | hist = (float *)SvPVX (histsv); |
1129 | float *hist = (float *)SvPVX (histsv); |
619 | |
1130 | |
620 | Zero (tmphist, sizeof (tmphist), char); |
1131 | struct feature f_h, f_s, f_v; |
|
|
1132 | feature_init (&f_h); |
|
|
1133 | feature_init (&f_s); |
|
|
1134 | feature_init (&f_v); |
621 | |
1135 | |
622 | for (j = len; j--; ) |
|
|
623 | { |
|
|
624 | unsigned int idx |
|
|
625 | = ((*buf & 0xc0) >> 2) |
|
|
626 | | ((*buf & 0x18) >> 1) |
|
|
627 | | (*buf & 0x03); |
|
|
628 | |
|
|
629 | ++tmphist[idx]; |
|
|
630 | ++buf; |
|
|
631 | } |
|
|
632 | |
1136 | { |
633 | for (j = 0; j < HISTSIZE; ++j) |
1137 | STRLEN len; |
634 | hist[j] = (float)tmphist[j] / (len + 1e-30); |
1138 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
1139 | while (len >= 3) |
|
|
1140 | { |
|
|
1141 | unsigned int r, g, b, h, s, v; |
|
|
1142 | r = *buf++; g = *buf++; b = *buf++; |
|
|
1143 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
1144 | |
|
|
1145 | feature_update_pass_1 (&f_h, h); |
|
|
1146 | feature_update_pass_1 (&f_s, s); |
|
|
1147 | feature_update_pass_1 (&f_v, v); |
|
|
1148 | |
|
|
1149 | len -= 3; |
|
|
1150 | } |
|
|
1151 | |
|
|
1152 | feature_finish_pass_1 (&f_h); |
|
|
1153 | feature_finish_pass_1 (&f_s); |
|
|
1154 | feature_finish_pass_1 (&f_v); |
|
|
1155 | } |
|
|
1156 | |
|
|
1157 | { |
|
|
1158 | STRLEN len; |
|
|
1159 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
1160 | while (len >= 3) |
|
|
1161 | { |
|
|
1162 | unsigned int r, g, b, h, s, v; |
|
|
1163 | r = *buf++; g = *buf++; b = *buf++; |
|
|
1164 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
1165 | |
|
|
1166 | feature_update_pass_2 (&f_h, h); |
|
|
1167 | feature_update_pass_2 (&f_s, s); |
|
|
1168 | feature_update_pass_2 (&f_v, v); |
|
|
1169 | |
|
|
1170 | len -= 3; |
|
|
1171 | } |
|
|
1172 | |
|
|
1173 | feature_finish_pass_2 (&f_h); |
|
|
1174 | feature_finish_pass_2 (&f_s); |
|
|
1175 | feature_finish_pass_2 (&f_v); |
|
|
1176 | } |
|
|
1177 | |
|
|
1178 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
1179 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
1180 | hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5; |
635 | |
1181 | |
636 | av_push (result, histsv); |
1182 | av_push (result, histsv); |
637 | } |
1183 | } |
638 | |
1184 | |
639 | RETVAL = newRV_noinc ((SV *)result); |
1185 | RETVAL = newRV_noinc ((SV *)result); |
640 | } |
1186 | } |
641 | OUTPUT: |
1187 | OUTPUT: |
642 | RETVAL |
1188 | RETVAL |
643 | |
1189 | |
644 | |
|
|