Changeset 2843 in Sophya for trunk/SophyaExt/FitsIOServer/fitsblkrw.h

Timestamp:

Nov 18, 2005, 6:44:57 PM (20 years ago)

Author:

ansari

Message:

debut modifs nouveau fits pour support I/O chaines de caracteres - Reza 18/11/2005

File:

• trunk/SophyaExt/FitsIOServer/fitsblkrw.h (modified) (1 diff)

trunk/SophyaExt/FitsIOServer/fitsblkrw.h

@@ -114 +114 @@
 };
 
+
+/*! 
+  Write character string data to binary/ascii HDU data in a fits file. 
+  See cfitsio function fits_write_col() for more detail.
+  \param colnum : table column number (starting from 1)
+  \param firstrow : the write operation starting row (starting from 1)
+  \param firstelem : the firstelem (for vector type columns)
+  \param d : pointer to string type array to be written
+  \param sz : number of data elements to be written
+  \param width : column width
+*/
+void WriteStringColumnData(FitsInOutFile& fios, int colnum, long firstrow,
+                           long firstelem, const string * d, size_t sz, 
+                           long width=0)
+{
+  int status = 0;
+  if (width < 1) width = 16;
+  for(size_t kk=0; kk<sz; kk++) {
+    char * cp = const_cast<char *>(d[kk].c_str());
+    status = 0;
+    fits_write_col(fios.FitsPtr(), FitsTypes::DataType(cp), colnum,
+                   firstrow+kk, firstelem, 1, &cp, &status);
+    if ( status ) {
+      fits_report_error(stderr, status);
+      char buff[32];
+      fits_get_errstatus(status, buff);
+      string msg = "WriteStringColumnData Error: " ;
+      msg += buff;
+      sprintf(buff," kk=%ld",kk);  msg += buff;
+      throw FitsIOException(msg);
+    }
+  }
+  return;
+}
+
+/*! 
+  Read character string data to binary/ascii HDU data in a fits file. 
+  See cfitsio function fits_read_col() for more detail.
+  \param colnum : table column number (starting from 1)
+  \param firstrow : the read operation starting point (row) (starting from 1)
+  \param firstelem : the firstelem (for vector type columns)
+  \param d : pointer to string type array to be read
+  \param sz : number of data elements to be read
+  \param width : column width
+*/
+void ReadStringColumnData(FitsInOutFile& fios, int colnum, long firstrow,
+                          long firstelem, string * d, size_t sz, long width)
+{
+  int status = 0;
+  int anynul = 0;
+  char buff[1024];
+  for(size_t kk=0; kk<sz; kk++) {
+    fits_read_col(fios.FitsPtr(), FitsTypes::DataType(buff), colnum,
+                 firstrow, firstelem+kk, 1, NULL, &buff, &anynul, &status);
+    d[kk] = buff;
+    if ( status ) {
+      fits_report_error(stderr, status);
+      char buff[32];
+      fits_get_errstatus(status, buff);
+      string msg = "ReadStringColumnData Error: " ;
+      msg += buff;
+      sprintf(buff," kk=%ld",kk);  msg += buff;
+      throw FitsIOException(msg);
+    }
+  }
+  return;
+}
+
 } // Fin du namespace