ViewVC Help
View File | Revision Log | Show Annotations | Download File
/cvs/CV/CV.xs
(Generate patch)

Comparing CV/CV.xs (file contents):
Revision 1.12 by root, Tue Feb 8 11:00:11 2005 UTC vs.
Revision 1.21 by root, Sat Jul 23 03:25:51 2005 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
8#include <jpeglib.h>
5#include <gdk-pixbuf/gdk-pixbuf.h> 9#include <gdk-pixbuf/gdk-pixbuf.h>
6 10
7#include <gperl.h> 11#include <gperl.h>
8#include <gtk2perl.h> 12#include <gtk2perl.h>
9 13
10#define IW 80 14#define IW 80 /* MUST match Schnauer.pm! */
15#define IH 60 /* MUST match Schnauer.pm! */
11 16
12#define RAND (seed = (seed + 7141) * 54773 % 134456) 17#define RAND (seed = (seed + 7141) * 54773 % 134456)
13 18
14#define LINELENGTH 240 19#define LINELENGTH 240
15 20
16#define ELLIPSIS "\xe2\x80\xa6" 21#define ELLIPSIS "\xe2\x80\xa6"
22
23struct jpg_err_mgr
24{
25 struct jpeg_error_mgr err;
26 jmp_buf setjmp_buffer;
27};
28
29static void
30cv_error_exit (j_common_ptr cinfo)
31{
32 longjmp (((struct jpg_err_mgr *)cinfo->err)->setjmp_buffer, 99);
33}
34
35static void
36cv_error_output (j_common_ptr cinfo)
37{
38 return;
39}
40
41static void
42rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
43 unsigned int *h, unsigned int *s, unsigned int *v)
44{
45 unsigned int mx = r; if (g > mx) mx = g; if (b > mx) mx = b;
46 unsigned int mn = r; if (g < mn) mn = g; if (b < mn) mn = b;
47 unsigned int delta = mx - mn;
48
49 *v = mx;
50
51 *s = mx ? delta * 255 / mx : 0;
52
53 if (delta == 0)
54 *h = 0;
55 else
56 {
57 if (r == mx)
58 *h = ((int)g - (int)b) * 255 / (int)(delta * 3);
59 else if (g == mx)
60 *h = ((int)b - (int)r) * 255 / (int)(delta * 3) + 52;
61 else if (b == mx)
62 *h = ((int)r - (int)g) * 255 / (int)(delta * 3) + 103;
63
64 *h &= 255;
65 }
66}
17 67
18static guint32 a85_val; 68static guint32 a85_val;
19static guint a85_cnt; 69static guint a85_cnt;
20static guchar a85_buf[LINELENGTH], *a85_ptr; 70static guchar a85_buf[LINELENGTH], *a85_ptr;
21 71
68 *a85_ptr++ = '\n'; 118 *a85_ptr++ = '\n';
69 119
70 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf); 120 PerlIO_write (fp, a85_buf, a85_ptr - a85_buf);
71} 121}
72 122
73static void
74rgb_to_hsv (unsigned int r, unsigned int g, unsigned int b,
75 unsigned int *h, unsigned int *s, unsigned int *v)
76{
77 unsigned int mx = r; if (g > mx) mx = g; if (b > mx) mx = b;
78 unsigned int mn = r; if (g < mn) mn = g; if (b < mn) mn = b;
79 unsigned int delta = mx - mn;
80
81 *v = mx;
82
83 *s = mx ? delta * 255 / mx : 0;
84
85 if (delta == 0)
86 *h = 0;
87 else
88 {
89 if (r == mx)
90 *h = ((int)g - (int)b) * 255 / (int)(delta * 3);
91 else if (g == mx)
92 *h = ((int)b - (int)r) * 255 / (int)(delta * 3) + 52;
93 else if (b == mx)
94 *h = ((int)r - (int)g) * 255 / (int)(delta * 3) + 103;
95
96 *h &= 255;
97 }
98}
99
100///////////////////////////////////////////////////////////////////////////// 123/////////////////////////////////////////////////////////////////////////////
101 124
102MODULE = Gtk2::CV PACKAGE = Gtk2::CV 125MODULE = Gtk2::CV PACKAGE = Gtk2::CV
103 126
104PROTOTYPES: ENABLE 127PROTOTYPES: ENABLE
149 dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i]; 172 dst[(w - 1 - x) * bpp + y * dstr + i] = src[x * bpp + y * sstr + i];
150} 173}
151 OUTPUT: 174 OUTPUT:
152 RETVAL 175 RETVAL
153 176
177GdkPixbuf_noinc *
178load_jpeg (char *path, int thumbnail=0)
179 CODE:
180{
181 struct jpeg_decompress_struct cinfo;
182 struct jpg_err_mgr jerr;
183 guchar *data;
184 int rs;
185 FILE *fp;
186 volatile GdkPixbuf *pb = 0;
187 RETVAL = 0;
188
189 if (!(fp = fopen (path, "rb")))
190 XSRETURN_UNDEF;
191
192 cinfo.err = jpeg_std_error (&jerr.err);
193
194 jerr.err.error_exit = cv_error_exit;
195 jerr.err.output_message = cv_error_output;
196
197 if ((rs = setjmp (jerr.setjmp_buffer)))
198 {
199 fclose (fp);
200 jpeg_destroy_decompress (&cinfo);
201
202 if (pb)
203 g_object_unref ((gpointer)pb);
204
205 XSRETURN_UNDEF;
206 }
207
208 jpeg_create_decompress (&cinfo);
209
210 jpeg_stdio_src (&cinfo, fp);
211 jpeg_read_header (&cinfo, TRUE);
212
213 cinfo.dct_method = JDCT_DEFAULT;
214 cinfo.do_fancy_upsampling = FALSE; /* worse quality, but nobody compained so far, and gdk-pixbuf does the same */
215 cinfo.do_block_smoothing = FALSE;
216 cinfo.out_color_space = JCS_RGB;
217 cinfo.quantize_colors = FALSE;
218
219 cinfo.scale_num = 1;
220 cinfo.scale_denom = 1;
221
222 jpeg_calc_output_dimensions (&cinfo);
223
224 if (thumbnail)
225 {
226 cinfo.dct_method = JDCT_FASTEST;
227 cinfo.do_fancy_upsampling = FALSE;
228
229 while (cinfo.scale_denom < 8
230 && cinfo.output_width >= IW*4
231 && cinfo.output_height >= IH*4)
232 {
233 cinfo.scale_denom <<= 1;
234 jpeg_calc_output_dimensions (&cinfo);
235 }
236 }
237
238 pb = RETVAL = gdk_pixbuf_new (GDK_COLORSPACE_RGB, 0, 8, cinfo.output_width, cinfo.output_height);
239 if (!RETVAL)
240 longjmp (jerr.setjmp_buffer, 2);
241
242 data = gdk_pixbuf_get_pixels (RETVAL);
243 rs = gdk_pixbuf_get_rowstride (RETVAL);
244
245 if (cinfo.output_components != 3)
246 longjmp (jerr.setjmp_buffer, 3);
247
248 jpeg_start_decompress (&cinfo);
249
250 while (cinfo.output_scanline < cinfo.output_height)
251 {
252 int remaining = cinfo.output_height - cinfo.output_scanline;
253 JSAMPROW rp[4];
254
255 rp [0] = data + cinfo.output_scanline * rs;
256 rp [1] = (guchar *)rp [0] + rs;
257 rp [2] = (guchar *)rp [1] + rs;
258 rp [3] = (guchar *)rp [2] + rs;
259
260 jpeg_read_scanlines (&cinfo, rp, remaining < 4 ? remaining : 4);
261 }
262
263 jpeg_finish_decompress (&cinfo);
264 fclose (fp);
265 jpeg_destroy_decompress (&cinfo);
266}
267 OUTPUT:
268 RETVAL
269
154############################################################################# 270#############################################################################
155 271
156MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer 272MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Schnauzer
273
274SV *
275foldcase (SV *pathsv)
276 PROTOTYPE: $
277 CODE:
278{
279 STRLEN plen;
280 U8 *path = SvPVutf8 (pathsv, plen);
281 U8 *pend = path + plen;
282 U8 dst [plen * 6 * 3], *dstp = dst;
283
284 while (path < pend)
285 {
286 U8 ch = *path;
287
288 if (ch >= 'a' && ch <= 'z')
289 *dstp++ = *path++;
290 else if (ch >= '0' && ch <= '9')
291 {
292 STRLEN el, nl = 0;
293 while (*path >= '0' && *path <= '9' && path < pend)
294 path++, nl++;
295
296 for (el = nl; el < 6; el++)
297 *dstp++ = '0';
298
299 memcpy (dstp, path - nl, nl);
300 dstp += nl;
301 }
302 else
303 {
304 STRLEN cl;
305 to_utf8_fold (path, dstp, &cl);
306 dstp += cl;
307 path += is_utf8_char (path);
308 }
309 }
310
311 RETVAL = newSVpvn (dst, dstp - dst);
312}
313 OUTPUT:
314 RETVAL
157 315
158GdkPixbuf_noinc * 316GdkPixbuf_noinc *
159p7_to_pb (int w, int h, guchar *src) 317p7_to_pb (int w, int h, guchar *src)
318 PROTOTYPE: @
160 CODE: 319 CODE:
161{ 320{
162 int x, y; 321 int x, y;
163 guchar *dst, *d; 322 guchar *dst, *d;
164 int dstr; 323 int dstr;
213 for (x = 0; x < w; x++) 372 for (x = 0; x < w; x++)
214 { 373 {
215 int r, g, b; 374 int r, g, b;
216 guchar *p = src + x * bpp + y * sstr; 375 guchar *p = src + x * bpp + y * sstr;
217 376
218 r = ((p[0] + er + Er[x]) * 7 + 128) / 255; 377 r = ((p[0] + er + Er[x]) * 7 + (RAND & 127) + 64) / 255;
219 g = ((p[1] + eg + Eg[x]) * 7 + 128) / 255; 378 g = ((p[1] + eg + Eg[x]) * 7 + (RAND & 127) + 64) / 255;
220 b = ((p[2] + eb + Eb[x]) * 3 + 128) / 255; 379 b = ((p[2] + eb + Eb[x]) * 3 + (RAND & 127) + 64) / 255;
221 380
222 r = r > 7 ? 7 : r < 0 ? 0 : r; 381 r = r > 7 ? 7 : r < 0 ? 0 : r;
223 g = g > 7 ? 7 : g < 0 ? 0 : g; 382 g = g > 7 ? 7 : g < 0 ? 0 : g;
224 b = b > 3 ? 3 : b < 0 ? 0 : b; 383 b = b > 3 ? 3 : b < 0 ? 0 : b;
225 384
226 er += p[0] - (r * 255 + 4) / 7; 385 er += p[0] - (r * 255 + 4) / 7;
227 eg += p[1] - (g * 255 + 4) / 7; 386 eg += p[1] - (g * 255 + 4) / 7;
228 eb += p[2] - (b * 255 + 2) / 3; 387 eb += p[2] - (b * 255 + 2) / 3;
229 388
230 Er[x] = er / 2; er -= er / 2 + RAND % 7 - 3; 389 Er[x] = er >> 1; er -= (er + 1) >> 1;
231 Eg[x] = eg / 2; eg -= eg / 2 + RAND % 7 - 3; 390 Eg[x] = eg >> 1; eg -= (eg + 1) >> 1;
232 Eb[x] = eb / 2; eb -= eb / 2 + RAND % 7 - 3; 391 Eb[x] = eb >> 1; eb -= (eb + 1) >> 1;
233 392
234 *dst++ = r << 5 | g << 2 | b; 393 *dst++ = r << 5 | g << 2 | b;
235 } 394 }
236 } 395 }
237} 396}
238 OUTPUT: 397 OUTPUT:
239 RETVAL 398 RETVAL
240
241SV *
242make_histogram (SV *ar)
243 CODE:
244{
245 int i;
246 AV *av, *result;
247
248 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
249 croak ("Not an array ref as first argument to make_histogram");
250
251 av = (AV *) SvRV (ar);
252 result = newAV ();
253
254 for (i = 0; i <= av_len (av); ++i)
255 {
256 int bigst, j;
257 SV *sv = *av_fetch (av, i, 1);
258 STRLEN len;
259 char *buf = SvPVbyte (sv, len);
260
261 int tmphist[256];
262 char *hist;
263
264 SV *histsv = newSV (257);
265 sv_upgrade (histsv, SVt_PV);
266 SvCUR_set (histsv, 256);
267 hist = SvPVX (histsv);
268
269 Zero (tmphist, 256, int);
270
271 while (len--)
272 ++tmphist[(unsigned char)*buf++];
273
274 bigst = 1;
275 for (j = 0; j < 256; ++j)
276 if (tmphist[j] > bigst)
277 bigst = tmphist[j];
278
279 for (j = 0; j < 256; ++j)
280 hist[j] = tmphist[j] * 255 / bigst;
281
282 av_push (result, histsv);
283 }
284
285 RETVAL = newRV_noinc ((SV *)result);
286}
287 OUTPUT:
288 RETVAL
289
290
291 399
292############################################################################# 400#############################################################################
293 401
294MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript 402MODULE = Gtk2::CV PACKAGE = Gtk2::CV::PostScript
295 403
419 } 527 }
420} 528}
421 OUTPUT: 529 OUTPUT:
422 RETVAL 530 RETVAL
423 531
532#############################################################################
424 533
534MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster
425 535
536SV *
537make_histograms (SV *ar)
538 CODE:
539{
540 int i;
541 AV *av, *result;
542
543 if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV)
544 croak ("Not an array ref as first argument to make_histogram");
545
546 av = (AV *) SvRV (ar);
547 result = newAV ();
548
549 for (i = 0; i <= av_len (av); ++i)
550 {
551 const int HISTSIZE = 64;
552
553 int j;
554 SV *sv = *av_fetch (av, i, 1);
555 STRLEN len;
556 char *buf = SvPVbyte (sv, len);
557
558 int tmphist[HISTSIZE];
559 float *hist;
560
561 SV *histsv = newSV (HISTSIZE * sizeof (float) + 1);
562 SvPOK_on (histsv);
563 SvCUR_set (histsv, HISTSIZE * sizeof (float));
564 hist = (float *)SvPVX (histsv);
565
566 Zero (tmphist, sizeof (tmphist), char);
567
568 for (j = len; j--; )
569 {
570 unsigned int idx
571 = ((*buf & 0xc0) >> 2)
572 | ((*buf & 0x18) >> 1)
573 | (*buf & 0x03);
574
575 ++tmphist[idx];
576 ++buf;
577 }
578
579 for (j = 0; j < HISTSIZE; ++j)
580 hist[j] = (float)tmphist[j] / (len + 1e-30);
581
582 av_push (result, histsv);
583 }
584
585 RETVAL = newRV_noinc ((SV *)result);
586}
587 OUTPUT:
588 RETVAL
589
590

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines