Changeset 3688 in Sophya for trunk/AddOn/TAcq/minifits.cc

Timestamp:

Nov 29, 2009, 12:13:54 AM (16 years ago)

Author:

ansari

Message:

Introduction de la classe FitsCubeWriter - Reza 28/11/2009

File:

• trunk/AddOn/TAcq/minifits.cc (modified) (3 diffs)

trunk/AddOn/TAcq/minifits.cc

@@ -220 +220 @@
   buff[strlen(buff)] = ' ';
   buff = header+320;
-  sprintf(buff, "NAXIS2  = %20ld / nb of rows = NumberOfPaquets", (long)nax2);
+  if (fgnax3) 
+    sprintf(buff, "NAXIS2  = %20ld / Number of fibers ", (long)nax2);
+  else 
+    sprintf(buff, "NAXIS2  = %20ld / nb of rows = NumberOfPaquets", (long)nax2);
   buff[strlen(buff)] = ' ';
   buff = header+400;
   if (fgnax3) 
-    sprintf(buff, "NAXIS3  = %20ld / nb of rows = NumberOfFibers", (long)nax3);
+    sprintf(buff, "NAXIS3  = %20ld / Number of paquets ", (long)nax3);
   else 
     strcpy(buff,"COMMENT  BAO-Radio / MiniFITSFile ");
@@ -311 +314 @@
 
 /* --Methode-- */
+string MiniFITSFile::GetKey(const char* key)
+{
+  char rs[80];
+  char cle[10];
+  strncpy(cle,key,8);
+  cle[8]='=';
+  for(int i=0;i<8;i++) 
+    if (cle[i]=='\0') cle[i]=' ';
+  cle[9]='\0';
+  for(int kh=80; kh<2800; kh+=80) {
+    char * buff = header+kh;
+    if (strncmp(buff, cle, 9) == 0)  {
+      strncpy(rs, buff+10, 70);
+      rs[70]='\0';
+      break;
+    }
+  }
+  return rs;
+}
+
+/* --Methode-- */
 void MiniFITSFile::DecodeHeader()
 {
@@ -322 +346 @@
   bool fgok3=false;
   for(int kh=80; kh<2800; kh+=80) {
-        buff = header+kh;
-        if (strncmp(buff, "NAXIS1  =", 9) == 0)  {
-          nax1 = atol(buff+10);
-          fgok1 = true;
-//        cout << " FOUND : NAXIS1= " << nax1 << endl;
-        } 
-        else if (strncmp(buff, "NAXIS2  =", 9) == 0)  {
-          nax2 = atol(buff+10);
-          fgok2 = true;
-//        cout << " FOUND : NAXIS2= " << nax2 << endl;            
-        } 
-        else if (strncmp(buff, "NAXIS3  =", 9) == 0)  {
-          nax3 = atol(buff+10);
-          fgok3 = true;
-//        cout << " FOUND : NAXIS2= " << nax2 << endl;            
-        } 
-        else if (strncmp(buff, "BITPIX  =", 9) == 0)  {
-          int bpix = atoi(buff+10);
-          fgokt = true;
-          if (bpix == 8)  dtype = MF_Byte;
-          else if (bpix == 16)  dtype = MF_Int16;
-          else if (bpix == -32)  dtype = MF_Float32;
-          else fgokt = false;
-//        cout << " FOUND : bpix= " << bpix << endl;              
-        } 
+    buff = header+kh;
+    if (strncmp(buff, "NAXIS1  =", 9) == 0)  {
+      nax1 = atol(buff+10);
+      fgok1 = true;
+      //          cout << " FOUND : NAXIS1= " << nax1 << endl;
+    } 
+    else if (strncmp(buff, "NAXIS2  =", 9) == 0)  {
+      nax2 = atol(buff+10);
+      fgok2 = true;
+      //          cout << " FOUND : NAXIS2= " << nax2 << endl;            
+    } 
+    else if (strncmp(buff, "NAXIS3  =", 9) == 0)  {
+      nax3 = atol(buff+10);
+      fgok3 = true;
+      //          cout << " FOUND : NAXIS2= " << nax2 << endl;            
+    } 
+    else if (strncmp(buff, "BITPIX  =", 9) == 0)  {
+      int bpix = atoi(buff+10);
+      fgokt = true;
+      if (bpix == 8)  dtype = MF_Byte;
+      else if (bpix == 16)  dtype = MF_Int16;
+      else if (bpix == -32)  dtype = MF_Float32;
+      else fgokt = false;
+      //          cout << " FOUND : bpix= " << bpix << endl;              
+    } 
   }
   if (!(fgok1&&fgok2&&fgokt))  
-     throw MiniFITSException("MiniFITSFile::DecodeHeader()/Error- Missing/wrong NAXIS1/2,BITPIX");
+    throw MiniFITSException("MiniFITSFile::DecodeHeader()/Error- Missing/wrong NAXIS1/2,BITPIX");
   if (fgok3)  fgnax3=true;
   return;