Changeset 31618 for branches/czw_branch/20110406/psLib/src/fits/psFitsImage.c

Timestamp:

Jun 8, 2011, 2:41:13 PM (15 years ago)

Author:

watersc1

Message:

compilable version of updated stack code (doesn't correctly reject) and asinh scaling (completely untested).

File:

• branches/czw_branch/20110406/psLib/src/fits/psFitsImage.c (modified) (10 diffs)

branches/czw_branch/20110406/psLib/src/fits/psFitsImage.c

@@ -55 +55 @@
     int psDatatype;                     // psLib data type
     bool is_logscaled;                  // is this image log scaled using BOFFSET?
+    bool is_asinh;                      // is this image asinh scaled using BSOFTEN?
 } p_psFitsReadInfo;
 
@@ -129 +130 @@
 
     // Check scale and zero
-    double bscale = 0.0, bzero = 0.0, boffset = NAN;    // Scale and zero point
+    double bscale = 0.0, bzero = 0.0, boffset = NAN, bsoften = NAN;    // Scale and zero point
     if (fits_read_key_dbl(fits->fd, "BSCALE", &bscale, NULL, &status) && status != KEY_NO_EXIST) {
         psFitsError(status, true, "Unable to read header.");
@@ -153 +154 @@
       info->is_logscaled = false;
     }
+    status = 0;
+    if (fits_read_key_dbl(fits->fd, "BSOFTEN", &bsoften, NULL, &status) && status != KEY_NO_EXIST) {
+        psFitsError(status, true, "Unable to read header.");
+        goto bad;
+    }
+    if (status == KEY_NO_EXIST) {
+      info->is_asinh = false;
+    }
+    else if (isfinite(bsoften)) {
+      info->is_asinh = true;
+    }
+    else {
+      info->is_asinh = false;
+    }
+
     status = 0;
 
@@ -262 +278 @@
                                           double *bzero, // Zero point applied
                                           double *boffset, // Log offset applied
+                                          double *bsoften, // asinh offset applied
                                           long *blank, // Blank value (integer data)
                                           psFitsFloat *floatType, // Type of custom floating-point
@@ -275 +292 @@
     psAssert(bzero, "impossible");
     psAssert(boffset, "impossible");
+    psAssert(bsoften, "impossible");
     psAssert(floatType, "impossible");
     psAssert(fits, "impossible");
@@ -322 +340 @@
         if (newScaleZero) {
             // Choose an appropriate BSCALE and BZERO
-          if (!psFitsScaleDetermine(bscale, bzero, boffset, blank, image, mask, maskVal, fits)) {
+          if (!psFitsScaleDetermine(bscale, bzero, boffset, bsoften, blank, image, mask, maskVal, fits)) {
                 // We can't have the write dying for this reason --- try to save it somehow!
                 psWarning("Unable to determine BSCALE and BZERO for image --- refusing to quantise.");
@@ -348 +366 @@
         }
 
-        return psFitsScaleForDisk(image, fits, *bscale, *bzero, *boffset, rng);
+        return psFitsScaleForDisk(image, fits, *bscale, *bzero, *boffset, *bsoften, rng);
     }
 
@@ -487 +505 @@
       status = 0;
       fits_read_key_dbl(fits->fd, "BOFFSET", &boffset, NULL, &status);
-      psImage *newImage = psFitsScaleFromDisk(outImage,boffset);
+      psImage *newImage = psFitsScaleFromDisk(outImage,boffset,0.0);
+      psFree (outImage);
+      outImage = newImage;
+    }
+    else if (info->is_asinh) {
+      double bsoften;
+      int status;
+      status = 0;
+      fits_read_key_dbl(fits->fd, "BSOFTEN", &bsoften, NULL, &status);
+      psImage *newImage = psFitsScaleFromDisk(outImage,0.0,bsoften);
       psFree (outImage);
       outImage = newImage;
@@ -602 +629 @@
     double bscale = NAN, bzero = NAN;   // Scale and zero point to put in header
     double boffset = NAN;               // Log offset to put into header.
+    double bsoften = NAN;               // Asinh softening parameter to put into header
     long blank = 0;                     // Blank (undefined) value for image
     psFitsFloat floatType;              // Custom floating-point convention type
-    psImage *diskImage = imageToDiskRepresentation(&bscale, &bzero, &boffset, &blank, &floatType, fits, image,
+    psImage *diskImage = imageToDiskRepresentation(&bscale, &bzero, &boffset, &bsoften, &blank, &floatType, fits, image,
                                                    mask, maskVal, NULL, true); // Image to write out
     if (!diskImage) {
@@ -831 +859 @@
     double bscale = NAN, bzero = NAN;   // Scale and zero point to put in header
     double boffset = NAN;               // Log offset to put in header
+    double bsoften = NAN;               // Asinh softening parameter to put in header
     long blank = 0;                     // Blank (undefined) value for image
     psFitsFloat floatType;              // Custom floating-point convention type
-    psImage *diskImage = imageToDiskRepresentation(&bscale, &bzero, &boffset, &blank, &floatType, fits, input,
+    psImage *diskImage = imageToDiskRepresentation(&bscale, &bzero, &boffset, &bsoften, &blank, &floatType, fits, input,
                                                    mask, maskVal, NULL, false); // Image to write out
     if (!diskImage) {