Changeset 6663 for branches/rel10_ifa/psModules/src/astrom/pmHDU.c

Timestamp:

Mar 21, 2006, 1:52:54 PM (20 years ago)

Author:

Paul Price

Message:

Now works with GPC1 dummy OTA FITS file from Sidik; some other changes to organise things.

File:

• branches/rel10_ifa/psModules/src/astrom/pmHDU.c (modified) (3 diffs)

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

@@ -114 +114 @@
     assert(fits);
 
+    psTrace(__func__, 7, "Writing HDU %s\n", hdu->extname);
+
     if (hdu->images && hdu->table && !hdu->header) {
         psError(PS_ERR_IO, true, "Both image and table data provided in HDU, but no header --- "
@@ -120 +122 @@
     }
 
+    // Preserve the extension name, if it's the PHU
+    char *extname = hdu->extname;       // The name of the extension
+    if (strcasecmp(extname, "PHU") == 0) {
+        bool mdok = true;               // Status of MD lookup
+        extname = psMetadataLookupStr(&mdok, hdu->header, "EXTNAME");
+        if (!mdok || strlen(extname) == 0) {
+            extname = "";
+        }
+    }
+
+    if (!hdu->images && !psFitsWriteHeader(hdu->header, fits)) {
+        psError(PS_ERR_IO, false, "Unable to write header for extension %s\n", extname);
+    }
+
     #ifdef FITS_TABLES
     if (hdu->images && (!hdu->table || psFitsIsImage(hdu->header))) {
-        #endif
-        psFitsWriteImageCube(fits, hdu->header, hdu->images, hdu->extname);
+        psFitsWriteImageCube(fits, hdu->header, hdu->images, extname);
         return true;
-        #ifdef FITS_TABLES
-
     }
 
@@ -136 +149 @@
     psError(PS_ERR_IO, true, "No idea what this HDU consists of!\n");
     return false;
+    #else
+
+    if (hdu->images) {
+        psTrace(__func__, 9, "Writing pixels for %s\n", hdu->extname);
+        if (!psFitsWriteImageCube(fits, hdu->header, hdu->images, extname)) {
+            psError(PS_ERR_IO, false, "Unable to write image to extension %s\n", hdu->extname);
+            return false;
+        }
+    }
+    return true;
     #endif
 }
 
-
-
-#if 0
-
-
-
-
-// Read the PHU of a file
-static bool readPHU(pmCameraFile *file  // File for which to read the PHU
-                   )
-{
-    if (! file->phu) {
-        if (! file->fits) {
-            file->fits = psFitsOpen(file->name, "r");
-        }
-        (void)psFitsMoveExtNum(file->fits, 0, false);
-        file->phu = psFitsReadHeader(NULL, file->fits);
-
-        return true;
-    }
-
-    return false;
-}
-
-bool pmCameraFileFormat(pmCameraFile *file, // File to inspect
-                        psMetadata *camera // Camera description
-                       )
-{
-    if (! camera) {
-        return false;
-    }
-
-    readPHU(file);
-
-    bool mdok = true;                   // Result of MD lookup
-    psMetadata *formats = psMetadataLookupMD(&mdok, camera, "FORMATS"); // List of known formats
-    if (! mdok || ! formats) {
-        psLogMsg(__func__, PS_LOG_WARN, "Unable to find list of FORMATS in camera configuration!\n");
-        return false;
-    }
-
-    bool result = false;                // Result of search
-    psMetadataIterator *formatsIter = psMetadataIteratorAlloc(formats, NULL); // Iterator for formats
-    psMetadataItem *formatsItem = NULL; // Item from formats
-    while ((formatsItem = psMetadataGetAndIncrement(formatsIter))) {
-        if (formatsItem->type != PS_DATA_STRING) {
-            psLogMsg(__func__, PS_LOG_WARN, "Entry not of type STR found in camera FORMATS list.\n");
-            continue;
-        }
-        psTrace(__func__, 7, "Trying format file: %s\n", formatsItem->name);
-        int numBadLines = 0;            // Number of bad lines in file
-        psMetadata *format = psMetadataConfigParse(NULL, &numBadLines, formatsItem->data.V, true);
-        if (numBadLines > 0) {
-            psLogMsg(__func__, PS_LOG_WARN, "%d bad lines found in %s.\n", formatsItem->name);
-        }
-        if (pmConfigValidateCamera(format, file->phu)) {
-            file->format = format;
-            result = true;
-            psLogMsg(__func__, PS_LOG_INFO + 2, "File %s matches camera format file %s\n", file->name,
-                     formatsItem->data.V);
-            break;
-        }
-        psFree(format);
-    }
-    psFree(formatsIter);
-
-    return result;
-}
-
-
-bool pmConfigValidateCamera(const psMetadata *camera, const psMetadata *header);
-
-pmFPA *pmFPASourceFile(pmFPA *fpa,      // FPA to which to add
-                       pmFile *file,    // File to add
-                       const char *name, // Symbolic name for file, e.g., "INPUT", "OUTPUT", "TESTFILE"
-                       pmConfig *config // Configuration set
-                      )
-{
-    assert(file);
-    assert(name);
-
-    if (! file->fits) {
-        file->fits = psFitsOpen(file->name, "r");
-    }
-    if (! file->phu) {
-        (void)psFitsMoveExtNum(file->fits, 0, false);
-        file->phu = psFitsReadHeader(NULL, file->fits);
-    }
-
-    psMetadata *fileConfig = NULL;      // The file configuration
-    if (fpa && fpa->camera) {
-        fpa
-
-        bool pmConfigValidateCamera(const psMetadata *camera, const psMetadata *header);
-        #endif