[ViewVC] Diff of: cvs/CV/CV.xs

Comparing CV/CV.xs (file contents):
Revision 1.28 by root, Mon Aug 22 00:30:58 2005 UTC vs.
Revision 1.64 by root, Wed Dec 1 03:32:40 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
		43	typedef char *octet_string;
		44
		45	static magic_t magic_cookie[2]; /* !mime, mime */
24		46
25	struct jpg_err_mgr	47	struct 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
31	static void	53	static void
32	cv_error_exit (j_common_ptr cinfo)	54	cv_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
37	static void	60	static void
38	cv_error_output (j_common_ptr cinfo)	61	cv_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
43	static void	71	static void
44	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,
…		…
65		93
66	*h &= 255;	94	*h &= 255;
67	}	95	}
68	}	96	}
69		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
70	static guint32 a85_val;	153	static guint32 a85_val;
71	static guint a85_cnt;	154	static guint a85_cnt;
72	static guchar a85_buf[LINELENGTH], *a85_ptr;	155	static guchar a85_buf[LINELENGTH], *a85_ptr;
73		156
74	static void	157	static 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
		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	/////////////////////////////////////////////////////////////////////////////
126		303
127	MODULE = Gtk2::CV PACKAGE = Gtk2::CV	304	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
128		305
129	PROTOTYPES: ENABLE	306	PROTOTYPES: ENABLE
130		307
		308	# calculate the common prefix length of two strings
		309	# missing function in perl. really :)
		310	int
		311	common_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
		326	const char *
		327	magic (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
		359	void
		360	gdk_window_clear_hints (GdkWindow *window)
		361	CODE:
		362	gdk_window_set_geometry_hints (window, 0, 0);
132		363
133	gboolean	364	gboolean
134	gdk_net_wm_supports (GdkAtom property)	365	gdk_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)
…		…
142	RETVAL	373	RETVAL
143		374
144	GdkPixbuf_noinc *	375	GdkPixbuf_noinc *
145	dealpha_expose (GdkPixbuf *pb)	376	dealpha_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;
…		…
160	for (x = 0; x < w; x++)	392	for (x = 0; x < w; x++)
161	for (y = 0; y < h; y++)	393	for (y = 0; y < h; y++)
162	for (i = 0; i < 3; i++)	394	for (i = 0; i < 3; i++)
163	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];
164	}	396	}
		397	perlinterp_acquire ();
165	OUTPUT:	398	OUTPUT:
166	RETVAL	399	RETVAL
167		400
168	GdkPixbuf_noinc *	401	GdkPixbuf_noinc *
169	rotate (GdkPixbuf *pb, int angle)	402	rotate (GdkPixbuf *pb, int angle)
170	CODE:	403	CODE:
		404	perlinterp_release ();
		405	if (angle < 0)
		406	angle += 360;
171	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
172	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE	408	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
173	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN	409	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
174	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE	410	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
175	: angle);	411	: angle);
		412	perlinterp_acquire ();
176	OUTPUT:	413	OUTPUT:
177	RETVAL	414	RETVAL
178		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
179	GdkPixbuf_noinc *	494	GdkPixbuf_noinc *
180	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
		596	RETVAL = 0;
		597
		598	perlinterp_release ();
		599
		600	dec = JxlDecoderCreate (0);
		601
		602	error = "JxlDecoderCreate failed";
		603	if (!dec)
		604	goto done;
		605
		606	runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
		607
		608	status = JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner);
		609	error = "JxlDecoderSetParallelRunner failed";
		610	if (status != JXL_DEC_SUCCESS)
		611	goto done;
		612
		613	error = "JxlDecoderSubscribeEvents failed";
		614	status = JxlDecoderSubscribeEvents (dec, JXL_DEC_BASIC_INFO \| JXL_DEC_FULL_IMAGE);
		615	if (status != JXL_DEC_SUCCESS)
		616	goto done;
		617
		618	status = JxlDecoderSetInput (dec, next_in, avail_in);
		619	error = "JxlDecoderSetInput failed";
		620	if (status != JXL_DEC_SUCCESS)
		621	goto done;
		622
		623	for (;;)
		624	{
		625	status = JxlDecoderProcessInput (dec);
		626
		627	switch (status)
		628	{
		629	case JXL_DEC_FULL_IMAGE:
		630	error = 0;
		631	goto done;
		632
		633	case JXL_DEC_ERROR:
		634	error = "JxlDecoderProcessInput failed";
		635	goto done;
		636
		637	case JXL_DEC_NEED_MORE_INPUT:
		638	error = "incomplete file";
		639	goto done;
		640
		641	case JXL_DEC_SUCCESS:
		642	error = "incomplete decode";
		643	goto done;
		644
		645	case JXL_DEC_BASIC_INFO:
		646	{
		647	status = JxlDecoderGetBasicInfo (dec, &info);
		648	error = "JxlDecoderGetBasicInfo failed";
		649	if (status != JXL_DEC_SUCCESS)
		650	goto done;
		651
		652	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
		653	error = "unable to allocate pixbuf";
		654	if (!RETVAL)
		655	goto done;
		656
		657	JxlPixelFormat format = {
		658	info.alpha_bits ? 4 : 3,
		659	JXL_TYPE_UINT8,
		660	JXL_NATIVE_ENDIAN,
		661	gdk_pixbuf_get_rowstride (RETVAL)
		662	};
		663
		664	// cannot use gdk_pixbuf_get_byte_length because that does
		665	// not return the size of the buffer, but the size of the buffer without
		666	// the last padding bytes. the internal buffer is rowstride * ysize,
		667	// and this is what the jxl decoder needs. none of this is documented
		668	// in either library, of course.
		669
		670	status = JxlDecoderSetImageOutBuffer (
		671	dec,
		672	&format,
		673	gdk_pixbuf_get_pixels (RETVAL),
		674	gdk_pixbuf_get_rowstride (RETVAL) * info.ysize
		675	);
		676	error = "JxlDecoderSetImageOutBuffer failed";
		677	if (status != JXL_DEC_SUCCESS)
		678	goto done;
		679	}
		680	break;
		681
		682	default:
		683	error = "unexpected event";
		684	goto done;
		685	}
		686	}
		687
		688	done:
		689	if (dec)
		690	JxlDecoderDestroy (dec);
		691
		692	if (runner)
		693	JxlThreadParallelRunnerDestroy (runner);
		694
		695	perlinterp_acquire ();
		696
		697	if (error)
		698	{
		699	if (RETVAL)
		700	g_object_unref (RETVAL);
		701
		702	croak ("load_jxl: %s (status %d)", error, status);
		703	}
		704	#else
		705	croak ("load_jxl: jpeg-xl not enabled at compile time");
		706	#endif
		707	}
		708	OUTPUT:
		709	RETVAL
		710
		711	GdkPixbuf_noinc *
		712	decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
181	CODE:	713	CODE:
182	{	714	{
183	struct jpeg_decompress_struct cinfo;	715	struct jpeg_decompress_struct cinfo;
184	struct jpg_err_mgr jerr;	716	struct jpg_err_mgr jerr;
185	guchar *data;
186	int rs;	717	int rs;
187	FILE *fp;
188	volatile GdkPixbuf *pb = 0;	718	volatile GdkPixbuf *pb = 0;
189	gchar *filename;	719	STRLEN data_len;
		720	guchar *data = SvPVbyte (image_data, data_len);
190		721
191	RETVAL = 0;	722	RETVAL = 0;
192		723
193	filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0);	724	perlinterp_release ();
194	fp = fopen (filename, "rb");
195	g_free (filename);
196
197	if (!fp)
198	XSRETURN_UNDEF;
199		725
200	cinfo.err = jpeg_std_error (&jerr.err);	726	cinfo.err = jpeg_std_error (&jerr.err);
201		727
202	jerr.err.error_exit = cv_error_exit;	728	jerr.err.error_exit = cv_error_exit;
203	jerr.err.output_message = cv_error_output;	729	jerr.err.output_message = cv_error_output;
204		730
205	if ((rs = setjmp (jerr.setjmp_buffer)))	731	if ((rs = setjmp (jerr.setjmp_buffer)))
206	{	732	{
207	fclose (fp);
208	jpeg_destroy_decompress (&cinfo);	733	jpeg_destroy_decompress (&cinfo);
209		734
210	if (pb)	735	if (pb)
211	g_object_unref ((gpointer)pb);	736	g_object_unref ((gpointer)pb);
212		737
		738	perlinterp_acquire ();
213	XSRETURN_UNDEF;	739	XSRETURN_UNDEF;
214	}	740	}
215		741
		742	if (!data_len)
		743	longjmp (jerr.setjmp_buffer, 4);
		744
216	jpeg_create_decompress (&cinfo);	745	jpeg_create_decompress (&cinfo);
		746	cv_jpeg_mem_src (&cinfo, data, data_len);
217		747
218	jpeg_stdio_src (&cinfo, fp);
219	jpeg_read_header (&cinfo, TRUE);	748	jpeg_read_header (&cinfo, TRUE);
220		749
221	cinfo.dct_method = JDCT_DEFAULT;	750	cinfo.dct_method = JDCT_DEFAULT;
222	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */	751	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
223	cinfo.do_block_smoothing = FALSE;	752	cinfo.do_block_smoothing = FALSE;
…		…
233	{	762	{
234	cinfo.dct_method = JDCT_FASTEST;	763	cinfo.dct_method = JDCT_FASTEST;
235	cinfo.do_fancy_upsampling = FALSE;	764	cinfo.do_fancy_upsampling = FALSE;
236		765
237	while (cinfo.scale_denom < 8	766	while (cinfo.scale_denom < 8
238	&& cinfo.output_width >= IW*4	767	&& cinfo.output_width >= iw*4
239	&& cinfo.output_height >= IH*4)	768	&& cinfo.output_height >= ih*4)
240	{	769	{
241	cinfo.scale_denom <<= 1;	770	cinfo.scale_denom <<= 1;
242	jpeg_calc_output_dimensions (&cinfo);	771	jpeg_calc_output_dimensions (&cinfo);
243	}	772	}
244	}	773	}
245		774
		775	if (cinfo.output_components != 3)
		776	longjmp (jerr.setjmp_buffer, 3);
		777
		778	if (cinfo.jpeg_color_space == JCS_YCCK \|\| cinfo.jpeg_color_space == JCS_CMYK)
		779	{
		780	cinfo.out_color_space = JCS_CMYK;
		781	cinfo.output_components = 4;
		782	}
		783
246	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height);	784	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
247	if (!RETVAL)	785	if (!RETVAL)
248	longjmp (jerr.setjmp_buffer, 2);	786	longjmp (jerr.setjmp_buffer, 2);
249		787
250	data = gdk_pixbuf_get_pixels (RETVAL);	788	data = gdk_pixbuf_get_pixels (RETVAL);
251	rs = gdk_pixbuf_get_rowstride (RETVAL);	789	rs = gdk_pixbuf_get_rowstride (RETVAL);
252
253	if (cinfo.output_components != 3)
254	longjmp (jerr.setjmp_buffer, 3);
255		790
256	jpeg_start_decompress (&cinfo);	791	jpeg_start_decompress (&cinfo);
257		792
258	while (cinfo.output_scanline < cinfo.output_height)	793	while (cinfo.output_scanline < cinfo.output_height)
259	{	794	{
…		…
266	rp [3] = (guchar *)rp [2] + rs;	801	rp [3] = (guchar *)rp [2] + rs;
267		802
268	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);	803	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
269	}	804	}
270		805
		806	if (cinfo.out_color_space == JCS_CMYK)
		807	{
		808	guchar end = data + cinfo.output_height rs;
		809
		810	while (data < end)
		811	{
		812	U32 c = data [0];
		813	U32 m = data [1];
		814	U32 y = data [2];
		815	U32 k = data [3];
		816
		817	if (0)
		818	if (cinfo.Adobe_transform == 2)
		819	{
		820	c ^= 0xff;
		821	m ^= 0xff;
		822	y ^= 0xff;
		823	k ^= 0xff;
		824	}
		825
		826	data [0] = (c * k + 0x80) / 0xff;
		827	data [1] = (m * k + 0x80) / 0xff;
		828	data [2] = (y * k + 0x80) / 0xff;
		829	data [3] = 0xff;
		830
		831	data += 4;
		832	}
		833	}
		834
271	jpeg_finish_decompress (&cinfo);	835	jpeg_finish_decompress (&cinfo);
272	fclose (fp);
273	jpeg_destroy_decompress (&cinfo);	836	jpeg_destroy_decompress (&cinfo);
		837	perlinterp_acquire ();
274	}	838	}
275	OUTPUT:	839	OUTPUT:
276	RETVAL	840	RETVAL
277		841
		842	void
		843	compare (GdkPixbuf a, GdkPixbuf b)
		844	PPCODE:
		845	perlinterp_release ();
		846	{
		847	int w = gdk_pixbuf_get_width (a);
		848	int h = gdk_pixbuf_get_height (a);
		849
		850	int sa = gdk_pixbuf_get_rowstride (a);
		851	int sb = gdk_pixbuf_get_rowstride (b);
		852
		853	guchar *pa = gdk_pixbuf_get_pixels (a);
		854	guchar *pb = gdk_pixbuf_get_pixels (b);
		855
		856	int x, y;
		857
		858	assert (w == gdk_pixbuf_get_width (b));
		859	assert (h == gdk_pixbuf_get_height (b));
		860
		861	assert (gdk_pixbuf_get_n_channels (a) == 3);
		862	assert (gdk_pixbuf_get_n_channels (b) == 3);
		863
		864	double diff = 0.;
		865	int peak = 0;
		866
		867	if (w && h)
		868	for (y = 0; y < h; y++)
		869	{
		870	guchar pa_ = pa + y sa;
		871	guchar pb_ = pb + y sb;
		872
		873	for (x = 0; x < w; x++)
		874	{
		875	int d;
		876
		877	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		878	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));
		880	}
		881	}
		882
		883	perlinterp_acquire ();
		884
		885	EXTEND (SP, 2);
		886	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		887	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		888	}
		889
278	#############################################################################	890	#############################################################################
279		891
280	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	892	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		893
		894	# currently only works for filenames (octet strings)
281		895
282	SV *	896	SV *
283	foldcase (SV *pathsv)	897	foldcase (SV *pathsv)
284	PROTOTYPE: $	898	PROTOTYPE: $
285	CODE:	899	CODE:
286	{	900	{
287	STRLEN plen;	901	STRLEN plen;
288	U8 path = (U8 )SvPVutf8 (pathsv, plen);	902	U8 path = (U8 )SvPV (pathsv, plen);
289	U8 *pend = path + plen;	903	U8 *pend = path + plen;
290	U8 dst [plen * 6 * 3], *dstp = dst;	904	U8 dst [plen * 8 * 3], *dstp = dst;
291		905
292	while (path < pend)	906	while (path < pend)
293	{	907	{
294	U8 ch = *path;	908	U8 ch = *path;
295		909
296	if (ch >= 'a' && ch <= 'z')	910	if (ch >= 'a' && ch <= 'z')
297	dstp++ = path++;	911	dstp++ = path++;
		912	else if (ch >= 'A' && ch <= 'Z')
		913	dstp++ = path++ + ('a' - 'A');
298	else if (ch >= '0' && ch <= '9')	914	else if (ch >= '0' && ch <= '9')
299	{	915	{
		916	/* version sort, up to 8 digits */
300	STRLEN el, nl = 0;	917	STRLEN el, nl = 0;
301	while (path >= '0' && path <= '9' && path < pend)	918	while (path >= '0' && path <= '9' && path < pend)
302	path++, nl++;	919	path++, nl++;
303		920
304	for (el = nl; el < 6; el++)	921	for (el = nl; el < 8; el++)
305	*dstp++ = '0';	922	*dstp++ = '0';
306		923
307	memcpy (dstp, path - nl, nl);	924	memcpy (dstp, path - nl, nl);
308	dstp += nl;	925	dstp += nl;
309	}	926	}
		927	else
		928	dstp++ = path++;
		929	#if 0
310	else	930	else
311	{	931	{
312	STRLEN cl;	932	STRLEN cl;
313	to_utf8_fold (path, dstp, &cl);	933	to_utf8_fold (path, dstp, &cl);
314	dstp += cl;	934	dstp += cl;
315	path += is_utf8_char (path);	935	path += is_utf8_char (path);
316	}	936	}
		937	#endif
317	}	938	}
318		939
319	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	940	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
320	}	941	}
321	OUTPUT:	942	OUTPUT:
…		…
483	#############################################################################	1104	#############################################################################
484		1105
485	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	1106	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
486		1107
487	SV *	1108	SV *
488	make_histograms (SV *ar)	1109	extract_features (SV *ar)
489	CODE:	1110	CODE:
490	{	1111	{
491	int i;	1112	int i;
492	AV av, result;	1113	AV av, result;
493		1114
494	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	1115	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
495	croak ("Not an array ref as first argument to make_histogram");	1116	croak ("Not an array ref as first argument to extract_features");
496		1117
497	av = (AV *) SvRV (ar);	1118	av = (AV *) SvRV (ar);
498	result = newAV ();	1119	result = newAV ();
499		1120
500	for (i = 0; i <= av_len (av); ++i)	1121	for (i = 0; i <= av_len (av); ++i)
501	{	1122	{
502	const int HISTSIZE = 64;
503
504	int j;
505	SV sv = av_fetch (av, i, 1);	1123	SV sv = av_fetch (av, i, 1);
506	STRLEN len;
507	char *buf = SvPVbyte (sv, len);
508
509	int tmphist[HISTSIZE];
510	float *hist;
511
512	SV histsv = newSV (HISTSIZE sizeof (float) + 1);	1124	SV histsv = newSV (9 sizeof (float) + 1);
		1125
513	SvPOK_on (histsv);	1126	SvPOK_on (histsv);
514	SvCUR_set (histsv, HISTSIZE * sizeof (float));	1127	SvCUR_set (histsv, 9 * sizeof (float));
515	hist = (float *)SvPVX (histsv);	1128	float hist = (float )SvPVX (histsv);
516		1129
517	Zero (tmphist, sizeof (tmphist), char);	1130	struct feature f_h, f_s, f_v;
		1131	feature_init (&f_h);
		1132	feature_init (&f_s);
		1133	feature_init (&f_v);
518		1134
519	for (j = len; j--; )
520	{
521	unsigned int idx
522	= ((*buf & 0xc0) >> 2)
523	\| ((*buf & 0x18) >> 1)
524	\| (*buf & 0x03);
525
526	++tmphist[idx];
527	++buf;
528	}
529		1135	{
530	for (j = 0; j < HISTSIZE; ++j)	1136	STRLEN len;
531	hist[j] = (float)tmphist[j] / (len + 1e-30);	1137	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1138	while (len >= 3)
		1139	{
		1140	unsigned int r, g, b, h, s, v;
		1141	r = buf++; g = buf++; b = *buf++;
		1142	rgb_to_hsv (r, g, b, &h, &s, &v);
		1143
		1144	feature_update_pass_1 (&f_h, h);
		1145	feature_update_pass_1 (&f_s, s);
		1146	feature_update_pass_1 (&f_v, v);
		1147
		1148	len -= 3;
		1149	}
		1150
		1151	feature_finish_pass_1 (&f_h);
		1152	feature_finish_pass_1 (&f_s);
		1153	feature_finish_pass_1 (&f_v);
		1154	}
		1155
		1156	{
		1157	STRLEN len;
		1158	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1159	while (len >= 3)
		1160	{
		1161	unsigned int r, g, b, h, s, v;
		1162	r = buf++; g = buf++; b = *buf++;
		1163	rgb_to_hsv (r, g, b, &h, &s, &v);
		1164
		1165	feature_update_pass_2 (&f_h, h);
		1166	feature_update_pass_2 (&f_s, s);
		1167	feature_update_pass_2 (&f_v, v);
		1168
		1169	len -= 3;
		1170	}
		1171
		1172	feature_finish_pass_2 (&f_h);
		1173	feature_finish_pass_2 (&f_s);
		1174	feature_finish_pass_2 (&f_v);
		1175	}
		1176
		1177	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		1178	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		1179	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
532		1180
533	av_push (result, histsv);	1181	av_push (result, histsv);
534	}	1182	}
535		1183
536	RETVAL = newRV_noinc ((SV *)result);	1184	RETVAL = newRV_noinc ((SV *)result);
537	}	1185	}
538	OUTPUT:	1186	OUTPUT:
539	RETVAL	1187	RETVAL
540		1188
541

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing CV/CV.xs (file contents): Revision 1.28 by root, Mon Aug 22 00:30:58 2005 UTC vs. Revision 1.64 by root, Wed Dec 1 03:32:40 2021 UTC

Diff Legend

Comparing CV/CV.xs (file contents):
Revision 1.28 by root, Mon Aug 22 00:30:58 2005 UTC vs.
Revision 1.64 by root, Wed Dec 1 03:32:40 2021 UTC