Changeset 19087 for trunk/psModules/src/objects/pmSourceIO.c

Timestamp:

Aug 16, 2008, 12:36:18 PM (18 years ago)

Author:

eugene

Message:

replace NAXIS1,2 in header with IMNAXIS1,2

File:

• trunk/psModules/src/objects/pmSourceIO.c (modified) (2 diffs)

trunk/psModules/src/objects/pmSourceIO.c

@@ -3 +3 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.63 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-07-17 22:38:15 $
+ *  @version $Revision: 1.64 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-08-16 22:36:18 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -424 +424 @@
                 pmConfigConformHeader(hdu->header, file->format);
 
+                // psFitsWriteBlank strips out the NAXISn keywords, forcing CFITSIO to take care of them
+                // save NAXIS1,NAXIS2 as IMNAXIS1,IMNAXIS2
+                int Nx = psMetadataLookupS32 (&status, hdu->header, "IMNAXIS1");
+                if (!status) {
+                    Nx = psMetadataLookupS32 (&status, hdu->header, "NAXIS1");
+                    if (status) {
+                        psMetadataAddS32 (hdu->header, PS_LIST_TAIL, "IMNAXIS1", PS_META_REPLACE, "original image size (x-dir)", Nx);
+                    } else {
+                        fprintf (stderr, "header is missing NAXIS1!");
+                    }
+                }
+
+                int Ny = psMetadataLookupS32 (&status, hdu->header, "IMNAXIS2");
+                if (!status) {
+                    Ny = psMetadataLookupS32 (&status, hdu->header, "NAXIS2");
+                    if (status) {
+                        psMetadataAddS32 (hdu->header, PS_LIST_TAIL, "IMNAXIS2", PS_META_REPLACE, "original image size (y-dir)", Ny);
+                    } else {
+                        fprintf (stderr, "header is missing NAXIS2!");
+                    }
+                }
+
                 psFitsWriteBlank (file->fits, hdu->header, headname);
                 psTrace ("pmFPAfile", 5, "wrote ext head %s (type: %d)\n", file->filename, file->type);