Changeset 30873 for trunk/psastro/src/psastroMaskUpdates.c

Timestamp:

Mar 10, 2011, 5:58:10 PM (15 years ago)

Author:

watersc1

Message:

diffraction spike width scaling

File:

• trunk/psastro/src/psastroMaskUpdates.c (modified) (2 diffs)

trunk/psastro/src/psastroMaskUpdates.c

@@ -137 +137 @@
     double REFSTAR_MASK_SATSPIKE_MAG_MAX   = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_SATSPIKE_MAG_MAX");
     double REFSTAR_MASK_SATSPIKE_WIDTH     = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_SATSPIKE_WIDTH");
+    double REFSTAR_MASK_SATSPIKE_L0        = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_SATSPIKE_L0");
+    double REFSTAR_MASK_SATSPIKE_WIDTH_SLOPE = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_SATSPIKE_WIDTH_SLOPE");
     double REFSTAR_MASK_SATSPIKE_OFFSET    = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_SATSPIKE_OFFSET");
     double REFSTAR_MASK_BLEED_MAG_MAX      = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_BLEED_MAG_MAX");
@@ -325 +327 @@
                           spikeLength = 0.0;
                         }
-                        // XXX we can make the width depend on the spike as well...
+                        // XXX we can make the width depend on the spike as well...
+                        // CZW: 2011-02-18 it does so now. It's a weak function (slope = 0.01), but it should help
+                        if (spikeLength > REFSTAR_MASK_SATSPIKE_L0) {
+                          spikeWidth += 0.5 * (spikeLength - REFSTAR_MASK_SATSPIKE_L0) * REFSTAR_MASK_SATSPIKE_WIDTH_SLOPE;
+                        }
+
                         // The length should also be a function of the image background level