11 #define ROI_X_FROM 261
12 #define ROI_X_TO (261+312)
14 #define ROI_Y_TO (42+246)
16 std::map<int, const char *> pix;
20 for (int i = FROM; i < TO; ++i) {
22 sprintf(filename, "%08d.ppm", i);
23 FILE *fp = fopen(filename, "rb");
29 fprintf(stderr, "Reading %s...\n", filename);
31 char *ptr = new char[WIDTH * HEIGHT * 3];
32 fseek(fp, HEADER, SEEK_SET);
34 if (fread(ptr, WIDTH * HEIGHT * 3, 1, fp) != 1) {
35 fprintf(stderr, "%s: Short read\n", filename);
42 for (int i = FROM; i < TO; ++i) {
43 for (int j = FROM; j < TO; ++j) {
44 const char *ptri = pix[i];
45 const char *ptrj = pix[j];
48 for (int y = ROI_Y_FROM; y < ROI_Y_TO; ++y) {
49 for (int x = ROI_X_FROM; x < ROI_X_TO; ++x) {
50 for (int c = 0; c < 3; ++c) {
51 sumdiff += abs(ptri[(y * WIDTH + x) * 3 + c] - ptrj[(y * WIDTH + x) * 3 + c]);
56 printf("%d %d %d\n", i, j, sumdiff);