My current scripting laguage of choice is python as its use is becoming increasingly common in astrophysics. It also has access to a host excellent third party modules. Below are a number of scripts I have written which you might find useful.
Photometry scripts
The python script 'polyphot.py' allows the user to display a FITS image (single plane only) and measure the integrated flux density of an irregular object by drawing polygons. Currently it assumes radio-data calibrated in Jy/beam and divides by the beam area.
Similarly, the script 'aperphot.py' draws a circular aperture and an annulus around a position in a Jy/beam calibrated FITS file in order to measure the net integrated flux.
USAGE: aperphot.py -x <deg> -y <deg> -s <src_rad_arcsec> -i <inner_rad_arcsec> -o <outer_rad_arcsec> <file.fits>
ATAC data reduction scripts
The TAR archive file 'atprocesspy.tar.gzip' contains scripts written in Python for the reduction of ATCA CABB continuum and spectral line data. These are primarily intended to be driven manually and allow the user to rapidly reduce their interferometry datasets. The file 'README_atprocesspy.txt' has information on installing and running the scripts. Full installs of miriad and python are required.
ASAP GUI for spectral line observations
The python script 'QuickData.py' is a graphical interface to the ATNF's ASAP spectral line reduction package. It is specifically configured for position-switched observations. This script is very much a work in progress!
Also available is a script to manually fourier-flag specific frequency componenst of spectra in an ASAP/SDF dataset in order to remove baseline ripples. 'ffilter.py' allows the user to fit and subtract multiple Gaussian lines before flagging the data in the fourier domain using one or more cosine-edge filters. The spectrum is then reconstructed.