PISA demonstration script
#+
# Name:
# pisa_demo
# Purpose:
# Script to demostrate PISA functions.
# Language:
# C shell
# Invocation:
# source pisa_demo
# Authors:
# PDRAPER: Peter Draper (STARLINK)
# {enter_new_authors_here}
# History:
# 7-MAY-1992 (PDRAPER):
# Original version.
# {enter_changes_here}
#-
#
# Find out which graphics device we're using.
echo -n "Device (image display):"
set device = $<
if ( $device == "" ) then
exit
endif
#
# initialise PISA
#
pisa
#
# display the image
#
echo " "
echo "......Displaying image......"
pisagrey noaxes device=$device in=/star/bin/pisa/frame drange=\[477,700\] \
accept reset
sleep 30
#
echo " "
echo "......Finding objects using PISAFIND......"
if ( -e pisafind.dat ) 'rm' pisafind.dat
if ( -e pisasize.dat ) 'rm' pisasize.dat
pisafind in=/star/bin/pisa/frame minpix=8 method=1 \
results=pisafind.dat sizes=pisasize.dat reset accept
#
echo " "
echo "......Plotting objects using PISAPLOT......"
pisaplot results=pisafind.dat device=$device palnum=2 reset accept
sleep 30
#
echo " "
echo "......Fitting stars using PISAFIT and PISA profiling function......"
pisafit in=/star/bin/pisa/frame positions=/star/bin/pisa/frame_stars.acc \
minmode=n radius=9 again=f weighted=f background=492.2 \
device=$device reset accept
sleep 30
#
echo " "
echo "......Generating model frame using PISAGEN and profile fit......"
pisagen input=/star/bin/pisa/frame positions=pisafind.dat \
output=model background=0 noise=g sigma=30 reset accept
#
# display the generated model data
echo " "
echo "......Displaying model frame......"
pisagrey model noaxes drange=\[0,450\] device=$device reset accept
sleep 30
#
# object classification
echo " "
echo "......Modifying PISA parameters using PISAPEAK......"
if ( -e pisapeak.dat ) 'rm' pisapeak.dat
pisapeak in=/star/bin/pisa/frame finddata=pisafind.dat \
results=pisapeak.dat reset accept
#
# apply a cut to the data in peakedness
echo " "
echo "......Applying cut to peakedness parameter using PISACUT......"
if ( -e stars1 ) 'rm' stars1
if ( -e gals1 ) 'rm' gals1
pisacut input=pisapeak.dat column=2 thresh=0.85 lower=stars1 \
higher=gals1 reset accept
#
# separate pisafind results
echo " "
echo "......Matching PISACUT classes with pisafind data using PISAMATCH......"
if ( -e stars11 ) 'rm' stars11
if ( -e gals11 ) 'rm' gals11
pisamatch one=stars1 two=pisafind.dat out=stars11 reset accept
pisamatch one=gals1 two=pisafind.dat out=gals11 reset accept
#
# display these
pisagrey /star/bin/pisa/frame drange=\[477,700\] \
device=$device noaxes reset accept
echo " "
echo "......Displaying stars (in red) using PISAPLOT......"
pisaplot overlay=t clear=f results=stars11 device=$device palnum=2 \
noannota thick=1.5 reset accept
#
echo " "
echo "......Displaying galaxies (in green) using PISAPLOT......"
pisaplot overlay=t clear=f results=gals11 device=$device palnum=3 \
noannota thick=1.5 reset accept
# end