相位插值

淩亂°似流年 2021-11-18 01:08 349阅读 0赞

1、原旋转中的

![ContractedBlock.gif][] ![ExpandedBlockStart.gif][]

double phase[33]={
         0};//16相位，包含端点存在33个距离
        for (int i=0;i<33;i++)
        {
            double i2=1.0*i;
            phase[i]=fs(i2/16);
        }
        int size;
        if(dbZoom>1)
        {
            size=iWidth*iHeight;
        }
        else
        {
            size=ceil(iWidth*dbZoom)*ceil(iHeight*dbZoom);
        }
        int newWidth=ceil(dbZoom*iWidth);
        //旋转中心为图像中心
        double rx0=ceil(dbZoom*iWidth)*0.5;  
        double ry0=ceil(dbZoom*iHeight)*0.5; 
        double srcx,srcy,u,v;
        int xOr,yOr;
        dbRotate=dbRotate*3.1415926/180.0;
        for (int y=0;y<ceil(dbZoom*iHeight);y++)
        {
            for (int x=0;x<ceil(dbZoom*iWidth);x++)
            {
                srcx=(double)((x-rx0)*cos(dbRotate) - (y-ry0)*sin(dbRotate) + rx0) ;
                srcy=(double)((x-rx0)*sin(dbRotate) + (y-ry0)*cos(dbRotate) + ry0) ;
                srcx=srcx*1/dbZoom;
                srcy=srcy*1/dbZoom;
                xOr = floor(srcx);
                yOr = floor(srcy);
                u=srcx-xOr;
                v=srcy-yOr;
                int phasex=floor(16*u+0.5);//16相位
                int phasey=floor(16*v+0.5);
                double A1,B1,C1,D1,A2,B2,C2,D2;
                A1=phase[16+phasex];
                B1=phase[phasex];
                C1=phase[16-phasex];
                D1=phase[32-phasex];
    
                A2=phase[16+phasey];
                B2=phase[phasey];
                C2=phase[16-phasey];
                D2=phase[32-phasey];
    
            /*    A1=fs(1+u);
                B1=fs(u);
                C1=fs(1-u);
                D1=fs(2-u);
    
                A2=fs(1+v);
                B2=fs(v);
                C2=fs(1-v);
                D2=fs(2-v);*/
    
                int newWidth=ceil(dbZoom*iWidth);
                if( !(xOr-1>=0 && xOr+2<=iWidth && yOr-1>=0 && yOr+2<=iHeight))
                    {
                        pbTag[y*newWidth+x]=0;
                    }
                else
                {   
                    double middle=
                        pbSrc[(yOr-1)*iWidth+(xOr-1)]*A1*A2+
                        pbSrc[(yOr)*iWidth+(xOr-1)]*A1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr-1)]*A1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr-1)]*A1*D2+
                        
                        pbSrc[(yOr-1)*iWidth+(xOr)]*B1*A2+
                        pbSrc[(yOr)*iWidth+(xOr)]*B1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr)]*B1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr)]*B1*D2+
    
                        pbSrc[(yOr-1)*iWidth+(xOr+1)]*C1*A2+
                        pbSrc[(yOr)*iWidth+(xOr+1)]*C1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr+1)]*C1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr+1)]*C1*D2+
    
                        pbSrc[(yOr-1)*iWidth+(xOr+2)]*D1*A2+
                        pbSrc[(yOr)*iWidth+(xOr+2)]*D1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr+2)]*D1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr+2)]*D1*D2;
    
                    if(middle<=255&&middle>=0)
                        pbTag[y*newWidth+x]=middle;
                    else if(middle>255)
                        pbTag[y*newWidth+x]=255;
                    else 
                        pbTag[y*newWidth+x]=0;
                }
            }
        }
        //memcpy(pbTag,pbSrc,size);
        return ERROR_SUCCESS;

2、插值后只显示中间1280\*1024。

![ContractedBlock.gif][] ![ExpandedBlockStart.gif][]

double phase[33]={
         0};//16相位，包含端点存在33个距离
        for (int i=0;i<33;i++)
        {
            double i2=1.0*i;
            phase[i]=fs(i2/16);
        }
    
        //旋转中心为图像中心
        double rx0=iWidth;  
        double ry0=iHeight; 
        double srcx,srcy,u,v;
        int xOr,yOr;
        int newWidth=ceil(dbZoom*iWidth);
        int newHeight=ceil(dbZoom*iHeight);
    
        for (int y=0;y<ceil(dbZoom*iHeight);y++)
        {
            for (int x=0;x<ceil(dbZoom*iWidth);x++)
            {
    
                srcx=(double)x;
                srcy=(double)y ;
                srcx=srcx*1/dbZoom;
                srcy=srcy*1/dbZoom;
                xOr = floor(srcx);
                yOr = floor(srcy);
                u=srcx-xOr;
                v=srcy-yOr;
    
                int phasex=floor(16*u+0.5);//16相位
                int phasey=floor(16*v+0.5);
                double A1,B1,C1,D1,A2,B2,C2,D2;
                A1=phase[16+phasex];
                B1=phase[phasex];
                C1=phase[16-phasex];
                D1=phase[32-phasex];
    
                A2=phase[16+phasey];
                B2=phase[phasey];
                C2=phase[16-phasey];
                D2=phase[32-phasey];
    
    
                if( !(srcx>=0 && srcx<=iWidth && srcy>=0 && srcy<=iHeight))
                    {
                        pbTag[y*2*iWidth+x]=0;//255
                    }
                else
                {   
                    double middle=
                        pbSrc[(yOr-1)*iWidth+(xOr-1)]*A1*A2+
                        pbSrc[(yOr)*iWidth+(xOr-1)]*A1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr-1)]*A1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr-1)]*A1*D2+
                        
                        pbSrc[(yOr-1)*iWidth+(xOr)]*B1*A2+
                        pbSrc[(yOr)*iWidth+(xOr)]*B1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr)]*B1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr)]*B1*D2+
    
                        pbSrc[(yOr-1)*iWidth+(xOr+1)]*C1*A2+
                        pbSrc[(yOr)*iWidth+(xOr+1)]*C1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr+1)]*C1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr+1)]*C1*D2+
    
                        pbSrc[(yOr-1)*iWidth+(xOr+2)]*D1*A2+
                        pbSrc[(yOr)*iWidth+(xOr+2)]*D1*B2+
                        pbSrc[(yOr+1)*iWidth+(xOr+2)]*D1*C2+
                        pbSrc[(yOr+2)*iWidth+(xOr+2)]*D1*D2;
    
                    if(middle<=255&&middle>=0)
                        pbTag[y*newWidth+x]=middle;
                    else if(middle>255)
                        pbTag[y*newWidth+x]=255;
                    else 
                        pbTag[y*newWidth+x]=0;
                }
            }
        }
    
        for (int y=0;y<1024;y++)
        {
            for (int x=0;x<1280;x++)
            {
                int y2=newHeight/2-512;
                int x2=newWidth/2-640;
                 pbFinal[y*1280+x]=pbTag[(y2+y)*newWidth+(x2+x)];
           }
        }
        return ERROR_SUCCESS;

双三。

转载于:https://www.cnblogs.com/wxl845235800/p/11152049.html

[ContractedBlock.gif]: https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif
[ExpandedBlockStart.gif]: /images/20211116/3180320133fd40d08d2f0eae65192a7c.png