ASIFT+OpenCV圖像特征匹配實戰

OpenCV包含頭文件：

1. #include "cv.h"
2. #include "highgui.h"
3. #include "cxcore.h"

1. if (!m_pImage1||!m_pImage2)
2. {
3. AfxMessageBox("please,select 2 images!");
4. return;
5. }
7. UpdateData(TRUE);
9. CvSize sz1 = cvSize(m_pImage1->width,m_pImage1->height);
10. CvSize sz2 = cvSize(m_pImage2->width,m_pImage2->height);
12. CvScalar s;
14. IplImage *gimg1 = cvCreateImage(sz1,IPL_DEPTH_8U,1);
15. cvCvtColor(m_pImage1,gimg1,CV_BGR2GRAY);
17. IplImage *gimg2 = cvCreateImage(sz2,IPL_DEPTH_8U,1);
18. cvCvtColor(m_pImage2,gimg2,CV_BGR2GRAY);
20. size_t w1, h1;
22. w1 = gimg1->width;
23. h1 = gimg1->height;
25. float * iarr1 = new float[w1*h1];
27. for(int i=0;i<h1;i++)
28. {
29. for(int j=0;j<w1;j++)
30. {
31. s=cvGet2D(gimg1,i,j);
32. iarr1[i*w1+j] = s.val[0];
33. }
34. }
36. vector<float> ipixels1(iarr1, iarr1 + w1 * h1);
38. delete [] iarr1;
40. size_t w2, h2;
42. w2 = gimg2->width;
43. h2 = gimg2->height;
45. float * iarr2 = new float[w2*h2];
47. for(int i=0;i<h2;i++)
48. {
49. for(int j=0;j<w2;j++)
50. {
51. s=cvGet2D(gimg2,i,j);
52. iarr2[i*w2+j] = s.val[0];
53. }
54. }
56. vector<float> ipixels2(iarr2, iarr2 + w2 * h2);
58. delete [] iarr2;
60. float wS = IM_X;
61. float hS = IM_Y;
63. float zoom1=0, zoom2=0;
64. int wS1=0, hS1=0, wS2=0, hS2=0;
65. vector<float> ipixels1_zoom, ipixels2_zoom;
67. if (!m_bOrininal)
68. {
69. if (m_lWidth==0 || m_lHeight == 0)
70. return;
72. wS = m_lWidth;
73. hS = m_lHeight;
75. float InitSigma_aa = 1.6;
77. float fproj_p, fproj_bg;
78. char fproj_i;
79. float *fproj_x4, *fproj_y4;
80. int fproj_o;
82. fproj_o = 3;
83. fproj_p = 0;
84. fproj_i = 0;
85. fproj_bg = 0;
86. fproj_x4 = 0;
87. fproj_y4 = 0;
89. float areaS = wS * hS;
91. // Resize image 1
92. float area1 = w1 * h1;
93. zoom1 = sqrt(area1/areaS);
95. wS1 = (int) (w1 / zoom1);
96. hS1 = (int) (h1 / zoom1);
98. int fproj_sx = wS1;
99. int fproj_sy = hS1;
101. float fproj_x1 = 0;
102. float fproj_y1 = 0;
103. float fproj_x2 = wS1;
104. float fproj_y2 = 0;
105. float fproj_x3 = 0;
106. float fproj_y3 = hS1;
108. /* Anti-aliasing filtering along vertical direction */
109. if ( zoom1 > 1 )
110. {
111. float sigma_aa = InitSigma_aa * zoom1 / 2;
112. GaussianBlur1D(ipixels1,w1,h1,sigma_aa,1);
113. GaussianBlur1D(ipixels1,w1,h1,sigma_aa,0);
114. }
116. // simulate a tilt: subsample the image along the vertical axis by a factor of t.
117. ipixels1_zoom.resize(wS1*hS1);
118. fproj (ipixels1, ipixels1_zoom, w1, h1, &fproj_sx, &fproj_sy, &fproj_bg, &fproj_o, &fproj_p,
119. &fproj_i , fproj_x1 , fproj_y1 , fproj_x2 , fproj_y2 , fproj_x3 , fproj_y3, fproj_x4, fproj_y4);
122. // Resize image 2
123. float area2 = w2 * h2;
124. zoom2 = sqrt(area2/areaS);
126. wS2 = (int) (w2 / zoom2);
127. hS2 = (int) (h2 / zoom2);
129. fproj_sx = wS2;
130. fproj_sy = hS2;
132. fproj_x2 = wS2;
133. fproj_y3 = hS2;
135. /* Anti-aliasing filtering along vertical direction */
136. if ( zoom1 > 1 )
137. {
138. float sigma_aa = InitSigma_aa * zoom2 / 2;
139. GaussianBlur1D(ipixels2,w2,h2,sigma_aa,1);
140. GaussianBlur1D(ipixels2,w2,h2,sigma_aa,0);
141. }
143. // simulate a tilt: subsample the image along the vertical axis by a factor of t.
144. ipixels2_zoom.resize(wS2*hS2);
145. fproj (ipixels2, ipixels2_zoom, w2, h2, &fproj_sx, &fproj_sy, &fproj_bg, &fproj_o, &fproj_p,
146. &fproj_i , fproj_x1 , fproj_y1 , fproj_x2 , fproj_y2 , fproj_x3 , fproj_y3, fproj_x4, fproj_y4);
147. }
148. else
149. {
150. ipixels1_zoom.resize(w1*h1);
151. ipixels1_zoom = ipixels1;
152. wS1 = w1;
153. hS1 = h1;
154. zoom1 = 1;
156. ipixels2_zoom.resize(w2*h2);
157. ipixels2_zoom = ipixels2;
158. wS2 = w2;
159. hS2 = h2;
160. zoom2 = 1;
161. }
163. int num_of_tilts1 = m_lTilts1;
164. int num_of_tilts2 = m_lTilts2;
166. int verb = 0;
167. // Define the SIFT parameters
168. siftPar siftparameters;
169. default_sift_parameters(siftparameters);
171. vector< vector< keypointslist > > keys1;
172. vector< vector< keypointslist > > keys2;
174. int num_keys1=0, num_keys2=0;
176. SetWindowText("Computing keypoints on the two images...");
178. CString str1,str2;
180. time_t tstart, tend1,tend2;
181. tstart = time(0);
182. DWORD dstart = GetTickCount();
184. num_keys1 = compute_asift_keypoints(ipixels1_zoom, wS1, hS1, num_of_tilts1, verb, keys1, siftparameters);
186. tend1 = time(0);
188. m_lKeyNum1 = num_keys1;
189. UpdateData(FALSE);
191. str1.Format("Img1 Keypoints computation accomplished in %f s",difftime(tend1, tstart));
192. SetWindowText(str1);
194. num_keys2 = compute_asift_keypoints(ipixels2_zoom, wS2, hS2, num_of_tilts2, verb, keys2, siftparameters);
196. tend2 = time(0);
198. m_lKeyNum2 = num_keys2;
199. UpdateData(FALSE);
201. str2.Format("Img2 Keypoints computation accomplished in %f s ,Matching the keypoints...",difftime(tend2, tstart));
202. SetWindowText(str2);
204. //// Match ASIFT keypoints
205. int num_matchings;
206. matchingslist matchings;
208. tstart = time(0);
209. num_matchings = compute_asift_matches(num_of_tilts1, num_of_tilts2, wS1, hS1, wS2,
210. hS2, verb, keys1, keys2, matchings, siftparameters);
211. tend1 = time(0);
212. DWORD dSpan = GetTickCount() - dstart;
214. cout << "Keypoints matching accomplished in " << difftime(tend1, tstart) << " seconds." << endl;
215. str2.Format("Keypoints matching accomplished in %f s",difftime(tend1, tstart));
216. SetWindowText(str2);
218. m_lMatches = num_matchings;
219. UpdateData(FALSE);
221. str1.Format("Total time used:%d ms",dSpan);
223. AfxMessageBox(str1);
225. cvRelease((void**)&gimg1);
226. cvRelease((void**)&gimg2);

運行界面：