1 /**************************************************************************
3 * Copyright 2009 Younes Manton.
6 * Permission is hereby granted, free of charge, to any person obtaining a
7 * copy of this software and associated documentation files (the
8 * "Software"), to deal in the Software without restriction, including
9 * without limitation the rights to use, copy, modify, merge, publish,
10 * distribute, sub license, and/or sell copies of the Software, and to
11 * permit persons to whom the Software is furnished to do so, subject to
12 * the following conditions:
14 * The above copyright notice and this permission notice (including the
15 * next paragraph) shall be included in all copies or substantial portions
18 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
19 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
20 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NON-INFRINGEMENT.
21 * IN NO EVENT SHALL TUNGSTEN GRAPHICS AND/OR ITS SUPPLIERS BE LIABLE FOR
22 * ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
23 * TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
24 * SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
26 **************************************************************************/
28 #include "util/u_math.h"
29 #include "util/u_debug.h"
34 * Color space conversion formulas
36 * To convert YCbCr to RGB,
41 * To calculate the color space conversion matrix csc with ProcAmp adjustments,
42 * mat44 csc, cstd, procamp, bias
43 * csc = cstd * (procamp * bias)
45 * Where cstd is a matrix corresponding to one of the color standards (BT.601, BT.709, etc)
46 * adjusted for the kind of YCbCr -> RGB mapping wanted (1:1, full),
47 * bias is a matrix corresponding to the kind of YCbCr -> RGB mapping wanted (1:1, full)
49 * To calculate procamp,
50 * mat44 procamp, hue, saturation, brightness, contrast
51 * procamp = brightness * (saturation * (contrast * hue))
53 * procamp = saturation * (brightness * (contrast * hue))
75 * [ 0, cos(h), sin(h), 0]
76 * [ 0, -sin(h), cos(h), 0]
81 * [ 0, c*s*cos(h), c*s*sin(h), 0]
82 * [ 0, -c*s*sin(h), c*s*cos(h), 0]
92 * [ c*cstd[ 0], c*cstd[ 1]*s*cos(h) - c*cstd[ 2]*s*sin(h), c*cstd[ 2]*s*cos(h) + c*cstd[ 1]*s*sin(h), cstd[ 3] + cstd[ 0]*(b + c*ybias) + cstd[ 1]*(c*cbbias*s*cos(h) + c*crbias*s*sin(h)) + cstd[ 2]*(c*crbias*s*cos(h) - c*cbbias*s*sin(h))]
93 * [ c*cstd[ 4], c*cstd[ 5]*s*cos(h) - c*cstd[ 6]*s*sin(h), c*cstd[ 6]*s*cos(h) + c*cstd[ 5]*s*sin(h), cstd[ 7] + cstd[ 4]*(b + c*ybias) + cstd[ 5]*(c*cbbias*s*cos(h) + c*crbias*s*sin(h)) + cstd[ 6]*(c*crbias*s*cos(h) - c*cbbias*s*sin(h))]
94 * [ c*cstd[ 8], c*cstd[ 9]*s*cos(h) - c*cstd[10]*s*sin(h), c*cstd[10]*s*cos(h) + c*cstd[ 9]*s*sin(h), cstd[11] + cstd[ 8]*(b + c*ybias) + cstd[ 9]*(c*cbbias*s*cos(h) + c*crbias*s*sin(h)) + cstd[10]*(c*crbias*s*cos(h) - c*cbbias*s*sin(h))]
95 * [ c*cstd[12], c*cstd[13]*s*cos(h) - c*cstd[14]*s*sin(h), c*cstd[14]*s*cos(h) + c*cstd[13]*s*sin(h), cstd[15] + cstd[12]*(b + c*ybias) + cstd[13]*(c*cbbias*s*cos(h) + c*crbias*s*sin(h)) + cstd[14]*(c*crbias*s*cos(h) - c*cbbias*s*sin(h))]
99 * Converts ITU-R BT.601 YCbCr pixels to RGB pixels where:
100 * Y is in [16,235], Cb and Cr are in [16,240]
101 * R, G, and B are in [16,235]
103 static const float bt_601
[16] =
105 1.0f
, 0.0f
, 1.371f
, 0.0f
,
106 1.0f
, -0.336f
, -0.698f
, 0.0f
,
107 1.0f
, 1.732f
, 0.0f
, 0.0f
,
108 0.0f
, 0.0f
, 0.0f
, 1.0f
112 * Converts ITU-R BT.601 YCbCr pixels to RGB pixels where:
113 * Y is in [16,235], Cb and Cr are in [16,240]
114 * R, G, and B are in [0,255]
116 static const float bt_601_full
[16] =
118 1.164f
, 0.0f
, 1.596f
, 0.0f
,
119 1.164f
, -0.391f
, -0.813f
, 0.0f
,
120 1.164f
, 2.018f
, 0.0f
, 0.0f
,
121 0.0f
, 0.0f
, 0.0f
, 1.0f
125 * Converts ITU-R BT.709 YCbCr pixels to RGB pixels where:
126 * Y is in [16,235], Cb and Cr are in [16,240]
127 * R, G, and B are in [16,235]
129 static const float bt_709
[16] =
131 1.0f
, 0.0f
, 1.540f
, 0.0f
,
132 1.0f
, -0.183f
, -0.459f
, 0.0f
,
133 1.0f
, 1.816f
, 0.0f
, 0.0f
,
134 0.0f
, 0.0f
, 0.0f
, 1.0f
138 * Converts ITU-R BT.709 YCbCr pixels to RGB pixels where:
139 * Y is in [16,235], Cb and Cr are in [16,240]
140 * R, G, and B are in [0,255]
142 static const float bt_709_full
[16] =
144 1.164f
, 0.0f
, 1.793f
, 0.0f
,
145 1.164f
, -0.213f
, -0.534f
, 0.0f
,
146 1.164f
, 2.115f
, 0.0f
, 0.0f
,
147 0.0f
, 0.0f
, 0.0f
, 1.0f
150 static const float smpte240m
[16] =
152 1.0f
, 0.0f
, 1.582f
, 0.0f
,
153 1.0f
, -0.228f
, -0.478f
, 0.0f
,
154 1.0f
, 1.833f
, 0.0f
, 0.0f
,
155 0.0f
, 0.0f
, 0.0f
, 1.0f
158 static const float smpte240m_full
[16] =
160 1.164f
, 0.0f
, 1.794f
, 0.0f
,
161 1.164f
, -0.258f
, -0.543f
, 0.0f
,
162 1.164f
, 2.079f
, 0.0f
, 0.0f
,
163 0.0f
, 0.0f
, 0.0f
, 1.0f
166 static const float identity
[16] =
168 1.0f
, 0.0f
, 0.0f
, 0.0f
,
169 0.0f
, 1.0f
, 0.0f
, 0.0f
,
170 0.0f
, 0.0f
, 1.0f
, 0.0f
,
171 0.0f
, 0.0f
, 0.0f
, 1.0f
174 const struct vl_procamp vl_default_procamp
= {
175 0.0f
, /* brightness */
177 1.0f
, /* saturation */
181 void vl_csc_get_matrix(enum VL_CSC_COLOR_STANDARD cs
,
182 struct vl_procamp
*procamp
,
186 float ybias
= full_range
? -16.0f
/255.0f
: 0.0f
;
187 float cbbias
= -128.0f
/255.0f
;
188 float crbias
= -128.0f
/255.0f
;
190 const struct vl_procamp
*p
= procamp
? procamp
: &vl_default_procamp
;
191 float c
= p
->contrast
;
192 float s
= p
->saturation
;
193 float b
= p
->brightness
;
201 case VL_CSC_COLOR_STANDARD_BT_601
:
202 cstd
= full_range
? &bt_601_full
[0] : &bt_601
[0];
204 case VL_CSC_COLOR_STANDARD_BT_709
:
205 cstd
= full_range
? &bt_709_full
[0] : &bt_709
[0];
207 case VL_CSC_COLOR_STANDARD_SMPTE_240M
:
208 cstd
= full_range
? &smpte240m_full
[0] : &smpte240m
[0];
210 case VL_CSC_COLOR_STANDARD_IDENTITY
:
212 assert(cs
== VL_CSC_COLOR_STANDARD_IDENTITY
);
213 memcpy(matrix
, &identity
[0], sizeof(float) * 16);
217 matrix
[ 0] = c
*cstd
[ 0];
218 matrix
[ 1] = c
*cstd
[ 1]*s
*cosf(h
) - c
*cstd
[ 2]*s
*sinf(h
);
219 matrix
[ 2] = c
*cstd
[ 2]*s
*cosf(h
) + c
*cstd
[ 1]*s
*sinf(h
);
220 matrix
[ 3] = cstd
[ 3] + cstd
[ 0]*(b
+ c
*ybias
) + cstd
[ 1]*(c
*cbbias
*s
*cosf(h
) + c
*crbias
*s
*sinf(h
)) + cstd
[ 2]*(c
*crbias
*s
*cosf(h
) - c
*cbbias
*s
*sinf(h
));
222 matrix
[ 4] = c
*cstd
[ 4];
223 matrix
[ 5] = c
*cstd
[ 5]*s
*cosf(h
) - c
*cstd
[ 6]*s
*sinf(h
);
224 matrix
[ 6] = c
*cstd
[ 6]*s
*cosf(h
) + c
*cstd
[ 5]*s
*sinf(h
);
225 matrix
[ 7] = cstd
[ 7] + cstd
[ 4]*(b
+ c
*ybias
) + cstd
[ 5]*(c
*cbbias
*s
*cosf(h
) + c
*crbias
*s
*sinf(h
)) + cstd
[ 6]*(c
*crbias
*s
*cosf(h
) - c
*cbbias
*s
*sinf(h
));
227 matrix
[ 8] = c
*cstd
[ 8];
228 matrix
[ 9] = c
*cstd
[ 9]*s
*cosf(h
) - c
*cstd
[10]*s
*sinf(h
);
229 matrix
[10] = c
*cstd
[10]*s
*cosf(h
) + c
*cstd
[ 9]*s
*sinf(h
);
230 matrix
[11] = cstd
[11] + cstd
[ 8]*(b
+ c
*ybias
) + cstd
[ 9]*(c
*cbbias
*s
*cosf(h
) + c
*crbias
*s
*sinf(h
)) + cstd
[10]*(c
*crbias
*s
*cosf(h
) - c
*cbbias
*s
*sinf(h
));
232 matrix
[12] = c
*cstd
[12];
233 matrix
[13] = c
*cstd
[13]*s
*cos(h
) - c
*cstd
[14]*s
*sin(h
);
234 matrix
[14] = c
*cstd
[14]*s
*cos(h
) + c
*cstd
[13]*s
*sin(h
);
235 matrix
[15] = cstd
[15] + cstd
[12]*(b
+ c
*ybias
) + cstd
[13]*(c
*cbbias
*s
*cos(h
) + c
*crbias
*s
*sin(h
)) + cstd
[14]*(c
*crbias
*s
*cos(h
) - c
*cbbias
*s
*sin(h
));