Caltech Astronomy

Caltech Optical Observatories

Optical/Infrared Instrumentation


Palomar LGS-AO Software Download Page

The IDL routines provided on this page is intended to assist observers in preparing LGS-AO programs, evaluating AO performance, and reducing AO science data. You are invited to download, use, and modify this software to suit your needs, but it is of course provided with no guarantee of support by the authors.

Many of the programs use IDL Astronomy User Library subroutines. Several also use the function optimization routines of Craig Markwardt. All include a header which details the syntax and use of the program.

LGS Observation Planning

July 08: Due to recent changes in the Vizier format, these routines require the most recent versions of QUERYVIZIER and WEBGET.

FINDTTREF.pro Find USNO-B1 or 2MASS tip/tilt stars near a list of LGS targets.
FINDLGSPSF.pro Find IR-bright PSF stars with similar target/TTref geometry.

AO performance analysis

PH_STREHL.pro Computes the Strehl of a single PHARO image. Uses PH_PSF.pro.
PH_PSF.pro Compute a diffraction-limited monochromatic PSF from a PHARO image header.

PHARO data reduction

READPHARO.pro Read a PHARO NxNx4 or NxNx8 FITS image cube and return the 2-D image.
PH_MKDARK.pro Process and save PHARO dark frames.
PH_MKFLAT.pro Process and save PHARO flat and bad pixel frames.
PH_REGISTER.pro Reduce sets of PHARO images taken either on/off source or using a dither pattern. Saves individual reduced images and an averaged stacked image. Best for faint sources (uses 2-stage registration by cross-correlation on a user-chosen object in the field).
PH_DITHER.pro Reduce sets of PHARO images taken either in a dither or mosaic pattern. Saves individual reduced images, mosaic, and mean image. Best for bright targets (registrers frames by cross-correlation over the overlap regions). Updated 4/10/09

Here is an example using the above routines to reduce one dither sequence and calibration files:
ddir = 'data/pharo_070425/'
pdir = 'proc/pharo_070425/'
PH_MKDARK, data=ddir, proc=pdir, frno=99+INDGEN(10) ; 19.8s darks
PH_MKFLAT, data=ddir, proc=pdir, frno=49+INDGEN(5) ; Ks flats (19.8s)
PH_DITHER, data=ddir, proc=pdir, obj=696+INDGEN(20) ; M20 Ks mosaic

Subroutines

READSTARLIST.pro Read Keck or Palomar format starlists into a structure array.
STARSTRING.pro Convert target structure into a Keck-format string.
RFITS.pro Read multiple FITS images into an image cube (use /PHARO keyword).
MEDIM.pro Compute the pixel-by-pixel median of a series of images (use /PHARO keyword). Updated 4/14/09
BPIXFIX.pro Replace bad pixels with the median of surrounding good pixels.
CRAYFIX.pro Find and replace anomalously high pixels.
PFFT.pro Returns a 2-D FFT, padded with zeros.
FSHIFT.pro Fourier shift function (sub-pixel shifts with no resolution loss.)


For more information please contact abouchez@astro.caltech.edu.