基于相关系数影像匹配实习报告

基于相关系数影像匹配实习报告一、实习内容与目的通过上机编程实现基于相关系数算法的影像匹配，从而更加熟练的掌握书本上关于影像匹配及相关系数算法的理论知识，并熟悉上机编程的操作。二、实习原理影像匹配实质上是在两幅或是多幅影像之间识别同名点，其中基于相关系数的算法是实现影像匹配的基本算法之一。相关系数是标准化的协方差函数，协方差函数除以两信号的方差即得相关系数。计算相关系数的公式如下：g(x,y)与g’(x’,y’)的相关系数为：),(),(),(qpCCqpCqpgggg其中：DyxggdxdyyxgEyxgC),(2)]},([),({DyxggdxdyqypxgEqypxgqpC),(2)]},([),({),(若(p0,q0)(p,q)，(p≠p0,q≠q0)，则p0,q0为搜索区影像相对于目标区影像的位移参数。由离散灰度数据对相关系数的估计为minjminjcrcjrijiminjcjrijiggggggggrc11112,,2,11,,)()())((),(minjcjrircgnmg11,,1，minjjignmg11,1考虑到计算工作量，相关系数的实用公式为minjminjminjcjricjriminjjijiminjminjcjriminjjicjrijignmggnmgggnmggrc1111211,,2211,,21111,11,,,])(1][)(1[)()(1)(),(三、实习步骤及相关代码在第一次点特征点实习的基础上我们组进行了第二次基于相关系数的影像匹配实习。所以程序实现的部分操作直接利用前次实习的中的程序。（1）读入左右影像FILE*pSrcFileLeft=NULL;pSrcFileLeft=fopen(srcFileLeft,rb);if(pSrcFileLeft==NULL){printf(原始左片影像打开失败\n);returnFAILURE;}FILE*pSrcFileRight=NULL;pSrcFileRight=fopen(srcFileRight,rb);if(pSrcFileRight==NULL){printf(原始右片影像打开失败\n);returnFAILURE;}（2）确定目标窗口大小通过老师课堂上的讲授和我们在编程过程中的不断尝试，我们默认使用的窗口的尺寸为11*11。（3）确定目标点的位置由前次的实习中得到左相片中的特征点作为这次实习的目标点。（4）预测右影像搜索范围本程序中也在左相片中计算出了特征点，然后将左右两相片中的特征点进行计算相关系数。（5）逐窗口计算相关系数并保存相应的储存单元将左相片中的目标点与右相片中的特征点进行相关系数的计算。（6）比较相关系数，取最大相关系数作为匹配点对于分别比较上一步骤中计算出的相关系数，选取出最大的一个相关系数作为匹配点。程序操作如下：在程序主界面中点击Extend中选取ImageMatch选项，然后进入参数设置界面。设置读取和保存影像路径和窗口大小等参数后，点击Progress进行运算。四、实习结果及分析作为数据源的左右影像如些图所示：左相片右相片经过计算后得到的目标点左右图像分别如下所示：左相片右相片显示匹配点号：通过老师得到的图片可以看出在左相片中得到的五个目标点中有四个在右相片中找到了同名点，由此可以看出本程序基本满足实习要求，并得到理想的结果。不足的地方就是选取的特征点的数量不够多，这样就无法判断本程序的“误判”的情况如何。五、实习体会在学习到相关系数的时候，一看到计算相关系数的公式的时候我就有种被折服的感觉，那个公式实在是太复杂，它使我感觉我怎么也不能记住这个复杂的公式。不过经过这次实习后，我觉得这个公式又不是当初认为的那样无法记忆。我想这次实习我最大的收获就是对于之前在学习过程中遇到不懂知识以及难以记忆的知识都得到了加强，方便了我对他们的理解与记忆。其次，在编程过程中更加锻炼自己动手编程的能力，使得我编程能力较之前有着一定的提高。每次编程实习都是一个积累经验的过程，然后通过一次次的实习后，我最终的能力就能够得到很好的体现。附件：voidMoravec(BYTE*pSrcBits,vectorFEATUREPOINT*FeaturePoint,intm_threshold,intm_window1,intm_window2,intwidth,intheight,int*Count1){。。。。。。。。此部分程序与点特征程序中代码一样，因此不再重复粘贴}vectorFEATUREPOINTRemoveReplicative(vectorFEATUREPOINTv){vectorFEATUREPOINTret;vectorFEATUREPOINT::iteratoriter=v.begin();ret.clear();ret.push_back(*iter);vectorFEATUREPOINT::iteratoriter2;BOOLb(FALSE);for(iter=v.begin();iter!=v.end();++iter){b=false;for(iter2=ret.begin();iter2!=ret.end();++iter2){if(iter-x==iter2-x&&iter-y==iter2-y){b=TRUE;//存在break;}}if(b==FALSE){ret.push_back(*iter);}}returnret;}voidSaveBand(intwidth,intheight,intbyteCount,intbiBitCount,LPBYTEpBits,CStringSavePath){BITMAPFILEHEADERbmfh;BITMAPINFOHEADERbmih;bmfh.bfType=0x4d42;//0x42=B0x4d=Mbmfh.bfReserved1=0;bmfh.bfReserved2=0;if(biBitCount==8){bmfh.bfOffBits=sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+256*4;bmfh.bfSize=sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+256*4+byteCount;}elseif(biBitCount==24){bmfh.bfOffBits=sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER);bmfh.bfSize=sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+byteCount;}else{MessageBox(NULL,8bitor24bitaccepted,错误,MB_ICONERROR);return;}bmih.biBitCount=biBitCount;bmih.biWidth=width;bmih.biHeight=height;bmih.biSize=40;bmih.biPlanes=1;bmih.biCompression=BI_RGB;bmih.biSizeImage=bmfh.bfSize;bmih.biClrImportant=0;bmih.biClrUsed=0;bmih.biXPelsPerMeter=0;bmih.biYPelsPerMeter=0;CFilef;//保存位图if(f.Open(SavePath,CFile::modeCreate|CFile::modeWrite)){f.Write(&bmfh,sizeof(bmfh));f.Write(&bmih,sizeof(bmih));if(biBitCount==8){RGBQUADrgb[256];inti(0);for(;i256;++i){rgb[i].rgbBlue=i;rgb[i].rgbGreen=i;rgb[i].rgbRed=i;rgb[i].rgbReserved=0;}//修改颜色表255对应的颜色rgb[255].rgbBlue=0;rgb[255].rgbGreen=0;rgb[255].rgbRed=255;f.Write(rgb,sizeof(RGBQUAD)*256);}f.Write(pBits,byteCount);f.Close();}else{MessageBox(NULL,保存图像失败,错误,MB_ICONERROR);return;}}voidDrawCross(BYTE*pBits,vectorFEATUREPOINTFeaturePoint,int_width,intbyteCount){vectorFEATUREPOINT::iteratoriter;inti(0),j(0);for(i=0;ibyteCount;++i){if(pBits[i]==255)pBits[i]=254;}for(iter=FeaturePoint.begin();iter!=FeaturePoint.end();++iter){i=iter-y;j=iter-x;pBits[i*_width+j]=255;pBits[i*_width+j+1]=255;pBits[i*_width+j-1]=255;pBits[(i+1)*_width+j]=255;pBits[(i-1)*_width+j]=255;pBits[i*_width+j+2]=255;pBits[i*_width+j-2]=255;pBits[(i+2)*_width+j]=255;pBits[(i-2)*_width+j]=255;}}voidDrawCross(BYTE*pBits1,BYTE*pBits2,vectorMATCHEDPOINTSMatchedPoints,int_width1,int_width2,intbyteCount1,intbyteCount2){vectorMATCHEDPOINTS::iteratoriter;inti(0),j(0);for(i=0;ibyteCount1;++i){if(pBits1[i]==255)pBits1[i]=254;}for(j=0;jbyteCount2;++j){if(pBits2[j]==255)pBits2[j]=254;}for(iter=MatchedPoints.begin();iter!=MatchedPoints.end();++iter){i=iter-y1;j=iter-x1;pBits1[i*_width1+j]=255;pBits1[i*_width1+j+1]=255;pBits1[i*_width1+j-1]=255;pBits1[(i+1)*_width1+j]=255;pBits1[(i-1)*_width1+j]=255;pBits1[i*_width1+j+2]=255;pBits1[i*_width1+j-2]=255;pBits1[(i+2)*_width1+j]=255;pBits1[(i-2)*_width1+j]=255;i=iter-y2;j=iter-x2;pBits2[i*_width2+j]=255;//pBits2[i*_width2+j+1]=255;pBits2[i*_width2+j-1]=255;pBits2[(i+1)*_width2+j]=255;pBits2[(i-1)*_width2+j]=255;pBits2[i*_width2+j+2]=255;pBits2[i*_width2+j-2]