单片机论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 1792|回复: 0
收起左侧

MATLAB程序 Harris角点检测算法

[复制链接]
51hei小小 发表于 2017-6-3 01:45 | 显示全部楼层 |阅读模式
  1. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  2. %特征点的匹配,主要应用 harris 角点的检测,match单向匹配  函数
  3. %适合于有白边的图像,因为,在加窗滤波的时候,没有限定范围的哈,尽量保证角点不在边上
  4. clc,clear all;
  5. filename1 = '..\.\photo9\cal1-L.txt';
  6. filename2 = '..\.\photo9\cal1-R.txt';
  7. b1 = imread('..\.\photo9\cal1-L.bmp');
  8. b2 = imread('..\.\photo9\cal1-R.bmp'); %double的作用很大的啊
  9. a1 = imadjust(b1);
  10. a2 = imadjust(b2);
  11. %subplot(1,2,1);
  12. %imshow(a1);
  13. %subplot(1,2,2);
  14. %imshow(a2);   
  15. %title('原来的图像');

  16. [result1,cnt1,c1,r1]=harris(a1);%角点检测,得到原始的焦点位置图result
  17. [result2,cnt2,c2,r2]=harris(a2);

  18. figure;
  19. subplot(1,2,1);
  20. imshow(a1);
  21. hold on;
  22. plot(r1,c1,'g.');
  23. subplot(1,2,2);
  24. imshow(a2);
  25. hold on;
  26. plot(r2,c2,'g.');
  27. title('图1,2的角点图');

  28. l=[r1 c1];
  29. r=[r2 c2];
  30. fid1=fopen(filename1,'wt');
  31. for i=1:cnt1
  32.     for j=1:2
  33.        if j==2
  34.          fprintf(fid1,'%g\n',l(i,j));
  35.       else
  36.         fprintf(fid1,'%g\t',l(i,j));
  37.        end
  38.     end
  39. end
  40. fclose(fid1);

  41. fid2=fopen(filename2,'wt');
  42. for i=1:cnt2
  43.     for j=1:2
  44.       if j==2
  45.          fprintf(fid2,'%g\n',r(i,j));
  46.       else
  47.         fprintf(fid2,'%g\t',r(i,j));
  48.        end
  49.     end
  50. end
  51. fclose(fid2);

  52. %res2=match(a1,cnt1,r1,c1,a2,cnt2,r2,c2);%从result1中开始搜索在result2中可能达到的
  53. %[r22,c22]=find(res2==1);
  54. %[m22,n22]=size(r22);
  55. %cnt22=m22;

  56. %res1=match(a2,cnt22,r22,c22,a1,cnt1,r1,c1);%反向搜索res2--result1
  57. %res1=and(res1,result1);   %保证反向匹配不会出现不可能的点
  58. %[r11,c11]=find(res1==1);
  59. %[m11,n11]=size(r11);
  60. %cnt11=m11;

  61. %res22=match(a1,cnt11,r11,c11,a2,cnt22,r22,c22);%从res1中开始搜索在res2中可能达到的
  62. ……………………

  63. …………限于本文篇幅 余下代码请从51黑下载附件…………
复制代码


完整版本下载:
http://www.51hei.com/bbs/dpj-86571-1.html

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|单片机论坛 |51黑电子论坛单片机 联系QQ:125739409;技术交流QQ群582644647

Powered by 单片机教程网

快速回复 返回顶部 返回列表