基于OpenCV的LBP算法(OpenCV1.0版本)

#include 
#include 
#include 

using namespace std;
using namespace cv;
//基于旧版本的opencv的LBP算法opencv1.0  
// 3 x 3 矩阵如下所示
// [ 1, 2, 3]
// [ 8, ij,4]
// [ 7, 6, 5]
void LBP(IplImage *src, IplImage *dst)
{
	// 进行处理的原图像为单通道图像
	assert(src != NULL && src->nChannels == 1);

	int tmp[8] = { 0 };

	int rows = src->height - 1;
	int cols = src->width - 1;
	
	for (int i = 1; i < rows; ++i)
	{
		for (int j = 1; j < cols; ++j)
		{
			int sum = 0;
			double val = cvGetReal2D(src, i, j);
			double tempVal = 0.0;

			tempVal = cvGetReal2D(src, i - 1, j - 1);		// 左上角
			tempVal > val ? tmp[0] = 255 : tmp[0] = 0;
			tempVal = cvGetReal2D(src, i, j - 1);			// 正上方
			tempVal > val ? tmp[1] = 255 : tmp[1] = 0;
			tempVal = cvGetReal2D(src, i + 1, j - 1);		// 右上角
			tempVal > val ? tmp[2] = 255 : tmp[2] = 0;
			tempVal = cvGetReal2D(src, i + 1, j);			// 右侧
			tempVal > val ? tmp[3] = 255 : tmp[3] = 0;		
			tempVal = cvGetReal2D(src, i + 1, j + 1);		// 右下角
			tempVal > val ? tmp[4] = 255 : tmp[4] = 0;
			tempVal = cvGetReal2D(src, i, j + 1);			// 正下方
			tempVal > val ? tmp[5] = 255 : tmp[5] = 0;
			tempVal = cvGetReal2D(src, i - 1, j + 1);		// 左下角
			tempVal > val ? tmp[6] = 255 : tmp[6] = 0;		
			tempVal = cvGetReal2D(src, i - 1, j);			// 左侧
			tempVal > val ? tmp[7] = 255 : tmp[7] = 0;

			// 计算 LBP 编码
			for (int k = 0; k < 8; ++k)
			{
				sum += tmp[k] * pow(2, k);
			}

			cvSetReal2D(dst, i, j, sum);
		}
	}
}

// 圆形 LBP 算子
void ELBP(IplImage* src, IplImage* dst, int radius, int neighbors)
{
	// 处理的图像为单通道图像
	assert(src->nChannels == 1);

	for (int i = 0; i < neighbors; ++i)
	{
		// 正弦弧度
		double sRadian = sin(2.0 * CV_PI * i / static_cast(neighbors));
		// 余弦弧度
		double cRadian = cos(2.0 * CV_PI * i / static_cast(neighbors));

		// 采样点的计算
		double x = static_cast(-radius * sRadian);
		double y = static_cast(radius * cRadian);

		// 下取整的值
		int fx = static_cast(floor(x));
		int fy = static_cast(floor(y));
		// 上取整的值
		int cx = static_cast(ceil(x));
		int cy = static_cast(ceil(y));
		// 小数部分
		double tx = x - fx;
		double ty = y - fy;

		// 设置插值权重
		double w1 = (1 - tx) * (1 - ty);
		double w2 = tx * (1 - ty);
		double w3 = (1 - tx) * ty;
		double w4 = tx * ty;

	
		// 循环处理图像数据
		for (int rows = radius; rows < src->height - radius; ++rows)
		{
			for (int cols = radius; cols < src->width - radius; ++cols)
			{
				// 计算插值
				double t1 = w1 * cvGetReal2D(src, rows + fy, cols + fx);
				double t2 = w2 * cvGetReal2D(src, rows + fy, cols + cx);
				double t3 = w3 * cvGetReal2D(src, rows + cy, cols + fx);
				double t4 = w4 * cvGetReal2D(src, rows + cy, cols + cx);
				double t = t1 + t2 + t3 + t4;

				double val = cvGetReal2D(src, rows, cols);
				double epsilon = std::numeric_limits::epsilon();

				uchar c = ((t > val) || abs(t - val) < epsilon);
				uchar tmp = c * pow(2, i);
				double v = cvGetReal2D(dst, rows - radius, cols - radius);
				v += tmp;
				cvSetReal2D(dst, rows - radius, cols - radius, v);
			}
		}

	}

}