Changeset 10815

Timestamp:

Dec 20, 2006, 6:44:45 AM (20 years ago)

Author:

eugene

Message:

mosaic astrometry work

Location:

trunk/psastro/src

Files:

• psastroMosaicAstrom.c (modified) (1 diff)
• psastroMosaicRescaleChips.c (modified) (2 diffs)

trunk/psastro/src/psastroMosaicAstrom.c

@@ -28 +28 @@
 
     // XXX no input distortion model yet; existing fpa distortion is identity: replace it
+    // XXX make this a test if these are NULL; the load step can NULL them if nothing is loaded
+    // XXX if we have an input distortion model, the gradient is measuring the error rel to that model
+    // XXX how does this couple to the individual chip fits?
     psFree (fpa->toTPA);
     psFree (fpa->fromTPA);

trunk/psastro/src/psastroMosaicRescaleChips.c

@@ -1 +1 @@
 # include "psastro.h"
+
+// XXX what is this doing? 
+// shouldn't it be doing toFPA -> fromFPA
 
 bool psastroMosaicRescaleChips (pmFPA *fpa) {
@@ -7 +10 @@
 
     // this loop selects the matched stars for all chips
+
+    // psRegion region = psMetadataLookupXXX (chip->concepts, "CHIP.TRIMSEC"); 
+    // psRegion region = psRegionSet (0, 4000, 0, 4000);
+    // chip->fromFPA = psPlaneTransformInvert(NULL, toFPA, region, 50);
+
     while ((chip = pmFPAviewNextChip (view, fpa, 1)) != NULL) {
         psTrace ("psastro", 4, "Chip %d: %x %x\n", view->chip, chip->file_exists, chip->process);
         if (!chip->process || !chip->file_exists) { continue; }
-        
-        psFree (chip->toFPA);
-        chip->toFPA = p_psPlaneTransformLinearInvert(chip->fromFPA);
+
+        // XXX make '50' be a fraction of the chip size?
+        chip->fromFPA = psPlaneTransformInvert(chip->fromFPA, chip->toFPA, region, 50);
     }
     psFree (view);