纯C实现的JPEG压缩算法
2010-09-03 22:19
190 查看
#include <stdio.h>
typedef unsigned int UINT32;
typedef int INT32;
typedef unsigned short UINT16;
typedef short INT16;
typedef unsigned char UINT8;
typedef char INT8;
#define CLIP(color) (unsigned char)(((color)>0xFF)?0xff:(((color)<0)?0:(color)))
#define BLOCK_SIZE 64
typedef struct JPEG_ENCODER_STRUCTURE
{
UINT16 mcu_width;
UINT16 mcu_height;
UINT16 horizontal_mcus;
UINT16 vertical_mcus;
UINT16 rows_in_bottom_mcus;
UINT16 cols_in_right_mcus;
UINT16 length_minus_mcu_width;
UINT16 length_minus_width;
UINT16 mcu_width_size;
UINT16 offset;
INT16 ldc1;
INT16 ldc2;
INT16 ldc3;
UINT16 rows;
UINT16 cols;
UINT16 incr;
} JPEG_ENCODER_STRUCTURE;
UINT8 Lqt[BLOCK_SIZE];
UINT8 Cqt[BLOCK_SIZE];
UINT16 ILqt[BLOCK_SIZE];
UINT16 ICqt[BLOCK_SIZE];
INT16 Y1[BLOCK_SIZE];
INT16 Y2[BLOCK_SIZE];
INT16 Y3[BLOCK_SIZE];
INT16 Y4[BLOCK_SIZE];
INT16 CB[BLOCK_SIZE];
INT16 CR[BLOCK_SIZE];
INT16 Temp[BLOCK_SIZE];
INT32 lcode;
UINT16 bitindex;
static UINT8 zigzag_table[] = {
0, 1, 5, 6, 14, 15, 27, 28, 2, 4, 7, 13, 16, 26, 29, 42, 3, 8, 12, 17,
25, 30, 41, 43, 9, 11, 18, 24, 31, 40, 44, 53, 10, 19, 23, 32, 39, 45, 52,
54, 20, 22, 33, 38, 46, 51, 55, 60, 21, 34, 37, 47, 50, 56, 59, 61, 35, 36,
48, 49, 57, 58, 62, 63
};
// Header for JPEG Encoder
static UINT16 markerdata[] = {
// dht
0xFFC4, 0x1A2, 0x00,
// luminance dc (2 - 16) + 1
0x0105, 0x0101, 0x00101, 0x0101, 0x0000, 0x00000, 00000, 00000,
// luminance dc (2 - 12) + 1
0x0102, 0x0304, 0x0506, 0x0708, 0x090A, 0x0B01,
// chrominance dc (1 - 16)
0x0003, 0x0101, 0x0101, 0x0101, 0x0101, 0x0100, 0x0000, 0x0000,
// chrominance dc (1 - 12)
0x0001, 0x00203, 0x0405, 0x0607, 0x0809, 0x00A0B,
// luminance ac 1 + (1 - 15)
0x1000, 0x0201, 0x0303, 0x0204, 0x0305, 0x0504, 0x0400, 0x0001,
// luminance ac 1 + (1 - 162) + 1
0x7D01, 0x0203, 0x0004, 0x1105, 0x1221, 0x3141, 0x0613, 0x5161, 0x0722,
0x7114, 0x3281, 0x91A1, 0x0823, 0x42B1, 0xC115, 0x52D1, 0xF024, 0x3362,
0x7282, 0x090A, 0x1617, 0x1819, 0x1A25, 0x2627, 0x2829, 0x2A34, 0x3536,
0x3738, 0x393A, 0x4344, 0x4546, 0x4748, 0x494A, 0x5354, 0x5556, 0x5758,
0x595A, 0x6364, 0x6566, 0x6768, 0x696A, 0x7374, 0x7576, 0x7778, 0x797A,
0x8384, 0x8586, 0x8788, 0x898A, 0x9293, 0x9495, 0x9697, 0x9899, 0x9AA2,
0xA3A4, 0xA5A6, 0xA7A8, 0xA9AA, 0xB2B3, 0xB4B5, 0xB6B7, 0xB8B9, 0xBAC2,
0xC3C4, 0xC5C6, 0xC7C8, 0xC9CA, 0xD2D3, 0xD4D5, 0xD6D7, 0xD8D9, 0xDAE1,
0xE2E3, 0xE4E5, 0xE6E7, 0xE8E9, 0xEAF1, 0xF2F3, 0xF4F5, 0xF6F7, 0xF8F9,
0xFA11,
// chrominance ac (1 - 16)
0x0002, 0x0102, 0x0404, 0x0304, 0x0705, 0x0404, 0x0001, 0x0277,
// chrominance ac (1 - 162)
0x0001, 0x0203, 0x1104, 0x0521, 0x3106, 0x1241, 0x5107, 0x6171, 0x1322,
0x3281, 0x0814, 0x4291, 0xA1B1, 0xC109, 0x2333, 0x52F0, 0x1562, 0x72D1,
0x0A16, 0x2434, 0xE125, 0xF117, 0x1819, 0x1A26, 0x2728, 0x292A, 0x3536,
0x3738, 0x393A, 0x4344, 0x4546, 0x4748, 0x494A, 0x5354, 0x5556, 0x5758,
0x595A, 0x6364, 0x6566, 0x6768, 0x696A, 0x7374, 0x7576, 0x7778, 0x797A,
0x8283, 0x8485, 0x8687, 0x8889, 0x8A92, 0x9394, 0x9596, 0x9798, 0x999A,
0xA2A3, 0xA4A5, 0xA6A7, 0xA8A9, 0xAAB2, 0xB3B4, 0xB5B6, 0xB7B8, 0xB9BA,
0xC2C3, 0xC4C5, 0xC6C7, 0xC8C9, 0xCAD2, 0xD3D4, 0xD5D6, 0xD7D8, 0xD9DA,
0xE2E3, 0xE4E5, 0xE6E7, 0xE8E9, 0xEAF2, 0xF3F4, 0xF5F6, 0xF7F8, 0xF9FA
};
UINT8 bitsize[] = {
0, 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6,
6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8
};
UINT16 luminance_dc_code_table[] = {
0x0000, 0x0002, 0x0003, 0x0004, 0x0005, 0x0006, 0x000E, 0x001E, 0x003E,
0x007E, 0x00FE, 0x01FE
};
UINT16 luminance_dc_size_table[] = {
0x0002, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0004, 0x0005, 0x0006,
0x0007, 0x0008, 0x0009
};
UINT16 chrominance_dc_code_table[] = {
0x0000, 0x0001, 0x0002, 0x0006, 0x000E, 0x001E, 0x003E, 0x007E, 0x00FE,
0x01FE, 0x03FE, 0x07FE
};
UINT16 chrominance_dc_size_table[] = {
0x0002, 0x0002, 0x0002, 0x0003, 0x0004, 0x0005, 0x0006, 0x0007, 0x0008,
0x0009, 0x000A, 0x000B
};
UINT16 luminance_ac_code_table[] = {
0x000A, 0x0000, 0x0001, 0x0004, 0x000B, 0x001A, 0x0078, 0x00F8, 0x03F6,
0xFF82, 0xFF83, 0x000C, 0x001B, 0x0079, 0x01F6, 0x07F6, 0xFF84, 0xFF85,
0xFF86, 0xFF87, 0xFF88, 0x001C, 0x00F9, 0x03F7, 0x0FF4, 0xFF89, 0xFF8A,
0xFF8b, 0xFF8C, 0xFF8D, 0xFF8E, 0x003A, 0x01F7, 0x0FF5, 0xFF8F, 0xFF90,
0xFF91, 0xFF92, 0xFF93, 0xFF94, 0xFF95, 0x003B, 0x03F8, 0xFF96, 0xFF97,
0xFF98, 0xFF99, 0xFF9A, 0xFF9B, 0xFF9C, 0xFF9D, 0x007A, 0x07F7, 0xFF9E,
0xFF9F, 0xFFA0, 0xFFA1, 0xFFA2, 0xFFA3, 0xFFA4, 0xFFA5, 0x007B, 0x0FF6,
0xFFA6, 0xFFA7, 0xFFA8, 0xFFA9, 0xFFAA, 0xFFAB, 0xFFAC, 0xFFAD, 0x00FA,
0x0FF7, 0xFFAE, 0xFFAF, 0xFFB0, 0xFFB1, 0xFFB2, 0xFFB3, 0xFFB4, 0xFFB5,
0x01F8, 0x7FC0, 0xFFB6, 0xFFB7, 0xFFB8, 0xFFB9, 0xFFBA, 0xFFBB, 0xFFBC,
0xFFBD, 0x01F9, 0xFFBE, 0xFFBF, 0xFFC0, 0xFFC1, 0xFFC2, 0xFFC3, 0xFFC4,
0xFFC5, 0xFFC6, 0x01FA, 0xFFC7, 0xFFC8, 0xFFC9, 0xFFCA, 0xFFCB, 0xFFCC,
0xFFCD, 0xFFCE, 0xFFCF, 0x03F9, 0xFFD0, 0xFFD1, 0xFFD2, 0xFFD3, 0xFFD4,
0xFFD5, 0xFFD6, 0xFFD7, 0xFFD8, 0x03FA, 0xFFD9, 0xFFDA, 0xFFDB, 0xFFDC,
0xFFDD, 0xFFDE, 0xFFDF, 0xFFE0, 0xFFE1, 0x07F8, 0xFFE2, 0xFFE3, 0xFFE4,
0xFFE5, 0xFFE6, 0xFFE7, 0xFFE8, 0xFFE9, 0xFFEA, 0xFFEB, 0xFFEC, 0xFFED,
0xFFEE, 0xFFEF, 0xFFF0, 0xFFF1, 0xFFF2, 0xFFF3, 0xFFF4, 0xFFF5, 0xFFF6,
0xFFF7, 0xFFF8, 0xFFF9, 0xFFFA, 0xFFFB, 0xFFFC, 0xFFFD, 0xFFFE,
0x07F9
};
UINT16 luminance_ac_size_table[] = {
0x0004, 0x0002, 0x0002, 0x0003, 0x0004, 0x0005, 0x0007, 0x0008, 0x000A,
0x0010, 0x0010, 0x0004, 0x0005, 0x0007, 0x0009, 0x000B, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0005, 0x0008, 0x000A, 0x000C, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x0009, 0x000C, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x000A, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007, 0x000B, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007, 0x000C,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0008,
0x000C, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0009, 0x000F, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x000A, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x000A, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000B, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x000B
};
UINT16 chrominance_ac_code_table[] = {
0x0000, 0x0001, 0x0004, 0x000A, 0x0018, 0x0019, 0x0038, 0x0078, 0x01F4,
0x03F6, 0x0FF4, 0x000B, 0x0039, 0x00F6, 0x01F5, 0x07F6, 0x0FF5, 0xFF88,
0xFF89, 0xFF8A, 0xFF8B, 0x001A, 0x00F7, 0x03F7, 0x0FF6, 0x7FC2, 0xFF8C,
0xFF8D, 0xFF8E, 0xFF8F, 0xFF90, 0x001B, 0x00F8, 0x03F8, 0x0FF7, 0xFF91,
0xFF92, 0xFF93, 0xFF94, 0xFF95, 0xFF96, 0x003A, 0x01F6, 0xFF97, 0xFF98,
0xFF99, 0xFF9A, 0xFF9B, 0xFF9C, 0xFF9D, 0xFF9E, 0x003B, 0x03F9, 0xFF9F,
0xFFA0, 0xFFA1, 0xFFA2, 0xFFA3, 0xFFA4, 0xFFA5, 0xFFA6, 0x0079, 0x07F7,
0xFFA7, 0xFFA8, 0xFFA9, 0xFFAA, 0xFFAB, 0xFFAC, 0xFFAD, 0xFFAE, 0x007A,
0x07F8, 0xFFAF, 0xFFB0, 0xFFB1, 0xFFB2, 0xFFB3, 0xFFB4, 0xFFB5, 0xFFB6,
0x00F9, 0xFFB7, 0xFFB8, 0xFFB9, 0xFFBA, 0xFFBB, 0xFFBC, 0xFFBD, 0xFFBE,
0xFFBF, 0x01F7, 0xFFC0, 0xFFC1, 0xFFC2, 0xFFC3, 0xFFC4, 0xFFC5, 0xFFC6,
0xFFC7, 0xFFC8, 0x01F8, 0xFFC9, 0xFFCA, 0xFFCB, 0xFFCC, 0xFFCD, 0xFFCE,
0xFFCF, 0xFFD0, 0xFFD1, 0x01F9, 0xFFD2, 0xFFD3, 0xFFD4, 0xFFD5, 0xFFD6,
0xFFD7, 0xFFD8, 0xFFD9, 0xFFDA, 0x01FA, 0xFFDB, 0xFFDC, 0xFFDD, 0xFFDE,
0xFFDF, 0xFFE0, 0xFFE1, 0xFFE2, 0xFFE3, 0x07F9, 0xFFE4, 0xFFE5, 0xFFE6,
0xFFE7, 0xFFE8, 0xFFE9, 0xFFEA, 0xFFEb, 0xFFEC, 0x3FE0, 0xFFED, 0xFFEE,
0xFFEF, 0xFFF0, 0xFFF1, 0xFFF2, 0xFFF3, 0xFFF4, 0xFFF5, 0x7FC3, 0xFFF6,
0xFFF7, 0xFFF8, 0xFFF9, 0xFFFA, 0xFFFB, 0xFFFC, 0xFFFD, 0xFFFE,
0x03FA
};
UINT16 chrominance_ac_size_table[] = {
0x0002, 0x0002, 0x0003, 0x0004, 0x0005, 0x0005, 0x0006, 0x0007, 0x0009,
0x000A, 0x000C, 0x0004, 0x0006, 0x0008, 0x0009, 0x000B, 0x000C, 0x0010,
0x0010, 0x0010, 0x0010, 0x0005, 0x0008, 0x000A, 0x000C, 0x000F, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0005, 0x0008, 0x000A, 0x000C, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x0009, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x000A, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007, 0x000B,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007,
0x000B, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0008, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000B, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000E, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000F, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x000A
};
#define FOUR_ZERO_ZERO 0 // Grey scale Y00 ...
#define FOUR_TWO_TWO 2 //Y00 Cb Y01 Cr
static void (*read_format) (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT8 * input_ptr);
static void RGB_2_400 (UINT8 * input_ptr, UINT8 * output_ptr, UINT32 image_width,UINT32 image_height)
{
UINT32 i, size;
UINT8 R, G, B;
INT32 Y;
UINT8 * inbuf = input_ptr;
size = image_width * image_height;
for (i = size; i > 0; i--)
{
B = inbuf[0];
G = inbuf[1];
R = inbuf[2];
inbuf += 3;
Y = CLIP ((77 * R + 150 * G + 29 * B) >> 8);
*input_ptr++ = (UINT8) Y;
}
}
static void read_400_format (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT8 * input_ptr)
{
INT32 i, j;
INT16 * Y1_Ptr = Y1;
UINT16 rows = jpeg_encoder_structure->rows;
UINT16 cols = jpeg_encoder_structure->cols;
UINT16 incr = jpeg_encoder_structure->incr;
for (i = rows; i > 0; i--)
{
for (j = cols; j > 0; j--)
*Y1_Ptr++ = *input_ptr++ - 128;
for (j = 8 - cols; j > 0; j--)
*Y1_Ptr++ = *(Y1_Ptr - 1);
input_ptr += incr;
}
for (i = 8 - rows; i > 0; i--)
{
for (j = 8; j > 0; j--)
*Y1_Ptr++ = *(Y1_Ptr - 8);
}
}
static void
RGB_2_422 (UINT8 * input_ptr, UINT8 * output_ptr, UINT32 image_width,
UINT32 image_height)
{
UINT32 i, size;
UINT8 R, G, B, R1, G1, B1;
INT32 Y, Yp, Cb, Cr;
UINT8 * inbuf = input_ptr;
size = image_width * image_height/2;
for (i = size; i > 0; i--)
{
B = inbuf[0];
G = inbuf[1];
R = inbuf[2];
B1 = inbuf[3];
G1 = inbuf[4];
R1 = inbuf[5];
inbuf += 6;
Y = CLIP ((77 * R + 150 * G + 29 * B) >> 8);
Yp = CLIP ((77 * R1 + 150 * G1 + 29 * B1) >> 8);
Cb = CLIP (((-43 * R - 85 * G + 128 * B) >> 8) + 128);
Cr = CLIP (((128 * R - 107 * G - 21 * B) >> 8) + 128);
*input_ptr++ = (UINT8) Y;
*input_ptr++ = (UINT8) Cb;
*input_ptr++ = (UINT8) Yp;
*input_ptr++ = (UINT8) Cr;
}
}
/* This function implements 16 Step division for Q.15 format data */
UINT16 DSP_Division (UINT32 numer, UINT32 denom)
{
UINT16 i;
denom <<= 15;
for (i = 16; i > 0; i--)
{
if (numer > denom)
{
numer -= denom;
numer <<= 1;
numer++;
}
else
numer <<= 1;
}
return (UINT16) numer;
}
UINT8 * write_markers (UINT8 * output_ptr, UINT32 image_format,UINT32 image_width, UINT32 image_height)
{
UINT16 i, header_length;
UINT8 number_of_components;
// Start of image marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xD8;
// Quantization table marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xDB;
// Quantization table length
*output_ptr++ = 0x00;
*output_ptr++ = 0x84;
// Pq, Tq
*output_ptr++ = 0x00;
// Lqt table
for (i = 0; i < 64; i++)
*output_ptr++ = Lqt[i];
// Pq, Tq
*output_ptr++ = 0x01;
// Cqt table
for (i = 0; i < 64; i++)
*output_ptr++ = Cqt[i];
// huffman table(DHT)
for (i = 0; i < 210; i++)
{
*output_ptr++ = (UINT8) (markerdata[i] >> 8);
*output_ptr++ = (UINT8) markerdata[i];
}
//FOUR_ZERO_ZERO
number_of_components = 3;
// Frame header(SOF)
// Start of frame marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xC0;
header_length = (UINT16) (8 + 3 * number_of_components);
// Frame header length
*output_ptr++ = (UINT8) (header_length >> 8);
*output_ptr++ = (UINT8) header_length;
// Precision (P)
*output_ptr++ = 0x08;
// image height
*output_ptr++ = (UINT8) (image_height >> 8);
*output_ptr++ = (UINT8) image_height;
// image width
*output_ptr++ = (UINT8) (image_width >> 8);
*output_ptr++ = (UINT8) image_width;
// Nf
*output_ptr++ = number_of_components;
// FOUR_ZERO_ZERO
*output_ptr++ = 0x01;
*output_ptr++ = 0x21;
*output_ptr++ = 0x00;
*output_ptr++ = 0x02;
*output_ptr++ = 0x11;
*output_ptr++ = 0x01;
*output_ptr++ = 0x03;
*output_ptr++ = 0x11;
*output_ptr++ = 0x01;
// Scan header(SOF)
// Start of scan marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xDA;
header_length = (UINT16) (6 + (number_of_components << 1));
// Scan header length
*output_ptr++ = (UINT8) (header_length >> 8);
*output_ptr++ = (UINT8) header_length;
// Ns
*output_ptr++ = number_of_components;
//FOUR_ZERO_ZERO
*output_ptr++ = 0x01;
*output_ptr++ = 0x00;
*output_ptr++ = 0x02;
*output_ptr++ = 0x11;
*output_ptr++ = 0x03;
*output_ptr++ = 0x11;
*output_ptr++ = 0x00;
*output_ptr++ = 0x3F;
*output_ptr++ = 0x00;
return output_ptr;
}
static void
read_422_format (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,
UINT8 * input_ptr)
{
INT32 i, j;
UINT16 Y1_cols, Y2_cols;
INT16 * Y1_Ptr = Y1;
INT16 * Y2_Ptr = Y2;
INT16 * CB_Ptr = CB;
INT16 * CR_Ptr = CR;
UINT16 rows = jpeg_encoder_structure->rows;
UINT16 cols = jpeg_encoder_structure->cols;
UINT16 incr = jpeg_encoder_structure->incr;
if (cols <= 8)
{
Y1_cols = cols;
Y2_cols = 0;
}
else
{
Y1_cols = 8;
Y2_cols = (UINT16) (cols - 8);
}
for (i = rows; i > 0; i--)
{
for (j = Y1_cols >> 1; j > 0; j--)
{
*Y1_Ptr++ = *input_ptr++ - 128;
*CB_Ptr++ = *input_ptr++ - 128;
*Y1_Ptr++ = *input_ptr++ - 128;
*CR_Ptr++ = *input_ptr++ - 128;
}
for (j = Y2_cols >> 1; j > 0; j--)
{
*Y2_Ptr++ = *input_ptr++ - 128;
*CB_Ptr++ = *input_ptr++ - 128;
*Y2_Ptr++ = *input_ptr++ - 128;
*CR_Ptr++ = *input_ptr++ - 128;
}
if (cols <= 8)
{
for (j = 8 - Y1_cols; j > 0; j--)
*Y1_Ptr++ = *(Y1_Ptr - 1);
for (j = 8 - Y2_cols; j > 0; j--)
*Y2_Ptr++ = *(Y1_Ptr - 1);
}
else
{
for (j = 8 - Y2_cols; j > 0; j--)
*Y2_Ptr++ = *(Y2_Ptr - 1);
}
for (j = (16 - cols) >> 1; j > 0; j--)
{
*CB_Ptr++ = *(CB_Ptr - 1);
*CR_Ptr++ = *(CR_Ptr - 1);
}
input_ptr += incr;
}
for (i = 8 - rows; i > 0; i--)
{
for (j = 8; j > 0; j--)
{
*Y1_Ptr++ = *(Y1_Ptr - 8);
*Y2_Ptr++ = *(Y2_Ptr - 8);
*CB_Ptr++ = *(CB_Ptr - 8);
*CR_Ptr++ = *(CR_Ptr - 8);
}
}
}
static void initialization (JPEG_ENCODER_STRUCTURE * jpeg, UINT32 image_format,UINT32 image_width, UINT32 image_height)
{
UINT16 mcu_width, mcu_height, bytes_per_pixel;
lcode = 0;
bitindex = 0;
jpeg->mcu_width = mcu_width = 16;
jpeg->horizontal_mcus = (UINT16) ((image_width + mcu_width - 1) >> 4);
jpeg->mcu_height = mcu_height = 8;
jpeg->vertical_mcus =
(UINT16) ((image_height + mcu_height - 1) >> 3);
bytes_per_pixel = 2;
read_format = read_422_format;
jpeg->rows_in_bottom_mcus =(UINT16) (image_height - (jpeg->vertical_mcus - 1) * mcu_height);
jpeg->cols_in_right_mcus =(UINT16) (image_width - (jpeg->horizontal_mcus - 1) * mcu_width);
jpeg->length_minus_mcu_width =(UINT16) ((image_width - mcu_width) * bytes_per_pixel);
jpeg->length_minus_width =(UINT16) ((image_width - jpeg->cols_in_right_mcus) * bytes_per_pixel);
jpeg->mcu_width_size = (UINT16) (mcu_width * bytes_per_pixel);
jpeg->offset =(UINT16) ((image_width * (mcu_height - 1) -(mcu_width - jpeg->cols_in_right_mcus)) * bytes_per_pixel);
jpeg->ldc1 = 0;
jpeg->ldc2 = 0;
jpeg->ldc3 = 0;
}
/* Multiply Quantization table with quality factor to get LQT and CQT */
void initialize_quantization_tables (UINT32 quality_factor)
{
UINT16 i, index;
UINT32 value;
UINT8 luminance_quant_table[] = {
16, 11, 10, 16, 24, 40, 51, 61, 12, 12, 14, 19, 26, 58, 60, 55, 14, 13,
16, 24, 40, 57, 69, 56, 14, 17, 22, 29, 51, 87, 80, 62, 18, 22, 37,
56, 68, 109, 103, 77, 24, 35, 55, 64, 81, 104, 113, 92, 49, 64, 78,
87, 103, 121, 120, 101, 72, 92, 95, 98, 112, 100, 103, 99
};
UINT8 chrominance_quant_table[] = {
17, 18, 24, 47, 99, 99, 99, 99, 18, 21, 26, 66, 99, 99, 99, 99, 24, 26,
56, 99, 99, 99, 99, 99, 47, 66, 99, 99, 99, 99, 99, 99, 99, 99, 99,
99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99,
99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99
};
for (i = 0; i < 64; i++)
{
index = zigzag_table[i];
/* luminance quantization table * quality factor */
value = luminance_quant_table[i] * quality_factor;
value = (value + 0x200) >> 10;
if (value == 0)
value = 1;
else if (value > 255)
value = 255;
Lqt[index] = (UINT8) value;
ILqt[i] = DSP_Division (0x8000, value);
/* chrominance quantization table * quality factor */
value = chrominance_quant_table[i] * quality_factor;
value = (value + 0x200) >> 10;
if (value == 0)
value = 1;
else if (value > 255)
value = 255;
Cqt[index] = (UINT8) value;
ICqt[i] = DSP_Division (0x8000, value);
}
}
UINT8 * close_bitstream (UINT8 * output_ptr)
{
UINT16 i, count;
UINT8 * ptr;
if (bitindex > 0)
{
lcode <<= (32 - bitindex);
count = (bitindex + 7) >> 3;
ptr = (UINT8 *) & lcode + 3;
for (i = count; i > 0; i--)
{
if ((*output_ptr++ = *ptr--) == 0xff)
*output_ptr++ = 0;
}
}
// End of image marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xD9;
return output_ptr;
}
/* DCT for One block(8x8) */
static void DCT (INT16 * data)
{
UINT16 i;
INT32 x0, x1, x2, x3, x4, x5, x6, x7, x8;
/* All values are shifted left by 10
and rounded off to nearest integer */
static const UINT16 c1 = 1420; /* cos PI/16 * root(2) */
static const UINT16 c2 = 1338; /* cos PI/8 * root(2) */
static const UINT16 c3 = 1204; /* cos 3PI/16 * root(2) */
static const UINT16 c5 = 805; /* cos 5PI/16 * root(2) */
static const UINT16 c6 = 554; /* cos 3PI/8 * root(2) */
static const UINT16 c7 = 283; /* cos 7PI/16 * root(2) */
static const UINT16 s1 = 3;
static const UINT16 s2 = 10;
static const UINT16 s3 = 13;
for (i = 8; i > 0; i--)
{
x8 = data[0] + data[7];
x0 = data[0] - data[7];
x7 = data[1] + data[6];
x1 = data[1] - data[6];
x6 = data[2] + data[5];
x2 = data[2] - data[5];
x5 = data[3] + data[4];
x3 = data[3] - data[4];
x4 = x8 + x5;
x8 -= x5;
x5 = x7 + x6;
x7 -= x6;
data[0] = (INT16) (x4 + x5);
data[4] = (INT16) (x4 - x5);
data[2] = (INT16) ((x8 * c2 + x7 * c6) >> s2);
data[6] = (INT16) ((x8 * c6 - x7 * c2) >> s2);
data[7] = (INT16) ((x0 * c7 - x1 * c5 + x2 * c3 - x3 * c1) >> s2);
data[5] = (INT16) ((x0 * c5 - x1 * c1 + x2 * c7 + x3 * c3) >> s2);
data[3] = (INT16) ((x0 * c3 - x1 * c7 - x2 * c1 - x3 * c5) >> s2);
data[1] = (INT16) ((x0 * c1 + x1 * c3 + x2 * c5 + x3 * c7) >> s2);
data += 8;
}
data -= 64;
for (i = 8; i > 0; i--)
{
x8 = data[0] + data[56];
x0 = data[0] - data[56];
x7 = data[8] + data[48];
x1 = data[8] - data[48];
x6 = data[16] + data[40];
x2 = data[16] - data[40];
x5 = data[24] + data[32];
x3 = data[24] - data[32];
x4 = x8 + x5;
x8 -= x5;
x5 = x7 + x6;
x7 -= x6;
data[0] = (INT16) ((x4 + x5) >> s1);
data[32] = (INT16) ((x4 - x5) >> s1);
data[16] = (INT16) ((x8 * c2 + x7 * c6) >> s3);
data[48] = (INT16) ((x8 * c6 - x7 * c2) >> s3);
data[56] = (INT16) ((x0 * c7 - x1 * c5 + x2 * c3 - x3 * c1) >> s3);
data[40] = (INT16) ((x0 * c5 - x1 * c1 + x2 * c7 + x3 * c3) >> s3);
data[24] = (INT16) ((x0 * c3 - x1 * c7 - x2 * c1 - x3 * c5) >> s3);
data[8] = (INT16) ((x0 * c1 + x1 * c3 + x2 * c5 + x3 * c7) >> s3);
data++;
}
}
/* multiply DCT Coefficients with Quantization table and store in ZigZag location */
void quantization (INT16 * const data, UINT16 * const quant_table_ptr)
{
INT16 i;
INT32 value;
for (i = 63; i >= 0; i--)
{
value = data[i] * quant_table_ptr[i];
value = (value + 0x4000) >> 15;
Temp[zigzag_table[i]] = (INT16) value;
}
}
UINT8 * huffman (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT16 component, UINT8 * output_ptr)
{
UINT16 i;
UINT16 * DcCodeTable, *DcSizeTable, *AcCodeTable, *AcSizeTable;
INT16 * Temp_Ptr, Coeff, LastDc;
UINT16 AbsCoeff, HuffCode, HuffSize, RunLength = 0, DataSize = 0, index;
INT16 bits_in_next_word;
UINT16 numbits;
UINT32 data;
Temp_Ptr = Temp;
Coeff = *Temp_Ptr++;
if (component == 1)
{
DcCodeTable = luminance_dc_code_table;
DcSizeTable = luminance_dc_size_table;
AcCodeTable = luminance_ac_code_table;
AcSizeTable = luminance_ac_size_table;
LastDc = jpeg_encoder_structure->ldc1;
jpeg_encoder_structure->ldc1 = Coeff;
}
else
{
DcCodeTable = chrominance_dc_code_table;
DcSizeTable = chrominance_dc_size_table;
AcCodeTable = chrominance_ac_code_table;
AcSizeTable = chrominance_ac_size_table;
if (component == 2)
{
LastDc = jpeg_encoder_structure->ldc2;
jpeg_encoder_structure->ldc2 = Coeff;
}
else
{
LastDc = jpeg_encoder_structure->ldc3;
jpeg_encoder_structure->ldc3 = Coeff;
}
}
Coeff -= LastDc;
AbsCoeff = (Coeff < 0) ? -Coeff-- : Coeff;
while (AbsCoeff != 0)
{
AbsCoeff >>= 1;
DataSize++;
}
HuffCode = DcCodeTable[DataSize];
HuffSize = DcSizeTable[DataSize];
Coeff &= (1 << DataSize) - 1;
data = (HuffCode << DataSize) | Coeff;
numbits = HuffSize + DataSize;
//PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
for (i = 63; i > 0; i--)
{
if ((Coeff = *Temp_Ptr++) != 0)
{
while (RunLength > 15)
{
RunLength -= 16;
data = AcCodeTable[161];
numbits = AcSizeTable[161];
//PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
}
AbsCoeff = (Coeff < 0) ? -Coeff-- : Coeff;
if (AbsCoeff >> 8 == 0)
DataSize = bitsize[AbsCoeff];
else
DataSize = bitsize[AbsCoeff >> 8] + 8;
index = RunLength * 10 + DataSize;
HuffCode = AcCodeTable[index];
HuffSize = AcSizeTable[index];
Coeff &= (1 << DataSize) - 1;
data = (HuffCode << DataSize) | Coeff;
numbits = HuffSize + DataSize;
// PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
RunLength = 0;
}
else
RunLength++;
}
if (RunLength != 0)
{
data = AcCodeTable[0];
numbits = AcSizeTable[0];
// PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
}
return output_ptr;
}
static UINT8 *encodeMCU (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT32 image_format, UINT8 * output_ptr)
{
DCT (Y1);
quantization (Y1, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
if (image_format == FOUR_ZERO_ZERO)
return output_ptr;
DCT (Y2);
quantization (Y2, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
if (image_format == FOUR_TWO_TWO)
goto chroma;
DCT (Y3);
quantization (Y3, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
DCT (Y4);
quantization (Y4, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
chroma:DCT (CB);
quantization (CB, ICqt);
output_ptr = huffman (jpeg_encoder_structure, 2, output_ptr);
DCT (CR);
quantization (CR, ICqt);
output_ptr = huffman (jpeg_encoder_structure, 3, output_ptr);
return output_ptr;
}
UINT32 encode_image (UINT8 * input_ptr, UINT8 * output_ptr,UINT32 quality_factor, UINT32 image_format,UINT32 image_width, UINT32 image_height)
{
UINT16 i, j;
UINT8 * output;
JPEG_ENCODER_STRUCTURE JpegStruct;
JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure = &JpegStruct;
output = output_ptr;
image_format = FOUR_TWO_TWO;
RGB_2_422 (input_ptr, output_ptr, image_width, image_height);
/* Initialization of JPEG control structure */
initialization (jpeg_encoder_structure, image_format, image_width,
image_height);
/* Quantization Table Initialization */
initialize_quantization_tables (quality_factor);
/* Writing Marker Data */
output_ptr =
write_markers (output_ptr, image_format, image_width, image_height);
for (i = 1; i <= jpeg_encoder_structure->vertical_mcus; i++)
{
if (i < jpeg_encoder_structure->vertical_mcus)
jpeg_encoder_structure->rows = jpeg_encoder_structure->mcu_height;
else
jpeg_encoder_structure->rows =
jpeg_encoder_structure->rows_in_bottom_mcus;
for (j = 1; j <= jpeg_encoder_structure->horizontal_mcus; j++)
{
if (j < jpeg_encoder_structure->horizontal_mcus)
{
jpeg_encoder_structure->cols =
jpeg_encoder_structure->mcu_width;
jpeg_encoder_structure->incr =
jpeg_encoder_structure->length_minus_mcu_width;
}
else
{
jpeg_encoder_structure->cols =
jpeg_encoder_structure->cols_in_right_mcus;
jpeg_encoder_structure->incr =
jpeg_encoder_structure->length_minus_width;
}
read_format (jpeg_encoder_structure, input_ptr);
/* Encode the data in MCU */
output_ptr =encodeMCU (jpeg_encoder_structure, image_format, output_ptr);
input_ptr += jpeg_encoder_structure->mcu_width_size;
}
input_ptr += jpeg_encoder_structure->offset;
}
/* Close Routine */
output_ptr = close_bitstream (output_ptr);
return (UINT32) (output_ptr - output);
}
void main()
{
int size=1280*720*3,sizeout=0,x,y;
unsigned char *outpict ;
outpict = (unsigned char *) malloc (size);
unsigned char *inpict ;
inpict = (unsigned char *) malloc (size);
for(x=0;x<1280;x++)//inpict初始化成RGB模式
for (y=0;y<720;y++)
{
*(inpict+x*3+y*1280*3)=y%2;
*(inpict+x*3+y*1280*3+1)=x%2;
*(inpict+x*3+y*1280*3+2)=y%2;
}
sizeout = encode_image (inpict,outpict,1,0,1280,720);
free(outpict);
free(inpict);
return 1;
}
typedef unsigned int UINT32;
typedef int INT32;
typedef unsigned short UINT16;
typedef short INT16;
typedef unsigned char UINT8;
typedef char INT8;
#define CLIP(color) (unsigned char)(((color)>0xFF)?0xff:(((color)<0)?0:(color)))
#define BLOCK_SIZE 64
typedef struct JPEG_ENCODER_STRUCTURE
{
UINT16 mcu_width;
UINT16 mcu_height;
UINT16 horizontal_mcus;
UINT16 vertical_mcus;
UINT16 rows_in_bottom_mcus;
UINT16 cols_in_right_mcus;
UINT16 length_minus_mcu_width;
UINT16 length_minus_width;
UINT16 mcu_width_size;
UINT16 offset;
INT16 ldc1;
INT16 ldc2;
INT16 ldc3;
UINT16 rows;
UINT16 cols;
UINT16 incr;
} JPEG_ENCODER_STRUCTURE;
UINT8 Lqt[BLOCK_SIZE];
UINT8 Cqt[BLOCK_SIZE];
UINT16 ILqt[BLOCK_SIZE];
UINT16 ICqt[BLOCK_SIZE];
INT16 Y1[BLOCK_SIZE];
INT16 Y2[BLOCK_SIZE];
INT16 Y3[BLOCK_SIZE];
INT16 Y4[BLOCK_SIZE];
INT16 CB[BLOCK_SIZE];
INT16 CR[BLOCK_SIZE];
INT16 Temp[BLOCK_SIZE];
INT32 lcode;
UINT16 bitindex;
static UINT8 zigzag_table[] = {
0, 1, 5, 6, 14, 15, 27, 28, 2, 4, 7, 13, 16, 26, 29, 42, 3, 8, 12, 17,
25, 30, 41, 43, 9, 11, 18, 24, 31, 40, 44, 53, 10, 19, 23, 32, 39, 45, 52,
54, 20, 22, 33, 38, 46, 51, 55, 60, 21, 34, 37, 47, 50, 56, 59, 61, 35, 36,
48, 49, 57, 58, 62, 63
};
// Header for JPEG Encoder
static UINT16 markerdata[] = {
// dht
0xFFC4, 0x1A2, 0x00,
// luminance dc (2 - 16) + 1
0x0105, 0x0101, 0x00101, 0x0101, 0x0000, 0x00000, 00000, 00000,
// luminance dc (2 - 12) + 1
0x0102, 0x0304, 0x0506, 0x0708, 0x090A, 0x0B01,
// chrominance dc (1 - 16)
0x0003, 0x0101, 0x0101, 0x0101, 0x0101, 0x0100, 0x0000, 0x0000,
// chrominance dc (1 - 12)
0x0001, 0x00203, 0x0405, 0x0607, 0x0809, 0x00A0B,
// luminance ac 1 + (1 - 15)
0x1000, 0x0201, 0x0303, 0x0204, 0x0305, 0x0504, 0x0400, 0x0001,
// luminance ac 1 + (1 - 162) + 1
0x7D01, 0x0203, 0x0004, 0x1105, 0x1221, 0x3141, 0x0613, 0x5161, 0x0722,
0x7114, 0x3281, 0x91A1, 0x0823, 0x42B1, 0xC115, 0x52D1, 0xF024, 0x3362,
0x7282, 0x090A, 0x1617, 0x1819, 0x1A25, 0x2627, 0x2829, 0x2A34, 0x3536,
0x3738, 0x393A, 0x4344, 0x4546, 0x4748, 0x494A, 0x5354, 0x5556, 0x5758,
0x595A, 0x6364, 0x6566, 0x6768, 0x696A, 0x7374, 0x7576, 0x7778, 0x797A,
0x8384, 0x8586, 0x8788, 0x898A, 0x9293, 0x9495, 0x9697, 0x9899, 0x9AA2,
0xA3A4, 0xA5A6, 0xA7A8, 0xA9AA, 0xB2B3, 0xB4B5, 0xB6B7, 0xB8B9, 0xBAC2,
0xC3C4, 0xC5C6, 0xC7C8, 0xC9CA, 0xD2D3, 0xD4D5, 0xD6D7, 0xD8D9, 0xDAE1,
0xE2E3, 0xE4E5, 0xE6E7, 0xE8E9, 0xEAF1, 0xF2F3, 0xF4F5, 0xF6F7, 0xF8F9,
0xFA11,
// chrominance ac (1 - 16)
0x0002, 0x0102, 0x0404, 0x0304, 0x0705, 0x0404, 0x0001, 0x0277,
// chrominance ac (1 - 162)
0x0001, 0x0203, 0x1104, 0x0521, 0x3106, 0x1241, 0x5107, 0x6171, 0x1322,
0x3281, 0x0814, 0x4291, 0xA1B1, 0xC109, 0x2333, 0x52F0, 0x1562, 0x72D1,
0x0A16, 0x2434, 0xE125, 0xF117, 0x1819, 0x1A26, 0x2728, 0x292A, 0x3536,
0x3738, 0x393A, 0x4344, 0x4546, 0x4748, 0x494A, 0x5354, 0x5556, 0x5758,
0x595A, 0x6364, 0x6566, 0x6768, 0x696A, 0x7374, 0x7576, 0x7778, 0x797A,
0x8283, 0x8485, 0x8687, 0x8889, 0x8A92, 0x9394, 0x9596, 0x9798, 0x999A,
0xA2A3, 0xA4A5, 0xA6A7, 0xA8A9, 0xAAB2, 0xB3B4, 0xB5B6, 0xB7B8, 0xB9BA,
0xC2C3, 0xC4C5, 0xC6C7, 0xC8C9, 0xCAD2, 0xD3D4, 0xD5D6, 0xD7D8, 0xD9DA,
0xE2E3, 0xE4E5, 0xE6E7, 0xE8E9, 0xEAF2, 0xF3F4, 0xF5F6, 0xF7F8, 0xF9FA
};
UINT8 bitsize[] = {
0, 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6,
6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
8, 8
};
UINT16 luminance_dc_code_table[] = {
0x0000, 0x0002, 0x0003, 0x0004, 0x0005, 0x0006, 0x000E, 0x001E, 0x003E,
0x007E, 0x00FE, 0x01FE
};
UINT16 luminance_dc_size_table[] = {
0x0002, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0004, 0x0005, 0x0006,
0x0007, 0x0008, 0x0009
};
UINT16 chrominance_dc_code_table[] = {
0x0000, 0x0001, 0x0002, 0x0006, 0x000E, 0x001E, 0x003E, 0x007E, 0x00FE,
0x01FE, 0x03FE, 0x07FE
};
UINT16 chrominance_dc_size_table[] = {
0x0002, 0x0002, 0x0002, 0x0003, 0x0004, 0x0005, 0x0006, 0x0007, 0x0008,
0x0009, 0x000A, 0x000B
};
UINT16 luminance_ac_code_table[] = {
0x000A, 0x0000, 0x0001, 0x0004, 0x000B, 0x001A, 0x0078, 0x00F8, 0x03F6,
0xFF82, 0xFF83, 0x000C, 0x001B, 0x0079, 0x01F6, 0x07F6, 0xFF84, 0xFF85,
0xFF86, 0xFF87, 0xFF88, 0x001C, 0x00F9, 0x03F7, 0x0FF4, 0xFF89, 0xFF8A,
0xFF8b, 0xFF8C, 0xFF8D, 0xFF8E, 0x003A, 0x01F7, 0x0FF5, 0xFF8F, 0xFF90,
0xFF91, 0xFF92, 0xFF93, 0xFF94, 0xFF95, 0x003B, 0x03F8, 0xFF96, 0xFF97,
0xFF98, 0xFF99, 0xFF9A, 0xFF9B, 0xFF9C, 0xFF9D, 0x007A, 0x07F7, 0xFF9E,
0xFF9F, 0xFFA0, 0xFFA1, 0xFFA2, 0xFFA3, 0xFFA4, 0xFFA5, 0x007B, 0x0FF6,
0xFFA6, 0xFFA7, 0xFFA8, 0xFFA9, 0xFFAA, 0xFFAB, 0xFFAC, 0xFFAD, 0x00FA,
0x0FF7, 0xFFAE, 0xFFAF, 0xFFB0, 0xFFB1, 0xFFB2, 0xFFB3, 0xFFB4, 0xFFB5,
0x01F8, 0x7FC0, 0xFFB6, 0xFFB7, 0xFFB8, 0xFFB9, 0xFFBA, 0xFFBB, 0xFFBC,
0xFFBD, 0x01F9, 0xFFBE, 0xFFBF, 0xFFC0, 0xFFC1, 0xFFC2, 0xFFC3, 0xFFC4,
0xFFC5, 0xFFC6, 0x01FA, 0xFFC7, 0xFFC8, 0xFFC9, 0xFFCA, 0xFFCB, 0xFFCC,
0xFFCD, 0xFFCE, 0xFFCF, 0x03F9, 0xFFD0, 0xFFD1, 0xFFD2, 0xFFD3, 0xFFD4,
0xFFD5, 0xFFD6, 0xFFD7, 0xFFD8, 0x03FA, 0xFFD9, 0xFFDA, 0xFFDB, 0xFFDC,
0xFFDD, 0xFFDE, 0xFFDF, 0xFFE0, 0xFFE1, 0x07F8, 0xFFE2, 0xFFE3, 0xFFE4,
0xFFE5, 0xFFE6, 0xFFE7, 0xFFE8, 0xFFE9, 0xFFEA, 0xFFEB, 0xFFEC, 0xFFED,
0xFFEE, 0xFFEF, 0xFFF0, 0xFFF1, 0xFFF2, 0xFFF3, 0xFFF4, 0xFFF5, 0xFFF6,
0xFFF7, 0xFFF8, 0xFFF9, 0xFFFA, 0xFFFB, 0xFFFC, 0xFFFD, 0xFFFE,
0x07F9
};
UINT16 luminance_ac_size_table[] = {
0x0004, 0x0002, 0x0002, 0x0003, 0x0004, 0x0005, 0x0007, 0x0008, 0x000A,
0x0010, 0x0010, 0x0004, 0x0005, 0x0007, 0x0009, 0x000B, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0005, 0x0008, 0x000A, 0x000C, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x0009, 0x000C, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x000A, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007, 0x000B, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007, 0x000C,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0008,
0x000C, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0009, 0x000F, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x000A, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x000A, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000B, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x000B
};
UINT16 chrominance_ac_code_table[] = {
0x0000, 0x0001, 0x0004, 0x000A, 0x0018, 0x0019, 0x0038, 0x0078, 0x01F4,
0x03F6, 0x0FF4, 0x000B, 0x0039, 0x00F6, 0x01F5, 0x07F6, 0x0FF5, 0xFF88,
0xFF89, 0xFF8A, 0xFF8B, 0x001A, 0x00F7, 0x03F7, 0x0FF6, 0x7FC2, 0xFF8C,
0xFF8D, 0xFF8E, 0xFF8F, 0xFF90, 0x001B, 0x00F8, 0x03F8, 0x0FF7, 0xFF91,
0xFF92, 0xFF93, 0xFF94, 0xFF95, 0xFF96, 0x003A, 0x01F6, 0xFF97, 0xFF98,
0xFF99, 0xFF9A, 0xFF9B, 0xFF9C, 0xFF9D, 0xFF9E, 0x003B, 0x03F9, 0xFF9F,
0xFFA0, 0xFFA1, 0xFFA2, 0xFFA3, 0xFFA4, 0xFFA5, 0xFFA6, 0x0079, 0x07F7,
0xFFA7, 0xFFA8, 0xFFA9, 0xFFAA, 0xFFAB, 0xFFAC, 0xFFAD, 0xFFAE, 0x007A,
0x07F8, 0xFFAF, 0xFFB0, 0xFFB1, 0xFFB2, 0xFFB3, 0xFFB4, 0xFFB5, 0xFFB6,
0x00F9, 0xFFB7, 0xFFB8, 0xFFB9, 0xFFBA, 0xFFBB, 0xFFBC, 0xFFBD, 0xFFBE,
0xFFBF, 0x01F7, 0xFFC0, 0xFFC1, 0xFFC2, 0xFFC3, 0xFFC4, 0xFFC5, 0xFFC6,
0xFFC7, 0xFFC8, 0x01F8, 0xFFC9, 0xFFCA, 0xFFCB, 0xFFCC, 0xFFCD, 0xFFCE,
0xFFCF, 0xFFD0, 0xFFD1, 0x01F9, 0xFFD2, 0xFFD3, 0xFFD4, 0xFFD5, 0xFFD6,
0xFFD7, 0xFFD8, 0xFFD9, 0xFFDA, 0x01FA, 0xFFDB, 0xFFDC, 0xFFDD, 0xFFDE,
0xFFDF, 0xFFE0, 0xFFE1, 0xFFE2, 0xFFE3, 0x07F9, 0xFFE4, 0xFFE5, 0xFFE6,
0xFFE7, 0xFFE8, 0xFFE9, 0xFFEA, 0xFFEb, 0xFFEC, 0x3FE0, 0xFFED, 0xFFEE,
0xFFEF, 0xFFF0, 0xFFF1, 0xFFF2, 0xFFF3, 0xFFF4, 0xFFF5, 0x7FC3, 0xFFF6,
0xFFF7, 0xFFF8, 0xFFF9, 0xFFFA, 0xFFFB, 0xFFFC, 0xFFFD, 0xFFFE,
0x03FA
};
UINT16 chrominance_ac_size_table[] = {
0x0002, 0x0002, 0x0003, 0x0004, 0x0005, 0x0005, 0x0006, 0x0007, 0x0009,
0x000A, 0x000C, 0x0004, 0x0006, 0x0008, 0x0009, 0x000B, 0x000C, 0x0010,
0x0010, 0x0010, 0x0010, 0x0005, 0x0008, 0x000A, 0x000C, 0x000F, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0005, 0x0008, 0x000A, 0x000C, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x0009, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0006, 0x000A, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007, 0x000B,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0007,
0x000B, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0008, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0009, 0x0010, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000B, 0x0010, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000E, 0x0010, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x000F, 0x0010,
0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010, 0x0010,
0x000A
};
#define FOUR_ZERO_ZERO 0 // Grey scale Y00 ...
#define FOUR_TWO_TWO 2 //Y00 Cb Y01 Cr
static void (*read_format) (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT8 * input_ptr);
static void RGB_2_400 (UINT8 * input_ptr, UINT8 * output_ptr, UINT32 image_width,UINT32 image_height)
{
UINT32 i, size;
UINT8 R, G, B;
INT32 Y;
UINT8 * inbuf = input_ptr;
size = image_width * image_height;
for (i = size; i > 0; i--)
{
B = inbuf[0];
G = inbuf[1];
R = inbuf[2];
inbuf += 3;
Y = CLIP ((77 * R + 150 * G + 29 * B) >> 8);
*input_ptr++ = (UINT8) Y;
}
}
static void read_400_format (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT8 * input_ptr)
{
INT32 i, j;
INT16 * Y1_Ptr = Y1;
UINT16 rows = jpeg_encoder_structure->rows;
UINT16 cols = jpeg_encoder_structure->cols;
UINT16 incr = jpeg_encoder_structure->incr;
for (i = rows; i > 0; i--)
{
for (j = cols; j > 0; j--)
*Y1_Ptr++ = *input_ptr++ - 128;
for (j = 8 - cols; j > 0; j--)
*Y1_Ptr++ = *(Y1_Ptr - 1);
input_ptr += incr;
}
for (i = 8 - rows; i > 0; i--)
{
for (j = 8; j > 0; j--)
*Y1_Ptr++ = *(Y1_Ptr - 8);
}
}
static void
RGB_2_422 (UINT8 * input_ptr, UINT8 * output_ptr, UINT32 image_width,
UINT32 image_height)
{
UINT32 i, size;
UINT8 R, G, B, R1, G1, B1;
INT32 Y, Yp, Cb, Cr;
UINT8 * inbuf = input_ptr;
size = image_width * image_height/2;
for (i = size; i > 0; i--)
{
B = inbuf[0];
G = inbuf[1];
R = inbuf[2];
B1 = inbuf[3];
G1 = inbuf[4];
R1 = inbuf[5];
inbuf += 6;
Y = CLIP ((77 * R + 150 * G + 29 * B) >> 8);
Yp = CLIP ((77 * R1 + 150 * G1 + 29 * B1) >> 8);
Cb = CLIP (((-43 * R - 85 * G + 128 * B) >> 8) + 128);
Cr = CLIP (((128 * R - 107 * G - 21 * B) >> 8) + 128);
*input_ptr++ = (UINT8) Y;
*input_ptr++ = (UINT8) Cb;
*input_ptr++ = (UINT8) Yp;
*input_ptr++ = (UINT8) Cr;
}
}
/* This function implements 16 Step division for Q.15 format data */
UINT16 DSP_Division (UINT32 numer, UINT32 denom)
{
UINT16 i;
denom <<= 15;
for (i = 16; i > 0; i--)
{
if (numer > denom)
{
numer -= denom;
numer <<= 1;
numer++;
}
else
numer <<= 1;
}
return (UINT16) numer;
}
UINT8 * write_markers (UINT8 * output_ptr, UINT32 image_format,UINT32 image_width, UINT32 image_height)
{
UINT16 i, header_length;
UINT8 number_of_components;
// Start of image marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xD8;
// Quantization table marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xDB;
// Quantization table length
*output_ptr++ = 0x00;
*output_ptr++ = 0x84;
// Pq, Tq
*output_ptr++ = 0x00;
// Lqt table
for (i = 0; i < 64; i++)
*output_ptr++ = Lqt[i];
// Pq, Tq
*output_ptr++ = 0x01;
// Cqt table
for (i = 0; i < 64; i++)
*output_ptr++ = Cqt[i];
// huffman table(DHT)
for (i = 0; i < 210; i++)
{
*output_ptr++ = (UINT8) (markerdata[i] >> 8);
*output_ptr++ = (UINT8) markerdata[i];
}
//FOUR_ZERO_ZERO
number_of_components = 3;
// Frame header(SOF)
// Start of frame marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xC0;
header_length = (UINT16) (8 + 3 * number_of_components);
// Frame header length
*output_ptr++ = (UINT8) (header_length >> 8);
*output_ptr++ = (UINT8) header_length;
// Precision (P)
*output_ptr++ = 0x08;
// image height
*output_ptr++ = (UINT8) (image_height >> 8);
*output_ptr++ = (UINT8) image_height;
// image width
*output_ptr++ = (UINT8) (image_width >> 8);
*output_ptr++ = (UINT8) image_width;
// Nf
*output_ptr++ = number_of_components;
// FOUR_ZERO_ZERO
*output_ptr++ = 0x01;
*output_ptr++ = 0x21;
*output_ptr++ = 0x00;
*output_ptr++ = 0x02;
*output_ptr++ = 0x11;
*output_ptr++ = 0x01;
*output_ptr++ = 0x03;
*output_ptr++ = 0x11;
*output_ptr++ = 0x01;
// Scan header(SOF)
// Start of scan marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xDA;
header_length = (UINT16) (6 + (number_of_components << 1));
// Scan header length
*output_ptr++ = (UINT8) (header_length >> 8);
*output_ptr++ = (UINT8) header_length;
// Ns
*output_ptr++ = number_of_components;
//FOUR_ZERO_ZERO
*output_ptr++ = 0x01;
*output_ptr++ = 0x00;
*output_ptr++ = 0x02;
*output_ptr++ = 0x11;
*output_ptr++ = 0x03;
*output_ptr++ = 0x11;
*output_ptr++ = 0x00;
*output_ptr++ = 0x3F;
*output_ptr++ = 0x00;
return output_ptr;
}
static void
read_422_format (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,
UINT8 * input_ptr)
{
INT32 i, j;
UINT16 Y1_cols, Y2_cols;
INT16 * Y1_Ptr = Y1;
INT16 * Y2_Ptr = Y2;
INT16 * CB_Ptr = CB;
INT16 * CR_Ptr = CR;
UINT16 rows = jpeg_encoder_structure->rows;
UINT16 cols = jpeg_encoder_structure->cols;
UINT16 incr = jpeg_encoder_structure->incr;
if (cols <= 8)
{
Y1_cols = cols;
Y2_cols = 0;
}
else
{
Y1_cols = 8;
Y2_cols = (UINT16) (cols - 8);
}
for (i = rows; i > 0; i--)
{
for (j = Y1_cols >> 1; j > 0; j--)
{
*Y1_Ptr++ = *input_ptr++ - 128;
*CB_Ptr++ = *input_ptr++ - 128;
*Y1_Ptr++ = *input_ptr++ - 128;
*CR_Ptr++ = *input_ptr++ - 128;
}
for (j = Y2_cols >> 1; j > 0; j--)
{
*Y2_Ptr++ = *input_ptr++ - 128;
*CB_Ptr++ = *input_ptr++ - 128;
*Y2_Ptr++ = *input_ptr++ - 128;
*CR_Ptr++ = *input_ptr++ - 128;
}
if (cols <= 8)
{
for (j = 8 - Y1_cols; j > 0; j--)
*Y1_Ptr++ = *(Y1_Ptr - 1);
for (j = 8 - Y2_cols; j > 0; j--)
*Y2_Ptr++ = *(Y1_Ptr - 1);
}
else
{
for (j = 8 - Y2_cols; j > 0; j--)
*Y2_Ptr++ = *(Y2_Ptr - 1);
}
for (j = (16 - cols) >> 1; j > 0; j--)
{
*CB_Ptr++ = *(CB_Ptr - 1);
*CR_Ptr++ = *(CR_Ptr - 1);
}
input_ptr += incr;
}
for (i = 8 - rows; i > 0; i--)
{
for (j = 8; j > 0; j--)
{
*Y1_Ptr++ = *(Y1_Ptr - 8);
*Y2_Ptr++ = *(Y2_Ptr - 8);
*CB_Ptr++ = *(CB_Ptr - 8);
*CR_Ptr++ = *(CR_Ptr - 8);
}
}
}
static void initialization (JPEG_ENCODER_STRUCTURE * jpeg, UINT32 image_format,UINT32 image_width, UINT32 image_height)
{
UINT16 mcu_width, mcu_height, bytes_per_pixel;
lcode = 0;
bitindex = 0;
jpeg->mcu_width = mcu_width = 16;
jpeg->horizontal_mcus = (UINT16) ((image_width + mcu_width - 1) >> 4);
jpeg->mcu_height = mcu_height = 8;
jpeg->vertical_mcus =
(UINT16) ((image_height + mcu_height - 1) >> 3);
bytes_per_pixel = 2;
read_format = read_422_format;
jpeg->rows_in_bottom_mcus =(UINT16) (image_height - (jpeg->vertical_mcus - 1) * mcu_height);
jpeg->cols_in_right_mcus =(UINT16) (image_width - (jpeg->horizontal_mcus - 1) * mcu_width);
jpeg->length_minus_mcu_width =(UINT16) ((image_width - mcu_width) * bytes_per_pixel);
jpeg->length_minus_width =(UINT16) ((image_width - jpeg->cols_in_right_mcus) * bytes_per_pixel);
jpeg->mcu_width_size = (UINT16) (mcu_width * bytes_per_pixel);
jpeg->offset =(UINT16) ((image_width * (mcu_height - 1) -(mcu_width - jpeg->cols_in_right_mcus)) * bytes_per_pixel);
jpeg->ldc1 = 0;
jpeg->ldc2 = 0;
jpeg->ldc3 = 0;
}
/* Multiply Quantization table with quality factor to get LQT and CQT */
void initialize_quantization_tables (UINT32 quality_factor)
{
UINT16 i, index;
UINT32 value;
UINT8 luminance_quant_table[] = {
16, 11, 10, 16, 24, 40, 51, 61, 12, 12, 14, 19, 26, 58, 60, 55, 14, 13,
16, 24, 40, 57, 69, 56, 14, 17, 22, 29, 51, 87, 80, 62, 18, 22, 37,
56, 68, 109, 103, 77, 24, 35, 55, 64, 81, 104, 113, 92, 49, 64, 78,
87, 103, 121, 120, 101, 72, 92, 95, 98, 112, 100, 103, 99
};
UINT8 chrominance_quant_table[] = {
17, 18, 24, 47, 99, 99, 99, 99, 18, 21, 26, 66, 99, 99, 99, 99, 24, 26,
56, 99, 99, 99, 99, 99, 47, 66, 99, 99, 99, 99, 99, 99, 99, 99, 99,
99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99,
99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99
};
for (i = 0; i < 64; i++)
{
index = zigzag_table[i];
/* luminance quantization table * quality factor */
value = luminance_quant_table[i] * quality_factor;
value = (value + 0x200) >> 10;
if (value == 0)
value = 1;
else if (value > 255)
value = 255;
Lqt[index] = (UINT8) value;
ILqt[i] = DSP_Division (0x8000, value);
/* chrominance quantization table * quality factor */
value = chrominance_quant_table[i] * quality_factor;
value = (value + 0x200) >> 10;
if (value == 0)
value = 1;
else if (value > 255)
value = 255;
Cqt[index] = (UINT8) value;
ICqt[i] = DSP_Division (0x8000, value);
}
}
UINT8 * close_bitstream (UINT8 * output_ptr)
{
UINT16 i, count;
UINT8 * ptr;
if (bitindex > 0)
{
lcode <<= (32 - bitindex);
count = (bitindex + 7) >> 3;
ptr = (UINT8 *) & lcode + 3;
for (i = count; i > 0; i--)
{
if ((*output_ptr++ = *ptr--) == 0xff)
*output_ptr++ = 0;
}
}
// End of image marker
*output_ptr++ = 0xFF;
*output_ptr++ = 0xD9;
return output_ptr;
}
/* DCT for One block(8x8) */
static void DCT (INT16 * data)
{
UINT16 i;
INT32 x0, x1, x2, x3, x4, x5, x6, x7, x8;
/* All values are shifted left by 10
and rounded off to nearest integer */
static const UINT16 c1 = 1420; /* cos PI/16 * root(2) */
static const UINT16 c2 = 1338; /* cos PI/8 * root(2) */
static const UINT16 c3 = 1204; /* cos 3PI/16 * root(2) */
static const UINT16 c5 = 805; /* cos 5PI/16 * root(2) */
static const UINT16 c6 = 554; /* cos 3PI/8 * root(2) */
static const UINT16 c7 = 283; /* cos 7PI/16 * root(2) */
static const UINT16 s1 = 3;
static const UINT16 s2 = 10;
static const UINT16 s3 = 13;
for (i = 8; i > 0; i--)
{
x8 = data[0] + data[7];
x0 = data[0] - data[7];
x7 = data[1] + data[6];
x1 = data[1] - data[6];
x6 = data[2] + data[5];
x2 = data[2] - data[5];
x5 = data[3] + data[4];
x3 = data[3] - data[4];
x4 = x8 + x5;
x8 -= x5;
x5 = x7 + x6;
x7 -= x6;
data[0] = (INT16) (x4 + x5);
data[4] = (INT16) (x4 - x5);
data[2] = (INT16) ((x8 * c2 + x7 * c6) >> s2);
data[6] = (INT16) ((x8 * c6 - x7 * c2) >> s2);
data[7] = (INT16) ((x0 * c7 - x1 * c5 + x2 * c3 - x3 * c1) >> s2);
data[5] = (INT16) ((x0 * c5 - x1 * c1 + x2 * c7 + x3 * c3) >> s2);
data[3] = (INT16) ((x0 * c3 - x1 * c7 - x2 * c1 - x3 * c5) >> s2);
data[1] = (INT16) ((x0 * c1 + x1 * c3 + x2 * c5 + x3 * c7) >> s2);
data += 8;
}
data -= 64;
for (i = 8; i > 0; i--)
{
x8 = data[0] + data[56];
x0 = data[0] - data[56];
x7 = data[8] + data[48];
x1 = data[8] - data[48];
x6 = data[16] + data[40];
x2 = data[16] - data[40];
x5 = data[24] + data[32];
x3 = data[24] - data[32];
x4 = x8 + x5;
x8 -= x5;
x5 = x7 + x6;
x7 -= x6;
data[0] = (INT16) ((x4 + x5) >> s1);
data[32] = (INT16) ((x4 - x5) >> s1);
data[16] = (INT16) ((x8 * c2 + x7 * c6) >> s3);
data[48] = (INT16) ((x8 * c6 - x7 * c2) >> s3);
data[56] = (INT16) ((x0 * c7 - x1 * c5 + x2 * c3 - x3 * c1) >> s3);
data[40] = (INT16) ((x0 * c5 - x1 * c1 + x2 * c7 + x3 * c3) >> s3);
data[24] = (INT16) ((x0 * c3 - x1 * c7 - x2 * c1 - x3 * c5) >> s3);
data[8] = (INT16) ((x0 * c1 + x1 * c3 + x2 * c5 + x3 * c7) >> s3);
data++;
}
}
/* multiply DCT Coefficients with Quantization table and store in ZigZag location */
void quantization (INT16 * const data, UINT16 * const quant_table_ptr)
{
INT16 i;
INT32 value;
for (i = 63; i >= 0; i--)
{
value = data[i] * quant_table_ptr[i];
value = (value + 0x4000) >> 15;
Temp[zigzag_table[i]] = (INT16) value;
}
}
UINT8 * huffman (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT16 component, UINT8 * output_ptr)
{
UINT16 i;
UINT16 * DcCodeTable, *DcSizeTable, *AcCodeTable, *AcSizeTable;
INT16 * Temp_Ptr, Coeff, LastDc;
UINT16 AbsCoeff, HuffCode, HuffSize, RunLength = 0, DataSize = 0, index;
INT16 bits_in_next_word;
UINT16 numbits;
UINT32 data;
Temp_Ptr = Temp;
Coeff = *Temp_Ptr++;
if (component == 1)
{
DcCodeTable = luminance_dc_code_table;
DcSizeTable = luminance_dc_size_table;
AcCodeTable = luminance_ac_code_table;
AcSizeTable = luminance_ac_size_table;
LastDc = jpeg_encoder_structure->ldc1;
jpeg_encoder_structure->ldc1 = Coeff;
}
else
{
DcCodeTable = chrominance_dc_code_table;
DcSizeTable = chrominance_dc_size_table;
AcCodeTable = chrominance_ac_code_table;
AcSizeTable = chrominance_ac_size_table;
if (component == 2)
{
LastDc = jpeg_encoder_structure->ldc2;
jpeg_encoder_structure->ldc2 = Coeff;
}
else
{
LastDc = jpeg_encoder_structure->ldc3;
jpeg_encoder_structure->ldc3 = Coeff;
}
}
Coeff -= LastDc;
AbsCoeff = (Coeff < 0) ? -Coeff-- : Coeff;
while (AbsCoeff != 0)
{
AbsCoeff >>= 1;
DataSize++;
}
HuffCode = DcCodeTable[DataSize];
HuffSize = DcSizeTable[DataSize];
Coeff &= (1 << DataSize) - 1;
data = (HuffCode << DataSize) | Coeff;
numbits = HuffSize + DataSize;
//PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
for (i = 63; i > 0; i--)
{
if ((Coeff = *Temp_Ptr++) != 0)
{
while (RunLength > 15)
{
RunLength -= 16;
data = AcCodeTable[161];
numbits = AcSizeTable[161];
//PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
}
AbsCoeff = (Coeff < 0) ? -Coeff-- : Coeff;
if (AbsCoeff >> 8 == 0)
DataSize = bitsize[AbsCoeff];
else
DataSize = bitsize[AbsCoeff >> 8] + 8;
index = RunLength * 10 + DataSize;
HuffCode = AcCodeTable[index];
HuffSize = AcSizeTable[index];
Coeff &= (1 << DataSize) - 1;
data = (HuffCode << DataSize) | Coeff;
numbits = HuffSize + DataSize;
// PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
RunLength = 0;
}
else
RunLength++;
}
if (RunLength != 0)
{
data = AcCodeTable[0];
numbits = AcSizeTable[0];
// PUTBITS
{
bits_in_next_word = (INT16) (bitindex + numbits - 32);
if (bits_in_next_word < 0)
{
lcode = (lcode << numbits) | data;
bitindex += numbits;
}
else
{
lcode = (lcode << (32 - bitindex)) | (data >> bits_in_next_word);
if ((*output_ptr++ = (UINT8) (lcode >> 24)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 16)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) (lcode >> 8)) == 0xff)
*output_ptr++ = 0;
if ((*output_ptr++ = (UINT8) lcode) == 0xff)
*output_ptr++ = 0;
lcode = data;
bitindex = bits_in_next_word;
}
}
}
return output_ptr;
}
static UINT8 *encodeMCU (JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure,UINT32 image_format, UINT8 * output_ptr)
{
DCT (Y1);
quantization (Y1, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
if (image_format == FOUR_ZERO_ZERO)
return output_ptr;
DCT (Y2);
quantization (Y2, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
if (image_format == FOUR_TWO_TWO)
goto chroma;
DCT (Y3);
quantization (Y3, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
DCT (Y4);
quantization (Y4, ILqt);
output_ptr = huffman (jpeg_encoder_structure, 1, output_ptr);
chroma:DCT (CB);
quantization (CB, ICqt);
output_ptr = huffman (jpeg_encoder_structure, 2, output_ptr);
DCT (CR);
quantization (CR, ICqt);
output_ptr = huffman (jpeg_encoder_structure, 3, output_ptr);
return output_ptr;
}
UINT32 encode_image (UINT8 * input_ptr, UINT8 * output_ptr,UINT32 quality_factor, UINT32 image_format,UINT32 image_width, UINT32 image_height)
{
UINT16 i, j;
UINT8 * output;
JPEG_ENCODER_STRUCTURE JpegStruct;
JPEG_ENCODER_STRUCTURE * jpeg_encoder_structure = &JpegStruct;
output = output_ptr;
image_format = FOUR_TWO_TWO;
RGB_2_422 (input_ptr, output_ptr, image_width, image_height);
/* Initialization of JPEG control structure */
initialization (jpeg_encoder_structure, image_format, image_width,
image_height);
/* Quantization Table Initialization */
initialize_quantization_tables (quality_factor);
/* Writing Marker Data */
output_ptr =
write_markers (output_ptr, image_format, image_width, image_height);
for (i = 1; i <= jpeg_encoder_structure->vertical_mcus; i++)
{
if (i < jpeg_encoder_structure->vertical_mcus)
jpeg_encoder_structure->rows = jpeg_encoder_structure->mcu_height;
else
jpeg_encoder_structure->rows =
jpeg_encoder_structure->rows_in_bottom_mcus;
for (j = 1; j <= jpeg_encoder_structure->horizontal_mcus; j++)
{
if (j < jpeg_encoder_structure->horizontal_mcus)
{
jpeg_encoder_structure->cols =
jpeg_encoder_structure->mcu_width;
jpeg_encoder_structure->incr =
jpeg_encoder_structure->length_minus_mcu_width;
}
else
{
jpeg_encoder_structure->cols =
jpeg_encoder_structure->cols_in_right_mcus;
jpeg_encoder_structure->incr =
jpeg_encoder_structure->length_minus_width;
}
read_format (jpeg_encoder_structure, input_ptr);
/* Encode the data in MCU */
output_ptr =encodeMCU (jpeg_encoder_structure, image_format, output_ptr);
input_ptr += jpeg_encoder_structure->mcu_width_size;
}
input_ptr += jpeg_encoder_structure->offset;
}
/* Close Routine */
output_ptr = close_bitstream (output_ptr);
return (UINT32) (output_ptr - output);
}
void main()
{
int size=1280*720*3,sizeout=0,x,y;
unsigned char *outpict ;
outpict = (unsigned char *) malloc (size);
unsigned char *inpict ;
inpict = (unsigned char *) malloc (size);
for(x=0;x<1280;x++)//inpict初始化成RGB模式
for (y=0;y<720;y++)
{
*(inpict+x*3+y*1280*3)=y%2;
*(inpict+x*3+y*1280*3+1)=x%2;
*(inpict+x*3+y*1280*3+2)=y%2;
}
sizeout = encode_image (inpict,outpict,1,0,1280,720);
free(outpict);
free(inpict);
return 1;
}
内容来自用户分享和网络整理，不保证内容的准确性，如有侵权内容，可联系管理员处理 
相关文章推荐
- JPEG图像压缩解压算法——C++实现
- 二进制图像压缩算法，2005IEEE一篇文章的实现及简单改进算法
- 关于用java io实现文件压缩与解压（不涉及压缩算法）
- Android图片压缩加密上传 - JPEG压缩算法解析
- java实现linux中gzip压缩解压缩算法：byte[]字节数组，文件，字符串，数据流的压缩解压缩
- c#实现linux中gzip压缩解压缩算法：byte[]字节数组，文件，字符串，数据流的压缩解压缩
- 数据的压缩存储与解压缩算法实现(C语言)
- SDT 旋转门压缩算法的C语言实现 压缩和解压缩
- 一种整数数据压缩存储的算法实现
- 轨迹压缩之Douglas-Peucker算法之C++实现
- JPEG压缩图像超分辨率重建算法
- JPEG图像压缩算法流程详解（转）
- URL短地址压缩算法 微博短地址原理解析 （Java实现）
- java笔试题：用java实现字符串压缩算法
- 数据的压缩存储与解压缩算法实现(C语言)
- Fast Compressive Tracking（快速压缩跟踪）算法的C++代码实现
- java实现的LZW 压缩算法源码
- union-find算法——使用路径压缩的加权quick-union算法实现
- 实现简易字符串压缩算法：由字母a-z或者A-Z组成，将其中连续出现2次以上（含2次）的字母转换为字母和出现次数，
- JPEG图像压缩算法流程详解