Changeset 15112 for branches/pap_branch_070920/psModules/src/camera/pmFPAfileDefine.c

Timestamp:

Sep 29, 2007, 12:01:01 PM (19 years ago)

Author:

Paul Price

Message:

Enabling compression when writing out. Introduces a new field in the FILERULES for output files (COMPRESSION), which provides a keyword which is looked up in COMPRESSION(METADATA) in the camera configuration. That then provides BITPIX,TILE.X,TILE.Y,TILE.Z,NOISE,HSCALE,HSMOOTH to set the compression and bits per pixel. Enabled compression, but not bitpix yet.

File:

• branches/pap_branch_070920/psModules/src/camera/pmFPAfileDefine.c (modified) (2 diffs)

branches/pap_branch_070920/psModules/src/camera/pmFPAfileDefine.c

@@ -78 +78 @@
 }
 
+// Parse an option from a metadata, returning the appropriate integer value
+static int parseOptionInt(const psMetadata *md, // Metadata containing the option
+                          const char *name, // Option name
+                          const char *source, // Description of source, for warning messages
+                          int defaultValue // Default value
+                          )
+{
+    psMetadataItem *item = psMetadataLookup(md, name); // Item with the value of interest
+    if (!item) {
+        psWarning("Unable to find value for %s in %s --- set to %d.", name, source, defaultValue);
+        return defaultValue;
+    }
+    int value = psMetadataItemParseS32(item); // Value of interst
+    if (value < 0) {
+        psWarning("Bad value for %s in %s (%d) --- set to %d.", name, source, value, defaultValue);
+        return defaultValue;
+    }
+    return value;
+}
+
 
 // define an input-type pmFPAfile, bind to the optional fpa if supplied
@@ -250 +270 @@
         file->fpa = pmFPAConstruct(file->camera);
     }
+
+    // Get compression scheme
+    const char *compName = psMetadataLookupStr(&status, data, "COMPRESSION"); // Name of compression scheme
+    if (compName && strcasecmp(compName, "NONE") != 0) {
+        psMetadata *compSchemes = psMetadataLookupMetadata(&status, camera, "COMPRESSION"); // Comp. schemes
+        if (!compSchemes) {
+            psWarning("Unable to find COMPRESSION in camera configuration --- compression disabled.");
+            goto COMPRESSION_DONE;
+        }
+        psMetadata *compression = psMetadataLookupMetadata(NULL, compSchemes, compName); // Compression info
+        if (!compression) {
+            psWarning("Unable to find %s in COMPRESSION in camera configuration --- compression disabled.",
+                      compName);
+            goto COMPRESSION_DONE;
+        }
+        const char *typeString = psMetadataLookupStr(NULL, compression, "COMPRESSION"); // Compression type
+        if (!typeString || strlen(typeString)) {
+            psWarning("Can't find COMPRESSION in compression scheme %s --- compression disabled.", compName);
+            goto COMPRESSION_DONE;
+        }
+        psFitsCompressionType type = psFitsCompressionTypeFromString(typeString); // Compression type enum
+
+        psString source = NULL;         // Source of options
+        psStringAppend(&source, "%s in COMPRESSION in camera %s", compName, cameraName);
+        psVector *tile = psVectorAlloc(3, PS_TYPE_S32); // Tile sizes
+        file->bitpix = parseOptionInt(compression, "BITPIX", source, 0); // Bits per pixel
+        tile->data.S32[0] = parseOptionInt(compression, "TILE.X", source, 0); // Tiling in x
+        tile->data.S32[1] = parseOptionInt(compression, "TILE.Y", source, 1); // Tiling in y
+        tile->data.S32[2] = parseOptionInt(compression, "TILE.Z", source, 1); // Tiling in z
+        int noise = parseOptionInt(compression, "NOISE", source, 16); // Noise bits
+        int hscale = parseOptionInt(compression, "HSCALE", source, 0); // Scaling for HCOMPRESS
+        int hsmooth = parseOptionInt(compression, "HSMOOTH", source, 0); // Smoothing for HCOMPRESS
+        psFree(source);
+
+        file->compression = psFitsCompressionAlloc(type, tile, noise, hscale, hsmooth);
+        psFree(tile);
+    }
+COMPRESSION_DONE:
 
     file->fileLevel = pmFPAPHULevel(format);