** Copyright 2004, Marcus Overhagen. All rights reserved.
** Distributed under the terms of the MIT license.
**
** This file generates the lookup tables used by image
** format conversion from planar YCbCr420 format to
** linear RGB32 format in file gfx_conv_c_lookup.cpp
*/
#include <strings.h>
#include <stdio.h>
double precision = 32768.0;
int shift = 15;
R = 1.164(Y - 16) + 1.596(Cr - 128)
G = 1.164(Y - 16) - 0.813(Cr - 128) - 0.392(Cb - 128)
B = 1.164(Y - 16) + 2.017(Cb - 128)
*/
int main()
{
printf("const int32 Cb_Gtab[256] = {\n\t");
for (int i = 0; i < 256; i++) {
int Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
int Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
int Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
int Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
printf("%d, ", Cb_G);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n\n");
printf("const int32 Cb_Btab[256] = {\n\t");
for (int i = 0; i < 256; i++) {
int Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
int Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
int Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
int Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
printf("%d, ", Cb_B);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n\n");
printf("const int32 Cr_Rtab[256] = {\n\t");
for (int i = 0; i < 256; i++) {
int Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
int Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
int Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
int Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
printf("%d, ", Cr_R);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n\n");
printf("const int32 Cr_Gtab[256] = {\n\t");
for (int i = 0; i < 256; i++) {
int Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
int Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
int Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
int Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
printf("%d, ", Cr_G);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n\n");
printf("const int32 Ytab[256] = {\n\t");
for (int i = 0; i < 256; i++) {
int Y = (i - 16) * 1.164 * precision + 0.5;
printf("%d, ", Y);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n\n");
int Ymin = (0 - 16) * 1.164 * precision + 0.5;
int Ymax = (255 - 16) * 1.164 * precision + 0.5;
int Cr_Rmin = (- 128 + 0) * 1.596 * precision + 0.5;
int Cr_Rmax = (- 128 + 255) * 1.596 * precision + 0.5;
int Cr_Gmin = (- 128 + 255) * -0.813 * precision + 0.5;
int Cr_Gmax = (- 128 + 0) * -0.813 * precision + 0.5;
int Cb_Gmin = (- 128 + 255) * -0.392 * precision + 0.5;
int Cb_Gmax = (- 128 + 0) * -0.392 * precision + 0.5;
int Cb_Bmin = (- 128 + 0) * 2.017 * precision + 0.5;
int Cb_Bmax = (- 128 + 255) * 2.017 * precision + 0.5;
printf("Cr_Rmin %d\n", Cr_Rmin);
printf("Cr_Rmax %d\n", Cr_Rmax);
printf("Cr_Gmin %d\n", Cr_Gmin);
printf("Cr_Gmax %d\n", Cr_Gmax);
printf("Cb_Gmin %d\n", Cb_Gmin);
printf("Cb_Gmax %d\n", Cb_Gmax);
printf("Cb_Bmin %d\n", Cb_Bmin);
printf("Cb_Bmax %d\n", Cb_Bmax);
int Rmax = (Ymax + Cr_Rmax) >> shift;
int Rmin = (Ymin + Cr_Rmin) >> shift;
int Gmax = (Ymax + Cr_Gmax + Cb_Gmax) >> shift;
int Gmin = (Ymin + Cr_Gmin + Cb_Gmin) >> shift;
int Bmax = (Ymax + Cb_Bmax) >> shift;
int Bmin = (Ymin + Cb_Bmin) >> shift;
printf("Rmax %d\n", Rmax);
printf("Rmin %d\n", Rmin);
printf("Gmax %d\n", Gmax);
printf("Gmin %d\n", Gmin);
printf("Bmax %d\n", Bmax);
printf("Bmin %d\n", Bmin);
int num_Rneg = (Rmin < 0) ? -Rmin : 0;
int num_Rpos = (Rmax > 255) ? Rmax - 255 : 0;
int num_Rent = num_Rneg + 256 + num_Rpos;
printf("num_Rneg %d\n", num_Rneg);
printf("num_Rpos %d\n", num_Rpos);
printf("num_Rent %d\n", num_Rent);
printf("const uint32 Rsat[%d] = {\n\t", num_Rent);
for (int i = 0; i < num_Rneg; i++) {
printf("0x000000, ");
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n\t");
for (int i = 0; i < 256; i++) {
printf("0x%06x, ", i << 16);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n\t");
for (int i = 0; i < num_Rpos; i++) {
printf("0x%06x, ", 255 << 16);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n");
printf("const uint32 *Rtab = &Rsat[%d];\n\n", num_Rneg);
int num_Gneg = (Gmin < 0) ? -Gmin : 0;
int num_Gpos = (Gmax > 255) ? Gmax - 255 : 0;
int num_Gent = num_Gneg + 256 + num_Gpos;
printf("num_Gneg %d\n", num_Gneg);
printf("num_Gpos %d\n", num_Gpos);
printf("num_Gent %d\n", num_Gent);
printf("const uint32 Gsat[%d] = {\n\t", num_Gent);
for (int i = 0; i < num_Gneg; i++) {
printf("0x000000, ");
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n\t");
for (int i = 0; i < 256; i++) {
printf("0x%06x, ", i << 8);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n\t");
for (int i = 0; i < num_Gpos; i++) {
printf("0x%06x, ", 255 << 8);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n");
printf("const uint32 *Gtab = &Gsat[%d];\n\n", num_Gneg);
int num_Bneg = (Bmin < 0) ? -Bmin : 0;
int num_Bpos = (Bmax > 255) ? Bmax - 255 : 0;
int num_Bent = num_Bneg + 256 + num_Bpos;
printf("num_Bneg %d\n", num_Bneg);
printf("num_Bpos %d\n", num_Bpos);
printf("num_Bent %d\n", num_Bent);
printf("const uint32 Bsat[%d] = {\n\t", num_Bent);
for (int i = 0; i < num_Bneg; i++) {
printf("0x000000, ");
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n\t");
for (int i = 0; i < 256; i++) {
printf("0x%06x, ", i << 0);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n\t");
for (int i = 0; i < num_Bpos; i++) {
printf("0x%06x, ", 255 << 0);
if ((i % 10) == 9)
printf("\n\t");
}
printf("\n};\n");
printf("const uint32 *Btab = &Bsat[%d];\n\n", num_Bneg);
return 0;
}
