Forstner算子提取

最新推荐文章于 2023-08-09 10:53:43 发布

OrangeMo0n

最新推荐文章于 2023-08-09 10:53:43 发布

阅读量3.6k

点赞数

分类专栏：摄影测量文章标签：摄影测量 C++ Forstner

本文链接：https://blog.csdn.net/OrangeMo0n/article/details/9149011

版权

本文介绍了使用Forstner算子进行兴趣点检测的过程，包括计算方差确定阈值，计算像素差值，窗口内的圆度和兴趣值评估，以及极值点判断。通过设定阈值和比较pw值，筛选出高质量的兴趣点。

摘要由CSDN通过智能技术生成

/**
 * @brief 提取Forstner算子
 * @param pBuffer 存放数据的指针
 * @param iCol,iRow,iWidth,iHeight 影像区块的起始行列号、宽高
 * @param vecForstnerPnts Forstner算子
 * @param iPointCount 提取的算子点数目
 */
bool Forstner_Operator(float *pBuffer, 
	int iCol, int iRow, int iWidth, int iHeight,
	vector<STPoint> &vecForstnerPnts, int iPointCount);

（1）计算数据的方差是否大于一定阈值（3.0），否则不存在兴趣点计算像素的和dSumPixel，像素的平方和dSquareSumPixel

方差dVar计算公式：Var(X)=E(X^2)-(E(X)^2)

float fVarThres = 3.0;
float fSum = 0.0; 
float fSquareSum = 0.0;
int iBufferSize = iWidth*iHeight;
for (int i=0; i<iBufferSize; i++)
{
	fSum += pBuffer[i];
	fSquareSum += pBuffer[i]*pBuffer[i];
}
float fVar = (fSquareSum - (fSum*fSum)/iBufferSize)/iBufferSize;
if (fVar < fVarThres)
{
	return false;
}

（2