Update for recent libpng
[svrt.git] / vision_problem_10.cc
1 /*
2  *  svrt is the ``Synthetic Visual Reasoning Test'', an image
3  *  generator for evaluating classification performance of machine
4  *  learning systems, humans and primates.
5  *
6  *  Copyright (c) 2009 Idiap Research Institute, http://www.idiap.ch/
7  *  Written by Francois Fleuret <francois.fleuret@idiap.ch>
8  *
9  *  This file is part of svrt.
10  *
11  *  svrt is free software: you can redistribute it and/or modify it
12  *  under the terms of the GNU General Public License version 3 as
13  *  published by the Free Software Foundation.
14  *
15  *  svrt is distributed in the hope that it will be useful, but
16  *  WITHOUT ANY WARRANTY; without even the implied warranty of
17  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
18  *  General Public License for more details.
19  *
20  *  You should have received a copy of the GNU General Public License
21  *  along with svrt.  If not, see <http://www.gnu.org/licenses/>.
22  *
23  */
24
25 #include "vision_problem_10.h"
26 #include "shape.h"
27
28 VisionProblem_10::VisionProblem_10() { }
29
30 void VisionProblem_10::generate(int label, Vignette *vignette) {
31   int nb_shapes = 4;
32   int xs[nb_shapes], ys[nb_shapes];
33
34   int error;
35   do {
36     do {
37       if(label == 1) {
38         scalar_t alpha = drand48() * 2 * M_PI;
39         scalar_t radius = drand48() * Vignette::width / 2;
40         scalar_t xc = drand48() * Vignette::width;
41         scalar_t yc = drand48() * Vignette::height;
42         for(int n = 0; n < nb_shapes; n++) {
43           xs[n] = int(xc + part_size/2 +
44                       radius * sin(alpha + n * 2 * M_PI / scalar_t(nb_shapes)));
45           ys[n] = int(yc + part_size/2 +
46                       radius * cos(alpha + n * 2 * M_PI / scalar_t(nb_shapes)));
47         }
48       } else {
49         for(int n = 0; n < nb_shapes; n++) {
50           xs[n] = int(drand48() * Vignette::width);
51           ys[n] = int(drand48() * Vignette::height);
52         }
53       }
54     } while(cluttered_shapes(part_size, nb_shapes, xs, ys));
55
56     vignette->clear();
57
58     Shape shape;
59     error = 0;
60     for(int n = 0; n < nb_shapes; n++) {
61       if(n == 0) {
62         shape.randomize(part_size / 2, part_hole_size / 2);
63       }
64       error |= shape.overwrites(vignette, xs[n], ys[n]);
65       shape.draw(vignette, xs[n], ys[n]);
66     }
67   } while(error);
68 }