Changeset 36835 for trunk/pswarp/src/pswarpUpdateMetadata.c

Timestamp:

Jun 7, 2014, 6:38:38 AM (12 years ago)

Author:

eugene

Message:

add pswarp backwards transformations so we can map warp pixels to chip pixels

File:

• trunk/pswarp/src/pswarpUpdateMetadata.c (modified) (1 diff)

trunk/pswarp/src/pswarpUpdateMetadata.c

@@ -78 +78 @@
                 psFree(md5string);
                 psFree(headerName);
+
+                bool status;
+                psString keyword = NULL;
+                psString mapstring = NULL;
+
+                // we (should) have an array of chipnames on the analysis
+                psArray *chipnames = psMetadataLookupPtr(&status, readout->analysis, PSWARP_ANALYSIS_CHIPNAMES); 
+                for (int i = 0; chipnames && (i < chipnames->n); i++) {
+                  psStringAppend (&keyword, "SRC_%04d", i);
+                  psMetadataAddStr(hdu->header, PS_LIST_TAIL, keyword, PS_META_REPLACE, "input image", chipnames->data[i]);
+                  psFree (keyword);
+                }
+
+                // we (should) have an array of chipnames on the analysis
+                psArray *chipRegions = psMetadataLookupPtr(&status, readout->analysis, PSWARP_ANALYSIS_CHIPREGIONS); 
+                for (int i = 0; chipRegions && (i < chipRegions->n); i++) {
+                  psStringAppend (&keyword, "SEC_%04d", i);
+                  psMetadataAddStr(hdu->header, PS_LIST_TAIL, keyword, PS_META_REPLACE, "input image", chipRegions->data[i]);
+                  psFree (keyword);
+                }
+
+                // we (should) also have an array of backwards maps (warp->chip) for the above chips
+                psArray *backmaps  = psMetadataLookupPtr(&status, readout->analysis, PSWARP_ANALYSIS_BACKMAPS); // 
+                for (int i = 0; backmaps && (i < backmaps->n); i++) {
+                  psPlaneTransform *map = backmaps->data[i];              
+                  psStringAppend (&keyword, "MPX_%04d", i);
+                  psStringAppend (&mapstring, "[%8.2f,%8.4f,%8.4f]", map->x->coeff[0][0], map->x->coeff[1][0], map->x->coeff[0][1]);
+                  psMetadataAddStr(hdu->header, PS_LIST_TAIL, keyword, PS_META_REPLACE, "warp to image map", mapstring);
+                  psFree (keyword);
+                  psFree (mapstring);
+                  psStringAppend (&keyword, "MPY_%04d", i);
+                  psStringAppend (&mapstring, "[%8.2f,%8.4f,%8.4f]", map->y->coeff[0][0], map->y->coeff[1][0], map->y->coeff[0][1]);
+                  psMetadataAddStr(hdu->header, PS_LIST_TAIL, keyword, PS_META_REPLACE, "warp to image map", mapstring);
+                  psFree (keyword);
+                  psFree (mapstring);
+                }
             }