libDataConversion.cpp File Reference

#include <stdio.h>
#include <iostream>
#include <string.h>
#include "TSystem.h"
#include "EdbLog.h"
#include "libDataConversion.h"
#include "dataIO/dataIO.cpp"

Include dependency graph for libDataConversion.cpp:

Functions
int	AddGrainsTXT (EdbRun *run, char *txtname)
int	AddMAP (EdbRun *run, char *mapname)
int	AddRWC (EdbRun *run, char *rwcname, int bAddRWD, const char *options)
int	AddRWD (EdbRun *run, char *rwdname, int fragID, const char *options)
double	FindConfig (IO_VS_Catalog *pCat, char *ConfigName, char *ConfigItem)

Function Documentation

AddGrainsTXT()

int AddGrainsTXT	(	EdbRun *	run,
		char *	txtname
	)

{
  printf("Testrun!\n");
  FILE* grfile;
  grfile=fopen(txtname,"r");
 
  int gs,gf,gTot,gn;                 //grains: view side frame Total index size
  float gX,gY,gZ,ga,gv;
 
  EdbRunHeader *Header;
  Header =  run->GetHeader();
  // customize run flags:
  //  eFlag[0] = 1  - UTS data
  //  eFlag[0] = 2  - SySal data
  //  eFlag[1] = 1  - Stage coordinates
  //  eFlag[1] = 2  - Absolute (marks) coordinates
  //  eFlag[2] = 1  - real (stage) coordinates for clusters
  //  eFlag[2] = 2  - pixels coordinates for clusters
  Header->SetFlag(0,2);  
  Header->SetFlag(1,1);
  Header->SetFlag(2,2);
 
  char instr[128];
  EdbView *View;
  View=run->GetView();
  View->Clear();
  float curv=0;
  int curs=1;
  int curf= -9 ;
  int curTot=0;
 
  int ngr=0;
  int nfr=0;
 
  if(fgets(instr,128,grfile)!=NULL)     //get initial side and viewID
    { 
      sscanf(instr,"%f %d %d %f %d %d %f %f %f",&gv,&gs,&gf,&gZ,&gTot,&gn,&gX,&gY,&ga);
    }
  curs = gs ;
  curv = gv ;
  fseek(grfile,0,0);
  printf("Reading grains from %s..\n",txtname);
  printf("Filling view %d side %d..",int(curv),curs);
 
  while(fgets(instr,128,grfile)!=NULL)
    {
      sscanf(instr,"%f %d %d %f %d %d %f %f %f",&gv,&gs,&gf,&gZ,&gTot,&gn,&gX,&gY,&ga);
      if ((gv!=curv)||(gs!=curs)) {     // new view found   
    // ADD VIEW
    View->GetHeader()->SetViewID((int) curv);
    if(curs)  View->GetHeader()->SetNframes(nfr, 0 );     //curs==1 top  curs==0 bottom
    else      View->GetHeader()->SetNframes(0  ,nfr);
    View->GetHeader()->SetNclusters(ngr);
    View->GetHeader()->SetNsegments(-1);//grain view
    run->AddView(View);
 
    View->Clear();
    curv = gv ;
    curs = gs ;
    printf("%d clusters.\nFilling view %d side %d..",ngr,int(curv),curs);
    ngr=0;
    nfr=0;
      }    
      if ( gf != curf) {  // new frame found
    View->AddFrame(curf, gZ, gTot);
    curf = gf ;
    curTot = gTot ;
    nfr++;
      }
      View->AddCluster(gX,gY,gZ,ga,gv, nfr-1 ,gs,-2);
      ngr++;
    }
 
  // ADD LAST VIEW
  View->GetHeader()->SetViewID((int) curv);
  if(curs)  View->GetHeader()->SetNframes(nfr, 0 );     //curs==1 top  curs==0 bottom
  else      View->GetHeader()->SetNframes(0  ,nfr);
  View->GetHeader()->SetNclusters(ngr);
  View->GetHeader()->SetNsegments(-1);  //grain view
  run->AddView(View);                       //the last view
 
  printf("%d clusters.\n",ngr);
  fclose(grfile);
 
  return true;
}

AddMAP()

int AddMAP	(	EdbRun *	run,
		char *	mapname
	)

{
#ifdef _USESYSAL    
  IO_Data* pMarks = 0;
  if (ReadMap((void**)&pMarks, (char*)mapname ) != 1)   return false;
 
  EdbMarksBox* stage = run->GetMarks()->GetStage() ;
  EdbMarksBox* abs   = run->GetMarks()->GetAbsolute();
    
  int nmarks = pMarks->MkMap.Map.Count ;
  for(int i=0;i<nmarks;i++)
    {
      float AbsX = pMarks->MkMap.Map.pMarks[i].Nominal.X;
      float AbsY = pMarks->MkMap.Map.pMarks[i].Nominal.Y;
      abs->AddMark(i,AbsX,AbsY);
      float StageX = pMarks->MkMap.Map.pMarks[i].Stage.X;
      float StageY = pMarks->MkMap.Map.pMarks[i].Stage.Y;
      stage->AddMark(i,StageX,StageY);
    }
  FreeMemory((void**)pMarks);
  delete pMarks ;
  return true;
#else
  cout << "IO_Data class not implemented for this op. sys. " << endl;
  return false;
#endif 
}

AddRWC()

int AddRWC	(	EdbRun *	run,
		char *	rwcname,
		int	bAddRWD,
		const char *	options
	)

{
  IO_VS_Catalog* pCat = 0;
 
  if( gSystem->AccessPathName(rwcname, kReadPermission) ) {
    Log(1,"AddRWC","ERROR! Can not access file %s",rwcname);
    return 0;
  }
 
  if( ReadCatalog((void**)&pCat, (char*)rwcname) != 1) return false; 
 
  int nFragments = pCat->Area.Fragments;
  cout<<"Hdr.ID :"
      <<pCat->Hdr.ID.Part[0] <<"\t"
      <<pCat->Hdr.ID.Part[1]<<"\t"
      <<pCat->Hdr.ID.Part[2]<<"\t"
      <<pCat->Hdr.ID.Part[3]<<endl;
  cout<<"Options:   "<<options<<endl;
  cout<<"Fragments: "<<nFragments<<endl;
 
  EdbRunHeader *Header;
  Header =  run->GetHeader();
  // customize run flags:
  //  eFlag[0] = 1  - UTS data
  //  eFlag[0] = 2  - SySal data
  //  eFlag[1] = 1  - Stage coordinates
  //  eFlag[1] = 2  - Absolute (marks) coordinates
  //  eFlag[2] = 1  - real (stage) coordinates for clusters
  //  eFlag[2] = 2  - pixels coordinates for clusters
  Header->SetFlag(0,2);  
  Header->SetFlag(1,1);
  Header->SetFlag(2,1);
 
  Header->SetLimits(pCat->Area.XMin,pCat->Area.XMax,
            pCat->Area.YMin,pCat->Area.YMax);
  // fiducial coordinates
  int lys;
  lys=(int)FindConfig(pCat,"Vertigo Scan","VLayers");
  Header->SetArea(pCat->Area.XViews*pCat->Area.YViews, 
          pCat->Area.XStep,pCat->Area.YStep, lys,lys,
          0); 
  //Header->SetNareas(pCat->Area.Fragments);
  Header->SetCCD((int) FindConfig(pCat,"Objective","Width"),
         (int) FindConfig(pCat,"Objective","Height"),
         -999,      // physical pixel size in microns along X
         -999,      // physical pixel size in microns along Y
         "",        // es. "Dalsa CAD4"
         "");       // es. "Bologna"
  Header->SetObjective(-999,        // magnification
               (float)FindConfig(pCat,"Objective","Width" )*
               fabs((float)FindConfig(pCat,"Objective","PixelToMicronX")),
               0,0,
               (float)FindConfig(pCat,"Objective","Height")*
               fabs((float)FindConfig(pCat,"Objective","PixelToMicronY")),
               0,0,
               "",                  // es. "Nikon CFI - oil"
               "");                // es, "50x"
  Header->SetPlate(-999,        // plate ID
           (float)FindConfig(pCat,"Vertigo Scan","VStep")*
           (float)FindConfig(pCat,"Vertigo Scan","Shrinkage"),
           (float)FindConfig(pCat,"Vertigo Scan","BaseThickness"),
           (float)FindConfig(pCat,"Vertigo Scan","VStep")*
           (float)FindConfig(pCat,"Vertigo Scan","Shrinkage"),
           (float)FindConfig(pCat,"Vertigo Scan","Shrinkage"),
           (float)FindConfig(pCat,"Vertigo Scan","Shrinkage"),
           "",                  // es. "Test Plate"
           "");             //
 
  // Store the catalog file as TObjString 
  TObjString objstr ;
  int  i, ch;
  FILE *stream;
  if( (stream = fopen( rwcname, "rb" )) == NULL ) exit( 0 );
  ch = fgetc( stream );
  for( i=0;  feof( stream ) == 0 ; i++ )
    {
      objstr.String().Append( ch ); 
      ch = fgetc( stream );
    }
 
  //    int NviewCols = (int) FindConfig(pCat,"Vertigo Scan","XFields")  ;
  //    int NviewRows = (int) FindConfig(pCat,"Vertigo Scan","YFields")  ;
  //    int NfragmCols = pCat->Area.XViews / NviewCols  ;
  //    int NfragmRows = pCat->Area.YViews / NviewRows  ;
 
  fclose( stream );
  objstr.Write("catalog"); 
 
  // Free memory
  FreeMemory((void**)pCat);
  delete pCat;
 
  run->Save() ;
 
  // loop on rwd files
  if(bAddRWD)
    {
      char rwdname[256], temp[256];
      for (int f = 1; f < nFragments+1; f++)
    {
      sprintf(temp, "%s", rwcname);
      sprintf(temp+strlen(temp)-1, "d");
      sprintf(rwdname, "%s.%08X", temp, f);
      cout <<"(tot. fragm.:"<<nFragments<<")  "; 
      if (! AddRWD(run, rwdname, f, options) )   continue;
    }
      cout <<endl;
    }
  return true;
}

AddRWD()

int AddRWD	(	EdbRun *	run,
		char *	rwdname,
		int	fragID,
		const char *	options
	)

{ 
  if( gSystem->AccessPathName(rwdname, kReadPermission) ) {
    Log(1,"AddRWD","ERROR! Can not access file %s",rwdname);
    return 0;
  }
 
   Bool_t addcl(true);
   Bool_t addclframe(false);
   Bool_t addareasum(false);
   // OPTIONS
   if (strstr(options,"NOCL") )    addcl=false;     // do not add clusters
   if (strstr(options,"CLFRAME") ) addclframe=true; // fill clusters.eFrame branch
   if (strstr(options,"SUM") )     addareasum=true; // encode the sum of cluster areas in the segment puls
   // ePuls = (sum of clust areas)*1000 + (number of clusters) 
   const char *key=0;
   int cutcl=0;
   if( (key = strstr(options,"CUTCL")) )  sscanf(key+5,"%d",&cutcl);
 
   EdbView*    edbView = run->GetView();
   EdbSegment* edbSegment = new EdbSegment(0,0,0,0,0,0,0,0);
 
   IO_VS_Fragment2* pFrag = 0;
   if (ReadFragment((void**)&pFrag, (char*)rwdname) != 1)   return false;
 
   int v, s, t, p;   // v=view, s=side, t=track, p=point
   int tracks;         // number of tracks in the fragment
   int fclusters;      // number of clusters in the fragment 
   int vclusters;      // number of clusters in the view
   float dz ;          // z-length of the track segment
   int tr_clusters; // number of cluster of the track
   int puls ;        // ePuls = (number of clusters ) OR  (sum of clust areas)*1000 + (number of clusters)
 
   Track2* rwdTrack;
   VS_View2* rwdView;
 
   tracks = 0;
   fclusters = 0;
   run->GetHeader()->SetNareas(run->GetHeader()->GetNareas()+1);
   for (s = 0; s < 2; s++) {
      for (v = 0; v < pFrag->Fragment.CountOfViews; v++) { 
      rwdView = &(pFrag->Fragment.pViews[v]);
      vclusters=0;
      tracks += rwdView->TCount[s];
 
      edbView->Clear();
      EdbViewHeader* edbViewHeader = edbView->GetHeader();
      edbViewHeader->SetAffine(rwdView->ImageMat[s][0][0],
                   rwdView->ImageMat[s][0][1],
                   rwdView->ImageMat[s][1][0],
                   rwdView->ImageMat[s][1][1],
                   rwdView->MapX[s], 
                   rwdView->MapY[s]); 
      edbViewHeader->SetAreaID(run->GetHeader()->GetNareas());
      edbViewHeader->SetCoordXY(rwdView->X[s], rwdView->Y[s]); 
      edbViewHeader->SetCoordZ(rwdView->RelevantZs.TopExt,rwdView->RelevantZs.TopInt,
                   rwdView->RelevantZs.BottomInt,rwdView->RelevantZs.BottomExt);
 
      if(s) edbViewHeader->SetNframes(0,rwdView->Layers[s].Count); // s==0 top , s==1 bottom
      else  edbViewHeader->SetNframes(rwdView->Layers[s].Count,0);
 
      edbViewHeader->SetNsegments(rwdView->TCount[s]);
      edbViewHeader->SetViewID(v);
      edbViewHeader->SetColRow(rwdView->TileX  ,rwdView->TileY );
 
      int nframes = rwdView->Layers[s].Count;
 
      for( int nlvl=0; nlvl<nframes; nlvl++ )
          edbView->AddFrame(nlvl,rwdView->Layers[s].pLayerInfo[nlvl].Z,
                    rwdView->Layers[s].pLayerInfo[nlvl].Clusters);
 
         for (t = 0; t < rwdView->TCount[s]; t++) {
            rwdTrack = &(rwdView->pTracks[s][t]);
            tr_clusters  = rwdTrack->Grains;
        if( cutcl>0 && tr_clusters<cutcl) continue;
            if (addareasum)   puls = rwdTrack->Grains + (2000000>rwdTrack->AreaSum?rwdTrack->AreaSum:2000000)*1000;
            else              puls = rwdTrack->Grains ;
 
            if ( rwdTrack->pGrains )      // fix the "null pointer to grains" bug
               dz = rwdTrack->pGrains[0].Z - rwdTrack->pGrains[tr_clusters-1].Z;
            else
               dz = 0 ;
            
            edbSegment->Set(  rwdTrack->Intercept.X,
                               rwdTrack->Intercept.Y,
                               rwdTrack->Intercept.Z,
                               rwdTrack->Slope.X,
                               rwdTrack->Slope.Y, 
                               dz, s , puls, t);
                          
            edbSegment->SetSigma(rwdTrack->Sigma,-999);
 
            // Add clusters 
            if (addcl && rwdTrack->pGrains ) {       // rwdTrack->pGrains -> fix the "null pointer to grains" bug
               float volume = 0;
               int frameId  = 0;
               for ( p=0; p<tr_clusters;p++)  {
                  if (addclframe) {
                     for(int nlvl=frameId;nlvl<nframes;nlvl++) {
                        if(rwdTrack->pGrains[p].Z == rwdView->Layers[s].pLayerInfo[nlvl].Z ) {
                           frameId=nlvl;
                           break;
                        }
                     }
                  }
                  edbView->AddCluster( rwdTrack->pGrains[p].X,
                                       rwdTrack->pGrains[p].Y,
                                       rwdTrack->pGrains[p].Z,
                                       rwdTrack->pGrains[p].Area,
                                       volume,
                                       frameId,
                                       s,t);                                        
               }
            }
            edbViewHeader->SetNclusters(vclusters);
            edbView->AddSegment(edbSegment) ;
            vclusters += tr_clusters;
         } // end of tracks (t)
         run->AddView(edbView);
         fclusters += vclusters;
      }; //end of views (v)
   };//end of sides (s)
 
  cout<<"Fragment:"<<fragID<<"\tmicrotracks: "
      <<tracks<<"\tclusters: "<<fclusters<<endl;
  cout << flush;
   
  FreeMemory((void**)pFrag);
#ifndef WIN32
  delete pFrag;    // WINDOWS VISTA CRASH but on linux - memory leak!!
#endif
  return true;
}

FindConfig()

double FindConfig	(	IO_VS_Catalog *	pCat,
		char *	ConfigName,
		char *	ConfigItem
	)

{
  int CountOfConfigs = pCat->Config.CountOfConfigs;
  VS_Config* pConfigs =pCat->Config.pConfigs;
 
  double dbl=0;
  for (int j=0;j<CountOfConfigs;j++)    {
    if(! strcmp(pConfigs[j].Config.Name, ConfigName) ) {
      for(int k=0;k<pConfigs[j].Config.CountOfItems;k++) {
        if(! strcmp(&pConfigs[j].Config.pItems[k*128],ConfigItem) ) {
      char* stopchar;
      dbl = strtod(&pConfigs[j].Config.pItems[k*128+64],&stopchar);
      break;
    }
      }
    }
  }
  return dbl;
}