⭐ 欢迎来到虫虫下载站! | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站

📄 compress.c

📁 图像处理入门书
💻 C
📖 第 1 页 / 共 3 页
字号:
	if((funcret=InitTag())!=FUNC_OK){

		GlobalUnlock(hJpegBuf);

		GlobalFree(hJpegBuf);

		showerror(funcret);

		return FALSE;

	}

	//create new bitmapfileheader and bitmapinfoheader

	memset((char *)&bf,0,sizeof(BITMAPFILEHEADER));	

	memset((char *)&bi,0,sizeof(BITMAPINFOHEADER));



	bi.biSize=(DWORD)sizeof(BITMAPINFOHEADER);

	bi.biWidth=(LONG)(ImgWidth);

	bi.biHeight=(LONG)(ImgHeight);

	bi.biPlanes=1;

	bi.biBitCount=24;

	bi.biClrUsed=0;

	bi.biClrImportant=0;

	bi.biCompression=BI_RGB;

	NumColors=0;

	LineBytes=(DWORD)WIDTHBYTES(bi.biWidth*bi.biBitCount);

	ImgSize=(DWORD)LineBytes*bi.biHeight;



	bf.bfType=0x4d42;

	bf.bfSize=sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+NumColors*sizeof(RGBQUAD)+ImgSize;

	bf.bfOffBits=(DWORD)(NumColors*sizeof(RGBQUAD)+sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER));

	BufSize=bf.bfSize-sizeof(BITMAPFILEHEADER);



	if((hImgData=GlobalAlloc(GHND,BufSize))==NULL){

		GlobalUnlock(hJpegBuf);

		GlobalFree(hJpegBuf);

		showerror(FUNC_MEMORY_ERROR);

		return FALSE;

	}

	lpImgData=(LPBITMAPINFOHEADER)GlobalLock(hImgData); 

	memcpy(lpImgData,(char *)&bi,sizeof(BITMAPINFOHEADER));

	lpPtr=(char *)lpImgData+sizeof(BITMAPINFOHEADER);



	if((SampRate_Y_H==0)||(SampRate_Y_V==0)){

		GlobalUnlock(hJpegBuf);

		GlobalFree(hJpegBuf);

		GlobalUnlock(hImgData);

		GlobalFree(hImgData);

		hImgData=NULL;

		showerror(FUNC_FORMAT_ERROR);

		return FALSE ;

	}



	funcret=Decode();

	if(funcret==FUNC_OK){

		hDc=GetDC(hWnd);

		hBitmap=CreateDIBitmap(hDc,	(LPBITMAPINFOHEADER)lpImgData, (LONG)CBM_INIT,

						(LPSTR)lpImgData+sizeof(BITMAPINFOHEADER) +NumColors*sizeof(RGBQUAD),

   						(LPBITMAPINFO)lpImgData, DIB_RGB_COLORS);



		hfbmp=_lcreat("c:\\jpeg2bmp.bmp",0);

		_lwrite(hfbmp,(LPSTR)&bf,sizeof(BITMAPFILEHEADER)); 

		_lwrite(hfbmp,(LPSTR)lpImgData,BufSize);

		_lclose(hfbmp);

		ReleaseDC(hWnd,hDc);

		GlobalUnlock(hJpegBuf);

		GlobalFree(hJpegBuf);

		GlobalUnlock(hImgData);

		return TRUE;

	}

	else{

		GlobalUnlock(hJpegBuf);

		GlobalFree(hJpegBuf);

		GlobalUnlock(hImgData);

		GlobalFree(hImgData);

		hImgData=NULL;

		showerror(funcret);

		return FALSE;

	}

}

/////////////////////////////////////////////////

void showerror(int funcret)

{

	switch(funcret){

	case FUNC_MEMORY_ERROR:

    	MessageBox(NULL,"Error alloc memory!","ErrorMessage",MB_OK|MB_ICONEXCLAMATION);

		break;

	case FUNC_FILE_ERROR:

    	MessageBox(NULL,"File not found!","ErrorMessage",MB_OK|MB_ICONEXCLAMATION);

		break;

	case FUNC_FORMAT_ERROR:

		MessageBox(NULL,"File format error!","Error Message",MB_OK|MB_ICONEXCLAMATION);

		break;

	}

}

////////////////////////////////////////////////////////////////////////////////

int InitTag()

{

	BOOL finish=FALSE;

	BYTE id;

	short  llength;

	short  i,j,k;

	short  huftab1,huftab2;

	short  huftabindex;

	BYTE hf_table_index;

	BYTE qt_table_index;

	BYTE comnum;



	unsigned char  *lptemp;

	short  ccount;



	lp=lpJpegBuf+2;



	while (!finish){

		id=*(lp+1);

		lp+=2;

		switch (id){

		case M_APP0:

			llength=MAKEWORD(*(lp+1),*lp);

			lp+=llength;

			break;

		case M_DQT:

			llength=MAKEWORD(*(lp+1),*lp);

			qt_table_index=(*(lp+2))&0x0f;

			lptemp=lp+3;

			if(llength<80){

				for(i=0;i<64;i++)

					qt_table[qt_table_index][i]=(short)*(lptemp++);

			}

			else{

				for(i=0;i<64;i++)

					qt_table[qt_table_index][i]=(short)*(lptemp++);

                qt_table_index=(*(lptemp++))&0x0f;

  				for(i=0;i<64;i++)

					qt_table[qt_table_index][i]=(short)*(lptemp++);

  			}

  			lp+=llength;		

			break;

		case M_SOF0:

	 		llength=MAKEWORD(*(lp+1),*lp);

	 		ImgHeight=MAKEWORD(*(lp+4),*(lp+3));

	 		ImgWidth=MAKEWORD(*(lp+6),*(lp+5));

            comp_num=*(lp+7);

		    if((comp_num!=1)&&(comp_num!=3))

	  			return FUNC_FORMAT_ERROR;

			if(comp_num==3){

				comp_index[0]=*(lp+8);

	  			SampRate_Y_H=(*(lp+9))>>4;

	  			SampRate_Y_V=(*(lp+9))&0x0f;

	  			YQtTable=(short *)qt_table[*(lp+10)];



				comp_index[1]=*(lp+11);

				SampRate_U_H=(*(lp+12))>>4;

	  			SampRate_U_V=(*(lp+12))&0x0f;

	  			UQtTable=(short *)qt_table[*(lp+13)];



	  			comp_index[2]=*(lp+14);

	  			SampRate_V_H=(*(lp+15))>>4;

	  			SampRate_V_V=(*(lp+15))&0x0f;

				VQtTable=(short *)qt_table[*(lp+16)];

	  		}

			else{

	  			comp_index[0]=*(lp+8);

				SampRate_Y_H=(*(lp+9))>>4;

	  			SampRate_Y_V=(*(lp+9))&0x0f;

	  			YQtTable=(short *)qt_table[*(lp+10)];



				comp_index[1]=*(lp+8);

	  			SampRate_U_H=1;

	  			SampRate_U_V=1;

	  			UQtTable=(short *)qt_table[*(lp+10)];



				comp_index[2]=*(lp+8);

				SampRate_V_H=1;

	  			SampRate_V_V=1;

	  			VQtTable=(short *)qt_table[*(lp+10)];

			}

  			lp+=llength;						    

			break;

		case M_DHT:             

			llength=MAKEWORD(*(lp+1),*lp);

			if (llength<0xd0){

				huftab1=(short)(*(lp+2))>>4;     //huftab1=0,1

		 		huftab2=(short)(*(lp+2))&0x0f;   //huftab2=0,1

				huftabindex=huftab1*2+huftab2;

		 		lptemp=lp+3;

				for (i=0; i<16; i++)

					code_len_table[huftabindex][i]=(short)(*(lptemp++));

				j=0;

				for (i=0; i<16; i++)

					if(code_len_table[huftabindex][i]!=0){

						k=0;

						while(k<code_len_table[huftabindex][i]){

							code_value_table[huftabindex][k+j]=(short)(*(lptemp++));

							k++;

						}

						j+=k;	

					}

				i=0;

				while (code_len_table[huftabindex][i]==0)

		 			i++;

				for (j=0;j<i;j++){

					huf_min_value[huftabindex][j]=0;

					huf_max_value[huftabindex][j]=0;

				}

				huf_min_value[huftabindex][i]=0;

				huf_max_value[huftabindex][i]=code_len_table[huftabindex][i]-1;

				for (j=i+1;j<16;j++){

					huf_min_value[huftabindex][j]=(huf_max_value[huftabindex][j-1]+1)<<1;

					huf_max_value[huftabindex][j]=huf_min_value[huftabindex][j]+code_len_table[huftabindex][j]-1;

				}

				code_pos_table[huftabindex][0]=0;

				for (j=1;j<16;j++)

		  			code_pos_table[huftabindex][j]=code_len_table[huftabindex][j-1]+code_pos_table[huftabindex][j-1];

		  		lp+=llength;

			}  //if

			else{

	 			hf_table_index=*(lp+2);

				lp+=2;

				while (hf_table_index!=0xff){

					huftab1=(short)hf_table_index>>4;     //huftab1=0,1

			 		huftab2=(short)hf_table_index&0x0f;   //huftab2=0,1

					huftabindex=huftab1*2+huftab2;

					lptemp=lp+1;

					ccount=0;

					for (i=0; i<16; i++){

						code_len_table[huftabindex][i]=(short)(*(lptemp++));

						ccount+=code_len_table[huftabindex][i];

					}

					ccount+=17;	

					j=0;

					for (i=0; i<16; i++)

						if(code_len_table[huftabindex][i]!=0){

							k=0;

							while(k<code_len_table[huftabindex][i])

							{

								code_value_table[huftabindex][k+j]=(short)(*(lptemp++));

								k++;

							}

							j+=k;

						}

					i=0;

					while (code_len_table[huftabindex][i]==0)

						i++;

					for (j=0;j<i;j++){

						huf_min_value[huftabindex][j]=0;

						huf_max_value[huftabindex][j]=0;

					}

					huf_min_value[huftabindex][i]=0;

					huf_max_value[huftabindex][i]=code_len_table[huftabindex][i]-1;

					for (j=i+1;j<16;j++){

						huf_min_value[huftabindex][j]=(huf_max_value[huftabindex][j-1]+1)<<1;

						huf_max_value[huftabindex][j]=huf_min_value[huftabindex][j]+code_len_table[huftabindex][j]-1;

					}

					code_pos_table[huftabindex][0]=0;

					for (j=1;j<16;j++)

						code_pos_table[huftabindex][j]=code_len_table[huftabindex][j-1]+code_pos_table[huftabindex][j-1];

					lp+=ccount;

					hf_table_index=*lp;

				}  //while

			}  //else

			break;

		case M_DRI:

			llength=MAKEWORD(*(lp+1),*lp);

			restart=MAKEWORD(*(lp+3),*(lp+2));

			lp+=llength;

			break;

		case M_SOS:

			llength=MAKEWORD(*(lp+1),*lp);

			comnum=*(lp+2);

			if(comnum!=comp_num)

				return FUNC_FORMAT_ERROR;

			lptemp=lp+3;

			for (i=0;i<comp_num;i++){

				if(*lptemp==comp_index[0]){

					YDcIndex=(*(lptemp+1))>>4;   //Y

					YAcIndex=((*(lptemp+1))&0x0f)+2;

				}

				else{

					UVDcIndex=(*(lptemp+1))>>4;   //U,V

					UVAcIndex=((*(lptemp+1))&0x0f)+2;

				}

				lptemp+=2;

			}

			lp+=llength;

			finish=TRUE;

			break;

		case M_EOI:    

			return FUNC_FORMAT_ERROR;

			break;

		default:

 			if ((id&0xf0)!=0xd0){

				llength=MAKEWORD(*(lp+1),*lp);

	 			lp+=llength;

			}

			else lp+=2;

			break;

  		}  //switch

	} //while

	return FUNC_OK;

}

/////////////////////////////////////////////////////////////////

void InitTable()

{

	short i,j;

	sizei=sizej=0;

	ImgWidth=ImgHeight=0;

	rrun=vvalue=0;

	BitPos=0;

	CurByte=0;

	IntervalFlag=FALSE;

	restart=0;

	for(i=0;i<3;i++)

		for(j=0;j<64;j++)

			qt_table[i][j]=0;

	comp_num=0;

	HufTabIndex=0;

	for(i=0;i<3;i++)

		comp_index[i]=0;

	for(i=0;i<4;i++)

		for(j=0;j<16;j++){

			code_len_table[i][j]=0;

			code_pos_table[i][j]=0;

			huf_max_value[i][j]=0;

			huf_min_value[i][j]=0;

		}

	for(i=0;i<4;i++)

		for(j=0;j<256;j++)

			code_value_table[i][j]=0;

	

	for(i=0;i<10*64;i++){

		MCUBuffer[i]=0;

		QtZzMCUBuffer[i]=0;

	}

	for(i=0;i<64;i++){

		Y[i]=0;

		U[i]=0;

		V[i]=0;

		BlockBuffer[i]=0;

	}

	ycoef=ucoef=vcoef=0;

}

/////////////////////////////////////////////////////////////////////////

int Decode()

{

	int funcret;



	Y_in_MCU=SampRate_Y_H*SampRate_Y_V;

	U_in_MCU=SampRate_U_H*SampRate_U_V;

	V_in_MCU=SampRate_V_H*SampRate_V_V;

	H_YtoU=SampRate_Y_H/SampRate_U_H;

	V_YtoU=SampRate_Y_V/SampRate_U_V;

	H_YtoV=SampRate_Y_H/SampRate_V_H;

	V_YtoV=SampRate_Y_V/SampRate_V_V;

	Initialize_Fast_IDCT();

	while((funcret=DecodeMCUBlock())==FUNC_OK){

		interval++;

		if((restart)&&(interval % restart==0))

			 IntervalFlag=TRUE;

		else

			IntervalFlag=FALSE;

		IQtIZzMCUComponent(0);

		IQtIZzMCUComponent(1);

		IQtIZzMCUComponent(2);

		GetYUV(0);

		GetYUV(1);

		GetYUV(2);

		StoreBuffer();

		sizej+=SampRate_Y_H*8;

		if(sizej>=ImgWidth){

			sizej=0;

			sizei+=SampRate_Y_V*8;

		}

		if ((sizej==0)&&(sizei>=ImgHeight))

			break;

		}

	return funcret;

}

/////////////////////////////////////////////////////////////////////////////////////////

void  GetYUV(short flag)

{

	short	H,VV;

	short	i,j,k,h;

	int		*buf;

	int		*pQtZzMCU;



	switch(flag){

	case 0:

		H=SampRate_Y_H;

		VV=SampRate_Y_V;

		buf=Y;

		pQtZzMCU=QtZzMCUBuffer;

		break;

	case 1:

		H=SampRate_U_H;

		VV=SampRate_U_V;

		buf=U;

		pQtZzMCU=QtZzMCUBuffer+Y_in_MCU*64;

		break;

	case 2:

		H=SampRate_V_H;

		VV=SampRate_V_V;

		buf=V;

		pQtZzMCU=QtZzMCUBuffer+(Y_in_MCU+U_in_MCU)*64;

		break;

	}

	for (i=0;i<VV;i++)

		for(j=0;j<H;j++)

			for(k=0;k<8;k++)

				for(h=0;h<8;h++)

					buf[(i*8+k)*SampRate_Y_H*8+j*8+h]=*pQtZzMCU++;

}

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -