- Gfxprim - ucw.cz

[repo.or.cz] gfxprim.git branch master updated: 9b6212e7e083c6a78dca9cf322553c27baea66b2
by metan 23 Jan '12

23 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via 9b6212e7e083c6a78dca9cf322553c27baea66b2 (commit) from d17743a7cfa2617c5f50d47f801f8722480553fc (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/9b6212e7e083c6a78dca9cf322553c27baea… commit 9b6212e7e083c6a78dca9cf322553c27baea66b2 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Mon Jan 23 22:10:20 2012 +0100 filters: Add fast Bilinear + low pass filter. diff --git a/demos/fbshow/fbshow.c b/demos/fbshow/fbshow.c index 4d9f6bc..f6187f8 100644 --- a/demos/fbshow/fbshow.c +++ b/demos/fbshow/fbshow.c @@ -166,17 +166,17 @@ static void *image_loader(void *ptr) GP_Context *ret; - if (rat < 1) { +/* if (rat < 1) { cpu_timer_start(&timer, "Blur"); callback.priv = "Blurring Image"; if (GP_FilterGaussianBlur(img, img, 0.5/rat, 0.5/rat, &callback) == NULL) return NULL; cpu_timer_stop(&timer); - } + } */ cpu_timer_start(&timer, "Resampling"); callback.priv = "Resampling Image"; - ret = GP_FilterResize(img, NULL, GP_INTERP_CUBIC_INT, img->w * rat, img->h * rat, &callback); + ret = GP_FilterResize(img, NULL, GP_INTERP_LINEAR_LF_INT, img->w * rat, img->h * rat, &callback); GP_ContextFree(img); cpu_timer_stop(&timer); diff --git a/include/filters/GP_Resize.h b/include/filters/GP_Resize.h index 53cca92..fe6afbf 100644 --- a/include/filters/GP_Resize.h +++ b/include/filters/GP_Resize.h @@ -36,6 +36,12 @@ Faster than Bicubic, but less precise. + Bilinear LF + ~~~~~~~~~~~ + + Bilinear with low-pass filter on downscaling, this is the best choice for + fast up and downscaling. + Bicubic ~~~~~~~ @@ -51,10 +57,11 @@ #include "GP_Filter.h" typedef enum GP_InterpolationType { - GP_INTERP_NN, /* Nearest Neighbour */ - GP_INTERP_LINEAR_INT, /* Bilinear - fixed point arithmetics */ - GP_INTERP_CUBIC, /* Bicubic */ - GP_INTERP_CUBIC_INT, /* Bicubic - fixed point arithmetics */ + GP_INTERP_NN, /* Nearest Neighbour */ + GP_INTERP_LINEAR_INT, /* Bilinear - fixed point arithmetics */ + GP_INTERP_LINEAR_LF_INT, /* Bilinear + low pass filter on downscaling */ + GP_INTERP_CUBIC, /* Bicubic */ + GP_INTERP_CUBIC_INT, /* Bicubic - fixed point arithmetics */ } GP_InterpolationType; /* diff --git a/libs/filters/GP_Resize.c b/libs/filters/GP_Resize.c index a812599..76d1621 100644 --- a/libs/filters/GP_Resize.c +++ b/libs/filters/GP_Resize.c @@ -391,6 +391,171 @@ int GP_FilterInterpolate_CubicInt(const GP_Context *src, GP_Context *dst, return 0; } +/* + * Sample row. + * + * The x and y are starting coordinates in source image. + * + * The xpix_dist is distance of two sampled pixels in source image coordinates. + * + * The xoff is offset of the first pixel. + * + * The r, g, b are used to store resulting values. + */ +static inline void linear_lp_sample_x(const GP_Context *src, + uint32_t x, uint32_t y, + uint32_t xpix_dist, uint32_t xoff, + uint32_t *r, uint32_t *g, uint32_t *b) +{ + GP_Pixel pix; + uint32_t i; + + pix = GP_GetPixel_Raw_24BPP(src, x, y); + + *r = (GP_Pixel_GET_R_RGB888(pix) * xoff) >> 9; + *g = (GP_Pixel_GET_G_RGB888(pix) * xoff) >> 9; + *b = (GP_Pixel_GET_B_RGB888(pix) * xoff) >> 9; + + for (i = (1<<14) - xoff; i > xpix_dist; i -= xpix_dist) { + if (x < src->w - 1) + x++; + + pix = GP_GetPixel_Raw_24BPP(src, x, y); + + *r += (GP_Pixel_GET_R_RGB888(pix) * xpix_dist) >> 9; + *g += (GP_Pixel_GET_G_RGB888(pix) * xpix_dist) >> 9; + *b += (GP_Pixel_GET_B_RGB888(pix) * xpix_dist) >> 9; + } + + if (i > 0) { + if (x < src->w - 1) + x++; + + pix = GP_GetPixel_Raw_24BPP(src, x, y); + + *r += (GP_Pixel_GET_R_RGB888(pix) * i) >> 9; + *g += (GP_Pixel_GET_G_RGB888(pix) * i) >> 9; + *b += (GP_Pixel_GET_B_RGB888(pix) * i) >> 9; + } +} + +/* + * Linear interpolation with low-pass filtering, used for fast downscaling + * on both x and y. + * + * Basically we do weighted arithmetic mean of the pixels: + * + * [x, y], [x + 1, y], [x + 2, y] ... [x + k, y] + * [x, y + 1] + * [x, y + 2] . + * . . . + * . . . + * . . + * [x, y + l] .... [x + k, y + l] + * + * + * The parameter k respectively l is determined by the distance between + * sampled coordinates on x respectively y. + * + * The pixels are weighted by how much they are 'hit' by the rectangle defined + * by the sampled pixel. + * + * The implementation is inspired by imlib2 downscaling algorithm. + */ +static int interpolate_linear_lp_xy(const GP_Context *src, GP_Context *dst, + GP_ProgressCallback *callback) +{ + uint32_t xmap[dst->w + 1]; + uint32_t ymap[dst->h + 1]; + uint16_t xoff[dst->w + 1]; + uint16_t yoff[dst->h + 1]; + uint32_t x, y; + uint32_t i, j; + + /* Pre-compute mapping for interpolation */ + uint32_t xstep = (src->w << 16) / dst->w; + uint32_t xpix_dist = (dst->w << 14) / src->w; + + for (i = 0; i < dst->w + 1; i++) { + uint32_t val = i * xstep; + xmap[i] = val >> 16; + xoff[i] = ((255 - ((val >> 8) & 0xff)) * xpix_dist)>>8; + } + + uint32_t ystep = (src->h << 16) / dst->h; + uint32_t ypix_dist = (dst->h << 14) / src->h; + + for (i = 0; i < dst->h + 1; i++) { + uint32_t val = i * ystep; + ymap[i] = val >> 16; + yoff[i] = ((255 - ((val >> 8) & 0xff)) * ypix_dist)>>8; + } + + /* Interpolate */ + for (y = 0; y < dst->h; y++) { + for (x = 0; x < dst->w; x++) { + uint32_t r, g, b; + uint32_t r1, g1, b1; + uint32_t x0, y0; + + x0 = xmap[x]; + y0 = ymap[y]; + + linear_lp_sample_x(src, x0, y0, + xpix_dist, xoff[x], + &r, &g, &b); + + r = (r * yoff[y]) >> 14; + g = (g * yoff[y]) >> 14; + b = (b * yoff[y]) >> 14; + + for (j = (1<<14) - yoff[y]; j > ypix_dist; j -= ypix_dist) { + + x0 = xmap[x]; + + if (y0 < src->h - 1) + y0++; + + linear_lp_sample_x(src, x0, y0, + xpix_dist, xoff[x], + &r1, &g1, &b1); + + r += (r1 * ypix_dist) >> 14; + g += (g1 * ypix_dist) >> 14; + b += (b1 * ypix_dist) >> 14; + } + + if (j > 0) { + x0 = xmap[x]; + + if (y0 < src->h - 1) + y0++; + + linear_lp_sample_x(src, x0, y0, + xpix_dist, xoff[x], + &r1, &g1, &b1); + + r += (r1 * j) >> 14; + g += (g1 * j) >> 14; + b += (b1 * j) >> 14; + } + + r = (r + (1<<4))>>5; + g = (g + (1<<4))>>5; + b = (b + (1<<4))>>5; + + GP_PutPixel_Raw_24BPP(dst, x, y, + GP_Pixel_CREATE_RGB888(r, g, b)); + } + + if (GP_ProgressCallbackReport(callback, y, dst->h, dst->w)) + return 1; + } + + GP_ProgressCallbackDone(callback); + return 0; +} + int GP_FilterInterpolate_LinearInt(const GP_Context *src, GP_Context *dst, GP_ProgressCallback *callback) { @@ -398,7 +563,7 @@ int GP_FilterInterpolate_LinearInt(const GP_Context *src, GP_Context *dst, uint32_t ymap[dst->h + 1]; uint8_t xoff[dst->w + 1]; uint8_t yoff[dst->h + 1]; - uint32_t i, j; + uint32_t x, y, i; GP_DEBUG(1, "Scaling image %ux%u -> %ux%u %2.2f %2.2f", src->w, src->h, dst->w, dst->h, @@ -422,20 +587,20 @@ int GP_FilterInterpolate_LinearInt(const GP_Context *src, GP_Context *dst, } /* Interpolate */ - for (i = 0; i < dst->h; i++) { - for (j = 0; j < dst->w; j++) { + for (y = 0; y < dst->h; y++) { + for (x = 0; x < dst->w; x++) { GP_Pixel pix00, pix01, pix10, pix11; uint32_t r0, r1, g0, g1, b0, b1; GP_Coord x0, x1, y0, y1; - x0 = xmap[j]; - x1 = xmap[j] + 1; + x0 = xmap[x]; + x1 = xmap[x] + 1; if (x1 >= (GP_Coord)src->w) x1 = src->w - 1; - y0 = ymap[i]; - y1 = ymap[i] + 1; + y0 = ymap[y]; + y1 = ymap[y] + 1; if (y1 >= (GP_Coord)src->h) y1 = src->h - 1; @@ -445,31 +610,31 @@ int GP_FilterInterpolate_LinearInt(const GP_Context *src, GP_Context *dst, pix01 = GP_GetPixel_Raw_24BPP(src, x0, y1); pix11 = GP_GetPixel_Raw_24BPP(src, x1, y1); - r0 = GP_Pixel_GET_R_RGB888(pix00) * (255 - xoff[j]); - g0 = GP_Pixel_GET_G_RGB888(pix00) * (255 - xoff[j]); - b0 = GP_Pixel_GET_B_RGB888(pix00) * (255 - xoff[j]); + r0 = GP_Pixel_GET_R_RGB888(pix00) * (255 - xoff[x]); + g0 = GP_Pixel_GET_G_RGB888(pix00) * (255 - xoff[x]); + b0 = GP_Pixel_GET_B_RGB888(pix00) * (255 - xoff[x]); - r0 += GP_Pixel_GET_R_RGB888(pix10) * xoff[j]; - g0 += GP_Pixel_GET_G_RGB888(pix10) * xoff[j]; - b0 += GP_Pixel_GET_B_RGB888(pix10) * xoff[j]; + r0 += GP_Pixel_GET_R_RGB888(pix10) * xoff[x]; + g0 += GP_Pixel_GET_G_RGB888(pix10) * xoff[x]; + b0 += GP_Pixel_GET_B_RGB888(pix10) * xoff[x]; - r1 = GP_Pixel_GET_R_RGB888(pix01) * (255 - xoff[j]); - g1 = GP_Pixel_GET_G_RGB888(pix01) * (255 - xoff[j]); - b1 = GP_Pixel_GET_B_RGB888(pix01) * (255 - xoff[j]); + r1 = GP_Pixel_GET_R_RGB888(pix01) * (255 - xoff[x]); + g1 = GP_Pixel_GET_G_RGB888(pix01) * (255 - xoff[x]); + b1 = GP_Pixel_GET_B_RGB888(pix01) * (255 - xoff[x]); - r1 += GP_Pixel_GET_R_RGB888(pix11) * xoff[j]; - g1 += GP_Pixel_GET_G_RGB888(pix11) * xoff[j]; - b1 += GP_Pixel_GET_B_RGB888(pix11) * xoff[j]; + r1 += GP_Pixel_GET_R_RGB888(pix11) * xoff[x]; + g1 += GP_Pixel_GET_G_RGB888(pix11) * xoff[x]; + b1 += GP_Pixel_GET_B_RGB888(pix11) * xoff[x]; - r0 = (r1 * yoff[i] + r0 * (255 - yoff[i])) >> 16; - g0 = (g1 * yoff[i] + g0 * (255 - yoff[i])) >> 16; - b0 = (b1 * yoff[i] + b0 * (255 - yoff[i])) >> 16; + r0 = (r1 * yoff[y] + r0 * (255 - yoff[y]) + (1<<15)) >> 16; + g0 = (g1 * yoff[y] + g0 * (255 - yoff[y]) + (1<<15)) >> 16; + b0 = (b1 * yoff[y] + b0 * (255 - yoff[y]) + (1<<15)) >> 16; - GP_PutPixel_Raw_24BPP(dst, j, i, + GP_PutPixel_Raw_24BPP(dst, x, y, GP_Pixel_CREATE_RGB888(r0, g0, b0)); } - if (GP_ProgressCallbackReport(callback, i, dst->h, dst->w)) + if (GP_ProgressCallbackReport(callback, y, dst->h, dst->w)) return 1; } @@ -477,6 +642,24 @@ int GP_FilterInterpolate_LinearInt(const GP_Context *src, GP_Context *dst, return 0; } +int GP_FilterInterpolate_LinearLFInt(const GP_Context *src, GP_Context *dst, + GP_ProgressCallback *callback) +{ + float x_rat = 1.00 * dst->w / src->w; + float y_rat = 1.00 * dst->h / src->h; + + if (x_rat < 1.00 && y_rat < 1.00) { + GP_DEBUG(1, "Downscaling image %ux%u -> %ux%u %2.2f %2.2f", + src->w, src->h, dst->w, dst->h, x_rat, y_rat); + return interpolate_linear_lp_xy(src, dst, callback); + } + + //TODO: x_rat > 1.00 && y_rat < 1.00 + //TODO: x_rat < 1.00 && y_rat > 1.00 + + return GP_FilterInterpolate_LinearInt(src, dst, callback); +} + int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, GP_InterpolationType type, GP_ProgressCallback *callback) @@ -486,6 +669,8 @@ int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, return GP_FilterInterpolate_NN(src, dst, callback); case GP_INTERP_LINEAR_INT: return GP_FilterInterpolate_LinearInt(src, dst, callback); + case GP_INTERP_LINEAR_LF_INT: + return GP_FilterInterpolate_LinearLFInt(src, dst, callback); case GP_INTERP_CUBIC: return GP_FilterInterpolate_Cubic(src, dst, callback); case GP_INTERP_CUBIC_INT: ----------------------------------------------------------------------- Summary of changes: demos/fbshow/fbshow.c | 6 +- include/filters/GP_Resize.h | 15 ++- libs/filters/GP_Resize.c | 233 ++++++++++++++++++++++++++++++++++++++----- 3 files changed, 223 insertions(+), 31 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: d17743a7cfa2617c5f50d47f801f8722480553fc
by metan 23 Jan '12

23 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via d17743a7cfa2617c5f50d47f801f8722480553fc (commit) from 66cbe1a3b25ba7d4c3361d420fa9de364c02557e (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/d17743a7cfa2617c5f50d47f801f87224805… commit d17743a7cfa2617c5f50d47f801f8722480553fc Author: Cyril Hrubis <metan(a)ucw.cz> Date: Mon Jan 23 16:57:24 2012 +0100 filters: Added Bilinear interpolation. diff --git a/include/filters/GP_Resize.h b/include/filters/GP_Resize.h index fcae054..53cca92 100644 --- a/include/filters/GP_Resize.h +++ b/include/filters/GP_Resize.h @@ -31,6 +31,10 @@ edges mostly consisting of big one color regions (eg doesn't blur the result on upscaling). + Bilinear + ~~~~~~~~ + + Faster than Bicubic, but less precise. Bicubic ~~~~~~~ @@ -47,9 +51,10 @@ #include "GP_Filter.h" typedef enum GP_InterpolationType { - GP_INTERP_NN, /* Nearest Neighbour */ - GP_INTERP_CUBIC, /* Bicubic */ - GP_INTERP_CUBIC_INT, /* Bicubic - integer arithmetics */ + GP_INTERP_NN, /* Nearest Neighbour */ + GP_INTERP_LINEAR_INT, /* Bilinear - fixed point arithmetics */ + GP_INTERP_CUBIC, /* Bicubic */ + GP_INTERP_CUBIC_INT, /* Bicubic - fixed point arithmetics */ } GP_InterpolationType; /* diff --git a/libs/filters/GP_Resize.c b/libs/filters/GP_Resize.c index a29ec0d..a812599 100644 --- a/libs/filters/GP_Resize.c +++ b/libs/filters/GP_Resize.c @@ -32,18 +32,30 @@ int GP_FilterInterpolate_NN(const GP_Context *src, GP_Context *dst, GP_ProgressCallback *callback) { + uint32_t xmap[dst->w]; + uint32_t ymap[dst->h]; + uint32_t i; GP_Coord x, y; GP_DEBUG(1, "Scaling image %ux%u -> %ux%u %2.2f %2.2f", src->w, src->h, dst->w, dst->h, 1.00 * dst->w / src->w, 1.00 * dst->h / src->h); + /* Pre-compute mapping for interpolation */ + uint32_t xstep = (src->w << 16) / dst->w; + + for (i = 0; i < dst->w + 1; i++) + xmap[i] = ((i * xstep) + (1<<15)) >> 16; + + uint32_t ystep = (src->h << 16) / dst->h; + + for (i = 0; i < dst->h + 1; i++) + ymap[i] = ((i * ystep) + (1<<15)) >> 16; + + /* Interpolate */ for (y = 0; y < (GP_Coord)dst->h; y++) { for (x = 0; x < (GP_Coord)dst->w; x++) { - GP_Coord xi = (1.00 * x / dst->w) * src->w; - GP_Coord yi = (1.00 * y / dst->h) * src->h; - - GP_Pixel pix = GP_GetPixel_Raw_24BPP(src, xi, yi); + GP_Pixel pix = GP_GetPixel_Raw_24BPP(src, xmap[x], ymap[y]); GP_PutPixel_Raw_24BPP(dst, x, y, pix); } @@ -379,6 +391,92 @@ int GP_FilterInterpolate_CubicInt(const GP_Context *src, GP_Context *dst, return 0; } +int GP_FilterInterpolate_LinearInt(const GP_Context *src, GP_Context *dst, + GP_ProgressCallback *callback) +{ + uint32_t xmap[dst->w + 1]; + uint32_t ymap[dst->h + 1]; + uint8_t xoff[dst->w + 1]; + uint8_t yoff[dst->h + 1]; + uint32_t i, j; + + GP_DEBUG(1, "Scaling image %ux%u -> %ux%u %2.2f %2.2f", + src->w, src->h, dst->w, dst->h, + 1.00 * dst->w / src->w, 1.00 * dst->h / src->h); + + /* Pre-compute mapping for interpolation */ + uint32_t xstep = (src->w << 16) / dst->w; + + for (i = 0; i < dst->w + 1; i++) { + uint32_t val = i * xstep; + xmap[i] = val >> 16; + xoff[i] = (val >> 8) & 0xff; + } + + uint32_t ystep = (src->h << 16) / dst->h; + + for (i = 0; i < dst->h + 1; i++) { + uint32_t val = i * ystep; + ymap[i] = val >> 16; + yoff[i] = (val >> 8) & 0xff; + } + + /* Interpolate */ + for (i = 0; i < dst->h; i++) { + for (j = 0; j < dst->w; j++) { + GP_Pixel pix00, pix01, pix10, pix11; + uint32_t r0, r1, g0, g1, b0, b1; + GP_Coord x0, x1, y0, y1; + + x0 = xmap[j]; + x1 = xmap[j] + 1; + + if (x1 >= (GP_Coord)src->w) + x1 = src->w - 1; + + y0 = ymap[i]; + y1 = ymap[i] + 1; + + if (y1 >= (GP_Coord)src->h) + y1 = src->h - 1; + + pix00 = GP_GetPixel_Raw_24BPP(src, x0, y0); + pix10 = GP_GetPixel_Raw_24BPP(src, x1, y0); + pix01 = GP_GetPixel_Raw_24BPP(src, x0, y1); + pix11 = GP_GetPixel_Raw_24BPP(src, x1, y1); + + r0 = GP_Pixel_GET_R_RGB888(pix00) * (255 - xoff[j]); + g0 = GP_Pixel_GET_G_RGB888(pix00) * (255 - xoff[j]); + b0 = GP_Pixel_GET_B_RGB888(pix00) * (255 - xoff[j]); + + r0 += GP_Pixel_GET_R_RGB888(pix10) * xoff[j]; + g0 += GP_Pixel_GET_G_RGB888(pix10) * xoff[j]; + b0 += GP_Pixel_GET_B_RGB888(pix10) * xoff[j]; + + r1 = GP_Pixel_GET_R_RGB888(pix01) * (255 - xoff[j]); + g1 = GP_Pixel_GET_G_RGB888(pix01) * (255 - xoff[j]); + b1 = GP_Pixel_GET_B_RGB888(pix01) * (255 - xoff[j]); + + r1 += GP_Pixel_GET_R_RGB888(pix11) * xoff[j]; + g1 += GP_Pixel_GET_G_RGB888(pix11) * xoff[j]; + b1 += GP_Pixel_GET_B_RGB888(pix11) * xoff[j]; + + r0 = (r1 * yoff[i] + r0 * (255 - yoff[i])) >> 16; + g0 = (g1 * yoff[i] + g0 * (255 - yoff[i])) >> 16; + b0 = (b1 * yoff[i] + b0 * (255 - yoff[i])) >> 16; + + GP_PutPixel_Raw_24BPP(dst, j, i, + GP_Pixel_CREATE_RGB888(r0, g0, b0)); + } + + if (GP_ProgressCallbackReport(callback, i, dst->h, dst->w)) + return 1; + } + + GP_ProgressCallbackDone(callback); + return 0; +} + int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, GP_InterpolationType type, GP_ProgressCallback *callback) @@ -386,6 +484,8 @@ int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, switch (type) { case GP_INTERP_NN: return GP_FilterInterpolate_NN(src, dst, callback); + case GP_INTERP_LINEAR_INT: + return GP_FilterInterpolate_LinearInt(src, dst, callback); case GP_INTERP_CUBIC: return GP_FilterInterpolate_Cubic(src, dst, callback); case GP_INTERP_CUBIC_INT: ----------------------------------------------------------------------- Summary of changes: include/filters/GP_Resize.h | 11 +++- libs/filters/GP_Resize.c | 108 +++++++++++++++++++++++++++++++++++++++++-- 2 files changed, 112 insertions(+), 7 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: 66cbe1a3b25ba7d4c3361d420fa9de364c02557e
by bluebear 22 Jan '12

22 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via 66cbe1a3b25ba7d4c3361d420fa9de364c02557e (commit) from 446922190513e6535202e8f45ab897bc622390ee (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/66cbe1a3b25ba7d4c3361d420fa9de364c02… commit 66cbe1a3b25ba7d4c3361d420fa9de364c02557e Author: Jiri BlueBear Dluhos <jiri.bluebear.dluhos(a)gmail.com> Date: Sun Jan 22 19:46:51 2012 +0100 First implementation of partial ellipse (plus some angle-related utils). diff --git a/include/core/GP_AngleUtils.h b/include/core/GP_AngleUtils.h new file mode 100644 index 0000000..f896b37 --- /dev/null +++ b/include/core/GP_AngleUtils.h @@ -0,0 +1,30 @@ +/***************************************************************************** + * This file is part of gfxprim library. * + * * + * Gfxprim is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Lesser General Public * + * License as published by the Free Software Foundation; either * + * version 2.1 of the License, or (at your option) any later version. * + * * + * Gfxprim is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * + * Lesser General Public License for more details. * + * * + * You should have received a copy of the GNU Lesser General Public * + * License along with gfxprim; if not, write to the Free Software * + * Foundation, Inc., 51 Franklin Street, Fifth Floor, * + * Boston, MA 02110-1301 USA * + * * + * Copyright (C) 2009-2011 Cyril Hrubis <metan(a)ucw.cz> * + * Copyright (c) 2012 Jiri Dluhos <jiri.bluebear.dluhos(a)gmail.com> * + * * + *****************************************************************************/ + +#ifndef GP_ANGLE_UTILS_H +#define GP_ANGLE_UTILS_H + +double GP_NormalizeAngle(double phi); +int GP_AngleInRange(double angle, double start, double end); + +#endif No newline at end of file diff --git a/libs/core/GP_AngleUtils.c b/libs/core/GP_AngleUtils.c new file mode 100644 index 0000000..f1c279b --- /dev/null +++ b/libs/core/GP_AngleUtils.c @@ -0,0 +1,49 @@ +/***************************************************************************** + * This file is part of gfxprim library. * + * * + * Gfxprim is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Lesser General Public * + * License as published by the Free Software Foundation; either * + * version 2.1 of the License, or (at your option) any later version. * + * * + * Gfxprim is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * + * Lesser General Public License for more details. * + * * + * You should have received a copy of the GNU Lesser General Public * + * License along with gfxprim; if not, write to the Free Software * + * Foundation, Inc., 51 Franklin Street, Fifth Floor, * + * Boston, MA 02110-1301 USA * + * * + * Copyright (C) 2011 Tomas Gavenciak <gavento(a)ucw.cz> * + * Copyright (C) 2012 Cyril Hrubis <metan(a)ucw.cz> * + * Copyright (C) 2012 Jiri Dluhos <jiri.bluebear.dluhos(a)gmail.com> * + * * + *****************************************************************************/ + +#include "GP_AngleUtils.h" + +#include <math.h> + +double GP_NormalizeAngle(double phi) +{ + // clamp angle to <-2*pi, 2*pi> + double phi2 = fmod(phi, 2*M_PI); + + // clamp angle to <0, 2*pi> + if (phi2 < 0) + phi2 += 2*M_PI; + + return phi2; +} + +int GP_AngleInRange(double angle, double start, double end) +{ + if (start < end) { + return (angle >= start && angle <= end); + } else { + return (angle >= start && angle <= 2*M_PI) + || (angle >= 0 && angle <= end); + } +} No newline at end of file diff --git a/libs/gfx/GP_PartialEllipse.c b/libs/gfx/GP_PartialEllipse.c new file mode 100644 index 0000000..5415a58 --- /dev/null +++ b/libs/gfx/GP_PartialEllipse.c @@ -0,0 +1,53 @@ +/***************************************************************************** + * This file is part of gfxprim library. * + * * + * Gfxprim is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Lesser General Public * + * License as published by the Free Software Foundation; either * + * version 2.1 of the License, or (at your option) any later version. * + * * + * Gfxprim is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * + * Lesser General Public License for more details. * + * * + * You should have received a copy of the GNU Lesser General Public * + * License along with gfxprim; if not, write to the Free Software * + * Foundation, Inc., 51 Franklin Street, Fifth Floor, * + * Boston, MA 02110-1301 USA * + * * + * Copyright (C) 2009-2012 Jiri "BlueBear" Dluhos * + * <jiri.bluebear.dluhos(a)gmail.com> * + * * + * Copyright (C) 2009-2012 Cyril Hrubis <metan(a)ucw.cz> * + * * + *****************************************************************************/ + +#include "GP_Gfx.h" +#include "core/GP_FnPerBpp.h" + +#include "algo/PartialEllipse.algo.h" + +/* Generate drawing functions for various bit depths. */ +GP_DEF_DRAW_FN_PER_BPP(GP_PartialEllipse_Raw, DEF_PARTIAL_ELLIPSE_FN) + +void GP_PartialEllipse_Raw(GP_Context *context, GP_Coord xcenter, GP_Coord ycenter, + GP_Size a, GP_Size b, int start, int end, GP_Pixel pixel) +{ + GP_CHECK_CONTEXT(context); + + GP_FN_PER_BPP_CONTEXT(GP_PartialEllipse_Raw, context, context, + xcenter, ycenter, a, b, start, end, pixel); +} + +void GP_PartialEllipse(GP_Context *context, GP_Coord xcenter, GP_Coord ycenter, + GP_Size a, GP_Size b, int start, int end, GP_Pixel pixel) +{ + GP_CHECK_CONTEXT(context); + + /* recalculate center point and swap a and b when axes are swapped */ + GP_TRANSFORM_POINT(context, xcenter, ycenter); + GP_TRANSFORM_SWAP(context, a, b); + + GP_PartialEllipse_Raw(context, xcenter, ycenter, a, b, start, end, pixel); +} No newline at end of file diff --git a/libs/gfx/algo/PartialEllipse.algo.h b/libs/gfx/algo/PartialEllipse.algo.h new file mode 100644 index 0000000..b87826e --- /dev/null +++ b/libs/gfx/algo/PartialEllipse.algo.h @@ -0,0 +1,68 @@ +/***************************************************************************** + * This file is part of gfxprim library. * + * * + * Gfxprim is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Lesser General Public * + * License as published by the Free Software Foundation; either * + * version 2.1 of the License, or (at your option) any later version. * + * * + * Gfxprim is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * + * Lesser General Public License for more details. * + * * + * You should have received a copy of the GNU Lesser General Public * + * License along with gfxprim; if not, write to the Free Software * + * Foundation, Inc., 51 Franklin Street, Fifth Floor, * + * Boston, MA 02110-1301 USA * + * * + * Copyright (C) 2009-2010 Jiri "BlueBear" Dluhos * + * <jiri.bluebear.dluhos(a)gmail.com> * + * * + * Copyright (C) 2009-2010 Cyril Hrubis <metan(a)ucw.cz> * + * * + *****************************************************************************/ + +#include "core/GP_AngleUtils.h" + +#include <math.h> + +/* + * This macro defines a partial ellipse drawing function. + * Arguments: + * CONTEXT_T - user-defined type of drawing context (passed to PUTPIXEL) + * PIXVAL_T - user-defined pixel value type (passed to PUTPIXEL) + * PUTPIXEL - a pixel drawing function f(context, x, y, pixval) + * FN_NAME - name of the function to be defined + */ +#define DEF_PARTIAL_ELLIPSE_FN(FN_NAME, CONTEXT_T, PIXVAL_T, PUTPIXEL) +void FN_NAME(CONTEXT_T context, int xcenter, int ycenter, int a, int b, + int start, int end, PIXVAL_T pixval) +{ + double startAngle = GP_NormalizeAngle(2*M_PI*(start / 360000.0)); + double endAngle = GP_NormalizeAngle(2*M_PI*(end / 360000.0)); ++ int x; + for (x = -a; x <= a; x++) { + double angle = acos(((double) x) / a); + double y = floor(b*sin(angle)); + if (GP_AngleInRange(angle, startAngle, endAngle)) { + PUTPIXEL(context, xcenter+x, ycenter-y, pixval); + } + if (GP_AngleInRange(2*M_PI - angle, startAngle, endAngle)) { + PUTPIXEL(context, xcenter+x, ycenter+y, pixval); + } + } ++ int y; + for (y = -b; y <= b; y++) { + double angle = asin(((double) y) / b); + double x = floor(a*cos(angle)); + if (GP_AngleInRange(angle, startAngle, endAngle)) { + PUTPIXEL(context, xcenter+x, ycenter-y, pixval); + } + if (GP_AngleInRange(M_PI - angle, startAngle, endAngle)) { + PUTPIXEL(context, x center-x, ycenter-y, pixval); + } + } +} No newline at end of file ----------------------------------------------------------------------- Summary of changes: .../GP_DefaultFont.h => core/GP_AngleUtils.h} | 14 ++-- libs/core/{GP_GetPutPixel.c => GP_AngleUtils.c} | 38 ++++++----- libs/gfx/{GP_Arc.c => GP_PartialEllipse.c} | 34 ++++------ .../algo/{Line.algo.h => PartialEllipse.algo.h} | 72 ++++++++------------ 4 files changed, 69 insertions(+), 89 deletions(-) copy include/{text/GP_DefaultFont.h => core/GP_AngleUtils.h} (87%) copy libs/core/{GP_GetPutPixel.c => GP_AngleUtils.c} (74%) copy libs/gfx/{GP_Arc.c => GP_PartialEllipse.c} (73%) copy libs/gfx/algo/{Line.algo.h => PartialEllipse.algo.h} (66%) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: 446922190513e6535202e8f45ab897bc622390ee
by metan 22 Jan '12

22 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via 446922190513e6535202e8f45ab897bc622390ee (commit) from 09241d96413a6376d399b4e81219fdf2d9fb1462 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/446922190513e6535202e8f45ab897bc6223… commit 446922190513e6535202e8f45ab897bc622390ee Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sun Jan 22 19:05:00 2012 +0100 filters: Blur fix row/column prefill. diff --git a/libs/filters/GP_Linear.c b/libs/filters/GP_Linear.c index 447e805..a55d1fe 100644 --- a/libs/filters/GP_Linear.c +++ b/libs/filters/GP_Linear.c @@ -154,7 +154,7 @@ int GP_FilterHLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, GP_Coord x, y; uint32_t i; int32_t ikernel[kw]; - uint32_t size = dst->w + kw; + uint32_t size = dst->w + kw - 1; for (i = 0; i < kw; i++) ikernel[i] = kernel[i] * MUL + 0.5; @@ -179,15 +179,15 @@ int GP_FilterHLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, B[i] = GP_Pixel_GET_B_RGB888(pix); } - for (i = kw/2; i < src->w; i++) { + for (i = 0; i < src->w; i++) { pix = GP_GetPixel_Raw_24BPP(src, i, y); - R[i] = GP_Pixel_GET_R_RGB888(pix); - G[i] = GP_Pixel_GET_G_RGB888(pix); - B[i] = GP_Pixel_GET_B_RGB888(pix); + R[i+kw/2] = GP_Pixel_GET_R_RGB888(pix); + G[i+kw/2] = GP_Pixel_GET_G_RGB888(pix); + B[i+kw/2] = GP_Pixel_GET_B_RGB888(pix); } - for (i = src->w; i < size; i++) { + for (i = src->w + kw/2; i < size; i++) { R[i] = GP_Pixel_GET_R_RGB888(pix); G[i] = GP_Pixel_GET_G_RGB888(pix); B[i] = GP_Pixel_GET_B_RGB888(pix); @@ -233,7 +233,7 @@ int GP_FilterVLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, GP_Coord x, y; uint32_t i; int32_t ikernel[kh]; - uint32_t size = dst->h + kh; + uint32_t size = dst->h + kh - 1; for (i = 0; i < kh; i++) ikernel[i] = kernel[i] * MUL + 0.5; @@ -258,15 +258,15 @@ int GP_FilterVLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, B[i] = GP_Pixel_GET_B_RGB888(pix); } - for (i = kh/2; i < src->h; i++) { + for (i = 0; i < src->h; i++) { pix = GP_GetPixel_Raw_24BPP(src, x, i); - R[i] = GP_Pixel_GET_R_RGB888(pix); - G[i] = GP_Pixel_GET_G_RGB888(pix); - B[i] = GP_Pixel_GET_B_RGB888(pix); + R[i+kh/2] = GP_Pixel_GET_R_RGB888(pix); + G[i+kh/2] = GP_Pixel_GET_G_RGB888(pix); + B[i+kh/2] = GP_Pixel_GET_B_RGB888(pix); } - for (i = src->h; i < size; i++) { + for (i = src->h + kh/2; i < size; i++) { R[i] = GP_Pixel_GET_R_RGB888(pix); G[i] = GP_Pixel_GET_G_RGB888(pix); B[i] = GP_Pixel_GET_B_RGB888(pix); ----------------------------------------------------------------------- Summary of changes: libs/filters/GP_Linear.c | 24 ++++++++++++------------ 1 files changed, 12 insertions(+), 12 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: 09241d96413a6376d399b4e81219fdf2d9fb1462
by metan 22 Jan '12

22 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via 09241d96413a6376d399b4e81219fdf2d9fb1462 (commit) from 1d97b3f64cb80fea14245d253f6c7db64b8495b2 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/09241d96413a6376d399b4e81219fdf2d9fb… commit 09241d96413a6376d399b4e81219fdf2d9fb1462 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sun Jan 22 18:09:25 2012 +0100 filters: BiCubic now the coordinate mapping is right. diff --git a/libs/filters/GP_Linear.c b/libs/filters/GP_Linear.c index 7101472..447e805 100644 --- a/libs/filters/GP_Linear.c +++ b/libs/filters/GP_Linear.c @@ -16,7 +16,7 @@ * Foundation, Inc., 51 Franklin Street, Fifth Floor, * * Boston, MA 02110-1301 USA * * * - * Copyright (C) 2009-2011 Cyril Hrubis <metan(a)ucw.cz> * + * Copyright (C) 2009-2012 Cyril Hrubis <metan(a)ucw.cz> * * * *****************************************************************************/ diff --git a/libs/filters/GP_Resize.c b/libs/filters/GP_Resize.c index 7b26eaa..a29ec0d 100644 --- a/libs/filters/GP_Resize.c +++ b/libs/filters/GP_Resize.c @@ -16,10 +16,12 @@ * Foundation, Inc., 51 Franklin Street, Fifth Floor, * * Boston, MA 02110-1301 USA * * * - * Copyright (C) 2009-2011 Cyril Hrubis <metan(a)ucw.cz> * + * Copyright (C) 2009-2012 Cyril Hrubis <metan(a)ucw.cz> * * * *****************************************************************************/ +#include <math.h> + #include <GP_Context.h> #include <GP_GetPutPixel.h> @@ -92,6 +94,13 @@ typedef union v4f { #define SUM_V4SF(a) ((a).f[0] + (a).f[1] + (a).f[2] + (a).f[3]) +#define CLAMP(val) do { + if (val < 0) + val = 0; + if (val > 255) + val = 255; +} while (0) + int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, GP_ProgressCallback *callback) { @@ -105,28 +114,36 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, for (i = 0; i < dst->w; i++) { float x = (1.00 * i / dst->w) * src->w; v4f cvx; - int xi = x - 1; + int xi[4]; + + xi[0] = floor(x - 1); + xi[1] = x; + xi[2] = x + 1; + xi[3] = x + 2; + + cvx.f[0] = cubic(xi[0] - x); + cvx.f[1] = cubic(xi[1] - x); + cvx.f[2] = cubic(xi[2] - x); + cvx.f[3] = cubic(xi[3] - x); - if (xi < 0) - xi = 0; + if (xi[0] < 0) + xi[0] = 0; - if (xi > (int)src->w - 4) - xi = src->w - 4; + if (xi[2] >= (int)src->w) + xi[2] = src->w - 1; - cvx.f[0] = cubic(xi - x); - cvx.f[1] = cubic(xi + 1 - x); - cvx.f[2] = cubic(xi + 2 - x); - cvx.f[3] = cubic(xi + 3 - x); + if (xi[3] >= (int)src->w) + xi[3] = src->w - 1; /* Generate interpolated column */ for (j = 0; j < src->h; j++) { v4f rv, gv, bv; GP_Pixel pix[4]; - pix[0] = GP_GetPixel_Raw_24BPP(src, xi, j); - pix[1] = GP_GetPixel_Raw_24BPP(src, xi + 1, j); - pix[2] = GP_GetPixel_Raw_24BPP(src, xi + 2, j); - pix[3] = GP_GetPixel_Raw_24BPP(src, xi + 3, j); + pix[0] = GP_GetPixel_Raw_24BPP(src, xi[0], j); + pix[1] = GP_GetPixel_Raw_24BPP(src, xi[1], j); + pix[2] = GP_GetPixel_Raw_24BPP(src, xi[2], j); + pix[3] = GP_GetPixel_Raw_24BPP(src, xi[3], j); rv.f[0] = GP_Pixel_GET_R_RGB888(pix[0]); rv.f[1] = GP_Pixel_GET_R_RGB888(pix[1]); @@ -157,33 +174,41 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, float y = (1.00 * j / dst->h) * src->h; v4f cvy, rv, gv, bv; float r, g, b; - int yi = y - 1; - - if (yi < 0) - yi = 0; - - if (yi > (int)src->h - 4) - yi = src->h - 4; + int yi[4]; + + yi[0] = floor(y - 1); + yi[1] = y; + yi[2] = y + 1; + yi[3] = y + 2; + + cvy.f[0] = cubic(yi[0] - y); + cvy.f[1] = cubic(yi[1] - y); + cvy.f[2] = cubic(yi[2] - y); + cvy.f[3] = cubic(yi[3] - y); + + if (yi[0] < 0) + yi[0] = 0; + + if (yi[2] >= (int)src->h) + yi[2] = src->h - 1; + + if (yi[3] >= (int)src->h) + yi[3] = src->h - 1; - cvy.f[0] = cubic(yi - y); - cvy.f[1] = cubic(yi + 1 - y); - cvy.f[2] = cubic(yi + 2 - y); - cvy.f[3] = cubic(yi + 3 - y); - - rv.f[0] = col_r[yi]; - rv.f[1] = col_r[yi + 1]; - rv.f[2] = col_r[yi + 2]; - rv.f[3] = col_r[yi + 3]; + rv.f[0] = col_r[yi[0]]; + rv.f[1] = col_r[yi[1]]; + rv.f[2] = col_r[yi[2]]; + rv.f[3] = col_r[yi[3]]; - gv.f[0] = col_g[yi]; - gv.f[1] = col_g[yi + 1]; - gv.f[2] = col_g[yi + 2]; - gv.f[3] = col_g[yi + 3]; + gv.f[0] = col_g[yi[0]]; + gv.f[1] = col_g[yi[1]]; + gv.f[2] = col_g[yi[2]]; + gv.f[3] = col_g[yi[3]]; - bv.f[0] = col_b[yi]; - bv.f[1] = col_b[yi + 1]; - bv.f[2] = col_b[yi + 2]; - bv.f[3] = col_b[yi + 3]; + bv.f[0] = col_b[yi[0]]; + bv.f[1] = col_b[yi[1]]; + bv.f[2] = col_b[yi[2]]; + bv.f[3] = col_b[yi[3]]; rv = MUL_V4SF(rv, cvy); gv = MUL_V4SF(gv, cvy); @@ -193,24 +218,10 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, g = SUM_V4SF(gv); b = SUM_V4SF(bv); - if (r > 255) - r = 255; - - if (r < 0) - r = 0; - - if (g > 255) - g = 255; - - if (g < 0) - g = 0; - - if (b > 255) - b = 255; - - if (b < 0) - b = 0; - + CLAMP(r); + CLAMP(g); + CLAMP(b); + GP_Pixel pix = GP_Pixel_CREATE_RGB888((uint8_t)r, (uint8_t)g, (uint8_t)b); GP_PutPixel_Raw_24BPP(dst, i, j, pix); } @@ -238,137 +249,6 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, int GP_FilterInterpolate_CubicInt(const GP_Context *src, GP_Context *dst, GP_ProgressCallback *callback) { - int32_t col_r[src->h], col_g[src->h], col_b[src->h]; - uint32_t i, j; - - GP_DEBUG(1, "Scaling image %ux%u -> %ux%u %2.2f %2.2f", - src->w, src->h, dst->w, dst->h, - 1.00 * dst->w / src->w, 1.00 * dst->h / src->h); - - for (i = 0; i < dst->w; i++) { - float x = (1.00 * i / dst->w) * src->w; - int32_t cvx[4]; - int xi = x - 1; - - if (xi < 0) - xi = 0; - - if (xi > (int)src->w - 4) - xi = src->w - 4; - - cvx[0] = cubic(xi - x) * MUL + 0.5; - cvx[1] = cubic(xi + 1 - x) * MUL + 0.5; - cvx[2] = cubic(xi + 2 - x) * MUL + 0.5; - cvx[3] = cubic(xi + 3 - x) * MUL + 0.5; - - /* Generate interpolated column */ - for (j = 0; j < src->h; j++) { - int32_t rv[4], gv[4], bv[4]; - GP_Pixel pix[4]; - - pix[0] = GP_GetPixel_Raw_24BPP(src, xi, j); - pix[1] = GP_GetPixel_Raw_24BPP(src, xi + 1, j); - pix[2] = GP_GetPixel_Raw_24BPP(src, xi + 2, j); - pix[3] = GP_GetPixel_Raw_24BPP(src, xi + 3, j); - - rv[0] = GP_Pixel_GET_R_RGB888(pix[0]); - rv[1] = GP_Pixel_GET_R_RGB888(pix[1]); - rv[2] = GP_Pixel_GET_R_RGB888(pix[2]); - rv[3] = GP_Pixel_GET_R_RGB888(pix[3]); - - gv[0] = GP_Pixel_GET_G_RGB888(pix[0]); - gv[1] = GP_Pixel_GET_G_RGB888(pix[1]); - gv[2] = GP_Pixel_GET_G_RGB888(pix[2]); - gv[3] = GP_Pixel_GET_G_RGB888(pix[3]); - - bv[0] = GP_Pixel_GET_B_RGB888(pix[0]); - bv[1] = GP_Pixel_GET_B_RGB888(pix[1]); - bv[2] = GP_Pixel_GET_B_RGB888(pix[2]); - bv[3] = GP_Pixel_GET_B_RGB888(pix[3]); - - MUL_I(rv, cvx); - MUL_I(gv, cvx); - MUL_I(bv, cvx); - - col_r[j] = SUM_I(rv); - col_g[j] = SUM_I(gv); - col_b[j] = SUM_I(bv); - } - - /* now interpolate column for new image */ - for (j = 0; j < dst->h; j++) { - float y = (1.00 * j / dst->h) * src->h; - int32_t cvy[4], rv[4], gv[4], bv[4]; - int32_t r, g, b; - int yi = y - 1; - - if (yi < 0) - yi = 0; - - if (yi > (int)src->h - 4) - yi = src->h - 4; - - cvy[0] = cubic(yi - y) * MUL + 0.5; - cvy[1] = cubic(yi + 1 - y) * MUL + 0.5; - cvy[2] = cubic(yi + 2 - y) * MUL + 0.5; - cvy[3] = cubic(yi + 3 - y) * MUL + 0.5; - - rv[0] = col_r[yi]; - rv[1] = col_r[yi + 1]; - rv[2] = col_r[yi + 2]; - rv[3] = col_r[yi + 3]; - - gv[0] = col_g[yi]; - gv[1] = col_g[yi + 1]; - gv[2] = col_g[yi + 2]; - gv[3] = col_g[yi + 3]; - - bv[0] = col_b[yi]; - bv[1] = col_b[yi + 1]; - bv[2] = col_b[yi + 2]; - bv[3] = col_b[yi + 3]; - - MUL_I(rv, cvy); - MUL_I(gv, cvy); - MUL_I(bv, cvy); - - r = (SUM_I(rv) + MUL*MUL/2) / MUL / MUL; - g = (SUM_I(gv) + MUL*MUL/2) / MUL / MUL; - b = (SUM_I(bv) + MUL*MUL/2) / MUL / MUL; - - if (r > 255) - r = 255; - - if (r < 0) - r = 0; - - if (g > 255) - g = 255; - - if (g < 0) - g = 0; - - if (b > 255) - b = 255; - - if (b < 0) - b = 0; - - GP_Pixel pix = GP_Pixel_CREATE_RGB888((uint8_t)r, (uint8_t)g, (uint8_t)b); - GP_PutPixel_Raw_24BPP(dst, i, j, pix); - } - - if (GP_ProgressCallbackReport(callback, i, dst->w, dst->h)) - return 1; - } - - GP_ProgressCallbackDone(callback); - return 0; -} - -int GP_FilterInterpolate_CubicInt2(const GP_Context *src, GP_Context *dst, - GP_ProgressCallback *callback) -{ int32_t col_r[src->w], col_g[src->w], col_b[src->w]; uint32_t i, j; @@ -379,28 +259,36 @@ int GP_FilterInterpolate_CubicInt2(const GP_Context *src, GP_Context *dst, for (i = 0; i < dst->h; i++) { float y = (1.00 * i / dst->h) * src->h; int32_t cvy[4]; - int yi = y - 1; + int yi[4]; - if (yi < 0) - yi = 0; - - if (yi > (int)src->h - 4) - yi = src->h - 4; + yi[0] = floor(y - 1); + yi[1] = y; + yi[2] = y + 1; + yi[3] = y + 2; - cvy[0] = cubic(yi - y) * MUL + 0.5; - cvy[1] = cubic(yi + 1 - y) * MUL + 0.5; - cvy[2] = cubic(yi + 2 - y) * MUL + 0.5; - cvy[3] = cubic(yi + 3 - y) * MUL + 0.5; + cvy[0] = cubic(yi[0] - y) * MUL + 0.5; + cvy[1] = cubic(yi[1] - y) * MUL + 0.5; + cvy[2] = cubic(yi[2] - y) * MUL + 0.5; + cvy[3] = cubic(yi[3] - y) * MUL + 0.5; + if (yi[0] < 0) + yi[0] = 0; + + if (yi[2] >= (int)src->h) + yi[2] = src->h - 1; + + if (yi[3] >= (int)src->h) + yi[3] = src->h - 1; + /* Generate interpolated row */ for (j = 0; j < src->w; j++) { int32_t rv[4], gv[4], bv[4]; GP_Pixel pix[4]; - pix[0] = GP_GetPixel_Raw_24BPP(src, j, yi); - pix[1] = GP_GetPixel_Raw_24BPP(src, j, yi + 1); - pix[2] = GP_GetPixel_Raw_24BPP(src, j, yi + 2); - pix[3] = GP_GetPixel_Raw_24BPP(src, j, yi + 3); + pix[0] = GP_GetPixel_Raw_24BPP(src, j, yi[0]); + pix[1] = GP_GetPixel_Raw_24BPP(src, j, yi[1]); + pix[2] = GP_GetPixel_Raw_24BPP(src, j, yi[2]); + pix[3] = GP_GetPixel_Raw_24BPP(src, j, yi[3]); rv[0] = GP_Pixel_GET_R_RGB888(pix[0]); rv[1] = GP_Pixel_GET_R_RGB888(pix[1]); @@ -416,7 +304,7 @@ int GP_FilterInterpolate_CubicInt2(const GP_Context *src, GP_Context *dst, bv[1] = GP_Pixel_GET_B_RGB888(pix[1]); bv[2] = GP_Pixel_GET_B_RGB888(pix[2]); bv[3] = GP_Pixel_GET_B_RGB888(pix[3]); - + MUL_I(rv, cvy); MUL_I(gv, cvy); MUL_I(bv, cvy); @@ -431,33 +319,41 @@ int GP_FilterInterpolate_CubicInt2(const GP_Context *src, GP_Context *dst, float x = (1.00 * j / dst->w) * src->w; int32_t cvx[4], rv[4], gv[4], bv[4]; int32_t r, g, b; - int xi = x - 1; + int xi[4]; + + xi[0] = floor(x - 1); + xi[1] = x; + xi[2] = x + 1; + xi[3] = x + 2; - if (xi < 0) - xi = 0; - - if (xi > (int)src->w - 4) - xi = src->w - 4; + cvx[0] = cubic(xi[0] - x) * MUL + 0.5; + cvx[1] = cubic(xi[1] - x) * MUL + 0.5; + cvx[2] = cubic(xi[2] - x) * MUL + 0.5; + cvx[3] = cubic(xi[3] - x) * MUL + 0.5; - cvx[0] = cubic(xi - x) * MUL + 0.5; - cvx[1] = cubic(xi + 1 - x) * MUL + 0.5; - cvx[2] = cubic(xi + 2 - x) * MUL + 0.5; - cvx[3] = cubic(xi + 3 - x) * MUL + 0.5; + if (xi[0] < 0) + xi[0] = 0; + + if (xi[2] >= (int)src->w) + xi[2] = src->w - 1; - rv[0] = col_r[xi]; - rv[1] = col_r[xi + 1]; - rv[2] = col_r[xi + 2]; - rv[3] = col_r[xi + 3]; + if (xi[3] >= (int)src->w) + xi[3] = src->w - 1; + + rv[0] = col_r[xi[0]]; + rv[1] = col_r[xi[1]]; + rv[2] = col_r[xi[2]]; + rv[3] = col_r[xi[3]]; - gv[0] = col_g[xi]; - gv[1] = col_g[xi + 1]; - gv[2] = col_g[xi + 2]; - gv[3] = col_g[xi + 3]; + gv[0] = col_g[xi[0]]; + gv[1] = col_g[xi[1]]; + gv[2] = col_g[xi[2]]; + gv[3] = col_g[xi[3]]; - bv[0] = col_b[xi]; - bv[1] = col_b[xi + 1]; - bv[2] = col_b[xi + 2]; - bv[3] = col_b[xi + 3]; + bv[0] = col_b[xi[0]]; + bv[1] = col_b[xi[1]]; + bv[2] = col_b[xi[2]]; + bv[3] = col_b[xi[3]]; MUL_I(rv, cvx); MUL_I(gv, cvx); @@ -467,23 +363,9 @@ int GP_FilterInterpolate_CubicInt2(const GP_Context *src, GP_Context *dst, g = (SUM_I(gv) + MUL*MUL/2) / MUL / MUL; b = (SUM_I(bv) + MUL*MUL/2) / MUL / MUL; - if (r > 255) - r = 255; - - if (r < 0) - r = 0; - - if (g > 255) - g = 255; - - if (g < 0) - g = 0; - - if (b > 255) - b = 255; - - if (b < 0) - b = 0; + CLAMP(r); + CLAMP(g); + CLAMP(b); GP_Pixel pix = GP_Pixel_CREATE_RGB888((uint8_t)r, (uint8_t)g, (uint8_t)b); GP_PutPixel_Raw_24BPP(dst, j, i, pix); @@ -507,7 +389,7 @@ int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, case GP_INTERP_CUBIC: return GP_FilterInterpolate_Cubic(src, dst, callback); case GP_INTERP_CUBIC_INT: - return GP_FilterInterpolate_CubicInt2(src, dst, callback); + return GP_FilterInterpolate_CubicInt(src, dst, callback); } return 1; ----------------------------------------------------------------------- Summary of changes: libs/filters/GP_Linear.c | 2 +- libs/filters/GP_Resize.c | 366 ++++++++++++++++------------------------------ 2 files changed, 125 insertions(+), 243 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: 1d97b3f64cb80fea14245d253f6c7db64b8495b2
by metan 22 Jan '12

22 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via 1d97b3f64cb80fea14245d253f6c7db64b8495b2 (commit) from f6d646faac5ebf7a3f0831af73a79709fae9d36a (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/1d97b3f64cb80fea14245d253f6c7db64b84… commit 1d97b3f64cb80fea14245d253f6c7db64b8495b2 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sun Jan 22 14:27:52 2012 +0100 filters: Bicubic interpolate first on x * That way we are using the caches better. diff --git a/libs/filters/GP_Resize.c b/libs/filters/GP_Resize.c index 05f2c3d..7b26eaa 100644 --- a/libs/filters/GP_Resize.c +++ b/libs/filters/GP_Resize.c @@ -366,6 +366,136 @@ int GP_FilterInterpolate_CubicInt(const GP_Context *src, GP_Context *dst, return 0; } +int GP_FilterInterpolate_CubicInt2(const GP_Context *src, GP_Context *dst, + GP_ProgressCallback *callback) +{ + int32_t col_r[src->w], col_g[src->w], col_b[src->w]; + uint32_t i, j; + + GP_DEBUG(1, "Scaling image %ux%u -> %ux%u %2.2f %2.2f", + src->w, src->h, dst->w, dst->h, + 1.00 * dst->w / src->w, 1.00 * dst->h / src->h); + + for (i = 0; i < dst->h; i++) { + float y = (1.00 * i / dst->h) * src->h; + int32_t cvy[4]; + int yi = y - 1; + + if (yi < 0) + yi = 0; + + if (yi > (int)src->h - 4) + yi = src->h - 4; + + cvy[0] = cubic(yi - y) * MUL + 0.5; + cvy[1] = cubic(yi + 1 - y) * MUL + 0.5; + cvy[2] = cubic(yi + 2 - y) * MUL + 0.5; + cvy[3] = cubic(yi + 3 - y) * MUL + 0.5; + + /* Generate interpolated row */ + for (j = 0; j < src->w; j++) { + int32_t rv[4], gv[4], bv[4]; + GP_Pixel pix[4]; + + pix[0] = GP_GetPixel_Raw_24BPP(src, j, yi); + pix[1] = GP_GetPixel_Raw_24BPP(src, j, yi + 1); + pix[2] = GP_GetPixel_Raw_24BPP(src, j, yi + 2); + pix[3] = GP_GetPixel_Raw_24BPP(src, j, yi + 3); + + rv[0] = GP_Pixel_GET_R_RGB888(pix[0]); + rv[1] = GP_Pixel_GET_R_RGB888(pix[1]); + rv[2] = GP_Pixel_GET_R_RGB888(pix[2]); + rv[3] = GP_Pixel_GET_R_RGB888(pix[3]); + + gv[0] = GP_Pixel_GET_G_RGB888(pix[0]); + gv[1] = GP_Pixel_GET_G_RGB888(pix[1]); + gv[2] = GP_Pixel_GET_G_RGB888(pix[2]); + gv[3] = GP_Pixel_GET_G_RGB888(pix[3]); + + bv[0] = GP_Pixel_GET_B_RGB888(pix[0]); + bv[1] = GP_Pixel_GET_B_RGB888(pix[1]); + bv[2] = GP_Pixel_GET_B_RGB888(pix[2]); + bv[3] = GP_Pixel_GET_B_RGB888(pix[3]); + + MUL_I(rv, cvy); + MUL_I(gv, cvy); + MUL_I(bv, cvy); + + col_r[j] = SUM_I(rv); + col_g[j] = SUM_I(gv); + col_b[j] = SUM_I(bv); + } + + /* now interpolate column for new image */ + for (j = 0; j < dst->w; j++) { + float x = (1.00 * j / dst->w) * src->w; + int32_t cvx[4], rv[4], gv[4], bv[4]; + int32_t r, g, b; + int xi = x - 1; + + if (xi < 0) + xi = 0; + + if (xi > (int)src->w - 4) + xi = src->w - 4; + + cvx[0] = cubic(xi - x) * MUL + 0.5; + cvx[1] = cubic(xi + 1 - x) * MUL + 0.5; + cvx[2] = cubic(xi + 2 - x) * MUL + 0.5; + cvx[3] = cubic(xi + 3 - x) * MUL + 0.5; + + rv[0] = col_r[xi]; + rv[1] = col_r[xi + 1]; + rv[2] = col_r[xi + 2]; + rv[3] = col_r[xi + 3]; + + gv[0] = col_g[xi]; + gv[1] = col_g[xi + 1]; + gv[2] = col_g[xi + 2]; + gv[3] = col_g[xi + 3]; + + bv[0] = col_b[xi]; + bv[1] = col_b[xi + 1]; + bv[2] = col_b[xi + 2]; + bv[3] = col_b[xi + 3]; + + MUL_I(rv, cvx); + MUL_I(gv, cvx); + MUL_I(bv, cvx); + + r = (SUM_I(rv) + MUL*MUL/2) / MUL / MUL; + g = (SUM_I(gv) + MUL*MUL/2) / MUL / MUL; + b = (SUM_I(bv) + MUL*MUL/2) / MUL / MUL; + + if (r > 255) + r = 255; + + if (r < 0) + r = 0; + + if (g > 255) + g = 255; + + if (g < 0) + g = 0; + + if (b > 255) + b = 255; + + if (b < 0) + b = 0; + + GP_Pixel pix = GP_Pixel_CREATE_RGB888((uint8_t)r, (uint8_t)g, (uint8_t)b); + GP_PutPixel_Raw_24BPP(dst, j, i, pix); + } + + if (GP_ProgressCallbackReport(callback, i, dst->h, dst->w)) + return 1; + } + + GP_ProgressCallbackDone(callback); + return 0; +} int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, GP_InterpolationType type, @@ -377,7 +507,7 @@ int GP_FilterResize_Raw(const GP_Context *src, GP_Context *dst, case GP_INTERP_CUBIC: return GP_FilterInterpolate_Cubic(src, dst, callback); case GP_INTERP_CUBIC_INT: - return GP_FilterInterpolate_CubicInt(src, dst, callback); + return GP_FilterInterpolate_CubicInt2(src, dst, callback); } return 1; ----------------------------------------------------------------------- Summary of changes: libs/filters/GP_Resize.c | 132 +++++++++++++++++++++++++++++++++++++++++++++- 1 files changed, 131 insertions(+), 1 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: f6d646faac5ebf7a3f0831af73a79709fae9d36a
by metan 22 Jan '12

22 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via f6d646faac5ebf7a3f0831af73a79709fae9d36a (commit) from 47b078113632ec8cf3e760af0dce2e0bf0c01c92 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/f6d646faac5ebf7a3f0831af73a79709fae9… commit f6d646faac5ebf7a3f0831af73a79709fae9d36a Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sun Jan 22 14:15:49 2012 +0100 filters: Bicubic resampling fix cubic coefs. * The cubic function is symetric, so this doesn't really matter, but this way we have it accordingly to the definition. diff --git a/libs/filters/GP_Resize.c b/libs/filters/GP_Resize.c index ac5ea0d..05f2c3d 100644 --- a/libs/filters/GP_Resize.c +++ b/libs/filters/GP_Resize.c @@ -113,10 +113,10 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, if (xi > (int)src->w - 4) xi = src->w - 4; - cvx.f[0] = cubic(x - xi); - cvx.f[1] = cubic(x - xi - 1); - cvx.f[2] = cubic(x - xi - 2); - cvx.f[3] = cubic(x - xi - 3); + cvx.f[0] = cubic(xi - x); + cvx.f[1] = cubic(xi + 1 - x); + cvx.f[2] = cubic(xi + 2 - x); + cvx.f[3] = cubic(xi + 3 - x); /* Generate interpolated column */ for (j = 0; j < src->h; j++) { @@ -157,7 +157,7 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, float y = (1.00 * j / dst->h) * src->h; v4f cvy, rv, gv, bv; float r, g, b; - int yi = y; + int yi = y - 1; if (yi < 0) yi = 0; @@ -165,10 +165,10 @@ int GP_FilterInterpolate_Cubic(const GP_Context *src, GP_Context *dst, if (yi > (int)src->h - 4) yi = src->h - 4; - cvy.f[0] = cubic(y - yi); - cvy.f[1] = cubic(y - yi - 1); - cvy.f[2] = cubic(y - yi - 2); - cvy.f[3] = cubic(y - yi - 3); + cvy.f[0] = cubic(yi - y); + cvy.f[1] = cubic(yi + 1 - y); + cvy.f[2] = cubic(yi + 2 - y); + cvy.f[3] = cubic(yi + 3 - y); rv.f[0] = col_r[yi]; rv.f[1] = col_r[yi + 1]; @@ -255,11 +255,11 @@ int GP_FilterInterpolate_CubicInt(const GP_Context *src, GP_Context *dst, if (xi > (int)src->w - 4) xi = src->w - 4; - - cvx[0] = cubic(x - xi) * MUL + 0.5; - cvx[1] = cubic(x - xi - 1) * MUL + 0.5; - cvx[2] = cubic(x - xi - 2) * MUL + 0.5; - cvx[3] = cubic(x - xi - 3) * MUL + 0.5; + + cvx[0] = cubic(xi - x) * MUL + 0.5; + cvx[1] = cubic(xi + 1 - x) * MUL + 0.5; + cvx[2] = cubic(xi + 2 - x) * MUL + 0.5; + cvx[3] = cubic(xi + 3 - x) * MUL + 0.5; /* Generate interpolated column */ for (j = 0; j < src->h; j++) { @@ -308,11 +308,11 @@ int GP_FilterInterpolate_CubicInt(const GP_Context *src, GP_Context *dst, if (yi > (int)src->h - 4) yi = src->h - 4; - cvy[0] = cubic(y - yi) * MUL + 0.5; - cvy[1] = cubic(y - yi - 1) * MUL + 0.5; - cvy[2] = cubic(y - yi - 2) * MUL + 0.5; - cvy[3] = cubic(y - yi - 3) * MUL + 0.5; - + cvy[0] = cubic(yi - y) * MUL + 0.5; + cvy[1] = cubic(yi + 1 - y) * MUL + 0.5; + cvy[2] = cubic(yi + 2 - y) * MUL + 0.5; + cvy[3] = cubic(yi + 3 - y) * MUL + 0.5; + rv[0] = col_r[yi]; rv[1] = col_r[yi + 1]; rv[2] = col_r[yi + 2]; ----------------------------------------------------------------------- Summary of changes: libs/filters/GP_Resize.c | 38 +++++++++++++++++++------------------- 1 files changed, 19 insertions(+), 19 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: 47b078113632ec8cf3e760af0dce2e0bf0c01c92
by metan 21 Jan '12

21 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via 47b078113632ec8cf3e760af0dce2e0bf0c01c92 (commit) via 7850ea0241dc2bef8cd4d45d3c22697206788cc9 (commit) from d093866437ca8adc100d6bc2a40011a26c26e813 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/47b078113632ec8cf3e760af0dce2e0bf0c0… commit 47b078113632ec8cf3e760af0dce2e0bf0c01c92 Merge: 7850ea0 d093866 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sat Jan 21 23:17:46 2012 +0100 Merge ssh://repo.or.cz/srv/git/gfxprim http://repo.or.cz/w/gfxprim.git/commit/7850ea0241dc2bef8cd4d45d3c2269720678… commit 7850ea0241dc2bef8cd4d45d3c22697206788cc9 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sat Jan 21 23:17:03 2012 +0100 filters: Some more fixes for blur. diff --git a/libs/filters/GP_Linear.c b/libs/filters/GP_Linear.c index 06145f6..7101472 100644 --- a/libs/filters/GP_Linear.c +++ b/libs/filters/GP_Linear.c @@ -90,7 +90,7 @@ int GP_FilterGaussianBlur_Raw(const GP_Context *src, GP_Context *dst, float kernel_x[size_x]; gaussian_kernel_init(sigma_x, kernel_x); - if (GP_FilterHLinearConvolutionInt_Raw(src, dst, kernel_x, size_x, + if (GP_FilterHLinearConvolution_Raw(src, dst, kernel_x, size_x, new_callback)) return 1; } @@ -103,7 +103,7 @@ int GP_FilterGaussianBlur_Raw(const GP_Context *src, GP_Context *dst, float kernel_y[size_y]; gaussian_kernel_init(sigma_y, kernel_y); - if (GP_FilterVLinearConvolutionInt_Raw(dst, dst, kernel_y, size_y, + if (GP_FilterVLinearConvolution_Raw(dst, dst, kernel_y, size_y, new_callback)) return 1; } @@ -137,213 +137,6 @@ GP_Context *GP_FilterGaussianBlur(const GP_Context *src, GP_Context *dst, return dst; } -int GP_FilterHLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, - float kernel[], uint32_t kw, - GP_ProgressCallback *callback) -{ - float kernel_sum = 0; - GP_Coord x, y; - uint32_t i; - GP_Pixel pix; - - GP_DEBUG(1, "Horizontal linear convolution kernel width %i image %ux%u", - kw, src->w, src->h); - - /* count kernel sum for normalization */ - for (i = 0; i < kw; i++) - kernel_sum += kernel[i]; - - /* do linear convolution */ - for (y = 0; y < (GP_Coord)dst->h; y++) { - uint32_t R[kw], G[kw], B[kw]; - - /* prefill the buffer on the start */ - for (i = 0; i < kw - 1; i++) { - int cx = i - kw/2; - - if (cx < 0) - cx = 0; - - pix = GP_GetPixel_Raw_24BPP(src, cx, y); - - R[i] = GP_Pixel_GET_R_RGB888(pix); - G[i] = GP_Pixel_GET_G_RGB888(pix); - B[i] = GP_Pixel_GET_B_RGB888(pix); - } - - int idx = kw - 1; - - for (x = 0; x < (GP_Coord)dst->w; x++) { - float r = 0, g = 0, b = 0; - - int cx = x + kw/2; - - if (cx >= (int)src->w) - cx = src->w - 1; - - pix = GP_GetPixel_Raw_24BPP(src, cx, y); - - R[idx] = GP_Pixel_GET_R_RGB888(pix); - G[idx] = GP_Pixel_GET_G_RGB888(pix); - B[idx] = GP_Pixel_GET_B_RGB888(pix); - - /* count the pixel value from neighbours weighted by kernel */ - for (i = 0; i < kw; i++) { - int k; - - if ((int)i < idx + 1) - k = kw - idx - 1 + i; - else - k = i - idx - 1; - - r += R[i] * kernel[k]; - g += G[i] * kernel[k]; - b += B[i] * kernel[k]; - } - - /* normalize the result */ - r /= kernel_sum; - g /= kernel_sum; - b /= kernel_sum; - - /* and clamp just to be extra sure */ - if (r > 255) - r = 255; - if (r < 0) - r = 0; - if (g > 255) - g = 255; - if (g < 0) - g = 0; - if (b > 255) - b = 255; - if (b < 0) - b = 0; - - pix = GP_Pixel_CREATE_RGB888((uint32_t)r, (uint32_t)g, (uint32_t)b); - - GP_PutPixel_Raw_24BPP(dst, x, y, pix); - - idx++; - - if (idx >= (int)kw) - idx = 0; - } - - if (GP_ProgressCallbackReport(callback, y, dst->h, dst->w)) - return 1; - } - - GP_ProgressCallbackDone(callback); - return 0; -} - -/* - * Special case for vertical only kernel (10-30% faster). - * - * Can be used in-place. - */ -int GP_FilterVLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, - float kernel[], uint32_t kh, - GP_ProgressCallback *callback) -{ - float kernel_sum = 0; - GP_Coord x, y; - uint32_t i; - - GP_DEBUG(1, "Vertical linear convolution kernel width %i image %ux%u", - kh, src->w, src->h); - - /* count kernel sum for normalization */ - for (i = 0; i < kh; i++) - kernel_sum += kernel[i]; - - /* do linear convolution */ - for (x = 0; x < (GP_Coord)dst->w; x++) { - GP_Pixel pix; - uint32_t R[kh], G[kh], B[kh]; - - /* prefill the buffer on the start */ - for (i = 0; i < kh - 1; i++) { - int cy = i - kh/2; - - if (cy < 0) - cy = 0; - - pix = GP_GetPixel_Raw_24BPP(src, x, cy); - - R[i] = GP_Pixel_GET_R_RGB888(pix); - G[i] = GP_Pixel_GET_G_RGB888(pix); - B[i] = GP_Pixel_GET_B_RGB888(pix); - } - - int idx = kh - 1; - - for (y = 0; y < (GP_Coord)dst->h; y++) { - float r = 0, g = 0, b = 0; - - int cy = y + kh/2; - - if (cy >= (int)src->h) - cy = src->h - 1; - - pix = GP_GetPixel_Raw_24BPP(src, x, cy); - - R[idx] = GP_Pixel_GET_R_RGB888(pix); - G[idx] = GP_Pixel_GET_G_RGB888(pix); - B[idx] = GP_Pixel_GET_B_RGB888(pix); - - /* count the pixel value from neighbours weighted by kernel */ - for (i = 0; i < kh; i++) { - int k; - - if ((int)i < idx + 1) - k = kh - idx - 1 + i; - else - k = i - idx - 1; - - r += R[i] * kernel[k]; - g += G[i] * kernel[k]; - b += B[i] * kernel[k]; - } - - /* normalize the result */ - r /= kernel_sum; - g /= kernel_sum; - b /= kernel_sum; - - /* and clamp just to be extra sure */ - if (r > 255) - r = 255; - if (r < 0) - r = 0; - if (g > 255) - g = 255; - if (g < 0) - g = 0; - if (b > 255) - b = 255; - if (b < 0) - b = 0; - - pix = GP_Pixel_CREATE_RGB888((uint32_t)r, (uint32_t)g, (uint32_t)b); - - GP_PutPixel_Raw_24BPP(dst, x, y, pix); - - idx++; - - if (idx >= (int)kh) - idx = 0; - } - - if (GP_ProgressCallbackReport(callback, x, dst->w, dst->h)) - return 1; - } - - GP_ProgressCallbackDone(callback); - return 0; -} - #define MUL 1024 #define CLAMP(val) do { @@ -353,9 +146,9 @@ int GP_FilterVLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, val = 0; } while (0) -int GP_FilterHLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, - float kernel[], uint32_t kw, - GP_ProgressCallback *callback) +int GP_FilterHLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, + float kernel[], uint32_t kw, + GP_ProgressCallback *callback) { int32_t kernel_sum = 0; GP_Coord x, y; @@ -378,25 +171,30 @@ int GP_FilterHLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, uint8_t R[size], G[size], B[size]; /* Fetch the whole row */ - for (i = 0; i < size; i++) { - GP_Pixel pix; - int cx = i - kw/2; + GP_Pixel pix = GP_GetPixel_Raw_24BPP(src, 0, y); - if (cx < 0) - cx = 0; - - if (cx >= (int)src->w) - cx = src->w - 1; - - pix = GP_GetPixel_Raw_24BPP(src, cx, y); + for (i = 0; i < kw/2; i++) { + R[i] = GP_Pixel_GET_R_RGB888(pix); + G[i] = GP_Pixel_GET_G_RGB888(pix); + B[i] = GP_Pixel_GET_B_RGB888(pix); + } + + for (i = kw/2; i < src->w; i++) { + pix = GP_GetPixel_Raw_24BPP(src, i, y); R[i] = GP_Pixel_GET_R_RGB888(pix); G[i] = GP_Pixel_GET_G_RGB888(pix); B[i] = GP_Pixel_GET_B_RGB888(pix); } + + for (i = src->w; i < size; i++) { + R[i] = GP_Pixel_GET_R_RGB888(pix); + G[i] = GP_Pixel_GET_G_RGB888(pix); + B[i] = GP_Pixel_GET_B_RGB888(pix); + } for (x = 0; x < (GP_Coord)dst->w; x++) { - int32_t r = 0, g = 0, b = 0; + int32_t r = MUL/2, g = MUL/2, b = MUL/2; /* count the pixel value from neighbours weighted by kernel */ for (i = 0; i < kw; i++) { @@ -406,9 +204,9 @@ int GP_FilterHLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, } /* normalize the result */ - r = (r + MUL/2) / kernel_sum; - g = (g + MUL/2) / kernel_sum; - b = (b + MUL/2) / kernel_sum; + r /= kernel_sum; + g /= kernel_sum; + b /= kernel_sum; /* and clamp just to be extra sure */ CLAMP(r); @@ -427,9 +225,9 @@ int GP_FilterHLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, return 0; } -int GP_FilterVLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, - float kernel[], uint32_t kh, - GP_ProgressCallback *callback) +int GP_FilterVLinearConvolution_Raw(const GP_Context *src, GP_Context *dst, + float kernel[], uint32_t kh, + GP_ProgressCallback *callback) { int32_t kernel_sum = 0; GP_Coord x, y; @@ -452,26 +250,30 @@ int GP_FilterVLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, uint8_t R[size], G[size], B[size]; /* Fetch the whole column */ - for (i = 0; i < size; i++) { - GP_Pixel pix; - - int cy = i - kh/2; + GP_Pixel pix = GP_GetPixel_Raw_24BPP(src, x, 0); - if (cy < 0) - cy = 0; + for (i = 0; i < kh/2; i++) { + R[i] = GP_Pixel_GET_R_RGB888(pix); + G[i] = GP_Pixel_GET_G_RGB888(pix); + B[i] = GP_Pixel_GET_B_RGB888(pix); + } - if (cy >= (int)src->h) - cy = src->h - 1; - - pix = GP_GetPixel_Raw_24BPP(src, x, cy); + for (i = kh/2; i < src->h; i++) { + pix = GP_GetPixel_Raw_24BPP(src, x, i); R[i] = GP_Pixel_GET_R_RGB888(pix); G[i] = GP_Pixel_GET_G_RGB888(pix); B[i] = GP_Pixel_GET_B_RGB888(pix); } + + for (i = src->h; i < size; i++) { + R[i] = GP_Pixel_GET_R_RGB888(pix); + G[i] = GP_Pixel_GET_G_RGB888(pix); + B[i] = GP_Pixel_GET_B_RGB888(pix); + } for (y = 0; y < (GP_Coord)dst->h; y++) { - int32_t r = 0, g = 0, b = 0; + int32_t r = MUL/2, g = MUL/2, b = MUL/2; /* count the pixel value from neighbours weighted by kernel */ for (i = 0; i < kh; i++) { @@ -481,9 +283,9 @@ int GP_FilterVLinearConvolutionInt_Raw(const GP_Context *src, GP_Context *dst, } /* normalize the result */ - r = (r + MUL/2) / kernel_sum; - g = (g + MUL/2) / kernel_sum; - b = (b + MUL/2) / kernel_sum; + r /= kernel_sum; + g /= kernel_sum; + b /= kernel_sum; /* and clamp just to be extra sure */ CLAMP(r); ----------------------------------------------------------------------- Summary of changes: libs/filters/GP_Linear.c | 288 +++++++-------------------------------------- 1 files changed, 45 insertions(+), 243 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: d093866437ca8adc100d6bc2a40011a26c26e813
by metan 21 Jan '12

21 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via d093866437ca8adc100d6bc2a40011a26c26e813 (commit) from c5d3e2a62fabfd142f86d47f51e313680d363ea3 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/d093866437ca8adc100d6bc2a40011a26c26… commit d093866437ca8adc100d6bc2a40011a26c26e813 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sat Jan 21 22:23:24 2012 +0100 demos: fbshow add missing -lrt. diff --git a/demos/fbshow/Makefile b/demos/fbshow/Makefile index 74d83fb..c5acfc4 100644 --- a/demos/fbshow/Makefile +++ b/demos/fbshow/Makefile @@ -4,7 +4,7 @@ CSOURCES=$(shell echo *.c) INCLUDE= CFLAGS+=-pthread -LDFLAGS+=-pthread +LDFLAGS+=-pthread -lrt LDLIBS+=-lGP -lGP_backends -lSDL -L$(TOPDIR)/build/ APPS=fbshow ----------------------------------------------------------------------- Summary of changes: demos/fbshow/Makefile | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0

[repo.or.cz] gfxprim.git branch master updated: c5d3e2a62fabfd142f86d47f51e313680d363ea3
by metan 21 Jan '12

21 Jan '12

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project gfxprim.git. The branch, master has been updated via c5d3e2a62fabfd142f86d47f51e313680d363ea3 (commit) from bccb0491d955bfbd1c0cb8d9ca7add8c886c50fe (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- http://repo.or.cz/w/gfxprim.git/commit/c5d3e2a62fabfd142f86d47f51e313680d36… commit c5d3e2a62fabfd142f86d47f51e313680d363ea3 Author: Cyril Hrubis <metan(a)ucw.cz> Date: Sat Jan 21 21:59:13 2012 +0100 demos: fbshow add more timers. diff --git a/demos/fbshow/Makefile b/demos/fbshow/Makefile index 6b04464..74d83fb 100644 --- a/demos/fbshow/Makefile +++ b/demos/fbshow/Makefile @@ -9,5 +9,7 @@ LDLIBS+=-lGP -lGP_backends -lSDL -L$(TOPDIR)/build/ APPS=fbshow +fbshow: cpu_timer.o + include $(TOPDIR)/include.mk include $(TOPDIR)/app.mk diff --git a/demos/fbshow/cpu_timer.c b/demos/fbshow/cpu_timer.c new file mode 100644 index 0000000..88455eb --- /dev/null +++ b/demos/fbshow/cpu_timer.c @@ -0,0 +1,48 @@ +/***************************************************************************** + * This file is part of gfxprim library. * + * * + * Gfxprim is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Lesser General Public * + * License as published by the Free Software Foundation; either * + * version 2.1 of the License, or (at your option) any later version. * + * * + * Gfxprim is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * + * Lesser General Public License for more details. * + * * + * You should have received a copy of the GNU Lesser General Public * + * License along with gfxprim; if not, write to the Free Software * + * Foundation, Inc., 51 Franklin Street, Fifth Floor, * + * Boston, MA 02110-1301 USA * + * * + * Copyright (C) 2009-2012 Cyril Hrubis <metan(a)ucw.cz> * + * * + *****************************************************************************/ + +#include <stdio.h> +#include "cpu_timer.h" + +void cpu_timer_start(struct cpu_timer *self, const char *name) +{ + self->name = name; + clock_gettime(CLOCK_THREAD_CPUTIME_ID, &self->t_start); +} + +void cpu_timer_stop(struct cpu_timer *self) +{ + clock_gettime(CLOCK_THREAD_CPUTIME_ID, &self->t_stop); + + int sec; + int nsec; + + if (self->t_stop.tv_nsec < self->t_start.tv_nsec) { + sec = self->t_stop.tv_sec - self->t_start.tv_sec - 1; + nsec = self->t_stop.tv_nsec + 1000000000 - self->t_start.tv_nsec; + } else { + sec = self->t_stop.tv_sec - self->t_start.tv_sec; + nsec = self->t_stop.tv_nsec - self->t_start.tv_nsec; + } + + printf("TIMER '%s' %i.%09i secn", self->name, sec, nsec); +} diff --git a/demos/fbshow/cpu_timer.h b/demos/fbshow/cpu_timer.h new file mode 100644 index 0000000..1b41144 --- /dev/null +++ b/demos/fbshow/cpu_timer.h @@ -0,0 +1,50 @@ +/***************************************************************************** + * This file is part of gfxprim library. * + * * + * Gfxprim is free software; you can redistribute it and/or * + * modify it under the terms of the GNU Lesser General Public * + * License as published by the Free Software Foundation; either * + * version 2.1 of the License, or (at your option) any later version. * + * * + * Gfxprim is distributed in the hope that it will be useful, * + * but WITHOUT ANY WARRANTY; without even the implied warranty of * + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * + * Lesser General Public License for more details. * + * * + * You should have received a copy of the GNU Lesser General Public * + * License along with gfxprim; if not, write to the Free Software * + * Foundation, Inc., 51 Franklin Street, Fifth Floor, * + * Boston, MA 02110-1301 USA * + * * + * Copyright (C) 2009-2012 Cyril Hrubis <metan(a)ucw.cz> * + * * + *****************************************************************************/ + + /* + + Simple timers to count cpu time. + + */ + +#ifndef __CPU_TIMER_H__ +#define __CPU_TIMER_H__ + +#include <time.h> + +struct cpu_timer { + struct timespec t_start; + struct timespec t_stop; + const char *name; +}; + +/* + * Inialize cpu timer. + */ +void cpu_timer_start(struct cpu_timer *self, const char *name); + +/* + * Stops cpu timer and prints result. + */ +void cpu_timer_stop(struct cpu_timer *self); + +#endif /* __CPU_TIMER_H__ */ diff --git a/demos/fbshow/fbshow.c b/demos/fbshow/fbshow.c index f7ab61f..4d9f6bc 100644 --- a/demos/fbshow/fbshow.c +++ b/demos/fbshow/fbshow.c @@ -35,6 +35,8 @@ #include <backends/GP_Backends.h> #include <input/GP_InputDriverLinux.h> +#include "cpu_timer.h" + static GP_Pixel black_pixel; static GP_Pixel white_pixel; @@ -106,40 +108,14 @@ static const char *img_name(const char *img_path) return NULL; } -#include <time.h> - -static struct timespec t_start; -static struct timespec t_stop; -static const char *t_name; - -static void timer_start(const char *name) -{ - t_name = name; - clock_gettime(CLOCK_THREAD_CPUTIME_ID, &t_start); -} - -static void timer_stop(void) -{ - clock_gettime(CLOCK_THREAD_CPUTIME_ID, &t_stop); - - int sec; - int nsec; - - if (t_stop.tv_nsec < t_start.tv_nsec) { - sec = t_stop.tv_sec - t_start.tv_sec - 1; - nsec = t_stop.tv_nsec + 1000000000 - t_start.tv_nsec; - } else { - sec = t_stop.tv_sec - t_start.tv_sec; - nsec = t_stop.tv_nsec - t_start.tv_nsec; - } - - printf("TIMER '%s' %i.%09i secn", t_name, sec, nsec); -} - static void *image_loader(void *ptr) { struct loader_params *params = ptr; GP_ProgressCallback callback = {.callback = image_loader_callback}; + struct cpu_timer timer; + struct cpu_timer sum_timer; + + cpu_timer_start(&sum_timer, "sum"); show_progress = params->show_progress || params->show_progress_once; params->show_progress_once = 0; @@ -149,7 +125,8 @@ static void *image_loader(void *ptr) GP_Context *img = NULL; callback.priv = "Loading image"; - + + cpu_timer_start(&timer, "Loading"); if (GP_LoadImage(params->img_path, &img, &callback) != 0) { GP_Fill(context, black_pixel); GP_Text(context, NULL, context->w/2, context->h/2, @@ -157,7 +134,8 @@ static void *image_loader(void *ptr) "Failed to load image :("); return NULL; } - + cpu_timer_stop(&timer); + GP_Size w, h; switch (rotate) { @@ -189,16 +167,18 @@ static void *image_loader(void *ptr) GP_Context *ret; if (rat < 1) { - timer_start("blur"); + cpu_timer_start(&timer, "Blur"); callback.priv = "Blurring Image"; if (GP_FilterGaussianBlur(img, img, 0.5/rat, 0.5/rat, &callback) == NULL) return NULL; - timer_stop(); + cpu_timer_stop(&timer); } + cpu_timer_start(&timer, "Resampling"); callback.priv = "Resampling Image"; ret = GP_FilterResize(img, NULL, GP_INTERP_CUBIC_INT, img->w * rat, img->h * rat, &callback); GP_ContextFree(img); + cpu_timer_stop(&timer); if (ret == NULL) return NULL; @@ -240,6 +220,8 @@ static void *image_loader(void *ptr) GP_FillRectXYWH(context, ret->w+cx, 0, cx, context->h, black_pixel); GP_FillRectXYWH(context, 0, ret->h+cy, context->w, cy, black_pixel); + cpu_timer_stop(&sum_timer); + if (!params->show_info) { GP_BackendFlip(backend); return NULL; ----------------------------------------------------------------------- Summary of changes: demos/fbshow/Makefile | 2 + .../GP_Backends.h => demos/fbshow/cpu_timer.c | 37 ++++++++------ .../GP_Backends.h => demos/fbshow/cpu_timer.h | 29 +++++++---- demos/fbshow/fbshow.c | 50 ++++++------------- 4 files changed, 57 insertions(+), 61 deletions(-) copy include/backends/GP_Backends.h => demos/fbshow/cpu_timer.c (71%) copy include/backends/GP_Backends.h => demos/fbshow/cpu_timer.h (80%) repo.or.cz automatic notification. Contact project admin jiri.bluebear.dluhos(a)gmail.com if you want to unsubscribe, or site admin admin(a)repo.or.cz if you receive no reply. -- gfxprim.git ("A simple 2D graphics library with emphasis on correctness and well-defined operation.")

1 0