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

Comparing CV/CV.xs (file contents):
Revision 1.6 by root, Wed Nov 12 00:59:01 2003 UTC vs.
Revision 1.64 by root, Wed Dec 1 03:32:40 2021 UTC

…		…
1	#include "EXTERN.h"	1	#include "EXTERN.h"
2	#include "perl.h"	2	#include "perl.h"
3	#include "XSUB.h"	3	#include "XSUB.h"
4		4
		5	#include <string.h>
		6	#include <setjmp.h>
		7	#include <math.h>
		8
		9	#include <magic.h>
		10
		11	#include <jpeglib.h>
		12	#include <jerror.h>
		13
		14	#include <glib.h>
		15	#include <gtk/gtk.h>
		16	#include <gdk/gdkx.h>
5	#include <gdk-pixbuf/gdk-pixbuf.h>	17	#include <gdk-pixbuf/gdk-pixbuf.h>
6		18
7	#include <gperl.h>	19	#include <gperl.h>
8	#include <gtk2perl.h>	20	#include <gtk2perl.h>
9		21
10	#define IW 80	22	#include <assert.h>
11		23
12	#define RAND (seed = (seed + 7141) * 54773 % 134456)	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
		36	#define IW 80 /* MUST match Schnauzer.pm! */
		37	#define IH 60 /* MUST match Schnauzer.pm! */
		38
		39	#define LINELENGTH 240
13		40
14	#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 */
		46
		47	struct jpg_err_mgr
		48	{
		49	struct jpeg_error_mgr err;
		50	jmp_buf setjmp_buffer;
		51	};
		52
		53	static void
		54	cv_error_exit (j_common_ptr cinfo)
		55	{
		56	cinfo->err->output_message (cinfo);
		57	longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
		58	}
		59
		60	static void
		61	cv_error_output (j_common_ptr cinfo)
		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);
		68	return;
		69	}
		70
		71	static void
		72	rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
		73	unsigned int h, unsigned int s, unsigned int *v)
		74	{
		75	unsigned int mx = r; if (g > mx) mx = g; if (b > mx) mx = b;
		76	unsigned int mn = r; if (g < mn) mn = g; if (b < mn) mn = b;
		77	unsigned int delta = mx - mn;
		78
		79	*v = mx;
		80
		81	s = mx ? delta 255 / mx : 0;
		82
		83	if (delta == 0)
		84	*h = 0;
		85	else
		86	{
		87	if (r == mx)
		88	h = ((int)g - (int)b) 255 / (int)(delta * 3);
		89	else if (g == mx)
		90	h = ((int)b - (int)r) 255 / (int)(delta * 3) + 52;
		91	else if (b == mx)
		92	h = ((int)r - (int)g) 255 / (int)(delta * 3) + 103;
		93
		94	*h &= 255;
		95	}
		96	}
		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	}
15		152
16	static guint32 a85_val;	153	static guint32 a85_val;
17	static guint a85_cnt;	154	static guint a85_cnt;
18	static guchar a85_buf[80], *a85_ptr;	155	static guchar a85_buf[LINELENGTH], *a85_ptr;
19		156
20	static void	157	static void
21	a85_init (void)	158	a85_init (void)
22	{	159	{
23	a85_cnt = 4;	160	a85_cnt = 4;
32	if (!--a85_cnt)	169	if (!--a85_cnt)
33	{	170	{
34	a85_cnt = 4;	171	a85_cnt = 4;
35	if (a85_val)	172	if (a85_val)
36	{	173	{
37	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;	174	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
38	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;	175	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
39	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;	176	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
40	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;	177	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
41	a85_ptr[0] = (a85_val ) + 33;	178	a85_ptr[0] = (a85_val ) + 33;
42		179
43	a85_ptr += 5;	180	a85_ptr += 5;
44	}	181	}
…		…
53	}	190	}
54	}	191	}
55		192
56	}	193	}
57		194
		195	static void
58	a85_finish (PerlIO *fp)	196	a85_finish (PerlIO *fp)
59	{	197	{
60	while (a85_cnt != 4)	198	while (a85_cnt != 4)
61	a85_push (fp, 0);	199	a85_push (fp, 0);
62		200
…		…
66		204
67	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);	205	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
68	}	206	}
69		207
70	/////////////////////////////////////////////////////////////////////////////	208	/////////////////////////////////////////////////////////////////////////////
		209	// memory source for libjpeg
71		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	/////////////////////////////////////////////////////////////////////////////
		303
72	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::ImageWindow	304	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
73		305
74	PROTOTYPES: ENABLE	306	PROTOTYPES: ENABLE
75		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
		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);
		363
		364	gboolean
		365	gdk_net_wm_supports (GdkAtom property)
		366	CODE:
		367	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
		368	RETVAL = gdk_net_wm_supports (property);
		369	#else
		370	RETVAL = 0;
		371	#endif
		372	OUTPUT:
		373	RETVAL
		374
76	GdkPixbuf_noinc *	375	GdkPixbuf_noinc *
77	transpose (GdkPixbuf *pb)	376	dealpha_expose (GdkPixbuf *pb)
78	CODE:	377	CODE:
		378	perlinterp_release ();
79	{	379	{
80	int w = gdk_pixbuf_get_width (pb);	380	int w = gdk_pixbuf_get_width (pb);
81	int h = gdk_pixbuf_get_height (pb);	381	int h = gdk_pixbuf_get_height (pb);
82	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;	382	int bpp = gdk_pixbuf_get_n_channels (pb);
83	int x, y, i;	383	int x, y, i;
84	guchar src = gdk_pixbuf_get_pixels (pb), dst;	384	guchar src = gdk_pixbuf_get_pixels (pb), dst;
85	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	385	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
86		386
87	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w);	387	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
88		388
89	dst = gdk_pixbuf_get_pixels (RETVAL);	389	dst = gdk_pixbuf_get_pixels (RETVAL);
90	dstr = gdk_pixbuf_get_rowstride (RETVAL);	390	dstr = gdk_pixbuf_get_rowstride (RETVAL);
91		391
92	for (y = 0; y < h; y++)
93	for (x = 0; x < w; x++)	392	for (x = 0; x < w; x++)
		393	for (y = 0; y < h; y++)
94	for (i = 0; i < bpp; i++)	394	for (i = 0; i < 3; i++)
95	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];
96	}	396	}
		397	perlinterp_acquire ();
97	OUTPUT:	398	OUTPUT:
98	RETVAL	399	RETVAL
99		400
100	GdkPixbuf_noinc *	401	GdkPixbuf_noinc *
101	flop (GdkPixbuf *pb)	402	rotate (GdkPixbuf *pb, int angle)
102	CODE:	403	CODE:
103	{	404	perlinterp_release ();
104	int w = gdk_pixbuf_get_width (pb);	405	if (angle < 0)
105	int h = gdk_pixbuf_get_height (pb);	406	angle += 360;
106	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;	407	RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
107	int x, y, i;	408	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
108	guchar src = gdk_pixbuf_get_pixels (pb), dst;	409	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
109	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	410	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
110		411	: angle);
111	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h);	412	perlinterp_acquire ();
112
113	dst = gdk_pixbuf_get_pixels (RETVAL);
114	dstr = gdk_pixbuf_get_rowstride (RETVAL);
115
116	for (y = 0; y < h; y++)
117	for (x = 0; x < w; x++)
118	for (i = 0; i < bpp; i++)
119	dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
120	}
121	OUTPUT:	413	OUTPUT:
122	RETVAL	414	RETVAL
123		415
		416	const char *
		417	filetype (SV *image_data)
		418	CODE:
		419	{
		420	STRLEN data_len;
		421	U8 *data = SvPVbyte (image_data, data_len);
		422	static const unsigned char jxl_header[] = {
		423	0, 0, 0, 0x0c, 0x4a, 0x58, 0x4c, 0x20, 0x0d, 0xa, 0x87, 0x0a
		424	};
		425
		426	if (data_len >= 20
		427	&& data[0] == 0xff
		428	&& data[1] == 0xd8
		429	&& data[2] == 0xff)
		430	RETVAL = "image/jpeg";
		431	else if (data_len >= 12
		432	&& data[ 0] == (U8)'R'
		433	&& data[ 1] == (U8)'I'
		434	&& data[ 2] == (U8)'F'
		435	&& data[ 3] == (U8)'F'
		436	&& data[ 8] == (U8)'W'
		437	&& data[ 9] == (U8)'E'
		438	&& data[10] == (U8)'B'
		439	&& data[11] == (U8)'P')
		440	RETVAL = "image/webp";
		441	else if (data_len >= 16
		442	&& data[ 0] == 0x89
		443	&& data[ 1] == (U8)'P'
		444	&& data[ 2] == (U8)'N'
		445	&& data[ 3] == (U8)'G'
		446	&& data[ 4] == 0x0d
		447	&& data[ 5] == 0x0a
		448	&& data[ 6] == 0x1a
		449	&& data[ 7] == 0x0a)
		450	RETVAL = "image/png";
		451	else if (data_len >= sizeof (jxl_header) && memcmp (data, jxl_header, sizeof (jxl_header)) == 0)
		452	RETVAL = "image/jxl"; // todo: might want to use JxlSignatureCheck
		453	else if (data_len >= 2
		454	&& data[0] == 0xff
		455	&& data[1] == 0x0a)
		456	RETVAL = "image/jxl";
		457	else if (data_len >= 13
		458	&& data[0] == 'G'
		459	&& data[1] == 'I'
		460	&& data[2] == 'F'
		461	&& data[3] == '8'
		462	//&& (data[4] == '7' \|\| data[4] == '9')
		463	&& data[5] == 'a')
		464	{
		465	RETVAL = "image/gif";
		466
		467	// now see if its animated - we require the netscape application header for this
		468	int ofs = 13;
		469
		470	if (data[10] & 0x80)
		471	ofs += (1 << ((data[10] & 7) + 1)) * 3;
		472
		473	if (data_len >= ofs + 2 + 1 + 11)
		474	{
		475
		476	// skip a graphic control extension block. we assume
		477	// there is at most one such block - while the NAB
		478	// has to come firstz, some files do not obey this
		479	if (data[ofs] == 0x21 && data[ofs + 1] == 0xf9)
		480	ofs += 3 + data[ofs + 2] + 1;
		481
		482	if (data_len >= ofs + 2 + 1 + 11)
		483	if (!memcmp (data + ofs, "\x21\xff\x0bNETSCAPE2.0", sizeof ("\x21\xff\x0bNETSCAPE2.0") - 1))
		484	RETVAL = "video/gif";
		485	}
		486	}
		487
		488	else
		489	XSRETURN_UNDEF;
		490	}
		491	OUTPUT:
		492	RETVAL
		493
		494	GdkPixbuf_noinc *
		495	decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		496	CODE:
		497	{
		498	#if WEBP
		499	STRLEN data_size;
		500	int alpha;
		501	WebPData data;
		502	WebPDemuxer *demux;
		503	WebPIterator iter;
		504	WebPDecoderConfig config;
		505	int inw, inh;
		506
		507	data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
		508	data.size = data_size;
		509
		510	perlinterp_release ();
		511
		512	RETVAL = 0;
		513
		514	if (!(demux = WebPDemux (&data)))
		515	goto err_demux;
		516
		517	if (!WebPDemuxGetFrame (demux, 1, &iter))
		518	goto err_iter;
		519
		520	if (!WebPInitDecoderConfig (&config))
		521	goto err_iter;
		522
		523	config.options.use_threads = 1;
		524
		525	if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
		526	goto err_iter;
		527
		528	inw = config.input.width;
		529	inh = config.input.height;
		530
		531	if (thumbnail)
		532	{
		533	if (inw * ih > inh * iw)
		534	ih = (iw * inh + inw - 1) / inw;
		535	else
		536	iw = (ih * inw + inh - 1) / inh;
		537
		538	config.options.bypass_filtering = 1;
		539	config.options.no_fancy_upsampling = 1;
		540
		541	config.options.use_scaling = 1;
		542	config.options.scaled_width = iw;
		543	config.options.scaled_height = ih;
		544	}
		545	else
		546	{
		547	iw = inw;
		548	ih = inh;
		549	}
		550
		551	alpha = !!config.input.has_alpha;
		552
		553	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
		554	if (!RETVAL)
		555	goto err_iter;
		556
		557	config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
		558	config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
		559	config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
		560	config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
		561	config.output.is_external_memory = 1;
		562
		563	if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
		564	{
		565	g_object_unref (RETVAL);
		566	RETVAL = 0;
		567	goto err_iter;
		568	}
		569
		570	err_iter:
		571	WebPDemuxReleaseIterator (&iter);
		572	err_demux:
		573	WebPDemuxDelete (demux);
		574
		575	perlinterp_acquire ();
		576	#else
		577	croak ("load_webp: webp not enabled at compile time");
		578	#endif
		579	}
		580	OUTPUT:
		581	RETVAL
		582
		583	GdkPixbuf_noinc *
		584	decode_jxl (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		585	CODE:
		586	{
		587	#if JXL
		588	JxlDecoder *dec;
		589	JxlBasicInfo info;
		590	const uint8_t next_in = (uint8_t )SvPVbyte_nolen (image_data);
		591	size_t avail_in = SvCUR (image_data);
		592	const char *error = 0;
		593	JxlDecoderStatus status;
		594	void *runner = 0;
		595
		596	RETVAL = 0;
		597
		598	perlinterp_release ();
		599
		600	dec = JxlDecoderCreate (0);
		601
		602	error = "JxlDecoderCreate failed";
		603	if (!dec)
		604	goto done;
		605
		606	runner = JxlThreadParallelRunnerCreate (0, JxlThreadParallelRunnerDefaultNumWorkerThreads ());
		607
		608	status = JxlDecoderSetParallelRunner (dec, JxlThreadParallelRunner, runner);
		609	error = "JxlDecoderSetParallelRunner failed";
		610	if (status != JXL_DEC_SUCCESS)
		611	goto done;
		612
		613	error = "JxlDecoderSubscribeEvents failed";
		614	status = JxlDecoderSubscribeEvents (dec, JXL_DEC_BASIC_INFO \| JXL_DEC_FULL_IMAGE);
		615	if (status != JXL_DEC_SUCCESS)
		616	goto done;
		617
		618	status = JxlDecoderSetInput (dec, next_in, avail_in);
		619	error = "JxlDecoderSetInput failed";
		620	if (status != JXL_DEC_SUCCESS)
		621	goto done;
		622
		623	for (;;)
		624	{
		625	status = JxlDecoderProcessInput (dec);
		626
		627	switch (status)
		628	{
		629	case JXL_DEC_FULL_IMAGE:
		630	error = 0;
		631	goto done;
		632
		633	case JXL_DEC_ERROR:
		634	error = "JxlDecoderProcessInput failed";
		635	goto done;
		636
		637	case JXL_DEC_NEED_MORE_INPUT:
		638	error = "incomplete file";
		639	goto done;
		640
		641	case JXL_DEC_SUCCESS:
		642	error = "incomplete decode";
		643	goto done;
		644
		645	case JXL_DEC_BASIC_INFO:
		646	{
		647	status = JxlDecoderGetBasicInfo (dec, &info);
		648	error = "JxlDecoderGetBasicInfo failed";
		649	if (status != JXL_DEC_SUCCESS)
		650	goto done;
		651
		652	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, !!info.alpha_bits, 8, info.xsize, info.ysize);
		653	error = "unable to allocate pixbuf";
		654	if (!RETVAL)
		655	goto done;
		656
		657	JxlPixelFormat format = {
		658	info.alpha_bits ? 4 : 3,
		659	JXL_TYPE_UINT8,
		660	JXL_NATIVE_ENDIAN,
		661	gdk_pixbuf_get_rowstride (RETVAL)
		662	};
		663
		664	// cannot use gdk_pixbuf_get_byte_length because that does
		665	// not return the size of the buffer, but the size of the buffer without
		666	// the last padding bytes. the internal buffer is rowstride * ysize,
		667	// and this is what the jxl decoder needs. none of this is documented
		668	// in either library, of course.
		669
		670	status = JxlDecoderSetImageOutBuffer (
		671	dec,
		672	&format,
		673	gdk_pixbuf_get_pixels (RETVAL),
		674	gdk_pixbuf_get_rowstride (RETVAL) * info.ysize
		675	);
		676	error = "JxlDecoderSetImageOutBuffer failed";
		677	if (status != JXL_DEC_SUCCESS)
		678	goto done;
		679	}
		680	break;
		681
		682	default:
		683	error = "unexpected event";
		684	goto done;
		685	}
		686	}
		687
		688	done:
		689	if (dec)
		690	JxlDecoderDestroy (dec);
		691
		692	if (runner)
		693	JxlThreadParallelRunnerDestroy (runner);
		694
		695	perlinterp_acquire ();
		696
		697	if (error)
		698	{
		699	if (RETVAL)
		700	g_object_unref (RETVAL);
		701
		702	croak ("load_jxl: %s (status %d)", error, status);
		703	}
		704	#else
		705	croak ("load_jxl: jpeg-xl not enabled at compile time");
		706	#endif
		707	}
		708	OUTPUT:
		709	RETVAL
		710
		711	GdkPixbuf_noinc *
		712	decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		713	CODE:
		714	{
		715	struct jpeg_decompress_struct cinfo;
		716	struct jpg_err_mgr jerr;
		717	int rs;
		718	volatile GdkPixbuf *pb = 0;
		719	STRLEN data_len;
		720	guchar *data = SvPVbyte (image_data, data_len);
		721
		722	RETVAL = 0;
		723
		724	perlinterp_release ();
		725
		726	cinfo.err = jpeg_std_error (&jerr.err);
		727
		728	jerr.err.error_exit = cv_error_exit;
		729	jerr.err.output_message = cv_error_output;
		730
		731	if ((rs = setjmp (jerr.setjmp_buffer)))
		732	{
		733	jpeg_destroy_decompress (&cinfo);
		734
		735	if (pb)
		736	g_object_unref ((gpointer)pb);
		737
		738	perlinterp_acquire ();
		739	XSRETURN_UNDEF;
		740	}
		741
		742	if (!data_len)
		743	longjmp (jerr.setjmp_buffer, 4);
		744
		745	jpeg_create_decompress (&cinfo);
		746	cv_jpeg_mem_src (&cinfo, data, data_len);
		747
		748	jpeg_read_header (&cinfo, TRUE);
		749
		750	cinfo.dct_method = JDCT_DEFAULT;
		751	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
		752	cinfo.do_block_smoothing = FALSE;
		753	cinfo.out_color_space = JCS_RGB;
		754	cinfo.quantize_colors = FALSE;
		755
		756	cinfo.scale_num = 1;
		757	cinfo.scale_denom = 1;
		758
		759	jpeg_calc_output_dimensions (&cinfo);
		760
		761	if (thumbnail)
		762	{
		763	cinfo.dct_method = JDCT_FASTEST;
		764	cinfo.do_fancy_upsampling = FALSE;
		765
		766	while (cinfo.scale_denom < 8
		767	&& cinfo.output_width >= iw*4
		768	&& cinfo.output_height >= ih*4)
		769	{
		770	cinfo.scale_denom <<= 1;
		771	jpeg_calc_output_dimensions (&cinfo);
		772	}
		773	}
		774
		775	if (cinfo.output_components != 3)
		776	longjmp (jerr.setjmp_buffer, 3);
		777
		778	if (cinfo.jpeg_color_space == JCS_YCCK \|\| cinfo.jpeg_color_space == JCS_CMYK)
		779	{
		780	cinfo.out_color_space = JCS_CMYK;
		781	cinfo.output_components = 4;
		782	}
		783
		784	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
		785	if (!RETVAL)
		786	longjmp (jerr.setjmp_buffer, 2);
		787
		788	data = gdk_pixbuf_get_pixels (RETVAL);
		789	rs = gdk_pixbuf_get_rowstride (RETVAL);
		790
		791	jpeg_start_decompress (&cinfo);
		792
		793	while (cinfo.output_scanline < cinfo.output_height)
		794	{
		795	int remaining = cinfo.output_height - cinfo.output_scanline;
		796	JSAMPROW rp[4];
		797
		798	rp [0] = data + cinfo.output_scanline * rs;
		799	rp [1] = (guchar *)rp [0] + rs;
		800	rp [2] = (guchar *)rp [1] + rs;
		801	rp [3] = (guchar *)rp [2] + rs;
		802
		803	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
		804	}
		805
		806	if (cinfo.out_color_space == JCS_CMYK)
		807	{
		808	guchar end = data + cinfo.output_height rs;
		809
		810	while (data < end)
		811	{
		812	U32 c = data [0];
		813	U32 m = data [1];
		814	U32 y = data [2];
		815	U32 k = data [3];
		816
		817	if (0)
		818	if (cinfo.Adobe_transform == 2)
		819	{
		820	c ^= 0xff;
		821	m ^= 0xff;
		822	y ^= 0xff;
		823	k ^= 0xff;
		824	}
		825
		826	data [0] = (c * k + 0x80) / 0xff;
		827	data [1] = (m * k + 0x80) / 0xff;
		828	data [2] = (y * k + 0x80) / 0xff;
		829	data [3] = 0xff;
		830
		831	data += 4;
		832	}
		833	}
		834
		835	jpeg_finish_decompress (&cinfo);
		836	jpeg_destroy_decompress (&cinfo);
		837	perlinterp_acquire ();
		838	}
		839	OUTPUT:
		840	RETVAL
		841
		842	void
		843	compare (GdkPixbuf a, GdkPixbuf b)
		844	PPCODE:
		845	perlinterp_release ();
		846	{
		847	int w = gdk_pixbuf_get_width (a);
		848	int h = gdk_pixbuf_get_height (a);
		849
		850	int sa = gdk_pixbuf_get_rowstride (a);
		851	int sb = gdk_pixbuf_get_rowstride (b);
		852
		853	guchar *pa = gdk_pixbuf_get_pixels (a);
		854	guchar *pb = gdk_pixbuf_get_pixels (b);
		855
		856	int x, y;
		857
		858	assert (w == gdk_pixbuf_get_width (b));
		859	assert (h == gdk_pixbuf_get_height (b));
		860
		861	assert (gdk_pixbuf_get_n_channels (a) == 3);
		862	assert (gdk_pixbuf_get_n_channels (b) == 3);
		863
		864	double diff = 0.;
		865	int peak = 0;
		866
		867	if (w && h)
		868	for (y = 0; y < h; y++)
		869	{
		870	guchar pa_ = pa + y sa;
		871	guchar pb_ = pb + y sb;
		872
		873	for (x = 0; x < w; x++)
		874	{
		875	int d;
		876
		877	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		878	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		879	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		880	}
		881	}
		882
		883	perlinterp_acquire ();
		884
		885	EXTEND (SP, 2);
		886	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		887	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		888	}
		889
124	#############################################################################	890	#############################################################################
125		891
126	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	892	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
127		893
		894	# currently only works for filenames (octet strings)
		895
		896	SV *
		897	foldcase (SV *pathsv)
		898	PROTOTYPE: $
		899	CODE:
		900	{
		901	STRLEN plen;
		902	U8 path = (U8 )SvPV (pathsv, plen);
		903	U8 *pend = path + plen;
		904	U8 dst [plen * 8 * 3], *dstp = dst;
		905
		906	while (path < pend)
		907	{
		908	U8 ch = *path;
		909
		910	if (ch >= 'a' && ch <= 'z')
		911	dstp++ = path++;
		912	else if (ch >= 'A' && ch <= 'Z')
		913	dstp++ = path++ + ('a' - 'A');
		914	else if (ch >= '0' && ch <= '9')
		915	{
		916	/* version sort, up to 8 digits */
		917	STRLEN el, nl = 0;
		918	while (path >= '0' && path <= '9' && path < pend)
		919	path++, nl++;
		920
		921	for (el = nl; el < 8; el++)
		922	*dstp++ = '0';
		923
		924	memcpy (dstp, path - nl, nl);
		925	dstp += nl;
		926	}
		927	else
		928	dstp++ = path++;
		929	#if 0
		930	else
		931	{
		932	STRLEN cl;
		933	to_utf8_fold (path, dstp, &cl);
		934	dstp += cl;
		935	path += is_utf8_char (path);
		936	}
		937	#endif
		938	}
		939
		940	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
		941	}
		942	OUTPUT:
		943	RETVAL
		944
128	GdkPixbuf_noinc *	945	GdkPixbuf_noinc *
129	p7_to_pb (int w, int h, guchar *src)	946	p7_to_pb (int w, int h, SV *src_sv)
		947	PROTOTYPE: @
130	CODE:	948	CODE:
131	{	949	{
132	int x, y;	950	int x, y;
133	guchar dst, d;	951	guchar dst, d;
134	int dstr;	952	int dstr;
		953	guchar src = (guchar )SvPVbyte_nolen (src_sv);
135		954
136	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);	955	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
137	dst = gdk_pixbuf_get_pixels (RETVAL);	956	dst = gdk_pixbuf_get_pixels (RETVAL);
138	dstr = gdk_pixbuf_get_rowstride (RETVAL);	957	dstr = gdk_pixbuf_get_rowstride (RETVAL);
139		958
…		…
148	}	967	}
149	}	968	}
150	OUTPUT:	969	OUTPUT:
151	RETVAL	970	RETVAL
152		971
		972	#############################################################################
		973
		974	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
		975
		976	void
		977	dump_ascii85 (PerlIO fp, GdkPixbuf pb)
		978	CODE:
		979	{
		980	int w = gdk_pixbuf_get_width (pb);
		981	int h = gdk_pixbuf_get_height (pb);
		982	int x, y, i;
		983	guchar *dst;
		984	int bpp = gdk_pixbuf_get_n_channels (pb);
		985	guchar *src = gdk_pixbuf_get_pixels (pb);
		986	int sstr = gdk_pixbuf_get_rowstride (pb);
		987
		988	a85_init ();
		989
		990	for (y = 0; y < h; y++)
		991	for (x = 0; x < w; x++)
		992	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
		993	a85_push (fp, src [x * bpp + y * sstr + i]);
		994
		995	a85_finish (fp);
		996	}
		997
		998	void
		999	dump_binary (PerlIO fp, GdkPixbuf pb)
		1000	CODE:
		1001	{
		1002	int w = gdk_pixbuf_get_width (pb);
		1003	int h = gdk_pixbuf_get_height (pb);
		1004	int x, y, i;
		1005	guchar *dst;
		1006	int bpp = gdk_pixbuf_get_n_channels (pb);
		1007	guchar *src = gdk_pixbuf_get_pixels (pb);
		1008	int sstr = gdk_pixbuf_get_rowstride (pb);
		1009
		1010	for (y = 0; y < h; y++)
		1011	for (x = 0; x < w; x++)
		1012	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
		1013	PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
		1014	}
		1015
		1016	#############################################################################
		1017
		1018	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
		1019
153	SV *	1020	SV *
154	pb_to_p7 (GdkPixbuf *pb)	1021	pb_to_hv84 (GdkPixbuf *pb)
155	CODE:	1022	CODE:
156	{	1023	{
157	int w = gdk_pixbuf_get_width (pb);	1024	int w = gdk_pixbuf_get_width (pb);
158	int h = gdk_pixbuf_get_height (pb);	1025	int h = gdk_pixbuf_get_height (pb);
159	int x, y;	1026	int x, y;
160	guchar *dst;	1027	guchar *dst;
161	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;	1028	int bpp = gdk_pixbuf_get_n_channels (pb);
162	guchar *src = gdk_pixbuf_get_pixels (pb);	1029	guchar *src = gdk_pixbuf_get_pixels (pb);
163	int sstr = gdk_pixbuf_get_rowstride (pb);	1030	int sstr = gdk_pixbuf_get_rowstride (pb);
164	int Er[IW], Eg[IW], Eb[IW];
165	int seed = 77;
166		1031
167	RETVAL = newSV (w * h);	1032	RETVAL = newSV (6 * 8 * 12 / 8);
168	SvPOK_only (RETVAL);	1033	SvPOK_only (RETVAL);
169	SvCUR_set (RETVAL, w * h);	1034	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
170		1035
171	dst = SvPVX (RETVAL);	1036	dst = (guchar *)SvPVX (RETVAL);
172
173	memset (Er, 0, sizeof (int) * IW);
174	memset (Eg, 0, sizeof (int) * IW);
175	memset (Eb, 0, sizeof (int) * IW);
176		1037
177	/* some primitive error distribution + random dithering */	1038	/* some primitive error distribution + random dithering */
178		1039
179	for (y = 0; y < h; y++)	1040	for (y = 0; y < h; y++)
180	{	1041	{
181	int er = 0, eg = 0, eb = 0;	1042	guchar p = src + y sstr;
182		1043
183	for (x = 0; x < w; x++)	1044	for (x = 0; x < w; x += 2)
184	{	1045	{
185	int r, g, b;	1046	unsigned int r, g, b, h, s, v, H, V1, V2;
186	guchar p = src + x bpp + y * sstr;
187		1047
188	r = ((p[0] + er + Er[x]) * 7 + 128) / 255;	1048	if (bpp == 3)
189	g = ((p[1] + eg + Eg[x]) * 7 + 128) / 255;	1049	r = p++, g = p++, b = *p++;
190	b = ((p[2] + eb + Eb[x]) * 3 + 128) / 255;	1050	else if (bpp == 1)
		1051	r = g = b = *p++;
		1052	else
		1053	abort ();
191		1054
192	r = r > 7 ? 7 : r < 0 ? 0 : r;	1055	rgb_to_hsv (r, g, b, &h, &s, &v);
193	g = g > 7 ? 7 : g < 0 ? 0 : g;
194	b = b > 3 ? 3 : b < 0 ? 0 : b;
195		1056
196	er += p[0] - (r * 255 + 4) / 7;	1057	H = (h * 15 / 255) << 4;
197	eg += p[1] - (g * 255 + 4) / 7;	1058	V1 = v;
198	eb += p[2] - (b * 255 + 2) / 3;
199		1059
200	Er[x] = er / 2; er -= er / 2 + RAND % 7 - 3;	1060	if (bpp == 3)
201	Eg[x] = eg / 2; eg -= eg / 2 + RAND % 7 - 3;	1061	r = p++, g = p++, b = *p++;
202	Eb[x] = eb / 2; eb -= eb / 2 + RAND % 7 - 3;	1062	else if (bpp == 1)
		1063	r = g = b = *p++;
		1064	else
		1065	abort ();
203		1066
204	*dst++ = r << 5 \| g << 2 \| b;	1067	rgb_to_hsv (r, g, b, &h, &s, &v);
		1068
		1069	H \|= h * 15 / 255;
		1070	V2 = v;
		1071
		1072	*dst++ = H;
		1073	*dst++ = V1;
		1074	*dst++ = V2;
205	}	1075	}
206	}	1076	}
207	}	1077	}
208	OUTPUT:	1078	OUTPUT:
209	RETVAL	1079	RETVAL
210		1080
		1081	SV *
		1082	hv84_to_av (unsigned char *hv84)
		1083	CODE:
		1084	{
		1085	int i = 72 / 3;
		1086	AV *av = newAV ();
		1087
		1088	RETVAL = (SV )newRV_noinc ((SV )av);
		1089	while (i--)
		1090	{
		1091	int h = *hv84++;
		1092	int v1 = *hv84++;
		1093	int v2 = *hv84++;
		1094
		1095	av_push (av, newSViv (v1));
		1096	av_push (av, newSViv ((h >> 4) * 255 / 15));
		1097	av_push (av, newSViv (v2));
		1098	av_push (av, newSViv ((h & 15) * 255 / 15));
		1099	}
		1100	}
		1101	OUTPUT:
		1102	RETVAL
		1103
211	#############################################################################	1104	#############################################################################
212		1105
213	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript	1106	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
214		1107
215	void	1108	SV *
216	dump_pb (PerlIO fp, GdkPixbuf pb)	1109	extract_features (SV *ar)
217	CODE:	1110	CODE:
218	{	1111	{
219	int w = gdk_pixbuf_get_width (pb);	1112	int i;
220	int h = gdk_pixbuf_get_height (pb);	1113	AV av, result;
221	int x, y, i;
222	guchar *dst;
223	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;
224	guchar *src = gdk_pixbuf_get_pixels (pb);
225	int sstr = gdk_pixbuf_get_rowstride (pb);
226		1114
227	a85_init ();	1115	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
		1116	croak ("Not an array ref as first argument to extract_features");
228		1117
229	for (y = 0; y < h; y++)	1118	av = (AV *) SvRV (ar);
230	for (x = 0; x < w; x++)	1119	result = newAV ();
231	for (i = 0; i < 3; i++)
232	a85_push (fp, src [x * bpp + y * sstr + i]);
233		1120
234	a85_finish (fp);	1121	for (i = 0; i <= av_len (av); ++i)
235	}	1122	{
		1123	SV sv = av_fetch (av, i, 1);
		1124	SV histsv = newSV (9 sizeof (float) + 1);
236		1125
		1126	SvPOK_on (histsv);
		1127	SvCUR_set (histsv, 9 * sizeof (float));
		1128	float hist = (float )SvPVX (histsv);
237		1129
		1130	struct feature f_h, f_s, f_v;
		1131	feature_init (&f_h);
		1132	feature_init (&f_s);
		1133	feature_init (&f_v);
238		1134
		1135	{
		1136	STRLEN len;
		1137	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1138	while (len >= 3)
		1139	{
		1140	unsigned int r, g, b, h, s, v;
		1141	r = buf++; g = buf++; b = *buf++;
		1142	rgb_to_hsv (r, g, b, &h, &s, &v);
239		1143
		1144	feature_update_pass_1 (&f_h, h);
		1145	feature_update_pass_1 (&f_s, s);
		1146	feature_update_pass_1 (&f_v, v);
		1147
		1148	len -= 3;
		1149	}
		1150
		1151	feature_finish_pass_1 (&f_h);
		1152	feature_finish_pass_1 (&f_s);
		1153	feature_finish_pass_1 (&f_v);
		1154	}
		1155
		1156	{
		1157	STRLEN len;
		1158	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1159	while (len >= 3)
		1160	{
		1161	unsigned int r, g, b, h, s, v;
		1162	r = buf++; g = buf++; b = *buf++;
		1163	rgb_to_hsv (r, g, b, &h, &s, &v);
		1164
		1165	feature_update_pass_2 (&f_h, h);
		1166	feature_update_pass_2 (&f_s, s);
		1167	feature_update_pass_2 (&f_v, v);
		1168
		1169	len -= 3;
		1170	}
		1171
		1172	feature_finish_pass_2 (&f_h);
		1173	feature_finish_pass_2 (&f_s);
		1174	feature_finish_pass_2 (&f_v);
		1175	}
		1176
		1177	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		1178	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		1179	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
		1180
		1181	av_push (result, histsv);
		1182	}
		1183
		1184	RETVAL = newRV_noinc ((SV *)result);
		1185	}
		1186	OUTPUT:
		1187	RETVAL
		1188

Diff Legend

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

Comparing CV/CV.xs (file contents): Revision 1.6 by root, Wed Nov 12 00:59:01 2003 UTC vs. Revision 1.64 by root, Wed Dec 1 03:32:40 2021 UTC

Diff Legend

Comparing CV/CV.xs (file contents):
Revision 1.6 by root, Wed Nov 12 00:59:01 2003 UTC vs.
Revision 1.64 by root, Wed Dec 1 03:32:40 2021 UTC