- 1、本文档共17页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
摄影测量学点特征提取模拟实践
辽宁工程技术大学模 拟 实 践教学单位 测绘学院专 业 测绘工程名 称 点特征提取模拟实践 班 级 测绘13-3 组 长 秦小茜 组 员 柳少冬 倪贺星 刘建琪指导教师 徐辛超 一、目的 理解影像中每个像素灰度值的概念;理解点特征在灰度方面的特点;掌握常用的点特征提取算子,及其对应的点特征提取方法,包括:Moravec算子、Forstner算子。二、要求 运用VC++、C#、MATLAB其中一种自己擅长的语言,实现两种点特征提取算法。三、方法与步骤 (一) Moravec算子提取过程a. 读取灰度影像,得到影像的灰度矩阵;b. 计算各像元的兴趣值 IV:c. 给定一经验阈值,将兴趣值()大于阈值的点作为候选点; d. 选取候选点中的极值点作为特征点;e. 将特征点在原始图像上进行标记、展示。 (二) Forstner算子提取过程a. 计算各像素的Robert’s梯度 ;b. 计算ll(如55或更大)窗口中灰度的协方差矩阵 ;其中:c. 计算兴趣值q与w ;式中:DetN代表矩阵N之行列式,trN代表矩阵N之迹。d. 确定待选点 ;当同时,该像元为待选点。e. 选取极值点。 即在一个适当窗口中选择最大的待选点。f. 将特征点在原始图像上进行标记、展示。、流程图 实验内容及成果1、Moravec算子程序相关内容(1)Moravec算子程序源代码:clcclear% I=imread(C:\Users\LBS\Desktop\Moravec MATLAB\girl.bmp);I=imread(C:\Users\Administrator\Desktop\点特征提取\Moravec MATLAB\left缩小.bmp);originalmap=I;% if isrgb(I)% I=rgb2gray(I); %将彩色图像转换为灰度图像% end I=double(I); %转换为双精度图像[row,column]=size(I);w1 = floor(5/2);%确定窗口大小w2 = 5;step=w1;ir=10000;%设置阈值num=0;%得到特征点的坐标xc=[];yc=[];tic;for i=w1+1:1:row-w1-1 if i+w2row-w1-1 continue; end for j=w1+1:1:column-w1-1 if j+w2column-w1-1 continue; end interest=0;%兴趣值 v1=0;v2=0;v3=0;v4=0; for k=-step:step-1 %计算纵向相邻像素灰度差的平方和 v1=v1+(I(i+k,j)-I(i+k+1,j))^2; %计算125度斜线方向相邻像素灰度差的平方和 v2=v2+(I(i+k,j+k)-I(i+k+1,j+k+1))^2; %计算横相邻像素灰度差的平方和 v3=v3+(I(i,j+k)-I(i,j+k+1))^2; %计算45度斜线方向相邻像素灰度差的平方和 v4=v4+(I(i+k,j-k)-I(i+k+1,j-k-1))^2; end %取v1,v2,v3,v4中最小者作为像素(c,r)的兴趣值 interest=min([v1,v2,v3,v4]); if interestir%//最大兴趣值与阈值比较 num=num+1; pointx(num)=j; pointy(num)=i; yc=[yc,i]; xc=[xc,j]; end endend %图像中显示特征点figure(4);imshow(originalmap);hold on;plot(xc,yc,R*);axis on;title(图像中显示特征点);xlabel(图像的列数);ylabel(图像的行数); toc;t=toc;disp([本程序的运行时间为,num2str(t),秒。]); 图一Moravec算子程序的相关截图 原始灰度影像 图二原始灰度影像 (3)程序运行结果 图三Moravec算子程序运行结果Forstner算子程序的相关内容Forstner算子程序源代码:clcclearI=imread(C:\Users\Administrator\Desktop\点特征提取\Moravec MATLAB\left缩小.bmp);% I=imread(C:\Users\LBS\Desktop\Lenna.bmp);originalmap=I;% if isrgb(I)% I=rgb2gray(I); %将彩色图像转换为灰度图像% end I=double(I); %转换为双精度图像[row,column]=size(I);w1 = 5;%确定窗口大小w2 = 5;k=floo
您可能关注的文档
- 总结材料(样本—简化)汪村互通立交.doc
- 急诊服务流程(急诊→医技检查→住院→手术→介入)与规范.doc
- 急诊科专科知识培训试题及答案.docx
- 慢性阻塞性肺疾病的家庭康复护理(修改版).doc
- 我们是被“他者”认同的符号所包裹着的僵尸.doc
- 恩格斯晚年的历史观对马克思主义唯物史观的发展.doc
- 我们都可以成为研究者.doc
- 思科协作解决方案.docx
- 建筑效果图后期ps处理.docx
- 我国信息化发展进入互联网时代答案.docx
- 应收账款重要实质性程序18zkzh lsnjc sj013.pdf
- 文本案例章概率1 acjc ii q13chapter 02.pdf
- 快速数据位移寄存器motorola摩托罗拉sn74ls95.pdf
- 十一章其他税种纳税申报审核提高练习全swsw lx1101.pdf
- 浅谈谨慎性原则在会计实务中的应用.docx
- 浅论企业会计信息质量的问题与对策-以淘宝网为例.doc
- 企业所得税纳税的会计处理现状1.doc
- 企业人力资源成本控制问题研究-以京东为例.doc
- 浅谈博维(杭州)知识产权代理有限公司成本管理.doc
- 宁波华尊汽车服务公司财务管理中的问题及对策——以宁波华尊汽车服务有限公司为例.doc
文档评论(0)