Changeset 41810

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

Location:

trunk/psModules/src/astrom

Files:

• pmAstrometryDistortion.c (modified) (1 diff)
• pmAstrometryModel.c (modified) (3 diffs)
• pmAstrometryWCS.c (modified) (4 diffs)

trunk/psModules/src/astrom/pmAstrometryDistortion.c

@@ -347 +347 @@
     // as of r40806, psPlaneTransformInvert supplies the extra order (if non-linear)
     psFree (fpa->fromTPA);
-    fpa->fromTPA = psPlaneTransformInvert(NULL, fpa->toTPA, *region, 50, 4);
+    fpa->fromTPA = psPlaneTransformInvert(NULL, fpa->toTPA, *region, 100, 6);
     psFree (region);
 

trunk/psModules/src/astrom/pmAstrometryModel.c

@@ -525 +525 @@
         // 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);
@@ -592 +592 @@
     // as of r40806, psPlaneTransformInvert supplies the extra order (if non-linear)
     psFree (file->fpa->fromTPA);
-    file->fpa->fromTPA = psPlaneTransformInvert(NULL, file->fpa->toTPA, *region, 50, 4);
+    file->fpa->fromTPA = psPlaneTransformInvert(NULL, file->fpa->toTPA, *region, 100, 6);
 
     psFree (model);
@@ -745 +745 @@
 
     psFree (file->fpa->fromTPA);
-    file->fpa->fromTPA = psPlaneTransformInvert(NULL, file->fpa->toTPA, *region, 50, 4);
+    file->fpa->fromTPA = psPlaneTransformInvert(NULL, file->fpa->toTPA, *region,100, 10);
 
     psFree (region);

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);