Moravec角點檢測算子

Moravec角點檢測算子

Moravec 在1981年提出Moravec角點檢測算子[1]，并將它應用于立體匹配。

首先, 計算每個像素點的興趣值, 即以該像素點為中心, 取一個w*w(如:5x5)的方形窗口, 計算0度、45度、90度、135度四個方向灰度差的平方和, 取其中的最小值作為該像素點的興趣值。

圖1- 1　以3x3為例　黑色窗口為Ｉ(x,y)　紅色窗口為I(x+u,y+v)

其中四種移位 (u,v) = (1,0), (1,1), (0,1), (-1, 1).w(x,y)為方形二值窗口,若像素點在窗口內，則取值為1, 否則為0。

其次,根據實際圖像設定一個閾值, 遍歷圖像以興趣值大于該閾值的點為候選點。

最后, 選一個一定大小的滑動窗口?, 讓該窗口遍歷灰度圖象, 在此過程中取窗口中興趣值最大的候選點為特征點, 算法結束。

?

?

圖1- 2?? Moravec角點檢測算子對簡單圖像的響應

Moravec角點檢測算子對斜邊緣的響應很強，因為只考慮了每隔45度的方向變化，而沒有在全部的方向上進行考慮；同時由于窗口函數是一個二值函數,不管像素點離中心點的距離，賦于一樣的權重，因此對噪聲響應也較強。最終對角點的定位也不是很準確.

[1]??? Moravec, H. 1981. Rover visual obstacle avoidance. In?International　Joint Conference on Artificial Intelligence, Vancouver, Canada,pp. 785–790.

?

opencv代碼:

#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#define M_YUZHI 400
using namespace cv;
using namespace std;void interestValue(int *a,int &v);//a為5*5窗口數組,v為最小值作為興趣點
void getmax(int *m_a,int &m_max,int &m_num);//m_a為滑動窗口數組,m_max為尋找的滑動窗口內最大值,m_num為最大值在數組中位置.
int  main(int argc,char **argv)
{Mat mat_src=imread(argv[1],0);Mat mat_c;cvtColor(mat_src,mat_c,CV_GRAY2BGR);int width=mat_src.cols;int height=mat_src.rows;int *temp=new int[width*height];int a[25];//5*5的模板int b[49];//7*7的模板int i,j;for(i=0;i<height;i++){for(j=0;j<width;j++){if(i<2||i>height-3||j<2||j>width-3)//處理邊界,設置為0{temp[i*width+j]=0;}else {a[0]=mat_src.at<uchar>(i-2,j-2);a[1]=mat_src.at<uchar>(i-2,j-1);a[2]=mat_src.at<uchar>(i-2,j);a[3]=mat_src.at<uchar>(i-2,j+1);a[4]=mat_src.at<uchar>(i-2,j+2);a[5]=mat_src.at<uchar>(i-1,j-2);a[6]=mat_src.at<uchar>(i-1,j-1);a[7]=mat_src.at<uchar>(i-1,j);a[8]=mat_src.at<uchar>(i-1,j+1);a[9]=mat_src.at<uchar>(i-1,j+2);a[10]=mat_src.at<uchar>(i,j-2);a[11]=mat_src.at<uchar>(i,j-1);a[12]=mat_src.at<uchar>(i,j);a[13]=mat_src.at<uchar>(i,j+1);a[14]=mat_src.at<uchar>(i,j+2);a[15]=mat_src.at<uchar>(i+1,j-2);a[16]=mat_src.at<uchar>(i+1,j-1);a[17]=mat_src.at<uchar>(i+1,j);a[18]=mat_src.at<uchar>(i+1,j+1);a[19]=mat_src.at<uchar>(i+1,j+2);a[20]=mat_src.at<uchar>(i+2,j-2);a[21]=mat_src.at<uchar>(i+2,j-1);a[22]=mat_src.at<uchar>(i+2,j);a[23]=mat_src.at<uchar>(i+2,j+1);a[24]=mat_src.at<uchar>(i+2,j+2);interestValue(a,temp[i*width+j]);}}}int m_w1=width/7;//滑動窗口寬int m_w2=width%7;int m_h1=height/7;int m_h2=height%7;int m_r=0;int m_c=0;for(i=0;i<m_h1;i++){for(j=0;j<m_w1;j++){for(int k=0;k<7;k++){b[7*k]=temp[i*width*7+j*7+k*width];b[7*k+1]=temp[i*width*7+j*7+1+k*width];b[7*k+2]=temp[i*width*7+j*7+2+k*width];b[7*k+3]=temp[i*width*7+j*7+3+k*width];b[7*k+4]=temp[i*width*7+j*7+4+k*width];b[7*k+5]=temp[i*width*7+j*7+5+k*width];b[7*k+6]=temp[i*width*7+j*7+6+k*width];}int m_tempnum=0;//獲取滑動區域的最大值int m_num=0;//位置getmax(b,m_tempnum,m_num); if(m_tempnum>M_YUZHI){m_r=m_num/7;m_c=m_num%7;m_r=7*i+m_r;m_c=7*j+m_c;circle(mat_c,Point(m_c,m_r),2,Scalar(0,0,255),1,8);}}}imshow("mat_c",mat_c);waitKey();
}void interestValue(int *a,int &v)//得到5*5區域模板的興趣值
{int v1,v2,v3,v4;//v1是水平方向,v2是135度方向,v3是垂直方向,v4是45度方向v1=v2=v3=v4=0;v1=(a[11]-a[10])*(a[11]-a[10])+(a[12]-a[11])*(a[12]-a[11])+(a[13]-a[12])*(a[13]-a[12])+(a[14]-a[13])*(a[14]-a[13]);v2=(a[6]-a[0])*(a[6]-a[0])+(a[12]-a[6])*(a[12]-a[6])+(a[18]-a[12])*(a[18]-a[12])+(a[24]-a[18])*(a[24]-a[18]);v3=(a[7]-a[2])*(a[7]-a[2])+(a[12]-a[7])*(a[12]-a[7])+(a[17]-a[12])*(a[17]-a[12])+(a[22]-a[17])*(a[22]-a[17]);v4=(a[8]-a[4])*(a[8]-a[4])+(a[12]-a[8])*(a[12]-a[8])+(a[16]-a[12])*(a[16]-a[12])+(a[20]-a[16])*(a[20]-a[16]);v=v1;if(v>v2)v=v2;if(v>v3)v=v3;if(v>v4)v=v4;
}
void getmax(int *m_a,  int &m_max,int &m_num)      //得到7*7區域的候選特征值
{m_max=m_a[0];for(int i=1;i<49;i++){if(m_max<m_a[i]){m_max=m_a[i];m_num=i;}}}

?