Changeset 18208

Timestamp:

Jun 19, 2008, 4:50:01 PM (18 years ago)

Author:

eugene

Message:

allow case where bscale/bzero cannot be determined

File:

• trunk/psLib/src/fits/psFitsImage.c (modified) (6 diffs)

trunk/psLib/src/fits/psFitsImage.c

@@ -7 +7 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.30 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-06-17 02:46:05 $
+ *  @version $Revision: 1.31 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-06-20 02:50:01 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -534 +534 @@
     }
 
+    bool useRequestedScale = true;
     if (!isfinite(bscale) || !isfinite(bzero)) {
         // Couldn't scale, so don't compress.  Save compression parameters for later
+        useRequestedScale = false;
         compress = psFitsCompressionGet(fits);
         if (!psFitsSetCompression(fits, PS_FITS_COMPRESS_NONE, NULL, 0, 0, 0)) {
@@ -556 +558 @@
         bzero = cfitsioBzero;
     }
+
     psFitsOptions *options = fits->options; // FITS I/O options
-    psAssert(!options || bitPix == options->bitpix || options->bitpix == 0, "impossible");
+    psAssert(!useRequestedScale || !options || bitPix == options->bitpix || options->bitpix == 0, "impossible");
 
     int naxis = 3;                      // Number of axes
@@ -603 +606 @@
     fits_set_bscale(fits->fd, 1.0, cfitsioBzero, &status);
 
-    if (bscale != 0.0) {
+    if (isfinite(bzero) && isfinite(bscale) && (bscale != 0.0)) {
         fits_write_key_dbl(fits->fd, "BZERO", bzero, 12,
                            "Scaling: TRUE = BZERO + BSCALE * DISK", &status);
@@ -693 +696 @@
     }
 
+    bool useRequestedScale = true;
     if (!isfinite(bscale) || !isfinite(bzero)) {
         // Couldn't scale, so don't compress.  Save compression parameters for later
+        useRequestedScale = false;
         compress = psFitsCompressionGet(fits);
         if (!psFitsSetCompression(fits, PS_FITS_COMPRESS_NONE, NULL, 0, 0, 0)) {
@@ -716 +721 @@
     }
     psFitsOptions *options = fits->options; // FITS I/O options
-    psAssert(!options || bitPix == options->bitpix || options->bitpix == 0, "impossible");
+    psAssert(!useRequestedScale || !options || bitPix == options->bitpix || options->bitpix == 0, "impossible");
 
     // Check to see if the HDU has the same datatype