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

Timestamp:

Apr 10, 2006, 10:22:42 AM (20 years ago)

Author:

magnier

Message:

adding concepts to copy from fpa to fpa

File:

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

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

@@ -7 +7 @@
 #include "pmFPAview.h"
 #include "pmFPAfile.h"
+#include "pmFPACopy.h"
 #include "pmFPARead.h"
 #include "pmFPAWrite.h"
@@ -18 +19 @@
 
     psFree (file->fpa);
+    psFree (file->readout);
     psFree (file->names);
 
@@ -46 +48 @@
 
     file->phu = NULL;
+    file->readout = NULL;
     file->header = NULL;
 
@@ -94 +97 @@
 
     pmFPAfile *file = pmFPAfileAlloc ();
+
+    // save the name of this pmFPAfile
+    file->name = psStringCopy (name);
+
     file->filerule = psMemIncrRefCounter(psMetadataLookupStr (&status, data, "FILENAME.RULE"));
     file->filextra = psMemIncrRefCounter(psMetadataLookupStr (&status, data, "FILENAME.XTRA"));
@@ -213 +220 @@
 
     if (file->mode == PM_FPA_MODE_NONE) {
+        return false;
+    }
+    if (file->mode == PM_FPA_MODE_INTERNAL) {
         return false;
     }
@@ -350 +360 @@
         pmFPAviewWriteObjects (view, file);
         psTrace ("pmFPAfile", 5, "wrote object %s (fpa: %p)\n", file->filename, file->fpa);
+        break;
+
+    default:
+        fprintf (stderr, "warning: type mismatch\n");
+        return false;
+    }
+    return true;
+}
+
+// create the data elements (headers, images) appropriate for this view
+bool pmFPAfileCreate (pmFPAfile *file, const pmFPAview *view)
+{
+    if (file->mode != PM_FPA_MODE_WRITE)
+        return false;
+
+    // get the current depth
+    pmFPAdepth depth = pmFPAviewDepth (view);
+
+    // do we need to write this file?
+    if (depth != file->dataDepth)
+        return false;
+
+    // XXX is this a sufficient check to avoid creating elements?
+    if (file->src == NULL)
+        return false;
+
+    switch (file->type) {
+    case PM_FPA_FILE_IMAGE:
+        /* create a PHU for thie file, if it does not exist */
+        pmFPAfileCopyView (file->src, file->fpa, view);
+        psTrace ("pmFPAfile", 5, "created fpa data elements for %s (fpa: %p)\n", file->filename, file->fpa);
+        break;
+
+    case PM_FPA_FILE_SX:
+    case PM_FPA_FILE_RAW:
+    case PM_FPA_FILE_OBJ:
+    case PM_FPA_FILE_CMP:
+    case PM_FPA_FILE_CMF:
         break;
 
@@ -409 +457 @@
         if (place == PM_FPA_BEFORE) {
             pmFPAfileRead (file, view);
+            pmFPAfileCreate (file, view);
         } else {
             pmFPAfileWrite (file, view);
@@ -418 +467 @@
 }
 
+// create a file with the given name, assign it type "INTERNAL", and supply it with an image
+// of the requested dimensions. (image only, mask and weight are ignored)
+pmReadout *pmFPAfileCreateInternal (psMetadata *files, char *name, int Nx, int Ny, int type)
+{
+    pmReadout *readout = pmReadoutAlloc (NULL);
+    readout->image = psImageAlloc (Nx, Ny, type);
+
+    // I want an image from the
+    pmFPAfile *file = pmFPAfileAlloc();
+    file->mode = PM_FPA_MODE_INTERNAL;
+
+    file->readout = readout;
+    psMetadataAddPtr (files, PS_LIST_TAIL, name, PS_DATA_UNKNOWN, "", file);
+    psFree (file);
+    // we free this copy of file, but 'files' still has a copy
+
+    return (readout);
+}
+
+bool pmFPAfileDropInternal (psMetadata *files, char *name)
+{
+    bool status;
+
+    pmFPAfile *file = psMetadataLookupPtr (&status, files, name);
+    if (file == NULL)
+        return false;
+
+    if (file->mode != PM_FPA_MODE_INTERNAL)
+        return false;
+
+    psMetadataRemoveKey (files, name);
+    return true;
+}
+
 // select the readout from the named pmFPAfile; if the named file does not exist,
-// create an image of the requested dimensions (save this image and file?)
-psImage *pmFPAfileReadoutImage (psMetadata *files, const pmFPAview *view, char *name, int Nx, int Ny, int type)
+pmReadout *pmFPAfileThisReadout (psMetadata *files, const pmFPAview *view, const char *name)
 {
     bool status;
-    psImage *image = NULL;
-
-    // I want an image from the
+
     pmFPAfile *file = psMetadataLookupPtr (&status, files, name);
-    if (file == NULL) {
-        // use Nx or Ny == 0 to avoid creating an image
-        if (Nx*Ny <= 0) {
-            return NULL;
-        }
-        image = psImageAlloc (Nx, Ny, type);
-        return (image);
-    }
-
-    // make an fpa consistent with this view
-    // save the image in the fpa
-    // save the fpa in the file
-
-    // inconsistent result: return NULL
+    if (file == NULL)
+        return NULL;
+
+    // internal files have the readout as a separate element:
+    if (file->mode == PM_FPA_MODE_INTERNAL) {
+        return file->readout;
+    }
+
     pmReadout *readout = pmFPAviewThisReadout (view, file->fpa);
-    if (readout == NULL) {
-        return NULL;
-    }
-
-    // resize the existing readout??? this is wrong...
-    // need to allocate it if it does not exist...
-    image = psImageRecycle (readout->image, Nx, Ny, type);
-    psMemIncrRefCounter (image);
-    return (image);
+    return readout;
 }
 
@@ -649 +715 @@
 }
 
+// XXX this this function through, then finish
+# if 0
 // look for the given name on the argument list.
 // returns the file (a view to the one saved on config->files)
@@ -734 +802 @@
     return file;
 }
+# endif
 
 // look for the given name on the argument list.
@@ -884 +953 @@
 }
 
+// given an already-opened fits file, write the components corresponding
+// to the specified view
+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) {
+        pmFPACopy (out, in);
+        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) {
+        pmChipCopy (outChip, inChip);
+        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) {
+        pmCellCopy (outCell, inCell);
+        return true;
+    }
+    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"