📄 jpeg-ls-decode.cpp
字号:
RowX++;
if(RowX>y)
{
RowX--;
flag1=1;
}
i++;
if(flag1==1)
break;
}
if((i==(1<<J[RUNindex]))&&RUNindex<31)
RUNindex++;
if(size<0&&flag1==0)
{
fread(&buffer,1,1,fp);
decode=buffer;
size=7;
}/*一个单元处理完且游程仍未结束，处理下一个单元*/
if(flag1==1)
break;
}
/***在已处理的位上置零***/
if(RowX==y&&flag1==1){/*情况1：在行末*/
if(size==-1)/*正好处理完一个单元且游程未结束*/
decode=0;/*处理器清零，完成游程解码*/
else/*对单元未处理完就已到行末*/
decode=(decode<<(32-size-1))>>(32-size-1);/*当被行末中断时，其剩下的rk位1前缀是1，所以也被前面的处理过size＋1是未处理的，所以32－size－1表示在当前处理范围内已处理的。此操作表示将已处理的位置零*/
size=size+1;/*待处理的位数是size＋1*/
flag1=0;
return;}/*游程解码结束*/
if(size==0)/*情况2：非行末，正好一个单元的最后一个是不同象素*/
decode2=0;
else
decode2=(decode<<(32-size))>>(32-size);/*剩下size个元素待处理*/
/***对长度小于rm的分段码字的处理***/
if(J[RUNindex]<=size)/*当有剩余比特出现时，且在当前处理单元内*/
{
output=(decode2>>(size-J[RUNindex]));/*将待处理的数据移至最右侧，即得到待处理的数据*/
if(J[RUNindex]==size)
decode=0;
else
decode=((decode<<(32-size+J[RUNindex]))>>(32-size+J[RUNindex]));/*将处理过的位置零*/
size=size-J[RUNindex];
}
else/*当前处理单元剩余的比特小于rk时*/
{
fread(&buffer,1,1,fp);
if(J[RUNindex]-size<=8)/*rk在下一个处理单元之内*/
decode1=decode2<<(J[RUNindex]-size);/*加入下一个单元中剩下部分的空间*/
while(J[RUNindex]-size>8)
{
decode1=decode1+(decode2<<(J[RUNindex]-size))+(buffer<<(J[RUNindex]-size-8));/*读入后面几个单元中rk部分.注意这里的加并不表示长度相加，而是在同一原点上叠加*/
J[RUNindex]=J[RUNindex]-8;
fread(&buffer,1,1,fp);
decode2=0;
}
output=decode1+(buffer>>(8+size-J[RUNindex]));/*右边第二项是当rk－size<8时，得到当前单元的剩余的rk，以接到n个满单元之后*/
if(J[RUNindex]-size==8)
decode=0;
else
decode=(buffer<<(32-8+J[RUNindex]-size))>>(32-8+J[RUNindex]-size);/*对前面已经处理过的（n个满单元及当前单元剩余rk）位置零*/
size=8-J[RUNindex]+size;
}
if(RUNindex>0)
RUNindex--;
for(j=0;j<output;j++){
*(f+LinX*(y+2)+RowX)=Ra;
RowX++;}
if(RowX>y){
RowX=RowX-y;
LinX++;
}
if(decode!=0){
while(((1<<(size-1))&decode)==0){
size--;
cnt++;
}
size--;
if(size==0)
decode=0;
else
decode=(decode<<(32-size))>>(32-size);
}
else{
while(decode==0)
{
fread(&buffer,1,1,fp);
decode=buffer;
size=size+8;
cnt=size;
}
int q=7;
while(((1<<q)&decode)==0)
q--;
cnt=cnt-8+7-q;
if(q==0)
decode=0;
else
decode=(decode<<(32-q))>>(32-q);
size=q;
}
Rb=*(f+(LinX-1)*(y+2)+RowX);
if(abs(Ra-Rb)<=NEAR)
RItype=1;
else
RItype=0;
if(RItype==1)
Px=Ra;
else
Px=Rb;
if((RItype==0)&&(Ra>Rb))
{
SIGN=-1;
}
else
SIGN=1;
if(RItype==0)
TEMP=A[365];
else
TEMP=A[366]+(N[366]>>1);
Q=RItype+365;
AQ=A[Q];
A[Q]=TEMP;
k=LG(Q);
if(cnt<(LIMIT-J[RUNindex]-qbpp-2))
{
if(k<=size)
{
output=decode>>(size-k);
if(size==k)
decode=0;
else
decode=((decode<<(32-size+k)))>>(32-size+k);
size=size-k;
}
else
{
fread(&buffer,1,1,fp);
decode1=0;
if(k-size<=8)
decode1=decode<<(k-size);
while(k-size>8)
{
decode1=decode1+(decode<<(k-size))+(buffer<<(k-size-8));
k=k-8;
fread(&buffer,1,1,fp);
decode=0;
}
output=decode1+(buffer>>(8+size-k));
if(8-k+size==0)
decode=0;
else
decode=(buffer<<(32-8+k-size))>>(32-8+k-size);
size=8-k+size;
}
EMErrval=(cnt<<k)+output;
}
else
{
if(qbpp<=size)
{
output=decode>>(size-qbpp);
if(size-qbpp==0)
decode=0;
else
decode=((decode<<(32-size+qbpp)))>>(32-size+qbpp);
size=size-qbpp;
}
else
{
fread(&buffer,1,1,fp);
decode1=0;
if(qbpp-size<=8)
decode1=decode<<(qbpp-size);
while(qbpp-size>8)
{
decode1=decode1+(decode<<(qbpp-size))+(buffer<<(qbpp-size-8));
qbpp=qbpp-8;
fread(&buffer,1,1,fp);
decode=0;
}
output=decode1+(buffer>>(8+size-qbpp));
if(8-qbpp+size==0)
decode=0;
else
decode=(buffer<<(32-8+qbpp-size))>>(32-8+qbpp-size);
size=8-qbpp+size;
}
EMErrval=output+1;
}
// **************************************
// * 由EMErrval求解Errval *
// **************************************
if(((EMErrval&1)==0)&&(RItype==1))
map=1;
else if(((EMErrval&1)==1)&&(RItype==1))
map=0;
else if(((EMErrval&1)==1)&&(RItype==0))
map=1;
else if(((EMErrval&1)==0)&&(RItype==0))
map=0;
if((map==0&&k!=0)||(map==0&&k==0&&2*Nn[Q-365]>=N[Q])||(map==1&&k==0&&2*Nn[Q-365]<N[Q]))
Errval=(EMErrval+RItype+map)>>1;
else
Errval=-((EMErrval+RItype+map)>>1);
flag=Errval;
Errval=Errval*(2*NEAR+1);
if(SIGN==-1)
Errval=-Errval;
Rx=Errval+Px;
if(Rx<(-NEAR))
Rx=Rx+RANGE*(2*NEAR+1);
if(Rx>(MAXVAL+NEAR))
Rx=Rx-RANGE*(2*NEAR+1);
if(Rx<0)
Rx=0;
else if(Rx>MAXVAL)
Rx=MAXVAL;
*(f+LinX*(y+2)+RowX)=Rx;
// **********************************
// * 更新其变量 *
// **********************************
if(flag<0)
Nn[Q-365]=Nn[Q-365]+1;
A[Q]=AQ;
A[Q]=A[Q]+((EMErrval+1-RItype)>>1);
if(N[Q]==RESET)
{
A[Q]=A[Q]>>1;
N[Q]=N[Q]>>1;
Nn[Q-365]=Nn[Q-365]>>1;
}
N[Q]=N[Q]+1;
}
// *************************************************
// * 主程序 *
// *************************************************
// i,j:循环变量计数；
// _______________________
// | | Rc| Rb| Rd| |
// -----------------------
// | | Ra| Ix| | |
// -----------------------
// D1,D2,D3:当前梯度
// x,y:图像的长度和宽度
void main()
{
int x,y,i,j;
int Ra,Rb,Rc,Rd,D1,D2,D3;
int BASIC_T1=3,BASIC_T2=7,BASIC_T3=21;
float T1,T2,T3;
int qbpp,filemode;
struct timeb start_time,end_time;
int second_d;
//printf("图象宽度和高度\n");
//scanf("%d,%d",&x,&y);
x=y=128;
//printf("选择文件格式(0:Adata5.dat,1:Adata5.raw)\n");
//scanf("%d",&filemode);
filemode=1;
ftime(&start_time);
// *********************************
// * 计算各参数的值 *
// *********************************
RANGE=((int)((MAXVAL+2*NEAR)/(2*NEAR+1)))+1;
qbpp=-(int)(floor(-log(RANGE)/log(2.0)));
bpp=__max(2,-(int)(floor(-log(MAXVAL+1)/log(2.0))));
LIMIT=2*(bpp+__max(8,bpp));
if(MAXVAL>=128){
T1=CLAMP_1((float)((int)(__min(MAXVAL,4095)+128)/256)*(BASIC_T1-2)+2+3*NEAR);
T2=CLAMP_2((float)((int)(__min(MAXVAL,4095)+128)/256)*(BASIC_T2-3)+3+5*NEAR,T1);
T3=CLAMP_3((float)((int)(__min(MAXVAL,4095)+128)/256)*(BASIC_T3-4)+4+7*NEAR,T2);}
else{
T1=CLAMP_1((float)__max(2,BASIC_T1/(int)(256/(MAXVAL+1))+3*NEAR));
T2=CLAMP_2((float)__max(3,BASIC_T2/(int)(256/(MAXVAL+1))+5*NEAR),T1);
T3=CLAMP_3((float)__max(4,BASIC_T3/(int)(256/(MAXVAL+1))+7*NEAR),T2);
}
f=(int *)calloc((x+1)*(y+2),sizeof(int));
if((fp=fopen("E:\\felicia\\Adata2","rb"))==NULL)
{
printf("cannot open the file\n");
return;
}
// ******************************
// * 对各数组的初始化 *
// ******************************
for(j=0;j<365;j++)
{
A[j]=__max(2,(RANGE+(1<<5))/(1<<6));
N[j]=1;
B[j]=0;
C[j]=0;
}
A[365]=A[0];
A[366]=A[0];
N[365]=1;
N[366]=1;
// ************************
// * 读入变量 *
// ************************
for(i=0;i<(x+1)*(y+2);i++)
*(f+i)=0;
while(RowX<=y&&LinX<=x)
{
Ra=(*(f+LinX*(y+2)+RowX-1));
Rb=(*(f+(LinX-1)*(y+2)+RowX));
Rc=(*(f+(LinX-1)*(y+2)+RowX-1));
Rd=(*(f+(LinX-1)*(y+2)+RowX+1));
D1=Rd-Rb;
D2=Rb-Rc;
D3=Rc-Ra;
// ******************************
// * 选择解码方式 *
// ******************************
if((abs(D1)<=NEAR)&&(abs(D2)<=NEAR)&&(abs(D3)<=NEAR))
RunModeProcessing(qbpp,Ra,Rb,y);
else
RegularModeProcessing(qbpp,Ra,Rb,Rc,Rd,D1,D2,D3,y,T1,T2,T3);
if(RowX==y){
*(f+LinX*(y+2)+y+1)=*(f+LinX*(y+2)+y);
*(f+(LinX+1)*(y+2))=*(f+LinX*(y+2)+1);}
RowX++;
if(RowX>y)
{
RowX=RowX-y;
LinX++;
}
}
// printf("max_k=%d\n",max_k);
fclose(fp);
// *************************************
// * 文件输出格式选择 *
// *************************************
if(filemode==0)
{
if((fp=fopen("E:\\felicia\\Adata5.dat","wt"))==NULL)
{
printf("cannot open the file\n");
exit(0);
}
for(i=1;i<x+1;i++)
{
for(j=1;j<y+1;j++)
fprintf(fp,"%5d",*(f+i*(y+2)+j));
fprintf(fp,"\n");
}
}
else
{
if((fp=fopen("E:\\felicia\\Adata5.raw","wb"))==NULL)
{
printf("cannot open the file\n");
exit(0);
}
for(i=1;i<x+1;i++)
for(j=1;j<y+1;j++)
fprintf(fp,"%c",(unsigned char)*(f+i*(y+2)+j));
//fwrite((f+i*(y+2)+j),2,1,fp);
}
fclose(fp);
// **************************
// * 计算编码时间 *
// **************************
ftime(&end_time);
second_d=end_time.millitm-start_time.millitm;
second_d=second_d+(end_time.time-start_time.time)*1000;
printf("The encoding costs:%.3f seconds.\n",(float)second_d/1000.0);
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -