Changeset 19270 for trunk/psLib/src/fits/psFitsImage.c

Timestamp:

Aug 28, 2008, 3:01:39 PM (18 years ago)

Author:

Paul Price

Message:

The bad pixels were coming through in ds9 as the maximum value. According to
​http://heasarc.gsfc.nasa.gov/docs/heasarc/fits/compress/compress_image.html,
the keyword is "BLANK" when BITPIX > 0, and "ZBLANK" when BITPIX < 0.
Since we always do our own quantisation, the correct keyword is always
"BLANK" instead of "ZBLANK".

File:

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

trunk/psLib/src/fits/psFitsImage.c

@@ -7 +7 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.37 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-08-02 02:50:41 $
+ *  @version $Revision: 1.38 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-08-29 01:01:39 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -624 +624 @@
 
     if (blank != 0 && bitPix > 0) {
-        // Some quantisation has taken place --- record the blank ("magic") pixel value
-        char *keyword = (psFitsCompressionGetType(fits) != PS_FITS_COMPRESS_NONE &&
-                         (!fits->options || fits->options->conventions.compression)) ?
-            "ZBLANK" : "BLANK";         // Keyword for recording blank pixel value
-        fits_write_key_lng(fits->fd, keyword, blank, "Value for undefined pixels", &status);
+        // Some quantisation has taken place --- record the blank ("magic") pixel value.
+        //
+        // According to http://heasarc.gsfc.nasa.gov/docs/heasarc/fits/compress/compress_image.html,
+        // the keyword is "BLANK" when BITPIX > 0, and "ZBLANK" when BITPIX < 0.  Since we do our
+        // own quantisation, the correct keyword is always "BLANK" instead of "ZBLANK".
+        fits_write_key_lng(fits->fd, "BLANK", blank, "Value for undefined pixels", &status);
         fits_set_imgnull(fits->fd, blank, &status);
         if (psFitsError(status, true, "Could not write BLANK header to file.")) {