Changeset 25710

Timestamp:

Oct 1, 2009, 9:32:14 AM (17 years ago)

Author:

bills

Message:

Put the RA and DEC of the center of the stamp into the header

File:

• trunk/pstamp/src/ppstampMakeStamp.c (modified) (7 diffs)

trunk/pstamp/src/ppstampMakeStamp.c

@@ -48 +48 @@
 }
 
-static bool copyMetadata(pmFPAfile *output, pmFPAfile *input, pmChip *inChip, ppstampOptions *options)
+static bool copyMetadata(pmFPAfile *output, pmFPAfile *input, pmChip *inChip, ppstampOptions *options, pmAstromObj *center)
 {
     pmChip    *outChip;
@@ -109 +109 @@
     }
 
+    psMetadataAddF64(outHDU->header, PS_LIST_TAIL, "RA_DEG", PS_META_REPLACE, "Right Ascension of stamp center", RAD_TO_DEG(center->sky->r));
+    psMetadataAddF64(outHDU->header, PS_LIST_TAIL, "DEC_DEG", PS_META_REPLACE, "Declination of stamp center", RAD_TO_DEG(center->sky->d));
+
     ppstampVersionMetadata(outHDU->header, options);
 
@@ -194 +197 @@
 
 static int makeStamp(pmConfig *config, ppstampOptions *options, pmFPAfile *input,
-                pmChip *inChip, pmFPAview *view)
+                pmChip *inChip, pmFPAview *view, pmAstromObj *center)
 {
     int status = false;
@@ -301 +304 @@
 
     if (status) {
-        status = copyMetadata(output, input, inChip, options);
+        status = copyMetadata(output, input, inChip, options, center);
     }
     return status ? PS_EXIT_SUCCESS : PS_EXIT_UNKNOWN_ERROR;
@@ -376 +379 @@
 {
     pmAstromObj *pt = pmAstromObjAlloc();
-
-    if (!options->roip.celestialCenter) {
-        // Center was specified in chip coordinates, transform to the sky and the TP
-        chipToSky(center, fpa, chip);
-    }
 
     // calculate the four corners of the bounding box in sky coordinates, translate them to
@@ -486 +484 @@
             center->chip->yErr = 0;
             onChip = true;
+            chipToSky(center, input->fpa, chip);
         }
     }
@@ -589 +588 @@
         case PPSTAMP_ON:
         case PPSTAMP_PARTIALLY_ON:
-            returnval = makeStamp(config, options, input, chip, view);
+            returnval = makeStamp(config, options, input, chip, view, center);
             allDone = true;
             foundOverlap = true;