20190710记录：去掉中转图，直接以1280*1024进行反坐标计算，填补pbFinal。

1、记录：去掉中转图，直接以1280*1024进行反坐标计算。pbFinal=1280*1024。

// Imagejoint.cpp : 定义控制台应用程序的入口点。

//

#include "stdafx.h"

#include "Imagejoint.h"

#include "math.h"

#include <afxwin.h>

#ifdef _DEBUG

#define new DEBUG_NEW

#endif

#include <atlimage.h>//CImage类

#include <locale.h>

// The one and only application object

CWinApp theApp;

using namespace std;

//双三次插值系数

double fs(double w)

{

        double a=-0.5;

        double fs;

        if (abs(w)<=)

            fs=(a+)*pow(abs(w),)-(a+)*pow(abs(w),)+;

        else if (abs(w)>&&abs(w)<=)

            fs=a*pow(abs(w),)-*a*pow(abs(w),)+*a*abs(w)-*a;

        else

            fs=;

        return fs;

}

HRESULT Imagejoint(PBYTE pbSrc,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal);

int _tmain(int argc, TCHAR* argv[], TCHAR* envp[])

{

    int nRetCode = ;//表示整数类型的函数返回码。n表示整数类型,Ret是Return的缩写,表示返回值,Code表示代码。

    setlocale(LC_ALL,"chs");

    HMODULE hModule = ::GetModuleHandle(NULL);

    if (hModule != NULL)

    {

        // initialize MFC and print and error on failure

        CImage cImage;

        HRESULT hResult;

        //初始化一些变量

        int        iWidth,iHeight,iBytePerPixel,iPitch;

        int        x,y;

        PBYTE    pbSrc=NULL,pbFinal=NULL;//源图、目标图

        PBYTE    pbImage=NULL;//load图像后存在这

        PDWORD    pdwImage=NULL;//用于保存图像

        double    dbZoom=2.25;

        CString str = "frame1.bmp";

        LPCTSTR filename = (LPCTSTR)str;

        do{

            if (!AfxWinInit(hModule, NULL, ::GetCommandLine(), ))

            {

                // TODO: change error code to suit your needs

                _tprintf(_T("Fatal Error: MFC initialization failed\n"));

                nRetCode = ;

                break;

            }

            //Load 图像到cImage对象中

            hResult=cImage.Load(filename);

            if(hResult!=ERROR_SUCCESS)

            {

                _tprintf(_T("源图像文件名错误！\n"));

                nRetCode= -;

                break;

            }

            iWidth=cImage.GetWidth();

            iHeight=cImage.GetHeight();

            //分配源图内存

            pbSrc = (PBYTE)malloc(iWidth*iHeight);

            pbFinal = (PBYTE)malloc(*);

        //while (dbZoom<3);

        //{

            //分配目标图内存

            if(pbSrc==NULL || pbFinal==NULL )

            {

                _tprintf(_T("内存申请错误！\n"));

                nRetCode= -;

                break;

            }

            //cImage数据存到pbImage，后再转换为源灰度图pbSrc

            iPitch=cImage.GetPitch();

            iBytePerPixel=(cImage.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage.Destroy();

            //TO DO:执行操作

            hResult=Imagejoint(pbSrc,iWidth,iHeight,dbZoom,pbFinal);

            if(hResult!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像处理错误！\n"));

                nRetCode= -;

                break;

            }

            //处理后保存图像

            iWidth=;

            iHeight=;

            cImage.Create(iWidth,-iHeight,);

            iPitch=cImage.GetPitch();

            for(y=;y<iHeight;y++)

            {

                pdwImage=(PDWORD)(PBYTE(cImage.GetBits())+iPitch*y);

                for(x=;x<iWidth;x++)

                {

                    pdwImage[x]=pbFinal[y*iWidth+x]*0x10101;

                }

            }

            //待保存图像文件名

            CString name1="target";

            CString name2;

            name2.Format(_T("%.2lf"),dbZoom);

            CString csTagName;

            csTagName=name1+name2;

            //CString csTagName="target";

            csTagName.Trim();

            csTagName.MakeUpper();

            if(csTagName.Right()!=_T(".BMP") ) csTagName.Append(_T(".BMP"));

            hResult=cImage.Save(csTagName);

            if(hResult!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像结果保存错误！\n"));

                nRetCode= -;

                break;

            }

            _tprintf(_T("图像处理成功！\n"));

            nRetCode= ERROR_SUCCESS;

            break;

            //dbZoom=dbZoom+0.5;

            //}

        }while();

        if(pbSrc) free(pbSrc);

        if(pbFinal) free(pbFinal);

    }

    else

    {

        _tprintf(_T("Fatal Error: GetModuleHandle failed\n"));

        nRetCode = ;

    }

    getchar();

    return nRetCode;

}

HRESULT Imagejoint(PBYTE pbSrc,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal)

{

    double phase[]={};//16相位，包含端点存在33个距离

    for (int i=;i<;i++)

    {

        double i2=1.0*i;

        phase[i]=fs(i2/);

    }

    //旋转中心为图像中心

    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=;y<;y++)

    {

        for (int x=;x<;x++)

        {

            srcx=(double)(newWidth/-+x)/dbZoom;

            srcy=(double)(newHeight/-+y)/dbZoom ;

            xOr = floor(srcx);

            yOr = floor(srcy);

            u=srcx-xOr;

            v=srcy-yOr;

            int phasex=floor(*u+0.5);//16相位

            int phasey=floor(*v+0.5);

            double A1,B1,C1,D1,A2,B2,C2,D2;

            A1=phase[+phasex];

            B1=phase[phasex];

            C1=phase[-phasex];

            D1=phase[-phasex];

            A2=phase[+phasey];

            B2=phase[phasey];

            C2=phase[-phasey];

            D2=phase[-phasey];

            if( !(srcx>= && srcx<=iWidth && srcy>= && srcy<=iHeight))

                {

                    pbFinal[y*+x]=;//

                }

            else

            {

                double middle=

                    pbSrc[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

            }

        }

    }

    return ERROR_SUCCESS;

}

2、外圈需拼接

目前dbZoom都是表示的，中间处的放大倍数。

由于目前测试图像是768*576的。

在dbZoom<1.75左右时，其外圈仍需补充拼接。

先进行简单测试：发现需补充拼接的点为，映射后落在原图外的部分。暂时仍采用同一图像进行外圈的拼接。

// Imagejoint.cpp : 定义控制台应用程序的入口点。

//

#include "stdafx.h"

#include "Imagejoint.h"

#include <afxwin.h>

#ifdef _DEBUG

#define new DEBUG_NEW

#endif

#include <atlimage.h>//CImage类

#include <locale.h>

#include "math.h"

using namespace std;

//双三次插值系数

double fs(double w)

{

        double a=-0.5;

        double fs;

        if (abs(w)<=)

            fs=(a+)*pow(abs(w),)-(a+)*pow(abs(w),)+;

        else if (abs(w)>&&abs(w)<=)

            fs=a*pow(abs(w),)-*a*pow(abs(w),)+*a*abs(w)-*a;

        else

            fs=;

        return fs;

}

HRESULT Imagejoint(PBYTE pbSrc,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal);

int _tmain(int argc, TCHAR* argv[], TCHAR* envp[])

{

    int nRetCode = ;//表示整数类型的函数返回码。n表示整数类型,Ret是Return的缩写,表示返回值,Code表示代码。

    setlocale(LC_ALL,"chs");

    HMODULE hModule = ::GetModuleHandle(NULL);

    if (hModule != NULL)

    {

        // initialize MFC and print and error on failure

        CImage cImage;

        HRESULT hResult;

        //初始化一些变量

        int        iWidth,iHeight,iBytePerPixel,iPitch;

        int        x,y;

        PBYTE    pbSrc=NULL,pbFinal=NULL;//源图、目标图

        PBYTE    pbImage=NULL;//load图像后存在这

        PDWORD    pdwImage=NULL;//用于保存图像

        double    dbZoom=1.25;

        CString str = "frame1.bmp";

        LPCTSTR filename = (LPCTSTR)str;

        do{

            if (!AfxWinInit(hModule, NULL, ::GetCommandLine(), ))

            {

                // TODO: change error code to suit your needs

                _tprintf(_T("Fatal Error: MFC initialization failed\n"));

                nRetCode = ;

                break;

            }

            //Load 图像到cImage对象中

            hResult=cImage.Load(filename);

            if(hResult!=ERROR_SUCCESS)

            {

                _tprintf(_T("源图像文件名错误！\n"));

                nRetCode= -;

                break;

            }

            iWidth=cImage.GetWidth();

            iHeight=cImage.GetHeight();

            //分配源图内存

            pbSrc = (PBYTE)malloc(iWidth*iHeight);

            pbFinal = (PBYTE)malloc(*);

        //while (dbZoom<3);

        //{

            //分配目标图内存

            if(pbSrc==NULL || pbFinal==NULL )

            {

                _tprintf(_T("内存申请错误！\n"));

                nRetCode= -;

                break;

            }

            //cImage数据存到pbImage，后再转换为源灰度图pbSrc

            iPitch=cImage.GetPitch();

            iBytePerPixel=(cImage.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage.Destroy();

            //TO DO:执行操作

            hResult=Imagejoint(pbSrc,iWidth,iHeight,dbZoom,pbFinal);

            if(hResult!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像处理错误！\n"));

                nRetCode= -;

                break;

            }

            //处理后保存图像

            iWidth=;

            iHeight=;

            cImage.Create(iWidth,-iHeight,);

            iPitch=cImage.GetPitch();

            for(y=;y<iHeight;y++)

            {

                pdwImage=(PDWORD)(PBYTE(cImage.GetBits())+iPitch*y);

                for(x=;x<iWidth;x++)

                {

                    pdwImage[x]=pbFinal[y*iWidth+x]*0x10101;

                }

            }

            //待保存图像文件名

            CString name1="target";

            CString name2;

            name2.Format(_T("%.2lf"),dbZoom);

            CString csTagName;

            csTagName=name1+name2;

            //CString csTagName="target";

            csTagName.Trim();

            csTagName.MakeUpper();

            if(csTagName.Right()!=_T(".BMP") ) csTagName.Append(_T(".BMP"));

            hResult=cImage.Save(csTagName);

            if(hResult!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像结果保存错误！\n"));

                nRetCode= -;

                break;

            }

            _tprintf(_T("图像处理成功！\n"));

            nRetCode= ERROR_SUCCESS;

            break;

            //dbZoom=dbZoom+0.5;

            //}

        }while();

        if(pbSrc) free(pbSrc);

        if(pbFinal) free(pbFinal);

    }

    else

    {

        _tprintf(_T("Fatal Error: GetModuleHandle failed\n"));

        nRetCode = ;

    }

    getchar();

    return nRetCode;

}

HRESULT Imagejoint(PBYTE pbSrc,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal)

{

    double phase[]={};//16相位，包含端点存在33个距离

    for (int i=;i<;i++)

    {

        double i2=1.0*i;

        phase[i]=fs(i2/);

    }

    //旋转中心为图像中心

    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=;y<;y++)

    {

        for (int x=;x<;x++)

        {

            srcx=(double)(newWidth/-+x)/dbZoom;

            srcy=(double)(newHeight/-+y)/dbZoom ;

            xOr = floor(srcx);

            yOr = floor(srcy);

            u=srcx-xOr;

            v=srcy-yOr;

            int phasex=floor(*u+0.5);//16相位

            int phasey=floor(*v+0.5);

            double A1,B1,C1,D1,A2,B2,C2,D2;

            A1=phase[+phasex];

            B1=phase[phasex];

            C1=phase[-phasex];

            D1=phase[-phasex];

            A2=phase[+phasey];

            B2=phase[phasey];

            C2=phase[-phasey];

            D2=phase[-phasey];

            if( !(srcx>= && srcx<=iWidth && srcy>= && srcy<=iHeight))//越界部分需拼接为外圈大视场图像，远景

                {

                    srcx=(double)(newWidth/-+x)/(*dbZoom);

                    srcy=(double)(newHeight/-+y)/(*dbZoom) ;

                    xOr = floor(srcx);

                    yOr = floor(srcy);

                    u=srcx-xOr;

                    v=srcy-yOr;

                    int phasex=floor(*u+0.5);//16相位

                    int phasey=floor(*v+0.5);

                    double A1,B1,C1,D1,A2,B2,C2,D2;

                    A1=phase[+phasex];

                    B1=phase[phasex];

                    C1=phase[-phasex];

                    D1=phase[-phasex];

                    A2=phase[+phasey];

                    B2=phase[phasey];

                    C2=phase[-phasey];

                    D2=phase[-phasey];

                    double middle=

                    pbSrc[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

                    //pbFinal[y*1280+x]=255;

                }

            else

            {

                double middle=

                    pbSrc[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

            }

        }

    }

    return ERROR_SUCCESS;

}

3、初步改为利用双图拼接而来，但是其中一些参数不对，拼接位置不对。注意：代码起始阶段的图像长宽和内存分配，默认是两图相同的。

// Imagejoint.cpp : 定义控制台应用程序的入口点。

//

#include "stdafx.h"

#include "Imagejoint.h"

#include <afxwin.h>

#ifdef _DEBUG

#define new DEBUG_NEW

#endif

#include <atlimage.h>//CImage类

#include <locale.h>

#include "math.h"

using namespace std;

//双三次插值系数

double fs(double w)

{

        double a=-0.5;

        double fs;

        if (abs(w)<=)

            fs=(a+)*pow(abs(w),)-(a+)*pow(abs(w),)+;

        else if (abs(w)>&&abs(w)<=)

            fs=a*pow(abs(w),)-*a*pow(abs(w),)+*a*abs(w)-*a;

        else

            fs=;

        return fs;

}

HRESULT Imagejoint(PBYTE pbSrc1,PBYTE pbSrc2,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal);

int _tmain(int argc, TCHAR* argv[], TCHAR* envp[])

{

    int nRetCode = ;//表示整数类型的函数返回码。n表示整数类型,Ret是Return的缩写,表示返回值,Code表示代码。

    setlocale(LC_ALL,"chs");

    HMODULE hModule = ::GetModuleHandle(NULL);

    if (hModule != NULL)

    {

        // initialize MFC and print and error on failure

        CImage cImage_far;

        CImage cImage_near;

        HRESULT hResult1,hResult2;

        //初始化一些变量

        int        iWidth,iHeight,iBytePerPixel,iPitch;

        int        x,y;

        PBYTE    pbSrc1=NULL,pbSrc2=NULL,pbFinal=NULL;//源图、目标图

        PBYTE    pbImage=NULL;//load图像后存在这

        PDWORD    pdwImage=NULL;//用于保存图像

        double    dbZoom=1.25;

        CString str = "far-frame1.bmp";

        CString str2 = "near-frame1.bmp";

        LPCTSTR filename1 = (LPCTSTR)str;

        LPCTSTR filename2 = (LPCTSTR)str2;

        do{

            if (!AfxWinInit(hModule, NULL, ::GetCommandLine(), ))

            {

                // TODO: change error code to suit your needs

                _tprintf(_T("Fatal Error: MFC initialization failed\n"));

                nRetCode = ;

                break;

            }

            //Load 图像到cImage对象中

            hResult1=cImage_far.Load(filename1);

            hResult2=cImage_near.Load(filename2);

            if(hResult1!=ERROR_SUCCESS||hResult2!=ERROR_SUCCESS)

            {

                _tprintf(_T("源图像文件名错误！\n"));

                nRetCode= -;

                break;

            }

            iWidth=cImage_far.GetWidth();

            iHeight=cImage_far.GetHeight();

            //分配源图内存

            pbSrc1 = (PBYTE)malloc(iWidth*iHeight);

            pbSrc2 = (PBYTE)malloc(iWidth*iHeight);

            pbFinal = (PBYTE)malloc(*);

        //while (dbZoom<3);

        //{

            //分配目标图内存

            if(pbSrc1==NULL||pbSrc2==NULL || pbFinal==NULL )

            {

                _tprintf(_T("内存申请错误！\n"));

                nRetCode= -;

                break;

            }

            //cImage数据存到pbImage，后再转换为源灰度图pbSrc

            iPitch=cImage_far.GetPitch();

            iBytePerPixel=(cImage_far.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage_far.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc1[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage_far.Destroy();

            iPitch=cImage_near.GetPitch();

            iBytePerPixel=(cImage_near.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage_near.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc2[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage_near.Destroy();

            //TO DO:执行操作

            hResult1=Imagejoint(pbSrc1,pbSrc2,iWidth,iHeight,dbZoom,pbFinal);

            if(hResult1!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像处理错误！\n"));

                nRetCode= -;

                break;

            }

            //处理后保存图像

            iWidth=;

            iHeight=;

            cImage_far.Create(iWidth,-iHeight,);

            iPitch=cImage_far.GetPitch();

            for(y=;y<iHeight;y++)

            {

                pdwImage=(PDWORD)(PBYTE(cImage_far.GetBits())+iPitch*y);

                for(x=;x<iWidth;x++)

                {

                    pdwImage[x]=pbFinal[y*iWidth+x]*0x10101;

                }

            }

            //待保存图像文件名

            CString name1="target";

            CString name2;

            name2.Format(_T("%.2lf"),dbZoom);

            CString csTagName;

            csTagName=name1+name2;

            //CString csTagName="target";

            csTagName.Trim();

            csTagName.MakeUpper();

            if(csTagName.Right()!=_T(".BMP") ) csTagName.Append(_T(".BMP"));

            hResult1=cImage_far.Save(csTagName);

            if(hResult1!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像结果保存错误！\n"));

                nRetCode= -;

                break;

            }

            _tprintf(_T("图像处理成功！\n"));

            nRetCode= ERROR_SUCCESS;

            break;

            //dbZoom=dbZoom+0.5;

            //}

        }while();

        if(pbSrc1) free(pbSrc1);

        if(pbSrc2) free(pbSrc2);

        if(pbFinal) free(pbFinal);

    }

    else

    {

        _tprintf(_T("Fatal Error: GetModuleHandle failed\n"));

        nRetCode = ;

    }

    getchar();

    return nRetCode;

}

HRESULT Imagejoint(PBYTE pbSrc1,PBYTE pbSrc2,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal)

{

    double phase[]={};//16相位，包含端点存在33个距离

    for (int i=;i<;i++)

    {

        double i2=1.0*i;

        phase[i]=fs(i2/);

    }

    //旋转中心为图像中心

    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=;y<;y++)

    {

        for (int x=;x<;x++)

        {

            srcx=(double)(newWidth/-+x)/dbZoom;

            srcy=(double)(newHeight/-+y)/dbZoom ;

            xOr = floor(srcx);

            yOr = floor(srcy);

            u=srcx-xOr;

            v=srcy-yOr;

            int phasex=floor(*u+0.5);//16相位

            int phasey=floor(*v+0.5);

            double A1,B1,C1,D1,A2,B2,C2,D2;

            A1=phase[+phasex];

            B1=phase[phasex];

            C1=phase[-phasex];

            D1=phase[-phasex];

            A2=phase[+phasey];

            B2=phase[phasey];

            C2=phase[-phasey];

            D2=phase[-phasey];

            if( !(srcx>= && srcx<=iWidth && srcy>= && srcy<=iHeight))//越界部分需拼接为外圈大视场图像，远景

                {

                    srcx=(double)(newWidth/-+x)/(*dbZoom);

                    srcy=(double)(newHeight/-+y)/(*dbZoom) ;

                    xOr = floor(srcx);

                    yOr = floor(srcy);

                    u=srcx-xOr;

                    v=srcy-yOr;

                    int phasex=floor(*u+0.5);//16相位

                    int phasey=floor(*v+0.5);

                    double A1,B1,C1,D1,A2,B2,C2,D2;

                    A1=phase[+phasex];

                    B1=phase[phasex];

                    C1=phase[-phasex];

                    D1=phase[-phasex];

                    A2=phase[+phasey];

                    B2=phase[phasey];

                    C2=phase[-phasey];

                    D2=phase[-phasey];

                    double middle=

                    pbSrc1[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

                    //pbFinal[y*1280+x]=255;

                }

            else

            {

                double middle=

                    pbSrc2[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

            }

        }

    }

    return ERROR_SUCCESS;

}

4、对拼接位置修正，x+6，y-3。右6，上3。

// Imagejoint.cpp : 定义控制台应用程序的入口点。

//

#include "stdafx.h"

#include "Imagejoint.h"

#include <afxwin.h>

#ifdef _DEBUG

#define new DEBUG_NEW

#endif

#include <atlimage.h>//CImage类

#include <locale.h>

#include "math.h"

using namespace std;

//双三次插值系数

double fs(double w)

{

        double a=-0.5;

        double fs;

        if (abs(w)<=)

            fs=(a+)*pow(abs(w),)-(a+)*pow(abs(w),)+;

        else if (abs(w)>&&abs(w)<=)

            fs=a*pow(abs(w),)-*a*pow(abs(w),)+*a*abs(w)-*a;

        else

            fs=;

        return fs;

}

HRESULT Imagejoint(PBYTE pbSrc1,PBYTE pbSrc2,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal);

int _tmain(int argc, TCHAR* argv[], TCHAR* envp[])

{

    int nRetCode = ;//表示整数类型的函数返回码。n表示整数类型,Ret是Return的缩写,表示返回值,Code表示代码。

    setlocale(LC_ALL,"chs");

    HMODULE hModule = ::GetModuleHandle(NULL);

    if (hModule != NULL)

    {

        // initialize MFC and print and error on failure

        CImage cImage_far;

        CImage cImage_near;

        HRESULT hResult1,hResult2;

        //初始化一些变量

        int        iWidth,iHeight,iBytePerPixel,iPitch;

        int        x,y;

        PBYTE    pbSrc1=NULL,pbSrc2=NULL,pbFinal=NULL;//源图、目标图

        PBYTE    pbImage=NULL;//load图像后存在这

        PDWORD    pdwImage=NULL;//用于保存图像

        double    dbZoom=1.5;

        CString str = "far-frame1.bmp";

        CString str2 = "near-frame1.bmp";

        LPCTSTR filename1 = (LPCTSTR)str;

        LPCTSTR filename2 = (LPCTSTR)str2;

        do{

            if (!AfxWinInit(hModule, NULL, ::GetCommandLine(), ))

            {

                // TODO: change error code to suit your needs

                _tprintf(_T("Fatal Error: MFC initialization failed\n"));

                nRetCode = ;

                break;

            }

            //Load 图像到cImage对象中

            hResult1=cImage_far.Load(filename1);

            hResult2=cImage_near.Load(filename2);

            if(hResult1!=ERROR_SUCCESS||hResult2!=ERROR_SUCCESS)

            {

                _tprintf(_T("源图像文件名错误！\n"));

                nRetCode= -;

                break;

            }

            iWidth=cImage_far.GetWidth();

            iHeight=cImage_far.GetHeight();

            //分配源图内存

            pbSrc1 = (PBYTE)malloc(iWidth*iHeight);

            pbSrc2 = (PBYTE)malloc(iWidth*iHeight);

            pbFinal = (PBYTE)malloc(*);

        //while (dbZoom<3);

        //{

            //分配目标图内存

            if(pbSrc1==NULL||pbSrc2==NULL || pbFinal==NULL )

            {

                _tprintf(_T("内存申请错误！\n"));

                nRetCode= -;

                break;

            }

            //cImage数据存到pbImage，后再转换为源灰度图pbSrc

            iPitch=cImage_far.GetPitch();

            iBytePerPixel=(cImage_far.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage_far.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc1[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage_far.Destroy();

            iPitch=cImage_near.GetPitch();

            iBytePerPixel=(cImage_near.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage_near.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc2[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage_near.Destroy();

            //TO DO:执行操作

            hResult1=Imagejoint(pbSrc1,pbSrc2,iWidth,iHeight,dbZoom,pbFinal);

            if(hResult1!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像处理错误！\n"));

                nRetCode= -;

                break;

            }

            //处理后保存图像

            iWidth=;

            iHeight=;

            cImage_far.Create(iWidth,-iHeight,);

            iPitch=cImage_far.GetPitch();

            for(y=;y<iHeight;y++)

            {

                pdwImage=(PDWORD)(PBYTE(cImage_far.GetBits())+iPitch*y);

                for(x=;x<iWidth;x++)

                {

                    pdwImage[x]=pbFinal[y*iWidth+x]*0x10101;

                }

            }

            //待保存图像文件名

            CString name1="target";

            CString name2;

            name2.Format(_T("%.2lf"),dbZoom);

            CString csTagName;

            csTagName=name1+name2;

            //CString csTagName="target";

            csTagName.Trim();

            csTagName.MakeUpper();

            if(csTagName.Right()!=_T(".BMP") ) csTagName.Append(_T(".BMP"));

            hResult1=cImage_far.Save(csTagName);

            if(hResult1!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像结果保存错误！\n"));

                nRetCode= -;

                break;

            }

            _tprintf(_T("图像处理成功！\n"));

            nRetCode= ERROR_SUCCESS;

            break;

            //dbZoom=dbZoom+0.5;

            //}

        }while();

        if(pbSrc1) free(pbSrc1);

        if(pbSrc2) free(pbSrc2);

        if(pbFinal) free(pbFinal);

    }

    else

    {

        _tprintf(_T("Fatal Error: GetModuleHandle failed\n"));

        nRetCode = ;

    }

    getchar();

    return nRetCode;

}

HRESULT Imagejoint(PBYTE pbSrc1,PBYTE pbSrc2,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal)

{

    double phase[]={};//16相位，包含端点存在33个距离

    for (int i=;i<;i++)

    {

        double i2=1.0*i;

        phase[i]=fs(i2/);

    }

    //旋转中心为图像中心

    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=;y<;y++)

    {

        for (int x=;x<;x++)

        {

            srcx=(double)(newWidth/-+x)/dbZoom;

            srcy=(double)(newHeight/-+y)/dbZoom ;

            xOr = floor(srcx);

            yOr = floor(srcy);

            u=srcx-xOr;

            v=srcy-yOr;

            int phasex=floor(*u+0.5);//16相位

            int phasey=floor(*v+0.5);

            double A1,B1,C1,D1,A2,B2,C2,D2;

            A1=phase[+phasex];

            B1=phase[phasex];

            C1=phase[-phasex];

            D1=phase[-phasex];

            A2=phase[+phasey];

            B2=phase[phasey];

            C2=phase[-phasey];

            D2=phase[-phasey];

            //Matlab中测试获取定值(34:541 -Y,41:728 -X)，需要微调。

            if( !(srcx>= && srcx<= && srcy>= && srcy<=))//越界部分需拼接为外圈大视场图像，远景

                {

                    srcx=(double)(newWidth-+x)/(*dbZoom)+;//补偏差

                    srcy=(double)(newHeight-+y)/(*dbZoom)-;

                    xOr = floor(srcx);

                    yOr = floor(srcy);

                    u=srcx-xOr;

                    v=srcy-yOr;

                    int phasex=floor(*u+0.5);//16相位

                    int phasey=floor(*v+0.5);

                    double A1,B1,C1,D1,A2,B2,C2,D2;

                    A1=phase[+phasex];

                    B1=phase[phasex];

                    C1=phase[-phasex];

                    D1=phase[-phasex];

                    A2=phase[+phasey];

                    B2=phase[phasey];

                    C2=phase[-phasey];

                    D2=phase[-phasey];

                    double middle=

                    pbSrc1[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

                    //pbFinal[y*1280+x]=255;

                }

            else

            {

                double middle=

                    pbSrc2[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

            }

        }

    }

    return ERROR_SUCCESS;

}

5、加入dbZoom自循环。

for(dbZoom=1.25;dbZoom<=2.25;dbZoom+=0.25)

但要注意！由于之后有:cImage_far.Destroy();和cImage_near.Destroy();

因此循环内应包含！CImage cImage_far;CImage cImage_near;

// Imagejoint.cpp : 定义控制台应用程序的入口点。

//

#include "stdafx.h"

#include "Imagejoint.h"

#include <afxwin.h>

#ifdef _DEBUG

#define new DEBUG_NEW

#endif

#include <atlimage.h>//CImage类

#include <locale.h>

#include "math.h"

using namespace std;

//双三次插值系数

double fs(double w)

{

        double a=-0.5;

        double fs;

        if (abs(w)<=)

            fs=(a+)*pow(abs(w),)-(a+)*pow(abs(w),)+;

        else if (abs(w)>&&abs(w)<=)

            fs=a*pow(abs(w),)-*a*pow(abs(w),)+*a*abs(w)-*a;

        else

            fs=;

        return fs;

}

HRESULT Imagejoint(PBYTE pbSrc1,PBYTE pbSrc2,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal);

int _tmain(int argc, TCHAR* argv[], TCHAR* envp[])

{

    int nRetCode = ;//表示整数类型的函数返回码。n表示整数类型,Ret是Return的缩写,表示返回值,Code表示代码。

    setlocale(LC_ALL,"chs");

    HMODULE hModule = ::GetModuleHandle(NULL);

    if (hModule != NULL)

    {

        HRESULT hResult1,hResult2;

        //初始化一些变量

        int        iWidth,iHeight,iBytePerPixel,iPitch;

        int        x,y;

        PBYTE    pbSrc1=NULL,pbSrc2=NULL,pbFinal=NULL;//源图、目标图

        PBYTE    pbImage=NULL;//load图像后存在这

        PDWORD    pdwImage=NULL;//用于保存图像

        double    dbZoom=1.25;

        for(dbZoom=1.25;dbZoom<=2.25;dbZoom+=0.25)

        {

            CImage cImage_far;

            CImage cImage_near;

            CString str = "far-frame1.bmp";

            CString str2 = "near-frame1.bmp";

            LPCTSTR filename1 = (LPCTSTR)str;

            LPCTSTR filename2 = (LPCTSTR)str2;

            if (!AfxWinInit(hModule, NULL, ::GetCommandLine(), ))

            {

                // TODO: change error code to suit your needs

                _tprintf(_T("Fatal Error: MFC initialization failed\n"));

                nRetCode = ;

            }

            //Load 图像到cImage对象中

            hResult1=cImage_far.Load(filename1);

            hResult2=cImage_near.Load(filename2);

            if(hResult1!=ERROR_SUCCESS||hResult2!=ERROR_SUCCESS)

            {

                _tprintf(_T("源图像文件名错误！\n"));

                nRetCode= -;

            }

            iWidth=cImage_far.GetWidth();

            iHeight=cImage_far.GetHeight();

            //分配源图内存

            pbSrc1 = (PBYTE)malloc(iWidth*iHeight);

            pbSrc2 = (PBYTE)malloc(iWidth*iHeight);

            //分配目标图内存

            pbFinal = (PBYTE)malloc(*);

            if(pbSrc1==NULL||pbSrc2==NULL || pbFinal==NULL )

            {

                _tprintf(_T("内存申请错误！\n"));

                nRetCode= -;

            }

            //cImage数据存到pbImage，后再转换为源灰度图pbSrc

            iPitch=cImage_far.GetPitch();

            iBytePerPixel=(cImage_far.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage_far.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc1[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage_far.Destroy();

            iPitch=cImage_near.GetPitch();

            iBytePerPixel=(cImage_near.GetBPP()+)/;

            if(iBytePerPixel==)

            {

                for(y=;y<iHeight;y++)

                {   //load的图像数据放到pbImage

                    pbImage=(PBYTE)(PBYTE(cImage_near.GetBits())+iPitch*y);//得到的是图像初始像素地址

                    for(x=;x<iWidth;x++)

                    {

                        //pbImage转换为灰度图pbSrc

                        pbSrc2[y*iWidth+x]=(pbImage[*x]*0.15+pbImage[*x+]*0.55+pbImage[*x+]*0.3);

                    }

                }

            }

            cImage_near.Destroy();

            //执行操作

            hResult1=Imagejoint(pbSrc1,pbSrc2,iWidth,iHeight,dbZoom,pbFinal);

            if(hResult1!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像处理错误！\n"));

                nRetCode= -;

            }

            //处理后保存图像

            iWidth=;

            iHeight=;

            cImage_far.Create(iWidth,-iHeight,);

            iPitch=cImage_far.GetPitch();

            for(y=;y<iHeight;y++)

            {

                pdwImage=(PDWORD)(PBYTE(cImage_far.GetBits())+iPitch*y);

                for(x=;x<iWidth;x++)

                {

                    pdwImage[x]=pbFinal[y*iWidth+x]*0x10101;

                }

            }

            //可预存待保存图像文件名

            CString name1="target";

            CString name2;

            name2.Format(_T("%.2lf"),dbZoom);

            CString csTagName;

            csTagName=name1+name2;

            csTagName.Trim();

            csTagName.MakeUpper();

            if(csTagName.Right()!=_T(".BMP") ) csTagName.Append(_T(".BMP"));

            hResult1=cImage_far.Save(csTagName);

            if(hResult1!=ERROR_SUCCESS)

            {

                _tprintf(_T("图像结果保存错误！\n"));

                nRetCode= -;

            }

            _tprintf(_T("图像处理成功！\n"));

            nRetCode= ERROR_SUCCESS;

            if(pbSrc1) free(pbSrc1);

            if(pbSrc2) free(pbSrc2);

            if(pbFinal) free(pbFinal);

        }//对应+=0.25倍数的for循环。

    }

    else

    {

        _tprintf(_T("Fatal Error: GetModuleHandle failed\n"));

        nRetCode = ;

    }

    getchar();

    return nRetCode;

}

HRESULT Imagejoint(PBYTE pbSrc1,PBYTE pbSrc2,int iWidth,int iHeight,double dbZoom,PBYTE pbFinal)

{

    double phase[]={};//16相位，包含端点存在33个距离

    for (int i=;i<;i++)

    {

        double i2=1.0*i;

        phase[i]=fs(i2/);

    }

    //旋转中心为图像中心

    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=;y<;y++)

    {

        for (int x=;x<;x++)

        {

            srcx=(double)(newWidth/-+x)/dbZoom;

            srcy=(double)(newHeight/-+y)/dbZoom ;

            xOr = floor(srcx);

            yOr = floor(srcy);

            u=srcx-xOr;

            v=srcy-yOr;

            int phasex=floor(*u+0.5);//16相位

            int phasey=floor(*v+0.5);

            double A1,B1,C1,D1,A2,B2,C2,D2;

            A1=phase[+phasex];

            B1=phase[phasex];

            C1=phase[-phasex];

            D1=phase[-phasex];

            A2=phase[+phasey];

            B2=phase[phasey];

            C2=phase[-phasey];

            D2=phase[-phasey];

            //Matlab中测试获取定值(34:541 -Y,41:728 -X)，需要微调。

            if( !(srcx>= && srcx<= && srcy>= && srcy<=))//越界部分需拼接为外圈大视场图像，远景

                {

                    srcx=(double)(newWidth-+x)/(*dbZoom)+;//补偏差

                    srcy=(double)(newHeight-+y)/(*dbZoom)-;

                    xOr = floor(srcx);

                    yOr = floor(srcy);

                    u=srcx-xOr;

                    v=srcy-yOr;

                    int phasex=floor(*u+0.5);//16相位

                    int phasey=floor(*v+0.5);

                    double A1,B1,C1,D1,A2,B2,C2,D2;

                    A1=phase[+phasex];

                    B1=phase[phasex];

                    C1=phase[-phasex];

                    D1=phase[-phasex];

                    A2=phase[+phasey];

                    B2=phase[phasey];

                    C2=phase[-phasey];

                    D2=phase[-phasey];

                    double middle=

                    pbSrc1[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc1[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc1[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc1[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

                    //pbFinal[y*1280+x]=255;

                }

            else

            {

                double middle=

                    pbSrc2[(yOr-)*iWidth+(xOr-)]*A1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr-)]*A1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr-)]*A1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr-)]*A1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr)]*B1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr)]*B1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr)]*B1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr)]*B1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr+)]*C1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr+)]*C1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*C1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*C1*D2+

                    pbSrc2[(yOr-)*iWidth+(xOr+)]*D1*A2+

                    pbSrc2[(yOr)*iWidth+(xOr+)]*D1*B2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*D1*C2+

                    pbSrc2[(yOr+)*iWidth+(xOr+)]*D1*D2;

                if(middle<=&&middle>=)

                    pbFinal[y*+x]=middle;

                else if(middle>)

                    pbFinal[y*+x]=;

                else

                    pbFinal[y*+x]=;

            }

        }

    }

    return ERROR_SUCCESS;

}

至此，基本完成效果仿真。

但目前只是针对图像的。下面考虑如何读取.avi格式的视频文件处理。

20190710记录：去掉中转图，直接以1280*1024进行反坐标计算，填补pbFinal。

20190710记录：去掉中转图，直接以1280*1024进行反坐标计算，填补pbFinal。的相关教程结束。

相关推荐

2022-09-03：n块石头放置在二维平面中的一些整数坐标点上每个坐标点上最多只能有一块石头如果一块石头的同行或者同列上有其他石头存在，那么就可以移除这块石头。给你一个长度为 n 的数组

Android自定义View研究--View中的原点坐标和XML中布局自定义View时View触摸原点问题

利用百度地图API，获取经纬度坐标

Android 性能优化（2）性能工具之「Hierarchy Viewer 」Optimizing Your UI：分析哪个view有性能问题，查看屏幕上某像素点的坐标，颜色等

php计算两个坐标（经度，纬度）之间距离的方法

在matlab中进行地理坐标和像素坐标的相互转换

基于均值坐标(Mean-Value Coordinates)的图像融合算法的优化实现

java 根据经纬度坐标计算两点的距离算法

20190710记录：去掉中转图，直接以1280*1024进行反坐标计算，填补pbFinal。

20190710记录：去掉中转图，直接以1280*1024进行反坐标计算，填补pbFinal。的相关教程结束。

相关推荐

2022-09-03：n块石头放置在二维平面中的一些整数坐标点上 每个坐标点上最多只能有一块石头 如果一块石头的 同行或者同列 上有其他石头存在，那么就可以移除这块石头。 给你一个长度为 n 的数组

Android自定义View研究--View中的原点坐标和XML中布局自定义View时View触摸原点问题

利用百度地图API，获取经纬度坐标

Android 性能优化（2）性能工具之「Hierarchy Viewer 」Optimizing Your UI：分析哪个view有性能问题，查看屏幕上某像素点的坐标，颜色等

php计算两个坐标（经度，纬度）之间距离的方法

在matlab中进行地理坐标和像素坐标的相互转换

基于均值坐标(Mean-Value Coordinates)的图像融合算法的优化实现

java 根据经纬度坐标计算两点的距离算法

2022-09-03：n块石头放置在二维平面中的一些整数坐标点上每个坐标点上最多只能有一块石头如果一块石头的同行或者同列上有其他石头存在，那么就可以移除这块石头。给你一个长度为 n 的数组