Changeset 41810 for trunk/psModules/src/astrom/pmAstrometryWCS.c

Timestamp:

Sep 14, 2021, 2:18:52 PM (5 years ago)

Author:

tdeboer

Message:

changes to order astrometry solution transformation, giving better reproduction of forward<->backward conversion

File:

• trunk/psModules/src/astrom/pmAstrometryWCS.c (modified) (4 diffs)

trunk/psModules/src/astrom/pmAstrometryWCS.c

@@ -608 +608 @@
     // as of r40806, psPlaneTransformInvert supplies the extra order (if non-linear)
     psFree (chip->fromFPA);
-    chip->fromFPA = psPlaneTransformInvert(NULL, chip->toFPA, *region, 50, 4);
+    chip->fromFPA = psPlaneTransformInvert(NULL, chip->toFPA, *region, 100, 6);
     psFree (region);
 
@@ -668 +668 @@
     // as of r40806, psPlaneTransformInvert supplies the extra order (if non-linear)
     psFree (chip->fromFPA);
-    chip->fromFPA = psPlaneTransformInvert(NULL, chip->toFPA, *region, 50, 4);
+    chip->fromFPA = psPlaneTransformInvert(NULL, chip->toFPA, *region, 100, 6);
     psFree (region);
 
@@ -706 +706 @@
       additional_orders = config_additional_orders; 
     }
-    fpa->fromTPA = psPlaneTransformInvert(NULL, fpa->toTPA, region, 50, additional_orders);
+    fpa->fromTPA = psPlaneTransformInvert(NULL, fpa->toTPA, region, 100, additional_orders);
     return true;
 }
@@ -1032 +1032 @@
 
     // NOTE: the extraOrders value (4) should be ignored since outToFPA is specified to be linear
-    psPlaneTransform *outFromFPA = psPlaneTransformInvert(NULL, outToFPA, *outputBounds, 50, 4);
+    psPlaneTransform *outFromFPA = psPlaneTransformInvert(NULL, outToFPA, *outputBounds, 100, 6);
     if (!outFromFPA) {
         psFree(outToFPA);