/*  
 *  Data file from http://yann.lecun.com/exdb/mnist/
 *  2006.10.26 A.Date 
 */

#include <stdio.h>
#include <stdlib.h>
#include <fcntl.h>

#define IMAGE_FILE "train-images-idx3-ubyte"
#define NUM_IMAGES 60000
#define SIZE 784   /* 28 x 28 */

static int num[10];
static unsigned char image[NUM_IMAGES][SIZE];

/* * Read the training images into memory */
void read_images() {

  int           i, j,k, fd;
  char          image_name[80];
  float f;
  unsigned char * ptr; 
  
  if ((fd=open(IMAGE_FILE,O_RDONLY))==-1){
		  printf("couldn't open image file");
		  exit(0);
  }
  
  /* skip headers */
  read(fd,num, 4*sizeof(int) ); 
  
  for (i=0; i<NUM_IMAGES; i++) { 
      read(fd,image[i], SIZE*sizeof(unsigned char) ); 
  }
  
  for (i=0; i<NUM_IMAGES; i++) { 
      printf("#####   image %d/%d\n",i, NUM_IMAGES); 
      for (j=0; j<SIZE; j++) { 
          printf("%3d", image[i][j]); 
          if ( (j+1) % 28 == 0 ){
              printf("\n"); 
          }
      }
      printf("\n\n"); 
  }
  close(fd);
  
}



int main(int argc, char **argv){

	read_images();

	return 0;
}