1*b2055c35SXin Li // Copyright 2010 Google Inc. All Rights Reserved.
2*b2055c35SXin Li //
3*b2055c35SXin Li // Use of this source code is governed by a BSD-style license
4*b2055c35SXin Li // that can be found in the COPYING file in the root of the source
5*b2055c35SXin Li // tree. An additional intellectual property rights grant can be found
6*b2055c35SXin Li // in the file PATENTS. All contributing project authors may
7*b2055c35SXin Li // be found in the AUTHORS file in the root of the source tree.
8*b2055c35SXin Li // -----------------------------------------------------------------------------
9*b2055c35SXin Li //
10*b2055c35SXin Li // YUV->RGB conversion functions
11*b2055c35SXin Li //
12*b2055c35SXin Li // Author: Skal ([email protected])
13*b2055c35SXin Li
14*b2055c35SXin Li #include "src/dsp/yuv.h"
15*b2055c35SXin Li
16*b2055c35SXin Li #include <assert.h>
17*b2055c35SXin Li #include <stdlib.h>
18*b2055c35SXin Li
19*b2055c35SXin Li //-----------------------------------------------------------------------------
20*b2055c35SXin Li // Plain-C version
21*b2055c35SXin Li
22*b2055c35SXin Li #define ROW_FUNC(FUNC_NAME, FUNC, XSTEP) \
23*b2055c35SXin Li static void FUNC_NAME(const uint8_t* y, \
24*b2055c35SXin Li const uint8_t* u, const uint8_t* v, \
25*b2055c35SXin Li uint8_t* dst, int len) { \
26*b2055c35SXin Li const uint8_t* const end = dst + (len & ~1) * (XSTEP); \
27*b2055c35SXin Li while (dst != end) { \
28*b2055c35SXin Li FUNC(y[0], u[0], v[0], dst); \
29*b2055c35SXin Li FUNC(y[1], u[0], v[0], dst + (XSTEP)); \
30*b2055c35SXin Li y += 2; \
31*b2055c35SXin Li ++u; \
32*b2055c35SXin Li ++v; \
33*b2055c35SXin Li dst += 2 * (XSTEP); \
34*b2055c35SXin Li } \
35*b2055c35SXin Li if (len & 1) { \
36*b2055c35SXin Li FUNC(y[0], u[0], v[0], dst); \
37*b2055c35SXin Li } \
38*b2055c35SXin Li } \
39*b2055c35SXin Li
40*b2055c35SXin Li // All variants implemented.
41*b2055c35SXin Li ROW_FUNC(YuvToRgbRow, VP8YuvToRgb, 3)
42*b2055c35SXin Li ROW_FUNC(YuvToBgrRow, VP8YuvToBgr, 3)
43*b2055c35SXin Li ROW_FUNC(YuvToRgbaRow, VP8YuvToRgba, 4)
44*b2055c35SXin Li ROW_FUNC(YuvToBgraRow, VP8YuvToBgra, 4)
45*b2055c35SXin Li ROW_FUNC(YuvToArgbRow, VP8YuvToArgb, 4)
46*b2055c35SXin Li ROW_FUNC(YuvToRgba4444Row, VP8YuvToRgba4444, 2)
47*b2055c35SXin Li ROW_FUNC(YuvToRgb565Row, VP8YuvToRgb565, 2)
48*b2055c35SXin Li
49*b2055c35SXin Li #undef ROW_FUNC
50*b2055c35SXin Li
51*b2055c35SXin Li // Main call for processing a plane with a WebPSamplerRowFunc function:
WebPSamplerProcessPlane(const uint8_t * y,int y_stride,const uint8_t * u,const uint8_t * v,int uv_stride,uint8_t * dst,int dst_stride,int width,int height,WebPSamplerRowFunc func)52*b2055c35SXin Li void WebPSamplerProcessPlane(const uint8_t* y, int y_stride,
53*b2055c35SXin Li const uint8_t* u, const uint8_t* v, int uv_stride,
54*b2055c35SXin Li uint8_t* dst, int dst_stride,
55*b2055c35SXin Li int width, int height, WebPSamplerRowFunc func) {
56*b2055c35SXin Li int j;
57*b2055c35SXin Li for (j = 0; j < height; ++j) {
58*b2055c35SXin Li func(y, u, v, dst, width);
59*b2055c35SXin Li y += y_stride;
60*b2055c35SXin Li if (j & 1) {
61*b2055c35SXin Li u += uv_stride;
62*b2055c35SXin Li v += uv_stride;
63*b2055c35SXin Li }
64*b2055c35SXin Li dst += dst_stride;
65*b2055c35SXin Li }
66*b2055c35SXin Li }
67*b2055c35SXin Li
68*b2055c35SXin Li //-----------------------------------------------------------------------------
69*b2055c35SXin Li // Main call
70*b2055c35SXin Li
71*b2055c35SXin Li WebPSamplerRowFunc WebPSamplers[MODE_LAST];
72*b2055c35SXin Li
73*b2055c35SXin Li extern VP8CPUInfo VP8GetCPUInfo;
74*b2055c35SXin Li extern void WebPInitSamplersSSE2(void);
75*b2055c35SXin Li extern void WebPInitSamplersSSE41(void);
76*b2055c35SXin Li extern void WebPInitSamplersMIPS32(void);
77*b2055c35SXin Li extern void WebPInitSamplersMIPSdspR2(void);
78*b2055c35SXin Li
WEBP_DSP_INIT_FUNC(WebPInitSamplers)79*b2055c35SXin Li WEBP_DSP_INIT_FUNC(WebPInitSamplers) {
80*b2055c35SXin Li WebPSamplers[MODE_RGB] = YuvToRgbRow;
81*b2055c35SXin Li WebPSamplers[MODE_RGBA] = YuvToRgbaRow;
82*b2055c35SXin Li WebPSamplers[MODE_BGR] = YuvToBgrRow;
83*b2055c35SXin Li WebPSamplers[MODE_BGRA] = YuvToBgraRow;
84*b2055c35SXin Li WebPSamplers[MODE_ARGB] = YuvToArgbRow;
85*b2055c35SXin Li WebPSamplers[MODE_RGBA_4444] = YuvToRgba4444Row;
86*b2055c35SXin Li WebPSamplers[MODE_RGB_565] = YuvToRgb565Row;
87*b2055c35SXin Li WebPSamplers[MODE_rgbA] = YuvToRgbaRow;
88*b2055c35SXin Li WebPSamplers[MODE_bgrA] = YuvToBgraRow;
89*b2055c35SXin Li WebPSamplers[MODE_Argb] = YuvToArgbRow;
90*b2055c35SXin Li WebPSamplers[MODE_rgbA_4444] = YuvToRgba4444Row;
91*b2055c35SXin Li
92*b2055c35SXin Li // If defined, use CPUInfo() to overwrite some pointers with faster versions.
93*b2055c35SXin Li if (VP8GetCPUInfo != NULL) {
94*b2055c35SXin Li #if defined(WEBP_HAVE_SSE2)
95*b2055c35SXin Li if (VP8GetCPUInfo(kSSE2)) {
96*b2055c35SXin Li WebPInitSamplersSSE2();
97*b2055c35SXin Li }
98*b2055c35SXin Li #endif // WEBP_HAVE_SSE2
99*b2055c35SXin Li #if defined(WEBP_HAVE_SSE41)
100*b2055c35SXin Li if (VP8GetCPUInfo(kSSE4_1)) {
101*b2055c35SXin Li WebPInitSamplersSSE41();
102*b2055c35SXin Li }
103*b2055c35SXin Li #endif // WEBP_HAVE_SSE41
104*b2055c35SXin Li #if defined(WEBP_USE_MIPS32)
105*b2055c35SXin Li if (VP8GetCPUInfo(kMIPS32)) {
106*b2055c35SXin Li WebPInitSamplersMIPS32();
107*b2055c35SXin Li }
108*b2055c35SXin Li #endif // WEBP_USE_MIPS32
109*b2055c35SXin Li #if defined(WEBP_USE_MIPS_DSP_R2)
110*b2055c35SXin Li if (VP8GetCPUInfo(kMIPSdspR2)) {
111*b2055c35SXin Li WebPInitSamplersMIPSdspR2();
112*b2055c35SXin Li }
113*b2055c35SXin Li #endif // WEBP_USE_MIPS_DSP_R2
114*b2055c35SXin Li }
115*b2055c35SXin Li }
116*b2055c35SXin Li
117*b2055c35SXin Li //-----------------------------------------------------------------------------
118*b2055c35SXin Li // ARGB -> YUV converters
119*b2055c35SXin Li
ConvertARGBToY_C(const uint32_t * argb,uint8_t * y,int width)120*b2055c35SXin Li static void ConvertARGBToY_C(const uint32_t* argb, uint8_t* y, int width) {
121*b2055c35SXin Li int i;
122*b2055c35SXin Li for (i = 0; i < width; ++i) {
123*b2055c35SXin Li const uint32_t p = argb[i];
124*b2055c35SXin Li y[i] = VP8RGBToY((p >> 16) & 0xff, (p >> 8) & 0xff, (p >> 0) & 0xff,
125*b2055c35SXin Li YUV_HALF);
126*b2055c35SXin Li }
127*b2055c35SXin Li }
128*b2055c35SXin Li
WebPConvertARGBToUV_C(const uint32_t * argb,uint8_t * u,uint8_t * v,int src_width,int do_store)129*b2055c35SXin Li void WebPConvertARGBToUV_C(const uint32_t* argb, uint8_t* u, uint8_t* v,
130*b2055c35SXin Li int src_width, int do_store) {
131*b2055c35SXin Li // No rounding. Last pixel is dealt with separately.
132*b2055c35SXin Li const int uv_width = src_width >> 1;
133*b2055c35SXin Li int i;
134*b2055c35SXin Li for (i = 0; i < uv_width; ++i) {
135*b2055c35SXin Li const uint32_t v0 = argb[2 * i + 0];
136*b2055c35SXin Li const uint32_t v1 = argb[2 * i + 1];
137*b2055c35SXin Li // VP8RGBToU/V expects four accumulated pixels. Hence we need to
138*b2055c35SXin Li // scale r/g/b value by a factor 2. We just shift v0/v1 one bit less.
139*b2055c35SXin Li const int r = ((v0 >> 15) & 0x1fe) + ((v1 >> 15) & 0x1fe);
140*b2055c35SXin Li const int g = ((v0 >> 7) & 0x1fe) + ((v1 >> 7) & 0x1fe);
141*b2055c35SXin Li const int b = ((v0 << 1) & 0x1fe) + ((v1 << 1) & 0x1fe);
142*b2055c35SXin Li const int tmp_u = VP8RGBToU(r, g, b, YUV_HALF << 2);
143*b2055c35SXin Li const int tmp_v = VP8RGBToV(r, g, b, YUV_HALF << 2);
144*b2055c35SXin Li if (do_store) {
145*b2055c35SXin Li u[i] = tmp_u;
146*b2055c35SXin Li v[i] = tmp_v;
147*b2055c35SXin Li } else {
148*b2055c35SXin Li // Approximated average-of-four. But it's an acceptable diff.
149*b2055c35SXin Li u[i] = (u[i] + tmp_u + 1) >> 1;
150*b2055c35SXin Li v[i] = (v[i] + tmp_v + 1) >> 1;
151*b2055c35SXin Li }
152*b2055c35SXin Li }
153*b2055c35SXin Li if (src_width & 1) { // last pixel
154*b2055c35SXin Li const uint32_t v0 = argb[2 * i + 0];
155*b2055c35SXin Li const int r = (v0 >> 14) & 0x3fc;
156*b2055c35SXin Li const int g = (v0 >> 6) & 0x3fc;
157*b2055c35SXin Li const int b = (v0 << 2) & 0x3fc;
158*b2055c35SXin Li const int tmp_u = VP8RGBToU(r, g, b, YUV_HALF << 2);
159*b2055c35SXin Li const int tmp_v = VP8RGBToV(r, g, b, YUV_HALF << 2);
160*b2055c35SXin Li if (do_store) {
161*b2055c35SXin Li u[i] = tmp_u;
162*b2055c35SXin Li v[i] = tmp_v;
163*b2055c35SXin Li } else {
164*b2055c35SXin Li u[i] = (u[i] + tmp_u + 1) >> 1;
165*b2055c35SXin Li v[i] = (v[i] + tmp_v + 1) >> 1;
166*b2055c35SXin Li }
167*b2055c35SXin Li }
168*b2055c35SXin Li }
169*b2055c35SXin Li
170*b2055c35SXin Li //-----------------------------------------------------------------------------
171*b2055c35SXin Li
ConvertRGB24ToY_C(const uint8_t * rgb,uint8_t * y,int width)172*b2055c35SXin Li static void ConvertRGB24ToY_C(const uint8_t* rgb, uint8_t* y, int width) {
173*b2055c35SXin Li int i;
174*b2055c35SXin Li for (i = 0; i < width; ++i, rgb += 3) {
175*b2055c35SXin Li y[i] = VP8RGBToY(rgb[0], rgb[1], rgb[2], YUV_HALF);
176*b2055c35SXin Li }
177*b2055c35SXin Li }
178*b2055c35SXin Li
ConvertBGR24ToY_C(const uint8_t * bgr,uint8_t * y,int width)179*b2055c35SXin Li static void ConvertBGR24ToY_C(const uint8_t* bgr, uint8_t* y, int width) {
180*b2055c35SXin Li int i;
181*b2055c35SXin Li for (i = 0; i < width; ++i, bgr += 3) {
182*b2055c35SXin Li y[i] = VP8RGBToY(bgr[2], bgr[1], bgr[0], YUV_HALF);
183*b2055c35SXin Li }
184*b2055c35SXin Li }
185*b2055c35SXin Li
WebPConvertRGBA32ToUV_C(const uint16_t * rgb,uint8_t * u,uint8_t * v,int width)186*b2055c35SXin Li void WebPConvertRGBA32ToUV_C(const uint16_t* rgb,
187*b2055c35SXin Li uint8_t* u, uint8_t* v, int width) {
188*b2055c35SXin Li int i;
189*b2055c35SXin Li for (i = 0; i < width; i += 1, rgb += 4) {
190*b2055c35SXin Li const int r = rgb[0], g = rgb[1], b = rgb[2];
191*b2055c35SXin Li u[i] = VP8RGBToU(r, g, b, YUV_HALF << 2);
192*b2055c35SXin Li v[i] = VP8RGBToV(r, g, b, YUV_HALF << 2);
193*b2055c35SXin Li }
194*b2055c35SXin Li }
195*b2055c35SXin Li
196*b2055c35SXin Li //-----------------------------------------------------------------------------
197*b2055c35SXin Li
198*b2055c35SXin Li void (*WebPConvertRGB24ToY)(const uint8_t* rgb, uint8_t* y, int width);
199*b2055c35SXin Li void (*WebPConvertBGR24ToY)(const uint8_t* bgr, uint8_t* y, int width);
200*b2055c35SXin Li void (*WebPConvertRGBA32ToUV)(const uint16_t* rgb,
201*b2055c35SXin Li uint8_t* u, uint8_t* v, int width);
202*b2055c35SXin Li
203*b2055c35SXin Li void (*WebPConvertARGBToY)(const uint32_t* argb, uint8_t* y, int width);
204*b2055c35SXin Li void (*WebPConvertARGBToUV)(const uint32_t* argb, uint8_t* u, uint8_t* v,
205*b2055c35SXin Li int src_width, int do_store);
206*b2055c35SXin Li
207*b2055c35SXin Li extern void WebPInitConvertARGBToYUVSSE2(void);
208*b2055c35SXin Li extern void WebPInitConvertARGBToYUVSSE41(void);
209*b2055c35SXin Li extern void WebPInitConvertARGBToYUVNEON(void);
210*b2055c35SXin Li
WEBP_DSP_INIT_FUNC(WebPInitConvertARGBToYUV)211*b2055c35SXin Li WEBP_DSP_INIT_FUNC(WebPInitConvertARGBToYUV) {
212*b2055c35SXin Li WebPConvertARGBToY = ConvertARGBToY_C;
213*b2055c35SXin Li WebPConvertARGBToUV = WebPConvertARGBToUV_C;
214*b2055c35SXin Li
215*b2055c35SXin Li WebPConvertRGB24ToY = ConvertRGB24ToY_C;
216*b2055c35SXin Li WebPConvertBGR24ToY = ConvertBGR24ToY_C;
217*b2055c35SXin Li
218*b2055c35SXin Li WebPConvertRGBA32ToUV = WebPConvertRGBA32ToUV_C;
219*b2055c35SXin Li
220*b2055c35SXin Li if (VP8GetCPUInfo != NULL) {
221*b2055c35SXin Li #if defined(WEBP_HAVE_SSE2)
222*b2055c35SXin Li if (VP8GetCPUInfo(kSSE2)) {
223*b2055c35SXin Li WebPInitConvertARGBToYUVSSE2();
224*b2055c35SXin Li }
225*b2055c35SXin Li #endif // WEBP_HAVE_SSE2
226*b2055c35SXin Li #if defined(WEBP_HAVE_SSE41)
227*b2055c35SXin Li if (VP8GetCPUInfo(kSSE4_1)) {
228*b2055c35SXin Li WebPInitConvertARGBToYUVSSE41();
229*b2055c35SXin Li }
230*b2055c35SXin Li #endif // WEBP_HAVE_SSE41
231*b2055c35SXin Li }
232*b2055c35SXin Li
233*b2055c35SXin Li #if defined(WEBP_HAVE_NEON)
234*b2055c35SXin Li if (WEBP_NEON_OMIT_C_CODE ||
235*b2055c35SXin Li (VP8GetCPUInfo != NULL && VP8GetCPUInfo(kNEON))) {
236*b2055c35SXin Li WebPInitConvertARGBToYUVNEON();
237*b2055c35SXin Li }
238*b2055c35SXin Li #endif // WEBP_HAVE_NEON
239*b2055c35SXin Li
240*b2055c35SXin Li assert(WebPConvertARGBToY != NULL);
241*b2055c35SXin Li assert(WebPConvertARGBToUV != NULL);
242*b2055c35SXin Li assert(WebPConvertRGB24ToY != NULL);
243*b2055c35SXin Li assert(WebPConvertBGR24ToY != NULL);
244*b2055c35SXin Li assert(WebPConvertRGBA32ToUV != NULL);
245*b2055c35SXin Li }
246