Changeset 3755 in Sophya for trunk/AddOn

Timestamp:

Mar 10, 2010, 3:32:34 PM (16 years ago)

Author:

ansari

Message:

Amelioration cosmetique de tstrdfits.cc et ajout fichier makef_auton pour compilation autonome de tstrdfits.cc avec minifits.cc brpaqu.cc, Reza 10/03/2010

Location:

trunk/AddOn/TAcq

Files:

• makef_auton (added)
• tstrdfits.cc (modified) (4 diffs)

trunk/AddOn/TAcq/tstrdfits.cc

@@ +1 @@
+/*    Projet BAORadio    
+    (C) LAL / Univ. Paris Sud ,   IRFU   2007-2010 
+    Programme de lecture des paquets Acquisition BAORadio 
+*/
 
 // include standard c/c++
@@ -8 +12 @@
 
 
-// include mini-fits lib 
+// include mini-fits lib et manipulation de paquets AcqBAORadio 
 #include "minifits.h"
+#include "brpaqu.h"
 
 int main(int narg, char* arg[])
@@ -19 +24 @@
   cout << " ---------- tstrdfits.cc Start ------------- " << endl;
   int rc = 0;
+  int maxprt = 7;
   try {
+    int paqsz = 0;
+    int nprt = 0;
+    BRPaquet paq; 
+    BRPaqChecker pcheck;
+    Byte* data=NULL;
+//---- Boucle sur les fichiers a lire ---- 
     for(int ifile=1; ifile<narg; ifile++) {
       string ffname = arg[ifile];
@@ -30 +42 @@
         cout << " PB : DataType!=MF_Byte --> skipping " << endl;
       }
-      size_t sx = mff.NAxis1(); 
-      size_t sy = mff.NAxis2();
-      Byte* data = new Byte[sx];
-      for(int j=0; j<sy; j++) {
-        mff.ReadB(data, sx, j*sx);
-        cout << " j=" << j << " data[0...3]=" << (int)data[0] << " , " 
-             << (int)data[1] << " , " <<  (int)data[2] << endl;
-      }
-      cout << "---- FIN lecture " << ffname << endl;
-    }   
+     if (paqsz == 0)  {  // premier passage, on fixe la taille de paquet et on alloue le buffer
+       paqsz = mff.NAxis1();
+       cout << " tstrdfits/ Allocating data , PaqSz=" << paqsz << endl;
+       data = new Byte[paqsz];
+       for(int ib=0; ib<paqsz; ib++) data[ib]=0;
+       paq.Set(data, paqsz); 
+     }
+     else {   // Sinon, on verifie que les differents fichiers ont la meme taille de paquet
+       if (paqsz != mff.NAxis1()) {
+         cout << " PB : paqsz=" << paqsz << " != mff.NAxis1()" << " --> skipping " << endl;
+         continue;
+       }
+     }
+  
+     size_t sx = mff.NAxis1(); 
+     size_t sy = mff.NAxis2();
+     // On construit un paquet sur 
+     // Boucle de lecture des paquets dans le fichier 
+     for(int j=0; j<sy; j++) {
+       mff.ReadB(data, sx, j*sx);    // Lecture du paquet 
+       pcheck.Check(paq);            // Verification, comptage des paquets OK 
+       if ((j%55==0)&&(nprt<maxprt)) {    // Impression optionnel du contenu du paquet 
+         paq.Print();    nprt++; 
+       }
+     }
+     cout << "---- FIN lecture fichier " << ffname << endl;
+    }
+    cout << "-------- FIN boucle de lecture des fichiers ----------" << endl;
+    cout << pcheck;
+    cout << "------------------------------------------------------" << endl;
+    delete[] data;
   }
   catch (MiniFITSException& exc) {