Changeset 39981 for trunk/psModules

Timestamp:

Feb 23, 2017, 4:22:35 PM (9 years ago)

Author:

watersc1

Message:

Add code to allow the number of additional polynomial terms in the WCS solution to be modified. This should allow updates to look like the original run.

File:

• trunk/psModules/src/astrom/pmAstrometryWCS.c (modified) (1 diff)

trunk/psModules/src/astrom/pmAstrometryWCS.c

@@ -695 +695 @@
     // the region defines the FPA pixels covered by the tranformation
     psFree (fpa->fromTPA);
-    psPlaneTransform *myPT = psPlaneTransformAlloc(fpa->toTPA->x->nX+4, fpa->toTPA->x->nY+4);
+
+    int additional_orders = 4;  // This is the number of orders that should be added.
+    bool status = false;
+    int config_additional_orders = psMetadataLookupS32(&status,fpa->analysis, "ADDITIONAL_WCS_ORDERS");
+    if (status) { // If this keyword is found, then use that number of orders instead
+      additional_orders = config_additional_orders; 
+    }
+    
+    psPlaneTransform *myPT = psPlaneTransformAlloc(fpa->toTPA->x->nX+additional_orders, fpa->toTPA->x->nY+additional_orders);
     fpa->fromTPA = psPlaneTransformInvert(myPT, fpa->toTPA, region, 50);
     psFree (myPT);