Initial commit

1 files changed, 139 insertions, 0 deletions

diff --git a/vo/votools/imgcat.cl b/vo/votools/imgcat.cl
new file mode 100644
index 00000000..7abdf74a
--- /dev/null
+++ b/vo/votools/imgcat.cl
@@ -0,0 +1,139 @@
+#{  IMGCAT -- Create a catalog of stellar objects in an image.
+
+procedure imgcat (image, catalog)
+
+string	image			{ prompt = "Image name" 	    	      }
+string  catalog 		{ prompt = "Output catalog name"	      }
+
+string	format    = "text"	{ prompt = "Output format"  		      }
+bool	verbose   = no		{ prompt = "Verbose?"	    		      }
+
+bool	dosplit   = yes		{ prompt = "Split detected objects?"	      }
+string	method    = "thresh"	{ prompt = "Detection method (ace|thresh)",
+				    enum = "ace|thresh"  		      }
+real	ellipse   = 0.33	{ prompt = "Ellipticity cutoff?"  	      }
+real	threshold = 5.		{ prompt = "Detected threshold?"	      }
+
+int	ndetected = 0		{ prompt = "Num detected objects"	      }
+int	nel_clip  = 0		{ prompt = "Num clipped for ellipticity"      }
+int	nobjs     = 0		{ prompt = "Num remaining objects"	      }
+  
+begin
+    string  img, fmt, catdef, cat, ccscript
+    string  tcat, tc1, tc2, tc3, tc4, tdef, expr, detcode
+    int     nstars, nx, ny
+    bool    verb, split
+    real    ellip, sigma, stddev, immean, lthresh
+
+
+    # Check for proper packages and reset environment.
+        ;
+    reset imtype  = "fits"
+    reset clobber = yes
+    #flpr 0
+
+
+    # Get the task parameters.
+    img	     = image
+    cat	     = catalog
+    verb     = verbose
+    fmt      = format
+    ellip    = ellipse
+    detcode  = method
+    split    = dosplit
+    sigma    = threshold
+
+
+    # If we have zero values, replace them with the image mean.
+    minmax (img, update-, >& "dev$null")
+    if (minmax.minval == 0) {
+        imstat (img, field="mean", format-) | scan (immean)
+        imreplace (img, immean, upper=0.0)
+    }
+
+
+    # Do the object detection.
+    if (detcode == "ace") {
+        # Initialize
+        catdef = mktemp ("/tmp/cd")
+
+        print ("WX\nWY\nPX\nPY\nFLUX\nELLIP\n", > catdef)
+	
+        # Remove the DATASEC keyword if present since it may not be valid for
+        # the image, but will be used by ACE anyway.
+	hselect (img, "DATASEC", yes) | scan (s1)
+	if (nscan() == 1)
+	   hedit (img, "DATASEC", delete+, update+, verify-, >& "dev$null")
+
+
+	# Run the DETECT task to locate the stars.
+	iferr {
+            detect (img, catdefs=catdef, catalog=cat, hsigma=sigma,
+	        doeval+, dosplit=split, dogrow+)
+return
+            detect (img, catdefs=catdef, catalog=cat, hsigma=sigma,
+	        doeval+, dosplit=split, dogrow+) | \
+	    match ("detected", "STDIN") | scan (nstars)
+            delete (catdef, verify-)
+	} then {
+	    error (0, "Error creating image catalog\n")
+	}
+
+        # Now filter the catalog to pull out the "stellar" objects identified
+        # by a low ellipticity.  Sort on FLUX.
+        tcat = mktemp ("/tmp/tc")
+        tinfo (cat, ttout-)
+        nstars = tinfo.nrows
+        fields (cat, "1-6", > tcat)
+        if (ellip > 0.0) {
+            expr = "c6 == INDEF || c6 < " // ellip
+            tselect (tcat, cat, expr)
+            tcopy (cat, tcat, >& "dev$null")
+        }
+        tinfo (tcat, ttout-)
+        tsort (tcat, col="c5", ascend-)
+
+
+        tcalc (tcat, "c1", "c1 * 15.0")
+        tdump (tcat, cd="dev$null", pf="dev$null", data="STDOUT",
+	    col="c1,c2,c3,c4,c5", row="1-", > cat)
+
+        # Sort the table and save the results to parameters.
+        tsort (cat, col="c4")
+
+        # Clean up.
+        if (access (tcat) == yes) 
+	    delete (tcat, verify-, >& "dev$null")
+
+    } else if (detcode == "thresh") {
+
+	# Get the stddev of the image
+	imstat (img, fields="stddev", format-) | scan (stddev)
+	lthresh = 1 * stddev
+
+	starfind (img, "image.obj", 1., lthresh, wcs="world", 
+	    wxf="%12.2H", wyf="%12.1h")
+
+	i = 0
+	printf ("|id            |ra            |dec           |\n", > cat)
+	list = "image.obj"
+	while (fscan (list, z, z, x, y) != EOF) {
+	    if (i > 18)
+	        printf (" %14d %14.8f %14.8f\n", (i-18), x, y, 	    >> cat)
+	    i = i + 1
+	}
+	list = ""
+	nstars = i - 18
+
+	delete ("image.obj", verify-, 	>& "dev$null")
+	ellip = 0.0				# disable ellipticity check
+    }
+
+    ndetected = nstars
+    nobjs     = tinfo.nrows
+
+return
+    # Run the cross-compare script.
+    ccscript = osfn ("votools$cross_comp.sh")
+    printf ("!%s %s 1.5 TWOMASS_PSC\n", ccscript, cat) | cl()
+end