2 24c02865 2005-01-04 devnull * rotate an image 180° in O(log Dx + log Dy) /dev/draw writes,
3 24c02865 2005-01-04 devnull * using an extra buffer same size as the image.
5 24c02865 2005-01-04 devnull * the basic concept is that you can invert an array by inverting
6 24c02865 2005-01-04 devnull * the top half, inverting the bottom half, and then swapping them.
7 24c02865 2005-01-04 devnull * the code does this slightly backwards to ensure O(log n) runtime.
8 24c02865 2005-01-04 devnull * (If you do it wrong, you can get O(log² n) runtime.)
10 24c02865 2005-01-04 devnull * This is usually overkill, but it speeds up slow remote
11 24c02865 2005-01-04 devnull * connections quite a bit.
14 24c02865 2005-01-04 devnull #include <u.h>
15 24c02865 2005-01-04 devnull #include <libc.h>
16 24c02865 2005-01-04 devnull #include <bio.h>
17 24c02865 2005-01-04 devnull #include <draw.h>
18 24c02865 2005-01-04 devnull #include <event.h>
19 24c02865 2005-01-04 devnull #include "page.h"
21 24c02865 2005-01-04 devnull int ndraw = 0;
23 24c02865 2005-01-04 devnull Xaxis = 0,
24 24c02865 2005-01-04 devnull Yaxis = 1,
27 24c02865 2005-01-04 devnull Image *mtmp;
30 24c02865 2005-01-04 devnull writefile(char *name, Image *im, int gran)
32 24c02865 2005-01-04 devnull static int c = 100;
34 24c02865 2005-01-04 devnull char buf[200];
36 24c02865 2005-01-04 devnull snprint(buf, sizeof buf, "%d%s%d", c++, name, gran);
37 24c02865 2005-01-04 devnull fd = create(buf, OWRITE, 0666);
38 24c02865 2005-01-04 devnull if(fd < 0)
40 24c02865 2005-01-04 devnull writeimage(fd, im, 0);
41 24c02865 2005-01-04 devnull close(fd);
45 24c02865 2005-01-04 devnull moveup(Image *im, Image *tmp, int a, int b, int c, int axis)
47 24c02865 2005-01-04 devnull Rectangle range;
48 24c02865 2005-01-04 devnull Rectangle dr0, dr1;
49 24c02865 2005-01-04 devnull Point p0, p1;
51 24c02865 2005-01-04 devnull if(a == b || b == c)
54 24c02865 2005-01-04 devnull drawop(tmp, tmp->r, im, nil, im->r.min, S);
56 24c02865 2005-01-04 devnull switch(axis){
57 24c02865 2005-01-04 devnull case Xaxis:
58 24c02865 2005-01-04 devnull range = Rect(a, im->r.min.y, c, im->r.max.y);
59 24c02865 2005-01-04 devnull dr0 = range;
60 24c02865 2005-01-04 devnull dr0.max.x = dr0.min.x+(c-b);
61 24c02865 2005-01-04 devnull p0 = Pt(b, im->r.min.y);
63 24c02865 2005-01-04 devnull dr1 = range;
64 24c02865 2005-01-04 devnull dr1.min.x = dr1.max.x-(b-a);
65 24c02865 2005-01-04 devnull p1 = Pt(a, im->r.min.y);
67 24c02865 2005-01-04 devnull case Yaxis:
68 24c02865 2005-01-04 devnull range = Rect(im->r.min.x, a, im->r.max.x, c);
69 24c02865 2005-01-04 devnull dr0 = range;
70 24c02865 2005-01-04 devnull dr0.max.y = dr0.min.y+(c-b);
71 24c02865 2005-01-04 devnull p0 = Pt(im->r.min.x, b);
73 24c02865 2005-01-04 devnull dr1 = range;
74 24c02865 2005-01-04 devnull dr1.min.y = dr1.max.y-(b-a);
75 24c02865 2005-01-04 devnull p1 = Pt(im->r.min.x, a);
78 24c02865 2005-01-04 devnull drawop(im, dr0, tmp, nil, p0, S);
79 24c02865 2005-01-04 devnull drawop(im, dr1, tmp, nil, p1, S);
83 24c02865 2005-01-04 devnull interlace(Image *im, Image *tmp, int axis, int n, Image *mask, int gran)
85 24c02865 2005-01-04 devnull Point p0, p1;
86 24c02865 2005-01-04 devnull Rectangle r0, r1;
88 24c02865 2005-01-04 devnull r0 = im->r;
89 24c02865 2005-01-04 devnull r1 = im->r;
90 24c02865 2005-01-04 devnull switch(axis) {
91 24c02865 2005-01-04 devnull case Xaxis:
92 24c02865 2005-01-04 devnull r0.max.x = n;
93 24c02865 2005-01-04 devnull r1.min.x = n;
94 24c02865 2005-01-04 devnull p0 = (Point){gran, 0};
95 24c02865 2005-01-04 devnull p1 = (Point){-gran, 0};
97 24c02865 2005-01-04 devnull case Yaxis:
98 24c02865 2005-01-04 devnull r0.max.y = n;
99 24c02865 2005-01-04 devnull r1.min.y = n;
100 24c02865 2005-01-04 devnull p0 = (Point){0, gran};
101 24c02865 2005-01-04 devnull p1 = (Point){0, -gran};
105 24c02865 2005-01-04 devnull drawop(tmp, im->r, im, display->opaque, im->r.min, S);
106 24c02865 2005-01-04 devnull gendrawop(im, r0, tmp, p0, mask, mask->r.min, S);
107 24c02865 2005-01-04 devnull gendrawop(im, r0, tmp, p1, mask, p1, S);
111 24c02865 2005-01-04 devnull * Halve the grating period in the mask.
112 24c02865 2005-01-04 devnull * The grating currently looks like
113 24c02865 2005-01-04 devnull * ####____####____####____####____
114 24c02865 2005-01-04 devnull * where #### is opacity.
116 24c02865 2005-01-04 devnull * We want
117 24c02865 2005-01-04 devnull * ##__##__##__##__##__##__##__##__
118 24c02865 2005-01-04 devnull * which is achieved by shifting the mask
119 24c02865 2005-01-04 devnull * and drawing on itself through itself.
120 24c02865 2005-01-04 devnull * Draw doesn't actually allow this, so
121 24c02865 2005-01-04 devnull * we have to copy it first.
123 24c02865 2005-01-04 devnull * ####____####____####____####____ (dst)
124 24c02865 2005-01-04 devnull * + ____####____####____####____#### (src)
125 24c02865 2005-01-04 devnull * in __####____####____####____####__ (mask)
126 24c02865 2005-01-04 devnull * ===========================================
127 24c02865 2005-01-04 devnull * ##__##__##__##__##__##__##__##__
130 24c02865 2005-01-04 devnull nextmask(Image *mask, int axis, int maskdim)
132 24c02865 2005-01-04 devnull Point delta;
134 24c02865 2005-01-04 devnull delta = axis==Xaxis ? Pt(maskdim,0) : Pt(0,maskdim);
135 24c02865 2005-01-04 devnull drawop(mtmp, mtmp->r, mask, nil, mask->r.min, S);
136 24c02865 2005-01-04 devnull gendrawop(mask, mask->r, mtmp, delta, mtmp, divpt(delta,-2), S);
137 24c02865 2005-01-04 devnull // writefile("mask", mask, maskdim/2);
138 24c02865 2005-01-04 devnull return maskdim/2;
142 24c02865 2005-01-04 devnull shuffle(Image *im, Image *tmp, int axis, int n, Image *mask, int gran,
143 24c02865 2005-01-04 devnull int lastnn)
145 24c02865 2005-01-04 devnull int nn, left;
147 24c02865 2005-01-04 devnull if(gran == 0)
149 24c02865 2005-01-04 devnull left = n%(2*gran);
150 24c02865 2005-01-04 devnull nn = n - left;
152 24c02865 2005-01-04 devnull interlace(im, tmp, axis, nn, mask, gran);
153 24c02865 2005-01-04 devnull // writefile("interlace", im, gran);
155 24c02865 2005-01-04 devnull gran = nextmask(mask, axis, gran);
156 24c02865 2005-01-04 devnull shuffle(im, tmp, axis, n, mask, gran, nn);
157 24c02865 2005-01-04 devnull // writefile("shuffle", im, gran);
158 24c02865 2005-01-04 devnull moveup(im, tmp, lastnn, nn, n, axis);
159 24c02865 2005-01-04 devnull // writefile("move", im, gran);
163 24c02865 2005-01-04 devnull rot180(Image *im)
165 24c02865 2005-01-04 devnull Image *tmp, *tmp0;
166 24c02865 2005-01-04 devnull Image *mask;
167 24c02865 2005-01-04 devnull Rectangle rmask;
168 24c02865 2005-01-04 devnull int gran;
170 24c02865 2005-01-04 devnull if(chantodepth(im->chan) < 8){
171 24c02865 2005-01-04 devnull /* this speeds things up dramatically; draw is too slow on sub-byte pixel sizes */
172 24c02865 2005-01-04 devnull tmp0 = xallocimage(display, im->r, CMAP8, 0, DNofill);
173 24c02865 2005-01-04 devnull drawop(tmp0, tmp0->r, im, nil, im->r.min, S);
175 24c02865 2005-01-04 devnull tmp0 = im;
177 24c02865 2005-01-04 devnull tmp = xallocimage(display, tmp0->r, tmp0->chan, 0, DNofill);
178 24c02865 2005-01-04 devnull if(tmp == nil){
179 24c02865 2005-01-04 devnull if(tmp0 != im)
180 24c02865 2005-01-04 devnull freeimage(tmp0);
183 24c02865 2005-01-04 devnull for(gran=1; gran<Dx(im->r); gran *= 2)
185 24c02865 2005-01-04 devnull gran /= 4;
187 24c02865 2005-01-04 devnull rmask.min = ZP;
188 24c02865 2005-01-04 devnull rmask.max = (Point){2*gran, 100};
190 24c02865 2005-01-04 devnull mask = xallocimage(display, rmask, GREY1, 1, DTransparent);
191 24c02865 2005-01-04 devnull mtmp = xallocimage(display, rmask, GREY1, 1, DTransparent);
192 24c02865 2005-01-04 devnull if(mask == nil || mtmp == nil) {
193 24c02865 2005-01-04 devnull fprint(2, "out of memory during rot180: %r\n");
194 24c02865 2005-01-04 devnull wexits("memory");
196 24c02865 2005-01-04 devnull rmask.max.x = gran;
197 24c02865 2005-01-04 devnull drawop(mask, rmask, display->opaque, nil, ZP, S);
198 24c02865 2005-01-04 devnull // writefile("mask", mask, gran);
199 24c02865 2005-01-04 devnull shuffle(im, tmp, Xaxis, Dx(im->r), mask, gran, 0);
200 24c02865 2005-01-04 devnull freeimage(mask);
201 24c02865 2005-01-04 devnull freeimage(mtmp);
203 24c02865 2005-01-04 devnull for(gran=1; gran<Dy(im->r); gran *= 2)
205 24c02865 2005-01-04 devnull gran /= 4;
206 24c02865 2005-01-04 devnull rmask.max = (Point){100, 2*gran};
207 24c02865 2005-01-04 devnull mask = xallocimage(display, rmask, GREY1, 1, DTransparent);
208 24c02865 2005-01-04 devnull mtmp = xallocimage(display, rmask, GREY1, 1, DTransparent);
209 24c02865 2005-01-04 devnull if(mask == nil || mtmp == nil) {
210 24c02865 2005-01-04 devnull fprint(2, "out of memory during rot180: %r\n");
211 24c02865 2005-01-04 devnull wexits("memory");
213 24c02865 2005-01-04 devnull rmask.max.y = gran;
214 24c02865 2005-01-04 devnull drawop(mask, rmask, display->opaque, nil, ZP, S);
215 24c02865 2005-01-04 devnull shuffle(im, tmp, Yaxis, Dy(im->r), mask, gran, 0);
216 24c02865 2005-01-04 devnull freeimage(mask);
217 24c02865 2005-01-04 devnull freeimage(mtmp);
218 24c02865 2005-01-04 devnull freeimage(tmp);
219 24c02865 2005-01-04 devnull if(tmp0 != im)
220 24c02865 2005-01-04 devnull freeimage(tmp0);
223 24c02865 2005-01-04 devnull /* rotates an image 90 degrees clockwise */
225 24c02865 2005-01-04 devnull rot90(Image *im)
227 24c02865 2005-01-04 devnull Image *tmp;
228 24c02865 2005-01-04 devnull int i, j, dx, dy;
230 24c02865 2005-01-04 devnull dx = Dx(im->r);
231 24c02865 2005-01-04 devnull dy = Dy(im->r);
232 24c02865 2005-01-04 devnull tmp = xallocimage(display, Rect(0, 0, dy, dx), im->chan, 0, DCyan);
233 24c02865 2005-01-04 devnull if(tmp == nil) {
234 24c02865 2005-01-04 devnull fprint(2, "out of memory during rot90: %r\n");
235 24c02865 2005-01-04 devnull wexits("memory");
238 24c02865 2005-01-04 devnull for(j = 0; j < dx; j++) {
239 24c02865 2005-01-04 devnull for(i = 0; i < dy; i++) {
240 24c02865 2005-01-04 devnull drawop(tmp, Rect(i, j, i+1, j+1), im, nil, Pt(j, dy-(i+1)), S);
243 24c02865 2005-01-04 devnull freeimage(im);
245 24c02865 2005-01-04 devnull return(tmp);
248 24c02865 2005-01-04 devnull /* from resample.c -- resize from → to using interpolation */
251 24c02865 2005-01-04 devnull #define K2 7 /* from -.7 to +.7 inclusive, meaning .2 into each adjacent pixel */
252 24c02865 2005-01-04 devnull #define NK (2*K2+1)
253 24c02865 2005-01-04 devnull double K[NK];
256 24c02865 2005-01-04 devnull fac(int L)
258 24c02865 2005-01-04 devnull int i, f;
261 24c02865 2005-01-04 devnull for(i=L; i>1; --i)
263 24c02865 2005-01-04 devnull return f;
267 24c02865 2005-01-04 devnull * i0(x) is the modified Bessel function, Σ (x/2)^2L / (L!)²
268 24c02865 2005-01-04 devnull * There are faster ways to calculate this, but we precompute
269 24c02865 2005-01-04 devnull * into a table so let's keep it simple.
272 24c02865 2005-01-04 devnull i0(double x)
274 24c02865 2005-01-04 devnull double v;
277 24c02865 2005-01-04 devnull v = 1.0;
278 24c02865 2005-01-04 devnull for(L=1; L<10; L++)
279 24c02865 2005-01-04 devnull v += pow(x/2., 2*L)/pow(fac(L), 2);
280 24c02865 2005-01-04 devnull return v;
284 24c02865 2005-01-04 devnull kaiser(double x, double tau, double alpha)
286 24c02865 2005-01-04 devnull if(fabs(x) > tau)
287 24c02865 2005-01-04 devnull return 0.;
288 24c02865 2005-01-04 devnull return i0(alpha*sqrt(1-(x*x/(tau*tau))))/i0(alpha);
292 24c02865 2005-01-04 devnull resamplex(uchar *in, int off, int d, int inx, uchar *out, int outx)
294 24c02865 2005-01-04 devnull int i, x, k;
295 24c02865 2005-01-04 devnull double X, xx, v, rat;
298 24c02865 2005-01-04 devnull rat = (double)inx/(double)outx;
299 24c02865 2005-01-04 devnull for(x=0; x<outx; x++){
300 24c02865 2005-01-04 devnull if(inx == outx){
301 24c02865 2005-01-04 devnull /* don't resample if size unchanged */
302 24c02865 2005-01-04 devnull out[off+x*d] = in[off+x*d];
303 24c02865 2005-01-04 devnull continue;
305 24c02865 2005-01-04 devnull v = 0.0;
306 24c02865 2005-01-04 devnull X = x*rat;
307 24c02865 2005-01-04 devnull for(k=-K2; k<=K2; k++){
308 24c02865 2005-01-04 devnull xx = X + rat*k/10.;
310 24c02865 2005-01-04 devnull if(i < 0)
312 24c02865 2005-01-04 devnull if(i >= inx)
313 24c02865 2005-01-04 devnull i = inx-1;
314 24c02865 2005-01-04 devnull v += in[off+i*d] * K[K2+k];
316 24c02865 2005-01-04 devnull out[off+x*d] = v;
321 24c02865 2005-01-04 devnull resampley(uchar **in, int off, int iny, uchar **out, int outy)
323 24c02865 2005-01-04 devnull int y, i, k;
324 24c02865 2005-01-04 devnull double Y, yy, v, rat;
326 24c02865 2005-01-04 devnull rat = (double)iny/(double)outy;
327 24c02865 2005-01-04 devnull for(y=0; y<outy; y++){
328 24c02865 2005-01-04 devnull if(iny == outy){
329 24c02865 2005-01-04 devnull /* don't resample if size unchanged */
330 24c02865 2005-01-04 devnull out[y][off] = in[y][off];
331 24c02865 2005-01-04 devnull continue;
333 24c02865 2005-01-04 devnull v = 0.0;
334 24c02865 2005-01-04 devnull Y = y*rat;
335 24c02865 2005-01-04 devnull for(k=-K2; k<=K2; k++){
336 24c02865 2005-01-04 devnull yy = Y + rat*k/10.;
338 24c02865 2005-01-04 devnull if(i < 0)
340 24c02865 2005-01-04 devnull if(i >= iny)
341 24c02865 2005-01-04 devnull i = iny-1;
342 24c02865 2005-01-04 devnull v += in[i][off] * K[K2+k];
344 24c02865 2005-01-04 devnull out[y][off] = v;
350 24c02865 2005-01-04 devnull resample(Image *from, Image *to)
352 24c02865 2005-01-04 devnull int i, j, bpl, nchan;
353 24c02865 2005-01-04 devnull uchar **oscan, **nscan;
354 24c02865 2005-01-04 devnull char tmp[20];
355 24c02865 2005-01-04 devnull int xsize, ysize;
356 24c02865 2005-01-04 devnull double v;
357 24c02865 2005-01-04 devnull Image *t1, *t2;
358 24c02865 2005-01-04 devnull ulong tchan;
360 24c02865 2005-01-04 devnull for(i=-K2; i<=K2; i++){
361 24c02865 2005-01-04 devnull K[K2+i] = kaiser(i/10., K2/10., 4.);
364 24c02865 2005-01-04 devnull /* normalize */
365 24c02865 2005-01-04 devnull v = 0.0;
366 24c02865 2005-01-04 devnull for(i=0; i<NK; i++)
367 24c02865 2005-01-04 devnull v += K[i];
368 24c02865 2005-01-04 devnull for(i=0; i<NK; i++)
369 24c02865 2005-01-04 devnull K[i] /= v;
371 24c02865 2005-01-04 devnull switch(from->chan){
372 24c02865 2005-01-04 devnull case GREY8:
373 24c02865 2005-01-04 devnull case RGB24:
374 24c02865 2005-01-04 devnull case RGBA32:
375 24c02865 2005-01-04 devnull case ARGB32:
376 24c02865 2005-01-04 devnull case XRGB32:
379 24c02865 2005-01-04 devnull case CMAP8:
380 24c02865 2005-01-04 devnull case RGB15:
381 24c02865 2005-01-04 devnull case RGB16:
382 24c02865 2005-01-04 devnull tchan = RGB24;
383 24c02865 2005-01-04 devnull goto Convert;
385 24c02865 2005-01-04 devnull case GREY1:
386 24c02865 2005-01-04 devnull case GREY2:
387 24c02865 2005-01-04 devnull case GREY4:
388 24c02865 2005-01-04 devnull tchan = GREY8;
389 24c02865 2005-01-04 devnull Convert:
390 24c02865 2005-01-04 devnull /* use library to convert to byte-per-chan form, then convert back */
391 24c02865 2005-01-04 devnull t1 = xallocimage(display, Rect(0, 0, Dx(from->r), Dy(from->r)), tchan, 0, DNofill);
392 24c02865 2005-01-04 devnull if(t1 == nil) {
393 24c02865 2005-01-04 devnull fprint(2, "out of memory for temp image 1 in resample: %r\n");
394 24c02865 2005-01-04 devnull wexits("memory");
396 24c02865 2005-01-04 devnull drawop(t1, t1->r, from, nil, ZP, S);
397 24c02865 2005-01-04 devnull t2 = xallocimage(display, to->r, tchan, 0, DNofill);
398 24c02865 2005-01-04 devnull if(t2 == nil) {
399 24c02865 2005-01-04 devnull fprint(2, "out of memory temp image 2 in resample: %r\n");
400 24c02865 2005-01-04 devnull wexits("memory");
402 24c02865 2005-01-04 devnull resample(t1, t2);
403 24c02865 2005-01-04 devnull drawop(to, to->r, t2, nil, ZP, S);
404 24c02865 2005-01-04 devnull freeimage(t1);
405 24c02865 2005-01-04 devnull freeimage(t2);
406 24c02865 2005-01-04 devnull return to;
408 24c02865 2005-01-04 devnull default:
409 24c02865 2005-01-04 devnull sysfatal("can't handle channel type %s", chantostr(tmp, from->chan));
412 24c02865 2005-01-04 devnull xsize = Dx(to->r);
413 24c02865 2005-01-04 devnull ysize = Dy(to->r);
414 24c02865 2005-01-04 devnull oscan = malloc(Dy(from->r)*sizeof(uchar*));
415 24c02865 2005-01-04 devnull nscan = malloc(max(ysize, Dy(from->r))*sizeof(uchar*));
416 24c02865 2005-01-04 devnull if(oscan == nil || nscan == nil)
417 24c02865 2005-01-04 devnull sysfatal("can't allocate: %r");
419 24c02865 2005-01-04 devnull /* unload original image into scan lines */
420 24c02865 2005-01-04 devnull bpl = bytesperline(from->r, from->depth);
421 24c02865 2005-01-04 devnull for(i=0; i<Dy(from->r); i++){
422 24c02865 2005-01-04 devnull oscan[i] = malloc(bpl);
423 24c02865 2005-01-04 devnull if(oscan[i] == nil)
424 24c02865 2005-01-04 devnull sysfatal("can't allocate: %r");
425 24c02865 2005-01-04 devnull j = unloadimage(from, Rect(from->r.min.x, from->r.min.y+i, from->r.max.x, from->r.min.y+i+1), oscan[i], bpl);
426 24c02865 2005-01-04 devnull if(j != bpl)
427 24c02865 2005-01-04 devnull sysfatal("unloadimage");
430 24c02865 2005-01-04 devnull /* allocate scan lines for destination. we do y first, so need at least Dy(from->r) lines */
431 24c02865 2005-01-04 devnull bpl = bytesperline(Rect(0, 0, xsize, Dy(from->r)), from->depth);
432 24c02865 2005-01-04 devnull for(i=0; i<max(ysize, Dy(from->r)); i++){
433 24c02865 2005-01-04 devnull nscan[i] = malloc(bpl);
434 24c02865 2005-01-04 devnull if(nscan[i] == nil)
435 24c02865 2005-01-04 devnull sysfatal("can't allocate: %r");
438 24c02865 2005-01-04 devnull /* resample in X */
439 24c02865 2005-01-04 devnull nchan = from->depth/8;
440 24c02865 2005-01-04 devnull for(i=0; i<Dy(from->r); i++){
441 24c02865 2005-01-04 devnull for(j=0; j<nchan; j++){
442 24c02865 2005-01-04 devnull if(j==0 && from->chan==XRGB32)
443 24c02865 2005-01-04 devnull continue;
444 24c02865 2005-01-04 devnull resamplex(oscan[i], j, nchan, Dx(from->r), nscan[i], xsize);
446 24c02865 2005-01-04 devnull free(oscan[i]);
447 24c02865 2005-01-04 devnull oscan[i] = nscan[i];
448 24c02865 2005-01-04 devnull nscan[i] = malloc(bpl);
449 24c02865 2005-01-04 devnull if(nscan[i] == nil)
450 24c02865 2005-01-04 devnull sysfatal("can't allocate: %r");
453 24c02865 2005-01-04 devnull /* resample in Y */
454 24c02865 2005-01-04 devnull for(i=0; i<xsize; i++)
455 24c02865 2005-01-04 devnull for(j=0; j<nchan; j++)
456 24c02865 2005-01-04 devnull resampley(oscan, nchan*i+j, Dy(from->r), nscan, ysize);
458 24c02865 2005-01-04 devnull /* pack data into destination */
459 24c02865 2005-01-04 devnull bpl = bytesperline(to->r, from->depth);
460 24c02865 2005-01-04 devnull for(i=0; i<ysize; i++){
461 24c02865 2005-01-04 devnull j = loadimage(to, Rect(0, i, xsize, i+1), nscan[i], bpl);
462 24c02865 2005-01-04 devnull if(j != bpl)
463 24c02865 2005-01-04 devnull sysfatal("loadimage: %r");
466 24c02865 2005-01-04 devnull for(i=0; i<Dy(from->r); i++){
467 24c02865 2005-01-04 devnull free(oscan[i]);
468 24c02865 2005-01-04 devnull free(nscan[i]);
470 24c02865 2005-01-04 devnull free(oscan);
471 24c02865 2005-01-04 devnull free(nscan);
473 24c02865 2005-01-04 devnull return to;
