Changeset 30334 for branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionMatch.c

Timestamp:

Jan 21, 2011, 11:32:03 AM (15 years ago)

Author:

eugene

Message:

adjust scale if radial moment and fwhm disagree too much

File:

• branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionMatch.c (modified) (1 diff)

branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionMatch.c

@@ -690 +690 @@
             }
 
+            // check on the kernel scaling -- if the kron-based radial moments are very different, adjust to match them
+            { 
+                float fwhm1;
+                float fwhm2;
+                pmSubtractionGetFWHMs(&fwhm1, &fwhm2);
+                psAssert(isfinite(fwhm1), "fwhm 1 not set");
+                psAssert(isfinite(fwhm2), "fwhm 2 not set");
+
+                // XXX this is BAD: depends on the relationship below:
+                // stamps->normWindow1 = 2.75*R1;
+                // stamps->normWindow2 = 2.75*R2;
+                float radMoment1 = stamps->normWindow1 / 2.75;
+                float radMoment2 = stamps->normWindow2 / 2.75;
+
+                float maxFWHM = PS_MAX(fwhm1, fwhm2);
+                float maxRadial = PS_MAX(radMoment1, radMoment2);
+            
+                if (fabs(2.0*(maxFWHM - maxRadial)/(maxFWHM + maxRadial)) > 0.25) {
+
+                    float scale = maxRadial / maxFWHM;
+                    psLogMsg ("psModules.imcombine", PS_LOG_INFO, "Kron and FWHM scales are quite different, re-scale by %f to use Kron", scale);
+                    
+                    for (int i = 0; i < isisWidths->n; i++) {
+                        isisWidths->data.F32[i] *= scale;
+                    }
+                }
+            }
+
             // Define kernel basis functions
             if (optimum && (type == PM_SUBTRACTION_KERNEL_ISIS || type == PM_SUBTRACTION_KERNEL_GUNK)) {