FaceDetection/Feature.cpp at master · limitforever/FaceDetection · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143

/*2016.4.13
特征选取和计算函数
给定图像矩阵，长宽以及坐标，利用积分图法快速计算给定5种特征的特征值*/

#include <string>
#include <iostream>
#include "Feature.h"


Feature::Feature(int t, int x, int y, int w, int h) {    //x,y表示位置，w,h表示大小，t表示类型
	type = t;
	xc = x;
	yc = y;
	width = w;
	height = h;
}

int Feature::getType() {
	return this->type;
}

int Feature::getWidth() {
	return this->width;
}

int Feature::getHeight() {
	return this->height;
}

int Feature::getxc() {
	return this->xc;
}

int Feature::getyc() {
	return this->yc;
}

/*利用积分图计算矩形特征值*/
float Feature::rectangleValue(float *ii, int iiwidth, int ix, int iy, int rx, int ry, int rw, int rh) {
	float value = ii[((iy + ry + rh - 1)*iiwidth) + (ix + rx + rw - 1)];         //A+B+C+D
	if ((ix + rx) > 0) value -= ii[((iy + ry + rh - 1)*iiwidth) + (ix + rx - 1)];  //B+D(x>0)
	if ((iy + ry) > 0) value -= ii[((iy + ry - 1)*iiwidth) + (ix + rx + rw - 1)];  //C+D(y>0)
	if ((ix + rx) > 0 && (iy + ry) > 0) value += ii[((iy + ry - 1)*iiwidth) + (ix + rx - 1)];   //D(x>0且y>0)
	return value;
}
/*设定特征类型*/
float Feature::getValue(float *ii, int iiwidth, int x, int y) {
	// 5 types of feature (A, B, C, C^t, D)
	float t = 0;
	switch (this->type) {
		//switch (type) {  //xc,yc为积分图起始位置的坐标（左上角）
	case 0:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc + (this->width / 2), this->yc, this->width / 2, this->height) - rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width / 2, this->height));
		//std::cout << "特征1的特征值为：" << t << std::endl;
		return t;
		break;
	case 1:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width, this->height / 2) - rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height / 2), this->width, this->height / 2));
		//std::cout << "特征2的特征值为：" << t << std::endl;
		return t;
		break;
	case 2:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width / 3, this->height) + rectangleValue(ii, iiwidth, x, y, this->xc + (width * 2 / 3), this->yc, this->width / 3, this->height) - 2 * rectangleValue(ii, iiwidth, x, y, this->xc + (width / 3), this->yc, this->width / 3, this->height)) / 2;
		//std::cout << "特征3的特征值为：" << t << std::endl;
		return t;
		break;
	case 3:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width, this->height / 3) + rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height * 2 / 3), this->width, this->height / 3) - 2 * rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height / 3), this->width, this->height / 3)) / 2;
		//std::cout << "特征4的特征值为：" << t << std::endl;
		return t;
		break;
	case 4:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width * 2 / 5, this->height) + rectangleValue(ii, iiwidth, x, y, this->xc + (width * 3 / 5), this->yc, this->width * 2 / 5, this->height) - 4 * rectangleValue(ii, iiwidth, x, y, this->xc + (width * 2 / 5), this->yc, this->width / 5, this->height)) / 4;
		//std::cout << "特征3的特征值为：" << t << std::endl;
		return t;
		break;
	case 5:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width, this->height * 2 / 5) - rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height * 3 / 5), this->width, this->height * 2 / 5) - 4 * rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height * 2 / 5), this->width, this->height / 5)) / 4;
		//std::cout << "特征4的特征值为：" << t << std::endl;
		return t;
		break;

	case 6:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width / 5, this->height) + rectangleValue(ii, iiwidth, x, y, this->xc + (width * 4 / 5), this->yc, this->width / 5, this->height) - 2 / 3 * rectangleValue(ii, iiwidth, x, y, this->xc + (width / 5), this->yc, this->width * 3 / 5, this->height))*1.5;
		//std::cout << "特征3的特征值为：" << t << std::endl;
		return t;
		break;
	case 7:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width, this->height / 5) + rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height * 4 / 5), this->width, this->height / 5) - 2 / 3 * rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (this->height / 5), this->width, this->height * 3 / 5))*1.5;
		//std::cout << "特征4的特征值为：" << t << std::endl;
		return t;
		break;

	case 8:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc + (width / 2), this->yc, this->width / 2, this->height / 2) + rectangleValue(ii, iiwidth, x, y, this->xc, this->yc + (height / 2), this->width / 2, this->height / 2) - rectangleValue(ii, iiwidth, x, y, this->xc + (width / 2), this->yc + (height / 2), this->width / 2, this->height / 2) - rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width / 2, this->height / 2)) / 2;
		//std::cout << "特征5的特征值为：" << t << std::endl;
		return t;
		break;
	case 9:
		t = (rectangleValue(ii, iiwidth, x, y, this->xc, this->yc, this->width, this->height) - 9 * rectangleValue(ii, iiwidth, x, y, this->xc + width / 3, this->yc + (height / 3), this->width / 3, this->height / 3)) / 8;
		//std::cout << "特征5的特征值为：" << t << std::endl;
		return t;
		break;
	default:
		std::cout << "错误的特征：" << this->type << " 并不存在！" << std::endl;
		exit(-1);
		break;
	}

}
/*特征比例*/
void Feature::scale(float s) {
	this->width = (this->width)*s;
	this->height = (this->height)*s;
	this->xc = (this->xc)*s;
	this->yc = (this->yc)*s;
}

/*备用快速积分法*/
void fastIntegral(unsigned char* inputMatrix, unsigned long* outputMatrix, int width, int height){
	unsigned long *columnSum = new unsigned long[width]; // sum of each column
	// calculate integral of the first line
	for (int i = 0; i<width; i++){
		columnSum[i] = inputMatrix[i];
		outputMatrix[i] = inputMatrix[i];
		if (i>0){
			outputMatrix[i] += outputMatrix[i - 1];
		}
	}
	for (int i = 1; i<height; i++){
		int offset = i*width;
		// first column of each line
		columnSum[0] += inputMatrix[offset];
		outputMatrix[offset] = columnSum[0];
		// other columns
		for (int j = 1; j<width; j++){
			columnSum[j] += inputMatrix[offset + j];
			outputMatrix[offset + j] = outputMatrix[offset + j - 1] + columnSum[j];
		}
	}
	return;
}