【正文】 inance_qt,scalefactor,)。writebyte()。for(i=0。for(i=0。i=161。i16。i++) writebyte([i])。i++) writebyte([i])。}void set_DHTinfo(){BYTE i。for(i=0。i=11。i16。i++) [i]=std_ac_luminance_values[i]。i++) [i]=std_dc_chrominance_nrcodes[i+1]。=0x11。for(i=0。writeword()。writebyte()。writebyte()。writeword(0xFFFE)。ilength。SBYTE posval。 mask[posval]) bytenew|=mask[bytepos]。writebyte(0)。}}}void pute_Huffman_table(BYTE *nrcodes,BYTE *std_table,bitstring *HT){BYTE k, j。pos_in_table=0。j=nrcodes[k]。codevalue++。pute_Huffman_table(std_ac_luminance_nrcodes,std_ac_luminance_values,YAC_HT)。}void set_numbers_category_and_bitcode(){SDWORD nr。if(category_alloc=NULL)exitmessage(Not enough memory.)。bitcode=bitcode_alloc+32767。cat=15。bitcode[nr].length=cat。nr++){category[nr]=cat。nrupper=1。R++) {YRtab[R]=(SDWORD)(65536*+)*R。G=255。}for(B=0。CrBtab[B]=(SDWORD)(65536*+)*B。N IDCT method, divisors are equal to quantization coefficients scaled by scalefactor[row]*scalefactor[col], where scalefactor[0] = 1scalefactor[k] = cos(k*PI/16) * sqrt(2) for k=l. .7 We apply a further scale factor of 8.What39。for (row=0。col++){fdtbl_Y[i]=(float)(((double)[zigzag[i]]*aanscalefactor[row]*aanscalefactor[col]*))。float tmp10,tmp11,tmp12,tmp13。float temp。i64。ctr=0。tmp6=dataptr[1]dataptr[6]。tmp4=dataptr[3]dataptr[4]。tmp12=tmp1tmp2。/* c4 */ dataptr[2]=tmp13+z1。tmp12=tmp6+tmp7。/*c2+c6 */ z3=tmp11*((float))。/* phase 6 */ dataptr[3]=z13z2。/* advance pointer to next row */}/* Pass 2: process columns. */dataptr=datafloat。tmp7=dataptr[0]dataptr[56]。tmp5=dataptr[16]dataptr[40]。 /* phase 2 */ tmp13=tmp0tmp3。 /* phase 3 */ dataptr[32]=tmp10tmp11。/* Odd part */tmp10=tmp4+tmp5。/* c6 */ z2=((float))*tmp10+z5。/* phase 5 */ z13=tmp7z3。dataptr[56]=z11z4。i++){/* Apply the quantization and scaling factor */temp=datafloat[i]*fdtbl[i]。BYTE i。BYTE nrmarker。i=63。//Encode DCif(Diff=0) writebits(HTDC[0])。(endOpos0)amp。//endOpos=first element in reverse order !=0 if(endOpos==0){writebits(EOB)。for(。i++)。nrmarker++)writebits(M16zeroes)。}if(endOpos!=63)writebits(EOB)。BYTE R,G,B。y++){for(x=0。B=RGB_buffer[location].B。location++。//DC coefficients used for differential encoding WORD xpos,ypos。xposXimage。process_DU(CbDU,fdtbl_Cb,amp。}}void load_bitmap(char *bitmap_name, WORD *Ximage_original,WORD *Yimage_original){WORD Xdiv8,Ydiv8。BYTE TMPBUF[256]。FILE *fp_bitmap=fopen(bitmap_name,rb)。B39。Ximage=(WORD)TMPBUF[19]*256^TMPBUF[18]。//Keep the old dimensions// of the imageif(Ximage%8!=0)Xdiv8=(Ximage/8)*8*8。// 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 plete X RGB_buffer=(colorRGB*) (malloc(3*Xdiv8*Ydiv8)) 。for(nrline=0。memcpy(amp。column++){memcpy(RGB_buffernrline*Xdiv8+column,amp。tmpline=(colorRGB *)malloc(dimline)。nrline_up,nrline_dn++){memcpy(tmpline,RGB_buffer+nrline_up*Ximage,dimline)。//Reorder in memory the inversed bitmapfor(nrline_up=Yimage1,nrline_dn=0。}}Ximage=Xdiv8。for(column=Ximage。nrline++){fread(RGB_buffer+nrline*Xdiv8,1,Ximage*3,fp_bitmap)。if(Ximage%4!=0) nr_fillingbytes=4(Ximage%4)。if(Yimage%8!=0) Ydiv8=(Yimage/8)*8+8。*Ximage_original=Ximage。M39。if(fread(TMPBUF,1,54,fp_bitmap)!=54)exitmessage(Need a truecolor BMP to encode.)。WORD dimline。//The number of the filling bytes in the BMP file// (the dimension in bytes of a BMP line on the disk is divisible by 4)colorRGB lastcolor。process_DU(CrDU,fdtbl_Cb,amp。process_DU(YDU,fdtbl_Y,amp。yposYimage。}location+=Ximage8。CbDU[pos]=Cb(R,G,B)。x++){R=RGB_buffer[location].R。for(y=0。BYTE pos=0。}writebits(HTAC[nrzeroes*16+category[DU[i]]])。if(nrzeroes=16){for(nrmarker=1。amp。}i=1。(DU[endOpos]==0)。 writebits(bitcode[Diff])。Diff=DU[0]*DC。fdct_and_quantization(ComponentDU,fdtbl,DU_DCT)。BYTE endOpos。}}void process_DU(SBYTE *ComponentDU,float *fdtbl,SWORD *DC,bitstring *HTDC,bitstring *HTAC){bitstring EOB=HTAC[0x00]。 /* advance pointer to next column */}// Quantize/descale the coefficients, and store into output array for(i=0。/* phase 6 */ dataptr[24]=z13z2。/* c2+c6 */ z3=tmp11*((float))。tmp12=tmp6+tmp7。/*c4*/dataptr[16]=tmp13+z1。tmp12=tmp1tmp2。tmp4=dataptr[24]dataptr[32]。tmp6=dataptr[8]dataptr[48]。ctr=0。dataptr[7]=z11z4。/* phase 5 */z13=tmp7z3。/* c6 */ z2=((float))*tmp10+z5。/* Odd part */tmp10=tmp4+tmp5。/* phase 3 */ dataptr[4]=tmp10tmp11。/*phase 2 */ tmp13=tmp0tmp3。tmp5=dataptr[2]dataptr[5]。tmp7=dataptr[0]dataptr[7]。/* Pass 1: process rows. */ dataptr=datafloat。BYTE i。float *dataptr。i++。row++){ for(col=0。BYTE row, col。 s C source:// Forward DCT routine idea taken from Independent JPEG Group39。B++) {YBtab[B]=(SDWORD)(65536*+)*B。CbGtab[G]=(SDWORD)(65536*+)*G。CrRtab[R]=(SDWORD)(32768)*R 。for(R=0。bitcode[nr].value=(WORD)(nrupper1+nr)。}//Negative numbersfor(nr=(nrupper1)。nrnrupper。nrupper=2。//allow negative subscripts bitcode_alloc=(bitstring *)malloc(65535*sizeof(bitstring))。BYTE cat,value。}void exitmessage(char *error_message){printf(%s\n,error_message)。}}void init_Huffman_tables(){pute_Huffman_table(std_dc_luminance_nrcodes,std_dc_luminance_values,YDC_HT)。HT[std_table[pos_in_table]].length=k。k=16。WORD codevalue。}bytepos=7。bytepos。posval=。}void writebits(bitstring bs)// A portable version。writeword(length+2)。writebyte()。writebyte()。writebyte()。i++) [i]=std_ac_chrominance_values[i]。i16。i=11。for(i=0。for (i=0。=0x10。i++) [i]=std_dc_luminance_nrcodes[i+1]。=0x01A2。i=161。for(i=0。for(i=0。writebyte()。i++) writebyte([i])。i++) writebyte([i])。i16。writeword()。=1。ll get, and the smaller // pression we39。if(temp255L) temp=255L。for(i=0。i64。i64。writeword()。w