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

Comparing CV/CV.xs (file contents):
Revision 1.32 by root, Wed Feb 15 08:36:55 2006 UTC vs.
Revision 1.41 by root, Sat Mar 1 04:11:26 2008 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>
11	#include <glib.h>	12	#include <glib.h>
…		…
13	#include <gdk-pixbuf/gdk-pixbuf.h>	14	#include <gdk-pixbuf/gdk-pixbuf.h>
14		15
15	#include <gperl.h>	16	#include <gperl.h>
16	#include <gtk2perl.h>	17	#include <gtk2perl.h>
17		18
18	#define IW 80 /* MUST match Schnauer.pm! */	19	#define IW 80 /* MUST match Schnauzer.pm! */
19	#define IH 60 /* MUST match Schnauer.pm! */	20	#define IH 60 /* MUST match Schnauzer.pm! */
20		21
21	#define RAND (seed = (seed + 7141) * 54773 % 134456)	22	#define RAND (seed = (seed + 7141) * 54773 % 134456)
22		23
23	#define LINELENGTH 240	24	#define LINELENGTH 240
24		25
…		…
67		68
68	*h &= 255;	69	*h &= 255;
69	}	70	}
70	}	71	}
71		72
		73	struct feature {
		74	float v1, v2, v3; // mean, square, cube
		75	int n;
		76	};
		77
		78	static void
		79	feature_init (struct feature *f)
		80	{
		81	f->v1 = 0.;
		82	f->v2 = 0.;
		83	f->v3 = 0.;
		84	f->n = 0;
		85	}
		86
		87	// didn't find an algorithm to neatly do mean, variance and skew in one pass.
		88	// elmex ist schuld.
		89	static void
		90	feature_update_pass_1 (struct feature *f, unsigned int v)
		91	{
		92	f->v1 += v;
		93	f->n += 1;
		94	}
		95
		96	static void
		97	feature_finish_pass_1 (struct feature *f)
		98	{
		99	if (f->n < 1)
		100	return;
		101
		102	f->v1 /= f->n;
		103	}
		104
		105	static void
		106	feature_update_pass_2 (struct feature *f, unsigned int v)
		107	{
		108	float d = v - f->v1;
		109
		110	f->v2 += d * d;
		111	f->v3 += d * d * d;
		112	}
		113
		114	static void
		115	feature_finish_pass_2 (struct feature *f)
		116	{
		117	if (f->n < 1)
		118	return;
		119
		120	f->v2 /= f->n;
		121	f->v3 /= f->n;
		122
		123	f->v1 /= 255.;
		124	f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2);
		125	f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.);
		126	}
		127
72	static guint32 a85_val;	128	static guint32 a85_val;
73	static guint a85_cnt;	129	static guint a85_cnt;
74	static guchar a85_buf[LINELENGTH], *a85_ptr;	130	static guchar a85_buf[LINELENGTH], *a85_ptr;
75		131
76	static void	132	static void
…		…
187	RETVAL = magic_file (cookie, path);	243	RETVAL = magic_file (cookie, path);
188	}	244	}
189	OUTPUT:	245	OUTPUT:
190	RETVAL	246	RETVAL
191		247
192	# missing in Gtk2 perl module	248	# missing/broken in Gtk2 perl module
		249
		250	void
		251	gdk_window_clear_hints (GdkWindow *window)
		252	CODE:
		253	gdk_window_set_geometry_hints (window, 0, 0);
193		254
194	gboolean	255	gboolean
195	gdk_net_wm_supports (GdkAtom property)	256	gdk_net_wm_supports (GdkAtom property)
196	CODE:	257	CODE:
197	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)	258	#if defined(GDK_WINDOWING_X11) && !defined(GDK_MULTIHEAD_SAFE)
…		…
331	jpeg_destroy_decompress (&cinfo);	392	jpeg_destroy_decompress (&cinfo);
332	}	393	}
333	OUTPUT:	394	OUTPUT:
334	RETVAL	395	RETVAL
335		396
		397	void
		398	compare (GdkPixbuf a, GdkPixbuf b)
		399	PPCODE:
		400	{
		401	int w = gdk_pixbuf_get_width (a);
		402	int h = gdk_pixbuf_get_height (a);
		403	int sa = gdk_pixbuf_get_rowstride (a);
		404	int sb = gdk_pixbuf_get_rowstride (b);
		405
		406	guchar *pa = gdk_pixbuf_get_pixels (a);
		407	guchar *pb = gdk_pixbuf_get_pixels (b);
		408
		409	int x, y;
		410
		411	assert (w == gdk_pixbuf_get_width (b));
		412	assert (h == gdk_pixbuf_get_height (b));
		413
		414	assert (gdk_pixbuf_get_n_channels (a) == 3);
		415	assert (gdk_pixbuf_get_n_channels (b) == 3);
		416
		417	double diff = 0.;
		418	int peak = 0;
		419
		420	if (w && h)
		421	for (y = 0; y < h; y++)
		422	{
		423	guchar pa_ = pa + y sa;
		424	guchar pb_ = pb + y sb;
		425
		426	for (x = 0; x < w; x++)
		427	{
		428	int d;
		429
		430	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		431	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		432	d = ((int)pa_++) - ((int)pb_++); diff += d*d; peak = MAX (peak, abs (d));
		433	}
		434	}
		435
		436	EXTEND (SP, 2);
		437	PUSHs (sv_2mortal (newSVnv (sqrt (diff / (w * h * 3. * 255. * 255.)))));
		438	PUSHs (sv_2mortal (newSVnv (peak / 255.)));
		439	}
		440
336	#############################################################################	441	#############################################################################
337		442
338	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer	443	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
		444
		445	# currently only works for filenames (octet strings)
339		446
340	SV *	447	SV *
341	foldcase (SV *pathsv)	448	foldcase (SV *pathsv)
342	PROTOTYPE: $	449	PROTOTYPE: $
343	CODE:	450	CODE:
344	{	451	{
345	STRLEN plen;	452	STRLEN plen;
346	U8 path = (U8 )SvPVutf8 (pathsv, plen);	453	U8 path = (U8 )SvPV (pathsv, plen);
347	U8 *pend = path + plen;	454	U8 *pend = path + plen;
348	U8 dst [plen * 6 * 3], *dstp = dst;	455	U8 dst [plen * 6 * 3], *dstp = dst;
349		456
350	while (path < pend)	457	while (path < pend)
351	{	458	{
352	U8 ch = *path;	459	U8 ch = *path;
353		460
354	if (ch >= 'a' && ch <= 'z')	461	if (ch >= 'a' && ch <= 'z')
355	dstp++ = path++;	462	dstp++ = path++;
		463	else if (ch >= 'A' && ch <= 'Z')
		464	dstp++ = path++ + ('a' - 'A');
356	else if (ch >= '0' && ch <= '9')	465	else if (ch >= '0' && ch <= '9')
357	{	466	{
358	STRLEN el, nl = 0;	467	STRLEN el, nl = 0;
359	while (path >= '0' && path <= '9' && path < pend)	468	while (path >= '0' && path <= '9' && path < pend)
360	path++, nl++;	469	path++, nl++;
…		…
363	*dstp++ = '0';	472	*dstp++ = '0';
364		473
365	memcpy (dstp, path - nl, nl);	474	memcpy (dstp, path - nl, nl);
366	dstp += nl;	475	dstp += nl;
367	}	476	}
		477	else
		478	dstp++ = path++;
		479	#if 0
368	else	480	else
369	{	481	{
370	STRLEN cl;	482	STRLEN cl;
371	to_utf8_fold (path, dstp, &cl);	483	to_utf8_fold (path, dstp, &cl);
372	dstp += cl;	484	dstp += cl;
373	path += is_utf8_char (path);	485	path += is_utf8_char (path);
374	}	486	}
		487	#endif
375	}	488	}
376		489
377	RETVAL = newSVpvn ((const char *)dst, dstp - dst);	490	RETVAL = newSVpvn ((const char *)dst, dstp - dst);
378	}	491	}
379	OUTPUT:	492	OUTPUT:
…		…
541	#############################################################################	654	#############################################################################
542		655
543	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster	656	MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
544		657
545	SV *	658	SV *
546	make_histograms (SV *ar)	659	extract_features (SV *ar)
547	CODE:	660	CODE:
548	{	661	{
549	int i;	662	int i;
550	AV av, result;	663	AV av, result;
551		664
552	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)	665	if (!SvROK (ar) \|\| SvTYPE (SvRV (ar)) != SVt_PVAV)
553	croak ("Not an array ref as first argument to make_histogram");	666	croak ("Not an array ref as first argument to extract_features");
554		667
555	av = (AV *) SvRV (ar);	668	av = (AV *) SvRV (ar);
556	result = newAV ();	669	result = newAV ();
557		670
558	for (i = 0; i <= av_len (av); ++i)	671	for (i = 0; i <= av_len (av); ++i)
559	{	672	{
560	const int HISTSIZE = 64;
561
562	int j;
563	SV sv = av_fetch (av, i, 1);	673	SV sv = av_fetch (av, i, 1);
564	STRLEN len;
565	char *buf = SvPVbyte (sv, len);
566
567	int tmphist[HISTSIZE];
568	float *hist;
569
570	SV histsv = newSV (HISTSIZE sizeof (float) + 1);	674	SV histsv = newSV (9 sizeof (float) + 1);
		675
571	SvPOK_on (histsv);	676	SvPOK_on (histsv);
572	SvCUR_set (histsv, HISTSIZE * sizeof (float));	677	SvCUR_set (histsv, 9 * sizeof (float));
573	hist = (float *)SvPVX (histsv);	678	float hist = (float )SvPVX (histsv);
574		679
575	Zero (tmphist, sizeof (tmphist), char);	680	struct feature f_h, f_s, f_v;
		681	feature_init (&f_h);
		682	feature_init (&f_s);
		683	feature_init (&f_v);
576		684
577	for (j = len; j--; )
578	{
579	unsigned int idx
580	= ((*buf & 0xc0) >> 2)
581	\| ((*buf & 0x18) >> 1)
582	\| (*buf & 0x03);
583
584	++tmphist[idx];
585	++buf;
586	}
587		685	{
588	for (j = 0; j < HISTSIZE; ++j)	686	STRLEN len;
589	hist[j] = (float)tmphist[j] / (len + 1e-30);	687	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		688	while (len >= 3)
		689	{
		690	unsigned int r, g, b, h, s, v;
		691	r = buf++; g = buf++; b = *buf++;
		692	rgb_to_hsv (r, g, b, &h, &s, &v);
		693
		694	feature_update_pass_1 (&f_h, h);
		695	feature_update_pass_1 (&f_s, s);
		696	feature_update_pass_1 (&f_v, v);
		697
		698	len -= 3;
		699	}
		700
		701	feature_finish_pass_1 (&f_h);
		702	feature_finish_pass_1 (&f_s);
		703	feature_finish_pass_1 (&f_v);
		704	}
		705
		706	{
		707	STRLEN len;
		708	unsigned char buf = (unsigned char )SvPVbyte (sv, len);
		709	while (len >= 3)
		710	{
		711	unsigned int r, g, b, h, s, v;
		712	r = buf++; g = buf++; b = *buf++;
		713	rgb_to_hsv (r, g, b, &h, &s, &v);
		714
		715	feature_update_pass_2 (&f_h, h);
		716	feature_update_pass_2 (&f_s, s);
		717	feature_update_pass_2 (&f_v, v);
		718
		719	len -= 3;
		720	}
		721
		722	feature_finish_pass_2 (&f_h);
		723	feature_finish_pass_2 (&f_s);
		724	feature_finish_pass_2 (&f_v);
		725	}
		726
		727	hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.;
		728	hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ;
		729	hist [6] = f_v.v1 * .5; hist [7] = f_v.v2 * .5; hist [8] = f_v.v3 * .5;
590		730
591	av_push (result, histsv);	731	av_push (result, histsv);
592	}	732	}
593		733
594	RETVAL = newRV_noinc ((SV *)result);	734	RETVAL = newRV_noinc ((SV *)result);
595	}	735	}
596	OUTPUT:	736	OUTPUT:
597	RETVAL	737	RETVAL
598		738
599

Diff Legend

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

Comparing CV/CV.xs (file contents): Revision 1.32 by root, Wed Feb 15 08:36:55 2006 UTC vs. Revision 1.41 by root, Sat Mar 1 04:11:26 2008 UTC

Diff Legend

Comparing CV/CV.xs (file contents):
Revision 1.32 by root, Wed Feb 15 08:36:55 2006 UTC vs.
Revision 1.41 by root, Sat Mar 1 04:11:26 2008 UTC