Changeset 6847 for branches/rel10_ifa/psModules/src/astrom/pmFPAfile.c

Timestamp:

Apr 12, 2006, 8:18:22 PM (20 years ago)

Author:

magnier

Message:

finishing up conversion to pmFPAfile paradigm

File:

• branches/rel10_ifa/psModules/src/astrom/pmFPAfile.c (modified) (13 diffs)

branches/rel10_ifa/psModules/src/astrom/pmFPAfile.c

@@ -2 +2 @@
 #include "pslib.h"
 #include "psAdditionals.h"
+#include "pmConfig.h"
 #include "pmHDU.h"
 #include "pmFPA.h"
@@ -22 +23 @@
     psFree (file->names);
 
+    psFree (file->format);
+    psFree (file->name);
+
     if (file->fits != NULL) {
         psFitsClose (file->fits);
@@ -54 +58 @@
     file->fits = NULL;
     file->names = psMetadataAlloc();
+
+    file->format = NULL;
+    file->name = NULL;
 
     file->filerule = NULL;
@@ -188 +195 @@
     if (fpa != NULL) {
         file->fpa = psMemIncrRefCounter(fpa);
+    }
+
+    // for WRITE type of data, the output format needs to be determined
+    // this is only needed if the output file is not identical in structure to the input
+    char *formatName = psMetadataLookupStr (&status, data, "FILE.FORMAT");
+    if (formatName && strcasecmp (formatName, "NONE")) {
+        psMetadata *formats = psMetadataLookupMD(&status, camera, "FORMATS"); // List of formats
+        char *formatFile = psMetadataLookupStr (&status, formats, formatName);
+        readConfig (&file->format, formatFile, formatName);
     }
 
@@ -332 +348 @@
 bool pmFPAfileWrite (pmFPAfile *file, const pmFPAview *view)
 {
+    // XXX test return false;
+
     if (file->mode != PM_FPA_MODE_WRITE)
         return false;
@@ -372 +390 @@
 bool pmFPAfileCreate (pmFPAfile *file, const pmFPAview *view)
 {
+    // XXX test return false;
     if (file->mode != PM_FPA_MODE_WRITE)
         return false;
@@ -389 +408 @@
     case PM_FPA_FILE_IMAGE:
         /* create a PHU for thie file, if it does not exist */
-        pmFPAfileCopyView (file->src, file->fpa, view);
+        pmFPAfileCopyStructureView (file->fpa, file->src, file->format, file->xBin, file->yBin, view);
         psTrace ("pmFPAfile", 5, "created fpa data elements for %s (fpa: %p)\n", file->filename, file->fpa);
         break;
@@ -698 +717 @@
 
         // set the view to the corresponding entry for this phu
-        pmFPAview *view = pmFPAAddSourceFromHeader(fpa, phu, format);
+        pmFPAview *view = pmFPAAddSourceFromHeader (fpa, phu, format);
 
         // XXX is this the correct psMD to save the filename?
@@ -786 +805 @@
 
         // set the view to the corresponding entry for this phu
-        pmFPAview *view = pmFPAAddSourceFromHeader(fpa, phu, format);
+        pmFPAview *view = pmFPAAddSourceFromHeader (fpa, phu, format);
 
         // XXX is this the correct psMD to save the filename?
@@ -803 +822 @@
 }
 # endif
+
+// create a new output pmFPAfile based on an existing FPA
+pmFPAfile *pmFPAfileFromFPA (pmConfig *config, pmFPA *src, int xBin, int yBin, char *filename)
+{
+    // XXX pmFPAConstruct has many leaks (6919)
+    pmFPA *fpa = pmFPAConstruct (config->camera);
+    pmFPAfile *file = pmFPAfileDefine (config->files, config->camera, fpa, filename);
+    file->src = src; // inherit output elements from this source pmFPA
+    file->xBin = xBin;
+    file->yBin = yBin;
+
+    return file;
+}
 
 // look for the given name on the argument list.
@@ -879 +911 @@
 
         // set the view to the corresponding entry for this phu
-        pmFPAview *view = pmFPAAddSourceFromHeader(file->fpa, phu, format);
+        pmFPAview *view = pmFPAAddSourceFromHeader (file->fpa, phu, format);
 
         // XXX is this the correct psMD to save the filename?
@@ -957 +989 @@
 bool pmFPAfileCopyView (pmFPA *out, pmFPA *in, const pmFPAview *view)
 {
-    // XXX how to we transmit the required binning factor?
     // pmFPAWrite takes care of all PHUs as needed
     if (view->chip == -1) {
@@ -985 +1016 @@
     return false;
 
-    // XXX disable readout write for now
-    # if (0)
-
-        if (view->readout >= cell->readouts->n) {
-            return false;
-        }
-    pmReadout *readout = cell->readouts->data[view->readout];
-
-    if (view->nRows == 0) {
-        pmReadoutWrite (readout, fits, NULL, NULL);
-    } else {
-        pmReadoutWriteSegment (readout, fits, view->nRows, view->iRows, NULL, NULL);
-    }
-    return true;
-    # endif
-}
-
-// need to distinguish between a "copy structure" and a "copy data"
-// need to allow for binning when performing the "copy structure"
+    // XXX add readout / segment equivalents
+}
+
+// given an already-opened fits file, write the components corresponding
+// to the specified view
+bool pmFPAfileCopyStructureView (pmFPA *out, pmFPA *in, psMetadata *format, int xBin, int yBin, const pmFPAview *view)
+{
+    // pmFPAWrite takes care of all PHUs as needed
+    if (view->chip == -1) {
+        pmFPAAddSourceFromView (out, view, format);
+        pmFPACopyStructure (out, in, xBin, yBin);
+        return true;
+    }
+    if (view->chip >= in->chips->n) {
+        return false;
+    }
+    pmChip *inChip = in->chips->data[view->chip];
+    pmChip *outChip = out->chips->data[view->chip];
+
+    if (view->cell == -1) {
+        pmFPAAddSourceFromView (out, view, format);
+        pmChipCopyStructure (outChip, inChip, xBin, yBin);
+        return true;
+    }
+    if (view->cell >= inChip->cells->n) {
+        return false;
+    }
+    pmCell *inCell = inChip->cells->data[view->cell];
+    pmCell *outCell = outChip->cells->data[view->cell];
+
+    if (view->readout == -1) {
+        pmFPAAddSourceFromView (out, view, format);
+        pmCellCopyStructure (outCell, inCell, xBin, yBin);
+        return true;
+    }
+    return false;
+
+    // XXX add readout / segment equivalents
+}