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

Comparing CV/CV.xs (file contents):
Revision 1.25 by root, Fri Aug 19 00:34:54 2005 UTC vs.
Revision 1.59 by root, Sat Apr 10 03:32:09 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)
…		…
140	#endif	371	#endif
141	OUTPUT:	372	OUTPUT:
142	RETVAL	373	RETVAL
143		374
144	GdkPixbuf_noinc *	375	GdkPixbuf_noinc *
145	transpose (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;
152	guchar src = gdk_pixbuf_get_pixels (pb), dst;	384	guchar src = gdk_pixbuf_get_pixels (pb), dst;
153	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	385	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
154		386
155	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w);	387	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
156		388
157	dst = gdk_pixbuf_get_pixels (RETVAL);	389	dst = gdk_pixbuf_get_pixels (RETVAL);
158	dstr = gdk_pixbuf_get_rowstride (RETVAL);	390	dstr = gdk_pixbuf_get_rowstride (RETVAL);
159		391
160	for (y = 0; y < h; y++)
161	for (x = 0; x < w; x++)	392	for (x = 0; x < w; x++)
		393	for (y = 0; y < h; y++)
162	for (i = 0; i < bpp; i++)	394	for (i = 0; i < 3; i++)
163	dst[y * bpp + x * 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	flop (GdkPixbuf *pb)	402	rotate (GdkPixbuf *pb, int angle)
170	CODE:	403	CODE:
171	{	404	perlinterp_release ();
172	int w = gdk_pixbuf_get_width (pb);	405	if (angle < 0)
173	int h = gdk_pixbuf_get_height (pb);	406	angle += 360;
174	int bpp = gdk_pixbuf_get_n_channels (pb);	407	RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
175	int x, y, i;	408	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
176	guchar src = gdk_pixbuf_get_pixels (pb), dst;	409	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
177	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	410	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
178		411	: angle);
179	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h);	412	perlinterp_acquire ();
180
181	dst = gdk_pixbuf_get_pixels (RETVAL);
182	dstr = gdk_pixbuf_get_rowstride (RETVAL);
183
184	for (y = 0; y < h; y++)
185	for (x = 0; x < w; x++)
186	for (i = 0; i < bpp; i++)
187	dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
188	}
189	OUTPUT:	413	OUTPUT:
190	RETVAL	414	RETVAL
191		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
		457	XSRETURN_UNDEF;
		458	}
		459	OUTPUT:
		460	RETVAL
		461
192	GdkPixbuf_noinc *	462	GdkPixbuf_noinc *
193	load_jpeg (SV *path, int thumbnail=0)	463	decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		464	CODE:
		465	{
		466	#if WEBP
		467	STRLEN data_size;
		468	int alpha;
		469	WebPData data;
		470	WebPDemuxer *demux;
		471	WebPIterator iter;
		472	WebPDecoderConfig config;
		473	int inw, inh;
		474
		475	data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
		476	data.size = data_size;
		477
		478	perlinterp_release ();
		479
		480	RETVAL = 0;
		481
		482	if (!(demux = WebPDemux (&data)))
		483	goto err_demux;
		484
		485	if (!WebPDemuxGetFrame (demux, 1, &iter))
		486	goto err_iter;
		487
		488	if (!WebPInitDecoderConfig (&config))
		489	goto err_iter;
		490
		491	config.options.use_threads = 1;
		492
		493	if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
		494	goto err_iter;
		495
		496	inw = config.input.width;
		497	inh = config.input.height;
		498
		499	if (thumbnail)
		500	{
		501	if (inw * ih > inh * iw)
		502	ih = (iw * inh + inw - 1) / inw;
		503	else
		504	iw = (ih * inw + inh - 1) / inh;
		505
		506	config.options.bypass_filtering = 1;
		507	config.options.no_fancy_upsampling = 1;
		508
		509	config.options.use_scaling = 1;
		510	config.options.scaled_width = iw;
		511	config.options.scaled_height = ih;
		512	}
		513	else
		514	{
		515	iw = inw;
		516	ih = inh;
		517	}
		518
		519	alpha = !!config.input.has_alpha;
		520
		521	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
		522	if (!RETVAL)
		523	goto err_iter;
		524
		525	config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
		526	config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
		527	config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
		528	config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
		529	config.output.is_external_memory = 1;
		530
		531	if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
		532	{
		533	g_object_unref (RETVAL);
		534	RETVAL = 0;
		535	goto err_iter;
		536	}
		537
		538	err_iter:
		539	WebPDemuxReleaseIterator (&iter);
		540	err_demux:
		541	WebPDemuxDelete (demux);
		542
		543	perlinterp_acquire ();
		544	#else
		545	croak ("load_webp: webp not enabled at compile time");
		546	#endif
		547	}
		548	OUTPUT:
		549	RETVAL
		550
		551	GdkPixbuf_noinc *
		552	decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		553	CODE:
		554	{
		555	#if JXL
		556	JxlDecoder *dec = JxlDecoderCreate (0);
		557	JxlBasicInfo info;
		558	const uint8_t next_in = (uint8_t )SvPVbyte_nolen (image_data);
		559	size_t avail_in = SvCUR (image_data);
		560	const char *error = 0;
		561	void *runner = 0;
		562	struct bmff_box box;
		563
		564	perlinterp_release ();
		565
		566	RETVAL = 0;
		567
		568	error = "JxlDecoderCreate failed";
		569	if (!dec)
		570	goto done;
		571
		572	runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
		573
		574	error = "JxlDecoderSetParallelRunner failed";
		575	if (JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner) != JXL_DEC_SUCCESS)
		576	goto done;
		577
		578	error = "JxlDecoderSubscribeEvents failed";
		579	if (JxlDecoderSubscribeEvents (dec, JXL_DEC_FULL_IMAGE \| JXL_DEC_BASIC_INFO) != JXL_DEC_SUCCESS)
		580	goto done;
		581
		582	error = "JxlDecoderSetInput failed";
		583	if (JxlDecoderSetInput (dec, next_in, avail_in) != JXL_DEC_SUCCESS)
		584	goto done;
		585
		586	for (;;)
		587	{
		588	JxlDecoderStatus status = JxlDecoderProcessInput (dec);
		589
		590	printf ("status %d\n",status);
		591
		592	switch (status)
		593	{
		594	case JXL_DEC_ERROR:
		595	error = "JxlDecoderProcessInput failed";
		596	goto done;
		597
		598	case JXL_DEC_NEED_MORE_INPUT:
		599	error = "incomplete file";
		600	goto done;
		601
		602	case JXL_DEC_SUCCESS:
		603	goto done;
		604
		605	case JXL_DEC_NEED_IMAGE_OUT_BUFFER:
		606	{
		607	error = "JxlDecoderGetBasicInfo failed";
		608	if (JxlDecoderGetBasicInfo (dec, &info) != JXL_DEC_SUCCESS)
		609	goto done;
		610
		611	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
		612	error = "unable to allocate pixbuf";
		613	if (!RETVAL)
		614	goto done;
		615
		616	JxlPixelFormat format = {
		617	info.alpha_bits ? 4 : 3,
		618	JXL_TYPE_UINT8,
		619	JXL_NATIVE_ENDIAN,
		620	gdk_pixbuf_get_rowstride (RETVAL)
		621	};
		622
		623	error = "JxlDecoderSetImageOutBuffer failed";
		624	if (JxlDecoderSetImageOutBuffer (
		625	dec,
		626	&format,
		627	gdk_pixbuf_get_pixels (RETVAL),
		628	gdk_pixbuf_get_byte_length (RETVAL)
		629	) != JXL_DEC_SUCCESS)
		630	goto done;
		631	}
		632	break;
		633
		634	default:
		635	error = "unexpected event";
		636	goto done;
		637	}
		638	}
		639
		640	done:
		641	if (dec)
		642	JxlDecoderDestroy (dec);
		643
		644	if (runner)
		645	JxlThreadParallelRunnerDestroy (runner);
		646
		647	perlinterp_acquire ();
		648
		649	if (error)
		650	{
		651	if (RETVAL)
		652	g_object_unref (RETVAL);
		653
		654	croak ("load_jxl: %s", error);
		655	}
		656	#else
		657	croak ("load_jxl: jpeg-xl not enabled at compile time");
		658	#endif
		659	}
		660	OUTPUT:
		661	RETVAL
		662
		663	GdkPixbuf_noinc *
		664	decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
194	CODE:	665	CODE:
195	{	666	{
196	struct jpeg_decompress_struct cinfo;	667	struct jpeg_decompress_struct cinfo;
197	struct jpg_err_mgr jerr;	668	struct jpg_err_mgr jerr;
198	guchar *data;
199	int rs;	669	int rs;
200	FILE *fp;
201	volatile GdkPixbuf *pb = 0;	670	volatile GdkPixbuf *pb = 0;
202	gchar *filename;	671	STRLEN data_len;
		672	guchar *data = SvPVbyte (image_data, data_len);
203		673
204	RETVAL = 0;	674	RETVAL = 0;
205		675
206	filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0);	676	perlinterp_release ();
207	fp = fopen (filename, "rb");
208	g_free (filename);
209
210	if (!fp)
211	XSRETURN_UNDEF;
212		677
213	cinfo.err = jpeg_std_error (&jerr.err);	678	cinfo.err = jpeg_std_error (&jerr.err);
214		679
215	jerr.err.error_exit = cv_error_exit;	680	jerr.err.error_exit = cv_error_exit;
216	jerr.err.output_message = cv_error_output;	681	jerr.err.output_message = cv_error_output;
217		682
218	if ((rs = setjmp (jerr.setjmp_buffer)))	683	if ((rs = setjmp (jerr.setjmp_buffer)))
219	{	684	{
220	fclose (fp);
221	jpeg_destroy_decompress (&cinfo);	685	jpeg_destroy_decompress (&cinfo);
222		686
223	if (pb)	687	if (pb)
224	g_object_unref ((gpointer)pb);	688	g_object_unref ((gpointer)pb);
225		689
		690	perlinterp_acquire ();
226	XSRETURN_UNDEF;	691	XSRETURN_UNDEF;
227	}	692	}
228		693
		694	if (!data_len)
		695	longjmp (jerr.setjmp_buffer, 4);
		696
229	jpeg_create_decompress (&cinfo);	697	jpeg_create_decompress (&cinfo);
		698	cv_jpeg_mem_src (&cinfo, data, data_len);
230		699
231	jpeg_stdio_src (&cinfo, fp);
232	jpeg_read_header (&cinfo, TRUE);	700	jpeg_read_header (&cinfo, TRUE);
233		701
234	cinfo.dct_method = JDCT_DEFAULT;	702	cinfo.dct_method = JDCT_DEFAULT;
235	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */	703	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
236	cinfo.do_block_smoothing = FALSE;	704	cinfo.do_block_smoothing = FALSE;
…		…
246	{	714	{
247	cinfo.dct_method = JDCT_FASTEST;	715	cinfo.dct_method = JDCT_FASTEST;
248	cinfo.do_fancy_upsampling = FALSE;	716	cinfo.do_fancy_upsampling = FALSE;
249		717
250	while (cinfo.scale_denom < 8	718	while (cinfo.scale_denom < 8
251	&& cinfo.output_width >= IW*4	719	&& cinfo.output_width >= iw*4
252	&& cinfo.output_height >= IH*4)	720	&& cinfo.output_height >= ih*4)
253	{	721	{
254	cinfo.scale_denom <<= 1;	722	cinfo.scale_denom <<= 1;
255	jpeg_calc_output_dimensions (&cinfo);	723	jpeg_calc_output_dimensions (&cinfo);
256	}	724	}
257	}	725	}
258		726
		727	if (cinfo.output_components != 3)
		728	longjmp (jerr.setjmp_buffer, 3);
		729
		730	if (cinfo.jpeg_color_space == JCS_YCCK \|\| cinfo.jpeg_color_space == JCS_CMYK)
		731	{
		732	cinfo.out_color_space = JCS_CMYK;
		733	cinfo.output_components = 4;
		734	}
		735
259	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height);	736	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
260	if (!RETVAL)	737	if (!RETVAL)
261	longjmp (jerr.setjmp_buffer, 2);	738	longjmp (jerr.setjmp_buffer, 2);
262		739
263	data = gdk_pixbuf_get_pixels (RETVAL);	740	data = gdk_pixbuf_get_pixels (RETVAL);
264	rs = gdk_pixbuf_get_rowstride (RETVAL);	741	rs = gdk_pixbuf_get_rowstride (RETVAL);
265
266	if (cinfo.output_components != 3)
267	longjmp (jerr.setjmp_buffer, 3);
268		742
269	jpeg_start_decompress (&cinfo);	743	jpeg_start_decompress (&cinfo);
270		744
271	while (cinfo.output_scanline < cinfo.output_height)	745	while (cinfo.output_scanline < cinfo.output_height)
272	{	746	{
…		…
279	rp [3] = (guchar *)rp [2] + rs;	753	rp [3] = (guchar *)rp [2] + rs;
280		754
281	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);	755	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
282	}	756	}
283		757
		758	if (cinfo.out_color_space == JCS_CMYK)
		759	{
		760	guchar end = data + cinfo.output_height rs;
		761
		762	while (data < end)
		763	{
		764	U32 c = data [0];
		765	U32 m = data [1];
		766	U32 y = data [2];
		767	U32 k = data [3];
		768
		769	if (0)
		770	if (cinfo.Adobe_transform == 2)
		771	{
		772	c ^= 0xff;
		773	m ^= 0xff;
		774	y ^= 0xff;
		775	k ^= 0xff;
		776	}
		777
		778	data [0] = (c * k + 0x80) / 0xff;
		779	data [1] = (m * k + 0x80) / 0xff;
		780	data [2] = (y * k + 0x80) / 0xff;
		781	data [3] = 0xff;
		782
		783	data += 4;
		784	}
		785	}
		786
284	jpeg_finish_decompress (&cinfo);	787	jpeg_finish_decompress (&cinfo);
285	fclose (fp);
286	jpeg_destroy_decompress (&cinfo);	788	jpeg_destroy_decompress (&cinfo);
		789	perlinterp_acquire ();
287	}	790	}
288	OUTPUT:	791	OUTPUT:
289	RETVAL	792	RETVAL
290		793
		794	void
		795	compare (GdkPixbuf a, GdkPixbuf b)
		796	PPCODE:
		797	perlinterp_release ();
		798	{
		799	int w = gdk_pixbuf_get_width (a);
		800	int h = gdk_pixbuf_get_height (a);
		801
		802	int sa = gdk_pixbuf_get_rowstride (a);
		803	int sb = gdk_pixbuf_get_rowstride (b);
		804
		805	guchar *pa = gdk_pixbuf_get_pixels (a);
		806	guchar *pb = gdk_pixbuf_get_pixels (b);
		807
		808	int x, y;
		809
		810	assert (w == gdk_pixbuf_get_width (b));
		811	assert (h == gdk_pixbuf_get_height (b));
		812
		813	assert (gdk_pixbuf_get_n_channels (a) == 3);
		814	assert (gdk_pixbuf_get_n_channels (b) == 3);
		815
		816	double diff = 0.;
		817	int peak = 0;
		818
		819	if (w && h)
		820	for (y = 0; y < h; y++)
		821	{
		822	guchar pa_ = pa + y sa;
		823	guchar pb_ = pb + y sb;
		824
		825	for (x = 0; x < w; x++)
		826	{
		827	int d;
		828
		829	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		830	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		831	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		832	}
		833	}
		834
		835	perlinterp_acquire ();
		836
		837	EXTEND (SP, 2);
		838	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		839	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		840	}
		841
291	#############################################################################	842	#############################################################################
292		843
293	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	844	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		845
		846	# currently only works for filenames (octet strings)
294		847
295	SV *	848	SV *
296	foldcase (SV *pathsv)	849	foldcase (SV *pathsv)
297	PROTOTYPE: $	850	PROTOTYPE: $
298	CODE:	851	CODE:
299	{	852	{
300	STRLEN plen;	853	STRLEN plen;
301	U8 path = (U8 )SvPVutf8 (pathsv, plen);	854	U8 path = (U8 )SvPV (pathsv, plen);
302	U8 *pend = path + plen;	855	U8 *pend = path + plen;
303	U8 dst [plen * 6 * 3], *dstp = dst;	856	U8 dst [plen * 6 * 3], *dstp = dst;
304		857
305	while (path < pend)	858	while (path < pend)
306	{	859	{
307	U8 ch = *path;	860	U8 ch = *path;
308		861
309	if (ch >= 'a' && ch <= 'z')	862	if (ch >= 'a' && ch <= 'z')
310	dstp++ = path++;	863	dstp++ = path++;
		864	else if (ch >= 'A' && ch <= 'Z')
		865	dstp++ = path++ + ('a' - 'A');
311	else if (ch >= '0' && ch <= '9')	866	else if (ch >= '0' && ch <= '9')
312	{	867	{
313	STRLEN el, nl = 0;	868	STRLEN el, nl = 0;
314	while (path >= '0' && path <= '9' && path < pend)	869	while (path >= '0' && path <= '9' && path < pend)
315	path++, nl++;	870	path++, nl++;
…		…
318	*dstp++ = '0';	873	*dstp++ = '0';
319		874
320	memcpy (dstp, path - nl, nl);	875	memcpy (dstp, path - nl, nl);
321	dstp += nl;	876	dstp += nl;
322	}	877	}
		878	else
		879	dstp++ = path++;
		880	#if 0
323	else	881	else
324	{	882	{
325	STRLEN cl;	883	STRLEN cl;
326	to_utf8_fold (path, dstp, &cl);	884	to_utf8_fold (path, dstp, &cl);
327	dstp += cl;	885	dstp += cl;
328	path += is_utf8_char (path);	886	path += is_utf8_char (path);
329	}	887	}
		888	#endif
330	}	889	}
331		890
332	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	891	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
333	}	892	}
334	OUTPUT:	893	OUTPUT:
…		…
359	}	918	}
360	}	919	}
361	OUTPUT:	920	OUTPUT:
362	RETVAL	921	RETVAL
363		922
		923	#############################################################################
		924
		925	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
		926
		927	void
		928	dump_ascii85 (PerlIO fp, GdkPixbuf pb)
		929	CODE:
		930	{
		931	int w = gdk_pixbuf_get_width (pb);
		932	int h = gdk_pixbuf_get_height (pb);
		933	int x, y, i;
		934	guchar *dst;
		935	int bpp = gdk_pixbuf_get_n_channels (pb);
		936	guchar *src = gdk_pixbuf_get_pixels (pb);
		937	int sstr = gdk_pixbuf_get_rowstride (pb);
		938
		939	a85_init ();
		940
		941	for (y = 0; y < h; y++)
		942	for (x = 0; x < w; x++)
		943	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
		944	a85_push (fp, src [x * bpp + y * sstr + i]);
		945
		946	a85_finish (fp);
		947	}
		948
		949	void
		950	dump_binary (PerlIO fp, GdkPixbuf pb)
		951	CODE:
		952	{
		953	int w = gdk_pixbuf_get_width (pb);
		954	int h = gdk_pixbuf_get_height (pb);
		955	int x, y, i;
		956	guchar *dst;
		957	int bpp = gdk_pixbuf_get_n_channels (pb);
		958	guchar *src = gdk_pixbuf_get_pixels (pb);
		959	int sstr = gdk_pixbuf_get_rowstride (pb);
		960
		961	for (y = 0; y < h; y++)
		962	for (x = 0; x < w; x++)
		963	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
		964	PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
		965	}
		966
		967	#############################################################################
		968
		969	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
		970
364	SV *	971	SV *
365	pb_to_p7 (GdkPixbuf *pb)	972	pb_to_hv84 (GdkPixbuf *pb)
366	CODE:	973	CODE:
367	{	974	{
368	int w = gdk_pixbuf_get_width (pb);	975	int w = gdk_pixbuf_get_width (pb);
369	int h = gdk_pixbuf_get_height (pb);	976	int h = gdk_pixbuf_get_height (pb);
370	int x, y;	977	int x, y;
371	guchar *dst;	978	guchar *dst;
372	int bpp = gdk_pixbuf_get_n_channels (pb);	979	int bpp = gdk_pixbuf_get_n_channels (pb);
373	guchar *src = gdk_pixbuf_get_pixels (pb);	980	guchar *src = gdk_pixbuf_get_pixels (pb);
374	int sstr = gdk_pixbuf_get_rowstride (pb);	981	int sstr = gdk_pixbuf_get_rowstride (pb);
375	int Er[IW], Eg[IW], Eb[IW];
376	int seed = 77;
377
378	RETVAL = newSV (w * h);
379	SvPOK_only (RETVAL);
380	SvCUR_set (RETVAL, w * h);
381
382	dst = (guchar *)SvPVX (RETVAL);
383
384	memset (Er, 0, sizeof (int) * IW);
385	memset (Eg, 0, sizeof (int) * IW);
386	memset (Eb, 0, sizeof (int) * IW);
387
388	/* some primitive error distribution + random dithering */
389
390	for (y = 0; y < h; y++)
391	{
392	int er = 0, eg = 0, eb = 0;
393
394	for (x = 0; x < w; x++)
395	{
396	int r, g, b;
397	guchar p = src + x bpp + y * sstr;
398
399	r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255;
400	g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255;
401	b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255;
402
403	r = r > 7 ? 7 : r < 0 ? 0 : r;
404	g = g > 7 ? 7 : g < 0 ? 0 : g;
405	b = b > 3 ? 3 : b < 0 ? 0 : b;
406
407	er += p[0] - (r * 255 + 4) / 7;
408	eg += p[1] - (g * 255 + 4) / 7;
409	eb += p[2] - (b * 255 + 2) / 3;
410
411	Er[x] = er >> 1; er -= (er + 1) >> 1;
412	Eg[x] = eg >> 1; eg -= (eg + 1) >> 1;
413	Eb[x] = eb >> 1; eb -= (eb + 1) >> 1;
414
415	*dst++ = r << 5 \| g << 2 \| b;
416	}
417	}
418	}
419	OUTPUT:
420	RETVAL
421
422	#############################################################################
423
424	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
425
426	void
427	dump_ascii85 (PerlIO fp, GdkPixbuf pb)
428	CODE:
429	{
430	int w = gdk_pixbuf_get_width (pb);
431	int h = gdk_pixbuf_get_height (pb);
432	int x, y, i;
433	guchar *dst;
434	int bpp = gdk_pixbuf_get_n_channels (pb);
435	guchar *src = gdk_pixbuf_get_pixels (pb);
436	int sstr = gdk_pixbuf_get_rowstride (pb);
437
438	a85_init ();
439
440	for (y = 0; y < h; y++)
441	for (x = 0; x < w; x++)
442	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
443	a85_push (fp, src [x * bpp + y * sstr + i]);
444
445	a85_finish (fp);
446	}
447
448	void
449	dump_binary (PerlIO fp, GdkPixbuf pb)
450	CODE:
451	{
452	int w = gdk_pixbuf_get_width (pb);
453	int h = gdk_pixbuf_get_height (pb);
454	int x, y, i;
455	guchar *dst;
456	int bpp = gdk_pixbuf_get_n_channels (pb);
457	guchar *src = gdk_pixbuf_get_pixels (pb);
458	int sstr = gdk_pixbuf_get_rowstride (pb);
459
460	for (y = 0; y < h; y++)
461	for (x = 0; x < w; x++)
462	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
463	PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
464	}
465
466	#############################################################################
467
468	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
469
470	SV *
471	pb_to_hv84 (GdkPixbuf *pb)
472	CODE:
473	{
474	int w = gdk_pixbuf_get_width (pb);
475	int h = gdk_pixbuf_get_height (pb);
476	int x, y;
477	guchar *dst;
478	int bpp = gdk_pixbuf_get_n_channels (pb);
479	guchar *src = gdk_pixbuf_get_pixels (pb);
480	int sstr = gdk_pixbuf_get_rowstride (pb);
481		982
482	RETVAL = newSV (6 * 8 * 12 / 8);	983	RETVAL = newSV (6 * 8 * 12 / 8);
483	SvPOK_only (RETVAL);	984	SvPOK_only (RETVAL);
484	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);	985	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
485		986
…		…
554	#############################################################################	1055	#############################################################################
555		1056
556	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	1057	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
557		1058
558	SV *	1059	SV *
559	make_histograms (SV *ar)	1060	extract_features (SV *ar)
560	CODE:	1061	CODE:
561	{	1062	{
562	int i;	1063	int i;
563	AV av, result;	1064	AV av, result;
564		1065
565	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	1066	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
566	croak ("Not an array ref as first argument to make_histogram");	1067	croak ("Not an array ref as first argument to extract_features");
567		1068
568	av = (AV *) SvRV (ar);	1069	av = (AV *) SvRV (ar);
569	result = newAV ();	1070	result = newAV ();
570		1071
571	for (i = 0; i <= av_len (av); ++i)	1072	for (i = 0; i <= av_len (av); ++i)
572	{	1073	{
573	const int HISTSIZE = 64;
574
575	int j;
576	SV sv = av_fetch (av, i, 1);	1074	SV sv = av_fetch (av, i, 1);
577	STRLEN len;
578	char *buf = SvPVbyte (sv, len);
579
580	int tmphist[HISTSIZE];
581	float *hist;
582
583	SV histsv = newSV (HISTSIZE sizeof (float) + 1);	1075	SV histsv = newSV (9 sizeof (float) + 1);
		1076
584	SvPOK_on (histsv);	1077	SvPOK_on (histsv);
585	SvCUR_set (histsv, HISTSIZE * sizeof (float));	1078	SvCUR_set (histsv, 9 * sizeof (float));
586	hist = (float *)SvPVX (histsv);	1079	float hist = (float )SvPVX (histsv);
587		1080
588	Zero (tmphist, sizeof (tmphist), char);	1081	struct feature f_h, f_s, f_v;
		1082	feature_init (&f_h);
		1083	feature_init (&f_s);
		1084	feature_init (&f_v);
589		1085
590	for (j = len; j--; )
591	{
592	unsigned int idx
593	= ((*buf & 0xc0) >> 2)
594	\| ((*buf & 0x18) >> 1)
595	\| (*buf & 0x03);
596
597	++tmphist[idx];
598	++buf;
599	}
600		1086	{
601	for (j = 0; j < HISTSIZE; ++j)	1087	STRLEN len;
602	hist[j] = (float)tmphist[j] / (len + 1e-30);	1088	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1089	while (len >= 3)
		1090	{
		1091	unsigned int r, g, b, h, s, v;
		1092	r = buf++; g = buf++; b = *buf++;
		1093	rgb_to_hsv (r, g, b, &h, &s, &v);
		1094
		1095	feature_update_pass_1 (&f_h, h);
		1096	feature_update_pass_1 (&f_s, s);
		1097	feature_update_pass_1 (&f_v, v);
		1098
		1099	len -= 3;
		1100	}
		1101
		1102	feature_finish_pass_1 (&f_h);
		1103	feature_finish_pass_1 (&f_s);
		1104	feature_finish_pass_1 (&f_v);
		1105	}
		1106
		1107	{
		1108	STRLEN len;
		1109	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1110	while (len >= 3)
		1111	{
		1112	unsigned int r, g, b, h, s, v;
		1113	r = buf++; g = buf++; b = *buf++;
		1114	rgb_to_hsv (r, g, b, &h, &s, &v);
		1115
		1116	feature_update_pass_2 (&f_h, h);
		1117	feature_update_pass_2 (&f_s, s);
		1118	feature_update_pass_2 (&f_v, v);
		1119
		1120	len -= 3;
		1121	}
		1122
		1123	feature_finish_pass_2 (&f_h);
		1124	feature_finish_pass_2 (&f_s);
		1125	feature_finish_pass_2 (&f_v);
		1126	}
		1127
		1128	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		1129	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		1130	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
603		1131
604	av_push (result, histsv);	1132	av_push (result, histsv);
605	}	1133	}
606		1134
607	RETVAL = newRV_noinc ((SV *)result);	1135	RETVAL = newRV_noinc ((SV *)result);
608	}	1136	}
609	OUTPUT:	1137	OUTPUT:
610	RETVAL	1138	RETVAL
611		1139
612

Diff Legend

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

Comparing CV/CV.xs (file contents): Revision 1.25 by root, Fri Aug 19 00:34:54 2005 UTC vs. Revision 1.59 by root, Sat Apr 10 03:32:09 2021 UTC

Diff Legend

Comparing CV/CV.xs (file contents):
Revision 1.25 by root, Fri Aug 19 00:34:54 2005 UTC vs.
Revision 1.59 by root, Sat Apr 10 03:32:09 2021 UTC