source: branches/kmk/Lucide/plugins/lujpeg/jpscale.cpp @ 306

Last change on this file since 306 was 306, checked in by dmik, 11 years ago

branches/kmk: Made all Lucide SOM plugins build with kmk/gcc.

File size: 11.7 KB
Line 
1/* ***** BEGIN LICENSE BLOCK *****
2 * Version: CDDL 1.0/LGPL 2.1
3 *
4 * The contents of this file are subject to the COMMON DEVELOPMENT AND
5 * DISTRIBUTION LICENSE (CDDL) Version 1.0 (the "License"); you may not use
6 * this file except in compliance with the License. You may obtain a copy of
7 * the License at http://www.sun.com/cddl/
8 *
9 * Software distributed under the License is distributed on an "AS IS" basis,
10 * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
11 * for the specific language governing rights and limitations under the
12 * License.
13 *
14 * The Initial Developer of the Original Code is
15 * Eugene Romanenko, netlabs.org.
16 * Portions created by the Initial Developer are Copyright (C) 2006
17 * the Initial Developer. All Rights Reserved.
18 *
19 * Contributor(s):
20 *
21 * Alternatively, the contents of this file may be used under the terms of
22 * the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
23 * in which case the provisions of the LGPL are applicable instead of those
24 * above. If you wish to allow use of your version of this file only under the
25 * terms of the LGPL, and not to allow others to use your version of this file
26 * under the terms of the CDDL, indicate your decision by deleting the
27 * provisions above and replace them with the notice and other provisions
28 * required by the LGPL. If you do not delete the provisions above, a recipient
29 * may use your version of this file under the terms of any one of the CDDL
30 * or the LGPL.
31 *
32 * ***** END LICENSE BLOCK ***** */
33
34// This file based on public domain code.
35// Original code by Dale Schumacher, public domain 1991
36// See _Graphics Gems III_ "General Filtered Image Rescaling", Dale A. Schumacher
37
38
39#include <stdint.h>
40#include <assert.h>
41#include <string.h>
42#include <math.h>
43#include <vector>
44#include <algorithm>
45
46#include <LuPixbuf.xh>
47
48
49#define ceilf ceil
50#define floorf floor
51
52inline int iclamp( int i, int min, int max ) {
53    assert( min <= max );
54    return std::max( min, std::min(i, max) );
55}
56
57inline int frnd( float f ) { return (int)( f + 0.5f ); }
58
59inline void *my_calloc( int count, int size )
60{
61    void *mem = (void *)new char[ count * size ];
62    memset( mem, 0, count * size );
63    return mem;
64}
65
66inline void my_cfree( void* mem )
67{
68    delete [] (char*)mem;
69}
70
71
72// Copy RGB data from the specified row into the given buffer.
73void get_row( Environment *ev, uint8_t* row, LuPixbuf *image, int x0, int xsize, int y )
74{
75    long m_width    = image->getWidth( ev );
76    long m_height   = image->getHeight( ev );
77    long m_pitch    = image->getRowSize( ev );
78    uint8_t *m_data = (uint8_t*)image->getDataPtr( ev );
79
80    y = iclamp( y, 0, m_height - 1 );
81    int x1 = x0 + xsize - 1;
82    if ( x1 >= m_width )
83    {
84        // clip, then extend.
85        int extra_pixels = x1 - m_width + 1;
86        uint8_t*  p = m_data + (y * m_pitch);
87        memcpy( row, p + x0 * 3, (3 * (m_width - x0)) );
88        // repeat last pixel
89        p = p + (m_width - 1) * 3;
90        uint8_t* q = row + (m_width - x0) * 3;
91        while (extra_pixels > 0)
92        {
93            *(q + 0) = *(p + 0);
94            *(q + 1) = *(p + 1);
95            *(q + 2) = *(p + 2);
96            q += 3;
97            extra_pixels--;
98        }
99    }
100    else {
101        memcpy( row, m_data + (y * m_pitch) + x0 * 3, (3 * xsize) );
102    }
103}
104
105
106
107// Copy RGB data from the specified column into the given buffer.
108void get_column( Environment *ev, uint8_t* column, LuPixbuf *image, int x )
109{
110    long m_width    = image->getWidth( ev );
111    long m_height   = image->getHeight( ev );
112    long m_pitch    = image->getRowSize( ev );
113    uint8_t *m_data = (uint8_t*)image->getDataPtr( ev );
114
115    if ((x < 0) || (x >= m_width)) {
116        assert(0);
117        x = iclamp(x, 0, m_width - 1);
118    }
119
120    int d = m_pitch;
121    uint8_t* p = m_data + x * 3;
122    for (int i = m_height; i-- > 0; p += d) {
123        *column++ = *p;
124        *column++ = *(p + 1);
125        *column++ = *(p + 2);
126    }
127}
128
129
130
131// Clamp {r, g, b} to [0,255], and write pixel data to the given image at (x, y).
132void put_pixel( Environment *ev, LuPixbuf *image, int x, int y, float r, float g, float b )
133{
134    long m_width    = image->getWidth( ev );
135    long m_height   = image->getHeight( ev );
136    long m_pitch    = image->getRowSize( ev );
137    uint8_t *m_data = (uint8_t*)image->getDataPtr( ev );
138
139    if ((x < 0) || (x >= m_width) || (y < 0) || (y >= m_height)) {
140        assert(0);
141        return;
142    }
143
144    uint8_t *p = m_data + (y * m_pitch);
145
146    p[x * 3 + 0] = iclamp(frnd(r), 0, 255);
147    p[x * 3 + 1] = iclamp(frnd(g), 0, 255);
148    p[x * 3 + 2] = iclamp(frnd(b), 0, 255);
149}
150
151
152
153/*
154 *   filter function definitions
155 */
156
157// MITCHELL
158
159#define SUPPORT (2.0f)
160
161#define B       (1.0f / 3.0f)
162#define C       (1.0f / 3.0f)
163
164float Mitchell_filter(float t)
165{
166    float tt = t * t;
167    if (t < 0.0f) t = -t;
168    if (t < 1.0f) {
169        t = (((12.0f - 9.0f * B - 6.0f * C) * (t * tt))
170             + ((-18.0f + 12.0f * B + 6.0f * C) * tt)
171             + (6.0f - 2.0f * B));
172        return(t / 6.0f);
173    } else if(t < 2.0f) {
174        t = (((-1.0f * B - 6.0f * C) * (t * tt))
175             + ((6.0f * B + 30.0f * C) * tt)
176             + ((-12.0f * B - 48.0f * C) * t)
177             + (8.0f * B + 24 * C));
178        return(t / 6.0f);
179    }
180    return(0.0f);
181}
182
183
184struct CONTRIB
185{
186    int   pixel;
187    float weight;
188
189    CONTRIB() : pixel(0), weight(0.f) {}
190
191    CONTRIB(int p, float w) : pixel(p), weight(w) {}
192};
193
194
195// Rescale the specified portion of the input image into the specified
196// portion of the output image.  Coordinates are *inclusive*.
197void resample( Environment *ev,
198               LuPixbuf *out, int out_x0, int out_y0, int out_x1, int out_y1,
199               LuPixbuf *in, float in_x0, float in_y0, float in_x1, float in_y1 )
200{
201    long out_m_width  = out->getWidth( ev );
202    long out_m_height = out->getHeight( ev );
203
204    assert(out_x0 <= out_x1);
205    assert(out_y0 <= out_y1);
206    assert(out_x0 >= 0 && out_x0 < out_m_width);
207    assert(out_x1 >= 0 && out_x1 < out_m_width);
208    assert(out_y0 >= 0 && out_y0 < out_m_height);
209    assert(out_y1 >= 0 && out_y1 < out_m_height);
210
211
212    int i, k;                   /* loop variables */
213    unsigned int j;                     /* loop variables */
214    int n;                              /* pixel number */
215    float center; int left, right;      /* filter calculation variables */
216    float width, fscale, weight;        /* filter calculation variables */
217    uint8_t*    raster;                 /* a row or column of pixels */
218
219    std::vector< std::vector<CONTRIB> > contrib;
220
221    int out_width = out_x1 - out_x0 + 1;
222    int out_height = out_y1 - out_y0 + 1;
223    assert(out_width > 0);
224    assert(out_height > 0);
225
226    float       in_width = in_x1 - in_x0;
227    float       in_height = in_y1 - in_y0;
228    assert(in_width > 0);
229    assert(in_height > 0);
230
231    int in_window_w = int(ceilf(in_x1) - floorf(in_x0) + 1);
232    int in_window_h = int(ceilf(in_y1) - floorf(in_y0) + 1);
233
234    /* create intermediate image to hold horizontal zoom */
235    LuPixbuf *tmp = new LuPixbuf( ev, out_width, in_window_h, 3 );
236    long tmp_m_width   = tmp->getWidth( ev );
237    long tmp_m_height  = tmp->getHeight( ev );
238
239    float xscale = (float) (out_width - 1) / in_width;
240    float yscale = (float) (out_height - 1) / in_height;
241
242    // xxxx protect against division by 0
243    if (yscale == 0) { yscale = 1.0f; }
244    if (xscale == 0) { xscale = 1.0f; }
245
246    /* pre-calculate filter contributions for a row */
247    contrib.resize(tmp_m_width);
248    if(xscale < 1.0f)
249    {
250        width = SUPPORT / xscale;
251        fscale = 1.0f / xscale;
252        for (i = 0; i < tmp_m_width; ++i)
253        {
254            contrib[i].resize(0);
255
256            center = (float) i / xscale;
257            left = int(ceilf(center - width));
258            right = int(floorf(center + width));
259            for (k = left; k <= right; ++k)
260            {
261                weight = center - (float) k;
262                weight = Mitchell_filter(weight / fscale) / fscale;
263                n = iclamp(k, 0, in_window_w - 1);
264                contrib[i].push_back(CONTRIB(n, weight));
265            }
266        }
267    }
268    else
269    {
270        for (i = 0; i < tmp_m_width; ++i)
271        {
272            contrib[i].resize(0);
273            center = (float) i / xscale;
274            left = int(ceilf(center - SUPPORT));
275            right = int(floorf(center + SUPPORT));
276            for(k = left; k <= right; ++k)
277            {
278                weight = center - (float) k;
279                weight = Mitchell_filter(weight);
280                n = iclamp(k, 0, in_window_w - 1);
281                contrib[i].push_back(CONTRIB(n, weight));
282            }
283        }
284    }
285
286    /* apply filter to zoom horizontally from src to tmp */
287    raster = (uint8_t*) my_calloc(in_window_w, 3);
288    for (k = 0; k < tmp_m_height; ++k)
289    {
290        get_row( ev, raster, in, int(floorf(in_x0)), in_window_w, k+in_y0 );
291        for (i = 0; i < tmp_m_width; ++i)
292        {
293            float       red = 0.0f;
294            float       green = 0.0f;
295            float       blue = 0.0f;
296            for(j = 0; j < contrib[i].size(); ++j)
297            {
298                int     pixel = contrib[i][j].pixel;
299                red     += raster[pixel * 3 + 0] * contrib[i][j].weight;
300                green   += raster[pixel * 3 + 1] * contrib[i][j].weight;
301                blue    += raster[pixel * 3 + 2] * contrib[i][j].weight;
302            }
303            put_pixel( ev, tmp, i, k, red, green, blue );
304        }
305    }
306    my_cfree(raster);
307
308    contrib.resize(out_height);
309
310    if (yscale < 1.0f)
311    {
312        width = SUPPORT / yscale;
313        fscale = 1.0f / yscale;
314        for (i = 0; i < out_height; ++i)
315        {
316            contrib[i].resize(0);
317
318            center = (float) i / yscale;
319            left = int(ceilf(center - width));
320            right = int(floorf(center + width));
321            for (k = left; k <= right; ++k)
322            {
323                weight = center - (float) k;
324                weight = Mitchell_filter(weight / fscale) / fscale;
325                n = iclamp(k, 0, tmp_m_height - 1);
326                contrib[i].push_back(CONTRIB(n, weight));
327            }
328        }
329    }
330    else
331    {
332        for (i = 0; i < out_height; ++i)
333        {
334            contrib[i].resize(0);
335            center = (float) i / yscale;
336            left = int(ceilf(center - SUPPORT));
337            right = int(floorf(center + SUPPORT));
338            for(k = left; k <= right; ++k)
339            {
340                weight = center - (float) k;
341                weight = Mitchell_filter(weight);
342                n = iclamp(k, 0, tmp_m_height - 1);
343                contrib[i].push_back(CONTRIB(n, weight));
344            }
345        }
346    }
347
348    /* apply filter to zoom vertically from tmp to dst */
349    raster = (uint8_t*) my_calloc(tmp_m_height, 3);
350    for (k = 0; k < tmp_m_width; ++k)
351    {
352        get_column( ev, raster, tmp, k );
353        for (i = 0; i < out_height; ++i)
354        {
355            float       red = 0.0f;
356            float       green = 0.0f;
357            float       blue = 0.0f;
358            for (j = 0; j < contrib[i].size(); ++j)
359            {
360                int     pixel = contrib[i][j].pixel;
361                red     += raster[pixel * 3 + 0] * contrib[i][j].weight;
362                green   += raster[pixel * 3 + 1] * contrib[i][j].weight;
363                blue    += raster[pixel * 3 + 2] * contrib[i][j].weight;
364            }
365            put_pixel( ev, out, k + out_x0, i + out_y0, red, green, blue );
366        }
367    }
368    my_cfree(raster);
369
370    contrib.resize(0);
371
372    delete tmp;
373}
374
375
Note: See TracBrowser for help on using the repository browser.