1 // Copyright 2014 Google Inc. All Rights Reserved. 2 // 3 // Use of this source code is governed by a BSD-style license 4 // that can be found in the COPYING file in the root of the source 5 // tree. An additional intellectual property rights grant can be found 6 // in the file PATENTS. All contributing project authors may 7 // be found in the AUTHORS file in the root of the source tree. 8 // ----------------------------------------------------------------------------- 9 // 10 // WebPPicture utils for colorspace conversion 11 // 12 // Author: Skal (pascal.massimino (at) gmail.com) 13 14 #include <assert.h> 15 #include <stdlib.h> 16 #include <math.h> 17 18 #include "./vp8enci.h" 19 #include "../utils/random.h" 20 #include "../dsp/yuv.h" 21 22 // Uncomment to disable gamma-compression during RGB->U/V averaging 23 #define USE_GAMMA_COMPRESSION 24 25 static const union { 26 uint32_t argb; 27 uint8_t bytes[4]; 28 } test_endian = { 0xff000000u }; 29 #define ALPHA_IS_LAST (test_endian.bytes[3] == 0xff) 30 31 static WEBP_INLINE uint32_t MakeARGB32(int a, int r, int g, int b) { 32 return (((uint32_t)a << 24) | (r << 16) | (g << 8) | b); 33 } 34 35 //------------------------------------------------------------------------------ 36 // Detection of non-trivial transparency 37 38 // Returns true if alpha[] has non-0xff values. 39 static int CheckNonOpaque(const uint8_t* alpha, int width, int height, 40 int x_step, int y_step) { 41 if (alpha == NULL) return 0; 42 while (height-- > 0) { 43 int x; 44 for (x = 0; x < width * x_step; x += x_step) { 45 if (alpha[x] != 0xff) return 1; // TODO(skal): check 4/8 bytes at a time. 46 } 47 alpha += y_step; 48 } 49 return 0; 50 } 51 52 // Checking for the presence of non-opaque alpha. 53 int WebPPictureHasTransparency(const WebPPicture* picture) { 54 if (picture == NULL) return 0; 55 if (!picture->use_argb) { 56 return CheckNonOpaque(picture->a, picture->width, picture->height, 57 1, picture->a_stride); 58 } else { 59 int x, y; 60 const uint32_t* argb = picture->argb; 61 if (argb == NULL) return 0; 62 for (y = 0; y < picture->height; ++y) { 63 for (x = 0; x < picture->width; ++x) { 64 if (argb[x] < 0xff000000u) return 1; // test any alpha values != 0xff 65 } 66 argb += picture->argb_stride; 67 } 68 } 69 return 0; 70 } 71 72 //------------------------------------------------------------------------------ 73 // RGB -> YUV conversion 74 75 static int RGBToY(int r, int g, int b, VP8Random* const rg) { 76 return VP8RGBToY(r, g, b, VP8RandomBits(rg, YUV_FIX)); 77 } 78 79 static int RGBToU(int r, int g, int b, VP8Random* const rg) { 80 return VP8RGBToU(r, g, b, VP8RandomBits(rg, YUV_FIX + 2)); 81 } 82 83 static int RGBToV(int r, int g, int b, VP8Random* const rg) { 84 return VP8RGBToV(r, g, b, VP8RandomBits(rg, YUV_FIX + 2)); 85 } 86 87 //------------------------------------------------------------------------------ 88 89 #if defined(USE_GAMMA_COMPRESSION) 90 91 // gamma-compensates loss of resolution during chroma subsampling 92 #define kGamma 0.80 93 #define kGammaFix 12 // fixed-point precision for linear values 94 #define kGammaScale ((1 << kGammaFix) - 1) 95 #define kGammaTabFix 7 // fixed-point fractional bits precision 96 #define kGammaTabScale (1 << kGammaTabFix) 97 #define kGammaTabRounder (kGammaTabScale >> 1) 98 #define kGammaTabSize (1 << (kGammaFix - kGammaTabFix)) 99 100 static int kLinearToGammaTab[kGammaTabSize + 1]; 101 static uint16_t kGammaToLinearTab[256]; 102 static int kGammaTablesOk = 0; 103 104 static void InitGammaTables(void) { 105 if (!kGammaTablesOk) { 106 int v; 107 const double scale = 1. / kGammaScale; 108 for (v = 0; v <= 255; ++v) { 109 kGammaToLinearTab[v] = 110 (uint16_t)(pow(v / 255., kGamma) * kGammaScale + .5); 111 } 112 for (v = 0; v <= kGammaTabSize; ++v) { 113 const double x = scale * (v << kGammaTabFix); 114 kLinearToGammaTab[v] = (int)(pow(x, 1. / kGamma) * 255. + .5); 115 } 116 kGammaTablesOk = 1; 117 } 118 } 119 120 static WEBP_INLINE uint32_t GammaToLinear(uint8_t v) { 121 return kGammaToLinearTab[v]; 122 } 123 124 // Convert a linear value 'v' to YUV_FIX+2 fixed-point precision 125 // U/V value, suitable for RGBToU/V calls. 126 static WEBP_INLINE int LinearToGamma(uint32_t base_value, int shift) { 127 const int v = base_value << shift; // final uplifted value 128 const int tab_pos = v >> (kGammaTabFix + 2); // integer part 129 const int x = v & ((kGammaTabScale << 2) - 1); // fractional part 130 const int v0 = kLinearToGammaTab[tab_pos]; 131 const int v1 = kLinearToGammaTab[tab_pos + 1]; 132 const int y = v1 * x + v0 * ((kGammaTabScale << 2) - x); // interpolate 133 return (y + kGammaTabRounder) >> kGammaTabFix; // descale 134 } 135 136 #else 137 138 static void InitGammaTables(void) {} 139 static WEBP_INLINE uint32_t GammaToLinear(uint8_t v) { return v; } 140 static WEBP_INLINE int LinearToGamma(uint32_t base_value, int shift) { 141 return (int)(base_value << shift); 142 } 143 144 #endif // USE_GAMMA_COMPRESSION 145 146 //------------------------------------------------------------------------------ 147 148 #define SUM4(ptr) LinearToGamma( \ 149 GammaToLinear((ptr)[0]) + \ 150 GammaToLinear((ptr)[step]) + \ 151 GammaToLinear((ptr)[rgb_stride]) + \ 152 GammaToLinear((ptr)[rgb_stride + step]), 0) \ 153 154 #define SUM2H(ptr) \ 155 LinearToGamma(GammaToLinear((ptr)[0]) + GammaToLinear((ptr)[step]), 1) 156 #define SUM2V(ptr) \ 157 LinearToGamma(GammaToLinear((ptr)[0]) + GammaToLinear((ptr)[rgb_stride]), 1) 158 #define SUM1(ptr) \ 159 LinearToGamma(GammaToLinear((ptr)[0]), 2) 160 161 #define RGB_TO_UV(x, y, SUM) { \ 162 const int src = (2 * (step * (x) + (y) * rgb_stride)); \ 163 const int dst = (x) + (y) * picture->uv_stride; \ 164 const int r = SUM(r_ptr + src); \ 165 const int g = SUM(g_ptr + src); \ 166 const int b = SUM(b_ptr + src); \ 167 picture->u[dst] = RGBToU(r, g, b, &rg); \ 168 picture->v[dst] = RGBToV(r, g, b, &rg); \ 169 } 170 171 static int ImportYUVAFromRGBA(const uint8_t* const r_ptr, 172 const uint8_t* const g_ptr, 173 const uint8_t* const b_ptr, 174 const uint8_t* const a_ptr, 175 int step, // bytes per pixel 176 int rgb_stride, // bytes per scanline 177 float dithering, 178 WebPPicture* const picture) { 179 int x, y; 180 const int width = picture->width; 181 const int height = picture->height; 182 const int has_alpha = CheckNonOpaque(a_ptr, width, height, step, rgb_stride); 183 VP8Random rg; 184 185 if (has_alpha) { 186 picture->colorspace |= WEBP_CSP_ALPHA_BIT; 187 } else { 188 picture->colorspace &= WEBP_CSP_UV_MASK; 189 } 190 picture->use_argb = 0; 191 192 if (!WebPPictureAllocYUVA(picture, width, height)) return 0; 193 194 VP8InitRandom(&rg, dithering); 195 InitGammaTables(); 196 197 // Import luma plane 198 for (y = 0; y < height; ++y) { 199 uint8_t* const dst = &picture->y[y * picture->y_stride]; 200 for (x = 0; x < width; ++x) { 201 const int offset = step * x + y * rgb_stride; 202 dst[x] = RGBToY(r_ptr[offset], g_ptr[offset], b_ptr[offset], &rg); 203 } 204 } 205 206 // Downsample U/V plane 207 for (y = 0; y < (height >> 1); ++y) { 208 for (x = 0; x < (width >> 1); ++x) { 209 RGB_TO_UV(x, y, SUM4); 210 } 211 if (width & 1) { 212 RGB_TO_UV(x, y, SUM2V); 213 } 214 } 215 if (height & 1) { 216 for (x = 0; x < (width >> 1); ++x) { 217 RGB_TO_UV(x, y, SUM2H); 218 } 219 if (width & 1) { 220 RGB_TO_UV(x, y, SUM1); 221 } 222 } 223 224 if (has_alpha) { 225 assert(step >= 4); 226 assert(picture->a != NULL); 227 for (y = 0; y < height; ++y) { 228 for (x = 0; x < width; ++x) { 229 picture->a[x + y * picture->a_stride] = 230 a_ptr[step * x + y * rgb_stride]; 231 } 232 } 233 } 234 return 1; 235 } 236 237 #undef SUM4 238 #undef SUM2V 239 #undef SUM2H 240 #undef SUM1 241 #undef RGB_TO_UV 242 243 //------------------------------------------------------------------------------ 244 // call for ARGB->YUVA conversion 245 246 int WebPPictureARGBToYUVADithered(WebPPicture* picture, WebPEncCSP colorspace, 247 float dithering) { 248 if (picture == NULL) return 0; 249 if (picture->argb == NULL) { 250 return WebPEncodingSetError(picture, VP8_ENC_ERROR_NULL_PARAMETER); 251 } else { 252 const uint8_t* const argb = (const uint8_t*)picture->argb; 253 const uint8_t* const r = ALPHA_IS_LAST ? argb + 2 : argb + 1; 254 const uint8_t* const g = ALPHA_IS_LAST ? argb + 1 : argb + 2; 255 const uint8_t* const b = ALPHA_IS_LAST ? argb + 0 : argb + 3; 256 const uint8_t* const a = ALPHA_IS_LAST ? argb + 3 : argb + 0; 257 258 picture->colorspace = colorspace; 259 return ImportYUVAFromRGBA(r, g, b, a, 4, 4 * picture->argb_stride, 260 dithering, picture); 261 } 262 } 263 264 int WebPPictureARGBToYUVA(WebPPicture* picture, WebPEncCSP colorspace) { 265 return WebPPictureARGBToYUVADithered(picture, colorspace, 0.f); 266 } 267 268 //------------------------------------------------------------------------------ 269 // call for YUVA -> ARGB conversion 270 271 int WebPPictureYUVAToARGB(WebPPicture* picture) { 272 if (picture == NULL) return 0; 273 if (picture->y == NULL || picture->u == NULL || picture->v == NULL) { 274 return WebPEncodingSetError(picture, VP8_ENC_ERROR_NULL_PARAMETER); 275 } 276 if ((picture->colorspace & WEBP_CSP_ALPHA_BIT) && picture->a == NULL) { 277 return WebPEncodingSetError(picture, VP8_ENC_ERROR_NULL_PARAMETER); 278 } 279 if ((picture->colorspace & WEBP_CSP_UV_MASK) != WEBP_YUV420) { 280 return WebPEncodingSetError(picture, VP8_ENC_ERROR_INVALID_CONFIGURATION); 281 } 282 // Allocate a new argb buffer (discarding the previous one). 283 if (!WebPPictureAllocARGB(picture, picture->width, picture->height)) return 0; 284 picture->use_argb = 1; 285 286 // Convert 287 { 288 int y; 289 const int width = picture->width; 290 const int height = picture->height; 291 const int argb_stride = 4 * picture->argb_stride; 292 uint8_t* dst = (uint8_t*)picture->argb; 293 const uint8_t *cur_u = picture->u, *cur_v = picture->v, *cur_y = picture->y; 294 WebPUpsampleLinePairFunc upsample = WebPGetLinePairConverter(ALPHA_IS_LAST); 295 296 // First row, with replicated top samples. 297 upsample(cur_y, NULL, cur_u, cur_v, cur_u, cur_v, dst, NULL, width); 298 cur_y += picture->y_stride; 299 dst += argb_stride; 300 // Center rows. 301 for (y = 1; y + 1 < height; y += 2) { 302 const uint8_t* const top_u = cur_u; 303 const uint8_t* const top_v = cur_v; 304 cur_u += picture->uv_stride; 305 cur_v += picture->uv_stride; 306 upsample(cur_y, cur_y + picture->y_stride, top_u, top_v, cur_u, cur_v, 307 dst, dst + argb_stride, width); 308 cur_y += 2 * picture->y_stride; 309 dst += 2 * argb_stride; 310 } 311 // Last row (if needed), with replicated bottom samples. 312 if (height > 1 && !(height & 1)) { 313 upsample(cur_y, NULL, cur_u, cur_v, cur_u, cur_v, dst, NULL, width); 314 } 315 // Insert alpha values if needed, in replacement for the default 0xff ones. 316 if (picture->colorspace & WEBP_CSP_ALPHA_BIT) { 317 for (y = 0; y < height; ++y) { 318 uint32_t* const argb_dst = picture->argb + y * picture->argb_stride; 319 const uint8_t* const src = picture->a + y * picture->a_stride; 320 int x; 321 for (x = 0; x < width; ++x) { 322 argb_dst[x] = (argb_dst[x] & 0x00ffffffu) | ((uint32_t)src[x] << 24); 323 } 324 } 325 } 326 } 327 return 1; 328 } 329 330 //------------------------------------------------------------------------------ 331 // automatic import / conversion 332 333 static int Import(WebPPicture* const picture, 334 const uint8_t* const rgb, int rgb_stride, 335 int step, int swap_rb, int import_alpha) { 336 int y; 337 const uint8_t* const r_ptr = rgb + (swap_rb ? 2 : 0); 338 const uint8_t* const g_ptr = rgb + 1; 339 const uint8_t* const b_ptr = rgb + (swap_rb ? 0 : 2); 340 const uint8_t* const a_ptr = import_alpha ? rgb + 3 : NULL; 341 const int width = picture->width; 342 const int height = picture->height; 343 344 if (!picture->use_argb) { 345 return ImportYUVAFromRGBA(r_ptr, g_ptr, b_ptr, a_ptr, step, rgb_stride, 346 0.f /* no dithering */, picture); 347 } 348 if (!WebPPictureAlloc(picture)) return 0; 349 350 assert(step >= (import_alpha ? 4 : 3)); 351 for (y = 0; y < height; ++y) { 352 uint32_t* const dst = &picture->argb[y * picture->argb_stride]; 353 int x; 354 for (x = 0; x < width; ++x) { 355 const int offset = step * x + y * rgb_stride; 356 dst[x] = MakeARGB32(import_alpha ? a_ptr[offset] : 0xff, 357 r_ptr[offset], g_ptr[offset], b_ptr[offset]); 358 } 359 } 360 return 1; 361 } 362 363 // Public API 364 365 int WebPPictureImportRGB(WebPPicture* picture, 366 const uint8_t* rgb, int rgb_stride) { 367 return (picture != NULL) ? Import(picture, rgb, rgb_stride, 3, 0, 0) : 0; 368 } 369 370 int WebPPictureImportBGR(WebPPicture* picture, 371 const uint8_t* rgb, int rgb_stride) { 372 return (picture != NULL) ? Import(picture, rgb, rgb_stride, 3, 1, 0) : 0; 373 } 374 375 int WebPPictureImportRGBA(WebPPicture* picture, 376 const uint8_t* rgba, int rgba_stride) { 377 return (picture != NULL) ? Import(picture, rgba, rgba_stride, 4, 0, 1) : 0; 378 } 379 380 int WebPPictureImportBGRA(WebPPicture* picture, 381 const uint8_t* rgba, int rgba_stride) { 382 return (picture != NULL) ? Import(picture, rgba, rgba_stride, 4, 1, 1) : 0; 383 } 384 385 int WebPPictureImportRGBX(WebPPicture* picture, 386 const uint8_t* rgba, int rgba_stride) { 387 return (picture != NULL) ? Import(picture, rgba, rgba_stride, 4, 0, 0) : 0; 388 } 389 390 int WebPPictureImportBGRX(WebPPicture* picture, 391 const uint8_t* rgba, int rgba_stride) { 392 return (picture != NULL) ? Import(picture, rgba, rgba_stride, 4, 1, 0) : 0; 393 } 394 395 //------------------------------------------------------------------------------ 396