2 ///////////////////////////////////////////////////////////////////////////
3 // This program is free software: you can redistribute it and/or modify //
4 // it under the terms of the version 3 of the GNU General Public License //
5 // as published by the Free Software Foundation. //
7 // This program is distributed in the hope that it will be useful, but //
8 // WITHOUT ANY WARRANTY; without even the implied warranty of //
9 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU //
10 // General Public License for more details. //
12 // You should have received a copy of the GNU General Public License //
13 // along with this program. If not, see <http://www.gnu.org/licenses/>. //
15 // Written by Francois Fleuret //
16 // (C) Idiap Research Institute //
18 // Contact <francois.fleuret@idiap.ch> for comments & bug reports //
19 ///////////////////////////////////////////////////////////////////////////
23 A RichImage is an Image with all the necessary pre-computation to
24 compute the feature responses: Edge detection and local variance
25 thresholding at multiple scale, with integral images. The heavy
26 machinery and ugly coding style is motivated by performance.
34 #include "rgb_image.h"
38 class RichImage : public Image {
40 static const int _image_size_min = 8;
44 int *_width_at_scale, *_height_at_scale;
46 unsigned int *_edge_map;
48 unsigned int **_line_edge_map;
49 unsigned int ***_tag_edge_map;
50 unsigned int ****_scale_edge_map;
54 void compute_one_scale_edge_maps(int width, int height,
55 unsigned int ***scale_edge_map,
56 unsigned char *pixel_map,
57 unsigned int *sum_pixel_map,
58 unsigned int *sum_sq_pixel_map);
62 // We have 8 edge orientations
63 static const int nb_edge_tags = 8;
64 static const int first_edge_tag = 0;
66 // The variance tag is 1 if the variance of the gray levels in the
67 // 16x16 square is greater than 10
68 static const int variance_tag = 8;
70 // We quantify the grayscales into 8 bins
71 static const int nb_gray_tags = 8;
72 static const int first_gray_tag = 9;
74 static const int _nb_tags = nb_edge_tags + 1 + nb_gray_tags;
76 inline int nb_tags() { return _nb_tags; }
77 inline int nb_scales() { return _nb_scales; }
79 inline int nb_tags_in_window(int scale, int tag, int xmin, int ymin, int xmax, int ymax) {
80 if(scale < 0 || scale >= _nb_scales) return 0;
81 if(xmin < 0) xmin = 0;
82 if(xmax >= _width_at_scale[scale]) xmax = _width_at_scale[scale] - 1;
83 if(ymin < 0) ymin = 0;
84 if(ymax >= _height_at_scale[scale]) ymax = _height_at_scale[scale] - 1;
85 if(xmin < xmax && ymin < ymax) {
86 unsigned int **tmp = _scale_edge_map[scale][tag];
87 return tmp[ymax][xmax] + tmp[ymin][xmin] - tmp[ymax][xmin] - tmp[ymin][xmax];
93 RichImage(int width, int height);
96 virtual void compute_rich_structure();
98 virtual void write(ostream *out);
100 // read does _NOT_ build the rich structure. One has to call
101 // compute_rich_structure() for that!
102 virtual void read(istream *in);
104 virtual void write_tag_png(const char *filename);