Changeset 26035 for trunk/psModules/src/imcombine/pmSubtractionStamps.h

Timestamp:

Nov 4, 2009, 4:52:05 PM (17 years ago)

Author:

Paul Price

Message:

Soften errors, both in the least-squares matrix/vector construction and in the rejection stage. Softening in the rejection stage seems to be very important so that bright stamps aren't rejected (after which the solution can go crazy). Now carrying around a weight image = 1/variance instead of the variance image. This can be softened, saving many divisions in the matrix/vector construction.

File:

• trunk/psModules/src/imcombine/pmSubtractionStamps.h (modified) (7 diffs)

trunk/psModules/src/imcombine/pmSubtractionStamps.h

@@ -24 +24 @@
     psArray *flux;                      ///< Fluxes for possible stamps (or NULL)
     int footprint;                      ///< Half-size of stamps
+    float sysErr;                       ///< Systematic error
 } pmSubtractionStampList;
 
@@ -31 +32 @@
                                                     const psRegion *region, // Region for stamps, or NULL
                                                     int footprint, // Half-size of stamps
-                                                    float spacing // Rough average spacing between stamps
+                                                    float spacing, // Rough average spacing between stamps
+                                                    float sysErr  // Relative systematic error or NAN
     );
 
@@ -68 +70 @@
     psKernel *image1;                   ///< Reference image postage stamp
     psKernel *image2;                   ///< Input image postage stamp
-    psKernel *variance;                 ///< Variance image postage stamp, or NULL
+    psKernel *weight;                   ///< Weight image (1/variance) postage stamp, or NULL
     psArray *convolutions1;             ///< Convolutions of image 1 for each kernel component, or NULL
     psArray *convolutions2;             ///< Convolutions of image 2 for each kernel component, or NULL
@@ -91 +93 @@
                                                 int footprint, ///< Half-size for stamps
                                                 float spacing, ///< Rough spacing for stamps
+                                                float sysErr,  ///< Relative systematic error in images
                                                 pmSubtractionMode mode ///< Mode for subtraction
     );
@@ -103 +106 @@
                                                int footprint, ///< Half-size for stamps
                                                float spacing, ///< Rough spacing for stamps
+                                               float sysErr,  ///< Systematic error in images
                                                pmSubtractionMode mode ///< Mode for subtraction
     );
@@ -115 +119 @@
     int footprint,                      ///< Half-size for stamps
     float spacing,                      ///< Rough spacing for stamps
+    float sysErr,                       ///< Systematic error in images
     pmSubtractionMode mode              ///< Mode for subtraction
     );
@@ -127 +132 @@
     int footprint,                      ///< Half-size for stamps
     float spacing,                      ///< Rough spacing for stamps
+    float sysErr,                       ///< Systematic error in images
     pmSubtractionMode mode              ///< Mode for subtraction
     );