基于相關(guān)系數(shù)影像匹配實(shí)習(xí)報(bào)告(共15頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上基于相關(guān)系數(shù)影像匹配實(shí)習(xí)報(bào)告1、 實(shí)習(xí)內(nèi)容與目的 通過(guò)上機(jī)編程實(shí)現(xiàn)基于相關(guān)系數(shù)算法的影像匹配，從而更加熟練的掌握書(shū)本上關(guān)于影像匹配及相關(guān)系數(shù)算法的理論知識(shí)，并熟悉上機(jī)編程的操作。2、 實(shí)習(xí)原理影像匹配實(shí)質(zhì)上是在兩幅或是多幅影像之間識(shí)別同名點(diǎn)，其中基于相關(guān)系數(shù)的算法是實(shí)現(xiàn)影像匹配的基本算法之一。相關(guān)系數(shù)是標(biāo)準(zhǔn)化的協(xié)方差函數(shù)，協(xié)方差函數(shù)除以?xún)尚盘?hào)的方差即得相關(guān)系數(shù)。計(jì)算相關(guān)系數(shù)的公式如下：g(x,y)與g(x,y)的相關(guān)系數(shù)為：其中：若(p0, q0)>(p, q)，(pp0,qq0)，則p0,q0為搜索區(qū)影像相對(duì)于目標(biāo)區(qū)影像的位移參數(shù)。由離散灰度數(shù)據(jù)對(duì)相關(guān)系數(shù)的估

2、計(jì)為，考慮到計(jì)算工作量，相關(guān)系數(shù)的實(shí)用公式為3、 實(shí)習(xí)步驟及相關(guān)代碼 在第一次點(diǎn)特征點(diǎn)實(shí)習(xí)的基礎(chǔ)上我們組進(jìn)行了第二次基于相關(guān)系數(shù)的影像匹配實(shí)習(xí)。所以程序?qū)崿F(xiàn)的部分操作直接利用前次實(shí)習(xí)的中的程序。（1） 讀入左右影像 FILE* pSrcFileLeft = NULL;pSrcFileLeft = fopen(srcFileLeft,"rb");if(pSrcFileLeft = NULL)printf("原始左片影像打開(kāi)失敗n");return FAILURE;FILE* pSrcFileRight = NULL;pSrcFileRight = fope

3、n(srcFileRight,"rb");if(pSrcFileRight = NULL)printf("原始右片影像打開(kāi)失敗n");return FAILURE;（2） 確定目標(biāo)窗口大小 通過(guò)老師課堂上的講授和我們?cè)诰幊踢^(guò)程中的不斷嘗試，我們默認(rèn)使用的窗口的尺寸為11*11。（3） 確定目標(biāo)點(diǎn)的位置 由前次的實(shí)習(xí)中得到左相片中的特征點(diǎn)作為這次實(shí)習(xí)的目標(biāo)點(diǎn)。（4） 預(yù)測(cè)右影像搜索范圍 本程序中也在左相片中計(jì)算出了特征點(diǎn)，然后將左右兩相片中的特征點(diǎn)進(jìn)行計(jì)算相關(guān)系數(shù)。（5） 逐窗口計(jì)算相關(guān)系數(shù)并保存相應(yīng)的儲(chǔ)存單元 將左相片中的目標(biāo)點(diǎn)與右相片中的特征點(diǎn)進(jìn)行相關(guān)

4、系數(shù)的計(jì)算。（6） 比較相關(guān)系數(shù)，取最大相關(guān)系數(shù)作為匹配點(diǎn) 對(duì)于分別比較上一步驟中計(jì)算出的相關(guān)系數(shù)，選取出最大的一個(gè)相關(guān)系數(shù)作為匹配點(diǎn)。 程序操作如下：在程序主界面中點(diǎn)擊Extend中選取ImageMatch選項(xiàng)，然后進(jìn)入?yún)?shù)設(shè)置界面。設(shè)置讀取和保存影像路徑和窗口大小等參數(shù)后，點(diǎn)擊Progress進(jìn)行運(yùn)算。4、 實(shí)習(xí)結(jié)果及分析 作為數(shù)據(jù)源的左右影像如些圖所示：左相片 右相片 經(jīng)過(guò)計(jì)算后得到的目標(biāo)點(diǎn)左右圖像分別如下所示：左相片 右相片 顯示匹配點(diǎn)號(hào)： 通過(guò)老師得到的圖片可以看出在左相片中得到的五個(gè)目標(biāo)點(diǎn)中有四個(gè)在右相片中找到了同名點(diǎn)，由此可以看出本程序基本滿(mǎn)足實(shí)習(xí)要求，并得到理想的結(jié)果。不足的地

5、方就是選取的特征點(diǎn)的數(shù)量不夠多，這樣就無(wú)法判斷本程序的“誤判”的情況如何。5、 實(shí)習(xí)體會(huì) 在學(xué)習(xí)到相關(guān)系數(shù)的時(shí)候，一看到計(jì)算相關(guān)系數(shù)的公式的時(shí)候我就有種被折服的感覺(jué)，那個(gè)公式實(shí)在是太復(fù)雜，它使我感覺(jué)我怎么也不能記住這個(gè)復(fù)雜的公式。不過(guò)經(jīng)過(guò)這次實(shí)習(xí)后，我覺(jué)得這個(gè)公式又不是當(dāng)初認(rèn)為的那樣無(wú)法記憶。我想這次實(shí)習(xí)我最大的收獲就是對(duì)于之前在學(xué)習(xí)過(guò)程中遇到不懂知識(shí)以及難以記憶的知識(shí)都得到了加強(qiáng)，方便了我對(duì)他們的理解與記憶。 其次，在編程過(guò)程中更加鍛煉自己動(dòng)手編程的能力，使得我編程能力較之前有著一定的提高。每次編程實(shí)習(xí)都是一個(gè)積累經(jīng)驗(yàn)的過(guò)程，然后通過(guò)一次次的實(shí)習(xí)后，我最終的能力就能夠得到很好的體現(xiàn)。附件：v

6、oid Moravec(BYTE* pSrcBits, vector<FEATUREPOINT> *FeaturePoint, int m_threshold, int m_window1, int m_window2, int width, int height, int* Count1)。此部分程序與點(diǎn)特征程序中代碼一樣，因此不再重復(fù)粘貼vector<FEATUREPOINT> RemoveReplicative(vector<FEATUREPOINT> v)vector<FEATUREPOINT> ret;vector<FEATURE

7、POINT>:iterator iter = v.begin();ret.clear();ret.push_back(*iter);vector<FEATUREPOINT>:iterator iter2;BOOL b(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

8、 = TRUE;/存在break;if (b = FALSE)ret.push_back(*iter);return ret;void SaveBand(int width, int height, int byteCount, int biBitCount, LPBYTE pBits, CString SavePath)BITMAPFILEHEADER bmfh;BITMAPINFOHEADER bmih;bmfh.bfType = 0x4d42; / 0x42 = "B" 0x4d = "M"bmfh.bfReserved1 = 0;bmfh.bfR

9、eserved2 = 0;if (biBitCount = 8)bmfh.bfOffBits = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER) + 256*4;bmfh.bfSize = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER) + 256*4 + byteCount;else if (biBitCount = 24)bmfh.bfOffBits = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER);bmfh.bfSize

10、 = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER) + byteCount;elseMessageBox(NULL, "8bit or 24bit accepted", "錯(cuò)誤", MB_ICONERROR);return;bmih.biBitCount = biBitCount;bmih.biWidth = width;bmih.biHeight = height;bmih.biSize = 40;bmih.biPlanes = 1;bmih.biCompression = BI_RGB;bmi

11、h.biSizeImage = bmfh.bfSize;bmih.biClrImportant = 0;bmih.biClrUsed = 0;bmih.biXPelsPerMeter = 0;bmih.biYPelsPerMeter = 0;CFile f;/保存位圖if (f.Open(SavePath, CFile:modeCreate | CFile:modeWrite)f.Write(&bmfh, sizeof(bmfh);f.Write(&bmih, sizeof(bmih);if (biBitCount = 8)RGBQUAD rgb256;int i(0);for

12、 (; i < 256; +i)rgbi.rgbBlue = i;rgbi.rgbGreen = i;rgbi.rgbRed = i;rgbi.rgbReserved = 0;/修改顏色表255對(duì)應(yīng)的顏色rgb255.rgbBlue = 0;rgb255.rgbGreen = 0;rgb255.rgbRed = 255;f.Write(rgb, sizeof(RGBQUAD)*256);f.Write(pBits, byteCount);f.Close();elseMessageBox(NULL, "保存圖像失敗", "錯(cuò)誤", MB_ICONER

13、ROR);return;void DrawCross(BYTE* pBits, vector<FEATUREPOINT> FeaturePoint, int _width, int byteCount)vector<FEATUREPOINT>:iterator iter;int i(0), j(0);for (i = 0; i < byteCount; +i)if (pBitsi = 255)pBitsi = 254;for (iter = FeaturePoint.begin(); iter != FeaturePoint.end(); +iter)i = it

14、er->y;j = iter->x;pBitsi*_width+j = 255;pBitsi*_width+j+1 = 255;pBitsi*_width+j-1 = 255;pBits(i+1)*_width+j = 255;pBits(i-1)*_width+j = 255;pBitsi*_width+j+2 = 255;pBitsi*_width+j-2 = 255;pBits(i+2)*_width+j = 255;pBits(i-2)*_width+j = 255;void DrawCross(BYTE* pBits1, BYTE* pBits2, vector<M

15、ATCHEDPOINTS> MatchedPoints, int _width1, int _width2, int byteCount1, int byteCount2)vector<MATCHEDPOINTS>:iterator iter;int i(0), j(0);for (i = 0; i < byteCount1; +i)if (pBits1i = 255)pBits1i = 254;for (j = 0; j < byteCount2; +j)if (pBits2j = 255)pBits2j = 254;for (iter = MatchedPoi

16、nts.begin(); iter != MatchedPoints.end(); +iter)i = iter->y1;j = iter->x1;pBits1i*_width1+j = 255;pBits1i*_width1+j+1 = 255;pBits1i*_width1+j-1 = 255;pBits1(i+1)*_width1+j = 255;pBits1(i-1)*_width1+j = 255;pBits1i*_width1+j+2 = 255;pBits1i*_width1+j-2 = 255;pBits1(i+2)*_width1+j = 255;pBits1(i

17、-2)*_width1+j = 255;i = iter->y2;j = iter->x2;pBits2i*_width2+j = 255;/pBits2i*_width2+j+1 = 255;pBits2i*_width2+j-1 = 255;pBits2(i+1)*_width2+j = 255;pBits2(i-1)*_width2+j = 255;pBits2i*_width2+j+2 = 255;pBits2i*_width2+j-2 = 255;pBits2(i+2)*_width2+j = 255;pBits2(i-2)*_width2+j = 255;void Sa

18、veList(vector<FEATUREPOINT> FeaturePoint, CString txtFile)FILE* pTxtFile= fopen(txtFile, "wt");if (!pTxtFile)return;fprintf(pTxtFile, "#t特征點(diǎn)數(shù): %dnn", FeaturePoint.size();fprintf(pTxtFile, "坐標(biāo)Xt坐標(biāo)Yt興趣值n");vector<FEATUREPOINT>:iterator iter;for (iter = Feature

19、Point.begin(); iter != FeaturePoint.end(); +iter)fprintf(pTxtFile, "%dt%dt%dn", iter->x, iter->y, (int)iter->IV);fclose(pTxtFile);void SaveList(vector<MATCHEDPOINTS> MatchedPoints, CString txtFile)FILE* pTxtFile= fopen(txtFile, "wt");if (!pTxtFile)return;fprintf(pT

20、xtFile, "#同名點(diǎn)對(duì)數(shù): %dnn",MatchedPoints.size();vector<MATCHEDPOINTS>:iterator iter;for (iter = MatchedPoints.begin(); iter != MatchedPoints.end(); +iter)fprintf(pTxtFile, "%dt%dt%dt%dt%fn",iter->x1, iter->y1, iter->x2, iter->y2, iter->coefficient);fclose(pTxtFil

21、e);int _tmain(int argc, TCHAR* argv, TCHAR* envp)if (argc != 15):MessageBox(NULL,"啟動(dòng)參數(shù)不正確,請(qǐng)從GUI.exe或命令行參數(shù)啟動(dòng)","",MB_ICONASTERISK);return FAILURE;system("Mode con: COLS=50 LINES=150");system("title Image_Match");int nRetCode = 0;if (!AfxWinInit(:GetModuleHandle(

22、NULL), NULL, :GetCommandLine(), 0)cerr << _T("Fatal Error: MFC initialization failed") << endl;nRetCode = 1;return nRetCode;CString srcFileLeft, srcFileRight, dstFileLeft, dstFileRight, txtFile;CString threshold_Left,window1_Left,window2_Left;CString threshold_Right,window1_Rig

23、ht,window2_Right;CString CoRelated,window3,flag1;/啟動(dòng)參數(shù)srcFileLeft = argv1;srcFileRight = argv2;dstFileLeft = argv3;dstFileRight = argv4;txtFile = argv5;threshold_Left = argv6;/閾值window1_Left = argv7;window2_Left = argv8;threshold_Right = argv9;window1_Right = argv10;window2_Right = argv11;CoRelated

24、= argv12;window3 = argv13;flag1 = argv14;/判斷是否在創(chuàng)建結(jié)果文件后顯示int m_threshold_L = atoi(threshold_Left);int m_window1_L = atoi(window1_Left);int m_window2_L = atoi(window2_Left);int m_threshold_R = atoi(threshold_Right);int m_window1_R = atoi(window1_Right);int m_window2_R = atoi(window2_Right);double m_Co

25、Related = (double)atoi(CoRelated) / 100;int m_window3 = atoi(window3);int m_flag1 = atoi(flag1);if (GetFileAttributes(srcFileLeft) = -1)CString err;err.Format("文件:%s不存在", srcFileLeft);:MessageBox(NULL,err, "Error", MB_ICONERROR);return FAILURE;if (GetFileAttributes(srcFileRight)

26、= -1)CString err;err.Format("文件:%s不存在", srcFileRight);:MessageBox(NULL,err, "Error", MB_ICONERROR);return FAILURE;/打開(kāi)2幅影像FILE* pSrcFileLeft = NULL;pSrcFileLeft = fopen(srcFileLeft,"rb");if(pSrcFileLeft = NULL)printf("原始左片影像打開(kāi)失敗n");return FAILURE;FILE* pSrcFile

27、Right = NULL;pSrcFileRight = fopen(srcFileRight,"rb");if(pSrcFileRight = NULL)printf("原始右片影像打開(kāi)失敗n");return FAILURE;/讀取左片數(shù)據(jù)BITMAPFILEHEADER bmfhLeft;BITMAPINFOHEADER bmihLeft;fread(&bmfhLeft, sizeof(BITMAPFILEHEADER), 1, pSrcFileLeft);fread(&bmihLeft, sizeof(BITMAPINFOHEAD

28、ER), 1, pSrcFileLeft);if(bmihLeft.biBitCount != 8)printf("只支持8位位圖: 左片n");return FAILURE;DWORD widthLeft = bmihLeft.biWidth;DWORD heightLeft = bmihLeft.biHeight;DWORD _widthLeft = (widthLeft*8 + 31)/32*4;/saved widthLeft in filesDWORD pixelCountLeft = _widthLeft*heightLeft;fseek(pSrcFileLef

29、t, 256*sizeof(RGBQUAD), SEEK_CUR);BYTE* pSrcBitsLeft = new BYTEpixelCountLeft;/new 1fread(pSrcBitsLeft, sizeof(BYTE), pixelCountLeft, pSrcFileLeft);fclose(pSrcFileLeft);pSrcFileLeft = NULL;/讀取右片數(shù)據(jù)BITMAPFILEHEADER bmfhRight;BITMAPINFOHEADER bmihRight;fread(&bmfhRight, sizeof(BITMAPFILEHEADER), 1,

30、 pSrcFileRight);fread(&bmihRight, sizeof(BITMAPINFOHEADER), 1, pSrcFileRight);if(bmihRight.biBitCount != 8)printf("只支持8位位圖: 右片n");return FAILURE;DWORD widthRight = bmihRight.biWidth;DWORD heightRight = bmihRight.biHeight;DWORD _widthRight = (widthRight*8 + 31)/32*4;/saved widthRight in

31、 filesDWORD pixelCountRight = _widthRight*heightRight;fseek(pSrcFileRight, 256*sizeof(RGBQUAD), SEEK_CUR);BYTE* pSrcBitsRight = new BYTEpixelCountRight;/new 2fread(pSrcBitsRight, sizeof(BYTE), pixelCountRight, pSrcFileRight);fclose(pSrcFileRight);pSrcFileRight = NULL;/開(kāi)始計(jì)算左片特征點(diǎn)/計(jì)算左片的興趣點(diǎn)存儲(chǔ)于該vector中ve

32、ctor<FEATUREPOINT> FeaturePointLeft;vector<FEATUREPOINT> FeaturePointRight;int Count1(0);int Count2(0);Moravec(pSrcBitsLeft, &FeaturePointLeft, m_threshold_L, m_window1_L, m_window2_L,widthLeft, heightLeft, &Count1);Moravec(pSrcBitsRight, &FeaturePointRight, m_threshold_R, m_

33、window1_R, m_window2_R,widthRight, heightRight, &Count2);cout << FeaturePointLeft.size() << endl;cout << FeaturePointRight.size() << endl;FeaturePointLeft = RemoveReplicative(FeaturePointLeft);FeaturePointRight = RemoveReplicative(FeaturePointRight);cout << "Fe

34、aturePointLeft.size() = " << FeaturePointLeft.size() << endl;cout << "FeaturePointRight.size() = " << FeaturePointRight.size() << endl;SaveList(FeaturePointLeft, "d:list1.txt");SaveList(FeaturePointRight, "d:list2.txt");/特征點(diǎn)計(jì)算完畢，存儲(chǔ)于vec

35、tor中/遍歷vector, 在右片中搜尋相關(guān)系數(shù)最大的窗口vector<FEATUREPOINT>:iterator iter1;vector<FEATUREPOINT>:iterator iter2;int k = m_window3 / 2;int i(0), j(0), m(0), n(0);BYTE* pTemplate = new BYTEm_window3*m_window3;/new 5ZeroMemory(pTemplate, sizeof(BYTE)*m_window3*m_window3);BYTE* pTarget = new BYTEm_win

36、dow3*m_window3;/new 6ZeroMemory(pTarget, sizeof(BYTE)*m_window3*m_window3);vector<MATCHEDPOINTS> MatchedPoints;float MaxR(0);float R(0);double A(0), B(0), C(0), D(0), E(0);/計(jì)算相關(guān)系數(shù)的中間量int lMaxWidth;int lMaxHeight;int Count3(0), Count4(0);for (iter1 = FeaturePointLeft.begin(); iter1 != FeaturePo

37、intLeft.end(); +iter1)cout << endl << "Left point ID: " << Count2+ << endl;/計(jì)算該點(diǎn)對(duì)應(yīng)的特征點(diǎn)B = 0.0;D = 0.0;MaxR = 0.0;R = 0.0;if (iter1->x < k | iter1->y < k | iter1->x > heightLeft - k | iter1->y > widthLeft - k)cout << "Left point abo

38、rt: " << Count3+ << endl;continue;/逐個(gè)搜索/將左片小窗口內(nèi)的數(shù)據(jù)拷貝/可先求出B,Dfor (i = 0; i < m_window3; +i)for(j = 0; j < m_window3 ; +j)/B = 小窗口所有像素灰度值求和/D = .的平方再求和B += pSrcBitsLeft(i+iter1->y)*_widthLeft + (j+iter1->x);D += pSrcBitsLeft(i+iter1->y)*_widthLeft + (j+iter1->x) * pS

39、rcBitsLeft(i+iter1->y)*_widthLeft + (j+iter1->x);Count4 = 0;/逐個(gè)小窗口計(jì)算Rfor (iter2 = FeaturePointRight.begin(); iter2 != FeaturePointRight.end(); +iter2)if (iter2->x < k | iter2->y < k | iter2->x > heightRight - k | iter2->y > widthRight - k)cout << "Right point abort: " << Count4+ << endl;continue;A = 0.0;C = 0.0;E = 0.0;j = iter2->x;/i,j為中心點(diǎn)i = iter2->y;j -= k;/i,j為頂點(diǎn)i -= k;int centerX = j + k;int centerY = i + k;/這里的centerX，centerY就是右片某個(gè)小窗口的中心/求A,C,Efor (m = 0; m < m_window3; +m)for