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

Comparing CV/CV.xs (file contents):
Revision 1.31 by root, Wed Feb 15 08:08:37 2006 UTC vs.
Revision 1.60 by root, Thu Jul 15 00:41:44 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>
7		8
8	#include <magic.h>	9	#include <magic.h>
9		10
10	#include <jpeglib.h>	11	#include <jpeglib.h>
		12	#include <jerror.h>
		13
11	#include <glib.h>	14	#include <glib.h>
12	#include <gtk/gtk.h>	15	#include <gtk/gtk.h>
		16	#include <gdk/gdkx.h>
13	#include <gdk-pixbuf/gdk-pixbuf.h>	17	#include <gdk-pixbuf/gdk-pixbuf.h>
14		18
15	#include <gperl.h>	19	#include <gperl.h>
16	#include <gtk2perl.h>	20	#include <gtk2perl.h>
17		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
18	#define IW 80 /* MUST match Schnauer.pm! */	36	#define IW 80 /* MUST match Schnauzer.pm! */
19	#define IH 60 /* MUST match Schnauer.pm! */	37	#define IH 60 /* MUST match Schnauzer.pm! */
20
21	#define RAND (seed = (seed + 7141) * 54773 % 134456)
22		38
23	#define LINELENGTH 240	39	#define LINELENGTH 240
24		40
25	#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 */
26		46
27	struct jpg_err_mgr	47	struct jpg_err_mgr
28	{	48	{
29	struct jpeg_error_mgr err;	49	struct jpeg_error_mgr err;
30	jmp_buf setjmp_buffer;	50	jmp_buf setjmp_buffer;
31	};	51	};
32		52
33	static void	53	static void
34	cv_error_exit (j_common_ptr cinfo)	54	cv_error_exit (j_common_ptr cinfo)
35	{	55	{
		56	cinfo->err->output_message (cinfo);
36	longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);	57	longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
37	}	58	}
38		59
39	static void	60	static void
40	cv_error_output (j_common_ptr cinfo)	61	cv_error_output (j_common_ptr cinfo)
41	{	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);
42	return;	68	return;
43	}	69	}
44		70
45	static void	71	static void
46	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,
…		…
67		93
68	*h &= 255;	94	*h &= 255;
69	}	95	}
70	}	96	}
71		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
72	static guint32 a85_val;	153	static guint32 a85_val;
73	static guint a85_cnt;	154	static guint a85_cnt;
74	static guchar a85_buf[LINELENGTH], *a85_ptr;	155	static guchar a85_buf[LINELENGTH], *a85_ptr;
75		156
76	static void	157	static void
…		…
88	if (!--a85_cnt)	169	if (!--a85_cnt)
89	{	170	{
90	a85_cnt = 4;	171	a85_cnt = 4;
91	if (a85_val)	172	if (a85_val)
92	{	173	{
93	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;	174	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
94	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;	175	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
95	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;	176	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
96	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;	177	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
97	a85_ptr[0] = (a85_val ) + 33;	178	a85_ptr[0] = (a85_val ) + 33;
98		179
99	a85_ptr += 5;	180	a85_ptr += 5;
100	}	181	}
…		…
123		204
124	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);	205	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
125	}	206	}
126		207
127	/////////////////////////////////////////////////////////////////////////////	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	/////////////////////////////////////////////////////////////////////////////
128		303
129	MODULE = Gtk2::CV PACKAGE = Gtk2::CV	304	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
130		305
131	PROTOTYPES: ENABLE	306	PROTOTYPES: ENABLE
132		307
		308	# calculate the common prefix length of two strings
133	# missing function in perl. really :)	309	# missing function in perl. really :)
134	int	310	int
135	common_prefix_length (a, b)	311	common_prefix_length (a, b)
136	unsigned char a = (unsigned char )SvPVutf8_nolen ($arg);	312	unsigned char a = (unsigned char )SvPVutf8_nolen ($arg);
137	unsigned char b = (unsigned char )SvPVutf8_nolen ($arg);	313	unsigned char b = (unsigned char )SvPVutf8_nolen ($arg);
…		…
146		322
147	OUTPUT:	323	OUTPUT:
148	RETVAL	324	RETVAL
149		325
150	const char *	326	const char *
151	magic (const char *path)	327	magic (SV *path_or_data)
		328	ALIAS:
		329	magic = 0
		330	magic_mime = 1
		331	magic_buffer = 2
		332	magic_buffer_mime = 3
152	CODE:	333	CODE:
153	{	334	{
154	static magic_t cookie;	335	STRLEN len;
		336	char *data = SvPVbyte (path_or_data, len);
155		337
156	if (!cookie)	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);
157	{	344	}
158	cookie = magic_open (MAGIC_MIME);
159		345
160	if (cookie)	346	perlinterp_release ();
161	magic_load (cookie, 0);
162	else
163	XSRETURN_UNDEF;
164	}
165		347
166	RETVAL = magic_file (cookie, path);	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 ();
167	}	353	}
168	OUTPUT:	354	OUTPUT:
169	RETVAL	355	RETVAL
170		356
171	# 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);
172		363
173	gboolean	364	gboolean
174	gdk_net_wm_supports (GdkAtom property)	365	gdk_net_wm_supports (GdkAtom property)
175	CODE:	366	CODE:
176	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)	367	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
…		…
182	RETVAL	373	RETVAL
183		374
184	GdkPixbuf_noinc *	375	GdkPixbuf_noinc *
185	dealpha_expose (GdkPixbuf *pb)	376	dealpha_expose (GdkPixbuf *pb)
186	CODE:	377	CODE:
		378	perlinterp_release ();
187	{	379	{
188	int w = gdk_pixbuf_get_width (pb);	380	int w = gdk_pixbuf_get_width (pb);
189	int h = gdk_pixbuf_get_height (pb);	381	int h = gdk_pixbuf_get_height (pb);
190	int bpp = gdk_pixbuf_get_n_channels (pb);	382	int bpp = gdk_pixbuf_get_n_channels (pb);
191	int x, y, i;	383	int x, y, i;
…		…
200	for (x = 0; x < w; x++)	392	for (x = 0; x < w; x++)
201	for (y = 0; y < h; y++)	393	for (y = 0; y < h; y++)
202	for (i = 0; i < 3; i++)	394	for (i = 0; i < 3; i++)
203	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];
204	}	396	}
		397	perlinterp_acquire ();
205	OUTPUT:	398	OUTPUT:
206	RETVAL	399	RETVAL
207		400
208	GdkPixbuf_noinc *	401	GdkPixbuf_noinc *
209	rotate (GdkPixbuf *pb, int angle)	402	rotate (GdkPixbuf *pb, int angle)
210	CODE:	403	CODE:
		404	perlinterp_release ();
		405	if (angle < 0)
		406	angle += 360;
211	RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE	407	RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
212	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE	408	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
213	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN	409	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
214	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE	410	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
215	: angle);	411	: angle);
		412	perlinterp_acquire ();
216	OUTPUT:	413	OUTPUT:
217	RETVAL	414	RETVAL
218		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
219	GdkPixbuf_noinc *	462	GdkPixbuf_noinc *
220	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)
221	CODE:	665	CODE:
222	{	666	{
223	struct jpeg_decompress_struct cinfo;	667	struct jpeg_decompress_struct cinfo;
224	struct jpg_err_mgr jerr;	668	struct jpg_err_mgr jerr;
225	guchar *data;
226	int rs;	669	int rs;
227	FILE *fp;
228	volatile GdkPixbuf *pb = 0;	670	volatile GdkPixbuf *pb = 0;
		671	STRLEN data_len;
		672	guchar *data = SvPVbyte (image_data, data_len);
229		673
230	RETVAL = 0;	674	RETVAL = 0;
231		675
232	fp = fopen (SvPVbyte_nolen (path), "rb");	676	perlinterp_release ();
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 * 8 * 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	{
		868	/* version sort, up to 8 digits */
337	STRLEN el, nl = 0;	869	STRLEN el, nl = 0;
338	while (path >= '0' && path <= '9' && path < pend)	870	while (path >= '0' && path <= '9' && path < pend)
339	path++, nl++;	871	path++, nl++;
340		872
341	for (el = nl; el < 6; el++)	873	for (el = nl; el < 8; el++)
342	*dstp++ = '0';	874	*dstp++ = '0';
343		875
344	memcpy (dstp, path - nl, nl);	876	memcpy (dstp, path - nl, nl);
345	dstp += nl;	877	dstp += nl;
346	}	878	}
		879	else
		880	dstp++ = path++;
		881	#if 0
347	else	882	else
348	{	883	{
349	STRLEN cl;	884	STRLEN cl;
350	to_utf8_fold (path, dstp, &cl);	885	to_utf8_fold (path, dstp, &cl);
351	dstp += cl;	886	dstp += cl;
352	path += is_utf8_char (path);	887	path += is_utf8_char (path);
353	}	888	}
		889	#endif
354	}	890	}
355		891
356	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	892	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
357	}	893	}
358	OUTPUT:	894	OUTPUT:
…		…
520	#############################################################################	1056	#############################################################################
521		1057
522	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	1058	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
523		1059
524	SV *	1060	SV *
525	make_histograms (SV *ar)	1061	extract_features (SV *ar)
526	CODE:	1062	CODE:
527	{	1063	{
528	int i;	1064	int i;
529	AV av, result;	1065	AV av, result;
530		1066
531	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	1067	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
532	croak ("Not an array ref as first argument to make_histogram");	1068	croak ("Not an array ref as first argument to extract_features");
533		1069
534	av = (AV *) SvRV (ar);	1070	av = (AV *) SvRV (ar);
535	result = newAV ();	1071	result = newAV ();
536		1072
537	for (i = 0; i <= av_len (av); ++i)	1073	for (i = 0; i <= av_len (av); ++i)
538	{	1074	{
539	const int HISTSIZE = 64;
540
541	int j;
542	SV sv = av_fetch (av, i, 1);	1075	SV sv = av_fetch (av, i, 1);
543	STRLEN len;
544	char *buf = SvPVbyte (sv, len);
545
546	int tmphist[HISTSIZE];
547	float *hist;
548
549	SV histsv = newSV (HISTSIZE sizeof (float) + 1);	1076	SV histsv = newSV (9 sizeof (float) + 1);
		1077
550	SvPOK_on (histsv);	1078	SvPOK_on (histsv);
551	SvCUR_set (histsv, HISTSIZE * sizeof (float));	1079	SvCUR_set (histsv, 9 * sizeof (float));
552	hist = (float *)SvPVX (histsv);	1080	float hist = (float )SvPVX (histsv);
553		1081
554	Zero (tmphist, sizeof (tmphist), char);	1082	struct feature f_h, f_s, f_v;
		1083	feature_init (&f_h);
		1084	feature_init (&f_s);
		1085	feature_init (&f_v);
555		1086
556	for (j = len; j--; )
557	{
558	unsigned int idx
559	= ((*buf & 0xc0) >> 2)
560	\| ((*buf & 0x18) >> 1)
561	\| (*buf & 0x03);
562
563	++tmphist[idx];
564	++buf;
565	}
566		1087	{
567	for (j = 0; j < HISTSIZE; ++j)	1088	STRLEN len;
568	hist[j] = (float)tmphist[j] / (len + 1e-30);	1089	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1090	while (len >= 3)
		1091	{
		1092	unsigned int r, g, b, h, s, v;
		1093	r = buf++; g = buf++; b = *buf++;
		1094	rgb_to_hsv (r, g, b, &h, &s, &v);
		1095
		1096	feature_update_pass_1 (&f_h, h);
		1097	feature_update_pass_1 (&f_s, s);
		1098	feature_update_pass_1 (&f_v, v);
		1099
		1100	len -= 3;
		1101	}
		1102
		1103	feature_finish_pass_1 (&f_h);
		1104	feature_finish_pass_1 (&f_s);
		1105	feature_finish_pass_1 (&f_v);
		1106	}
		1107
		1108	{
		1109	STRLEN len;
		1110	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		1111	while (len >= 3)
		1112	{
		1113	unsigned int r, g, b, h, s, v;
		1114	r = buf++; g = buf++; b = *buf++;
		1115	rgb_to_hsv (r, g, b, &h, &s, &v);
		1116
		1117	feature_update_pass_2 (&f_h, h);
		1118	feature_update_pass_2 (&f_s, s);
		1119	feature_update_pass_2 (&f_v, v);
		1120
		1121	len -= 3;
		1122	}
		1123
		1124	feature_finish_pass_2 (&f_h);
		1125	feature_finish_pass_2 (&f_s);
		1126	feature_finish_pass_2 (&f_v);
		1127	}
		1128
		1129	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		1130	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		1131	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
569		1132
570	av_push (result, histsv);	1133	av_push (result, histsv);
571	}	1134	}
572		1135
573	RETVAL = newRV_noinc ((SV *)result);	1136	RETVAL = newRV_noinc ((SV *)result);
574	}	1137	}
575	OUTPUT:	1138	OUTPUT:
576	RETVAL	1139	RETVAL
577		1140
578

Diff Legend

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

Comparing CV/CV.xs (file contents): Revision 1.31 by root, Wed Feb 15 08:08:37 2006 UTC vs. Revision 1.60 by root, Thu Jul 15 00:41:44 2021 UTC

Diff Legend

Comparing CV/CV.xs (file contents):
Revision 1.31 by root, Wed Feb 15 08:08:37 2006 UTC vs.
Revision 1.60 by root, Thu Jul 15 00:41:44 2021 UTC