[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.57 by root, Mon Jul 30 22:55:46 2018 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>
		23
		24	#if WEBP
		25	#include <webp/demux.h>
		26	#include <webp/decode.h>
		27	#endif
		28
		29	#include "perlmulticore.h"
		30
		31	#define IW 80 /* MUST match Schnauzer.pm! */
		32	#define IH 60 /* MUST match Schnauzer.pm! */
11		33
12	#define RAND (seed = (seed + 7141) * 54773 % 134456)	34	#define RAND (seed = (seed + 7141) * 54773 % 134456)
13		35
		36	#define LINELENGTH 240
		37
14	#define ELLIPSIS "\xe2\x80\xa6"	38	#define ELLIPSIS "\xe2\x80\xa6"
		39
		40	typedef char *octet_string;
		41
		42	static magic_t magic_cookie[2]; /* !mime, mime */
		43
		44	struct jpg_err_mgr
		45	{
		46	struct jpeg_error_mgr err;
		47	jmp_buf setjmp_buffer;
		48	};
		49
		50	static void
		51	cv_error_exit (j_common_ptr cinfo)
		52	{
		53	cinfo->err->output_message (cinfo);
		54	longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
		55	}
		56
		57	static void
		58	cv_error_output (j_common_ptr cinfo)
		59	{
		60	char msg[JMSG_LENGTH_MAX];
		61
		62	cinfo->err->format_message (cinfo, msg);
		63
		64	fprintf (stderr, "JPEG decoding error: %s\n", msg);
		65	return;
		66	}
		67
		68	static void
		69	rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
		70	unsigned int h, unsigned int s, unsigned int *v)
		71	{
		72	unsigned int mx = r; if (g > mx) mx = g; if (b > mx) mx = b;
		73	unsigned int mn = r; if (g < mn) mn = g; if (b < mn) mn = b;
		74	unsigned int delta = mx - mn;
		75
		76	*v = mx;
		77
		78	s = mx ? delta 255 / mx : 0;
		79
		80	if (delta == 0)
		81	*h = 0;
		82	else
		83	{
		84	if (r == mx)
		85	h = ((int)g - (int)b) 255 / (int)(delta * 3);
		86	else if (g == mx)
		87	h = ((int)b - (int)r) 255 / (int)(delta * 3) + 52;
		88	else if (b == mx)
		89	h = ((int)r - (int)g) 255 / (int)(delta * 3) + 103;
		90
		91	*h &= 255;
		92	}
		93	}
		94
		95	struct feature {
		96	float v1, v2, v3; // mean, square, cube
		97	int n;
		98	};
		99
		100	static void
		101	feature_init (struct feature *f)
		102	{
		103	f->v1 = 0.;
		104	f->v2 = 0.;
		105	f->v3 = 0.;
		106	f->n = 0;
		107	}
		108
		109	// didn't find an algorithm to neatly do mean, variance and skew in one pass.
		110	// elmex ist schuld.
		111	static void
		112	feature_update_pass_1 (struct feature *f, unsigned int v)
		113	{
		114	f->v1 += v;
		115	f->n += 1;
		116	}
		117
		118	static void
		119	feature_finish_pass_1 (struct feature *f)
		120	{
		121	if (f->n < 1)
		122	return;
		123
		124	f->v1 /= f->n;
		125	}
		126
		127	static void
		128	feature_update_pass_2 (struct feature *f, unsigned int v)
		129	{
		130	float d = v - f->v1;
		131
		132	f->v2 += d * d;
		133	f->v3 += d * d * d;
		134	}
		135
		136	static void
		137	feature_finish_pass_2 (struct feature *f)
		138	{
		139	if (f->n < 1)
		140	return;
		141
		142	f->v2 /= f->n;
		143	f->v3 /= f->n;
		144
		145	f->v1 /= 255.;
		146	f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2);
		147	f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.);
		148	}
15		149
16	static guint32 a85_val;	150	static guint32 a85_val;
17	static guint a85_cnt;	151	static guint a85_cnt;
18	static guchar a85_buf[80], *a85_ptr;	152	static guchar a85_buf[LINELENGTH], *a85_ptr;
19		153
20	static void	154	static void
21	a85_init (void)	155	a85_init (void)
22	{	156	{
23	a85_cnt = 4;	157	a85_cnt = 4;
32	if (!--a85_cnt)	166	if (!--a85_cnt)
33	{	167	{
34	a85_cnt = 4;	168	a85_cnt = 4;
35	if (a85_val)	169	if (a85_val)
36	{	170	{
37	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;	171	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
38	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;	172	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
39	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;	173	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
40	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;	174	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
41	a85_ptr[0] = (a85_val ) + 33;	175	a85_ptr[0] = (a85_val ) + 33;
42		176
43	a85_ptr += 5;	177	a85_ptr += 5;
44	}	178	}
…		…
53	}	187	}
54	}	188	}
55		189
56	}	190	}
57		191
		192	static void
58	a85_finish (PerlIO *fp)	193	a85_finish (PerlIO *fp)
59	{	194	{
60	while (a85_cnt != 4)	195	while (a85_cnt != 4)
61	a85_push (fp, 0);	196	a85_push (fp, 0);
62		197
…		…
66		201
67	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);	202	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
68	}	203	}
69		204
70	/////////////////////////////////////////////////////////////////////////////	205	/////////////////////////////////////////////////////////////////////////////
		206	// memory source for libjpeg
71		207
		208	static void cv_ms_init (j_decompress_ptr cinfo)
		209	{
		210	}
		211
		212	static void cv_ms_term (j_decompress_ptr cinfo)
		213	{
		214	}
		215
		216	static boolean cv_ms_fill (j_decompress_ptr cinfo)
		217	{
		218	// unexpected EOF, warn and generate fake EOI marker
		219
		220	WARNMS (cinfo, JWRN_JPEG_EOF);
		221
		222	struct jpeg_source_mgr src = (struct jpeg_source_mgr )cinfo->src;
		223
		224	static const JOCTET eoi[] = { 0xFF, JPEG_EOI };
		225
		226	src->next_input_byte = eoi;
		227	src->bytes_in_buffer = sizeof (eoi);
		228
		229	return TRUE;
		230	}
		231
		232	static void cv_ms_skip (j_decompress_ptr cinfo, long num_bytes)
		233	{
		234	struct jpeg_source_mgr src = (struct jpeg_source_mgr )cinfo->src;
		235
		236	src->next_input_byte += num_bytes;
		237	src->bytes_in_buffer -= num_bytes;
		238	}
		239
		240	static void cv_jpeg_mem_src (j_decompress_ptr cinfo, void *buf, size_t buflen)
		241	{
		242	struct jpeg_source_mgr *src;
		243
		244	if (!cinfo->src)
		245	cinfo->src = (struct jpeg_source_mgr *)
		246	(*cinfo->mem->alloc_small) (
		247	(j_common_ptr) cinfo, JPOOL_PERMANENT, sizeof (struct jpeg_source_mgr)
		248	);
		249
		250	src = (struct jpeg_source_mgr *)cinfo->src;
		251	src->init_source = cv_ms_init;
		252	src->fill_input_buffer = cv_ms_fill;
		253	src->skip_input_data = cv_ms_skip;
		254	src->resync_to_restart = jpeg_resync_to_restart;
		255	src->term_source = cv_ms_term;
		256	src->next_input_byte = (JOCTET *)buf;
		257	src->bytes_in_buffer = buflen;
		258	}
		259
		260	/////////////////////////////////////////////////////////////////////////////
		261
72	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::ImageWindow	262	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
73		263
74	PROTOTYPES: ENABLE	264	PROTOTYPES: ENABLE
75		265
		266	# calculate the common prefix length of two strings
		267	# missing function in perl. really :)
		268	int
		269	common_prefix_length (a, b)
		270	unsigned char a = (unsigned char )SvPVutf8_nolen ($arg);
		271	unsigned char b = (unsigned char )SvPVutf8_nolen ($arg);
		272	CODE:
		273	RETVAL = 0;
		274
		275	while (a == b && *a)
		276	{
		277	RETVAL += (*a & 0xc0) != 0x80;
		278	a++, b++;
		279	}
		280
		281	OUTPUT:
		282	RETVAL
		283
		284	const char *
		285	magic (SV *path_or_data)
		286	ALIAS:
		287	magic = 0
		288	magic_mime = 1
		289	magic_buffer = 2
		290	magic_buffer_mime = 3
		291	CODE:
		292	{
		293	STRLEN len;
		294	char *data = SvPVbyte (path_or_data, len);
		295
		296	if (!magic_cookie[0])
		297	{
		298	magic_cookie[0] = magic_open (MAGIC_SYMLINK);
		299	magic_cookie[1] = magic_open (MAGIC_SYMLINK \| MAGIC_MIME_TYPE);
		300	magic_load (magic_cookie[0], 0);
		301	magic_load (magic_cookie[1], 0);
		302	}
		303
		304	perlinterp_release ();
		305
		306	RETVAL = ix & 2
		307	? magic_buffer (magic_cookie[ix & 1], data, len)
		308	: magic_file (magic_cookie[ix & 1], data);
		309
		310	perlinterp_acquire ();
		311	}
		312	OUTPUT:
		313	RETVAL
		314
		315	# missing/broken in Gtk2 perl module
		316
		317	void
		318	gdk_window_clear_hints (GdkWindow *window)
		319	CODE:
		320	gdk_window_set_geometry_hints (window, 0, 0);
		321
		322	gboolean
		323	gdk_net_wm_supports (GdkAtom property)
		324	CODE:
		325	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
		326	RETVAL = gdk_net_wm_supports (property);
		327	#else
		328	RETVAL = 0;
		329	#endif
		330	OUTPUT:
		331	RETVAL
		332
76	GdkPixbuf_noinc *	333	GdkPixbuf_noinc *
77	transpose (GdkPixbuf *pb)	334	dealpha_expose (GdkPixbuf *pb)
78	CODE:	335	CODE:
		336	perlinterp_release ();
79	{	337	{
80	int w = gdk_pixbuf_get_width (pb);	338	int w = gdk_pixbuf_get_width (pb);
81	int h = gdk_pixbuf_get_height (pb);	339	int h = gdk_pixbuf_get_height (pb);
82	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;	340	int bpp = gdk_pixbuf_get_n_channels (pb);
83	int x, y, i;	341	int x, y, i;
84	guchar src = gdk_pixbuf_get_pixels (pb), dst;	342	guchar src = gdk_pixbuf_get_pixels (pb), dst;
85	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	343	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
86		344
87	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w);	345	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
88		346
89	dst = gdk_pixbuf_get_pixels (RETVAL);	347	dst = gdk_pixbuf_get_pixels (RETVAL);
90	dstr = gdk_pixbuf_get_rowstride (RETVAL);	348	dstr = gdk_pixbuf_get_rowstride (RETVAL);
91		349
92	for (y = 0; y < h; y++)
93	for (x = 0; x < w; x++)	350	for (x = 0; x < w; x++)
		351	for (y = 0; y < h; y++)
94	for (i = 0; i < bpp; i++)	352	for (i = 0; i < 3; i++)
95	dst[y * bpp + x * dstr + i] = src[x * bpp + y * sstr + i];	353	dst[x * 3 + y * dstr + i] = src[x * bpp + y * sstr + i];
96	}	354	}
		355	perlinterp_acquire ();
97	OUTPUT:	356	OUTPUT:
98	RETVAL	357	RETVAL
99		358
100	GdkPixbuf_noinc *	359	GdkPixbuf_noinc *
101	flop (GdkPixbuf *pb)	360	rotate (GdkPixbuf *pb, int angle)
102	CODE:	361	CODE:
103	{	362	perlinterp_release ();
104	int w = gdk_pixbuf_get_width (pb);	363	if (angle < 0)
105	int h = gdk_pixbuf_get_height (pb);	364	angle += 360;
106	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;	365	RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
107	int x, y, i;	366	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
108	guchar src = gdk_pixbuf_get_pixels (pb), dst;	367	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
109	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	368	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
110		369	: angle);
111	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h);	370	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:	371	OUTPUT:
122	RETVAL	372	RETVAL
123		373
		374	const char *
		375	filetype (SV *image_data)
		376	CODE:
		377	{
		378	STRLEN data_len;
		379	U8 *data = SvPVbyte (image_data, data_len);
		380
		381	if (data_len >= 20
		382	&& data[0] == 0xff
		383	&& data[1] == 0xd8
		384	&& data[2] == 0xff)
		385	RETVAL = "jpg";
		386	else if (data_len >= 12
		387	&& data[ 0] == (U8)'R'
		388	&& data[ 1] == (U8)'I'
		389	&& data[ 2] == (U8)'F'
		390	&& data[ 3] == (U8)'F'
		391	&& data[ 8] == (U8)'W'
		392	&& data[ 9] == (U8)'E'
		393	&& data[10] == (U8)'B'
		394	&& data[11] == (U8)'P')
		395	RETVAL = "webp";
		396	else if (data_len >= 16
		397	&& data[ 0] == 0x89
		398	&& data[ 1] == (U8)'P'
		399	&& data[ 2] == (U8)'N'
		400	&& data[ 3] == (U8)'G'
		401	&& data[ 4] == 0x0d
		402	&& data[ 5] == 0x0a
		403	&& data[ 6] == 0x1a
		404	&& data[ 7] == 0x0a)
		405	RETVAL = "png";
		406	else
		407	RETVAL = "other";
		408	}
		409	OUTPUT:
		410	RETVAL
		411
		412	GdkPixbuf_noinc *
		413	decode_webp (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		414	CODE:
		415	{
		416	#if WEBP
		417	STRLEN data_size;
		418	int alpha;
		419	WebPData data;
		420	WebPDemuxer *demux;
		421	WebPIterator iter;
		422	WebPDecoderConfig config;
		423	int inw, inh;
		424
		425	data.bytes = (uint8_t *)SvPVbyte (image_data, data_size);
		426	data.size = data_size;
		427
		428	perlinterp_release ();
		429
		430	RETVAL = 0;
		431
		432	if (!(demux = WebPDemux (&data)))
		433	goto err_demux;
		434
		435	if (!WebPDemuxGetFrame (demux, 1, &iter))
		436	goto err_iter;
		437
		438	if (!WebPInitDecoderConfig (&config))
		439	goto err_iter;
		440
		441	config.options.use_threads = 1;
		442
		443	if (WebPGetFeatures (iter.fragment.bytes, iter.fragment.size, &config.input) != VP8_STATUS_OK)
		444	goto err_iter;
		445
		446	inw = config.input.width;
		447	inh = config.input.height;
		448
		449	if (thumbnail)
		450	{
		451	if (inw * ih > inh * iw)
		452	ih = (iw * inh + inw - 1) / inw;
		453	else
		454	iw = (ih * inw + inh - 1) / inh;
		455
		456	config.options.bypass_filtering = 1;
		457	config.options.no_fancy_upsampling = 1;
		458
		459	config.options.use_scaling = 1;
		460	config.options.scaled_width = iw;
		461	config.options.scaled_height = ih;
		462	}
		463	else
		464	{
		465	iw = inw;
		466	ih = inh;
		467	}
		468
		469	alpha = !!config.input.has_alpha;
		470
		471	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, alpha, 8, iw, ih);
		472	if (!RETVAL)
		473	goto err_iter;
		474
		475	config.output.colorspace = alpha ? MODE_RGBA : MODE_RGB;
		476	config.output.u.RGBA.rgba = gdk_pixbuf_get_pixels (RETVAL);
		477	config.output.u.RGBA.stride = gdk_pixbuf_get_rowstride (RETVAL);
		478	config.output.u.RGBA.size = gdk_pixbuf_get_byte_length (RETVAL);
		479	config.output.is_external_memory = 1;
		480
		481	if (WebPDecode (iter.fragment.bytes, iter.fragment.size, &config) != VP8_STATUS_OK)
		482	{
		483	g_object_unref (RETVAL);
		484	RETVAL = 0;
		485	goto err_iter;
		486	}
		487
		488	err_iter:
		489	WebPDemuxReleaseIterator (&iter);
		490	err_demux:
		491	WebPDemuxDelete (demux);
		492
		493	perlinterp_acquire ();
		494	#else
		495	croak ("load_webp: webp not enabled at compile time");
		496	#endif
		497	}
		498	OUTPUT:
		499	RETVAL
		500
		501	GdkPixbuf_noinc *
		502	decode_jpeg (SV *image_data, int thumbnail = 0, int iw = 0, int ih = 0)
		503	CODE:
		504	{
		505	struct jpeg_decompress_struct cinfo;
		506	struct jpg_err_mgr jerr;
		507	int rs;
		508	volatile GdkPixbuf *pb = 0;
		509	STRLEN data_len;
		510	guchar *data = SvPVbyte (image_data, data_len);
		511
		512	RETVAL = 0;
		513
		514	perlinterp_release ();
		515
		516	cinfo.err = jpeg_std_error (&jerr.err);
		517
		518	jerr.err.error_exit = cv_error_exit;
		519	jerr.err.output_message = cv_error_output;
		520
		521	if ((rs = setjmp (jerr.setjmp_buffer)))
		522	{
		523	jpeg_destroy_decompress (&cinfo);
		524
		525	if (pb)
		526	g_object_unref ((gpointer)pb);
		527
		528	perlinterp_acquire ();
		529	XSRETURN_UNDEF;
		530	}
		531
		532	if (!data_len)
		533	longjmp (jerr.setjmp_buffer, 4);
		534
		535	jpeg_create_decompress (&cinfo);
		536	cv_jpeg_mem_src (&cinfo, data, data_len);
		537
		538	jpeg_read_header (&cinfo, TRUE);
		539
		540	cinfo.dct_method = JDCT_DEFAULT;
		541	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
		542	cinfo.do_block_smoothing = FALSE;
		543	cinfo.out_color_space = JCS_RGB;
		544	cinfo.quantize_colors = FALSE;
		545
		546	cinfo.scale_num = 1;
		547	cinfo.scale_denom = 1;
		548
		549	jpeg_calc_output_dimensions (&cinfo);
		550
		551	if (thumbnail)
		552	{
		553	cinfo.dct_method = JDCT_FASTEST;
		554	cinfo.do_fancy_upsampling = FALSE;
		555
		556	while (cinfo.scale_denom < 8
		557	&& cinfo.output_width >= iw*4
		558	&& cinfo.output_height >= ih*4)
		559	{
		560	cinfo.scale_denom <<= 1;
		561	jpeg_calc_output_dimensions (&cinfo);
		562	}
		563	}
		564
		565	if (cinfo.output_components != 3)
		566	longjmp (jerr.setjmp_buffer, 3);
		567
		568	if (cinfo.jpeg_color_space == JCS_YCCK \|\| cinfo.jpeg_color_space == JCS_CMYK)
		569	{
		570	cinfo.out_color_space = JCS_CMYK;
		571	cinfo.output_components = 4;
		572	}
		573
		574	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
		575	if (!RETVAL)
		576	longjmp (jerr.setjmp_buffer, 2);
		577
		578	data = gdk_pixbuf_get_pixels (RETVAL);
		579	rs = gdk_pixbuf_get_rowstride (RETVAL);
		580
		581	jpeg_start_decompress (&cinfo);
		582
		583	while (cinfo.output_scanline < cinfo.output_height)
		584	{
		585	int remaining = cinfo.output_height - cinfo.output_scanline;
		586	JSAMPROW rp[4];
		587
		588	rp [0] = data + cinfo.output_scanline * rs;
		589	rp [1] = (guchar *)rp [0] + rs;
		590	rp [2] = (guchar *)rp [1] + rs;
		591	rp [3] = (guchar *)rp [2] + rs;
		592
		593	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
		594	}
		595
		596	if (cinfo.out_color_space == JCS_CMYK)
		597	{
		598	guchar end = data + cinfo.output_height rs;
		599
		600	while (data < end)
		601	{
		602	U32 c = data [0];
		603	U32 m = data [1];
		604	U32 y = data [2];
		605	U32 k = data [3];
		606
		607	if (0)
		608	if (cinfo.Adobe_transform == 2)
		609	{
		610	c ^= 0xff;
		611	m ^= 0xff;
		612	y ^= 0xff;
		613	k ^= 0xff;
		614	}
		615
		616	data [0] = (c * k + 0x80) / 0xff;
		617	data [1] = (m * k + 0x80) / 0xff;
		618	data [2] = (y * k + 0x80) / 0xff;
		619	data [3] = 0xff;
		620
		621	data += 4;
		622	}
		623	}
		624
		625	jpeg_finish_decompress (&cinfo);
		626	jpeg_destroy_decompress (&cinfo);
		627	perlinterp_acquire ();
		628	}
		629	OUTPUT:
		630	RETVAL
		631
		632	void
		633	compare (GdkPixbuf a, GdkPixbuf b)
		634	PPCODE:
		635	perlinterp_release ();
		636	{
		637	int w = gdk_pixbuf_get_width (a);
		638	int h = gdk_pixbuf_get_height (a);
		639
		640	int sa = gdk_pixbuf_get_rowstride (a);
		641	int sb = gdk_pixbuf_get_rowstride (b);
		642
		643	guchar *pa = gdk_pixbuf_get_pixels (a);
		644	guchar *pb = gdk_pixbuf_get_pixels (b);
		645
		646	int x, y;
		647
		648	assert (w == gdk_pixbuf_get_width (b));
		649	assert (h == gdk_pixbuf_get_height (b));
		650
		651	assert (gdk_pixbuf_get_n_channels (a) == 3);
		652	assert (gdk_pixbuf_get_n_channels (b) == 3);
		653
		654	double diff = 0.;
		655	int peak = 0;
		656
		657	if (w && h)
		658	for (y = 0; y < h; y++)
		659	{
		660	guchar pa_ = pa + y sa;
		661	guchar pb_ = pb + y sb;
		662
		663	for (x = 0; x < w; x++)
		664	{
		665	int d;
		666
		667	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		668	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		669	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		670	}
		671	}
		672
		673	perlinterp_acquire ();
		674
		675	EXTEND (SP, 2);
		676	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		677	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		678	}
		679
124	#############################################################################	680	#############################################################################
125		681
126	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	682	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
127		683
		684	# currently only works for filenames (octet strings)
		685
		686	SV *
		687	foldcase (SV *pathsv)
		688	PROTOTYPE: $
		689	CODE:
		690	{
		691	STRLEN plen;
		692	U8 path = (U8 )SvPV (pathsv, plen);
		693	U8 *pend = path + plen;
		694	U8 dst [plen * 6 * 3], *dstp = dst;
		695
		696	while (path < pend)
		697	{
		698	U8 ch = *path;
		699
		700	if (ch >= 'a' && ch <= 'z')
		701	dstp++ = path++;
		702	else if (ch >= 'A' && ch <= 'Z')
		703	dstp++ = path++ + ('a' - 'A');
		704	else if (ch >= '0' && ch <= '9')
		705	{
		706	STRLEN el, nl = 0;
		707	while (path >= '0' && path <= '9' && path < pend)
		708	path++, nl++;
		709
		710	for (el = nl; el < 6; el++)
		711	*dstp++ = '0';
		712
		713	memcpy (dstp, path - nl, nl);
		714	dstp += nl;
		715	}
		716	else
		717	dstp++ = path++;
		718	#if 0
		719	else
		720	{
		721	STRLEN cl;
		722	to_utf8_fold (path, dstp, &cl);
		723	dstp += cl;
		724	path += is_utf8_char (path);
		725	}
		726	#endif
		727	}
		728
		729	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
		730	}
		731	OUTPUT:
		732	RETVAL
		733
128	GdkPixbuf_noinc *	734	GdkPixbuf_noinc *
129	p7_to_pb (int w, int h, guchar *src)	735	p7_to_pb (int w, int h, SV *src_sv)
		736	PROTOTYPE: @
130	CODE:	737	CODE:
131	{	738	{
132	int x, y;	739	int x, y;
133	guchar dst, d;	740	guchar dst, d;
134	int dstr;	741	int dstr;
		742	guchar src = (guchar )SvPVbyte_nolen (src_sv);
135		743
136	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);	744	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
137	dst = gdk_pixbuf_get_pixels (RETVAL);	745	dst = gdk_pixbuf_get_pixels (RETVAL);
138	dstr = gdk_pixbuf_get_rowstride (RETVAL);	746	dstr = gdk_pixbuf_get_rowstride (RETVAL);
139		747
…		…
148	}	756	}
149	}	757	}
150	OUTPUT:	758	OUTPUT:
151	RETVAL	759	RETVAL
152		760
		761	#############################################################################
		762
		763	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
		764
		765	void
		766	dump_ascii85 (PerlIO fp, GdkPixbuf pb)
		767	CODE:
		768	{
		769	int w = gdk_pixbuf_get_width (pb);
		770	int h = gdk_pixbuf_get_height (pb);
		771	int x, y, i;
		772	guchar *dst;
		773	int bpp = gdk_pixbuf_get_n_channels (pb);
		774	guchar *src = gdk_pixbuf_get_pixels (pb);
		775	int sstr = gdk_pixbuf_get_rowstride (pb);
		776
		777	a85_init ();
		778
		779	for (y = 0; y < h; y++)
		780	for (x = 0; x < w; x++)
		781	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
		782	a85_push (fp, src [x * bpp + y * sstr + i]);
		783
		784	a85_finish (fp);
		785	}
		786
		787	void
		788	dump_binary (PerlIO fp, GdkPixbuf pb)
		789	CODE:
		790	{
		791	int w = gdk_pixbuf_get_width (pb);
		792	int h = gdk_pixbuf_get_height (pb);
		793	int x, y, i;
		794	guchar *dst;
		795	int bpp = gdk_pixbuf_get_n_channels (pb);
		796	guchar *src = gdk_pixbuf_get_pixels (pb);
		797	int sstr = gdk_pixbuf_get_rowstride (pb);
		798
		799	for (y = 0; y < h; y++)
		800	for (x = 0; x < w; x++)
		801	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
		802	PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
		803	}
		804
		805	#############################################################################
		806
		807	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
		808
153	SV *	809	SV *
154	pb_to_p7 (GdkPixbuf *pb)	810	pb_to_hv84 (GdkPixbuf *pb)
155	CODE:	811	CODE:
156	{	812	{
157	int w = gdk_pixbuf_get_width (pb);	813	int w = gdk_pixbuf_get_width (pb);
158	int h = gdk_pixbuf_get_height (pb);	814	int h = gdk_pixbuf_get_height (pb);
159	int x, y;	815	int x, y;
160	guchar *dst;	816	guchar *dst;
161	int bpp = gdk_pixbuf_get_has_alpha (pb) ? 4 : 3;	817	int bpp = gdk_pixbuf_get_n_channels (pb);
162	guchar *src = gdk_pixbuf_get_pixels (pb);	818	guchar *src = gdk_pixbuf_get_pixels (pb);
163	int sstr = gdk_pixbuf_get_rowstride (pb);	819	int sstr = gdk_pixbuf_get_rowstride (pb);
164	int Er[IW], Eg[IW], Eb[IW];
165	int seed = 77;
166		820
167	RETVAL = newSV (w * h);	821	RETVAL = newSV (6 * 8 * 12 / 8);
168	SvPOK_only (RETVAL);	822	SvPOK_only (RETVAL);
169	SvCUR_set (RETVAL, w * h);	823	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
170		824
171	dst = SvPVX (RETVAL);	825	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		826
177	/* some primitive error distribution + random dithering */	827	/* some primitive error distribution + random dithering */
178		828
179	for (y = 0; y < h; y++)	829	for (y = 0; y < h; y++)
180	{	830	{
181	int er = 0, eg = 0, eb = 0;	831	guchar p = src + y sstr;
182		832
183	for (x = 0; x < w; x++)	833	for (x = 0; x < w; x += 2)
184	{	834	{
185	int r, g, b;	835	unsigned int r, g, b, h, s, v, H, V1, V2;
186	guchar p = src + x bpp + y * sstr;
187		836
188	r = ((p[0] + er + Er[x]) * 7 + 128) / 255;	837	if (bpp == 3)
189	g = ((p[1] + eg + Eg[x]) * 7 + 128) / 255;	838	r = p++, g = p++, b = *p++;
190	b = ((p[2] + eb + Eb[x]) * 3 + 128) / 255;	839	else if (bpp == 1)
		840	r = g = b = *p++;
		841	else
		842	abort ();
191		843
192	r = r > 7 ? 7 : r < 0 ? 0 : r;	844	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		845
196	er += p[0] - (r * 255 + 4) / 7;	846	H = (h * 15 / 255) << 4;
197	eg += p[1] - (g * 255 + 4) / 7;	847	V1 = v;
198	eb += p[2] - (b * 255 + 2) / 3;
199		848
200	Er[x] = er / 2; er -= er / 2 + RAND % 7 - 3;	849	if (bpp == 3)
201	Eg[x] = eg / 2; eg -= eg / 2 + RAND % 7 - 3;	850	r = p++, g = p++, b = *p++;
202	Eb[x] = eb / 2; eb -= eb / 2 + RAND % 7 - 3;	851	else if (bpp == 1)
		852	r = g = b = *p++;
		853	else
		854	abort ();
203		855
204	*dst++ = r << 5 \| g << 2 \| b;	856	rgb_to_hsv (r, g, b, &h, &s, &v);
		857
		858	H \|= h * 15 / 255;
		859	V2 = v;
		860
		861	*dst++ = H;
		862	*dst++ = V1;
		863	*dst++ = V2;
205	}	864	}
206	}	865	}
207	}	866	}
208	OUTPUT:	867	OUTPUT:
209	RETVAL	868	RETVAL
210		869
		870	SV *
		871	hv84_to_av (unsigned char *hv84)
		872	CODE:
		873	{
		874	int i = 72 / 3;
		875	AV *av = newAV ();
		876
		877	RETVAL = (SV )newRV_noinc ((SV )av);
		878	while (i--)
		879	{
		880	int h = *hv84++;
		881	int v1 = *hv84++;
		882	int v2 = *hv84++;
		883
		884	av_push (av, newSViv (v1));
		885	av_push (av, newSViv ((h >> 4) * 255 / 15));
		886	av_push (av, newSViv (v2));
		887	av_push (av, newSViv ((h & 15) * 255 / 15));
		888	}
		889	}
		890	OUTPUT:
		891	RETVAL
		892
211	#############################################################################	893	#############################################################################
212		894
213	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript	895	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
214		896
215	void	897	SV *
216	dump_pb (PerlIO fp, GdkPixbuf pb)	898	extract_features (SV *ar)
217	CODE:	899	CODE:
218	{	900	{
219	int w = gdk_pixbuf_get_width (pb);	901	int i;
220	int h = gdk_pixbuf_get_height (pb);	902	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		903
227	a85_init ();	904	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
		905	croak ("Not an array ref as first argument to extract_features");
228		906
229	for (y = 0; y < h; y++)	907	av = (AV *) SvRV (ar);
230	for (x = 0; x < w; x++)	908	result = newAV ();
231	for (i = 0; i < 3; i++)
232	a85_push (fp, src [x * bpp + y * sstr + i]);
233		909
234	a85_finish (fp);	910	for (i = 0; i <= av_len (av); ++i)
235	}	911	{
		912	SV sv = av_fetch (av, i, 1);
		913	SV histsv = newSV (9 sizeof (float) + 1);
236		914
		915	SvPOK_on (histsv);
		916	SvCUR_set (histsv, 9 * sizeof (float));
		917	float hist = (float )SvPVX (histsv);
237		918
		919	struct feature f_h, f_s, f_v;
		920	feature_init (&f_h);
		921	feature_init (&f_s);
		922	feature_init (&f_v);
238		923
		924	{
		925	STRLEN len;
		926	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		927	while (len >= 3)
		928	{
		929	unsigned int r, g, b, h, s, v;
		930	r = buf++; g = buf++; b = *buf++;
		931	rgb_to_hsv (r, g, b, &h, &s, &v);
239		932
		933	feature_update_pass_1 (&f_h, h);
		934	feature_update_pass_1 (&f_s, s);
		935	feature_update_pass_1 (&f_v, v);
		936
		937	len -= 3;
		938	}
		939
		940	feature_finish_pass_1 (&f_h);
		941	feature_finish_pass_1 (&f_s);
		942	feature_finish_pass_1 (&f_v);
		943	}
		944
		945	{
		946	STRLEN len;
		947	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		948	while (len >= 3)
		949	{
		950	unsigned int r, g, b, h, s, v;
		951	r = buf++; g = buf++; b = *buf++;
		952	rgb_to_hsv (r, g, b, &h, &s, &v);
		953
		954	feature_update_pass_2 (&f_h, h);
		955	feature_update_pass_2 (&f_s, s);
		956	feature_update_pass_2 (&f_v, v);
		957
		958	len -= 3;
		959	}
		960
		961	feature_finish_pass_2 (&f_h);
		962	feature_finish_pass_2 (&f_s);
		963	feature_finish_pass_2 (&f_v);
		964	}
		965
		966	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		967	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		968	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
		969
		970	av_push (result, histsv);
		971	}
		972
		973	RETVAL = newRV_noinc ((SV *)result);
		974	}
		975	OUTPUT:
		976	RETVAL
		977

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.57 by root, Mon Jul 30 22:55:46 2018 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.57 by root, Mon Jul 30 22:55:46 2018 UTC