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

Comparing CV/CV.xs (file contents):
Revision 1.30 by root, Wed Dec 28 20:56:20 2005 UTC vs.
Revision 1.62 by root, Sun Jul 25 11:32:46 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
131	# missing function in perl. really :)	309	# missing function in perl. really :)
132	int	310	int
133	common_prefix_length (a, b)	311	common_prefix_length (a, b)
134	unsigned char a = (unsigned char )SvPVutf8_nolen ($arg);	312	unsigned char a = (unsigned char )SvPVutf8_nolen ($arg);
135	unsigned char b = (unsigned char )SvPVutf8_nolen ($arg);	313	unsigned char b = (unsigned char )SvPVutf8_nolen ($arg);
…		…
143	}	321	}
144		322
145	OUTPUT:	323	OUTPUT:
146	RETVAL	324	RETVAL
147		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
148	# 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);
149		363
150	gboolean	364	gboolean
151	gdk_net_wm_supports (GdkAtom property)	365	gdk_net_wm_supports (GdkAtom property)
152	CODE:	366	CODE:
153	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)	367	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
…		…
159	RETVAL	373	RETVAL
160		374
161	GdkPixbuf_noinc *	375	GdkPixbuf_noinc *
162	dealpha_expose (GdkPixbuf *pb)	376	dealpha_expose (GdkPixbuf *pb)
163	CODE:	377	CODE:
		378	perlinterp_release ();
164	{	379	{
165	int w = gdk_pixbuf_get_width (pb);	380	int w = gdk_pixbuf_get_width (pb);
166	int h = gdk_pixbuf_get_height (pb);	381	int h = gdk_pixbuf_get_height (pb);
167	int bpp = gdk_pixbuf_get_n_channels (pb);	382	int bpp = gdk_pixbuf_get_n_channels (pb);
168	int x, y, i;	383	int x, y, i;
…		…
177	for (x = 0; x < w; x++)	392	for (x = 0; x < w; x++)
178	for (y = 0; y < h; y++)	393	for (y = 0; y < h; y++)
179	for (i = 0; i < 3; i++)	394	for (i = 0; i < 3; i++)
180	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];
181	}	396	}
		397	perlinterp_acquire ();
182	OUTPUT:	398	OUTPUT:
183	RETVAL	399	RETVAL
184		400
185	GdkPixbuf_noinc *	401	GdkPixbuf_noinc *
186	rotate (GdkPixbuf *pb, int angle)	402	rotate (GdkPixbuf *pb, int angle)
187	CODE:	403	CODE:
		404	perlinterp_release ();
		405	if (angle < 0)
		406	angle += 360;
188	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
189	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE	408	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
190	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN	409	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
191	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE	410	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
192	: angle);	411	: angle);
		412	perlinterp_acquire ();
193	OUTPUT:	413	OUTPUT:
194	RETVAL	414	RETVAL
195		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
		475	// skip a graphic control extension block. we assume
		476	// there is at most one such block - while the NAB
		477	// has to come firstz, some files do not obey this
		478	if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9)
		479	ofs += 3 + data[ofs + 2] + 1;
		480
		481	if (data_len >= ofs + 2 + 1 + 11)
		482	if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1))
		483	RETVAL = "video/gif";
		484	}
		485	}
		486
		487	else
		488	XSRETURN_UNDEF;
		489	}
		490	OUTPUT:
		491	RETVAL
		492
196	GdkPixbuf_noinc *	493	GdkPixbuf_noinc *
197	load_jpeg (SV *path, int thumbnail=0)	494	decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		495	CODE:
		496	{
		497	#if WEBP
		498	STRLEN data_size;
		499	int alpha;
		500	WebPData data;
		501	WebPDemuxer *demux;
		502	WebPIterator iter;
		503	WebPDecoderConfig config;
		504	int inw, inh;
		505
		506	data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
		507	data.size = data_size;
		508
		509	perlinterp_release ();
		510
		511	RETVAL = 0;
		512
		513	if (!(demux = WebPDemux (&data)))
		514	goto err_demux;
		515
		516	if (!WebPDemuxGetFrame (demux, 1, &iter))
		517	goto err_iter;
		518
		519	if (!WebPInitDecoderConfig (&config))
		520	goto err_iter;
		521
		522	config.options.use_threads = 1;
		523
		524	if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
		525	goto err_iter;
		526
		527	inw = config.input.width;
		528	inh = config.input.height;
		529
		530	if (thumbnail)
		531	{
		532	if (inw * ih > inh * iw)
		533	ih = (iw * inh + inw - 1) / inw;
		534	else
		535	iw = (ih * inw + inh - 1) / inh;
		536
		537	config.options.bypass_filtering = 1;
		538	config.options.no_fancy_upsampling = 1;
		539
		540	config.options.use_scaling = 1;
		541	config.options.scaled_width = iw;
		542	config.options.scaled_height = ih;
		543	}
		544	else
		545	{
		546	iw = inw;
		547	ih = inh;
		548	}
		549
		550	alpha = !!config.input.has_alpha;
		551
		552	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
		553	if (!RETVAL)
		554	goto err_iter;
		555
		556	config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
		557	config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
		558	config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
		559	config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
		560	config.output.is_external_memory = 1;
		561
		562	if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
		563	{
		564	g_object_unref (RETVAL);
		565	RETVAL = 0;
		566	goto err_iter;
		567	}
		568
		569	err_iter:
		570	WebPDemuxReleaseIterator (&iter);
		571	err_demux:
		572	WebPDemuxDelete (demux);
		573
		574	perlinterp_acquire ();
		575	#else
		576	croak ("load_webp: webp not enabled at compile time");
		577	#endif
		578	}
		579	OUTPUT:
		580	RETVAL
		581
		582	GdkPixbuf_noinc *
		583	decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		584	CODE:
		585	{
		586	#if JXL
		587	JxlDecoder *dec = JxlDecoderCreate (0);
		588	JxlBasicInfo info;
		589	const uint8_t next_in = (uint8_t )SvPVbyte_nolen (image_data);
		590	size_t avail_in = SvCUR (image_data);
		591	const char *error = 0;
		592	void *runner = 0;
		593	struct bmff_box box;
		594
		595	perlinterp_release ();
		596
		597	RETVAL = 0;
		598
		599	error = "JxlDecoderCreate failed";
		600	if (!dec)
		601	goto done;
		602
		603	runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
		604
		605	error = "JxlDecoderSetParallelRunner failed";
		606	if (JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner) != JXL_DEC_SUCCESS)
		607	goto done;
		608
		609	error = "JxlDecoderSubscribeEvents failed";
		610	if (JxlDecoderSubscribeEvents (dec, JXL_DEC_FULL_IMAGE \| JXL_DEC_BASIC_INFO) != JXL_DEC_SUCCESS)
		611	goto done;
		612
		613	error = "JxlDecoderSetInput failed";
		614	if (JxlDecoderSetInput (dec, next_in, avail_in) != JXL_DEC_SUCCESS)
		615	goto done;
		616
		617	for (;;)
		618	{
		619	JxlDecoderStatus status = JxlDecoderProcessInput (dec);
		620
		621	printf ("status %d\n",status);
		622
		623	switch (status)
		624	{
		625	case JXL_DEC_ERROR:
		626	error = "JxlDecoderProcessInput failed";
		627	goto done;
		628
		629	case JXL_DEC_NEED_MORE_INPUT:
		630	error = "incomplete file";
		631	goto done;
		632
		633	case JXL_DEC_SUCCESS:
		634	goto done;
		635
		636	case JXL_DEC_NEED_IMAGE_OUT_BUFFER:
		637	{
		638	error = "JxlDecoderGetBasicInfo failed";
		639	if (JxlDecoderGetBasicInfo (dec, &info) != JXL_DEC_SUCCESS)
		640	goto done;
		641
		642	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
		643	error = "unable to allocate pixbuf";
		644	if (!RETVAL)
		645	goto done;
		646
		647	JxlPixelFormat format = {
		648	info.alpha_bits ? 4 : 3,
		649	JXL_TYPE_UINT8,
		650	JXL_NATIVE_ENDIAN,
		651	gdk_pixbuf_get_rowstride (RETVAL)
		652	};
		653
		654	error = "JxlDecoderSetImageOutBuffer failed";
		655	if (JxlDecoderSetImageOutBuffer (
		656	dec,
		657	&format,
		658	gdk_pixbuf_get_pixels (RETVAL),
		659	gdk_pixbuf_get_byte_length (RETVAL)
		660	) != JXL_DEC_SUCCESS)
		661	goto done;
		662	}
		663	break;
		664
		665	default:
		666	error = "unexpected event";
		667	goto done;
		668	}
		669	}
		670
		671	done:
		672	if (dec)
		673	JxlDecoderDestroy (dec);
		674
		675	if (runner)
		676	JxlThreadParallelRunnerDestroy (runner);
		677
		678	perlinterp_acquire ();
		679
		680	if (error)
		681	{
		682	if (RETVAL)
		683	g_object_unref (RETVAL);
		684
		685	croak ("load_jxl: %s", error);
		686	}
		687	#else
		688	croak ("load_jxl: jpeg-xl not enabled at compile time");
		689	#endif
		690	}
		691	OUTPUT:
		692	RETVAL
		693
		694	GdkPixbuf_noinc *
		695	decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
198	CODE:	696	CODE:
199	{	697	{
200	struct jpeg_decompress_struct cinfo;	698	struct jpeg_decompress_struct cinfo;
201	struct jpg_err_mgr jerr;	699	struct jpg_err_mgr jerr;
202	guchar *data;
203	int rs;	700	int rs;
204	FILE *fp;
205	volatile GdkPixbuf *pb = 0;	701	volatile GdkPixbuf *pb = 0;
		702	STRLEN data_len;
		703	guchar *data = SvPVbyte (image_data, data_len);
206		704
207	RETVAL = 0;	705	RETVAL = 0;
208		706
209	fp = fopen (SvPVbyte_nolen (path), "rb");	707	perlinterp_release ();
210
211	if (!fp)
212	XSRETURN_UNDEF;
213		708
214	cinfo.err = jpeg_std_error (&jerr.err);	709	cinfo.err = jpeg_std_error (&jerr.err);
215		710
216	jerr.err.error_exit = cv_error_exit;	711	jerr.err.error_exit = cv_error_exit;
217	jerr.err.output_message = cv_error_output;	712	jerr.err.output_message = cv_error_output;
218		713
219	if ((rs = setjmp (jerr.setjmp_buffer)))	714	if ((rs = setjmp (jerr.setjmp_buffer)))
220	{	715	{
221	fclose (fp);
222	jpeg_destroy_decompress (&cinfo);	716	jpeg_destroy_decompress (&cinfo);
223		717
224	if (pb)	718	if (pb)
225	g_object_unref ((gpointer)pb);	719	g_object_unref ((gpointer)pb);
226		720
		721	perlinterp_acquire ();
227	XSRETURN_UNDEF;	722	XSRETURN_UNDEF;
228	}	723	}
229		724
		725	if (!data_len)
		726	longjmp (jerr.setjmp_buffer, 4);
		727
230	jpeg_create_decompress (&cinfo);	728	jpeg_create_decompress (&cinfo);
		729	cv_jpeg_mem_src (&cinfo, data, data_len);
231		730
232	jpeg_stdio_src (&cinfo, fp);
233	jpeg_read_header (&cinfo, TRUE);	731	jpeg_read_header (&cinfo, TRUE);
234		732
235	cinfo.dct_method = JDCT_DEFAULT;	733	cinfo.dct_method = JDCT_DEFAULT;
236	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */	734	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
237	cinfo.do_block_smoothing = FALSE;	735	cinfo.do_block_smoothing = FALSE;
…		…
247	{	745	{
248	cinfo.dct_method = JDCT_FASTEST;	746	cinfo.dct_method = JDCT_FASTEST;
249	cinfo.do_fancy_upsampling = FALSE;	747	cinfo.do_fancy_upsampling = FALSE;
250		748
251	while (cinfo.scale_denom < 8	749	while (cinfo.scale_denom < 8
252	&& cinfo.output_width >= IW*4	750	&& cinfo.output_width >= iw*4
253	&& cinfo.output_height >= IH*4)	751	&& cinfo.output_height >= ih*4)
254	{	752	{
255	cinfo.scale_denom <<= 1;	753	cinfo.scale_denom <<= 1;
256	jpeg_calc_output_dimensions (&cinfo);	754	jpeg_calc_output_dimensions (&cinfo);
257	}	755	}
258	}	756	}
259		757
		758	if (cinfo.output_components != 3)
		759	longjmp (jerr.setjmp_buffer, 3);
		760
		761	if (cinfo.jpeg_color_space == JCS_YCCK \|\| cinfo.jpeg_color_space == JCS_CMYK)
		762	{
		763	cinfo.out_color_space = JCS_CMYK;
		764	cinfo.output_components = 4;
		765	}
		766
260	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height);	767	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
261	if (!RETVAL)	768	if (!RETVAL)
262	longjmp (jerr.setjmp_buffer, 2);	769	longjmp (jerr.setjmp_buffer, 2);
263		770
264	data = gdk_pixbuf_get_pixels (RETVAL);	771	data = gdk_pixbuf_get_pixels (RETVAL);
265	rs = gdk_pixbuf_get_rowstride (RETVAL);	772	rs = gdk_pixbuf_get_rowstride (RETVAL);
266
267	if (cinfo.output_components != 3)
268	longjmp (jerr.setjmp_buffer, 3);
269		773
270	jpeg_start_decompress (&cinfo);	774	jpeg_start_decompress (&cinfo);
271		775
272	while (cinfo.output_scanline < cinfo.output_height)	776	while (cinfo.output_scanline < cinfo.output_height)
273	{	777	{
…		…
280	rp [3] = (guchar *)rp [2] + rs;	784	rp [3] = (guchar *)rp [2] + rs;
281		785
282	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);	786	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
283	}	787	}
284		788
		789	if (cinfo.out_color_space == JCS_CMYK)
		790	{
		791	guchar end = data + cinfo.output_height rs;
		792
		793	while (data < end)
		794	{
		795	U32 c = data [0];
		796	U32 m = data [1];
		797	U32 y = data [2];
		798	U32 k = data [3];
		799
		800	if (0)
		801	if (cinfo.Adobe_transform == 2)
		802	{
		803	c ^= 0xff;
		804	m ^= 0xff;
		805	y ^= 0xff;
		806	k ^= 0xff;
		807	}
		808
		809	data [0] = (c * k + 0x80) / 0xff;
		810	data [1] = (m * k + 0x80) / 0xff;
		811	data [2] = (y * k + 0x80) / 0xff;
		812	data [3] = 0xff;
		813
		814	data += 4;
		815	}
		816	}
		817
285	jpeg_finish_decompress (&cinfo);	818	jpeg_finish_decompress (&cinfo);
286	fclose (fp);
287	jpeg_destroy_decompress (&cinfo);	819	jpeg_destroy_decompress (&cinfo);
		820	perlinterp_acquire ();
288	}	821	}
289	OUTPUT:	822	OUTPUT:
290	RETVAL	823	RETVAL
291		824
		825	void
		826	compare (GdkPixbuf a, GdkPixbuf b)
		827	PPCODE:
		828	perlinterp_release ();
		829	{
		830	int w = gdk_pixbuf_get_width (a);
		831	int h = gdk_pixbuf_get_height (a);
		832
		833	int sa = gdk_pixbuf_get_rowstride (a);
		834	int sb = gdk_pixbuf_get_rowstride (b);
		835
		836	guchar *pa = gdk_pixbuf_get_pixels (a);
		837	guchar *pb = gdk_pixbuf_get_pixels (b);
		838
		839	int x, y;
		840
		841	assert (w == gdk_pixbuf_get_width (b));
		842	assert (h == gdk_pixbuf_get_height (b));
		843
		844	assert (gdk_pixbuf_get_n_channels (a) == 3);
		845	assert (gdk_pixbuf_get_n_channels (b) == 3);
		846
		847	double diff = 0.;
		848	int peak = 0;
		849
		850	if (w && h)
		851	for (y = 0; y < h; y++)
		852	{
		853	guchar pa_ = pa + y sa;
		854	guchar pb_ = pb + y sb;
		855
		856	for (x = 0; x < w; x++)
		857	{
		858	int d;
		859
		860	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		861	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		862	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		863	}
		864	}
		865
		866	perlinterp_acquire ();
		867
		868	EXTEND (SP, 2);
		869	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		870	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		871	}
		872
292	#############################################################################	873	#############################################################################
293		874
294	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	875	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		876
		877	# currently only works for filenames (octet strings)
295		878
296	SV *	879	SV *
297	foldcase (SV *pathsv)	880	foldcase (SV *pathsv)
298	PROTOTYPE: $	881	PROTOTYPE: $
299	CODE:	882	CODE:
300	{	883	{
301	STRLEN plen;	884	STRLEN plen;
302	U8 path = (U8 )SvPVutf8 (pathsv, plen);	885	U8 path = (U8 )SvPV (pathsv, plen);
303	U8 *pend = path + plen;	886	U8 *pend = path + plen;
304	U8 dst [plen * 6 * 3], *dstp = dst;	887	U8 dst [plen * 8 * 3], *dstp = dst;
305		888
306	while (path < pend)	889	while (path < pend)
307	{	890	{
308	U8 ch = *path;	891	U8 ch = *path;
309		892
310	if (ch >= 'a' && ch <= 'z')	893	if (ch >= 'a' && ch <= 'z')
311	dstp++ = path++;	894	dstp++ = path++;
		895	else if (ch >= 'A' && ch <= 'Z')
		896	dstp++ = path++ + ('a' - 'A');
312	else if (ch >= '0' && ch <= '9')	897	else if (ch >= '0' && ch <= '9')
313	{	898	{
		899	/* version sort, up to 8 digits */
314	STRLEN el, nl = 0;	900	STRLEN el, nl = 0;
315	while (path >= '0' && path <= '9' && path < pend)	901	while (path >= '0' && path <= '9' && path < pend)
316	path++, nl++;	902	path++, nl++;
317		903
318	for (el = nl; el < 6; el++)	904	for (el = nl; el < 8; el++)
319	*dstp++ = '0';	905	*dstp++ = '0';
320		906
321	memcpy (dstp, path - nl, nl);	907	memcpy (dstp, path - nl, nl);
322	dstp += nl;	908	dstp += nl;
323	}	909	}
		910	else
		911	dstp++ = path++;
		912	#if 0
324	else	913	else
325	{	914	{
326	STRLEN cl;	915	STRLEN cl;
327	to_utf8_fold (path, dstp, &cl);	916	to_utf8_fold (path, dstp, &cl);
328	dstp += cl;	917	dstp += cl;
329	path += is_utf8_char (path);	918	path += is_utf8_char (path);
330	}	919	}
		920	#endif
331	}	921	}
332		922
333	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	923	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
334	}	924	}
335	OUTPUT:	925	OUTPUT:
…		…
497	#############################################################################	1087	#############################################################################
498		1088
499	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	1089	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
500		1090
501	SV *	1091	SV *
502	make_histograms (SV *ar)	1092	extract_features (SV *ar)
503	CODE:	1093	CODE:
504	{	1094	{
505	int i;	1095	int i;
506	AV av, result;	1096	AV av, result;
507		1097
508	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	1098	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
509	croak ("Not an array ref as first argument to make_histogram");	1099	croak ("Not an array ref as first argument to extract_features");
510		1100
511	av = (AV *) SvRV (ar);	1101	av = (AV *) SvRV (ar);
512	result = newAV ();	1102	result = newAV ();
513		1103
514	for (i = 0; i <= av_len (av); ++i)	1104	for (i = 0; i <= av_len (av); ++i)
515	{	1105	{
516	const int HISTSIZE = 64;
517
518	int j;
519	SV sv = av_fetch (av, i, 1);	1106	SV sv = av_fetch (av, i, 1);
520	STRLEN len;
521	char *buf = SvPVbyte (sv, len);
522
523	int tmphist[HISTSIZE];
524	float *hist;
525
526	SV histsv = newSV (HISTSIZE sizeof (float) + 1);	1107	SV histsv = newSV (9 sizeof (float) + 1);
		1108
527	SvPOK_on (histsv);	1109	SvPOK_on (histsv);
528	SvCUR_set (histsv, HISTSIZE * sizeof (float));	1110	SvCUR_set (histsv, 9 * sizeof (float));
529	hist = (float *)SvPVX (histsv);	1111	float hist = (float )SvPVX (histsv);
530		1112
531	Zero (tmphist, sizeof (tmphist), char);	1113	struct feature f_h, f_s, f_v;
		1114	feature_init (&f_h);
		1115	feature_init (&f_s);
		1116	feature_init (&f_v);
532		1117
533	for (j = len; j--; )
534	{
535	unsigned int idx
536	= ((*buf & 0xc0) >> 2)
537	\| ((*buf & 0x18) >> 1)
538	\| (*buf & 0x03);
539
540	++tmphist[idx];
541	++buf;
542	}
543		1118	{
544	for (j = 0; j < HISTSIZE; ++j)	1119	STRLEN len;
545	hist[j] = (float)tmphist[j] / (len + 1e-30);	1120	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1121	while (len >= 3)
		1122	{
		1123	unsigned int r, g, b, h, s, v;
		1124	r = buf++; g = buf++; b = *buf++;
		1125	rgb_to_hsv (r, g, b, &h, &s, &v);
		1126
		1127	feature_update_pass_1 (&f_h, h);
		1128	feature_update_pass_1 (&f_s, s);
		1129	feature_update_pass_1 (&f_v, v);
		1130
		1131	len -= 3;
		1132	}
		1133
		1134	feature_finish_pass_1 (&f_h);
		1135	feature_finish_pass_1 (&f_s);
		1136	feature_finish_pass_1 (&f_v);
		1137	}
		1138
		1139	{
		1140	STRLEN len;
		1141	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1142	while (len >= 3)
		1143	{
		1144	unsigned int r, g, b, h, s, v;
		1145	r = buf++; g = buf++; b = *buf++;
		1146	rgb_to_hsv (r, g, b, &h, &s, &v);
		1147
		1148	feature_update_pass_2 (&f_h, h);
		1149	feature_update_pass_2 (&f_s, s);
		1150	feature_update_pass_2 (&f_v, v);
		1151
		1152	len -= 3;
		1153	}
		1154
		1155	feature_finish_pass_2 (&f_h);
		1156	feature_finish_pass_2 (&f_s);
		1157	feature_finish_pass_2 (&f_v);
		1158	}
		1159
		1160	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		1161	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		1162	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
546		1163
547	av_push (result, histsv);	1164	av_push (result, histsv);
548	}	1165	}
549		1166
550	RETVAL = newRV_noinc ((SV *)result);	1167	RETVAL = newRV_noinc ((SV *)result);
551	}	1168	}
552	OUTPUT:	1169	OUTPUT:
553	RETVAL	1170	RETVAL
554		1171
555

Diff Legend

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

Comparing CV/CV.xs (file contents): Revision 1.30 by root, Wed Dec 28 20:56:20 2005 UTC vs. Revision 1.62 by root, Sun Jul 25 11:32:46 2021 UTC

Diff Legend

Comparing CV/CV.xs (file contents):
Revision 1.30 by root, Wed Dec 28 20:56:20 2005 UTC vs.
Revision 1.62 by root, Sun Jul 25 11:32:46 2021 UTC