#!/bin/sh
# imsex: run SExtractor on a FITS image
# By Doug Mink, Smithsonian Astrophysical Observatory
# September 15, 2000

if test $# = 0
    then
    echo "imsex: Run SExtractor on a FITS image and return stars in DAOFIND x y mag"
    echo "  usage: imsex file.fits returns image coordinates in file.sex"
    echo "  -v  note progress"
    echo "  -m  extract stars from all images in multi-extension file"
    echo "  -n  number of brightest stars to return (default is 100)"
    exit
fi
if test $1 = 'help'
    then
    echo "Run SExtractor on a FITS image and return brightest stars in DAOFIND x y mag"
    echo "-v  note progress"
    echo "-n  number of brightest stars to return (default is 100)"
    exit
fi
if test $1 = "-v"
    then
    v=1
    shift
else
    v=
fi
if test $1 = "-m"
    then
    m=1
    shift
else
    m=
fi
if test $1 = "-n"
    then
    s=$2
    shift
    shift
else
    s=100
fi

f=$1
g=`echo $f | grep ,`
if test $g
    then
    r=`fileroot $1`
    echo $f '->' $r
elif test $m
    then
    n=1
    r=$f
    f=$r,$n
else
    r=$f
fi
while test -r $r
    do
    if test $m
	then
	if test $v
	    then
	    e=`imrot -v $f`
	else
	    e=`imrot $f`
	fi
    elif test $g
	then
	if test $v
	    then
	    e=`imrot -v $f`
	else
	    e=`imrot $f`
	fi
    else
	e=$f
    fi
    if test $e
	then
	echo $f '->' $e
	if test $v
	    then
	    echo Running SExtractor on $e
	    fi
	sex -c default.sex $e

	f=`filename $e`
	c=`fileroot $f`.sex
	if test $v
	    then
	    echo Saving $s brightest stars in $c
	    fi
	sort -n +2.0 test.cat | head -$s > $c
	/bin/rm test.cat
	if test $m
	    then
	    n=`expr $n + 1`
	    f=$r,$n
	elif test $# = 0
	    then
	    exit
	else
	    shift
	    g=`echo $1 | grep ,`
	    if test $g
		then
		f=`fileroot $1`
	    else
		f=$1
	    fi
	    echo $f
	fi
    else
	exit
    fi
done

# May  2 2000	Use fileroot and filename instead of root and fname
# Sep 12 2000	Fix help if no arguments
# Sep 15 2000	Add extraction of single image from multi-extension file