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

Comparing CV/CV.xs (file contents):
Revision 1.27 by root, Sun Aug 21 02:23:52 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)
…		…
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	transpose (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, h, w);	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[y * bpp + x * 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	flop (GdkPixbuf *pb)	463	decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
194	CODE:	464	CODE:
195	{	465	{
196	int w = gdk_pixbuf_get_width (pb);	466	#if WEBP
197	int h = gdk_pixbuf_get_height (pb);	467	STRLEN data_size;
198	int bpp = gdk_pixbuf_get_n_channels (pb);	468	int alpha;
		469	WebPData data;
		470	WebPDemuxer *demux;
		471	WebPIterator iter;
		472	WebPDecoderConfig config;
199	int x, y, i;	473	int inw, inh;
200	guchar src = gdk_pixbuf_get_pixels (pb), dst;
201	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
202		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
203	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h);	521	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
		522	if (!RETVAL)
		523	goto err_iter;
204		524
205	dst = gdk_pixbuf_get_pixels (RETVAL);	525	config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
		526	config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
206	dstr = gdk_pixbuf_get_rowstride (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;
207		530
208	for (y = 0; y < h; y++)	531	if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
209	for (x = 0; x < w; x++)	532	{
210	for (i = 0; i < bpp; i++)	533	g_object_unref (RETVAL);
211	dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];	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
212	}	547	}
213	OUTPUT:	548	OUTPUT:
214	RETVAL	549	RETVAL
215		550
216	GdkPixbuf_noinc *	551	GdkPixbuf_noinc *
217	load_jpeg (SV *path, int thumbnail=0)	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)
218	CODE:	665	CODE:
219	{	666	{
220	struct jpeg_decompress_struct cinfo;	667	struct jpeg_decompress_struct cinfo;
221	struct jpg_err_mgr jerr;	668	struct jpg_err_mgr jerr;
222	guchar *data;
223	int rs;	669	int rs;
224	FILE *fp;
225	volatile GdkPixbuf *pb = 0;	670	volatile GdkPixbuf *pb = 0;
226	gchar *filename;	671	STRLEN data_len;
		672	guchar *data = SvPVbyte (image_data, data_len);
227		673
228	RETVAL = 0;	674	RETVAL = 0;
229		675
230	filename = g_filename_from_utf8 (SvPVutf8_nolen (path), -1, 0, 0, 0);	676	perlinterp_release ();
231	fp = fopen (filename, "rb");
232	g_free (filename);
233
234	if (!fp)
235	XSRETURN_UNDEF;
236		677
237	cinfo.err = jpeg_std_error (&jerr.err);	678	cinfo.err = jpeg_std_error (&jerr.err);
238		679
239	jerr.err.error_exit = cv_error_exit;	680	jerr.err.error_exit = cv_error_exit;
240	jerr.err.output_message = cv_error_output;	681	jerr.err.output_message = cv_error_output;
241		682
242	if ((rs = setjmp (jerr.setjmp_buffer)))	683	if ((rs = setjmp (jerr.setjmp_buffer)))
243	{	684	{
244	fclose (fp);
245	jpeg_destroy_decompress (&cinfo);	685	jpeg_destroy_decompress (&cinfo);
246		686
247	if (pb)	687	if (pb)
248	g_object_unref ((gpointer)pb);	688	g_object_unref ((gpointer)pb);
249		689
		690	perlinterp_acquire ();
250	XSRETURN_UNDEF;	691	XSRETURN_UNDEF;
251	}	692	}
252		693
		694	if (!data_len)
		695	longjmp (jerr.setjmp_buffer, 4);
		696
253	jpeg_create_decompress (&cinfo);	697	jpeg_create_decompress (&cinfo);
		698	cv_jpeg_mem_src (&cinfo, data, data_len);
254		699
255	jpeg_stdio_src (&cinfo, fp);
256	jpeg_read_header (&cinfo, TRUE);	700	jpeg_read_header (&cinfo, TRUE);
257		701
258	cinfo.dct_method = JDCT_DEFAULT;	702	cinfo.dct_method = JDCT_DEFAULT;
259	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 */
260	cinfo.do_block_smoothing = FALSE;	704	cinfo.do_block_smoothing = FALSE;
…		…
270	{	714	{
271	cinfo.dct_method = JDCT_FASTEST;	715	cinfo.dct_method = JDCT_FASTEST;
272	cinfo.do_fancy_upsampling = FALSE;	716	cinfo.do_fancy_upsampling = FALSE;
273		717
274	while (cinfo.scale_denom < 8	718	while (cinfo.scale_denom < 8
275	&& cinfo.output_width >= IW*4	719	&& cinfo.output_width >= iw*4
276	&& cinfo.output_height >= IH*4)	720	&& cinfo.output_height >= ih*4)
277	{	721	{
278	cinfo.scale_denom <<= 1;	722	cinfo.scale_denom <<= 1;
279	jpeg_calc_output_dimensions (&cinfo);	723	jpeg_calc_output_dimensions (&cinfo);
280	}	724	}
281	}	725	}
282		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
283	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);
284	if (!RETVAL)	737	if (!RETVAL)
285	longjmp (jerr.setjmp_buffer, 2);	738	longjmp (jerr.setjmp_buffer, 2);
286		739
287	data = gdk_pixbuf_get_pixels (RETVAL);	740	data = gdk_pixbuf_get_pixels (RETVAL);
288	rs = gdk_pixbuf_get_rowstride (RETVAL);	741	rs = gdk_pixbuf_get_rowstride (RETVAL);
289
290	if (cinfo.output_components != 3)
291	longjmp (jerr.setjmp_buffer, 3);
292		742
293	jpeg_start_decompress (&cinfo);	743	jpeg_start_decompress (&cinfo);
294		744
295	while (cinfo.output_scanline < cinfo.output_height)	745	while (cinfo.output_scanline < cinfo.output_height)
296	{	746	{
…		…
303	rp [3] = (guchar *)rp [2] + rs;	753	rp [3] = (guchar *)rp [2] + rs;
304		754
305	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);	755	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
306	}	756	}
307		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
308	jpeg_finish_decompress (&cinfo);	787	jpeg_finish_decompress (&cinfo);
309	fclose (fp);
310	jpeg_destroy_decompress (&cinfo);	788	jpeg_destroy_decompress (&cinfo);
		789	perlinterp_acquire ();
311	}	790	}
312	OUTPUT:	791	OUTPUT:
313	RETVAL	792	RETVAL
314		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
315	#############################################################################	842	#############################################################################
316		843
317	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	844	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		845
		846	# currently only works for filenames (octet strings)
318		847
319	SV *	848	SV *
320	foldcase (SV *pathsv)	849	foldcase (SV *pathsv)
321	PROTOTYPE: $	850	PROTOTYPE: $
322	CODE:	851	CODE:
323	{	852	{
324	STRLEN plen;	853	STRLEN plen;
325	U8 path = (U8 )SvPVutf8 (pathsv, plen);	854	U8 path = (U8 )SvPV (pathsv, plen);
326	U8 *pend = path + plen;	855	U8 *pend = path + plen;
327	U8 dst [plen * 6 * 3], *dstp = dst;	856	U8 dst [plen * 6 * 3], *dstp = dst;
328		857
329	while (path < pend)	858	while (path < pend)
330	{	859	{
331	U8 ch = *path;	860	U8 ch = *path;
332		861
333	if (ch >= 'a' && ch <= 'z')	862	if (ch >= 'a' && ch <= 'z')
334	dstp++ = path++;	863	dstp++ = path++;
		864	else if (ch >= 'A' && ch <= 'Z')
		865	dstp++ = path++ + ('a' - 'A');
335	else if (ch >= '0' && ch <= '9')	866	else if (ch >= '0' && ch <= '9')
336	{	867	{
337	STRLEN el, nl = 0;	868	STRLEN el, nl = 0;
338	while (path >= '0' && path <= '9' && path < pend)	869	while (path >= '0' && path <= '9' && path < pend)
339	path++, nl++;	870	path++, nl++;
…		…
342	*dstp++ = '0';	873	*dstp++ = '0';
343		874
344	memcpy (dstp, path - nl, nl);	875	memcpy (dstp, path - nl, nl);
345	dstp += nl;	876	dstp += nl;
346	}	877	}
		878	else
		879	dstp++ = path++;
		880	#if 0
347	else	881	else
348	{	882	{
349	STRLEN cl;	883	STRLEN cl;
350	to_utf8_fold (path, dstp, &cl);	884	to_utf8_fold (path, dstp, &cl);
351	dstp += cl;	885	dstp += cl;
352	path += is_utf8_char (path);	886	path += is_utf8_char (path);
353	}	887	}
		888	#endif
354	}	889	}
355		890
356	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	891	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
357	}	892	}
358	OUTPUT:	893	OUTPUT:
…		…
383	}	918	}
384	}	919	}
385	OUTPUT:	920	OUTPUT:
386	RETVAL	921	RETVAL
387		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
388	SV *	971	SV *
389	pb_to_p7 (GdkPixbuf *pb)	972	pb_to_hv84 (GdkPixbuf *pb)
390	CODE:	973	CODE:
391	{	974	{
392	int w = gdk_pixbuf_get_width (pb);	975	int w = gdk_pixbuf_get_width (pb);
393	int h = gdk_pixbuf_get_height (pb);	976	int h = gdk_pixbuf_get_height (pb);
394	int x, y;	977	int x, y;
395	guchar *dst;	978	guchar *dst;
396	int bpp = gdk_pixbuf_get_n_channels (pb);	979	int bpp = gdk_pixbuf_get_n_channels (pb);
397	guchar *src = gdk_pixbuf_get_pixels (pb);	980	guchar *src = gdk_pixbuf_get_pixels (pb);
398	int sstr = gdk_pixbuf_get_rowstride (pb);	981	int sstr = gdk_pixbuf_get_rowstride (pb);
399	int Er[IW], Eg[IW], Eb[IW];
400	int seed = 77;
401
402	RETVAL = newSV (w * h);
403	SvPOK_only (RETVAL);
404	SvCUR_set (RETVAL, w * h);
405
406	dst = (guchar *)SvPVX (RETVAL);
407
408	memset (Er, 0, sizeof (int) * IW);
409	memset (Eg, 0, sizeof (int) * IW);
410	memset (Eb, 0, sizeof (int) * IW);
411
412	/* some primitive error distribution + random dithering */
413
414	for (y = 0; y < h; y++)
415	{
416	int er = 0, eg = 0, eb = 0;
417
418	for (x = 0; x < w; x++)
419	{
420	int r, g, b;
421	guchar p = src + x bpp + y * sstr;
422
423	r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255;
424	g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255;
425	b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255;
426
427	r = r > 7 ? 7 : r < 0 ? 0 : r;
428	g = g > 7 ? 7 : g < 0 ? 0 : g;
429	b = b > 3 ? 3 : b < 0 ? 0 : b;
430
431	er += p[0] - (r * 255 + 4) / 7;
432	eg += p[1] - (g * 255 + 4) / 7;
433	eb += p[2] - (b * 255 + 2) / 3;
434
435	Er[x] = er >> 1; er -= (er + 1) >> 1;
436	Eg[x] = eg >> 1; eg -= (eg + 1) >> 1;
437	Eb[x] = eb >> 1; eb -= (eb + 1) >> 1;
438
439	*dst++ = r << 5 \| g << 2 \| b;
440	}
441	}
442	}
443	OUTPUT:
444	RETVAL
445
446	#############################################################################
447
448	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
449
450	void
451	dump_ascii85 (PerlIO fp, GdkPixbuf pb)
452	CODE:
453	{
454	int w = gdk_pixbuf_get_width (pb);
455	int h = gdk_pixbuf_get_height (pb);
456	int x, y, i;
457	guchar *dst;
458	int bpp = gdk_pixbuf_get_n_channels (pb);
459	guchar *src = gdk_pixbuf_get_pixels (pb);
460	int sstr = gdk_pixbuf_get_rowstride (pb);
461
462	a85_init ();
463
464	for (y = 0; y < h; y++)
465	for (x = 0; x < w; x++)
466	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
467	a85_push (fp, src [x * bpp + y * sstr + i]);
468
469	a85_finish (fp);
470	}
471
472	void
473	dump_binary (PerlIO fp, GdkPixbuf pb)
474	CODE:
475	{
476	int w = gdk_pixbuf_get_width (pb);
477	int h = gdk_pixbuf_get_height (pb);
478	int x, y, i;
479	guchar *dst;
480	int bpp = gdk_pixbuf_get_n_channels (pb);
481	guchar *src = gdk_pixbuf_get_pixels (pb);
482	int sstr = gdk_pixbuf_get_rowstride (pb);
483
484	for (y = 0; y < h; y++)
485	for (x = 0; x < w; x++)
486	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
487	PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
488	}
489
490	#############################################################################
491
492	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
493
494	SV *
495	pb_to_hv84 (GdkPixbuf *pb)
496	CODE:
497	{
498	int w = gdk_pixbuf_get_width (pb);
499	int h = gdk_pixbuf_get_height (pb);
500	int x, y;
501	guchar *dst;
502	int bpp = gdk_pixbuf_get_n_channels (pb);
503	guchar *src = gdk_pixbuf_get_pixels (pb);
504	int sstr = gdk_pixbuf_get_rowstride (pb);
505		982
506	RETVAL = newSV (6 * 8 * 12 / 8);	983	RETVAL = newSV (6 * 8 * 12 / 8);
507	SvPOK_only (RETVAL);	984	SvPOK_only (RETVAL);
508	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);	985	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
509		986
…		…
578	#############################################################################	1055	#############################################################################
579		1056
580	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	1057	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
581		1058
582	SV *	1059	SV *
583	make_histograms (SV *ar)	1060	extract_features (SV *ar)
584	CODE:	1061	CODE:
585	{	1062	{
586	int i;	1063	int i;
587	AV av, result;	1064	AV av, result;
588		1065
589	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	1066	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
590	croak ("Not an array ref as first argument to make_histogram");	1067	croak ("Not an array ref as first argument to extract_features");
591		1068
592	av = (AV *) SvRV (ar);	1069	av = (AV *) SvRV (ar);
593	result = newAV ();	1070	result = newAV ();
594		1071
595	for (i = 0; i <= av_len (av); ++i)	1072	for (i = 0; i <= av_len (av); ++i)
596	{	1073	{
597	const int HISTSIZE = 64;
598
599	int j;
600	SV sv = av_fetch (av, i, 1);	1074	SV sv = av_fetch (av, i, 1);
601	STRLEN len;
602	char *buf = SvPVbyte (sv, len);
603
604	int tmphist[HISTSIZE];
605	float *hist;
606
607	SV histsv = newSV (HISTSIZE sizeof (float) + 1);	1075	SV histsv = newSV (9 sizeof (float) + 1);
		1076
608	SvPOK_on (histsv);	1077	SvPOK_on (histsv);
609	SvCUR_set (histsv, HISTSIZE * sizeof (float));	1078	SvCUR_set (histsv, 9 * sizeof (float));
610	hist = (float *)SvPVX (histsv);	1079	float hist = (float )SvPVX (histsv);
611		1080
612	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);
613		1085
614	for (j = len; j--; )
615	{
616	unsigned int idx
617	= ((*buf & 0xc0) >> 2)
618	\| ((*buf & 0x18) >> 1)
619	\| (*buf & 0x03);
620
621	++tmphist[idx];
622	++buf;
623	}
624		1086	{
625	for (j = 0; j < HISTSIZE; ++j)	1087	STRLEN len;
626	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;
627		1131
628	av_push (result, histsv);	1132	av_push (result, histsv);
629	}	1133	}
630		1134
631	RETVAL = newRV_noinc ((SV *)result);	1135	RETVAL = newRV_noinc ((SV *)result);
632	}	1136	}
633	OUTPUT:	1137	OUTPUT:
634	RETVAL	1138	RETVAL
635		1139
636

Diff Legend

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

Comparing CV/CV.xs (file contents): Revision 1.27 by root, Sun Aug 21 02:23:52 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.27 by root, Sun Aug 21 02:23:52 2005 UTC vs.
Revision 1.59 by root, Sat Apr 10 03:32:09 2021 UTC