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

Timestamp:

Mar 13, 2006, 4:21:07 PM (20 years ago)

Author:

magnier

Message:

updates to pmFPAview, pmFPAfile

File:

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

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

@@ -3 +3 @@
 #include "pmHDU.h"
 #include "pmFPA.h"
+#include "pmFPAConstruct.h"
 #include "pmFPAview.h"
 #include "pmFPAfile.h"
+#include "pmFPARead.h"
+#include "pmFPAWrite.h"
 #include "pmFPAviewObjectsIO.h"
 
@@ -40 +43 @@
 
 // XXX is a file allowed to be both READ and WRITE?
-pmFPAfile *pmFPAfileDefine (pmFPAview *view, psMetadata *camera, char *name)
+pmFPAfile *pmFPAfileDefine (psMetadata *files, psMetadata *camera, pmFPA *fpa, char *name)
 {
 
@@ -49 +52 @@
     psMetadata *data = psMetadataLookupPtr (&status, camera, name);
     if (data == NULL) {
-        return false;
+        return NULL;
     }
 
@@ -114 +117 @@
     }
 
+    if (fpa == NULL) {
+        fpa = pmFPAConstruct (camera);
+    }
+    file->fpa = fpa;
+
+    psMetadataAddPtr (files, PS_LIST_TAIL, name, PS_META_REPLACE, "", file);
     return (file);
 }
 
-bool pmFPAfileOpen (pmFPAfile *file, pmFPAview *view, psMetadata *camera)
+bool pmFPAfileOpen (pmFPAfile *file, pmFPAview *view)
 {
 
@@ -130 +139 @@
 
     // determine the file name
-    file->filename = pmFPAviewNameFromRule (file->filerule, view);
+    file->filename = pmFPAviewNameFromRule (file->filerule, view, file->fpa);
 
     // indirect filenames
     if (!strcasecmp (file->filename, "@FPAIO")) {
-        extra = pmFPAviewNameFromRule (file->filextra, view);
-        file->filename = psMetadataLookupStr (&status, view->IO, extra);
+        extra = pmFPAviewNameFromRule (file->filextra, view, file->fpa);
+        file->filename = psMetadataLookupStr (&status, file->names, extra);
     }
     if (!strcasecmp (file->filename, "@DETDB")) {
-        extra = pmFPAviewNameFromRule (file->filextra, view);
+        extra = pmFPAviewNameFromRule (file->filextra, view, file->fpa);
         // file->filename = pmDetrendSelect (extra);
     }
@@ -165 +174 @@
     switch (file->type) {
     case PM_FPA_FILE_IMAGE:
-        pmFPAviewReadFitsImage (view, file->fits);
+        pmFPAviewReadFitsImage (view, file);
         break;
 
@@ -216 +225 @@
 }
 
+bool pmFPAfileClose (pmFPAfile *file, pmFPAview *view)
+{
+    // is current depth == open depth?
+    pmFPAdepth depth = pmFPAviewDepth (view);
+    if (file->openDepth != depth) {
+        return false;
+    }
+
+    // check if we are actually open
+    psFitsClose (file->fits);
+    return true;
+}
+
 // attempt open and read for all pmFPAfiles in view->IO
-bool pmFPAfileIOCheck (pmFPAview *view, psMetadata *camera)
-{
-
+bool pmFPAfileReadChecks (psMetadata *files, pmFPAview *view)
+{
     // recipe override values (command-line options):
     psMetadataItem *item = NULL;
-    psMetadataIterator *iter = psMetadataIteratorAlloc (view->IO, PS_LIST_HEAD, NULL);
+    psMetadataIterator *iter = psMetadataIteratorAlloc (files, PS_LIST_HEAD, NULL);
     while ((item = psMetadataGetAndIncrement (iter)) != NULL) {
         pmFPAfile *file = item->data.V;
 
-        pmFPAfileOpen (file, view, camera);
+        pmFPAfileOpen (file, view);
         pmFPAfileRead (file, view);
     }
@@ -232 +253 @@
     return true;
 }
+
+// attempt open and read for all pmFPAfiles in view->IO
+bool pmFPAfileWriteChecks (psMetadata *files, pmFPAview *view)
+{
+    // recipe override values (command-line options):
+    psMetadataItem *item = NULL;
+    psMetadataIterator *iter = psMetadataIteratorAlloc (files, PS_LIST_HEAD, NULL);
+    while ((item = psMetadataGetAndIncrement (iter)) != NULL) {
+        pmFPAfile *file = item->data.V;
+
+        pmFPAfileWrite (file, view);
+        pmFPAfileClose (file, view);
+    }
+    psFree (iter);
+    return true;
+}
+
+psImage *pmFPAfileReadoutImage (psMetadata *files, pmFPAview *view, char *name, int Nx, int Ny, int type)
+{
+    bool status;
+    psImage *image = NULL;
+
+    // I want an image from the
+    pmFPA *fpa = psMetadataLookupPtr (&status, files, name);
+    if (fpa == NULL) {
+        // use Nx or Ny == 0 to avoid creating an image
+        if (Nx*Ny <= 0) {
+            return NULL;
+        }
+        image = psImageAlloc (Nx, Ny, type);
+        return (image);
+    }
+
+    // inconsistent result: return NULL
+    pmReadout *readout = pmFPAviewThisReadout (view, 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);
+}
+
+# if (0)
+    // return the named pmFPAfile
+    pmFPAfile *pmFPAfileSelect (pmFPAview *view, char *name)
+{
+
+    // select the name from the camera config data
+    pmFPAfile *file = psMetadataLookupPtr (&status, view->IO, name);
+    if (file == NULL) {
+        return NULL;
+    }
+    return file;
+}
+
+bool pmFPAfileWriteChecks (pmFPAview *view, char *name)
+{
+
+    pmFPAfile *file = pmFPAfileSelect (view, name);
+    pmFPAfileWrite (file, view);
+    return true;
+}
+#endif
+
+// given an already-opened fits file, read the components corresponding
+// to the specified view
+bool pmFPAviewReadFitsImage (pmFPAview *view, pmFPAfile *file)
+{
+    pmFPA *fpa = file->fpa;
+    psFits *fits = file->fits;
+
+    if (view->chip == -1) {
+        pmFPARead (fpa, fits, NULL);
+        return true;
+    }
+
+    if (view->chip >= fpa->chips->n) {
+        return false;
+    }
+    pmChip *chip = fpa->chips->data[view->chip];
+
+    if (view->cell == -1) {
+        pmChipRead (chip, fits, NULL);
+        return true;
+    }
+
+    if (view->cell >= chip->cells->n) {
+        return false;
+    }
+    pmCell *cell = chip->cells->data[view->cell];
+
+    if (view->readout == -1) {
+        pmCellRead (cell, fits, NULL);
+        return true;
+    }
+
+    // XXX pmReadoutRead, pmReadoutReadSegement disabled for now
+    return false;
+
+    # if (0)
+
+        if (view->readout >= cell->readouts->n) {
+            return false;
+        }
+    pmReadout *readout = cell->readouts->data[view->readout];
+
+    if (view->nRows == 0) {
+        pmReadoutRead (readout, fits, NULL);
+    } else {
+        pmReadoutReadSegment (readout, fits, view->nRows, view->iRows, NULL, NULL);
+    }
+    return true;
+    # endif
+}
+
+// XXX image writes disabled for now
+# if (0)
+    // given an already-opened fits file, write the components corresponding
+    // to the specified view
+    bool pmFPAviewWriteFitsImage (pmFPAview *view, pmFPAfile *file)
+{
+
+    pmFPA *fpa = file->fpa;
+    psFits *fits = file->fits;
+
+    if (view->chip == -1) {
+        pmFPAWrite (fpa, fits, NULL, NULL);
+        return true;
+    }
+
+    if (view->chip >= fpa->chips->n) {
+        return false;
+    }
+    pmChip *chip = fpa->chips->data[view->chip];
+
+    if (view->cell == -1) {
+        pmChipWrite (chip, fits, NULL, NULL);
+        return true;
+    }
+
+    if (view->cell >= chip->cells->n) {
+        return false;
+    }
+    pmCell *cell = chip->cells->data[view->cell];
+
+    if (view->readout == -1) {
+        pmCellWrite (cell, fits, NULL, NULL);
+        return true;
+    }
+
+    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