sparrow.c

   1 /* Copyright (C) <2010> Douglas Bagnall <douglas@halo.gen.nz>
   2  *
   3  * This library is free software; you can redistribute it and/or
   4  * modify it under the terms of the GNU Library General Public
   5  * License as published by the Free Software Foundation; either
   6  * version 2 of the License, or (at your option) any later version.
   7  *
   8  * This library is distributed in the hope that it will be useful,
   9  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11  * Library General Public License for more details.
  12  *
  13  * You should have received a copy of the GNU Library General Public
  14  * License along with this library; if not, write to the
  15  * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  16  * Boston, MA 02111-1307, USA.
  17  */
  18
  19 #include "sparrow.h"
  20 #include "gstsparrow.h"
  21
  22 #include <string.h>
  23 #include <math.h>
  24
  25 /* static functions (via `make cproto`) */
  26 static void change_state(GstSparrow *sparrow, sparrow_state state);
  27
  28
  29 /* set up whatever debugging methods are enabled */
  30 static void
  31 init_debug(GstSparrow *sparrow){
  32   if (!sparrow->debug_frame){
  33     sparrow->debug_frame = malloc_aligned_or_die(MAX(sparrow->in.size, sparrow->out.size));
  34   }
  35 }
  36
  37 /*RNG code */
  38 /*seed with -1 for automatic seed choice */
  39 static void rng_init(GstSparrow *sparrow, guint32 seed){
  40   GST_DEBUG("in RNG init\n");
  41   if (seed == (guint32)-1){
  42     /* XXX should really use /dev/urandom */
  43     seed = rand() ^ (int)&(sparrow);
  44     GST_DEBUG("Real seed %u\n", seed);
  45   }
  46   if (seed == 0)
  47     seed = 12345;
  48
  49   dsfmt_init_gen_rand(sparrow->dsfmt, seed);
  50   dsfmt_gen_rand_all(sparrow->dsfmt);
  51   GST_DEBUG("RNG seeded with %u\n", seed);
  52 }
  53
  54 /** debugging: write frames out somewhere. **/
  55
  56 /*spit out the frame as a ppm image */
  57 static void
  58 ppm_dump(sparrow_format *rgb, guint8 *data, guint32 width, guint32 height, char *name)
  59 {
  60   guint i;
  61   FILE *fh = fopen(name, "w");
  62   guint32 size = width * height;
  63   fprintf(fh, "P6\n%u %u\n255\n", width, height);
  64   /* 4 cases: xBGR xRGB BGRx RGBx
  65      need to convert to 24bit R G B
  66      XX maybe could optimise some cases?
  67   */
  68   guint32 *p = (guint32 *)data;
  69   for (i = 0; i < size; i++){
  70     putc((*p >> rgb->rshift) & 255, fh);
  71     putc((*p >> rgb->gshift) & 255, fh);
  72     putc((*p >> rgb->bshift) & 255, fh);
  73     p++;
  74   }
  75   fflush(fh);
  76   fclose(fh);
  77 }
  78
  79 /*pgm for greyscale */
  80 static void
  81 pgm_dump(guint8 *data, guint32 width, guint32 height, char *name)
  82 {
  83   FILE *fh = fopen(name, "w");
  84   size_t size = width * height;
  85   fprintf(fh, "P5\n%u %u\n255\n", width, height);
  86   size_t wrote = fwrite(data, 1, size, fh);
  87   if (wrote != size){
  88     GST_DEBUG("wanted to write %u bytes; fwrite said %u\n", size, wrote);
  89   }
  90   fflush(fh);
  91   fclose(fh);
  92 }
  93
  94
  95 #define PPM_FILENAME_TEMPLATE "/tmp/sparrow_%05d.%s"
  96 #define PPM_FILENAME_LENGTH (sizeof(PPM_FILENAME_TEMPLATE) + 15)
  97
  98 void INVISIBLE
  99 debug_frame(GstSparrow *sparrow, guint8 *data, guint32 width, guint32 height, int pixsize){
 100   char name[PPM_FILENAME_LENGTH];
 101   int res;
 102   if (pixsize == 4){ /*rgb*/
 103     res = snprintf(name, PPM_FILENAME_LENGTH, PPM_FILENAME_TEMPLATE, sparrow->frame_count, "ppm");
 104     if (res > 0){
 105       ppm_dump(&(sparrow->in), data, width, height, name);
 106     }
 107   }
 108   else {/*greyscale*/
 109     res = snprintf(name, PPM_FILENAME_LENGTH, PPM_FILENAME_TEMPLATE, sparrow->frame_count, "pgm");
 110     if (res > 0){
 111       pgm_dump(data, width, height, name);
 112     }
 113   }
 114 }
 115
 116 /** interpret gst attributes **/
 117
 118 /* Extract a colour (R,G,B) bitmask from gobject  */
 119 static guint32 get_mask(GstStructure *s, char *mask_name){
 120   gint32 mask;
 121   int res = gst_structure_get_int(s, mask_name, &mask);
 122   if (!res){
 123     GST_WARNING("No mask for '%s' !\n", mask_name);
 124   }
 125   return (guint32)mask;
 126 }
 127
 128 static void
 129 extract_caps(sparrow_format *im, GstCaps *caps)
 130 {
 131   GstStructure *s = gst_caps_get_structure (caps, 0);
 132   gst_structure_get_int(s, "width", &(im->width));
 133   gst_structure_get_int(s, "height", &(im->height));
 134   im->rshift = mask_to_shift(get_mask(s, "red_mask"));
 135   im->gshift = mask_to_shift(get_mask(s, "green_mask"));
 136   im->bshift = mask_to_shift(get_mask(s, "blue_mask"));
 137   /* recalculate shifts as little-endian */
 138   im->rmask = 0xff << im->rshift;
 139   im->gmask = 0xff << im->gshift;
 140   im->bmask = 0xff << im->bshift;
 141
 142   im->pixcount = im->width * im->height;
 143   im->size = im->pixcount * PIXSIZE;
 144   im->colours[SPARROW_WHITE] = im->rmask | im->gmask | im->bmask;
 145   im->colours[SPARROW_GREEN] = im->gmask;
 146   im->colours[SPARROW_MAGENTA] = im->rmask | im->bmask;
 147
 148   GST_DEBUG("\ncaps:\n%" GST_PTR_FORMAT, caps);
 149   GST_DEBUG("shifts: r %u g %u b %u\n", im->rshift, im->gshift, im->bshift);
 150   GST_DEBUG("dimensions: w %u h %u pix %u size %u\n", im->width, im->height,
 151       im->pixcount, im->size);
 152 }
 153
 154
 155
 156 /*Most functions below here are called from gstsparrow.c and are NOT static */
 157
 158 void INVISIBLE
 159 sparrow_rotate_history(GstSparrow *sparrow, GstBuffer *inbuf){
 160   if (sparrow->in_buffer){
 161     gst_buffer_unref(sparrow->prev_buffer);
 162     sparrow->prev_buffer = sparrow->in_buffer;
 163     sparrow->prev_frame = sparrow->in_frame;
 164   }
 165   gst_buffer_ref(inbuf);
 166   sparrow->in_buffer = inbuf;
 167
 168   sparrow->in_frame = GST_BUFFER_DATA(inbuf);
 169 }
 170
 171 /* called by gst_sparrow_init(). The source/sink capabilities (and commandline
 172    arguments[?]) are unknown at this stage, so there isn't much useful to do
 173    here.*/
 174 void INVISIBLE
 175 sparrow_pre_init(GstSparrow *sparrow){
 176 }
 177
 178 /* called by gst_sparrow_set_caps(). This sets up everything after gstreamer
 179    has worked out what the pipeline will look like.
 180  */
 181 gboolean INVISIBLE
 182 sparrow_init(GstSparrow *sparrow, GstCaps *incaps, GstCaps *outcaps){
 183   change_state(sparrow, SPARROW_INIT);
 184   extract_caps(&(sparrow->in), incaps);
 185   extract_caps(&(sparrow->out), outcaps);
 186   sparrow_format *in = &(sparrow->in);
 187
 188   sparrow->work_frame = zalloc_aligned_or_die(in->size);
 189   sparrow->dsfmt = zalloc_aligned_or_die(sizeof(dsfmt_t));
 190   sparrow->screenmask = malloc_aligned_or_die(in->pixcount);
 191
 192 #if USE_SPARSE_MAP
 193   size_t point_memsize = (sizeof(sparrow_map_point_t) * sparrow->out.pixcount / LINE_PERIOD) + 1;
 194   size_t row_memsize = sizeof(sparrow_map_row_t) * sparrow->out.height + 1;
 195   sparrow->map.point_mem = malloc_aligned_or_die(point_memsize);
 196   sparrow->map.rows = zalloc_aligned_or_die(row_memsize);
 197 #else
 198   size_t lutsize = sizeof(sparrow_map_lut_t) * sparrow->out.pixcount;
 199   sparrow->map_lut = malloc_aligned_or_die(lutsize);
 200 #endif
 201
 202   sparrow->prev_buffer = gst_buffer_new_and_alloc(in->size);
 203   sparrow->prev_frame  = GST_BUFFER_DATA(sparrow->prev_buffer);
 204   memset(sparrow->prev_frame, 0, in->size);
 205
 206   sparrow->timer_start.tv_sec = 0;
 207   sparrow->timer_stop.tv_sec = 0;
 208
 209
 210   rng_init(sparrow, sparrow->rng_seed);
 211
 212   if (sparrow->debug){
 213     init_debug(sparrow);
 214   }
 215
 216   if (sparrow->rng_seed & 2){/*XXX need better test */
 217     sparrow->colour = sparrow->out.colours[SPARROW_GREEN];
 218   }
 219   else {
 220     sparrow->colour = sparrow->out.colours[SPARROW_MAGENTA];
 221   }
 222
 223   sparrow->timer_log = (sparrow->use_timer) ? fopen(TIMER_LOG_FILE, "w") : NULL;
 224   change_state(sparrow, SPARROW_NEXT_STATE);
 225   return TRUE;
 226 }
 227
 228 void INVISIBLE
 229 sparrow_finalise(GstSparrow *sparrow)
 230 {
 231   free(sparrow->work_frame);
 232   free(sparrow->dsfmt);
 233   free(sparrow->screenmask);
 234 #if USE_SPARSE_MAP
 235   free(sparrow->map.point_mem);
 236   free(sparrow->map.rows);
 237 #else
 238   free(sparrow->map_lut);
 239 #endif
 240
 241
 242   if (sparrow->timer_log){
 243     fclose(sparrow->timer_log);
 244   }
 245   //free everything
 246   //cvReleaseImageHeader(IplImage** image)
 247 }
 248
 249
 250 /* initialisation functions and sparrow_transform() use this to set up a new
 251    state. */
 252 static void
 253 change_state(GstSparrow *sparrow, sparrow_state state)
 254 {
 255   GST_DEBUG("state is %d, sparrow->state is %d\n", state, sparrow->state);
 256   switch(sparrow->state){
 257   case SPARROW_FIND_SELF:
 258     finalise_find_self(sparrow);
 259     break;
 260   case SPARROW_FIND_SCREEN:
 261     finalise_find_screen(sparrow);
 262     break;
 263   case SPARROW_FIND_EDGES:
 264     finalise_find_edges(sparrow);
 265     break;
 266   case SPARROW_PLAY:
 267     finalise_process_frame(sparrow);
 268     break;
 269   case SPARROW_INIT:
 270     break;
 271   default:
 272     GST_DEBUG("change_state got unknown state: %d\n", state);
 273   }
 274   if (state == SPARROW_NEXT_STATE){
 275     state = sparrow->state + 1;
 276   }
 277   switch(state){
 278   case SPARROW_FIND_SELF:
 279     init_find_self(sparrow);
 280     break;
 281   case SPARROW_FIND_SCREEN:
 282     init_find_screen(sparrow);
 283     break;
 284   case SPARROW_FIND_EDGES:
 285     init_find_edges(sparrow);
 286     break;
 287   case SPARROW_INIT:
 288   case SPARROW_PLAY:
 289     break;
 290   default:
 291     GST_DEBUG("change_state got unknown state: %d\n", state);
 292   }
 293   sparrow->state = state;
 294 }
 295
 296
 297 /*called by gst_sparrow_transform_ip every frame.
 298
 299   decide what to do based on sparrow->state. All the processing is done in a
 300   "mode_*" function, which returns a state or SPARROW_STATUS_QUO.  If a state
 301   is returned, then change_state() is called to initialise the state, even if
 302   it is the current state (so states can use this to reset).
 303 */
 304 void INVISIBLE
 305 sparrow_transform(GstSparrow *sparrow, guint8 *in, guint8 *out)
 306 {
 307   sparrow_state new_state;
 308 #if TIME_TRANSFORM
 309   TIMER_START(sparrow);
 310 #endif
 311   switch(sparrow->state){
 312   case SPARROW_FIND_SELF:
 313     new_state = mode_find_self(sparrow, in, out);
 314     break;
 315   case SPARROW_FIND_SCREEN:
 316     new_state = mode_find_screen(sparrow, in, out);
 317     break;
 318   case SPARROW_FIND_EDGES:
 319     new_state = mode_find_edges(sparrow, in, out);
 320     break;
 321   default:
 322     new_state = mode_process_frame(sparrow, in, out);
 323   }
 324   sparrow->frame_count++;
 325   if (new_state != SPARROW_STATUS_QUO){
 326     change_state(sparrow, new_state);
 327   }
 328 #if TIME_TRANSFORM
 329   TIMER_STOP(sparrow);
 330 #endif
 331 }
 332