Changeset 1444 for trunk/src

Timestamp:

10/28/08 21:09:24 (16 years ago)

Author:

Malte Marquarding

Message:

Ticket #110: I have added a fix for the FITS header now, This modifies CRVAL1, CTYPE1 to be read by CLASS

Location:

trunk/src

Files:

• STFITSImageWriter.cpp (modified) (2 diffs)
• STFITSImageWriter.h (modified) (1 diff)

trunk/src/STFITSImageWriter.cpp

@@ -30 +30 @@
 //#---------------------------------------------------------------------------
 
+#include <fitsio.h>
 #include <images/Images/TempImage.h>
 
@@ -217 +218 @@
 
 void STFITSImageWriter::classHackHeader(const String& filename) {
-  cout << "Not yet implemented" << endl;
+  int status = 0;
+  fitsfile *fptr;     
+  cout << "filename" << endl;
+  if( fits_open_file(&fptr, filename.c_str(), READWRITE, &status) ) 
+    throw AipsError("FCoudn't open fits file for CLASS modification");
+
+  if ( fits_update_key(fptr, TSTRING, "CTYPE1", (char *)"FREQ",
+                       NULL, &status) )
+    throw AipsError("Couldn't modify CTYPE1.");
+  float restf,refval,newfreq;
+  fits_read_key(fptr, TFLOAT, "CRVAL1", 
+                &refval, NULL, &status);
+  fits_read_key(fptr, TFLOAT, "RESTFREQ", 
+                &restf, NULL, &status);
+  newfreq = refval - restf;
+  if ( fits_update_key(fptr, TFLOAT, "CRVAL1", &newfreq,  NULL, &status) )
+    throw AipsError("Couldn't modify CRVAL1");
+  fits_close_file(fptr, &status);
+
 }

trunk/src/STFITSImageWriter.h

@@ -61 +61 @@
                                                    casa::Double lon, 
                                                    casa::Double lat);
+  // Apply hacks to fits file so it can be read by class
+  // modifications are
+  // CTYPE1 :  FREQ-XYZ -> FREQ
+  // CRVAL1 :  -> CRVAL1-RESTFREQ
   void classHackHeader(const casa::String& filename);