-
Notifications
You must be signed in to change notification settings - Fork 0
/
moravec.c
47 lines (37 loc) · 923 Bytes
/
moravec.c
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
#include<stdio.h>
#define min (a, b, c, d)
int main() {
Mat img = read
int S01 = 0;
int S0_1 = 0;
int S10 = 0;
int S_10 = 0;
int x;
int y;
int C = 0;
//S01 = 4
for (y = 4; y <= 6; y++) {
for (x = 2; x <= 4; x++) {
S01 += ((img.at<uchar>(y, x+1) - img.ar<uchar>(y, x)) * (img.at<uchar>(y, x+1) - img.ar<uchar>(y, x)));
}
}
//S0-1 = 3
for (y = 4; y <= 6; y++) {
for (x = 2; x <= 4; x++) {
S0_1 += ((img.ar<uchar>(y, x-1) - img.at<uchar>(y, x)) * (img.ar<uchar>(y, x-1) - img.at<uchar>(y, x)));
}
}
//S10
for (y = 4; y <= 6; y++) {
for (x = 2; x <= 4; x++) {
S10 += ((img.ar<uchar>(y+1, x) - img.at<uchar>(y, x)) * (img.ar<uchar>(y+1, x) - img.at<uchar>(y, x)));
}
}
//S-10
for (y = 4; y <= 6; y++) {
for (x = 2; x <= 4; x++) {
S10 += ((img.ar<uchar>(y-1, x) - img.at<uchar>(y, x)) * (img.ar<uchar>(y-1, x) - img.at<uchar>(y, x)));
}
}
C = min(S01, S0_1, S10, S_10);
}