Changeset 4954 for trunk/psphot/src/psLibUtils.c

Timestamp:

Sep 6, 2005, 5:33:01 PM (21 years ago)

Author:

eugene

Message:

psphot reorganization: merged in pmObjects code from psModule, modified to work with psLib 0.7.0

File:

• trunk/psphot/src/psLibUtils.c (modified) (15 diffs)

trunk/psphot/src/psLibUtils.c

@@ -1 @@
-# include "psphot.h"
+# include <strings.h>  // for strncasecmp
+# include <pslib.h>
+
+// XXX EAM : this is NOT included in the C99 headers ??
+FILE *fdopen(int fildes, const char *mode);
 
 // XXX EAM : these utility functions should be added back into PSLib
@@ -29 +33 @@
     if (timers == NULL) timers = psHashAlloc (16);
 
-    start = psTimeGetTime (PS_TIME_UTC);
+    start = psTimeGetNow (PS_TIME_UTC);
     psHashAdd (timers, name, start);
     psFree (start); 
@@ -47 +51 @@
     if (start == NULL) return (0);
 
-    mark = psTimeGetTime (PS_TIME_UTC);
+    mark = psTimeGetNow (PS_TIME_UTC);
     delta = psTimeDelta (mark, start);
     psFree (mark);
@@ -53 +57 @@
 
     return (delta);
+}
+
+// find the location of the specified argument
+int psArgumentGet (int argc, char **argv, char *arg) {
+
+    int i;
+
+    for (i = 0; i < argc; i++) {
+        if (!strcmp(argv[i], arg))
+            return (i);
+    }
+  
+    return ((int)NULL);
+}
+
+// remove the specified argument (by location)
+int psArgumentRemove (int N, int *argc, char **argv) {
+
+    int i;
+
+    if ((N != (int)NULL) && (N != 0)) {
+        (*argc)--;
+        for (i = N; i < *argc; i++) {
+            argv[i] = argv[i+1];
+        }
+    }
+
+    return (N);
 }
 
@@ -114 +146 @@
 }
 
-// find the location of the specified argument
-int psArgumentGet (int argc, char **argv, char *arg) {
-
-    int i;
-
-    for (i = 0; i < argc; i++) {
-        if (!strcmp(argv[i], arg))
-            return (i);
-    }
-  
-    return ((int)NULL);
-}
-
-// remove the specified argument (by location)
-int psArgumentRemove (int N, int *argc, char **argv) {
-
-    int i;
-
-    if ((N != (int)NULL) && (N != 0)) {
-        (*argc)--;
-        for (i = N; i < *argc; i++) {
-            argv[i] = argv[i+1];
-        }
-    }
-
-    return (N);
-}
-
 // alternate implementation of this function from pmObjects.c
+// XXX EAM : move this to pmObjects_EAM.c
 psVector *psGetRowVectorFromImage(psImage *image,
                                   psU32 row)
 {
-    PS_IMAGE_CHECK_NULL(image, NULL);
-    PS_IMAGE_CHECK_TYPE(image, PS_TYPE_F32, NULL);
+    PS_ASSERT_IMAGE_NON_NULL(image, NULL);
+    PS_ASSERT_IMAGE_TYPE(image, PS_TYPE_F32, NULL);
 
     psVector *tmpVector = psVectorAlloc(image->numCols, PS_TYPE_F32);
@@ -201 +206 @@
 }
 
-void psImageSmooth (psImage *image, float sigma, float Nsigma) {
+// XXX EAM : this is now in psLib
+void psImageSmooth_EAM (psImage *image, float sigma, float Nsigma) {
 
     int Nx, Ny, Npixel, Nrange;
@@ -314 +320 @@
 }           
 
+// count number of pixels with given mask value
+int psImageCountPixelMask (psImage *mask, psU8 value) 
+{
+    int Npixels = 0;
+  
+    for (int i = 0; i < mask->numRows; i++) {
+        for (int j = 0; j < mask->numCols; j++) {
+            if (mask->data.U8[i][j] & value) {
+                Npixels ++;
+            }
+        }
+    }
+    return (Npixels);
+}
+
 // define a square region centered on the given coordinate
+// XXX EAM : this is now in psLib
+# if (0)
 psRegion *psRegionSquare (psF32 x, psF32 y, psF32 radius) {
     psRegion *region;
@@ -321 +344 @@
     return (region);
 }
+# endif
 
 // set actual region based on image parameters:
@@ -327 +351 @@
 //     frame, which means the negative values should subtract from Nx,Ny of
 //     the parent, not the child.  but, we don't carry the dimensions of the 
-//     parent in the psImage container.  for now, us the child Nx,Ny
-// force range to be on this subimage 
-// XXX EAM : this needs to be changes to use psRegion rather than psRegion*
+//     parent in the psImage container.  for now, use the child Nx,Ny
+//     force range to be on this subimage 
+// XXX EAM : this needs to be changed to use psRegion rather than psRegion*
+// XXX EAM : this is now in psLib
+# if (0)
 psRegion *psRegionForImage (psRegion *out, psImage *image, psRegion *in) {
     
@@ -356 +382 @@
     return (out);
 }
+# endif
 
 // mask the area contained by the region
 // the region is defined wrt the parent image
+// XXX EAM : this is now in psLib
+# if (0)
 void psImageMaskRegion (psImage *image, psRegion *region, bool logical_and, int maskValue) {
 
@@ -375 +404 @@
     }
 }
+# endif
 
 // mask the area not contained by the region
 // the region is defined wrt the parent image
+// XXX EAM : this is now in psLib
+# if (0)
 void psImageKeepRegion (psImage *image, psRegion *region, bool logical_and, int maskValue) {
 
@@ -396 +428 @@
     }
 }
+# endif
 
 // mask the area contained by the region
 // the region is defined wrt the parent image
+// XXX EAM : this is now in psLib
+# if (0)
 void psImageMaskCircle (psImage *image, double x, double y, double radius, bool logical_and, int maskValue) {
 
@@ -419 +454 @@
     }
 }
+# endif
 
 // mask the area contained by the region
 // the region is defined wrt the parent image
+// XXX EAM : this is now in psLib
+# if (0)
 void psImageKeepCircle (psImage *image, double x, double y, double radius, bool logical_and, int maskValue) {
 
@@ -442 +480 @@
     }
 }
+# endif
 
 psVector *psVectorCreate (double lower, double upper, double delta, psElemType type) {
@@ -455 +494 @@
   return (out);
 }
+
+// XXX EAM a utility function
+bool p_psVectorPrintRow (int fd, psVector *a, char *name)
+{
+
+    FILE *f;
+    f = fdopen(fd, "a+");
+    fprintf (f, "vector: %s\n", name);
+
+    for (int i = 0; i < a[0].n; i++) {
+        fprintf (f, "%f  ", p_psVectorGetElementF64(a, i));
+    }
+    fprintf (f, "\n");
+    fclose(f);
+    return (true);
+}
+// XXX EAM is the use of fdopen here a good way to do this?