Changeset 4216 for trunk/psphot/src/psphot-utils.c

Timestamp:

Jun 12, 2005, 3:49:21 PM (21 years ago)

Author:

eugene

Message:

rework with masks and noise array

File:

• trunk/psphot/src/psphot-utils.c (modified) (4 diffs)

trunk/psphot/src/psphot-utils.c

@@ -129 +129 @@
 }
 
+bool DumpImage (psImage *image, char *filename) {
+
+    unlink (filename);
+    psFits *fits = psFitsAlloc (filename);
+    psFitsWriteImage (fits, NULL, image, 0, NULL);
+    psFree (fits);
+    return;
+}
+
 void pmSourceMaskSaturated (psSource *source, float saturate) {
 
@@ -142 +151 @@
 }
 
-void pmSourceMaskRegion (psSource *source, psRegion *region) {
-
-    psImage *image = source->pixels;
-    psImage *mask  = source->mask;
+// mask the area contained by the region
+// the region is defined wrt the parent image
+void psImageMaskRegion (psImage *image, psRegion *region, bool logical_and, int maskValue) {
+
+    for (int iy = 0; iy < image->numRows; iy++) {
+        for (int ix = 0; ix < image->numCols; ix++) {
+            if (ix + image->col0 <  region->x0) continue;
+            if (ix + image->col0 >= region->x1) continue;
+            if (iy + image->row0 <  region->y0) continue;
+            if (iy + image->row0 >= region->y1) continue;
+            if (logical_and) {
+                image->data.U8[iy][ix] &= maskValue;
+            } else {
+                image->data.U8[iy][ix] |= maskValue;
+            }
+        }
+    }
+}
+
+// mask the area not contained by the region
+// the region is defined wrt the parent image
+void psImageKeepRegion (psImage *image, psRegion *region, bool logical_and, int maskValue) {
 
     for (int iy = 0; iy < image->numRows; iy++) {
@@ -155 +182 @@
             continue;
         maskit:
-            mask->data.U8[iy][ix] |= 0x02;
+            if (logical_and) {
+                image->data.U8[iy][ix] &= maskValue;
+            } else {
+                image->data.U8[iy][ix] |= maskValue;
+            }
         }
     }
@@ -161 +192 @@
 
 // mask the area contained by the region
-void psImageMaskRegion (psImage *image, psRegion *region, int maskValue) {
-
-    for (int iy = 0; iy < image->numRows; iy++) {
-        for (int ix = 0; ix < image->numCols; ix++) {
-            if (ix + image->col0 <  region->x0) continue;
-            if (ix + image->col0 >= region->x1) continue;
-            if (iy + image->row0 <  region->y0) continue;
-            if (iy + image->row0 >= region->y1) continue;
-            image->data.U8[iy][ix] |= maskValue;
-        }
-    }
-}
-
-// mask the area not contained by the region
-void psImageKeepRegion (psImage *image, psRegion *region, int maskValue) {
-
-    for (int iy = 0; iy < image->numRows; iy++) {
-        for (int ix = 0; ix < image->numCols; ix++) {
-            if (ix + image->col0 <  region->x0) goto maskit;
-            if (ix + image->col0 >= region->x1) goto maskit;
-            if (iy + image->row0 <  region->y0) goto maskit;
-            if (iy + image->row0 >= region->y1) goto maskit;
-            continue;
-        maskit:
-            image->data.U8[iy][ix] |= maskValue;
+// the region is defined wrt the parent image
+void psImageMaskCircle (psImage *image, double x, double y, double radius, bool logical_and, int maskValue) {
+
+    double R2 = PS_SQR(radius);
+
+    for (int iy = 0; iy < image->numRows; iy++) {
+        for (int ix = 0; ix < image->numCols; ix++) {
+            dx = ix + image->col0 - x;
+            dy = iy + image->row0 - y;
+            r2 = PS_SQR(dx) + PS_SQR(dy);
+            if (r2 > R2) continue;
+            if (logical_and) {
+                image->data.U8[iy][ix] &= maskValue;
+            } else {
+                image->data.U8[iy][ix] |= maskValue;
+            }
+        }
+    }
+}
+
+// mask the area contained by the region
+// the region is defined wrt the parent image
+void psImageKeepCircle (psImage *image, double x, double y, double radius, bool logical_and, int maskValue) {
+
+    double R2 = PS_SQR(radius);
+
+    for (int iy = 0; iy < image->numRows; iy++) {
+        for (int ix = 0; ix < image->numCols; ix++) {
+            dx = ix + image->col0 - x;
+            dy = iy + image->row0 - y;
+            r2 = PS_SQR(dx) + PS_SQR(dy);
+            if (r2 < R2) continue;
+            if (logical_and) {
+                image->data.U8[iy][ix] &= maskValue;
+            } else {
+                image->data.U8[iy][ix] |= maskValue;
+            }
         }
     }