Changeset 23666

Timestamp:

Apr 1, 2009, 4:47:30 PM (17 years ago)

Author:

eugene

Message:

add INST_MAG to output header, fix ghost count

File:

• trunk/psastro/src/psastroExtractStars.c (modified) (5 diffs)

trunk/psastro/src/psastroExtractStars.c

@@ -76 +76 @@
 
     int nExt = 0;
+    int nExtGhost = 0;
 
     // open the output file handle: we are just saving a series of extensions to this file
@@ -157 +158 @@
                     psMetadataAddF32 (header, PS_LIST_TAIL, "FPA_Y",   0, "focal-plane coordinate", ref->FP->y);
                     psMetadataAddF32 (header, PS_LIST_TAIL, "MAG",     0, "magnitude",              ref->Mag);
+                    psMetadataAddF32 (header, PS_LIST_TAIL, "INST_MAG",0, "instrumental magnitude", ref->Mag - MagOffset);
 
                     snprintf (extname, 80, "extname.%05d", nExt);
@@ -175 +177 @@
 
                     pmReadout *inReadout = pmFPAviewThisReadout (view, input->fpa);
-                    psImage *subraster = psastroExtractStar (inReadout->image, ghost->chip->x, ghost->chip->y, 200.0, 200.0);
+                    psImage *subraster = psastroExtractStar (inReadout->image, ghost->chip->x, ghost->chip->y, 400.0, 400.0);
                     if (!subraster) continue;
                     
@@ -187 +189 @@
                     psMetadataAddF32 (header, PS_LIST_TAIL, "FPA_Y0",  0, "star focal-plane coordinate", ghost->TP->y); // note : over-loaded value 
                     psMetadataAddF32 (header, PS_LIST_TAIL, "MAG",     0, "magnitude",              ghost->Mag);
-
-                    snprintf (extname, 80, "extname.%05d", nExt);
-                    psFitsWriteImage (outStars, header, subraster, 0, extname);
-                    nExt ++;
+                    psMetadataAddF32 (header, PS_LIST_TAIL, "INST_MAG",0, "instrumental magnitude", ghost->Mag - MagOffset);
+
+                    snprintf (extname, 80, "extname.%05d", nExtGhost);
+                    psFitsWriteImage (outGhosts, header, subraster, 0, extname);
+                    nExtGhost ++;
                     
                     psFree (header);
@@ -204 +207 @@
 
     psFitsClose (outStars);
+    psFitsClose (outGhosts);
     psFree (view);
     return true;