Changeset 2898 in Sophya for trunk/SophyaExt/FitsIOServer/fitsmanager.cc

Timestamp:

Jan 13, 2006, 7:26:40 PM (20 years ago)

Author:

ansari

Message:

Ajout methode FitsManager::ScanFile() et corrections diverses - Reza 13/01/2006

File:

• trunk/SophyaExt/FitsIOServer/fitsmanager.cc (modified) (3 diffs)

trunk/SophyaExt/FitsIOServer/fitsmanager.cc

@@ -75 +75 @@
   for(it = hlistp->begin(); it != hlistp->end(); it++) {
     hfg = (*it).fhi->CheckHandling(o);
+    if (hfg < 1) continue;
     if ( ( hfg > bhfg ) || ( (hfg == bhfg) && ((*it).glev > blev) ) ) {
       fhi = (*it).fhi;  bhfg = hfg;  blev = (*it).glev;
@@ -114 +115 @@
   for(it = hlistp->begin(); it != hlistp->end(); it++) {
     hfg = (*it).fhi->CheckReadability(is);
+    if (hfg < 1) continue;
     if ( ( hfg > bhfg ) || ( (hfg == bhfg) && ((*it).glev > blev) ) ) {
       fhi = (*it).fhi;  bhfg = hfg;  blev = (*it).glev;
@@ -136 +138 @@
   return fhi2;
 }
+
+/*!
+  \param filename : FITS file name to be scanned
+  \param os : infomation will be sent to formatted stream os
+  \param slev : scan level , bit 0 (1/3) print HDU keywords,
+  bit 2 (2,3) try to read HDU data using the appropraite handler
+  \param Rc : return number of scanned HDU's
+ */
+int FitsManager::ScanFile(string filename, ostream& os, int slev)
+{
+  FitsInOutFile is(filename, FitsInOutFile::Fits_RO);
+  os << "=== FitsManager::ScanFile( " << filename << " ) NbHDUs= " 
+     << is.NbHDUs() << endl;
+  int rc = 0;
+  for(int k=0; k<is.NbHDUs(); k++) {
+    os << " ------ HDU No " << is.CurrentHDU() << " Type= " 
+         << is.CurrentHDUTypeStr() << endl;
+    int hdutyp = is.CurrentHDUType();
+    if (hdutyp == IMAGE_HDU) {
+      long naxes[5] = {0,0,0,0,0};
+      int naxis=5;
+      int imgtyp = is.GetImageHDUInfo(naxis, naxes);
+      os << ">> IMAGE_HDU:  naxis= " << naxis << " : ";
+      for(int i=0; i<naxis; i++) {
+        if (i>0) os << " x " ;
+        os << naxes[i];
+      }
+      os << endl;
+    }
+    else {
+      vector<string> colnames;
+      vector<int> coltypes;
+      vector<long> repcnt;
+      vector<long> width;
+      int ncols = is.GetColInfo(colnames, coltypes, repcnt, width);
+      if (hdutyp == BINARY_TBL) os << ">> BINARY_TBL :  NRows= " << is.GetNbRows();
+      else os << ">> ASCII_TBL :  NRows= " << is.GetNbRows();
+      os << " x NCols= " << ncols << endl;
+      for(int kk=0; kk<colnames.size(); kk++) {
+        os << "Col[" << kk+1 << "]  Name= " << colnames[kk] 
+             << " Type= " << FitsTypes::DataTypeToTypeString(coltypes[kk])
+             << " Repeat= " << repcnt[kk]
+             << " W= " << width[kk] << endl;
+      }
+    }
+    // Fin the appropriate handler :
+    ChkHLP(); 
+    FitsHandlerInterface * fhi = NULL;
+    HandlerList::iterator it;
+    string hdesc;
+    int hfg = 0;
+    int bhfg = 0;
+    int blev = 0;
+    for(it = hlistp->begin(); it != hlistp->end(); it++) {
+      hfg = (*it).fhi->CheckReadability(is);
+      if (hfg < 1) continue;
+      if ( ( hfg > bhfg ) || ( (hfg == bhfg) && ((*it).glev > blev) ) ) {
+        fhi = (*it).fhi;  bhfg = hfg;  blev = (*it).glev; hdesc = (*it).desc; 
+      }
+    }
+    if (fhi == NULL) 
+      os << ">>> Warning : No handler found for this HDU ... " << endl;
+    else 
+      os << ">>> Reader/handler: " <<  hdesc << " : " 
+         << typeid(*fhi).name() << " HandLevel= " << blev << ", "  << bhfg << endl;
+    if (fhi && (slev >= 2)) {
+      os << ">>> Trying to read HDU data using the handler ..." << endl;
+      FitsHandlerInterface* fhic = fhi->Clone();
+      fhic->Read(is);
+      os << " FitsHandler.Read() OK " << endl;
+    }
+    if ( (slev == 1) || (slev == 3) ) {
+      os << ">>>> HDU keywords list :  " <<  endl;
+      DVList dvl;
+      is.GetHeaderRecords(dvl);
+      os << dvl;
+    }
+    os << "               --------------------- " << endl;
+    is.MoveToNextHDU();
+    rc++;
+  }
+  os << "===================================================" << endl; 
+  return rc;
+}