/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Name: Subdivsion
 Copyright: 
 Authors: COFFIGNAL Thomas
 Date: 24/02/09
 Description: Subdivision de l'image en paves
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/

#include <stdio.h>
#include <limits.h>
#include <float.h>
#include <math.h>
#include <string.h>
#include <stdlib.h>	
#define MAX 256
#define imageParPointeur(i,j,k) (image_ptr + (long int)(i) + largeur*(long int)(j) + largeur*hauteur*(long int)(k))
#define petitCubeParPointeur(i,j,k) (petitCube_ptr + (long int)(i) + largeurPetitCube*(long int)(j) + largeurPetitCube*hauteurPetitCube*(long int)(k))


main()
{   
	
	//Declaration de variables
	FILE *entree=NULL;
    FILE *sortie;
	
	char  nomimage[MAX],nomsortie[MAX];
	
	int largeur, hauteur, profondeur ;
    int largeurPetitCube, hauteurPetitCube, profondeurPetitCube ;
    int i,j,k;    
    int nbCubesDansLargeur, nbCubesDansHauteur, nbCubesDansProfondeur ;
    int iDecalage, jDecalage, kDecalage ;
    int iCube, jCube, kCube ;
    int iDansImage,jDansImage,kDansImage ;
   	int xmin,ymin,zmin,xmax,ymax,zmax;
    
    long int taille, taillePetitCube ;
    
    unsigned char *image_ptr,*petitCube_ptr ;
	
	//lecture du nom du fichier a traiter
	scanf("%s\n",&nomimage);   
	entree = fopen(nomimage,"rb");
	
	//Effacement des petits cubes precedents
	system("del /Q Resultats\\cubes_decoupes\\*.*");
	
	//Recuperation des donnees du fichier inputPCE 
	scanf("%d %d %d\n",&xmin,&ymin,&zmin);
	scanf("%d %d %d\n",&xmax,&ymax,&zmax);	
	scanf("%d %d %d\n",&largeurPetitCube,&hauteurPetitCube,&profondeurPetitCube);  
	
	//Parametrage dimension de l'image
	largeur=xmax-xmin+1;
	hauteur=ymax-ymin+1;
	profondeur=zmax-zmin+1;
	taille = largeur*hauteur*profondeur ;
	
	//Parametrage dimension des pave
    taillePetitCube = largeurPetitCube*hauteurPetitCube*profondeurPetitCube; 
	
	//Decoupage	
	nbCubesDansLargeur = largeur/largeurPetitCube;
    nbCubesDansHauteur =  hauteur/hauteurPetitCube;
    nbCubesDansProfondeur = profondeur/profondeurPetitCube;
	
	//Affichage des donnees de decoupe
	printf("\til y aura %d cubes selon x\n",nbCubesDansLargeur);
	printf("\til y aura %d cubes selon y\n",nbCubesDansHauteur);
	printf("\til y aura %d cubes selon z\n",nbCubesDansProfondeur);
    printf("\tLargeur du petit cube = %d\n",largeurPetitCube);
	printf("\tHauteur du petit cube = %d\n",hauteurPetitCube);
	printf("\tProfondeur du petit cube =%d\n",profondeurPetitCube); 
	
	//Allocation du pointeur de l'image
	image_ptr = (unsigned char *)malloc(taille*sizeof(unsigned char));
	
	//Lecture du pointeur de l'image
	fread(image_ptr,sizeof(unsigned char),taille,entree);
    
    //Fermeture du fichier
	fclose(entree);
    
    //Allocation du pointeur du petit cube
    petitCube_ptr = (unsigned char *)malloc(taillePetitCube*sizeof(unsigned char));
		  
    //Debut de la subdivision
    for (iCube=0;iCube<nbCubesDansLargeur;iCube++)
        { 
          iDecalage = iCube*largeurPetitCube ;
          for (jCube=0;jCube < nbCubesDansHauteur;jCube++)
              {
              jDecalage = jCube*hauteurPetitCube ;
              for (kCube=0;kCube < nbCubesDansProfondeur;kCube++)
                  {
                  kDecalage = kCube*profondeurPetitCube ;
                  for ( i=0; i < largeurPetitCube; i++ )
                      {
                      iDansImage = i+iDecalage ;
                      for ( j=0; j < hauteurPetitCube; j++ )
                          {
                          jDansImage = j+jDecalage ;             
                          for ( k=0; k < profondeurPetitCube;k++ )
                              {
                              kDansImage = k+kDecalage ;           
                              *petitCubeParPointeur(i,j,k)=*imageParPointeur(iDansImage,jDansImage,kDansImage);
                              }
                          }
                      }
                      //Creation des paves_iCube_jCube_kCube.mat
   	                  sprintf(nomsortie,"%s%s%d%s%d%s%d%s","Resultats\\cubes_decoupes\\paves","_",iCube,"_",jCube,"_",kCube,".mat");
       	              
                      //Ouverture fichier binaire 
                      sortie=fopen(nomsortie,"ab");
                      
                      //Ecriture de l'en-tete
                      fprintf(sortie,"%3d\t%3d\t%3d\n",largeurPetitCube,hauteurPetitCube,profondeurPetitCube);
       	              fflush(sortie);
       	                     	              
                      //Affichage ecran
       	              printf("\t\npaves_%d_%d_%d.mat",iCube, jCube, kCube);      
                      
                      //Ecriture du fichier binaire                       	                      	                      
                      fwrite(petitCube_ptr,sizeof(unsigned char),taillePetitCube,sortie); 
                      
                      //Fermeture fichier binaire
                      fclose(sortie);                              
                  }
              }
         }       
}