svrt.c

   1
   2 /*
   3  *  svrt is the ``Synthetic Visual Reasoning Test'', an image
   4  *  generator for evaluating classification performance of machine
   5  *  learning systems, humans and primates.
   6  *
   7  *  Copyright (c) 2017 Idiap Research Institute, http://www.idiap.ch/
   8  *  Written by Francois Fleuret <francois.fleuret@idiap.ch>
   9  *
  10  *  This file is part of svrt.
  11  *
  12  *  svrt is free software: you can redistribute it and/or modify it
  13  *  under the terms of the GNU General Public License version 3 as
  14  *  published by the Free Software Foundation.
  15  *
  16  *  svrt is distributed in the hope that it will be useful, but
  17  *  WITHOUT ANY WARRANTY; without even the implied warranty of
  18  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19  *  General Public License for more details.
  20  *
  21  *  You should have received a copy of the GNU General Public License
  22  *  along with selector.  If not, see <http://www.gnu.org/licenses/>.
  23  *
  24  */
  25
  26 #include <TH/TH.h>
  27
  28 #include "svrt_generator.h"
  29
  30 THByteStorage *compress(THByteStorage *x) {
  31   long k, g, n;
  32
  33   k = 0; n = 0;
  34   while(k < x->size) {
  35     g = 0;
  36     while(k < x->size && x->data[k] == 255 && g < 255) { g++; k++; }
  37     n++;
  38     if(k < x->size && g < 255) { k++; }
  39   }
  40
  41   if(x->data[k-1] == 0) {
  42     n++;
  43   }
  44
  45   THByteStorage *result = THByteStorage_newWithSize(n);
  46
  47   k = 0; n = 0;
  48   while(k < x->size) {
  49     g = 0;
  50     while(k < x->size && x->data[k] == 255 && g < 255) { g++; k++; }
  51     result->data[n++] = g;
  52     if(k < x->size && g < 255) { k++; }
  53   }
  54   if(x->data[k-1] == 0) {
  55     result->data[n++] = 0;
  56   }
  57
  58   return result;
  59 }
  60
  61 THByteStorage *uncompress(THByteStorage *x) {
  62   long k, g, n;
  63
  64   k = 0;
  65   for(n = 0; n < x->size - 1; n++) {
  66     k = k + x->data[n];
  67     if(x->data[n] < 255) { k++; }
  68   }
  69   k = k + x->data[n];
  70
  71   THByteStorage *result = THByteStorage_newWithSize(k);
  72
  73   k = 0;
  74   for(n = 0; n < x->size - 1; n++) {
  75     for(g = 0; g < x->data[n]; g++) {
  76       result->data[k++] = 255;
  77     }
  78     if(x->data[n] < 255) {
  79       result->data[k++] = 0;
  80     }
  81   }
  82   for(g = 0; g < x->data[n]; g++) {
  83     result->data[k++] = 255;
  84   }
  85
  86   return result;
  87 }
  88
  89 THByteTensor *generate_vignettes(long n_problem, THLongTensor *labels) {
  90   struct VignetteSet vs;
  91   long nb_vignettes;
  92   long st0, st1, st2;
  93   long v, i, j;
  94   long *m, *l;
  95   unsigned char *a, *b;
  96
  97   if(THLongTensor_nDimension(labels) != 1) {
  98     printf("Label tensor has to be of dimension 1.\n");
  99     exit(1);
 100   }
 101
 102   nb_vignettes = THLongTensor_size(labels, 0);
 103   m = THLongTensor_storage(labels)->data + THLongTensor_storageOffset(labels);
 104   st0 = THLongTensor_stride(labels, 0);
 105   l = (long *) malloc(sizeof(long) * nb_vignettes);
 106   for(v = 0; v < nb_vignettes; v++) {
 107     l[v] = *m;
 108     m += st0;
 109   }
 110
 111   svrt_generate_vignettes(n_problem, nb_vignettes, l, &vs);
 112   free(l);
 113
 114   THLongStorage *size = THLongStorage_newWithSize(3);
 115   size->data[0] = vs.nb_vignettes;
 116   size->data[1] = vs.height;
 117   size->data[2] = vs.width;
 118
 119   THByteTensor *result = THByteTensor_newWithSize(size, NULL);
 120   THLongStorage_free(size);
 121
 122   st0 = THByteTensor_stride(result, 0);
 123   st1 = THByteTensor_stride(result, 1);
 124   st2 = THByteTensor_stride(result, 2);
 125
 126   unsigned char *r = vs.data;
 127   for(v = 0; v < vs.nb_vignettes; v++) {
 128     a = THByteTensor_storage(result)->data + THByteTensor_storageOffset(result) + v * st0;
 129     for(i = 0; i < vs.height; i++) {
 130       b = a + i * st1;
 131       for(j = 0; j < vs.width; j++) {
 132         *b = (unsigned char) (*r);
 133         r++;
 134         b += st2;
 135       }
 136     }
 137   }
 138
 139   free(vs.data);
 140
 141   return result;
 142 }