[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.61 by root, Sun Jul 25 11:24:48 2021 UTC

…		…
2	#include "perl.h"	2	#include "perl.h"
3	#include "XSUB.h"	3	#include "XSUB.h"
4		4
5	#include <string.h>	5	#include <string.h>
6	#include <setjmp.h>	6	#include <setjmp.h>
		7	#include <math.h>
		8
		9	#include <magic.h>
7		10
8	#include <jpeglib.h>	11	#include <jpeglib.h>
		12	#include <jerror.h>
		13
9	#include <glib.h>	14	#include <glib.h>
10	#include <gtk/gtk.h>	15	#include <gtk/gtk.h>
		16	#include <gdk/gdkx.h>
11	#include <gdk-pixbuf/gdk-pixbuf.h>	17	#include <gdk-pixbuf/gdk-pixbuf.h>
12		18
13	#include <gperl.h>	19	#include <gperl.h>
14	#include <gtk2perl.h>	20	#include <gtk2perl.h>
15		21
		22	#include <assert.h>
		23
		24	#if WEBP
		25	#include <webp/demux.h>
		26	#include <webp/decode.h>
		27	#endif
		28
		29	#if JXL
		30	#include <jxl/decode.h>
		31	#include "jxl/thread_parallel_runner.h"
		32	#endif
		33
		34	#include "perlmulticore.h"
		35
16	#define IW 80 /* MUST match Schnauer.pm! */	36	#define IW 80 /* MUST match Schnauzer.pm! */
17	#define IH 60 /* MUST match Schnauer.pm! */	37	#define IH 60 /* MUST match Schnauzer.pm! */
18
19	#define RAND (seed = (seed + 7141) * 54773 % 134456)
20		38
21	#define LINELENGTH 240	39	#define LINELENGTH 240
22		40
23	#define ELLIPSIS "\xe2\x80\xa6"	41	#define ELLIPSIS "\xe2\x80\xa6"
		42
		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,
		424	0x0d, 0xa, 0x87, 0x0a, 0, 0, 0, 0x14,
		425	0x66, 0x74, 0x79, 0x70, 0x6a, 0x78, 0x6c, 0x20,
		426	0, 0, 0, 0, 0x6a, 0x78, 0x6c, 0x20,
		427	};
		428
		429	if (data_len >= 20
		430	&& data[0] == 0xff
		431	&& data[1] == 0xd8
		432	&& data[2] == 0xff)
		433	RETVAL = "image/jpeg";
		434	else if (data_len >= 12
		435	&& data[ 0] == (U8)'R'
		436	&& data[ 1] == (U8)'I'
		437	&& data[ 2] == (U8)'F'
		438	&& data[ 3] == (U8)'F'
		439	&& data[ 8] == (U8)'W'
		440	&& data[ 9] == (U8)'E'
		441	&& data[10] == (U8)'B'
		442	&& data[11] == (U8)'P')
		443	RETVAL = "image/webp";
		444	else if (data_len >= 16
		445	&& data[ 0] == 0x89
		446	&& data[ 1] == (U8)'P'
		447	&& data[ 2] == (U8)'N'
		448	&& data[ 3] == (U8)'G'
		449	&& data[ 4] == 0x0d
		450	&& data[ 5] == 0x0a
		451	&& data[ 6] == 0x1a
		452	&& data[ 7] == 0x0a)
		453	RETVAL = "image/png";
		454	else if (data_len >= sizeof (jxl_header) && memcmp (data, jxl_header, sizeof (jxl_header)) == 0)
		455	RETVAL = "image/jxl";
		456	else if (data_len >= 13
		457	&& data[0] == 'G'
		458	&& data[1] == 'I'
		459	&& data[2] == 'F'
		460	&& data[3] == '8'
		461	//&& (data[4] == '7' \|\| data[4] == '9')
		462	&& data[5] == 'a')
		463	{
		464	RETVAL = "image/gif";
		465
		466	// now see if its animated - we require the netscape application header for this
		467	int ofs = 13;
		468
		469	if (data[10] & 0x80)
		470	ofs += (1 << ((data[10] & 7) + 1)) * 3;
		471
		472	if (data_len >= ofs + 2 + 1 + 11)
		473	{
		474	// skip a rare first graphic control extension block
		475	// we assume there is exactly one block
		476	if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9)
		477	ofs += 3 + data[ofs + 2] + 1;
		478
		479	if (data_len >= ofs + 2 + 1 + 11)
		480	if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1))
		481	RETVAL = "video/gif";
		482	}
		483	}
		484
		485	else
		486	XSRETURN_UNDEF;
		487	}
		488	OUTPUT:
		489	RETVAL
		490
179	GdkPixbuf_noinc *	491	GdkPixbuf_noinc *
180	load_jpeg (SV *path, int thumbnail=0)	492	decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		493	CODE:
		494	{
		495	#if WEBP
		496	STRLEN data_size;
		497	int alpha;
		498	WebPData data;
		499	WebPDemuxer *demux;
		500	WebPIterator iter;
		501	WebPDecoderConfig config;
		502	int inw, inh;
		503
		504	data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
		505	data.size = data_size;
		506
		507	perlinterp_release ();
		508
		509	RETVAL = 0;
		510
		511	if (!(demux = WebPDemux (&data)))
		512	goto err_demux;
		513
		514	if (!WebPDemuxGetFrame (demux, 1, &iter))
		515	goto err_iter;
		516
		517	if (!WebPInitDecoderConfig (&config))
		518	goto err_iter;
		519
		520	config.options.use_threads = 1;
		521
		522	if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
		523	goto err_iter;
		524
		525	inw = config.input.width;
		526	inh = config.input.height;
		527
		528	if (thumbnail)
		529	{
		530	if (inw * ih > inh * iw)
		531	ih = (iw * inh + inw - 1) / inw;
		532	else
		533	iw = (ih * inw + inh - 1) / inh;
		534
		535	config.options.bypass_filtering = 1;
		536	config.options.no_fancy_upsampling = 1;
		537
		538	config.options.use_scaling = 1;
		539	config.options.scaled_width = iw;
		540	config.options.scaled_height = ih;
		541	}
		542	else
		543	{
		544	iw = inw;
		545	ih = inh;
		546	}
		547
		548	alpha = !!config.input.has_alpha;
		549
		550	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
		551	if (!RETVAL)
		552	goto err_iter;
		553
		554	config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
		555	config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
		556	config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
		557	config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
		558	config.output.is_external_memory = 1;
		559
		560	if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
		561	{
		562	g_object_unref (RETVAL);
		563	RETVAL = 0;
		564	goto err_iter;
		565	}
		566
		567	err_iter:
		568	WebPDemuxReleaseIterator (&iter);
		569	err_demux:
		570	WebPDemuxDelete (demux);
		571
		572	perlinterp_acquire ();
		573	#else
		574	croak ("load_webp: webp not enabled at compile time");
		575	#endif
		576	}
		577	OUTPUT:
		578	RETVAL
		579
		580	GdkPixbuf_noinc *
		581	decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		582	CODE:
		583	{
		584	#if JXL
		585	JxlDecoder *dec = JxlDecoderCreate (0);
		586	JxlBasicInfo info;
		587	const uint8_t next_in = (uint8_t )SvPVbyte_nolen (image_data);
		588	size_t avail_in = SvCUR (image_data);
		589	const char *error = 0;
		590	void *runner = 0;
		591	struct bmff_box box;
		592
		593	perlinterp_release ();
		594
		595	RETVAL = 0;
		596
		597	error = "JxlDecoderCreate failed";
		598	if (!dec)
		599	goto done;
		600
		601	runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
		602
		603	error = "JxlDecoderSetParallelRunner failed";
		604	if (JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner) != JXL_DEC_SUCCESS)
		605	goto done;
		606
		607	error = "JxlDecoderSubscribeEvents failed";
		608	if (JxlDecoderSubscribeEvents (dec, JXL_DEC_FULL_IMAGE \| JXL_DEC_BASIC_INFO) != JXL_DEC_SUCCESS)
		609	goto done;
		610
		611	error = "JxlDecoderSetInput failed";
		612	if (JxlDecoderSetInput (dec, next_in, avail_in) != JXL_DEC_SUCCESS)
		613	goto done;
		614
		615	for (;;)
		616	{
		617	JxlDecoderStatus status = JxlDecoderProcessInput (dec);
		618
		619	printf ("status %d\n",status);
		620
		621	switch (status)
		622	{
		623	case JXL_DEC_ERROR:
		624	error = "JxlDecoderProcessInput failed";
		625	goto done;
		626
		627	case JXL_DEC_NEED_MORE_INPUT:
		628	error = "incomplete file";
		629	goto done;
		630
		631	case JXL_DEC_SUCCESS:
		632	goto done;
		633
		634	case JXL_DEC_NEED_IMAGE_OUT_BUFFER:
		635	{
		636	error = "JxlDecoderGetBasicInfo failed";
		637	if (JxlDecoderGetBasicInfo (dec, &info) != JXL_DEC_SUCCESS)
		638	goto done;
		639
		640	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
		641	error = "unable to allocate pixbuf";
		642	if (!RETVAL)
		643	goto done;
		644
		645	JxlPixelFormat format = {
		646	info.alpha_bits ? 4 : 3,
		647	JXL_TYPE_UINT8,
		648	JXL_NATIVE_ENDIAN,
		649	gdk_pixbuf_get_rowstride (RETVAL)
		650	};
		651
		652	error = "JxlDecoderSetImageOutBuffer failed";
		653	if (JxlDecoderSetImageOutBuffer (
		654	dec,
		655	&format,
		656	gdk_pixbuf_get_pixels (RETVAL),
		657	gdk_pixbuf_get_byte_length (RETVAL)
		658	) != JXL_DEC_SUCCESS)
		659	goto done;
		660	}
		661	break;
		662
		663	default:
		664	error = "unexpected event";
		665	goto done;
		666	}
		667	}
		668
		669	done:
		670	if (dec)
		671	JxlDecoderDestroy (dec);
		672
		673	if (runner)
		674	JxlThreadParallelRunnerDestroy (runner);
		675
		676	perlinterp_acquire ();
		677
		678	if (error)
		679	{
		680	if (RETVAL)
		681	g_object_unref (RETVAL);
		682
		683	croak ("load_jxl: %s", error);
		684	}
		685	#else
		686	croak ("load_jxl: jpeg-xl not enabled at compile time");
		687	#endif
		688	}
		689	OUTPUT:
		690	RETVAL
		691
		692	GdkPixbuf_noinc *
		693	decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
181	CODE:	694	CODE:
182	{	695	{
183	struct jpeg_decompress_struct cinfo;	696	struct jpeg_decompress_struct cinfo;
184	struct jpg_err_mgr jerr;	697	struct jpg_err_mgr jerr;
185	guchar *data;
186	int rs;	698	int rs;
187	FILE *fp;
188	volatile GdkPixbuf *pb = 0;	699	volatile GdkPixbuf *pb = 0;
189	gchar *filename;	700	STRLEN data_len;
		701	guchar *data = SvPVbyte (image_data, data_len);
190		702
191	RETVAL = 0;	703	RETVAL = 0;
192		704
193	filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0);	705	perlinterp_release ();
194	fp = fopen (filename, "rb");
195	g_free (filename);
196
197	if (!fp)
198	XSRETURN_UNDEF;
199		706
200	cinfo.err = jpeg_std_error (&jerr.err);	707	cinfo.err = jpeg_std_error (&jerr.err);
201		708
202	jerr.err.error_exit = cv_error_exit;	709	jerr.err.error_exit = cv_error_exit;
203	jerr.err.output_message = cv_error_output;	710	jerr.err.output_message = cv_error_output;
204		711
205	if ((rs = setjmp (jerr.setjmp_buffer)))	712	if ((rs = setjmp (jerr.setjmp_buffer)))
206	{	713	{
207	fclose (fp);
208	jpeg_destroy_decompress (&cinfo);	714	jpeg_destroy_decompress (&cinfo);
209		715
210	if (pb)	716	if (pb)
211	g_object_unref ((gpointer)pb);	717	g_object_unref ((gpointer)pb);
212		718
		719	perlinterp_acquire ();
213	XSRETURN_UNDEF;	720	XSRETURN_UNDEF;
214	}	721	}
215		722
		723	if (!data_len)
		724	longjmp (jerr.setjmp_buffer, 4);
		725
216	jpeg_create_decompress (&cinfo);	726	jpeg_create_decompress (&cinfo);
		727	cv_jpeg_mem_src (&cinfo, data, data_len);
217		728
218	jpeg_stdio_src (&cinfo, fp);
219	jpeg_read_header (&cinfo, TRUE);	729	jpeg_read_header (&cinfo, TRUE);
220		730
221	cinfo.dct_method = JDCT_DEFAULT;	731	cinfo.dct_method = JDCT_DEFAULT;
222	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */	732	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
223	cinfo.do_block_smoothing = FALSE;	733	cinfo.do_block_smoothing = FALSE;
…		…
233	{	743	{
234	cinfo.dct_method = JDCT_FASTEST;	744	cinfo.dct_method = JDCT_FASTEST;
235	cinfo.do_fancy_upsampling = FALSE;	745	cinfo.do_fancy_upsampling = FALSE;
236		746
237	while (cinfo.scale_denom < 8	747	while (cinfo.scale_denom < 8
238	&& cinfo.output_width >= IW*4	748	&& cinfo.output_width >= iw*4
239	&& cinfo.output_height >= IH*4)	749	&& cinfo.output_height >= ih*4)
240	{	750	{
241	cinfo.scale_denom <<= 1;	751	cinfo.scale_denom <<= 1;
242	jpeg_calc_output_dimensions (&cinfo);	752	jpeg_calc_output_dimensions (&cinfo);
243	}	753	}
244	}	754	}
245		755
		756	if (cinfo.output_components != 3)
		757	longjmp (jerr.setjmp_buffer, 3);
		758
		759	if (cinfo.jpeg_color_space == JCS_YCCK \|\| cinfo.jpeg_color_space == JCS_CMYK)
		760	{
		761	cinfo.out_color_space = JCS_CMYK;
		762	cinfo.output_components = 4;
		763	}
		764
246	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height);	765	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
247	if (!RETVAL)	766	if (!RETVAL)
248	longjmp (jerr.setjmp_buffer, 2);	767	longjmp (jerr.setjmp_buffer, 2);
249		768
250	data = gdk_pixbuf_get_pixels (RETVAL);	769	data = gdk_pixbuf_get_pixels (RETVAL);
251	rs = gdk_pixbuf_get_rowstride (RETVAL);	770	rs = gdk_pixbuf_get_rowstride (RETVAL);
252
253	if (cinfo.output_components != 3)
254	longjmp (jerr.setjmp_buffer, 3);
255		771
256	jpeg_start_decompress (&cinfo);	772	jpeg_start_decompress (&cinfo);
257		773
258	while (cinfo.output_scanline < cinfo.output_height)	774	while (cinfo.output_scanline < cinfo.output_height)
259	{	775	{
…		…
266	rp [3] = (guchar *)rp [2] + rs;	782	rp [3] = (guchar *)rp [2] + rs;
267		783
268	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);	784	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
269	}	785	}
270		786
		787	if (cinfo.out_color_space == JCS_CMYK)
		788	{
		789	guchar end = data + cinfo.output_height rs;
		790
		791	while (data < end)
		792	{
		793	U32 c = data [0];
		794	U32 m = data [1];
		795	U32 y = data [2];
		796	U32 k = data [3];
		797
		798	if (0)
		799	if (cinfo.Adobe_transform == 2)
		800	{
		801	c ^= 0xff;
		802	m ^= 0xff;
		803	y ^= 0xff;
		804	k ^= 0xff;
		805	}
		806
		807	data [0] = (c * k + 0x80) / 0xff;
		808	data [1] = (m * k + 0x80) / 0xff;
		809	data [2] = (y * k + 0x80) / 0xff;
		810	data [3] = 0xff;
		811
		812	data += 4;
		813	}
		814	}
		815
271	jpeg_finish_decompress (&cinfo);	816	jpeg_finish_decompress (&cinfo);
272	fclose (fp);
273	jpeg_destroy_decompress (&cinfo);	817	jpeg_destroy_decompress (&cinfo);
		818	perlinterp_acquire ();
274	}	819	}
275	OUTPUT:	820	OUTPUT:
276	RETVAL	821	RETVAL
277		822
		823	void
		824	compare (GdkPixbuf a, GdkPixbuf b)
		825	PPCODE:
		826	perlinterp_release ();
		827	{
		828	int w = gdk_pixbuf_get_width (a);
		829	int h = gdk_pixbuf_get_height (a);
		830
		831	int sa = gdk_pixbuf_get_rowstride (a);
		832	int sb = gdk_pixbuf_get_rowstride (b);
		833
		834	guchar *pa = gdk_pixbuf_get_pixels (a);
		835	guchar *pb = gdk_pixbuf_get_pixels (b);
		836
		837	int x, y;
		838
		839	assert (w == gdk_pixbuf_get_width (b));
		840	assert (h == gdk_pixbuf_get_height (b));
		841
		842	assert (gdk_pixbuf_get_n_channels (a) == 3);
		843	assert (gdk_pixbuf_get_n_channels (b) == 3);
		844
		845	double diff = 0.;
		846	int peak = 0;
		847
		848	if (w && h)
		849	for (y = 0; y < h; y++)
		850	{
		851	guchar pa_ = pa + y sa;
		852	guchar pb_ = pb + y sb;
		853
		854	for (x = 0; x < w; x++)
		855	{
		856	int d;
		857
		858	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		859	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		860	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		861	}
		862	}
		863
		864	perlinterp_acquire ();
		865
		866	EXTEND (SP, 2);
		867	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		868	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		869	}
		870
278	#############################################################################	871	#############################################################################
279		872
280	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	873	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		874
		875	# currently only works for filenames (octet strings)
281		876
282	SV *	877	SV *
283	foldcase (SV *pathsv)	878	foldcase (SV *pathsv)
284	PROTOTYPE: $	879	PROTOTYPE: $
285	CODE:	880	CODE:
286	{	881	{
287	STRLEN plen;	882	STRLEN plen;
288	U8 path = (U8 )SvPVutf8 (pathsv, plen);	883	U8 path = (U8 )SvPV (pathsv, plen);
289	U8 *pend = path + plen;	884	U8 *pend = path + plen;
290	U8 dst [plen * 6 * 3], *dstp = dst;	885	U8 dst [plen * 8 * 3], *dstp = dst;
291		886
292	while (path < pend)	887	while (path < pend)
293	{	888	{
294	U8 ch = *path;	889	U8 ch = *path;
295		890
296	if (ch >= 'a' && ch <= 'z')	891	if (ch >= 'a' && ch <= 'z')
297	dstp++ = path++;	892	dstp++ = path++;
		893	else if (ch >= 'A' && ch <= 'Z')
		894	dstp++ = path++ + ('a' - 'A');
298	else if (ch >= '0' && ch <= '9')	895	else if (ch >= '0' && ch <= '9')
299	{	896	{
		897	/* version sort, up to 8 digits */
300	STRLEN el, nl = 0;	898	STRLEN el, nl = 0;
301	while (path >= '0' && path <= '9' && path < pend)	899	while (path >= '0' && path <= '9' && path < pend)
302	path++, nl++;	900	path++, nl++;
303		901
304	for (el = nl; el < 6; el++)	902	for (el = nl; el < 8; el++)
305	*dstp++ = '0';	903	*dstp++ = '0';
306		904
307	memcpy (dstp, path - nl, nl);	905	memcpy (dstp, path - nl, nl);
308	dstp += nl;	906	dstp += nl;
309	}	907	}
		908	else
		909	dstp++ = path++;
		910	#if 0
310	else	911	else
311	{	912	{
312	STRLEN cl;	913	STRLEN cl;
313	to_utf8_fold (path, dstp, &cl);	914	to_utf8_fold (path, dstp, &cl);
314	dstp += cl;	915	dstp += cl;
315	path += is_utf8_char (path);	916	path += is_utf8_char (path);
316	}	917	}
		918	#endif
317	}	919	}
318		920
319	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	921	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
320	}	922	}
321	OUTPUT:	923	OUTPUT:
…		…
483	#############################################################################	1085	#############################################################################
484		1086
485	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	1087	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
486		1088
487	SV *	1089	SV *
488	make_histograms (SV *ar)	1090	extract_features (SV *ar)
489	CODE:	1091	CODE:
490	{	1092	{
491	int i;	1093	int i;
492	AV av, result;	1094	AV av, result;
493		1095
494	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	1096	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
495	croak ("Not an array ref as first argument to make_histogram");	1097	croak ("Not an array ref as first argument to extract_features");
496		1098
497	av = (AV *) SvRV (ar);	1099	av = (AV *) SvRV (ar);
498	result = newAV ();	1100	result = newAV ();
499		1101
500	for (i = 0; i <= av_len (av); ++i)	1102	for (i = 0; i <= av_len (av); ++i)
501	{	1103	{
502	const int HISTSIZE = 64;
503
504	int j;
505	SV sv = av_fetch (av, i, 1);	1104	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);	1105	SV histsv = newSV (9 sizeof (float) + 1);
		1106
513	SvPOK_on (histsv);	1107	SvPOK_on (histsv);
514	SvCUR_set (histsv, HISTSIZE * sizeof (float));	1108	SvCUR_set (histsv, 9 * sizeof (float));
515	hist = (float *)SvPVX (histsv);	1109	float hist = (float )SvPVX (histsv);
516		1110
517	Zero (tmphist, sizeof (tmphist), char);	1111	struct feature f_h, f_s, f_v;
		1112	feature_init (&f_h);
		1113	feature_init (&f_s);
		1114	feature_init (&f_v);
518		1115
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		1116	{
530	for (j = 0; j < HISTSIZE; ++j)	1117	STRLEN len;
531	hist[j] = (float)tmphist[j] / (len + 1e-30);	1118	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1119	while (len >= 3)
		1120	{
		1121	unsigned int r, g, b, h, s, v;
		1122	r = buf++; g = buf++; b = *buf++;
		1123	rgb_to_hsv (r, g, b, &h, &s, &v);
		1124
		1125	feature_update_pass_1 (&f_h, h);
		1126	feature_update_pass_1 (&f_s, s);
		1127	feature_update_pass_1 (&f_v, v);
		1128
		1129	len -= 3;
		1130	}
		1131
		1132	feature_finish_pass_1 (&f_h);
		1133	feature_finish_pass_1 (&f_s);
		1134	feature_finish_pass_1 (&f_v);
		1135	}
		1136
		1137	{
		1138	STRLEN len;
		1139	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1140	while (len >= 3)
		1141	{
		1142	unsigned int r, g, b, h, s, v;
		1143	r = buf++; g = buf++; b = *buf++;
		1144	rgb_to_hsv (r, g, b, &h, &s, &v);
		1145
		1146	feature_update_pass_2 (&f_h, h);
		1147	feature_update_pass_2 (&f_s, s);
		1148	feature_update_pass_2 (&f_v, v);
		1149
		1150	len -= 3;
		1151	}
		1152
		1153	feature_finish_pass_2 (&f_h);
		1154	feature_finish_pass_2 (&f_s);
		1155	feature_finish_pass_2 (&f_v);
		1156	}
		1157
		1158	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		1159	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		1160	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
532		1161
533	av_push (result, histsv);	1162	av_push (result, histsv);
534	}	1163	}
535		1164
536	RETVAL = newRV_noinc ((SV *)result);	1165	RETVAL = newRV_noinc ((SV *)result);
537	}	1166	}
538	OUTPUT:	1167	OUTPUT:
539	RETVAL	1168	RETVAL
540		1169
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.61 by root, Sun Jul 25 11:24:48 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.61 by root, Sun Jul 25 11:24:48 2021 UTC