遥感提取特征点Word格式.docx

1、imgn=zeros（h,w）;n=4;for y=1+n:h-n for x=1+n:w-n sq=img（y-n:y+n,x-n:x+n）; V=zeros（1,4）; for i=2:2*n+1 %垂直，水平，对角，反对角四个方向领域灰度差的平方和 V（1）=V（1）+（sq（i,n+1）-sq（i-1,n+1）2; V（2）=V（2）+（sq（n+1,i）-sq（n+1,i-1）2; V（3）=V（3）+（sq（i,i）-sq（i-1,i-1）2; V（4）=V（4）+（sq（i,（2*n+1）-（i-1）-sq（i-1,（2*n+1）-（i-2）2; end pix=min（V）;

2、%四个方向中选最小值 imgn（y,x）=pix;endT=mean（imgn（: %设阈值，小于均值置零ind=find（imgnImage（i,j）-t & Image（i-1,j-1） Image（i-1,j） Image（i-1,j+1） Image（i,j-1） Image（i,j+1） Image（i+1,j-1） Image（i+1,j） Image（i+1,j+1）=2 & nlike CRFmax CRFmax = CRF（i,j）;%CRFmaxcount = 0; % 用来记录角点的个数t=0.01;% 下面通过一个3*3的窗口来判断当前位置是否为角点 if Corner

3、（i,j）=1 %只关注候选点的八邻域 t*CRFmax & CRF（i,j） CRF（i-1,j-1） . & CRF（i-1,j） & CRF（i-1,j+1） . CRF（i,j-1） & CRF（i,j+1） . CRF（i+1,j-1） & CRF（i+1,j）. CRF（i+1,j+1） count=count+1;%这个是角点，count加1 else % 如果当前位置（i,j）不是角点，则在Corner（i,j）中删除对该候选角点的记录 Corner（i,j） = 0;% disp（角点个数% disp（count）figure,imshow（Image）; % display

4、 Intensity Imagehold on;% toc（t1）for j=boundary: column_ave=0; row_ave=0; k=0; if Corner（i,j）=1 for x=i-3:i+3 %7*7邻域 for y=j-3:j+3 if Corner（x,y）=1% 用算数平均数作为角点坐标，如果改用几何平均数求点的平均坐标，对角点的提取意义不大 row_ave=row_ave+x; column_ave=column_ave+y; k=k+1; if k0 %周围不止一个角点 plot（ column_ave/k,row_ave/k ,g.三 运行结果 1001

5、特征点 1002特征点二、 （二）（三）Susan算子一 算法公式（1）借助图3-1来解释Susan检测的原理，其中图片是白色背景，有一个颜色比较暗淡的矩形。用一个圆形模板在图像上移动，若模板内的像素灰度与模板中心的像素（被称为核Nucleus）灰度值小于一定的阈值，则认为该点与核Nucleus具有相同的灰度，满足该条件的像素组成的区域就称为USAN。在图片上有5个圆形区域。圆形区域表示的是掩码区域。把圆形区域内的每一个位置的像素值与圆心处的像素值相比较，那么圆中的的像素可以分为两类，一类是像素值与圆心处的像素值相近的，另一类是像素值与圆心的处的像素值相差比较大的。 图3-1 图3-2如果将模

6、板中各个像素的灰度都与模板中心的核像素的灰度进行比较，那么就会发现总有一部分模板区域和灰度与核像素的灰度相同或相似，这部分区域可以称为USAN。USAN区域包含很多与图像结构有关的信息。利用这种区域的尺寸、重心、二阶矩的分析，可以得到图像中的角点，边缘等信息。从上图所示，当核像素处在图像中的灰度一致区域时，USAN的面积会达到最大。第e个模板就是属于这种情况。（2）Susan进行角点检测时，遵循了常规的思路：使用一个窗口在图像上逐点滑动，在每一个位置上计算出一个角点量，再进行局部极大值抑制，得到最终的角点。我们这里使用的窗口是圆形窗口，最小的窗口是33的，此次使用的是37个像素的圆形窗口，如图

7、3-2。（3）在角点检测中，有两种类型的阈值，一种用来约束角点的数量，另一种用来约束角点的质量。当然，一个阈值不能完全做到只影响质量或数量，只是会有一个侧重点。阈值g是角点质量的，尽管也会影响数量，但是相对来说更侧重于影响质量（角点的形状）。例如，g值减小，那么Susan会更加侧重于检测到更加“尖锐”的角点，所以，可以更加自己的实际需求来确定阈值g；而阈值t，是角点的数量，当t减小时，会检测到更多的角点，所以，阈值t可以在不影响角点质量的情况下，控制检测到的角点的数量，如果图像的对比度比较低，可以修改t值以适应变化。 下面简单叙述下利用Susan算子检测角点的步骤：1:利用圆形模板遍历图像，计

8、算每点处的USAN值。设置一阈值g，一般取值为1/2（Max（n）， 也即取值为USAN最大值的一半，进行阈值化，得到角点响应。使用非极大值抑制来寻找角点。通过上面的方式得到的角点，存在很大伪角点。为了去除伪角点，Susan算子可以由以下方法实现：计算USAN区域的重心，然后计算重心和模板中心的距离，如果距离较小则不是正确的角点；判断USAN区域的重心和模板中心的连线所经过的像素都是否属于USAN区域的像素，如果属于那么这个模板中心的点就是角点。clc;img=imread（img=rgb2gray（img）;imshow（img）;m n=size（img）;img=double（img）;

9、t=45; %模板中心像素灰度和周围灰度差别的阈值，自己设置usan=; %当前像素和周围在像素差别在t以下的个数%这里用了37个像素的模板for i=4:m-3 %没有在外围扩展图像，最终图像会缩小 for j=4:n-3 tmp=img（i-3:i+3,j-3:j+3）; %先构造7*7的模板，49个像素 c=0; for p=1:7 for q=1: if （p-4）2+（q-4）2=12 %在其中筛选，最终模板类似一个圆形 % usan（k）=usan（k）+exp（-（img（i,j）-tmp（p,q）/t）6）; if abs（img（i,j）-tmp（p,q）t %判断灰度是否相

10、近，t是自己设置的 c=c+1; usan=usan c;g=2*max（usan）/3; %确定角点提取的数量，值比较高时会提取出边缘，自己设置for i=1:length（usan） if usan（i）max（max（imgn（i-1,j-1:j+1） imgn（i,j-1） imgn（i,j+1） max（imgn（i+1,j-1:j+1）; re（i,j）=1; re（i,j）=0;imshow（re=1）;（四）三种方法的比较：（1）Moravec算子对斜边缘的响应很强，因为只考虑了每隔45度的方向变化，而没有在全部的方向上进行考虑；同时由于窗口函数是一个二值函数,不管像素点离中心

11、点的距离，赋于一样的权重，因此对噪声响应也一般，最终对角点的定位较准确。（2）Harris算子是一种有效的点特征提取算子,其优点总结起来有:计算简单:Harris算子中只用到灰度的一阶差分以及滤波,操作简单。2:提取的点特征均匀而且合理:Harris算子对图像中的每个点都计算其兴趣值,然后在邻域中选择最优点。实验表明,在纹理信息丰富的区域,Harris算子可以提取出大量有用的特征点,而在纹理信息少的区域,提取的特征点则较少。3:稳定:Harris算子的计算公式中只涉及到一阶导数，因此对图像旋转、灰度变化、噪声影响和视点变换不敏感,它也是比较稳定的一种点特征提取算子。Harris算子的局限性有：

12、它对尺度很敏感，不具有尺度不变性。提取的角点是像素级的。（3）Susan算子的特点有：在用Susan算子对边缘和角点进行检测时不需要计算微分，这使得Susan算子对噪声更加鲁棒。Susan检测算子能提供不依赖于模板尺寸的边缘精度。换句话说，最小USAN区域面积的计算是个相对的概念，与模版尺寸无关，所以Susan边缘算子的性能不受模版尺寸影响。控制参数的选择很简单，且任意性小，容易实现自动化选取。（5）心得体会通过本次实验，我对特征点提取方法的计算原理及实现过程有了深刻的认识，在本次实验中，我遇到了很多困难，但是在同学们的帮助下，我们互相商讨，这些问题都一一得到了解决，在解决困难的过程中的编程能力得到了提高，对其所涉及到的知识的印象也得到了加深。总之，感谢老师给了我这次锻炼自己的机会，之后还要继续学习研究MATLAB，提升自己的编程能力。