相关：SIFT原理与源码解析
SURF原理与源码解析
在实时的视频流处理中，需要对每一帧特征提取，对算法处理速度上有很高的要求，传统的SIFT,Harris等特征点提取很难满足。由此提出Fast（Features from Accelerated Segment Test），由于不涉及尺度，梯度，等复杂运算，Fast检测器速度非常快。它使用一定邻域内像元的灰度值与中心点比较大小去判断是否为一个角点。但它的缺点是不具有方向性,尺度不变性。
转载请注明出处：
Fast角点提取步骤（以Fast-12-16为例）：
1.以固定半径为圆的边上取16个像素点（图中白色框出的位置），与中心点像素值Ip做差。
2.若边上存在连续的12（N>12,若为Fast-9,只需要N>9）个点满足 ( I(x)-I(p) )>threshold 或者( I(x)-I(p) ) < -threshold。(其中I(x)表示边上的像素值，I(p)为中心点像素值，threshold为设定的阈值。)则此点作为一个候选角点。如图上的虚线连接的位置。通常为了加速计算，直接比较1,5,9,13位置的差值，超过三个即视为一个候选点（存在连续的12个像元的必要条件），否则直接排除。
3.非极大值抑制，排除不稳定角点。采用强度响应函数：
即一个角点强度值定义为中心点与边上的12个角点像素差值的绝对值累加和。
opencv源码解析：
角点强度计算方法不采用上面的公式所描述，而是采用最小的差值(见代码分析)作为其角点强度值。
#include <stdio.h>#include <iostream>#include "cv.h"#include "opencv2/highgui/highgui.hpp"#include "opencv2/core/core.hpp"#include "opencv2/features2d/features2d.hpp"using namespace std;using namespace cv;int main( int argc, char** argv ){Mat img_1 = imread( "F:\\Picture\\hotel.jpg");if( !img_1.data ){ return -1; }FastFeatureDetector detector(50,true); //第一个参数为阈值，第二个采用非最大值抑制std::vector<KeyPoint> keypoints_1;detector.detect( img_1, keypoints_1 );//调用FAST_t函数检测，见下面详细解析drawKeypoints(img_1,keypoints_1,img_1,Scalar::all(255));imshow("HOTEL",img_1);waitKey(0);return 0;}void FAST_t(InputArray _img, std::vector<KeyPoint>& keypoints, int threshold, bool nonmax_suppression){Mat img = _img.getMat();const int K = patternSize/2, N = patternSize + K + 1;int i, j, k, pixel[25];makeOffsets(pixel, (int)img.step, patternSize);keypoints.clear();threshold = std::min(std::max(threshold, 0), 255);//保证阈值在0-255之间。uchar threshold_tab[512];for( i = -255; i <= 255; i++ )threshold_tab[i+255] = (uchar)(i < -threshold ? 1 : i > threshold ? 2 : 0); //分类成为darker、similar、brighter三种AutoBuffer<uchar> _buf((img.cols+16)*3*(sizeof(int) + sizeof(uchar)) + 128);uchar* buf[3];buf[0] = _buf; buf[1] = buf[0] + img.cols; buf[2] = buf[1] + img.cols;//保存对应角点强度值，否则为0int* cpbuf[3];cpbuf[0] = (int*)alignPtr(buf[2] + img.cols, sizeof(int)) + 1;//保存角点位置，+1为了存储这一行的角点总数cpbuf[1] = cpbuf[0] + img.cols + 1;cpbuf[2] = cpbuf[1] + img.cols + 1;memset(buf[0], 0, img.cols*3);for(i = 3; i < img.rows-2; i++){const uchar* ptr = img.ptr<uchar>(i) + 3;uchar* curr = buf[(i – 3)%3];int* cornerpos = cpbuf[(i – 3)%3];memset(curr, 0, img.cols);int ncorners = 0;if( i < img.rows – 3 ){j = 3;/*采用9点分割测试，加快检测速度检测任意一个直径两端的像素点，若同时与中心点相似，必定不是角点因为至少要占一半的数量*/for( ; j < img.cols – 3; j++, ptr++ ){int v = ptr[0];const uchar* tab = &threshold_tab[0] – v + 255;int d = tab[ptr[pixel[0]]] | tab[ptr[pixel[8]]];//if( d == 0 )// 加快检测速度[0]与[8]两个点都与中心点灰度值相近，，排除这个点continue;d &= tab[ptr[pixel[2]]] | tab[ptr[pixel[10]]];//直径两端两个点都相近，则为0d &= tab[ptr[pixel[4]]] | tab[ptr[pixel[12]]];//d &= tab[ptr[pixel[6]]] | tab[ptr[pixel[14]]];//每隔45度选取一个点if( d == 0 ) //continue;d &= tab[ptr[pixel[1]]] | tab[ptr[pixel[9]]];d &= tab[ptr[pixel[3]]] | tab[ptr[pixel[11]]];d &= tab[ptr[pixel[5]]] | tab[ptr[pixel[13]]];d &= tab[ptr[pixel[7]]] | tab[ptr[pixel[15]]];if( d & 1 ) //darker 中心值大，周围小的情况{int vt = v – threshold, count = 0;for( k = 0; k < N; k++ ) //且连续一半的像素点灰度差值( v-x > threshold )大于阈值{int x = ptr[pixel[k]];if(x < vt){if( ++count > K ){cornerpos[ncorners++] = j;if(nonmax_suppression)//非最大值抑制curr[j] = (uchar)cornerScore<patternSize>(ptr, pixel, threshold);//计算角点的强度响应值,最小的差值（绝对值）break;}}elsecount = 0;}}if( d & 2 )//brighter 中心值小，周围值大的情况{int vt = v + threshold, count = 0;for( k = 0; k < N; k++ ) //连续一半的像素点灰度差值( x-v < threshold )大于阈值{int x = ptr[pixel[k]];if(x > vt){if( ++count > K ){cornerpos[ncorners++] = j;if(nonmax_suppression)curr[j] = (uchar)cornerScore<patternSize>(ptr, pixel, threshold);//计算角点的强度响应值,最小的差值（绝对值）break;}}elsecount = 0;}}}}cornerpos[-1] = ncorners;//存储第i行上的角点总数量if( i == 3 )continue;/*与邻域的8个角点响应值做比较，非角点的响应值为0*/const uchar* prev = buf[(i – 4 + 3)%3]; //相邻的两行const uchar* pprev = buf[(i – 5 + 3)%3];//cornerpos = cpbuf[(i – 4 + 3)%3];//存储角点的列位置ncorners = cornerpos[-1]; //存储第i行上的角点总数量for( k = 0; k < ncorners; k++ ){j = cornerpos[k];int score = prev[j];if( !nonmax_suppression || //非极大值抑制,用角点强度值比较周围8个强度响应值(score > prev[j+1] && score > prev[j-1] &&score > pprev[j-1] && score > pprev[j] && score > pprev[j+1] &&score > curr[j-1] && score > curr[j] && score > curr[j+1]) ){keypoints.push_back(KeyPoint((float)j, (float)(i-1), 7.f, -1, (float)score));}}}}学习不是人生的全部，但学习都征服不了，那我还能做什么?
