基于C语言的JPEG编码代码详解

最新推荐文章于 2024-07-22 23:58:15 发布
最新推荐文章于 2024-07-22 23:58:15 发布
阅读量5k 收藏 57
点赞数 3
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_37633745/article/details/107960668
版权

一、基于C语言的JPEG编码代码详解

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#pragma warning(disable : 4996)

static struct APP0infotype {
	unsigned short int marker;           // = 0xFFE0
	unsigned short int length;           // = 16 for usual JPEG, no thumbnail
	unsigned char      JFIFsignature[5]; // = "JFIF",'\0'
	unsigned char      versionhi;        // 1
	unsigned char      versionlo;        // 1
	unsigned char      xyunits;          // 0 = no units, normal density
	unsigned short int xdensity;         // 1
	unsigned short int ydensity;         // 1
	unsigned char      thumbnwidth;      // 0
	unsigned char      thumbnheight;     // 0
} APP0info = { 0xFFE0,16,'J','F','I','F',0,1,1,0,1,1,0,0 };

static struct  SOF0infotype {
	unsigned short int marker;         // = 0xFFC0
	unsigned short int length;         // = 17 for a truecolor YCbCr JPG
	unsigned char      precision;     // Should be 8: 8 bits/sample
	unsigned short int height;
	unsigned short int width;
	unsigned char      nrofcomponents; //Should be 3: We encode a truecolor JPG
	unsigned char      IdY;            // = 1
	unsigned char      HVY;            // sampling factors for Y
									   // (bit 0-3 vert., 4-7 hor.)
	unsigned char      QTY;            // Quantization Table number for Y = 0
	unsigned char      IdCb;           // = 2
	unsigned char      HVCb;
	unsigned char      QTCb;           // 1
	unsigned char      IdCr;           // = 3
	unsigned char      HVCr;
	unsigned char      QTCr;           // Normally equal to QTCb = 1
} SOF0info = { 0xFFC0,17,8,0,0,3,1,0x22,0,2,0x11,1,3,0x11,1 };

// Default sampling factors are 1,1 for every image component: No downsampling

static struct DQTinfotype {
	unsigned short int marker;      // = 0xFFDB（标记代码）
	unsigned short int length;      // = 132（除标记代码后的数据总长度）

	unsigned char      QTYinfo;     //（高四位量化表精度0：8位，1：16位；低四位量化表ID：0~3） 
	unsigned char      Ytable[64];

	unsigned char      QTCbinfo;    // = 1 (quantization table for Cb,Cr}
	unsigned char      Cbtable[64];
} DQTinfo;

// Ytable from DQTinfo should be equal to a scaled and zizag reordered version
// of the table which can be found in "tables.h": std_luminance_qt
// Cbtable , similar = std_chrominance_qt
// We'll init them in the program using set_DQTinfo function

static struct DHTinfotype {
	unsigned short int marker;           // = 0xFFC4（两字节的标记代码）
	unsigned short int length;           //0x01A2（除标记代码后的数据总长度）

	unsigned char      HTYDCinfo;        //bit 0..3: number of HT (0..3), for Y =0（表ID和表类型：高四位表类型0直流1AC交流；低四位huffman表ID）
										 //bit 4  :type of HT, 0 = DC table,1 = AC table
										 //bit 5..7: not used, must be 0
	unsigned char      YDC_nrcodes[16];  //at index i = nr of codes with length i（不同位数的码字数量）
	unsigned char      YDC_values[12];	 //

	unsigned char      HTYACinfo;        // = 0x10
	unsigned char      YAC_nrcodes[16];
	unsigned char      YAC_values[162];  //we'll use the standard Huffman tables

	unsigned char      HTCbDCinfo;       // = 1
	unsigned char      CbDC_nrcodes[16];
	unsigned char      CbDC_values[12];

	unsigned char      HTCbACinfo;       // = 0x11
	unsigned char      CbAC_nrcodes[16];
	unsigned char      CbAC_values[162];
} DHTinfo;

static struct SOSinfotype {
	unsigned short int marker;    		// = 0xFFDA
	unsigned short int length;    		// = 12
	unsigned char      nrofcomponents; 	// Should be 3: truecolor JPG

	unsigned char      IdY;            	//1
	unsigned char      HTY;            	//0 // bits 0..3: AC table (0..3)
										  // bits 4..7: DC table (0..3)

	unsigned char      IdCb;           	//2
	unsigned char      HTCb;           	//0x11

	unsigned char      IdCr;           	//3
	unsigned char      HTCr;           	//0x11

	unsigned char      Ss, Se, Bf;       	// not interesting, they should be 0,63,0
} SOSinfo = { 0xFFDA,12,3,1,0,2,0x11,3,0x11,0,0x3F,0 };

typedef struct { unsigned char length; unsigned short int value; } bitstring;

#define  Y(R,G,B) ((unsigned char)( (YRtab[(R)] + YGtab[(G)] + YBtab[(B)]) >> 16 ) - 128)
#define Cb(R,G,B) ((unsigned char)( (CbRtab[(R)] + CbGtab[(G)] + CbBtab[(B)]) >> 16 ))
#define Cr(R,G,B) ((unsigned char)( (CrRtab[(R)] + CrGtab[(G)] + CrBtab[(B)]) >> 16 ))

//fputc((b),fp_jpeg_stream) 把参数 b 指定的字符（一个无符号字符）写入到指定的流 fp_jpeg_stream 中，并把位置标识符往前移动。
#define writebyte(b) fputc((b),fp_jpeg_stream)
#define writeword(w) writebyte((w)/256);writebyte((w)%256);

static unsigned char zigzag[64] = {
  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
};

/*
These are the sample quantization tables given in JPEG spec section K.1.
The spec says that the values given produce "good" quality, and
when divided by 2, "very good" quality.
*/

static unsigned char DQT_Y[64] = {
  0x08, 0x06, 0x06, 0x07, 0x06, 0x05, 0x08, 0x07,
  0x07, 0x07, 0x09, 0x09, 0x08, 0x0a, 0x0c, 0x14,
  0x0d, 0x0c, 0x0b, 0x0b, 0x0c, 0x19, 0x12, 0x13,
  0x0f, 0x14, 0x1d, 0x1a, 0x1f, 0x1e, 0x1d, 0x1a,
  0x1c, 0x1c, 0x20, 0x24, 0x2e, 0x27, 0x20, 0x22,
  0x2c, 0x23, 0x1c, 0x1c, 0x28, 0x37, 0x29, 0x2c,
  0x30, 0x31, 0x34, 0x34, 0x34, 0x1f, 0x27, 0x39,
  0x3d, 0x38, 0x32, 0x3c, 0x2e, 0x33, 0x34, 0x32
};

static unsigned char DQT_C[64] = {
  0x09, 0x09, 0x09, 0x0c, 0x0b, 0x0c, 0x18, 0x0d,
  0x0d, 0x18, 0x32, 0x21, 0x1c, 0x21, 0x32, 0x32,
  0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32,
  0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32,
  0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32,
  0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32,
  0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32,
  0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32, 0x32
};


static unsigned char std_luminance_qt[64] = {
  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
};
static unsigned char std_chrominance_qt[64] = {
  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
};
// Standard Huffman tables (cf. JPEG standard section K.3) */

static unsigned char std_dc_luminance_nrcodes[17] = { 0,0,1,5,1,1,1,1,1,1,0,0,0,0,0,0,0 };
static unsigned char std_dc_luminance_values[12] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 };

static unsigned char std_dc_chrominance_nrcodes[17] = { 0,0,3,1,1,1,1,1,1,1,1,1,0,0,0,0,0 };
static unsigned char std_dc_chrominance_values[12] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 };

static unsigned char std_ac_luminance_nrcodes[17] = { 0,0,2,1,3,3,2,4,3,5,5,4,4,0,0,1,0x7d };
static unsigned char std_ac_luminance_values[162] = {
  0x01, 0x02, 0x03, 0x00, 0x04, 0x11, 0x05, 0x12,
  0x21, 0x31, 0x41, 0x06, 0x13, 0x51, 0x61, 0x07,
  0x22, 0x71, 0x14, 0x32, 0x81, 0x91, 0xa1, 0x08,
  0x23, 0x42, 0xb1, 0xc1, 0x15, 0x52, 0xd1, 0xf0,
  0x24, 0x33, 0x62, 0x72, 0x82, 0x09, 0x0a, 0x16,
  0x17, 0x18, 0x19, 0x1a, 0x25, 0x26, 0x27, 0x28,
  0x29, 0x2a, 0x34, 0x35, 0x36, 0x37, 0x38, 0x39,
  0x3a, 0x43, 0x44, 0x45, 0x46, 0x47, 0x48, 0x49,
  0x4a, 0x53, 0x54, 0x55, 0x56, 0x57, 0x58, 0x59,
  0x5a, 0x63, 0x64, 0x65, 0x66, 0x67, 0x68, 0x69,
  0x6a, 0x73, 0x74, 0x75, 0x76, 0x77, 0x78, 0x79,
  0x7a, 0x83, 0x84, 0x85, 0x86, 0x87, 0x88, 0x89,
  0x8a, 0x92, 0x93, 0x94, 0x95, 0x96, 0x97, 0x98,
  0x99, 0x9a, 0xa2, 0xa3, 0xa4, 0xa5, 0xa6, 0xa7,
  0xa8, 0xa9, 0xaa, 0xb2, 0xb3, 0xb4, 0xb5, 0xb6,
  0xb7, 0xb8, 0xb9, 0xba, 0xc2, 0xc3, 0xc4, 0xc5,
  0xc6, 0xc7, 0xc8, 0xc9, 0xca, 0xd2, 0xd3, 0xd4,
  0xd5, 0xd6, 0xd7, 0xd8, 0xd9, 0xda, 0xe1, 0xe2,
  0xe3, 0xe4, 0xe5, 0xe6, 0xe7, 0xe8, 0xe9, 0xea,
  0xf1, 0xf2, 0xf3, 0xf4, 0xf5, 0xf6, 0xf7, 0xf8,
  0xf9, 0xfa };

static unsigned char std_ac_chrominance_nrcodes[17] = { 0,0,2,1,2,4,4,3,4,7,5,4,4,0,1,2,0x77 };
static unsigned char std_ac_chrominance_values[162] = {
  0x00, 0x01, 0x02, 0x03, 0x11, 0x04, 0x05, 0x21,
  0x31, 0x06, 0x12, 0x41, 0x51, 0x07, 0x61, 0x71,
  0x13, 0x22, 0x32, 0x81, 0x08, 0x14, 0x42, 0x91,
  0xa1, 0xb1, 0xc1, 0x09, 0x23, 0x33, 0x52, 0xf0,
  0x15, 0x62, 0x72, 0xd1, 0x0a, 0x16, 0x24, 0x34,
  0xe1, 0x25, 0xf1, 0x17, 0x18, 0x19, 0x1a, 0x26,
  0x27, 0x28, 0x29, 0x2a, 0x35, 0x36, 0x37, 0x38,
  0x39, 0x3a, 0x43, 0x44, 0x45, 0x46, 0x47, 0x48,
  0x49, 0x4a, 0x53, 0x54, 0x55, 0x56, 0x57, 0x58,
  0x59, 0x5a, 0x63, 0x64, 0x65, 0x66, 0x67, 0x68,
  0x69, 0x6a, 0x73, 0x74, 0x75, 0x76, 0x77, 0x78,
  0x79, 0x7a, 0x82, 0x83, 0x84, 0x85, 0x86, 0x87,
  0x88, 0x89, 0x8a, 0x92, 0x93, 0x94, 0x95, 0x96,
  0x97, 0x98, 0x99, 0x9a, 0xa2, 0xa3, 0xa4, 0xa5,
  0xa6, 0xa7, 0xa8, 0xa9, 0xaa, 0xb2, 0xb3, 0xb4,
  0xb5, 0xb6, 0xb7, 0xb8, 0xb9, 0xba, 0xc2, 0xc3,
  0xc4, 0xc5, 0xc6, 0xc7, 0xc8, 0xc9, 0xca, 0xd2,
  0xd3, 0xd4, 0xd5, 0xd6, 0xd7, 0xd8, 0xd9, 0xda,
  0xe2, 0xe3, 0xe4, 0xe5, 0xe6, 0xe7, 0xe8, 0xe9,
  0xea, 0xf2, 0xf3, 0xf4, 0xf5, 0xf6, 0xf7, 0xf8,
  0xf9, 0xfa };

static unsigned char bytenew = 0; 	//The byte that will be written in the JPG file
static signed char bytepos = 7; 	//bit position in the byte we write (bytenew)
									//should be<=7 and >=0
static unsigned short int mask[16] = { 1,2,4,8,16,32,64,128,256,512,1024,2048,4096,8192,16384,32768 };

// The Huffman tables we'll use:
static bitstring YDC_HT[12];
static bitstring CbDC_HT[12];
static bitstring YAC_HT[256];
static bitstring CbAC_HT[256];

static unsigned char *category_alloc;
static unsigned char *category; 		//Here we'll keep the category of the numbers in range: -32767..32767
static bitstring *bitcode_alloc;
static bitstring *bitcode; 				//their bitcoded representation

static unsigned char *buffer; 			//image to be encoded
unsigned short int Ximage, Yimage;		//image dimensions divisible by 8

//Data Unit YCbCr
static signed char DU_Y[256];
static signed char DU_Cb[64];
static signed char DU_Cr[64];

static signed short int DU_DCT[64]; // Current DU (after DCT and quantization) which we'll zigzag

static signed short int DU[64]; 	//zigzag reordered DU which will be Huffman coded

FILE *fp_jpeg_stream;

void write_APP0info()
{
	//Nothing to overwrite for APP0info
	writeword(APP0info.marker);						//标记代码
	writeword(APP0info.length);						//数据长度
	writebyte('J'); writebyte('F'); writebyte('I'); writebyte('F'); writebyte(0);
	writebyte(APP0info.versionhi);
	writebyte(APP0info.versionlo);
	writebyte(APP0info.xyunits);
	writeword(APP0info.xdensity);
	writeword(APP0info.ydensity);
	writebyte(APP0info.thumbnwidth);
	writebyte(APP0info.thumbnheight);
}

void write_SOF0info() {
	// We should overwrite width and height
	writeword(SOF0info.marker);
	writeword(SOF0info.length);
	writebyte(SOF0info.precision);
	writeword(SOF0info.height);
	writeword(SOF0info.width);
	writebyte(SOF0info.nrofcomponents);
	writebyte(SOF0info.IdY);
	writebyte(SOF0info.HVY);
	writebyte(SOF0info.QTY);
	writebyte(SOF0info.IdCb);
	writebyte(SOF0info.HVCb);
	writebyte(SOF0info.QTCb);
	writebyte(SOF0info.IdCr);
	writebyte(SOF0info.HVCr);
	writebyte(SOF0info.QTCr);
}

void write_DQTinfo() {
	unsigned char i;
	writeword(DQTinfo.marker);
	writeword(DQTinfo.length);
	writebyte(DQTinfo.QTYinfo);
	for (i = 0; i < 64; i++) writebyte(DQTinfo.Ytable[i]);
	writebyte(DQTinfo.QTCbinfo);
	for (i = 0; i < 64; i++) writebyte(DQTinfo.Cbtable[i]);
}

void set_DQTinfo() {
	int i;
	DQTinfo.marker = 0xFFDB;	//标记段0xFFDB
	DQTinfo.length = 132;		//数据长度固定为132字节
	DQTinfo.QTYinfo = 0;		//亮度分量的精度（0），ID（0）
	DQTinfo.QTCbinfo = 1;		//色度分量的精度（0），ID（1）
	for (i = 0; i < 64; i++) 
		DQTinfo.Ytable[i] = DQT_Y[i];
	for (i = 0; i < 64; i++) 
		DQTinfo.Cbtable[i] = DQT_C[i];
}

void write_DHTinfo() {
	unsigned char i;
	writeword(DHTinfo.marker);
	writeword(DHTinfo.length);
	writebyte(DHTinfo.HTYDCinfo);
	for (i = 0; i < 16; i++)  writebyte(DHTinfo.YDC_nrcodes[i]);
	for (i = 0; i <= 11; i++) writebyte(DHTinfo.YDC_values[i]);
	writebyte(DHTinfo.HTYACinfo);
	for (i = 0; i < 16; i++)  writebyte(DHTinfo.YAC_nrcodes[i]);
	for (i = 0; i <= 161; i++) writebyte(DHTinfo.YAC_values[i]);
	writebyte(DHTinfo.HTCbDCinfo);
	for (i = 0; i < 16; i++)  writebyte(DHTinfo.CbDC_nrcodes[i]);
	for (i = 0; i <= 11; i++)  writebyte(DHTinfo.CbDC_values[i]);
	writebyte(DHTinfo.HTCbACinfo);
	for (i = 0; i < 16; i++)  writebyte(DHTinfo.CbAC_nrcodes[i]);
	for (i = 0; i <= 161; i++) writebyte(DHTinfo.CbAC_values[i]);
}

void set_DHTinfo() {
	unsigned char i;
	DHTinfo.marker = 0xFFC4;
	DHTinfo.length = 0x01A2;
	DHTinfo.HTYDCinfo = 0;												//亮度直流哈夫曼表
	for (i = 0; i < 16; i++)											//不同位数的码字数量
		DHTinfo.YDC_nrcodes[i] = std_dc_luminance_nrcodes[i + 1];
	for (i = 0; i <= 11; i++)											//编码内容
		DHTinfo.YDC_values[i] = std_dc_luminance_values[i];

	DHTinfo.HTYACinfo = 0x10;											
	for (i = 0; i < 16; i++)  
		DHTinfo.YAC_nrcodes[i] = std_ac_luminance_nrcodes[i + 1];
	for (i = 0; i <= 161; i++) 
		DHTinfo.YAC_values[i] = std_ac_luminance_values[i];

	DHTinfo.HTCbDCinfo = 1;
	for (i = 0; i < 16; i++)  
		DHTinfo.CbDC_nrcodes[i] = std_dc_chrominance_nrcodes[i + 1];
	for (i = 0; i <= 11; i++)  
		DHTinfo.CbDC_values[i] = std_dc_chrominance_values[i];

	DHTinfo.HTCbACinfo = 0x11;
	for (i = 0; i < 16; i++)  
		DHTinfo.CbAC_nrcodes[i] = std_ac_chrominance_nrcodes[i + 1];
	for (i = 0; i <= 161; i++) 
		DHTinfo.CbAC_values[i] = std_ac_chrominance_values[i];
}

void write_SOSinfo() {
	//Nothing to overwrite for SOSinfo
	writeword(SOSinfo.marker);
	writeword(SOSinfo.length);
	writebyte(SOSinfo.nrofcomponents);
	writebyte(SOSinfo.IdY);
	writebyte(SOSinfo.HTY);
	writebyte(SOSinfo.IdCb);
	writebyte(SOSinfo.HTCb);
	writebyte(SOSinfo.IdCr);
	writebyte(SOSinfo.HTCr);
	writebyte(SOSinfo.Ss);
	writebyte(SOSinfo.Se);
	writebyte(SOSinfo.Bf);
}

void writebits(bitstring bs)
{
	unsigned short int value;
	signed char posval;
	value = bs.value;				//获得码字		
	posval = bs.length - 1;			//编码的长度（1-16）
	while (posval >= 0)
	{
		if (value & mask[posval]) 
			bytenew |= mask[bytepos];
		posval--;
		bytepos--;
		if (bytepos < 0)
		{
			if (bytenew == 0xFF)
			{
				writebyte(0xFF);	//数据若为0xFF需要再写入一个0x00以作区分
				writebyte(0);
			}
			else
			{
				writebyte(bytenew);
			}
			bytepos = 7;
			bytenew = 0;
		}
	}
}

void compute_Huffman_table(unsigned char *nrcodes, unsigned char *std_table, bitstring *HT)
{
	unsigned char k, j;
	unsigned char pos_in_table;
	unsigned short int codevalue;
	codevalue = 0; pos_in_table = 0;
	for (k = 1; k <= 16; k++)
	{
		for (j = 1; j <= nrcodes[k]; j++)					//不同位数的编码数量
		{
			HT[std_table[pos_in_table]].value = codevalue;	//存储码字
			HT[std_table[pos_in_table]].length = k;			//编码的长度
			pos_in_table++;									//位置加加
			codevalue++;									//值加加
		}
		codevalue *= 2;										//长度加1
	}
}
void init_Huffman_tables()
{
	compute_Huffman_table(std_dc_luminance_nrcodes, std_dc_luminance_values, YDC_HT);			//获得亮度直流Huffman表
	compute_Huffman_table(std_dc_chrominance_nrcodes, std_dc_chrominance_values, CbDC_HT);		//获得色度直流Huffman表
	compute_Huffman_table(std_ac_luminance_nrcodes, std_ac_luminance_values, YAC_HT);			//获得交流亮度Huffman表	
	compute_Huffman_table(std_ac_chrominance_nrcodes, std_ac_chrominance_values, CbAC_HT);		// Huffman表
}

void set_numbers_category_and_bitcode()
{
	signed long int nr;
	signed long int nrlower, nrupper;
	unsigned char cat, value;

	category_alloc = (unsigned char *)malloc(65535 * sizeof(unsigned char));//对类别分配内存
	category = category_alloc + 32767; 										//allow negative subscripts（允许复数下标）
	bitcode_alloc = (bitstring *)malloc(65535 * sizeof(bitstring));			//对码值分配内存
	bitcode = bitcode_alloc + 32767;		//允许数组负数下标
	nrlower = 1;
	nrupper = 2;
	for (cat = 1; cat <= 15; cat++)
	{
		//Positive numbers（设置正数的类别及编码值）
		for (nr = nrlower; nr < nrupper; nr++)
		{
			category[nr] = cat;											//nr便是实际的数值，cat表示编码的长度
			bitcode[nr].length = cat;									//码字的长度
			bitcode[nr].value = (unsigned short int)nr;					//值对应的编码
		}	
		//Negative numbers（设置负数的类别及编码值）
		for (nr = -(nrupper - 1); nr <= -nrlower; nr++)
		{
			category[nr] = cat;								
			bitcode[nr].length = cat;
			bitcode[nr].value = (unsigned short int)(nrupper - 1 + nr);
		}
		nrlower <<= 1;
		nrupper <<= 1;
	}
}

void DCT(signed char *data, unsigned char *fdtbl, signed short int *outdata)
{
//aanscalefactor[0] = 1,aanscalefactor[k] = cos(k*PI/16) * sqrt(2) for k = 1..7
//量化表系数=1/(原量化系数*行aa&n因子*列aa&n因子*8)
	double aanscalefactor[8] = { 1.0, 1.387039845, 1.306562965, 1.175875602,1.0, 0.785694958, 0.541196100, 0.275899379 };
	float tmp0,  tmp1,  tmp2,  tmp3,  tmp4,  tmp5,  tmp6,  tmp7;
	float tmp10, tmp11, tmp12, tmp13, tmp14, tmp15, tmp16, tmp17, tmp18, tmp19;
	float z1, z2, z3, z4, z5, z11, z13;
	float *dataptr;
	float datafloat[64];
	float temp;
	signed char ctr;
	unsigned char i;
	for (i = 0; i < 64; i++)
		datafloat[i] = data[i];

	/* Pass 1: process rows. */
	dataptr = datafloat;
	for (ctr = 7; ctr >= 0; ctr--)
	{
		tmp0 = dataptr[0] + dataptr[7];
		tmp1 = dataptr[1] + dataptr[6];
		tmp2 = dataptr[2] + dataptr[5];
		tmp3 = dataptr[3] + dataptr[4];
		
		tmp4 = dataptr[3] - dataptr[4];		
		tmp5 = dataptr[2] - dataptr[5];		
		tmp6 = dataptr[1] - dataptr[6];		
		tmp7 = dataptr[0] - dataptr[7];			

		// Phase 2
		tmp10 = tmp0 + tmp3;		
		tmp11 = tmp1 + tmp2;		
		tmp12 = tmp1 - tmp2; 		
		tmp13 = tmp0 - tmp3;		
		
		
		
		tmp14 = (tmp4 + tmp5) * ((float) 0.923879533);		//PI/8
		tmp15 = (tmp4 + tmp5) * ((float) 0.382683432);		
		tmp16 = (tmp6 + tmp7) * ((float) 0.382683432);
		tmp17 = (tmp6 + tmp7) * ((float) 0.923879533);
		tmp18 = tmp7;
		tmp19 = (tmp5 + tmp6) * ((float) 0.707106781);		//sqrt(2)/2

		// Phase 3
		dataptr[0] = (tmp10  + tmp11) / aanscalefactor[0];
		dataptr[4] = (tmp10  - tmp11) / aanscalefactor[4];
		z1 = (tmp12 + tmp13) * ((float) 0.707106781);		//sqrt(2)/2
		z2 = tmp14  - tmp16;
		z4 = tmp15  + tmp17;
		z11 = tmp18 + tmp19;
		z13 = tmp18 - tmp19;

		// Phase 4
		dataptr[2] = (tmp13 + z1)   / aanscalefactor[2];
		dataptr[6] = (tmp13 - z1)   / aanscalefactor[6];
		dataptr[5] = (z13   + z2)   / aanscalefactor[5];
		dataptr[3] = (z13   - z2)   / aanscalefactor[3];
		dataptr[1] = (z11   + z4)   / aanscalefactor[1];
		dataptr[7] = (z11   - z4)   / aanscalefactor[7];

		dataptr += 8;
	}

	/* Pass 2: process columns. */

	dataptr -= 64;
	for (ctr = 7; ctr >= 0; ctr--)
	{
		tmp0 = dataptr[0]  + dataptr[56];
		tmp7 = dataptr[0]  - dataptr[56];
		tmp1 = dataptr[8]  + dataptr[48];
		tmp6 = dataptr[8]  - dataptr[48];
		tmp2 = dataptr[16] + dataptr[40];
		tmp5 = dataptr[16] - dataptr[40];
		tmp3 = dataptr[24] + dataptr[32];
		tmp4 = dataptr[24] - dataptr[32];

		// Phase 2
		tmp10 = tmp0  + tmp3;
		tmp13 = tmp0  - tmp3;
		tmp11 = tmp1  + tmp2;
		tmp12 = tmp1  - tmp2;
		tmp14 = (tmp4 + tmp5) * ((float) 0.923879533);
		tmp15 = (tmp4 + tmp5) * ((float) 0.382683432);
		tmp16 = (tmp6 + tmp7) * ((float) 0.382683432);
		tmp17 = (tmp6 + tmp7) * ((float) 0.923879533);
		tmp18 = tmp7;
		tmp19 = (tmp5 + tmp6) * ((float) 0.707106781);

		// Phase 3
		dataptr[0] = (tmp10 + tmp11)  / aanscalefactor[0];
		dataptr[32] = (tmp10 - tmp11) / aanscalefactor[1];
		z1 = (tmp12 + tmp13) * ((float) 0.707106781);
		z2 = tmp14  - tmp16;
		z4 = tmp15  + tmp17;
		z11 = tmp18 + tmp19;
		z13 = tmp18 - tmp19;

		// Phase 4
		dataptr[16] = (tmp13 + z1) / aanscalefactor[2];
		dataptr[48] = (tmp13 - z1) / aanscalefactor[6];
		dataptr[40] = (z13 + z2)   / aanscalefactor[5];
		dataptr[24] = (z13 - z2)   / aanscalefactor[3];
		dataptr[8]  = (z11 + z4)   / aanscalefactor[1];
		dataptr[56] = (z11 - z4)   / aanscalefactor[7];

		dataptr++;
	}

	// DQT
	for (i = 0; i < 64; i++)
	{
		temp = datafloat[i] / (float)((double)fdtbl[zigzag[i]] * 8.0);					//量化
		//Round to nearest integer.
		//Since C does not specify the direction of rounding for negative
		//quotients, we have to force the dividend positive for portability.
		//The maximum coefficient size is + -16K(for 12 - bit data), so this
		//code should work for either 16 - bit or 32 - bit ints.
		outdata[i] = (signed short int) ((signed short int)(temp + 16384.5) - 16384);
	}
}

void process_DU(signed char *ComponentDU, unsigned char *fdtbl, signed short int *DC, bitstring *HTDC, bitstring *HTAC)
{
	unsigned char i;
	unsigned char startpos;
	unsigned char end0pos;
	unsigned char nrzeroes;
	unsigned char nrmarker;
	signed short int Diff;

	DCT(ComponentDU, fdtbl, DU_DCT);					//DCT和量化后的输出矩阵	

	// 之字型扫描
	for (i = 0; i <= 63; i++) 
		DU[zigzag[i]] = DU_DCT[i];

	Diff = DU[0] - *DC;									//直流分量采用DPCM编码
	*DC = DU[0];

	if (Diff == 0)
		writebits(HTDC[0]); 							//Diff might be 0
	else
	{
		writebits(HTDC[category[Diff]]);				//category[Diff]，求得码字的长度即为huffman编码的值，HTDC[category[Diff]]表示huffman编码的值对应的码字和码字长度
		writebits(bitcode[Diff]);						//再写数据
	}
	for (end0pos = 63; (end0pos > 0) && (DU[end0pos] == 0); end0pos--);	//检测到非0数据跳出循环
	if (end0pos == 0)
	{
		writebits(HTAC[0x00]);							//EOB
		return;
	}

	i = 1;
	while (i <= end0pos) 
	{
		startpos = i;
		for (;(DU[i] == 0) && (i <= end0pos); i++);
		nrzeroes = i - startpos;										//获得前置0的个数
		if (nrzeroes >= 16) 
		{
			
			for (nrmarker = 1; nrmarker <= nrzeroes / 16; nrmarker++) 
				writebits(HTAC[0xF0]);									//说明中间有16个连0
			nrzeroes = nrzeroes % 16;
		}
		writebits(HTAC[nrzeroes * 16 + category[DU[i]]]);				//写码字
		writebits(bitcode[DU[i]]);										//写值
		i++;
	}
	if (end0pos != 63) 													//写入EOB表示一个块结束
		writebits(HTAC[0x00]); 
}

void RGB2YCbCr(unsigned short int xpos, unsigned short int ypos)
{
	unsigned char x, y;
	unsigned char R, G, B;
	for (y = 0; y < 16; y++)
	{
		for (x = 0; x < 16; x++)
		{
			R = buffer[ypos * Ximage * 3 + xpos * 3 + y * Ximage * 3 + x * 3 + 2];
			G = buffer[ypos * Ximage * 3 + xpos * 3 + y * Ximage * 3 + x * 3 + 1];
			B = buffer[ypos * Ximage * 3 + xpos * 3 + y * Ximage * 3 + x * 3 + 0];
			DU_Y[y * 16 + x] = (unsigned char)(0.299 *R + 0.587 *G + 0.114 *B) - 128;				//亮度分量 
			if ((x % 2 == 0) & (y % 2 == 0)) 
			{
				DU_Cb[y / 2 * 8 + x / 2] = (unsigned char)(0.5 *B - 0.33126 *G - 0.16874 *R);		//色差分量Cb
				DU_Cr[y / 2 * 8 + x / 2] = (unsigned char)(0.5 *R - 0.41869 *G - 0.08131 *B);		//色差分量Cr
			}
		}
	}
}

void main_encoder()
{
	signed short int DCY = 0, DCCb = 0, DCCr = 0;						//DC系数用来做差分编码
	unsigned short int xpos, ypos;
	signed char DU[64];
	int i;
	for (ypos = 0; ypos < Yimage; ypos += 16) 
	{
		for (xpos = 0; xpos < Ximage; xpos += 16) 
		{
			RGB2YCbCr(xpos, ypos);
			for (i = 0; i < 64; i++) 									//获得MCU左上角的Y分量矩阵
				DU[i] = DU_Y[(i / 8) * 16 + i % 8 + 0];
			process_DU(DU, DQTinfo.Ytable, &DCY, YDC_HT, YAC_HT);		//第一个参数为数据值，第二个参数为亮度的量化表，第四个参数为亮度直流分量的huffman表，第四个参数为亮度交流分量的huffman表
			for (i = 0; i < 64; i++) 
				DU[i] = DU_Y[(i / 8) * 16 + i % 8 + 8];
			process_DU(DU, DQTinfo.Ytable, &DCY, YDC_HT, YAC_HT);
			for (i = 0; i < 64; i++) 
				DU[i] = DU_Y[(i / 8) * 16 + i % 8 + 128];
			process_DU(DU, DQTinfo.Ytable, &DCY, YDC_HT, YAC_HT);
			for (i = 0; i < 64; i++) 
				DU[i] = DU_Y[(i / 8) * 16 + i % 8 + 136];
			process_DU(DU, DQTinfo.Ytable, &DCY, YDC_HT, YAC_HT);

			process_DU(DU_Cb, DQTinfo.Cbtable, &DCCb, CbDC_HT, CbAC_HT);
			process_DU(DU_Cr, DQTinfo.Cbtable, &DCCr, CbDC_HT, CbAC_HT);
		}
	}
}

void load_bitmap(char *bitmap_name, unsigned short int *Ximage_original, unsigned short int *Yimage_original)
{
	unsigned short int Xdiv8, Ydiv8;
	unsigned char nr_fillingbytes;								//The number of the filling bytes in the BMP file（bmp以4字节对齐，因此有填充）
																//(the dimension in bytes of a BMP line on the disk is divisible by 4)
	unsigned char lastcolor[3];
	unsigned short int column;
	unsigned char TMPBUF[256];
	unsigned short int nrline_up, nrline_dn, nrline;
	unsigned short int dimline;
	unsigned char *tmpline;
	FILE *fp_bitmap = fopen(bitmap_name, "rb");					//打开bmp文件
	fread(TMPBUF, 1, 54, fp_bitmap);							//读入位图的bmp文件头(bmp file header)、位图信息头(bitmap information)
	Ximage = (unsigned short int)TMPBUF[19] * 256 + TMPBUF[18];	//获得图像的宽度
	Yimage = (unsigned short int)TMPBUF[23] * 256 + TMPBUF[22];	//获得图像的高度
	*Ximage_original = Ximage;
	*Yimage_original = Yimage; 									//Keep the old dimensions of the image
	if (Ximage % 16 != 0)										//保证长宽均为MCU宽度的整数倍
		Xdiv8 = (Ximage / 16) * 16 + 16;
	else
		Xdiv8 = Ximage;
	if (Yimage % 16 != 0)
		Ydiv8 = (Yimage / 16) * 16 + 16;
	else
		Ydiv8 = Yimage;

	// The image we encode shall be filled with the last line and the last column
	// from the original bitmap, until Ximage and Yimage are divisible by 8
	// Load BMP image from disk and complete X
	buffer = (unsigned char *)(malloc(3 * Xdiv8 * Ydiv8));			//24位位图一个像素占24个字节
	if (Ximage % 4 != 0) 											//读取文件数一行字节数必须是4的倍数，如果不够四字节要补齐四字节
		nr_fillingbytes = 4 - (Ximage % 4);							//填充的字节数
	else
		nr_fillingbytes = 0;
	for (nrline = 0; nrline < Yimage; nrline++)						//读出每一行数据并在其后填充以满足与16求余为0
	{
		fread(buffer + nrline * Xdiv8 * 3, 1, Ximage * 3, fp_bitmap);			//读取图像一行的数据
		fread(TMPBUF, 1, nr_fillingbytes, fp_bitmap);							//读取填充的字节数
		memcpy(&lastcolor, buffer + nrline * Xdiv8 * 3 + Ximage * 3 - 3, 3);	//读出上一像素的RGB值用于扩展
		for (column = Ximage; column < Xdiv8; column++)
		{
			memcpy(buffer + nrline * Xdiv8 * 3 + column * 3, &lastcolor, 3);	//扩展列
		}
	}
	Ximage = Xdiv8;																//扩展后的图像宽度
	dimline = Ximage * 3;														//扩展后每一行所包含的字节数
	tmpline = (unsigned char *)malloc(dimline);

	//Reorder in memory the inversed bitmap（位图信息头的图像高度是正数，所以位图数据在文件中的排列顺序是从左下角到右上角的，以行为主序排列的）
	for (nrline_up = Yimage - 1, nrline_dn = 0; nrline_up > nrline_dn; nrline_up--, nrline_dn++)
	{
		memcpy(tmpline, buffer + nrline_up * Ximage * 3, dimline);
		memcpy(buffer + nrline_up * Ximage * 3, buffer + nrline_dn * Ximage * 3, dimline);
		memcpy(buffer + nrline_dn * Ximage * 3, tmpline, dimline);
	}

	//Y completion:
	memcpy(tmpline, buffer + (Yimage - 1) * Ximage * 3, dimline);				//行扩展为Ydiv8行
	for (nrline = Yimage; nrline < Ydiv8; nrline++)
	{
		memcpy(buffer + nrline * Ximage * 3, tmpline, dimline);
	}

	Yimage = Ydiv8;				//扩展后的图像长宽
	free(tmpline);				//释放内存空间
	fclose(fp_bitmap);			//关闭bmp图像文件
}

void init_all()
{
	set_DQTinfo();						//设置量化表信息
	set_DHTinfo();						//设置哈夫曼表信息
	init_Huffman_tables();				//初始化哈夫曼表
	set_numbers_category_and_bitcode();	//设置数值的类别（编码位数及编码值）
}

/*
int main(int argc,char *argv[],char **env)
第一个参数，int型的argc，为整型，用来统计程序运行时发送给main函数的命令行参数的个数
第二个参数，char*型的argv[]，为字符串数组，用来存放指向的字符串参数的指针数组，每一个元素指向一个参数。
argv[0] 指向程序运行的全路径名
argv[1] 指向在DOS命令行中执行程序名后的第一个字符串
argv[2] 指向执行程序名后的第二个字符串
argv[3] 指向执行程序名后的第三个字符串
argv[argc] 为NULL
第三个参数，char**型的env，为字符串数组。env[]的每一个元素都包含ENVVAR=value形式的字符串，其中ENVVAR为环境变量，value为其对应的值。平时使用到的比较少。
*/


int main(int argc, char *argv[])
{
	char BMP_filename[64];									//保存文件名（bmp文件名）
	char JPG_filename[64];									//保存Jpeg压缩后的文件名
	unsigned short int Ximage_original, Yimage_original; 	//the original image dimensions,
															//before we made them divisible by 8
	unsigned char len_filename;
	bitstring fillbits; 									//filling bitstring for the bit alignment of the EOI marker
	if (argc > 1) 											//大于1说明有其他参数输入
	{
		strcpy(BMP_filename, argv[1]);
		if (argc > 2)
			strcpy(JPG_filename, argv[2]);
		else
		{
			strcpy(JPG_filename, BMP_filename);
			len_filename = strlen(BMP_filename);
			strcpy(JPG_filename + (len_filename - 3), "jpg");
		}
	}
	load_bitmap(BMP_filename, &Ximage_original, &Yimage_original);		//将图像获取到buffer所指向的内存空间中
	fp_jpeg_stream = fopen(JPG_filename, "wb");							//新建一个二进制文件，只允许写，用于存放压缩后的Jpeg数据
	init_all();												//初始化一些标记头文件，以及开辟一个存储空间用于存放正负数的编码（正数使用源码，负数使用反码）
	SOF0info.width = Ximage_original;
	SOF0info.height = Yimage_original;

	writeword(0xFFD8);	//写入SOI图像开始标志

	write_APP0info();	//写入JFIF应用数据块
	write_DQTinfo();	//写入量化表
	write_SOF0info();	//帧图像分开始
	write_DHTinfo();	//写入哈夫曼表
	write_SOSinfo();	//扫描开始

	bytenew = 0;														//The byte that will be written in the JPG file
	bytepos = 7;														//bit position in the byte we write (bytenew)
																		//should be <= 7 and >=0	
	main_encoder();		//编码子函数

	//Do the bit alignment of the EOI marker
	if (bytepos >= 0)
	{
		fillbits.length = bytepos + 1;
		fillbits.value = (1 << (bytepos + 1)) - 1;
		writebits(fillbits);
	}
	writeword(0xFFD9); 			//EOI图像结束标志

	free(buffer);				//释放内存
	free(category_alloc);
	free(bitcode_alloc);
	fclose(fp_jpeg_stream);
}
Da Liu
关注
JPEG解码_C语言实现
刘东的博客
06-01 1万+
一、实验原理 1.JPEG编解码原理 JPEG（Joint Photographic Experts Group）采用有损压缩方式去除冗余的图像和彩色数据，在获得极高压缩率的同时能展现十分丰富生动的图像。 编码原理图如下: ① Level Offset（零偏置） 为了使使像素的绝对值出现3位10进制的概率大大减少，对于灰度级是2的n次方的像素，通过减去2的n-1次方，将无符号的整数值变成
JPEG编解码的c语言实现
02-25
JPEG编解码的实现，包括哈夫曼编码，解码。还有JPG转BMP格式程序.
C语言编写的JPEG编码输出.JPG文件（已在Cortex-A9内核上验证通过）
05-26
C语言编写的JPEG编码输出.JPG文件（已在Cortex-A9内核上验证通过）
cjpeg配置
最新发布
这是一个c++热爱者的博客哟
07-22 106
【代码】cjpeg配置。
数字图像处理(1)：JPEG分析 & 实践（C语言）
ZhaoDongyu_AK47的博客
06-26 2848
文章目录1. JPEG格式2. JPEG原理&压缩流程2.1 步骤一：图像分割2.2 步骤二：颜色空间转换RGB->YCbCr2.3 步骤三：离散余弦变换2.4 步骤四：数据量化2.5 步骤五：哈弗曼编码3. 下载实现代码并调通运行3.1 测试平台3.2 测试图片3.3 压缩后输出图片4. 压缩率附：代码 1. JPEG格式 JPEG是Joint Photographic Experts Group(联合图像专家组)的缩写，文件后辍名为"．jpg"或"．jpeg"，是最常用的图像文件格式，由一
JPEG.zip_C jpeg 源代码
09-24
在"JPEG.zip"压缩包中，包含了第2章的源代码，这部分代码可能详细阐述了JPEG编码流程中的关键步骤，包括DCT计算、量化、哈夫曼编码和文件写入等。通过阅读和理解这些源代码，开发者可以更深入地了解JPEG压缩的实现...
基于C语言的Motion JPEG2000文件的构造
09-25
2. **识别和提取JPEG2000编码流**：使用C语言编写代码，从数据文件中查找并提取JPEG2000编码流。这一步骤至关重要，因为编码流的正确性和完整性直接影响到最终文件的播放效果。 3. **合并编码流与文件头和文件尾**...
JPEG格式编码示例代码
08-03
总之，这个“JPEG格式编码示例代码”提供了一个基础的JPEG编码框架，可以帮助开发者深入理解JPEG编码流程，并动手实现自己的JPEG编码器。通过研究和实践，不仅可以提升图像处理技术，还能增强对数据压缩原理的认识。
JPEG编解码的C语言实现
07-01
JPEG编解码的实现，包括哈夫曼编码，解码。还有JPG转BMP格式程序.
jpeg压缩纯C语言实现
05-07
利用visual studio 2013平台，实现了对灰度图像的JPEG压缩，纯C语言实现，可直接移植到DSP，单片机等嵌入式平台当中
JPEG图像编码
03-12
经典JPEG编码 可实现灰色图像、彩色图像编码 调试可通过
JPEG压缩编码代码
04-19
包含DCT变换频谱显示、8*8子块的DCT变换图像压缩的代码等。
二、基于C语言的JPEG解码代码详解
m0_37633745的博客
08-13 3510
// JPEG解码 // 输入JPEG后输出Bitmap #include <stdio.h> #include <stdlib.h> #pragma warning(disable : 4996) unsigned int BuffIndex; // JPEG数据的位置 unsigned int BuffSize; // JPEG数据的大小 unsigned int BuffX; // 图像的横向尺寸 unsigned int BuffY; //
JPEG
rosenberg的专栏
05-14 4784
引用:http://www.cnblogs.com/laughterwym/archive/2006/06/23/433793.html一、JPEG算法概要　　JPEG（Joint Photographic Experts Group）是一个由ISO和IEC两个组织机构联合组成的一个专家组，负责制定静态的数字图像数据压缩编码标准，这个专家组开发的算法称为JPEG算法，并且成为国际上通用的标准，
C语言实现的BMP和JPEG图片的解码
YANGSENG1987的专栏
08-05 8646
/* 这是read_picture.c 文件的代码 ，主要是对输入的图片的文件（BMP和JPEG类型）进行解码，转换成在LCD中显示的数据； 使用该程序时，要先把jpeglib库文件和头文件放到对应的标谁库和标准头文件中；否则编译出错的 */ #include #include #include #include #include #include #include
评论 5
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Da Liu

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值