source: libpipi/trunk/pipi/resize.c @ 2902

Last change on this file since 2902 was 2902, checked in by Sam Hocevar, 11 years ago

Support C99 types on Win32 through the same hacks as in libcaca.

File size: 3.2 KB
Line 
1/*
2 *  libpipi       Pathetic image processing interface library
3 *  Copyright (c) 2004-2008 Sam Hocevar <sam@zoy.org>
4 *                All Rights Reserved
5 *
6 *  $Id$
7 *
8 *  This library is free software. It comes without any warranty, to
9 *  the extent permitted by applicable law. You can redistribute it
10 *  and/or modify it under the terms of the Do What The Fuck You Want
11 *  To Public License, Version 2, as published by Sam Hocevar. See
12 *  http://sam.zoy.org/wtfpl/COPYING for more details.
13 */
14
15/*
16 * resize.c: image resizing functions
17 */
18
19#include "config.h"
20
21#include <stdlib.h>
22#include <string.h>
23
24#include "pipi.h"
25#include "pipi_internals.h"
26
27pipi_image_t *pipi_resize(pipi_image_t *src, int w, int h)
28{
29    float *srcdata, *dstdata, *aline, *line;
30    pipi_image_t *dst;
31    pipi_pixels_t *srcp, *dstp;
32    int x, y, x0, y0, sw, dw, sh, dh, remy;
33
34    srcp = pipi_getpixels(src, PIPI_PIXELS_RGBA_F);
35    srcdata = (float *)srcp->pixels;
36
37    dst = pipi_new(w, h);
38    dstp = pipi_getpixels(dst, PIPI_PIXELS_RGBA_F);
39    dstdata = (float *)dstp->pixels;
40
41    sw = src->w; sh = src->h;
42    dw = dst->w; dh = dst->h;
43
44    aline = malloc(3 * dw * sizeof(float));
45    line = malloc(3 * dw * sizeof(float));
46
47    memset(line, 0, 3 * dw * sizeof(float));
48    remy = 0;
49
50    for(y = 0, y0 = 0; y < dh; y++)
51    {
52        int toty = 0, ny;
53
54        memset(aline, 0, 3 * dw * sizeof(float));
55
56        while(toty < sh)
57        {
58            if(remy == 0)
59            {
60                float r = 0, g = 0, b = 0;
61                int remx = 0;
62
63                for(x = 0, x0 = 0; x < dw; x++)
64                {
65                    float ar = 0, ag = 0, ab = 0;
66                    int totx = 0, nx;
67
68                    while(totx < sw)
69                    {
70                        if(remx == 0)
71                        {
72                            r = srcdata[(y0 * sw + x0) * 4];
73                            g = srcdata[(y0 * sw + x0) * 4 + 1];
74                            b = srcdata[(y0 * sw + x0) * 4 + 2];
75                            x0++;
76                            remx = dw;
77                        }
78
79                        nx = (totx + remx <= sw) ? remx : sw - totx;
80                        ar += nx * r; ag += nx * g; ab += nx * b;
81                        totx += nx;
82                        remx -= nx;
83                    }
84
85                    line[3 * x] = ar;
86                    line[3 * x + 1] = ag;
87                    line[3 * x + 2] = ab;
88                }
89
90                y0++;
91                remy = dh;
92            }
93
94            ny = (toty + remy <= sh) ? remy : sh - toty;
95            for(x = 0; x < dw; x++)
96            {
97                aline[3 * x] += ny * line[3 * x];
98                aline[3 * x + 1] += ny * line[3 * x + 1];
99                aline[3 * x + 2] += ny * line[3 * x + 2];
100            }
101            toty += ny;
102            remy -= ny;
103        }
104
105        for(x = 0; x < dw; x++)
106        {
107            dstdata[(y * dw + x) * 4] = aline[3 * x] / (sw * sh);
108            dstdata[(y * dw + x) * 4 + 1] = aline[3 * x + 1] / (sw * sh);
109            dstdata[(y * dw + x) * 4 + 2] = aline[3 * x + 2] / (sw * sh);
110        }
111    }
112
113    free(aline);
114    free(line);
115
116    return dst;
117}
118
Note: See TracBrowser for help on using the repository browser.