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

Comparing CV/CV.xs (file contents):
Revision 1.18 by root, Wed Jul 20 09:06:37 2005 UTC vs.
Revision 1.57 by root, Mon Jul 30 22:55:46 2018 UTC

…		…
2	#include "perl.h"	2	#include "perl.h"
3	#include "XSUB.h"	3	#include "XSUB.h"
4		4
5	#include <string.h>	5	#include <string.h>
6	#include <setjmp.h>	6	#include <setjmp.h>
		7	#include <math.h>
		8
		9	#include <magic.h>
7		10
8	#include <jpeglib.h>	11	#include <jpeglib.h>
		12	#include <jerror.h>
		13
		14	#include <glib.h>
		15	#include <gtk/gtk.h>
		16	#include <gdk/gdkx.h>
9	#include <gdk-pixbuf/gdk-pixbuf.h>	17	#include <gdk-pixbuf/gdk-pixbuf.h>
10		18
11	#include <gperl.h>	19	#include <gperl.h>
12	#include <gtk2perl.h>	20	#include <gtk2perl.h>
13		21
		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
14	#define IW 80 /* MUST match Schnauer.pm! */	31	#define IW 80 /* MUST match Schnauzer.pm! */
15	#define IH 60 /* MUST match Schnauer.pm! */	32	#define IH 60 /* MUST match Schnauzer.pm! */
16		33
17	#define RAND (seed = (seed + 7141) * 54773 % 134456)	34	#define RAND (seed = (seed + 7141) * 54773 % 134456)
18		35
19	#define LINELENGTH 240	36	#define LINELENGTH 240
20		37
21	#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 */
22		43
23	struct jpg_err_mgr	44	struct jpg_err_mgr
24	{	45	{
25	struct jpeg_error_mgr err;	46	struct jpeg_error_mgr err;
26	jmp_buf setjmp_buffer;	47	jmp_buf setjmp_buffer;
27	};	48	};
28		49
29	static void	50	static void
30	cv_error_exit (j_common_ptr cinfo)	51	cv_error_exit (j_common_ptr cinfo)
31	{	52	{
		53	cinfo->err->output_message (cinfo);
32	longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);	54	longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
33	}	55	}
34		56
35	static void	57	static void
36	cv_error_output (j_common_ptr cinfo)	58	cv_error_output (j_common_ptr cinfo)
37	{	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);
38	return;	65	return;
39	}	66	}
40		67
41	static void	68	static void
42	rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,	69	rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
…		…
63		90
64	*h &= 255;	91	*h &= 255;
65	}	92	}
66	}	93	}
67		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	}
		149
68	static guint32 a85_val;	150	static guint32 a85_val;
69	static guint a85_cnt;	151	static guint a85_cnt;
70	static guchar a85_buf[LINELENGTH], *a85_ptr;	152	static guchar a85_buf[LINELENGTH], *a85_ptr;
71		153
72	static void	154	static void
…		…
84	if (!--a85_cnt)	166	if (!--a85_cnt)
85	{	167	{
86	a85_cnt = 4;	168	a85_cnt = 4;
87	if (a85_val)	169	if (a85_val)
88	{	170	{
89	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;	171	a85_ptr[4] = (a85_val % 85) + 33; a85_val /= 85;
90	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;	172	a85_ptr[3] = (a85_val % 85) + 33; a85_val /= 85;
91	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;	173	a85_ptr[2] = (a85_val % 85) + 33; a85_val /= 85;
92	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;	174	a85_ptr[1] = (a85_val % 85) + 33; a85_val /= 85;
93	a85_ptr[0] = (a85_val ) + 33;	175	a85_ptr[0] = (a85_val ) + 33;
94		176
95	a85_ptr += 5;	177	a85_ptr += 5;
96	}	178	}
…		…
119		201
120	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);	202	PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
121	}	203	}
122		204
123	/////////////////////////////////////////////////////////////////////////////	205	/////////////////////////////////////////////////////////////////////////////
		206	// memory source for libjpeg
		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	/////////////////////////////////////////////////////////////////////////////
124		261
125	MODULE = Gtk2::CV PACKAGE = Gtk2::CV	262	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
126		263
127	PROTOTYPES: ENABLE	264	PROTOTYPES: ENABLE
128		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
129	GdkPixbuf_noinc *	333	GdkPixbuf_noinc *
130	transpose (GdkPixbuf *pb)	334	dealpha_expose (GdkPixbuf *pb)
131	CODE:	335	CODE:
		336	perlinterp_release ();
132	{	337	{
133	int w = gdk_pixbuf_get_width (pb);	338	int w = gdk_pixbuf_get_width (pb);
134	int h = gdk_pixbuf_get_height (pb);	339	int h = gdk_pixbuf_get_height (pb);
135	int bpp = gdk_pixbuf_get_n_channels (pb);	340	int bpp = gdk_pixbuf_get_n_channels (pb);
136	int x, y, i;	341	int x, y, i;
137	guchar src = gdk_pixbuf_get_pixels (pb), dst;	342	guchar src = gdk_pixbuf_get_pixels (pb), dst;
138	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	343	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;
139		344
140	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, h, w);	345	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
141		346
142	dst = gdk_pixbuf_get_pixels (RETVAL);	347	dst = gdk_pixbuf_get_pixels (RETVAL);
143	dstr = gdk_pixbuf_get_rowstride (RETVAL);	348	dstr = gdk_pixbuf_get_rowstride (RETVAL);
144		349
145	for (y = 0; y < h; y++)
146	for (x = 0; x < w; x++)	350	for (x = 0; x < w; x++)
		351	for (y = 0; y < h; y++)
147	for (i = 0; i < bpp; i++)	352	for (i = 0; i < 3; i++)
148	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];
149	}	354	}
		355	perlinterp_acquire ();
150	OUTPUT:	356	OUTPUT:
151	RETVAL	357	RETVAL
152		358
153	GdkPixbuf_noinc *	359	GdkPixbuf_noinc *
154	flop (GdkPixbuf *pb)	360	rotate (GdkPixbuf *pb, int angle)
155	CODE:	361	CODE:
156	{	362	perlinterp_release ();
157	int w = gdk_pixbuf_get_width (pb);	363	if (angle < 0)
158	int h = gdk_pixbuf_get_height (pb);	364	angle += 360;
159	int bpp = gdk_pixbuf_get_n_channels (pb);	365	RETVAL = gdk_pixbuf_rotate_simple (pb, angle == 0 ? GDK_PIXBUF_ROTATE_NONE
160	int x, y, i;	366	: angle == 90 ? GDK_PIXBUF_ROTATE_COUNTERCLOCKWISE
161	guchar src = gdk_pixbuf_get_pixels (pb), dst;	367	: angle == 180 ? GDK_PIXBUF_ROTATE_UPSIDEDOWN
162	int sstr = gdk_pixbuf_get_rowstride (pb), dstr;	368	: angle == 270 ? GDK_PIXBUF_ROTATE_CLOCKWISE
163		369	: angle);
164	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, bpp == 4, 8, w, h);	370	perlinterp_acquire ();
165
166	dst = gdk_pixbuf_get_pixels (RETVAL);
167	dstr = gdk_pixbuf_get_rowstride (RETVAL);
168
169	for (y = 0; y < h; y++)
170	for (x = 0; x < w; x++)
171	for (i = 0; i < bpp; i++)
172	dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
173	}
174	OUTPUT:	371	OUTPUT:
175	RETVAL	372	RETVAL
176		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
177	GdkPixbuf_noinc *	412	GdkPixbuf_noinc *
178	load_jpeg (char *path, int thumbnail=0)	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)
179	CODE:	503	CODE:
180	{	504	{
181	struct jpeg_decompress_struct cinfo;	505	struct jpeg_decompress_struct cinfo;
182	struct jpg_err_mgr jerr;	506	struct jpg_err_mgr jerr;
183	guchar *data;
184	int rs;	507	int rs;
185	FILE *fp;
186	volatile GdkPixbuf *pb = 0;	508	volatile GdkPixbuf *pb = 0;
		509	STRLEN data_len;
		510	guchar *data = SvPVbyte (image_data, data_len);
		511
187	RETVAL = 0;	512	RETVAL = 0;
188		513
189	if (!(fp = fopen (path, "rb")))	514	perlinterp_release ();
190	XSRETURN_UNDEF;
191		515
192	cinfo.err = jpeg_std_error (&jerr.err);	516	cinfo.err = jpeg_std_error (&jerr.err);
193		517
194	jerr.err.error_exit = cv_error_exit;	518	jerr.err.error_exit = cv_error_exit;
195	jerr.err.output_message = cv_error_output;	519	jerr.err.output_message = cv_error_output;
196		520
197	if ((rs = setjmp (jerr.setjmp_buffer)))	521	if ((rs = setjmp (jerr.setjmp_buffer)))
198	{	522	{
199	fclose (fp);
200	jpeg_destroy_decompress (&cinfo);	523	jpeg_destroy_decompress (&cinfo);
201		524
202	if (pb)	525	if (pb)
203	g_object_unref ((gpointer)pb);	526	g_object_unref ((gpointer)pb);
204		527
		528	perlinterp_acquire ();
205	XSRETURN_UNDEF;	529	XSRETURN_UNDEF;
206	}	530	}
207		531
		532	if (!data_len)
		533	longjmp (jerr.setjmp_buffer, 4);
		534
208	jpeg_create_decompress (&cinfo);	535	jpeg_create_decompress (&cinfo);
		536	cv_jpeg_mem_src (&cinfo, data, data_len);
209		537
210	jpeg_stdio_src (&cinfo, fp);
211	jpeg_read_header (&cinfo, TRUE);	538	jpeg_read_header (&cinfo, TRUE);
212		539
213	cinfo.dct_method = JDCT_DEFAULT;	540	cinfo.dct_method = JDCT_DEFAULT;
214	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */	541	cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
215	cinfo.do_block_smoothing = FALSE;	542	cinfo.do_block_smoothing = FALSE;
…		…
225	{	552	{
226	cinfo.dct_method = JDCT_FASTEST;	553	cinfo.dct_method = JDCT_FASTEST;
227	cinfo.do_fancy_upsampling = FALSE;	554	cinfo.do_fancy_upsampling = FALSE;
228		555
229	while (cinfo.scale_denom < 8	556	while (cinfo.scale_denom < 8
230	&& (cinfo.output_width >> 1) >= IW	557	&& cinfo.output_width >= iw*4
231	&& (cinfo.output_height >> 1) >= IH)	558	&& cinfo.output_height >= ih*4)
232	{	559	{
233	cinfo.scale_denom <<= 1;	560	cinfo.scale_denom <<= 1;
234	jpeg_calc_output_dimensions (&cinfo);	561	jpeg_calc_output_dimensions (&cinfo);
235	}	562	}
236	}	563	}
237		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
238	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height);	574	pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, cinfo.output_components == 4, 8, cinfo.output_width, cinfo.output_height);
239	if (!RETVAL)	575	if (!RETVAL)
240	longjmp (jerr.setjmp_buffer, 2);	576	longjmp (jerr.setjmp_buffer, 2);
241		577
242	data = gdk_pixbuf_get_pixels (RETVAL);	578	data = gdk_pixbuf_get_pixels (RETVAL);
243	rs = gdk_pixbuf_get_rowstride (RETVAL);	579	rs = gdk_pixbuf_get_rowstride (RETVAL);
244
245	if (cinfo.output_components != 3)
246	longjmp (jerr.setjmp_buffer, 3);
247		580
248	jpeg_start_decompress (&cinfo);	581	jpeg_start_decompress (&cinfo);
249		582
250	while (cinfo.output_scanline < cinfo.output_height)	583	while (cinfo.output_scanline < cinfo.output_height)
251	{	584	{
…		…
258	rp [3] = (guchar *)rp [2] + rs;	591	rp [3] = (guchar *)rp [2] + rs;
259		592
260	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);	593	jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
261	}	594	}
262		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
263	jpeg_finish_decompress (&cinfo);	625	jpeg_finish_decompress (&cinfo);
264	fclose (fp);
265	jpeg_destroy_decompress (&cinfo);	626	jpeg_destroy_decompress (&cinfo);
		627	perlinterp_acquire ();
266	}	628	}
267	OUTPUT:	629	OUTPUT:
268	RETVAL	630	RETVAL
269		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
270	#############################################################################	680	#############################################################################
271		681
272	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	682	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		683
		684	# currently only works for filenames (octet strings)
273		685
274	SV *	686	SV *
275	foldcase (SV *pathsv)	687	foldcase (SV *pathsv)
276	PROTOTYPE: $	688	PROTOTYPE: $
277	CODE:	689	CODE:
278	{	690	{
279	STRLEN plen;	691	STRLEN plen;
280	U8 *path = SvPVutf8 (pathsv, plen);	692	U8 path = (U8 )SvPV (pathsv, plen);
281	U8 *pend = path + plen;	693	U8 *pend = path + plen;
282	U8 dst [plen * 6 * 3], *dstp = dst;	694	U8 dst [plen * 6 * 3], *dstp = dst;
283		695
284	while (path < pend)	696	while (path < pend)
285	{	697	{
		698	U8 ch = *path;
		699
286	if (path >= '0' && path <= '9')	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')
287	{	705	{
288	STRLEN el, nl = 0;	706	STRLEN el, nl = 0;
289	while (path >= '0' && path <= '9')	707	while (path >= '0' && path <= '9' && path < pend)
290	path++, nl++;	708	path++, nl++;
291		709
292	for (el = nl; el < 6; el++)	710	for (el = nl; el < 6; el++)
293	*dstp++ = '0';	711	*dstp++ = '0';
294		712
295	memcpy (dstp, path - nl, nl);	713	memcpy (dstp, path - nl, nl);
296	dstp += nl;	714	dstp += nl;
297	}	715	}
		716	else
		717	dstp++ = path++;
		718	#if 0
298	else	719	else
299	{	720	{
300	STRLEN cl;	721	STRLEN cl;
301	to_utf8_fold (path, dstp, &cl);	722	to_utf8_fold (path, dstp, &cl);
302	dstp += cl;	723	dstp += cl;
303	path += is_utf8_char (path);	724	path += is_utf8_char (path);
304	}	725	}
		726	#endif
305	}	727	}
306		728
307	RETVAL = newSVpvn (dst, dstp - dst);	729	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
308	}	730	}
309	OUTPUT:	731	OUTPUT:
310	RETVAL	732	RETVAL
311		733
312	GdkPixbuf_noinc *	734	GdkPixbuf_noinc *
313	p7_to_pb (int w, int h, guchar *src)	735	p7_to_pb (int w, int h, SV *src_sv)
		736	PROTOTYPE: @
314	CODE:	737	CODE:
315	{	738	{
316	int x, y;	739	int x, y;
317	guchar dst, d;	740	guchar dst, d;
318	int dstr;	741	int dstr;
		742	guchar src = (guchar )SvPVbyte_nolen (src_sv);
319		743
320	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);	744	RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, w, h);
321	dst = gdk_pixbuf_get_pixels (RETVAL);	745	dst = gdk_pixbuf_get_pixels (RETVAL);
322	dstr = gdk_pixbuf_get_rowstride (RETVAL);	746	dstr = gdk_pixbuf_get_rowstride (RETVAL);
323		747
…		…
332	}	756	}
333	}	757	}
334	OUTPUT:	758	OUTPUT:
335	RETVAL	759	RETVAL
336		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
337	SV *	809	SV *
338	pb_to_p7 (GdkPixbuf *pb)	810	pb_to_hv84 (GdkPixbuf *pb)
339	CODE:	811	CODE:
340	{	812	{
341	int w = gdk_pixbuf_get_width (pb);	813	int w = gdk_pixbuf_get_width (pb);
342	int h = gdk_pixbuf_get_height (pb);	814	int h = gdk_pixbuf_get_height (pb);
343	int x, y;	815	int x, y;
344	guchar *dst;	816	guchar *dst;
345	int bpp = gdk_pixbuf_get_n_channels (pb);	817	int bpp = gdk_pixbuf_get_n_channels (pb);
346	guchar *src = gdk_pixbuf_get_pixels (pb);	818	guchar *src = gdk_pixbuf_get_pixels (pb);
347	int sstr = gdk_pixbuf_get_rowstride (pb);	819	int sstr = gdk_pixbuf_get_rowstride (pb);
348	int Er[IW], Eg[IW], Eb[IW];
349	int seed = 77;
350
351	RETVAL = newSV (w * h);
352	SvPOK_only (RETVAL);
353	SvCUR_set (RETVAL, w * h);
354
355	dst = SvPVX (RETVAL);
356
357	memset (Er, 0, sizeof (int) * IW);
358	memset (Eg, 0, sizeof (int) * IW);
359	memset (Eb, 0, sizeof (int) * IW);
360
361	/* some primitive error distribution + random dithering */
362
363	for (y = 0; y < h; y++)
364	{
365	int er = 0, eg = 0, eb = 0;
366
367	for (x = 0; x < w; x++)
368	{
369	int r, g, b;
370	guchar p = src + x bpp + y * sstr;
371
372	r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255;
373	g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255;
374	b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255;
375
376	r = r > 7 ? 7 : r < 0 ? 0 : r;
377	g = g > 7 ? 7 : g < 0 ? 0 : g;
378	b = b > 3 ? 3 : b < 0 ? 0 : b;
379
380	er += p[0] - (r * 255 + 4) / 7;
381	eg += p[1] - (g * 255 + 4) / 7;
382	eb += p[2] - (b * 255 + 2) / 3;
383
384	Er[x] = er >> 1; er -= (er + 1) >> 1;
385	Eg[x] = eg >> 1; eg -= (eg + 1) >> 1;
386	Eb[x] = eb >> 1; eb -= (eb + 1) >> 1;
387
388	*dst++ = r << 5 \| g << 2 \| b;
389	}
390	}
391	}
392	OUTPUT:
393	RETVAL
394
395	SV *
396	make_histogram (SV *ar)
397	CODE:
398	{
399	int i;
400	AV av, result;
401
402	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
403	croak ("Not an array ref as first argument to make_histogram");
404
405	av = (AV *) SvRV (ar);
406	result = newAV ();
407
408	for (i = 0; i <= av_len (av); ++i)
409	{
410	const int HISTSIZE = 64;
411
412	int j;
413	SV sv = av_fetch (av, i, 1);
414	STRLEN len;
415	char *buf = SvPVbyte (sv, len);
416
417	int tmphist[HISTSIZE];
418	float *hist;
419
420	SV histsv = newSV (HISTSIZE sizeof (float) + 1);
421	SvPOK_on (histsv);
422	SvCUR_set (histsv, HISTSIZE * sizeof (float));
423	hist = (float *)SvPVX (histsv);
424
425	Zero (tmphist, sizeof (tmphist), char);
426
427	for (j = len; j--; )
428	{
429	unsigned int idx
430	= ((*buf & 0xc0) >> 2)
431	\| ((*buf & 0x18) >> 1)
432	\| (*buf & 0x03);
433
434	++tmphist[idx];
435	++buf;
436	}
437
438	for (j = 0; j < HISTSIZE; ++j)
439	hist[j] = (float)tmphist[j] / (len + 1e-30);
440
441	av_push (result, histsv);
442	}
443
444	RETVAL = newRV_noinc ((SV *)result);
445	}
446	OUTPUT:
447	RETVAL
448
449	#############################################################################
450
451	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
452
453	void
454	dump_ascii85 (PerlIO fp, GdkPixbuf pb)
455	CODE:
456	{
457	int w = gdk_pixbuf_get_width (pb);
458	int h = gdk_pixbuf_get_height (pb);
459	int x, y, i;
460	guchar *dst;
461	int bpp = gdk_pixbuf_get_n_channels (pb);
462	guchar *src = gdk_pixbuf_get_pixels (pb);
463	int sstr = gdk_pixbuf_get_rowstride (pb);
464
465	a85_init ();
466
467	for (y = 0; y < h; y++)
468	for (x = 0; x < w; x++)
469	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
470	a85_push (fp, src [x * bpp + y * sstr + i]);
471
472	a85_finish (fp);
473	}
474
475	void
476	dump_binary (PerlIO fp, GdkPixbuf pb)
477	CODE:
478	{
479	int w = gdk_pixbuf_get_width (pb);
480	int h = gdk_pixbuf_get_height (pb);
481	int x, y, i;
482	guchar *dst;
483	int bpp = gdk_pixbuf_get_n_channels (pb);
484	guchar *src = gdk_pixbuf_get_pixels (pb);
485	int sstr = gdk_pixbuf_get_rowstride (pb);
486
487	for (y = 0; y < h; y++)
488	for (x = 0; x < w; x++)
489	for (i = 0; i < (bpp < 3 ? 1 : 3); i++)
490	PerlIO_putc (fp, src [x * bpp + y * sstr + i]);
491	}
492
493	#############################################################################
494
495	MODULE = Gtk2::CV PACKAGE = Gtk2::CV
496
497	SV *
498	pb_to_hv84 (GdkPixbuf *pb)
499	CODE:
500	{
501	int w = gdk_pixbuf_get_width (pb);
502	int h = gdk_pixbuf_get_height (pb);
503	int x, y;
504	guchar *dst;
505	int bpp = gdk_pixbuf_get_n_channels (pb);
506	guchar *src = gdk_pixbuf_get_pixels (pb);
507	int sstr = gdk_pixbuf_get_rowstride (pb);
508		820
509	RETVAL = newSV (6 * 8 * 12 / 8);	821	RETVAL = newSV (6 * 8 * 12 / 8);
510	SvPOK_only (RETVAL);	822	SvPOK_only (RETVAL);
511	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);	823	SvCUR_set (RETVAL, 6 * 8 * 12 / 8);
512		824
513	dst = SvPVX (RETVAL);	825	dst = (guchar *)SvPVX (RETVAL);
514		826
515	/* some primitive error distribution + random dithering */	827	/* some primitive error distribution + random dithering */
516		828
517	for (y = 0; y < h; y++)	829	for (y = 0; y < h; y++)
518	{	830	{
…		…
576	}	888	}
577	}	889	}
578	OUTPUT:	890	OUTPUT:
579	RETVAL	891	RETVAL
580		892
		893	#############################################################################
581		894
		895	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
582		896
		897	SV *
		898	extract_features (SV *ar)
		899	CODE:
		900	{
		901	int i;
		902	AV av, result;
		903
		904	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
		905	croak ("Not an array ref as first argument to extract_features");
		906
		907	av = (AV *) SvRV (ar);
		908	result = newAV ();
		909
		910	for (i = 0; i <= av_len (av); ++i)
		911	{
		912	SV sv = av_fetch (av, i, 1);
		913	SV histsv = newSV (9 sizeof (float) + 1);
		914
		915	SvPOK_on (histsv);
		916	SvCUR_set (histsv, 9 * sizeof (float));
		917	float hist = (float )SvPVX (histsv);
		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);
		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);
		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.18 by root, Wed Jul 20 09:06:37 2005 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.18 by root, Wed Jul 20 09:06:37 2005 UTC vs.
Revision 1.57 by root, Mon Jul 30 22:55:46 2018 UTC