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

Last change on this file since 3342 was 3342, checked in by Sam Hocevar, 12 years ago

Change _C pixel format suffixes into _U8 for more clarity.

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_get_pixels(src, PIPI_PIXELS_RGBA_F32);
35    srcdata = (float *)srcp->pixels;
36
37    dst = pipi_new(w, h);
38    dstp = pipi_get_pixels(dst, PIPI_PIXELS_RGBA_F32);
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.