#include<errno.h>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#include "neuzeichnen.h"
#include "mygraph.h"
#include "ReadData.h"

void usage(const char *message){
  printf("%s\n",message);
  printf("Usage: Contor_mpg [-magnify newx newy xofs yofs]\n"
	 "                  [-center]\n"
         "                  [-nocuts nocuts]\n"
         "                  [-noadjust min max cut1 ... cutn]\n"
	 "                  filename1 ... filenameN\n"
	 "where the datatypes are given by:\n"
	 "     int newx newy\n"
	 "     double xofs yofs min max cut1 ... cutN\n"
	 "-nocuts needs to be defined before -noadjust\n"
	 "If the filenames should be obtained automatically there is also the\n"
	 "option '-name formatstring increment LargestNumber' which will \n"
	 "automatically write the appropriate filename\n"
	 "(c) Alexander Wagner 2003 (requires mjpegtools)\n");
  exit(1);
}

int main(int argc, char *argv[])
{
  int number=200,i,j,k,Xsize,Ysize,dummy,
    compress=0,magnify=0,center=0,newx,newy,
    nocuts=1,adjust;
  char **names=NULL,*spointer,*command,addstr[100],
    ppmstr[500];
  int ppmx,ppmy;
  double *datap=NULL,*cut=NULL,dmin,dmax,xofs,yofs;
  FILE *ppmfile;


  if (argc>1) { /* filenames are given as parameters to the program */
    i=1;
    errno=0;
    while ((i<argc)&&(argv[i][0]=='-')){
      if (0==strcmp(argv[i],"-magnify")){
	newx=strtod(argv[i+1],&spointer);
	newy=strtod(argv[i+2],&spointer);
	xofs=strtod(argv[i+3],&spointer);
	yofs=strtod(argv[i+4],&spointer);
	magnify=1;
	i+=5;
      } 
      else if (0==strcmp(argv[i],"-center")){
	center=0;
	i++;
      }
      else if (0==strcmp(argv[i],"-nocuts")){
	nocuts=strtod(argv[i+1],&spointer);
	cut=(double *) malloc(nocuts*sizeof(double));
	i+=2;
      }
      else if (0==strcmp(argv[i],"-noadjust")){
	adjust=1;
	dmin=strtod(argv[i+1],&spointer);
	dmax=strtod(argv[i+2],&spointer);
	i+=3;
	for (j=0;j<nocuts;j++){
	  cut[j]=strtod(argv[i+j],&spointer);
	}
	i+=nocuts;
      }
      else if (0==strcmp(argv[i],"-name")){
	number=strtod(argv[i+3],&spointer)/strtod(argv[i+2],&spointer);
	names=(char**) malloc(number*sizeof(char*));
	for (j=strtod(argv[i+2],&spointer),k=0;
	     j<=strtod(argv[i+3],&spointer);
	     j+=strtod(argv[i+2],&spointer),k++){
	  names[k]=malloc(256*sizeof(char));
	  sprintf(names[k],argv[i+1],j);
	}
	i+=4;
      }
      else {
	printf("Option %s not recognized\n",argv[i]);
	i++;
	errno=1;
      }
      if (errno!=0) usage("Error in number conversion!");
      if ((i>=argc)&&(names==NULL)) usage("No files to process!");
    }
   
    if (cut==NULL) cut=(double *) malloc(nocuts*sizeof(double));
    if (names==NULL){
      names = &(argv[i]);
      number= argc-i;
    }
  } else usage("No Arguments!");

  InitGraphics(dummy);
  command = (char *) malloc(256*number*sizeof(char));

  for (i=0;i<number;i++){
    compress = 0;
    spointer=strstr(names[i],".gz");
    if (spointer!=NULL){
      sprintf(command,"gunzip %s",names[i]);
      printf("%s\n",command);
      system(command);
      compress = 1;
      *spointer='\0';
    }
    ReadData(names[i],&datap,&Xsize,&Ysize);
    if (compress == 1){
      sprintf(command,"gzip %s",names[i]);
      printf("%s\n",command);
      system(command);
    }
    Contour_print("epsfile.eps",1,1,0,0,nocuts,cut,datap,Ysize,Xsize,
		  0,"comment",magnify,center,newx,newy,xofs,yofs,
		  0,1,-1,1,1);
    sprintf(command,"gs -sDEVICE=ppm -r72x72 -q -dNOPAUSE -sOutputFile='-' -dBATCH epsfile.eps | pnmcrop > pic.ppm");
    system(command);
    /* There is an undocumented restriction on ppmtoy4m which requires
       the dimensions of the ppm image to be powers of two. I could not
       find any documentation of this problem, but this is supposed to be
       a workaround. First we determine the size of the ppm image and then,
       using pnmcut, we reduce its size as required.*/
    ppmfile=fopen("pic.ppm","r");
    fscanf(ppmfile,"%s %i %i",&(ppmstr[0]),&ppmx,&ppmy);
    fclose(ppmfile);
    printf("Analyse pic.ppm: %s %i %i\n",ppmstr,ppmx,ppmy);
    if (((ppmx&1)==1) || ((ppmy&1)==1)){
      sprintf(command,
	      "cat pic.ppm | pnmcut -width %i -height %i >pic_%05i.ppm",
	      (ppmx>>1)<<1,(ppmy>>1)<<1,i);
      system(command);
    } else {
      sprintf(command,"mv pic.ppm pic_%05i.ppm",i);
      system(command);
    }
  }

  sprintf(command,"cat ");
  for (i=0;i<number; i++) {
    sprintf(addstr,"pic_%05i.ppm ",i);
    strcat(command,addstr);
  }
  strcat(command," | ppmtoy4m -F 24:1 | mpeg2enc -o anim.m1v");
  printf("%s\n",command);
  system(command);
 
  CloseGraphics(dummy);
  return 0;
}