Changeset 21087 for branches/eam_branch_20081230/psphot/src/psphotSourceSize.c

Timestamp:

Jan 6, 2009, 7:13:34 AM (18 years ago)

Author:

eugene

Message:

update to psImageMaskType and associated (psphot, psastro, ppImage, ppMerge, ppSim, ppStats, ppArith)

File:

• branches/eam_branch_20081230/psphot/src/psphotSourceSize.c (modified) (11 diffs)

branches/eam_branch_20081230/psphot/src/psphotSourceSize.c

@@ -3 +3 @@
 
 static float psphotModelContour(const psImage *image, const psImage *weight, const psImage *mask,
-                                psMaskType maskVal, const pmModel *model, float Ro);
-
-bool psphotMaskCosmicRay_Old (pmSource *source, psMaskType maskVal, psMaskType crMask);
-bool psphotMaskCosmicRay_New (psImage *mask, pmSource *source, psMaskType maskVal, psMaskType crMask);
+                                psImageMaskType maskVal, const pmModel *model, float Ro);
+
+bool psphotMaskCosmicRay_Old (pmSource *source, psImageMaskType maskVal, psImageMaskType crMask);
+bool psphotMaskCosmicRay_New (psImage *mask, pmSource *source, psImageMaskType maskVal, psImageMaskType crMask);
 
 // we need to call this function after sources have been fitted to the PSF model and
@@ -21 +21 @@
 
     // user-defined masks to test for good/bad pixels (build from recipe list if not yet set)
-    psMaskType maskVal = psMetadataLookupU8(&status, recipe, "MASK.PSPHOT"); // Mask value for bad pixels
+    psImageMaskType maskVal = psMetadataLookupImageMask(&status, recipe, "MASK.PSPHOT"); // Mask value for bad pixels
     assert (maskVal);
 
     // bit to mask the cosmic-ray pixels
-    psMaskType crMask  = pmConfigMaskGet("CR", config); // Mask value for cosmic rays
+    psImageMaskType crMask  = pmConfigMaskGet("CR", config); // Mask value for cosmic rays
 
     float CR_NSIGMA_LIMIT = psMetadataLookupF32 (&status, recipe, "PSPHOT.CR.NSIGMA.LIMIT");
@@ -63 +63 @@
         psF32 **resid  = source->pixels->data.F32;
         psF32 **weight = source->weight->data.F32;
-        psU8 **mask    = source->maskObj->data.U8;
+        psImageMaskType **mask    = source->maskObj->data.PS_TYPE_IMAGE_MASK_DATA;
 
         // check for extendedness: measure the delta flux significance at the 1 sigma contour
@@ -191 +191 @@
 // deviations?
 static float psphotModelContour(const psImage *image, const psImage *weight, const psImage *mask,
-                                psMaskType maskVal, const pmModel *model, float Ro)
+                                psImageMaskType maskVal, const pmModel *model, float Ro)
 {
     psF32 *PAR = model->params->data.F32; // Model parameters
@@ -260 +260 @@
 }
 
-bool psphotMaskCosmicRay_New (psImage *mask, pmSource *source, psMaskType maskVal, psMaskType crMask) {
+bool psphotMaskCosmicRay_New (psImage *mask, pmSource *source, psImageMaskType maskVal, psImageMaskType crMask) {
 
     // replace the source flux
@@ -294 +294 @@
         
         for (int ix = xs; ix < xe; ix++) {
-            mask->data.U8[iy][ix] |= crMask;
+            mask->data.PS_TYPE_IMAGE_MASK_DATA[iy][ix] |= crMask;
         }
     }
@@ -300 +300 @@
 }
 
-bool psphotMaskCosmicRay_Old (pmSource *source, psMaskType maskVal, psMaskType crMask) {
+bool psphotMaskCosmicRay_Old (pmSource *source, psImageMaskType maskVal, psImageMaskType crMask) {
 
     source->mode |= PM_SOURCE_MODE_CR_LIMIT;
@@ -320 +320 @@
         float SN = pixels->data.F32[yo][ix] / sqrt(weight->data.F32[yo][ix]);
         if (SN > SN_LIMIT) {
-            mask->data.U8[yo][ix] |= crMask;
+            mask->data.PS_TYPE_IMAGE_MASK_DATA[yo][ix] |= crMask;
         }
     }
@@ -326 +326 @@
         float SN = pixels->data.F32[yo][ix] / sqrt(weight->data.F32[yo][ix]);
         if (SN > SN_LIMIT) {
-            mask->data.U8[yo][ix] |= crMask;
+            mask->data.PS_TYPE_IMAGE_MASK_DATA[yo][ix] |= crMask;
         }
     }
@@ -339 +339 @@
 
             bool valid = false;
-            valid |= (mask->data.U8[iy+1][ix] & crMask);
-            valid |= (ix > 0) ? (mask->data.U8[iy+1][ix-1] & crMask) : 0;
-            valid |= (ix <= mask->numCols) ? (mask->data.U8[iy+1][ix+1] & crMask) : 0;
+            valid |= (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy+1][ix] & crMask);
+            valid |= (ix > 0) ? (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy+1][ix-1] & crMask) : 0;
+            valid |= (ix <= mask->numCols) ? (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy+1][ix+1] & crMask) : 0;
 
             if (!valid) continue;
-            mask->data.U8[iy][ix] |= crMask;
+            mask->data.PS_TYPE_IMAGE_MASK_DATA[iy][ix] |= crMask;
         }
     }
@@ -355 +355 @@
 
             bool valid = false;
-            valid |= (mask->data.U8[iy-1][ix] & crMask);
-            valid |= (ix > 0) ? (mask->data.U8[iy-1][ix-1] & crMask) : 0;
-            valid |= (ix <= mask->numCols) ? (mask->data.U8[iy-1][ix+1] & crMask) : 0;
+            valid |= (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy-1][ix] & crMask);
+            valid |= (ix > 0) ? (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy-1][ix-1] & crMask) : 0;
+            valid |= (ix <= mask->numCols) ? (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy-1][ix+1] & crMask) : 0;
 
             if (!valid) continue;
-            mask->data.U8[iy][ix] |= crMask;
+            mask->data.PS_TYPE_IMAGE_MASK_DATA[iy][ix] |= crMask;
         }
     }