二次线性插值的缩放函数

/


void ResizeImage(unsigned char *src,int nSrcWidth, int nSrcHeight, unsigned char *dst,int nDstWidth, int nDstHeight,int nChannel)
{


int step1= nSrcWidth*nChannel;
int step2= nDstWidth*nChannel;
int j;
int i;
int k;


unsigned char* data1= (unsigned char*)src;
unsigned char* data2= (unsigned char*)dst;


long xrIntFloat_16=((nSrcWidth)<<16)/nDstWidth+1; 
long yrIntFloat_16=((nSrcHeight)<<16)/nDstHeight+1;
const long csDErrorX=-(1<<15)+(xrIntFloat_16>>1);
const long csDErrorY=-(1<<15)+(yrIntFloat_16>>1);


long srcy_16=csDErrorY;


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


long srcx_16=csDErrorX;
unsigned long v_8=(srcy_16 & 0xFFFF)>>8;
long kq=srcy_16>>16;
kq = max_int(0,kq);
kq = min_int(kq,nSrcHeight-2);


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


unsigned long u_8=(srcx_16 & 0xFFFF)>>8;
unsigned long pm3_16=(u_8*v_8);
unsigned long pm2_16=(u_8*(unsigned long)(256-v_8));
unsigned long pm1_16=(v_8*(unsigned long)(256-u_8));
unsigned long pm0_16=((256-u_8)*(256-v_8));


long kp=srcx_16>>16;
kp = max_int(0,kp);
kp = min_int(kp,nSrcWidth-2);


for (k =0;k<3;k++)
{
data2[j*step2+i*3+k] =(pm0_16*data1[kq*step1+kp*3+k] + pm1_16*data1[(kq+1)*step1+kp*3+k] +
pm2_16*data1[kq*step1+(kp+1)*3+k] +pm3_16*data1[(kq+1)*step1+(kp+1)*3+k])>>16 ;


}
srcx_16+=xrIntFloat_16;


}
srcy_16+=yrIntFloat_16;


}


}