图形图像处理－之－高质量的快速的图像缩放中篇二次线性插值和三次卷积插值

墨蓝 2022-06-12 11:47 115阅读 0赞

转自 [http://blog.csdn.net/housisong/article/details/1452249][http_blog.csdn.net_housisong_article_details_1452249]

图形图像处理－之－高质量的快速的图像缩放 中篇 二次线性插值和三次卷积插值  
                     [HouSisong@GMail.com][HouSisong_GMail.com]   2006255.12.13

(2009.03.07  可以到这里下载缩放算法的完整的可以编译的项目源代码:  [http://blog.csdn.net/housisong/archive/2009/03/07/3967270.aspx][http_blog.csdn.net_housisong_archive_2009_03_07_3967270.aspx]  )

(2007.11.12 替换了二次线性插值的实现(以前偷懒使用了一个近似公式），改进后在图片边缘的插值效果更好(包括三次卷积插值的边界也更精确)；   
(2007.09.14 修正三次卷积的MMX版本中表的精度太低(7bit),造成卷积结果误差较大的问题,该版本提高了插值质量,并且速度加快12-25%)  
(2007.09.07 PicZoom\_ThreeOrder2和PicZoom\_ThreeOrder\_MMX在缩放的图片宽或高  
小于3个像素的时候有一个Bug(边界计算错误);将unsigned long xrIntFloat\_16,   
yrIntFloat\_16的定义改成long xrIntFloat\_16,yrIntFloat\_16就可以了)  
(2007.07.02 ThreeOrder2\_Fast一点小的改进,加快14%)  
(2007.06.18 优化PicZoom\_BilInear\_MMX的实现(由138.5fps提高到147.9fps),  
            并添加更快的两路展开的实现版本BilInear\_MMX\_expand2函数；  
            补充新的SSE2的实现PicZoom\_BilInear\_SSE2函数)  
(2007.06.06 更新测试数据,编译器由vc6改为vc2005,CPU由赛扬2G改为AMD64x2 4200+(2.1G) )  
(2007.03.06 更新)

tag:图像缩放,速度优化,定点数优化,近邻取样插值,二次线性插值,三次线性插值,  
   MipMap链,三次卷积插值,MMX,SSE,SSE2,CPU缓存优化

摘要:首先给出一个基本的图像缩放算法，然后一步一步的优化其速度和缩放质量；

高质量的快速的图像缩放 全文 分为:  
     上篇 近邻取样插值和其速度优化  
     中篇 二次线性插值和三次卷积插值  
     下篇 三次线性插值和MipMap链  
     补充 使用SSE2优化

正文：  
  为了便于讨论，这里只处理32bit的ARGB颜色；  
  代码使用C++;涉及到汇编优化的时候假定为x86平台;使用的编译器为vc2005;  
  为了代码的可读性,没有加入异常处理代码;  
   测试使用的CPU为AMD64x2 4200+(2.37G)  和 Intel Core2 4400(2.00G);

速度测试说明:  
  只测试内存数据到内存数据的缩放  
  测试图片都是800\*600缩放到1024\*768; fps表示每秒钟的帧数,值越大表示函数越快

A:近邻取样插值、二次线性插值、三次卷积插值 缩放**效果对比**

![zaka0.JPG][]             ![zaka1.JPG][]               ![zaka2.JPG][]    
       原图         近邻取样缩放到0.6倍     近邻取样缩放到1.6倍

![zoom2_0.6.JPG][]               ![zoom2_1.6.JPG][]   
                二次线性插值缩放到0.6倍   二次线性插值缩放到1.6倍

![zoom3_0.6.JPG][]                ![zoom3_1.6.JPG][]  
               三次卷积插值缩放到0.6倍   三次卷积插值缩放到1.6倍

![test0.GIF][] ![zoom0_8.GIF][] ![zoom2_true_8.JPG][] ![zoom3_8.JPG][] ![zoom2_8.JPG][]  
 原图 近邻取样缩放到8倍 二次线性插值缩放到8倍 三次卷积插值缩放到8倍 二次线性插值(近似公式)

近邻取样插值缩放简单、速度快，但很多时候缩放出的图片质量比较差(特别是对于人物、景色等),  
图片的缩放有比较明显的锯齿；使用二次或更高次插值有利于改善缩放效果；

B: 首先定义图像数据结构:

\#define asm \_\_asm  
  
typedef unsigned char TUInt8; // \[0..255\]  
struct TARGB32      //32 bit color  
\{  
    TUInt8  b,g,r,a;          //a is alpha  
\};  
  
struct TPicRegion  //一块颜色数据区的描述，便于参数传递  
\{  
    TARGB32\*    pdata;         //颜色数据首地址  
    long        byte\_width;    //一行数据的物理宽度(字节宽度)；  
                //abs(byte\_width)有可能大于等于width\*sizeof(TARGB32);  
    long        width;         //像素宽度  
    long        height;        //像素高度  
\};  
  
//那么访问一个点的函数可以写为：  
inline TARGB32& Pixels(const TPicRegion& pic,const long x,const long y)  
\{  
    return ( (TARGB32\*)((TUInt8\*)pic.pdata+pic.byte\_width\*y) )\[x\];  
\}

**二次线性插值缩放：**

C: 二次线性插值缩放原理和公式图示:

![zoom2_map0.gif][]

缩放后图片                 原图片  
            (宽DW,高DH)              (宽SW,高SH)

缩放映射原理:  
  (Sx-0)/(SW-0)=(Dx-0)/(DW-0)   (Sy-0)/(SH-0)=(Dy-0)/(DH-0)  
 =>   Sx=Dx\*SW/DW                    Sy=Dy\*SH/DH

聚焦看看(Sx,Sy)坐标点(Sx,Sy为浮点数)附近的情况；

![zoom2_map1.gif][]

对于近邻取样插值的缩放算法，直接取Color0颜色作为缩放后点的颜色；  
二次线性插值需要考虑(Sx,Sy)坐标点周围的4个颜色值Color0/Color1/Color2/Color3，  
把(Sx,Sy)到A/B/C/D坐标点的距离作为系数来把4个颜色混合出缩放后点的颜色；  
( u=Sx-floor(Sx); v=Sy-floor(Sy); 说明:floor函数的返回值为小于等于参数的最大整数 )    
  二次线性插值公式为:  
 tmpColor0=Color0\*(1-u) + Color2\*u;  
 tmpColor1=Color1\*(1-u) + Color3\*u;  
        DstColor =tmpColor0\*(1-v) + tmpColor2\*v;

展开公式为:  
        pm0=(1-u)\*(1-v);  
        pm1=v\*(1-u);  
        pm2=u\*(1-v);  
        pm3=u\*v;  
  则颜色混合公式为:  
        DstColor = Color0\*pm0 + Color1\*pm1 + Color2\*pm2 + Color3\*pm3;

参数函数图示：

![Bilinear.GIF][]

二次线性插值函数图示

对于上面的公式，它将图片向右下各移动了半个像素，需要对此做一个修正；  
  =>   Sx=(Dx+0.5)\*SW/DW-0.5; Sy=(Dy+0.5)\*SH/DH-0.5;  
而实际的程序,还需要考虑到边界(访问源图片可能超界)对于算法的影响,边界的处理可能有各种  
方案(不处理边界或边界回绕或边界饱和或边界映射或用背景颜色混合等;文章中默认使用边界饱和来处理超界);  
比如:边界饱和函数:

//访问一个点的函数，(x,y)坐标可能超出图片边界； //边界处理模式:边界饱和  
inline TARGB32 Pixels\_Bound(const TPicRegion& pic,long x,long y)  
\{  
    //assert((pic.width>0)&&(pic.height>0));  
    bool IsInPic=true;  
    if (x<0) \{x=0; IsInPic=false; \} else if (x>=pic.width ) \{x=pic.width -1; IsInPic=false; \}  
    if (y<0) \{y=0; IsInPic=false; \} else if (y>=pic.height) \{y=pic.height-1; IsInPic=false; \}  
    TARGB32 result=Pixels(pic,x,y);  
    if (!IsInPic) result.a=0;  
    return result;  
\}

D: 二次线性插值缩放算法的一个参考实现：PicZoom\_BilInear0  
  该函数并没有做什么优化，只是一个简单的浮点实现版本；

inline void Bilinear0(const TPicRegion& pic,float fx,float fy,TARGB32\* result)  
    \{  
        long x=(long)fx; if (x>fx) --x; //x=floor(fx);      
        long y=(long)fy; if (y>fy) --y; //y=floor(fy);  
          
        TARGB32 Color0=Pixels\_Bound(pic,x,y);  
        TARGB32 Color2=Pixels\_Bound(pic,x+1,y);  
        TARGB32 Color1=Pixels\_Bound(pic,x,y+1);  
        TARGB32 Color3=Pixels\_Bound(pic,x+1,y+1);  
  
        float u=fx-x;  
        float v=fy-y;  
        float pm3=u\*v;  
        float pm2=u\*(1-v);  
        float pm1=v\*(1-u);  
        float pm0=(1-u)\*(1-v);  
  
        result->a=(pm0\*Color0.a+pm1\*Color1.a+pm2\*Color2.a+pm3\*Color3.a);  
        result->r=(pm0\*Color0.r+pm1\*Color1.r+pm2\*Color2.r+pm3\*Color3.r);  
        result->g=(pm0\*Color0.g+pm1\*Color1.g+pm2\*Color2.g+pm3\*Color3.g);  
        result->b=(pm0\*Color0.b+pm1\*Color1.b+pm2\*Color2.b+pm3\*Color3.b);  
    \}  
  
void PicZoom\_Bilinear0(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    unsigned long dst\_width=Dst.width;  
    TARGB32\* pDstLine=Dst.pdata;  
    for (unsigned long y=0;y<Dst.height;++y)  
    \{  
        float srcy=(y+0.4999999)\*Src.height/Dst.height-0.5;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            float srcx=(x+0.4999999)\*Src.width/Dst.width-0.5;  
            Bilinear0(Src,srcx,srcy,&pDstLine\[x\]);  
        \}  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_BilInear0      8.3 fps

E: 把PicZoom\_BilInear0的浮点计算改写为定点数实现：PicZoom\_BilInear1

inline void Bilinear1(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        long x=x\_16>>16;  
        long y=y\_16>>16;  
        TARGB32 Color0=Pixels\_Bound(pic,x,y);  
        TARGB32 Color2=Pixels\_Bound(pic,x+1,y);  
        TARGB32 Color1=Pixels\_Bound(pic,x,y+1);  
        TARGB32 Color3=Pixels\_Bound(pic,x+1,y+1);  
  
        unsigned long u\_8=(x\_16 & 0xFFFF)>>8;  
        unsigned long v\_8=(y\_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));  
  
        result->a=((pm0\_16\*Color0.a+pm1\_16\*Color1.a+pm2\_16\*Color2.a+pm3\_16\*Color3.a)>>16);  
        result->r=((pm0\_16\*Color0.r+pm1\_16\*Color1.r+pm2\_16\*Color2.r+pm3\_16\*Color3.r)>>16);  
        result->g=((pm0\_16\*Color0.g+pm1\_16\*Color1.g+pm2\_16\*Color2.g+pm3\_16\*Color3.g)>>16);  
        result->b=((pm0\_16\*Color0.b+pm1\_16\*Color1.b+pm2\_16\*Color2.b+pm3\_16\*Color3.b)>>16);  
    \}  
  
void PicZoom\_Bilinear1(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear1(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_BilInear1     17.7 fps

F: 二次线性插值需要考略边界访问超界的问题，我们可以将边界区域和内部区域分开处理，这样就可以优化内部的插值实现函数了：比如不需要判断访问超界、减少颜色数据复制、减少一些不必要的重复坐标计算等等

inline void Bilinear2\_Fast(TARGB32\* PColor0,TARGB32\* PColor1,unsigned long u\_8,unsigned long v\_8,TARGB32\* result)  
    \{  
        unsigned long pm3\_16=u\_8\*v\_8;  
        unsigned long pm2\_16=(u\_8<<8)-pm3\_16;  
        unsigned long pm1\_16=(v\_8<<8)-pm3\_16;  
        unsigned long pm0\_16=(1<<16)-pm1\_16-pm2\_16-pm3\_16;  
     
        result->a=((pm0\_16\*PColor0\[0\].a+pm2\_16\*PColor0\[1\].a+pm1\_16\*PColor1\[0\].a+pm3\_16\*PColor1\[1\].a)>>16);  
        result->r=((pm0\_16\*PColor0\[0\].r+pm2\_16\*PColor0\[1\].r+pm1\_16\*PColor1\[0\].r+pm3\_16\*PColor1\[1\].r)>>16);  
        result->g=((pm0\_16\*PColor0\[0\].g+pm2\_16\*PColor0\[1\].g+pm1\_16\*PColor1\[0\].g+pm3\_16\*PColor1\[1\].g)>>16);  
        result->b=((pm0\_16\*PColor0\[0\].b+pm2\_16\*PColor0\[1\].b+pm1\_16\*PColor1\[0\].b+pm3\_16\*PColor1\[1\].b)>>16);  
    \}  
  
    inline void Bilinear2\_Border(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        long x=(x\_16>>16);  
        long y=(y\_16>>16);  
        unsigned long u\_16=((unsigned short)(x\_16));  
        unsigned long v\_16=((unsigned short)(y\_16));  
  
        TARGB32 pixel\[4\];  
        pixel\[0\]=Pixels\_Bound(pic,x,y);  
        pixel\[1\]=Pixels\_Bound(pic,x+1,y);  
        pixel\[2\]=Pixels\_Bound(pic,x,y+1);  
        pixel\[3\]=Pixels\_Bound(pic,x+1,y+1);  
          
        Bilinear2\_Fast(&pixel\[0\],&pixel\[2\],u\_16>>8,v\_16>>8,result);  
    \}  
  
void PicZoom\_Bilinear2(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    //计算出需要特殊处理的边界  
    long border\_y0=-csDErrorY/yrIntFloat\_16+1;              //y0+y\*yr>=0; y0=csDErrorY => y>=-csDErrorY/yr  
    if (border\_y0>=Dst.height) border\_y0=Dst.height;  
    long border\_x0=-csDErrorX/xrIntFloat\_16+1;       
    if (border\_x0>=Dst.width ) border\_x0=Dst.width;   
    long border\_y1=(((Src.height-2)<<16)-csDErrorY)/yrIntFloat\_16+1; //y0+y\*yr<=(height-2) => y<=(height-2-csDErrorY)/yr  
    if (border\_y1<border\_y0) border\_y1=border\_y0;  
    long border\_x1=(((Src.width-2)<<16)-csDErrorX)/xrIntFloat\_16+1;   
    if (border\_x1<border\_x0) border\_x1=border\_x0;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long Src\_byte\_width=Src.byte\_width;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<border\_y0;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear2\_Border(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y0;y<border\_y1;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        long x;  
        for (x=0;x<border\_x0;++x)  
        \{  
            Bilinear2\_Border(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
  
        \{  
            unsigned long v\_8=(srcy\_16 & 0xFFFF)>>8;  
            TARGB32\* PSrcLineColor= (TARGB32\*)((TUInt8\*)(Src.pdata)+Src\_byte\_width\*(srcy\_16>>16)) ;  
            for (unsigned long x=border\_x0;x<border\_x1;++x)  
            \{  
                TARGB32\* PColor0=&PSrcLineColor\[srcx\_16>>16\];  
                TARGB32\* PColor1=(TARGB32\*)((TUInt8\*)(PColor0)+Src\_byte\_width);  
                Bilinear2\_Fast(PColor0,PColor1,(srcx\_16 & 0xFFFF)>>8,v\_8,&pDstLine\[x\]);  
                srcx\_16+=xrIntFloat\_16;  
            \}  
        \}  
  
        for (x=border\_x1;x<dst\_width;++x)  
        \{  
            Bilinear2\_Border(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y1;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear2\_Border(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_BilInear2     43.4 fps

(F'补充:  
  如果不想处理边界访问超界问题，可以考虑扩大源图片的尺寸，加一个边框 (“哨兵”优化)；  
这样插值算法就不用考虑边界问题了，程序写起来也简单很多!   
  如果对缩放结果的边界像素级精度要求不是太高，我还有一个方案，一个稍微改变的缩放公式:  
  Sx=Dx\*(SW-1)/DW； Sy=Dy\*(SH-1)/DH；  (源图片宽和高:SW>=2;SH>=2)  
  证明这个公式不会造成内存访问超界:   
   要求Dx=DW-1时: sx+1=int( (dw-1)/dw\*(dw-1) ) +1 <= (sw-1)  
      有:  int( (sw-1)\*(dw-1)/dw ) <=sw-2  
             (sw-1)\*(dw-1)/dw <(sw-1)  
             (dw-1) /dw<1  
             (dw-1) <dw  
  比如，按这个公式的一个简单实现: (缩放效果见前面的"二次线性插值(近似公式)"图示)

void PicZoom\_ftBilinear\_Common(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(2>Src.width)||(2>Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width-1)<<16)/Dst.width;   
    long yrIntFloat\_16=((Src.height-1)<<16)/Dst.height;  
  
    unsigned long dst\_width=Dst.width;  
    long Src\_byte\_width=Src.byte\_width;  
    TARGB32\* pDstLine=Dst.pdata;  
    long srcy\_16=0;  
    for (unsigned long y=0;y<Dst.height;++y)  
    \{  
        unsigned long v\_8=(srcy\_16 & 0xFFFF)>>8;  
        TARGB32\* PSrcLineColor= (TARGB32\*)((TUInt8\*)(Src.pdata)+Src\_byte\_width\*(srcy\_16>>16)) ;  
        long srcx\_16=0;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            TARGB32\* PColor0=&PSrcLineColor\[srcx\_16>>16\];  
            Bilinear\_Fast\_Common(PColor0,(TARGB32\*)((TUInt8\*)(PColor0)+Src\_byte\_width),(srcx\_16 & 0xFFFF)>>8,v\_8,&pDstLine\[x\]);  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

)

G:利用单指令多数据处理的MMX指令一般都可以加快颜色的运算；在使用MMX改写之前，利用  
32bit寄存器(或变量)来模拟单指令多数据处理；  
数据储存原理:一个颜色数据分量只有一个字节,用2个字节来储存单个颜色分量的计算结果，  
对于很多颜色计算来说精度就够了；那么一个32bit寄存器(或变量)就可以储存2个计算出的  
临时颜色分量；从而达到了单个指令两路数据处理的目的;  
单个指令两路数据处理的计算:   
  乘法： ((0x00AA\*a)<<16) | (0x00BB\*a) = 0x00AA00BB \* a   
    可见只要保证0x00AA\*a和0x00BB\*a都小于(1<<16)那么乘法可以直接使用无符号数乘法了  
  加法: ((0x00AA+0x00CC)<<16) | (0x00BB+0x00DD) = 0x00AA00BB + 0x00CC00DD   
    可见只要0x00AA+0x00CC和0x00BB+0x00DD小于(1<<16)那么加法可以直接使用无符号数加法了  
  (移位、减法等稍微复杂一点,因为这里没有用到就不推倒运算公式了)

inline void Bilinear\_Fast\_Common(TARGB32\* PColor0,TARGB32\* PColor1,unsigned long u\_8,unsigned long v\_8,TARGB32\* result)  
    \{  
        unsigned long pm3\_8=(u\_8\*v\_8)>>8;  
        unsigned long pm2\_8=u\_8-pm3\_8;  
        unsigned long pm1\_8=v\_8-pm3\_8;  
        unsigned long pm0\_8=256-pm1\_8-pm2\_8-pm3\_8;  
  
        unsigned long Color=\*(unsigned long\*)(PColor0);  
        unsigned long BR=(Color & 0x00FF00FF)\*pm0\_8;  
        unsigned long GA=((Color & 0xFF00FF00)>>8)\*pm0\_8;  
                      Color=((unsigned long\*)(PColor0))\[1\];  
                      GA+=((Color & 0xFF00FF00)>>8)\*pm2\_8;  
                      BR+=(Color & 0x00FF00FF)\*pm2\_8;  
                      Color=\*(unsigned long\*)(PColor1);  
                      GA+=((Color & 0xFF00FF00)>>8)\*pm1\_8;  
                      BR+=(Color & 0x00FF00FF)\*pm1\_8;  
                      Color=((unsigned long\*)(PColor1))\[1\];  
                      GA+=((Color & 0xFF00FF00)>>8)\*pm3\_8;  
                      BR+=(Color & 0x00FF00FF)\*pm3\_8;  
  
        \*(unsigned long\*)(result)=(GA & 0xFF00FF00)|((BR & 0xFF00FF00)>>8);  
    \}  
  
    inline void Bilinear\_Border\_Common(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        long x=(x\_16>>16);  
        long y=(y\_16>>16);  
        unsigned long u\_16=((unsigned short)(x\_16));  
        unsigned long v\_16=((unsigned short)(y\_16));  
  
        TARGB32 pixel\[4\];  
        pixel\[0\]=Pixels\_Bound(pic,x,y);  
        pixel\[1\]=Pixels\_Bound(pic,x+1,y);  
        pixel\[2\]=Pixels\_Bound(pic,x,y+1);  
        pixel\[3\]=Pixels\_Bound(pic,x+1,y+1);  
          
        Bilinear\_Fast\_Common(&pixel\[0\],&pixel\[2\],u\_16>>8,v\_16>>8,result);  
    \}  
  
void PicZoom\_Bilinear\_Common(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    //计算出需要特殊处理的边界  
    long border\_y0=-csDErrorY/yrIntFloat\_16+1;              //y0+y\*yr>=0; y0=csDErrorY => y>=-csDErrorY/yr  
    if (border\_y0>=Dst.height) border\_y0=Dst.height;  
    long border\_x0=-csDErrorX/xrIntFloat\_16+1;       
    if (border\_x0>=Dst.width ) border\_x0=Dst.width;   
    long border\_y1=(((Src.height-2)<<16)-csDErrorY)/yrIntFloat\_16+1; //y0+y\*yr<=(height-2) => y<=(height-2-csDErrorY)/yr  
    if (border\_y1<border\_y0) border\_y1=border\_y0;  
    long border\_x1=(((Src.width-2)<<16)-csDErrorX)/xrIntFloat\_16+1;   
    if (border\_x1<border\_x0) border\_x1=border\_x0;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long Src\_byte\_width=Src.byte\_width;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<border\_y0;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y0;y<border\_y1;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        long x;  
        for (x=0;x<border\_x0;++x)  
        \{  
            Bilinear\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
  
        \{  
            unsigned long v\_8=(srcy\_16 & 0xFFFF)>>8;  
            TARGB32\* PSrcLineColor= (TARGB32\*)((TUInt8\*)(Src.pdata)+Src\_byte\_width\*(srcy\_16>>16)) ;  
            for (unsigned long x=border\_x0;x<border\_x1;++x)  
            \{  
                TARGB32\* PColor0=&PSrcLineColor\[srcx\_16>>16\];  
                TARGB32\* PColor1=(TARGB32\*)((TUInt8\*)(PColor0)+Src\_byte\_width);  
                Bilinear\_Fast\_Common(PColor0,PColor1,(srcx\_16 & 0xFFFF)>>8,v\_8,&pDstLine\[x\]);  
                srcx\_16+=xrIntFloat\_16;  
            \}  
        \}  
  
        for (x=border\_x1;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y1;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_BilInear\_Common   65.3 fps

H:使用MMX指令改写:PicZoom\_Bilinear\_MMX

inline void  Bilinear\_Fast\_MMX(TARGB32\* PColor0,TARGB32\* PColor1,unsigned long u\_8,unsigned long v\_8,TARGB32\* result)  
    \{  
        asm  
        \{      
              MOVD      MM6,v\_8  
              MOVD      MM5,u\_8  
              mov       edx,PColor0  
              mov       eax,PColor1  
              PXOR      mm7,mm7  
  
              MOVD         MM2,dword ptr \[eax\]    
              MOVD         MM0,dword ptr \[eax+4\]  
              PUNPCKLWD    MM5,MM5  
              PUNPCKLWD    MM6,MM6  
              MOVD         MM3,dword ptr \[edx\]    
              MOVD         MM1,dword ptr \[edx+4\]  
              PUNPCKLDQ    MM5,MM5   
              PUNPCKLBW    MM0,MM7  
              PUNPCKLBW    MM1,MM7  
              PUNPCKLBW    MM2,MM7  
              PUNPCKLBW    MM3,MM7  
              PSUBw        MM0,MM2  
              PSUBw        MM1,MM3  
              PSLLw        MM2,8  
              PSLLw        MM3,8  
              PMULlw       MM0,MM5  
              PMULlw       MM1,MM5  
              PUNPCKLDQ    MM6,MM6   
              PADDw        MM0,MM2  
              PADDw        MM1,MM3  
  
              PSRLw        MM0,8  
              PSRLw        MM1,8  
              PSUBw        MM0,MM1  
              PSLLw        MM1,8  
              PMULlw       MM0,MM6  
              mov       eax,result  
              PADDw        MM0,MM1  
  
              PSRLw        MM0,8  
              PACKUSwb     MM0,MM7  
              movd      \[eax\],MM0   
              //emms  
        \}  
    \}  
  
    void Bilinear\_Border\_MMX(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        long x=(x\_16>>16);  
        long y=(y\_16>>16);  
        unsigned long u\_16=((unsigned short)(x\_16));  
        unsigned long v\_16=((unsigned short)(y\_16));  
  
        TARGB32 pixel\[4\];  
        pixel\[0\]=Pixels\_Bound(pic,x,y);  
        pixel\[1\]=Pixels\_Bound(pic,x+1,y);  
        pixel\[2\]=Pixels\_Bound(pic,x,y+1);  
        pixel\[3\]=Pixels\_Bound(pic,x+1,y+1);  
          
        Bilinear\_Fast\_MMX(&pixel\[0\],&pixel\[2\],u\_16>>8,v\_16>>8,result);  
    \}  
  
void PicZoom\_Bilinear\_MMX(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    //计算出需要特殊处理的边界  
    long border\_y0=-csDErrorY/yrIntFloat\_16+1;              //y0+y\*yr>=0; y0=csDErrorY => y>=-csDErrorY/yr  
    if (border\_y0>=Dst.height) border\_y0=Dst.height;  
    long border\_x0=-csDErrorX/xrIntFloat\_16+1;       
    if (border\_x0>=Dst.width ) border\_x0=Dst.width;   
    long border\_y1=(((Src.height-2)<<16)-csDErrorY)/yrIntFloat\_16+1; //y0+y\*yr<=(height-2) => y<=(height-2-csDErrorY)/yr  
    if (border\_y1<border\_y0) border\_y1=border\_y0;  
    long border\_x1=(((Src.width-2)<<16)-csDErrorX)/xrIntFloat\_16+1;   
    if (border\_x1<border\_x0) border\_x1=border\_x0;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long Src\_byte\_width=Src.byte\_width;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<border\_y0;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y0;y<border\_y1;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        long x;  
        for (x=0;x<border\_x0;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
  
        \{  
            unsigned long v\_8=(srcy\_16 & 0xFFFF)>>8;  
            TARGB32\* PSrcLineColor= (TARGB32\*)((TUInt8\*)(Src.pdata)+Src\_byte\_width\*(srcy\_16>>16)) ;  
            for (unsigned long x=border\_x0;x<border\_x1;++x)  
            \{  
                TARGB32\* PColor0=&PSrcLineColor\[srcx\_16>>16\];  
                TARGB32\* PColor1=(TARGB32\*)((TUInt8\*)(PColor0)+Src\_byte\_width);  
                Bilinear\_Fast\_MMX(PColor0,PColor1,(srcx\_16 & 0xFFFF)>>8,v\_8,&pDstLine\[x\]);  
                srcx\_16+=xrIntFloat\_16;  
            \}  
        \}  
  
        for (x=border\_x1;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y1;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    asm emms  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_BilInear\_MMX 132.9 fps

H' 对BilInear\_MMX简单改进:PicZoom\_Bilinear\_MMX\_Ex

void PicZoom\_Bilinear\_MMX\_Ex(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    //计算出需要特殊处理的边界  
    long border\_y0=-csDErrorY/yrIntFloat\_16+1;              //y0+y\*yr>=0; y0=csDErrorY => y>=-csDErrorY/yr  
    if (border\_y0>=Dst.height) border\_y0=Dst.height;  
    long border\_x0=-csDErrorX/xrIntFloat\_16+1;       
    if (border\_x0>=Dst.width ) border\_x0=Dst.width;   
    long border\_y1=(((Src.height-2)<<16)-csDErrorY)/yrIntFloat\_16+1; //y0+y\*yr<=(height-2) => y<=(height-2-csDErrorY)/yr  
    if (border\_y1<border\_y0) border\_y1=border\_y0;  
    long border\_x1=(((Src.width-2)<<16)-csDErrorX)/xrIntFloat\_16+1;   
    if (border\_x1<border\_x0) border\_x1=border\_x0;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long Src\_byte\_width=Src.byte\_width;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<border\_y0;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
  
    for (y=border\_y0;y<border\_y1;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        long x;  
        for (x=0;x<border\_x0;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
  
        \{  
            long dst\_width\_fast=border\_x1-border\_x0;  
            if (dst\_width\_fast>0)  
            \{  
                unsigned long v\_8=(srcy\_16 & 0xFFFF)>>8;  
                TARGB32\* PSrcLineColor= (TARGB32\*)((TUInt8\*)(Src.pdata)+Src\_byte\_width\*(srcy\_16>>16)) ;  
                TARGB32\* PSrcLineColorNext= (TARGB32\*)((TUInt8\*)(PSrcLineColor)+Src\_byte\_width) ;  
                TARGB32\* pDstLine\_Fast=&pDstLine\[border\_x0\];  
                asm  
                \{  
                      movd         mm6,v\_8  
                      pxor         mm7,mm7 //mm7=0  
                      PUNPCKLWD    MM6,MM6  
                      PUNPCKLDQ    MM6,MM6//mm6=v\_8  
                      
                      mov       esi,PSrcLineColor  
                      mov       ecx,PSrcLineColorNext  
                      mov       edx,srcx\_16  
                      mov       ebx,dst\_width\_fast  
                      mov       edi,pDstLine\_Fast  
                      lea       edi,\[edi+ebx\*4\]  
                      push      ebp  
                      mov       ebp,xrIntFloat\_16  
                      neg       ebx  
  
                loop\_start:  
  
                          mov       eax,edx  
                          shl       eax,16  
                          shr       eax,24  
                          //== movzx       eax,dh  //eax=u\_8  
                          MOVD      MM5,eax  
                          mov       eax,edx  
                          shr       eax,16     //srcx\_16>>16  
  
                          MOVD         MM2,dword ptr \[ecx+eax\*4\]    
                          MOVD         MM0,dword ptr \[ecx+eax\*4+4\]  
                          PUNPCKLWD    MM5,MM5  
                          MOVD         MM3,dword ptr \[esi+eax\*4\]    
                          MOVD         MM1,dword ptr \[esi+eax\*4+4\]  
                          PUNPCKLDQ    MM5,MM5 //mm5=u\_8  
                          PUNPCKLBW    MM0,MM7  
                          PUNPCKLBW    MM1,MM7  
                          PUNPCKLBW    MM2,MM7  
                          PUNPCKLBW    MM3,MM7  
                          PSUBw        MM0,MM2  
                          PSUBw        MM1,MM3  
                          PSLLw        MM2,8  
                          PSLLw        MM3,8  
                          PMULlw       MM0,MM5  
                          PMULlw       MM1,MM5  
                          PADDw        MM0,MM2  
                          PADDw        MM1,MM3  
  
                          PSRLw        MM0,8  
                          PSRLw        MM1,8  
                          PSUBw        MM0,MM1  
                          PSLLw        MM1,8  
                          PMULlw       MM0,MM6  
                          PADDw        MM0,MM1  
  
                          PSRLw     MM0,8  
                          PACKUSwb  MM0,MM7  
                          MOVd   dword ptr    \[edi+ebx\*4\],MM0 //write DstColor  
                                        
                          add       edx,ebp //srcx\_16+=xrIntFloat\_16  
                          inc       ebx  
                          jnz       loop\_start  
  
                      pop       ebp  
                      mov       srcx\_16,edx  
                \}  
            \}  
        \}  
  
        for (x=border\_x1;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y1;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            Bilinear\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    asm emms  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_Bilinear\_MMX\_Ex 157.0 fps

I: 把测试成绩放在一起：

//CPU: AMD64x2 4200+(2.37G)  zoom 800\*600 to 1024\*768  
//==============================================================================  
// StretchBlt                   232.7 fps     
// PicZoom3\_SSE                 711.7 fps   
//   
// PicZoom\_BilInear0              8.3 fps  
// PicZoom\_BilInear1             17.7 fps  
// PicZoom\_BilInear2             43.4 fps  
// PicZoom\_BilInear\_Common       65.3 fps  
// PicZoom\_BilInear\_MMX         132.9 fps  
// PicZoom\_BilInear\_MMX\_Ex      157.0 fps

补充Intel Core2 4400上的测试成绩:

//CPU: Intel Core2 4400(2.00G)  zoom 800\*600 to 1024\*768  
//==============================================================================  
// PicZoom3\_SSE                1099.7 fps    
//   
// PicZoom\_BilInear0             10.7 fps  
// PicZoom\_BilInear1             24.2 fps  
// PicZoom\_BilInear2             54.3 fps  
// PicZoom\_BilInear\_Common       59.8 fps  
// PicZoom\_BilInear\_MMX         118.4 fps  
// PicZoom\_BilInear\_MMX\_Ex      142.9 fps

**三次卷积插值:**

J: 二次线性插值缩放出的图片很多时候让人感觉变得模糊(术语叫低通滤波)，特别是在放大  
的时候；使用三次卷积插值来改善插值结果；三次卷积插值考虑映射点周围16个点(4x4)的颜色来  
计算最终的混合颜色，如图；

![zoom4_map.gif][]  
         P(0,0)所在像素为映射的点，加上它周围的15个点，按一定系数混合得到最终输出结果；

混合公式参见PicZoom\_ThreeOrder0的实现；

插值曲线公式sin(x\*PI)/(x\*PI),如图:

![sinxdivx2.GIF][]  
                             三次卷积插值曲线sin(x\*PI)/(x\*PI) （其中PI=3.1415926...）

K:三次卷积插值缩放算法的一个参考实现：PicZoom\_ThreeOrder0  
  该函数并没有做过多的优化，只是一个简单的浮点实现版本;

inline double SinXDivX(double x)   
        \{  
            //该函数计算插值曲线sin(x\*PI)/(x\*PI)的值 //PI=3.1415926535897932385;   
            //下面是它的近似拟合表达式  
            const float a = -1; //a还可以取 a=-2,-1,-0.75,-0.5等等，起到调节锐化或模糊程度的作用  
  
            if (x<0) x=-x; //x=abs(x);  
            double x2=x\*x;  
            double x3=x2\*x;  
            if (x<=1)  
              return (a+2)\*x3 - (a+3)\*x2 + 1;  
            else if (x<=2)   
              return a\*x3 - (5\*a)\*x2 + (8\*a)\*x - (4\*a);  
            else  
              return 0;  
        \}   
  
        inline TUInt8 border\_color(long Color)  
        \{  
            if (Color<=0)  
                return 0;  
            else if (Color>=255)  
                return 255;  
            else  
                return Color;  
        \}  
          
    void ThreeOrder0(const TPicRegion& pic,const float fx,const float fy,TARGB32\* result)  
    \{  
        long x0=(long)fx; if (x0>fx) --x0; //x0=floor(fx);      
        long y0=(long)fy; if (y0>fy) --y0; //y0=floor(fy);  
        float fu=fx-x0;  
        float fv=fy-y0;  
  
        TARGB32 pixel\[16\];  
        long i,j;  
  
        for (i=0;i<4;++i)  
        \{  
            for (j=0;j<4;++j)  
            \{  
                long x=x0-1+j;  
                long y=y0-1+i;  
                pixel\[i\*4+j\]=Pixels\_Bound(pic,x,y);  
            \}  
        \}  
  
        float afu\[4\],afv\[4\];  
        //  
        afu\[0\]=SinXDivX(1+fu);  
        afu\[1\]=SinXDivX(fu);  
        afu\[2\]=SinXDivX(1-fu);  
        afu\[3\]=SinXDivX(2-fu);  
        afv\[0\]=SinXDivX(1+fv);  
        afv\[1\]=SinXDivX(fv);  
        afv\[2\]=SinXDivX(1-fv);  
        afv\[3\]=SinXDivX(2-fv);  
  
        float sR=0,sG=0,sB=0,sA=0;  
        for (i=0;i<4;++i)  
        \{  
            float aR=0,aG=0,aB=0,aA=0;  
            for (long j=0;j<4;++j)  
            \{  
                aA+=afu\[j\]\*pixel\[i\*4+j\].a;  
                aR+=afu\[j\]\*pixel\[i\*4+j\].r;  
                aG+=afu\[j\]\*pixel\[i\*4+j\].g;  
                aB+=afu\[j\]\*pixel\[i\*4+j\].b;  
            \}  
            sA+=aA\*afv\[i\];  
            sR+=aR\*afv\[i\];  
            sG+=aG\*afv\[i\];  
            sB+=aB\*afv\[i\];  
        \}  
  
        result->a=border\_color((long)(sA+0.5));  
        result->r=border\_color((long)(sR+0.5));  
        result->g=border\_color((long)(sG+0.5));  
        result->b=border\_color((long)(sB+0.5));  
    \}  
  
void PicZoom\_ThreeOrder0(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
  
    unsigned long dst\_width=Dst.width;  
    TARGB32\* pDstLine=Dst.pdata;  
    for (unsigned long y=0;y<Dst.height;++y)  
    \{  
        float srcy=(y+0.4999999)\*Src.height/Dst.height-0.5;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            float srcx=(x+0.4999999)\*Src.width/Dst.width-0.5;  
            ThreeOrder0(Src,srcx,srcy,&pDstLine\[x\]);  
        \}  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_ThreeOrder0    3.6 fps

L: 使用定点数来优化缩放函数;边界和内部分开处理；对SinXDivX做一个查找表;对border\_color做一个查找表；

static long SinXDivX\_Table\_8\[(2<<8)+1\];  
    class \_CAutoInti\_SinXDivX\_Table \{  
    private:   
        void \_Inti\_SinXDivX\_Table()  
        \{  
            for (long i=0;i<=(2<<8);++i)  
                SinXDivX\_Table\_8\[i\]=long(0.5+256\*SinXDivX(i\*(1.0/(256))))\*1;  
        \};  
    public:  
        \_CAutoInti\_SinXDivX\_Table() \{ \_Inti\_SinXDivX\_Table(); \}  
    \};  
    static \_CAutoInti\_SinXDivX\_Table \_\_tmp\_CAutoInti\_SinXDivX\_Table;  
  
  
    //颜色查表  
    static TUInt8 \_color\_table\[256\*3\];  
    static const TUInt8\* color\_table=&\_color\_table\[256\];  
    class \_CAuto\_inti\_color\_table  
    \{  
    public:  
        \_CAuto\_inti\_color\_table() \{  
            for (int i=0;i<256\*3;++i)  
                \_color\_table\[i\]=border\_color(i-256);  
        \}  
    \};  
    static \_CAuto\_inti\_color\_table \_Auto\_inti\_color\_table;  
  
    void ThreeOrder\_Fast\_Common(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        unsigned long u\_8=(unsigned char)((x\_16)>>8);  
        unsigned long v\_8=(unsigned char)((y\_16)>>8);  
        const TARGB32\* pixel=&Pixels(pic,(x\_16>>16)-1,(y\_16>>16)-1);  
        long pic\_byte\_width=pic.byte\_width;  
  
        long au\_8\[4\],av\_8\[4\];  
        //  
        au\_8\[0\]=SinXDivX\_Table\_8\[(1<<8)+u\_8\];  
        au\_8\[1\]=SinXDivX\_Table\_8\[u\_8\];  
        au\_8\[2\]=SinXDivX\_Table\_8\[(1<<8)-u\_8\];  
        au\_8\[3\]=SinXDivX\_Table\_8\[(2<<8)-u\_8\];  
        av\_8\[0\]=SinXDivX\_Table\_8\[(1<<8)+v\_8\];  
        av\_8\[1\]=SinXDivX\_Table\_8\[v\_8\];  
        av\_8\[2\]=SinXDivX\_Table\_8\[(1<<8)-v\_8\];  
        av\_8\[3\]=SinXDivX\_Table\_8\[(2<<8)-v\_8\];  
  
        long sR=0,sG=0,sB=0,sA=0;  
        for (long i=0;i<4;++i)  
        \{  
            long aA=au\_8\[0\]\*pixel\[0\].a + au\_8\[1\]\*pixel\[1\].a + au\_8\[2\]\*pixel\[2\].a + au\_8\[3\]\*pixel\[3\].a;  
            long aR=au\_8\[0\]\*pixel\[0\].r + au\_8\[1\]\*pixel\[1\].r + au\_8\[2\]\*pixel\[2\].r + au\_8\[3\]\*pixel\[3\].r;  
            long aG=au\_8\[0\]\*pixel\[0\].g + au\_8\[1\]\*pixel\[1\].g + au\_8\[2\]\*pixel\[2\].g + au\_8\[3\]\*pixel\[3\].g;  
            long aB=au\_8\[0\]\*pixel\[0\].b + au\_8\[1\]\*pixel\[1\].b + au\_8\[2\]\*pixel\[2\].b + au\_8\[3\]\*pixel\[3\].b;  
            sA+=aA\*av\_8\[i\];  
            sR+=aR\*av\_8\[i\];  
            sG+=aG\*av\_8\[i\];  
            sB+=aB\*av\_8\[i\];  
            ((TUInt8\*&)pixel)+=pic\_byte\_width;  
        \}  
  
        result->a=color\_table\[sA>>16\];  
        result->r=color\_table\[sR>>16\];  
        result->g=color\_table\[sG>>16\];  
        result->b=color\_table\[sB>>16\];  
    \}  
  
    void ThreeOrder\_Border\_Common(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        long x0\_sub1=(x\_16>>16)-1;  
        long y0\_sub1=(y\_16>>16)-1;  
        unsigned long u\_16\_add1=((unsigned short)(x\_16))+(1<<16);  
        unsigned long v\_16\_add1=((unsigned short)(y\_16))+(1<<16);  
  
        TARGB32 pixel\[16\];  
        long i;  
  
        for (i=0;i<4;++i)  
        \{  
            long y=y0\_sub1+i;  
            pixel\[i\*4+0\]=Pixels\_Bound(pic,x0\_sub1+0,y);  
            pixel\[i\*4+1\]=Pixels\_Bound(pic,x0\_sub1+1,y);  
            pixel\[i\*4+2\]=Pixels\_Bound(pic,x0\_sub1+2,y);  
            pixel\[i\*4+3\]=Pixels\_Bound(pic,x0\_sub1+3,y);  
        \}  
          
        TPicRegion npic;  
        npic.pdata     =&pixel\[0\];  
        npic.byte\_width=4\*sizeof(TARGB32);  
        //npic.width     =4;  
        //npic.height    =4;  
        ThreeOrder\_Fast\_Common(npic,u\_16\_add1,v\_16\_add1,result);  
    \}  
  
void PicZoom\_ThreeOrder\_Common(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    //计算出需要特殊处理的边界  
    long border\_y0=((1<<16)-csDErrorY)/yrIntFloat\_16+1;//y0+y\*yr>=1; y0=csDErrorY => y>=(1-csDErrorY)/yr  
    if (border\_y0>=Dst.height) border\_y0=Dst.height;  
    long border\_x0=((1<<16)-csDErrorX)/xrIntFloat\_16+1;  
    if (border\_x0>=Dst.width ) border\_x0=Dst.width;  
    long border\_y1=(((Src.height-3)<<16)-csDErrorY)/yrIntFloat\_16+1; //y0+y\*yr<=(height-3) => y<=(height-3-csDErrorY)/yr  
    if (border\_y1<border\_y0) border\_y1=border\_y0;  
    long border\_x1=(((Src.width-3)<<16)-csDErrorX)/xrIntFloat\_16+1;;   
    if (border\_x1<border\_x0) border\_x1=border\_x0;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<border\_y0;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            ThreeOrder\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y0;y<border\_y1;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        long x;  
        for (x=0;x<border\_x0;++x)  
        \{  
            ThreeOrder\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        for (x=border\_x0;x<border\_x1;++x)  
        \{  
            ThreeOrder\_Fast\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//fast  !  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        for (x=border\_x1;x<dst\_width;++x)  
        \{  
            ThreeOrder\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y1;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            ThreeOrder\_Border\_Common(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_ThreeOrder\_Common    16.9 fps

M: 用MMX来优化ThreeOrder\_Common函数:ThreeOrder\_MMX

typedef   unsigned long TMMXData32;  
    static TMMXData32 SinXDivX\_Table\_MMX\[(2<<8)+1\];  
    class \_CAutoInti\_SinXDivX\_Table\_MMX \{  
    private:   
        void \_Inti\_SinXDivX\_Table\_MMX()  
        \{  
            for (long i=0;i<=(2<<8);++i)  
            \{  
                unsigned short t=long(0.5+(1<<14)\*SinXDivX(i\*(1.0/(256))));  
                unsigned long tl=t | (((unsigned long)t)<<16);  
                SinXDivX\_Table\_MMX\[i\]=tl;  
            \}  
        \};  
    public:  
        \_CAutoInti\_SinXDivX\_Table\_MMX() \{ \_Inti\_SinXDivX\_Table\_MMX(); \}  
    \};  
    static \_CAutoInti\_SinXDivX\_Table\_MMX \_\_tmp\_CAutoInti\_SinXDivX\_Table\_MMX;  
  
  
    void \_\_declspec(naked) \_private\_ThreeOrder\_Fast\_MMX()  
    \{  
        asm  
        \{  
            movd        mm1,dword ptr \[edx\]  
            movd        mm2,dword ptr \[edx+4\]  
            movd        mm3,dword ptr \[edx+8\]  
            movd        mm4,dword ptr \[edx+12\]  
            movd        mm5,dword ptr \[(offset SinXDivX\_Table\_MMX)+256\*4+eax\*4\]  
            movd        mm6,dword ptr \[(offset SinXDivX\_Table\_MMX)+eax\*4\]  
            punpcklbw   mm1,mm7  
            punpcklbw   mm2,mm7  
            punpcklwd   mm5,mm5  
            punpcklwd   mm6,mm6  
            psllw       mm1,7  
            psllw       mm2,7  
            pmulhw      mm1,mm5  
            pmulhw      mm2,mm6  
            punpcklbw   mm3,mm7  
            punpcklbw   mm4,mm7  
            movd        mm5,dword ptr \[(offset SinXDivX\_Table\_MMX)+256\*4+ecx\*4\]  
            movd        mm6,dword ptr \[(offset SinXDivX\_Table\_MMX)+512\*4+ecx\*4\]  
            punpcklwd   mm5,mm5  
            punpcklwd   mm6,mm6  
            psllw       mm3,7  
            psllw       mm4,7  
            pmulhw      mm3,mm5  
            pmulhw      mm4,mm6  
            paddsw      mm1,mm2  
            paddsw      mm3,mm4  
            movd        mm6,dword ptr \[ebx\] //v  
            paddsw      mm1,mm3  
            punpcklwd   mm6,mm6  
  
            pmulhw      mm1,mm6  
            add     edx,esi  //\+pic.byte\_width  
            paddsw      mm0,mm1  
  
            ret  
        \}  
    \}  
  
    inline void ThreeOrder\_Fast\_MMX(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        asm  
        \{  
            mov     ecx,pic  
            mov     eax,y\_16  
            mov     ebx,x\_16  
            movzx   edi,ah //v\_8  
            mov     edx,\[ecx+TPicRegion::pdata\]  
            shr     eax,16  
            mov     esi,\[ecx+TPicRegion::byte\_width\]  
            dec     eax  
            movzx   ecx,bh //u\_8  
            shr     ebx,16  
            imul    eax,esi  
            lea     edx,\[edx+ebx\*4-4\]  
            add     edx,eax //pixel  
  
            mov     eax,ecx  
            neg     ecx  
  
            pxor    mm7,mm7  //0  
            //mov     edx,pixel  
            pxor    mm0,mm0  //result=0  
            //lea     eax,auv\_7  
  
            lea    ebx,\[(offset SinXDivX\_Table\_MMX)+256\*4+edi\*4\]  
            call  \_private\_ThreeOrder\_Fast\_MMX  
            lea    ebx,\[(offset SinXDivX\_Table\_MMX)+edi\*4\]  
            call  \_private\_ThreeOrder\_Fast\_MMX  
            neg    edi  
            lea    ebx,\[(offset SinXDivX\_Table\_MMX)+256\*4+edi\*4\]  
            call  \_private\_ThreeOrder\_Fast\_MMX  
            lea    ebx,\[(offset SinXDivX\_Table\_MMX)+512\*4+edi\*4\]  
            call  \_private\_ThreeOrder\_Fast\_MMX  
  
            psraw     mm0,3  
            mov       eax,result  
            packuswb  mm0,mm7  
            movd      \[eax\],mm0  
            //emms  
        \}  
    \}  
  
    void ThreeOrder\_Border\_MMX(const TPicRegion& pic,const long x\_16,const long y\_16,TARGB32\* result)  
    \{  
        unsigned long x0\_sub1=(x\_16>>16)-1;  
        unsigned long y0\_sub1=(y\_16>>16)-1;  
        long u\_16\_add1=((unsigned short)(x\_16))+(1<<16);  
        long v\_16\_add1=((unsigned short)(y\_16))+(1<<16);  
  
        TARGB32 pixel\[16\];  
  
        for (long i=0;i<4;++i)  
        \{  
            long y=y0\_sub1+i;  
            pixel\[i\*4+0\]=Pixels\_Bound(pic,x0\_sub1  ,y);  
            pixel\[i\*4+1\]=Pixels\_Bound(pic,x0\_sub1+1,y);  
            pixel\[i\*4+2\]=Pixels\_Bound(pic,x0\_sub1+2,y);  
            pixel\[i\*4+3\]=Pixels\_Bound(pic,x0\_sub1+3,y);  
        \}  
          
        TPicRegion npic;  
        npic.pdata     =&pixel\[0\];  
        npic.byte\_width=4\*sizeof(TARGB32);  
        //npic.width     =4;  
        //npic.height    =4;  
        ThreeOrder\_Fast\_MMX(npic,u\_16\_add1,v\_16\_add1,result);  
    \}  
  
void PicZoom\_ThreeOrder\_MMX(const TPicRegion& Dst,const TPicRegion& Src)  
\{  
    if (  (0==Dst.width)||(0==Dst.height)  
        ||(0==Src.width)||(0==Src.height)) return;  
  
    long xrIntFloat\_16=((Src.width)<<16)/Dst.width+1;   
    long yrIntFloat\_16=((Src.height)<<16)/Dst.height+1;  
    const long csDErrorX=-(1<<15)+(xrIntFloat\_16>>1);  
    const long csDErrorY=-(1<<15)+(yrIntFloat\_16>>1);  
  
    unsigned long dst\_width=Dst.width;  
  
    //计算出需要特殊处理的边界  
    long border\_y0=((1<<16)-csDErrorY)/yrIntFloat\_16+1;//y0+y\*yr>=1; y0=csDErrorY => y>=(1-csDErrorY)/yr  
    if (border\_y0>=Dst.height) border\_y0=Dst.height;  
    long border\_x0=((1<<16)-csDErrorX)/xrIntFloat\_16+1;  
    if (border\_x0>=Dst.width ) border\_x0=Dst.width;  
    long border\_y1=(((Src.height-3)<<16)-csDErrorY)/yrIntFloat\_16+1; //y0+y\*yr<=(height-3) => y<=(height-3-csDErrorY)/yr  
    if (border\_y1<border\_y0) border\_y1=border\_y0;  
    long border\_x1=(((Src.width-3)<<16)-csDErrorX)/xrIntFloat\_16+1;;   
    if (border\_x1<border\_x0) border\_x1=border\_x0;  
  
    TARGB32\* pDstLine=Dst.pdata;  
    long srcy\_16=csDErrorY;  
    long y;  
    for (y=0;y<border\_y0;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            ThreeOrder\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y0;y<border\_y1;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        long x;  
        for (x=0;x<border\_x0;++x)  
        \{  
            ThreeOrder\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        for (x=border\_x0;x<border\_x1;++x)  
        \{  
            ThreeOrder\_Fast\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//fast MMX !  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        for (x=border\_x1;x<dst\_width;++x)  
        \{  
            ThreeOrder\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]);//border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    for (y=border\_y1;y<Dst.height;++y)  
    \{  
        long srcx\_16=csDErrorX;  
        for (unsigned long x=0;x<dst\_width;++x)  
        \{  
            ThreeOrder\_Border\_MMX(Src,srcx\_16,srcy\_16,&pDstLine\[x\]); //border  
            srcx\_16+=xrIntFloat\_16;  
        \}  
        srcy\_16+=yrIntFloat\_16;  
        ((TUInt8\*&)pDstLine)+=Dst.byte\_width;  
    \}  
    asm emms  
\}

//速度测试:  
//==============================================================================  
// PicZoom\_ThreeOrder\_MMX   34.3 fps

N：将测试结果放到一起：

//CPU: AMD64x2 4200+(2.37G)  zoom 800\*600 to 1024\*768  
//==============================================================================  
// StretchBlt                   232.7 fps     
// PicZoom3\_SSE                 711.7 fps    
// PicZoom\_BilInear\_MMX\_Ex      157.0 fps  
//   
// PicZoom\_ThreeOrder0            3.6 fps  
// PicZoom\_ThreeOrder\_Common     16.9 fps  
// PicZoom\_ThreeOrder\_MMX        34.3 fps

补充Intel Core2 4400上的测试成绩:

//CPU: Intel Core2 4400(2.00G)  zoom 800\*600 to 1024\*768  
//==============================================================================  
// PicZoom3\_SSE                1099.7 fps    
// PicZoom\_BilInear\_MMX\_Ex      142.9 fps  
//   
// PicZoom\_ThreeOrder0            4.2 fps  
// PicZoom\_ThreeOrder\_Common     17.6 fps  
// PicZoom\_ThreeOrder\_MMX        34.4 fps

[http_blog.csdn.net_housisong_article_details_1452249]: http://blog.csdn.net/housisong/article/details/1452249
[HouSisong_GMail.com]: mailto:HouSisong@263.net
[http_blog.csdn.net_housisong_archive_2009_03_07_3967270.aspx]: http://blog.csdn.net/housisong/archive/2009/03/07/3967270.aspx
[zaka0.JPG]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zaka0.JPG
[zaka1.JPG]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zaka1.JPG
[zaka2.JPG]: /images/20220612/59bce493868a4f31b9e0d3d4f5aaf146.png
[zoom2_0.6.JPG]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zoom2_0.6.JPG
[zoom2_1.6.JPG]: /images/20220612/6be0980caa464a659a25ef94741164e8.png
[zoom3_0.6.JPG]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zoom3_0.6.JPG
[zoom3_1.6.JPG]: /images/20220612/6efe00525365480f8adc465ff2342bb7.png
[test0.GIF]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/test0.GIF
[zoom0_8.GIF]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zoom0_8.GIF
[zoom2_true_8.JPG]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zoom2_true_8.JPG
[zoom3_8.JPG]: https://p-blog.csdn.net/images/p_blog_csdn_net/housisong/zoom3_8.JPG
[zoom2_8.JPG]: /images/20220612/d9424dafc7c34c85aa800dd2b4171fe2.png
[zoom2_map0.gif]: /images/20220612/0573df6b52304a2bac783dd9bf939a1c.png
[zoom2_map1.gif]: /images/20220612/b0281c3bdc21497ebe6df4b8e465d407.png
[Bilinear.GIF]: /images/20220612/d8385e1fe6ec42eabfd033e6cdea2bce.png
[zoom4_map.gif]: /images/20220612/55266d80f98f426aaf92362b7f098826.png
[sinxdivx2.GIF]: /images/20220612/aaceb9d8ee91462f8e180714cbb9af98.png