Changeset 6712 for branches/rel10_ifa/psModules/src/objects/pmFPAviewWriteObjects.c

Timestamp:

Mar 27, 2006, 4:14:56 PM (20 years ago)

Author:

magnier

Message:

substantial work on the pmFPAfile I/O concepts

File:

• branches/rel10_ifa/psModules/src/objects/pmFPAviewWriteObjects.c (modified) (4 diffs)

branches/rel10_ifa/psModules/src/objects/pmFPAviewWriteObjects.c

@@ -3 +3 @@
  *  @author EAM, IfA; GLG, MHPCC
  *
- *  @version $Revision: 1.1.2.3 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-03-14 02:21:07 $
+ *  @version $Revision: 1.1.2.4 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-28 02:14:56 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -27 +27 @@
 #include "pmFPAviewObjectsIO.h"
 
+bool psFitsWriteHeaderNotImage (psFits *fits, psMetadata *header)
+{
+
+    int status = 0;
+    char fitsErr[1024];
+
+    psMetadataItem *mdi = NULL;
+    mdi = psMetadataLookup (header, "NAXIS");
+    if (mdi != NULL) {
+        mdi->data.S32 = 0;
+        mdi->type = PS_DATA_S32;
+    }
+
+    fits_create_img (fits->fd, SHORT_IMG, 0, NULL, &status);
+
+    if (status != 0) {
+        fits_get_errstatus (status, fitsErr);
+        psError (PS_ERR_IO, true, PS_ERRORTEXT_psFits_WRITE_FAILED, fitsErr);
+        exit (2);
+    }
+
+    psFitsWriteHeader (header, fits);
+
+    return true;
+}
+
 // Given a FITS file pointer, read the table of object data
 bool pmFPAviewWriteObjects (pmFPAview *view, pmFPAfile *file)
 {
 
-    pmFPA *fpa = view->fpa;
+    pmFPA *fpa = file->fpa;
 
     if (view->chip == -1) {
@@ -112 +138 @@
     char *headname;
     pmHDU *hdu;
+    pmHDU *phu;
 
     psArray *sources = psMetadataLookupPtr (&status, readout->analysis, "PSPHOT.SOURCES");
+    return true;
 
     switch (file->type) {
+    case PM_FPA_FILE_RAW:
+        filename = pmFPAfileNameFromRule (file->filerule, file, view);
+        pmSourcesWriteRAW (sources, filename);
+        psFree (filename);
+        break;
+
     case PM_FPA_FILE_OBJ:
-        filename = pmFPAviewNameFromRule (file->filerule, view, file->fpa);
+        filename = pmFPAfileNameFromRule (file->filerule, file, view);
         pmSourcesWriteOBJ (sources, filename);
+        psFree (filename);
         break;
 
     case PM_FPA_FILE_SX:
-        filename = pmFPAviewNameFromRule (file->filerule, view, file->fpa);
+        filename = pmFPAfileNameFromRule (file->filerule, file, view);
         pmSourcesWriteSX (sources, filename);
+        psFree (filename);
         break;
 
@@ -129 +165 @@
         // a SPLIT format : only one header and object table per file
         hdu = pmFPAviewThisHDU (view, file->fpa);
-        filename = pmFPAviewNameFromRule (file->filerule, view, file->fpa);
+        filename = pmFPAfileNameFromRule (file->filerule, file, view);
         pmSourcesWriteCMP (sources, filename, hdu->header);
+        psFree (filename);
         break;
 
     case PM_FPA_FILE_CMF:
-        // write header, if not yet written?
+        // write a PHU?
+        // write a header?
+        // write the data
+
+        // get the current header
         hdu = pmFPAviewThisHDU (view, file->fpa);
 
-        // XXX is this needed? is it automatically added?
-        // psMetadataAdd (header, PS_LIST_TAIL, "EXTEND", PS_DATA_BOOL | PS_META_REPLACE, "this file has extensions", true);
-
-        // I need a test to see if the Cell header has been written out
-        if (hdu->header) {
-            psMetadataItem *mdi = NULL;
-            headname = pmFPAviewNameFromRule (file->extxtra, view, file->fpa);
-
-            // set NAXIS to 0 (we don't write out the data array)
-            // XXX part of the header cleanup?
-            mdi = psMetadataLookup (hdu->header, "NAXIS");
-            mdi->data.S32 = 0;
-            mdi->type = PS_DATA_S32;
-
-            if (strcasecmp (headname, "PHU")) {
-                psMetadataAdd (hdu->header, PS_LIST_HEAD, "EXTNAME", PS_DATA_STRING, "extension name", headname);
-            }
-            psFitsWriteHeader (hdu->header, file->fits);
+        // if file does not yet have a PHU, write it to disk
+        if (file->phu == NULL) {
+            // get the corresponding phu
+            phu = pmFPAviewThisPHU (view, file->fpa);
+
+            // CMF always has extensions
+            psMetadataAdd (phu->header, PS_LIST_TAIL, "EXTEND", PS_DATA_BOOL | PS_META_REPLACE, "this file has extensions", true);
+            psFitsWriteHeaderNotImage (file->fits, phu->header);
+            file->phu = phu->header;
+            psTrace ("pmFPAfile", 5, "wrote phu %s (type: %d)\n", file->filename, file->type);
+
+            file->phu = phu->header;
         }
-        dataname = pmFPAviewNameFromRule (file->extrule, view, file->fpa);
+
+        // this this header block is new, write it to disk,
+        if (hdu->header != file->header) {
+            // determine name for header extension
+            headname = pmFPAfileNameFromRule (file->extxtra, file, view);
+            psMetadataAddStr (hdu->header, PS_LIST_TAIL, "EXTNAME", PS_META_REPLACE, "extension name", headname);
+            psFitsWriteHeaderNotImage (file->fits, hdu->header);
+            psTrace ("pmFPAfile", 5, "wrote ext head %s (type: %d)\n", file->filename, file->type);
+            file->header = hdu->header;
+            psFree (headname);
+        }
+
+        // write this table to disk
+        dataname = pmFPAfileNameFromRule (file->extrule, file, view);
         pmSourcesWriteCMF (file->fits, sources, hdu->header, dataname);
+        psFree (dataname);
         break;