CTIO
Published on CTIO (http://www.ctio.noao.edu/noao)

CTIO Home > Astronomers > Instruments > IR Instruments > ISPI

ISPI

The CTIO Blanco 4.0-m Infrared Side Port Imager

NOTE: Not currently in use

The Infrared Side Port Imager ISPI (eye-spy) is a facility infrared camera at the CTIO Blanco 4-m telescope, serving a broad range of science programs through the following capabilities:

1-2.4 micron imaging with 2K x 2K HgCdTe HAWAII-2 array

0.3 arcsec/pixel sampling matched to f/8 IR image quality of ~0.6 arcsec

10.25 x 10.25 arcmin field of view,

broad band J,H and Ks, as well as a set of narrow band filters.

 

ISPI Overview

The Infrared Side Port Imager ISPI (eye-spy) is a facility infrared camera at the CTIO Blanco 4-m telescope, serving a broad range of science programs through the following capabilities:

  • 1-2.4 micron imaging with 2K x 2K HgCdTe HAWAII-2 array
  • 0.3 arcsec/pixel sampling matched to f/8 IR image quality of ~0.6 arcsec
  • 10.25 x 10.25 arcmin field of view,
  • broad band J,H and Ks, as well as a set of narrow band filters.

ISPI is one component of a fixed instrument complement for future Blanco operations; the others are Prime Focus MOSAIC for optical imaging, and Hydra for multiobject optical spectroscopy. This fixed complement will lower operations costs on the Blanco as CTIO meets its committments for SOAR commissioning and operation, as well as Gemini operations support.

While ISPI's basic capabilities derive from its science mission, the project also faced realities of cost and schedule, space envelope at the 4-m Cass focus, and the future mix of instruments available to the community on the Blanco, SOAR, and Gemini South telescopes. The result is a conservative design which employs existing elements wherever possible (e.g. the optical design), uses a single large array, and has no spatial filtering (occulting masks) or spectroscopic capability. The ISPI project delivered a highly capable instrument quickly and at low cost. It complements future IR capability for high spatial resolution imaging and spectroscopy on SOAR and Gemini South, and a very-wide-field IR imager to be shared between the NOAO 4-m telescopes (NEWFIRM).

Block Diagram [1]

ISPI Filters [2]

ISPI Mechanical Design

ISPI is side-looking so that it may be mounted on the 4-m simultaneously with Hydra II. The space envelope is constrained by the mirror cell above and Hydra below. The refractive optical system is enclosed in a cylindrical, straight-through design with a strong heritage from the facility IR camera CIRIM. Liquid cryogens are used for instrument cooling.

The design separates into two dewars with separate cryogen tanks which surround the optical path. The fore dewar has the entrance window, collimator, and filter wheel assemblies. The aft dewar holds the camera and detector assemblies. The cylindrical construction provides a very stiff assembly to meet flexure requirements. The diameter is 14 inches and the total length of the dewar assembly is 45 inches, the weight for the complete assembly is expected to be ~220lbs.

The filter wheels provide 16 positions including darks and opens. There is a single fixed Lyot stop. Alignment to the telescope is accomplished by adjustments at the telescope interface.

 

March 22, 2004

ISPI Optical Design

The ISPI optical design originated with a design produced by Charles Harmer (NOAO) for the U. Florida IR camera FLAMINGOS [3]. This design was a close match to our science-defined specifications, was already optimized for the Mayall 4-m, and allowed cost savings by joint procurement of optics. Hence we adopted it for ISPI, with a slight reoptimization of element spacings for the Blanco telescope by Harmer and Maxime Boccas (CTIO).

This is a refractive, collimator-camera system with an intermediate cold pupil image. The f/3 camera employs one aspheric element; all other surfaces are spheres. Design performance is shown graphically below with spot diagrams (box size is 2x2pixels, ie. 36x36 microns)
The camera optics coupled to the Blanco 4-m give the following performance for 80% encircled energy diameter:

 

Band Best in Field Worst in Field Spot Diagram
J 0.36" 0.48" J Spot Diagram [4]
H 0.36" 0.48" H Spot Diagram [5]
K 0.48" 0.55" K Spot Diagram [6]

 

Note: The best and worst are not obtained respectively on axis and in the corner

The optics have been produced by Janos Technology, Inc [7]. The FLAMINGOS implementation has had first light [8] confimation of excellent image quality on Kitt Peak telescopes.

 

Ray Trace

ISPI Ray Trace

March 22, 2004

ISPI GUI

ISPI is operated through a Graphical User Interface, allowing the user to control the array and the filter mechanism. Through the GUI a limited set of telescope control functions is accessible, enabling the user to perform small off sets and to set the telescope focus. Telemetry of ISPI is also made available through the GUI. The ISPI GUI is an ARCview application, which allows the user to control both the instrument (array control & filter mechanisms) as well as some limited telescope , such as small off sets and telescope focus.

[9]

March 22, 2003

ISPI Performance

ISPI throughput

based on data taken on Dec. 21, 2004

 


Band Background Flux
(electrons/sec)
per pixel
Integrated Flux
(electrons/sec)
for m=15 star
J 165 3000
H 1125 4565
K_short 1925 3085

The exposure calculator [10] on the CTIO Infrared Instrument web page [11] has been updated accordingly.

 

March 26, 2005

ISPI Quick reference guide

General information

The Infrared Side Port Imager ISPI (eye-spy) is a facility imager at the CTIO Blanco 4-m telescope, serving a broad range of science programs through the following capabilities:

  • 1-2.4 micron imaging with 2K x 2K HgCdTe HAWAII-2 array
  • 0.3 arcsec/pixel sampling matched to f/8 IR image quality of ~0.6 arcsec
  • 10.25 x 10.25 arcmin field of view,
  • broad band J,H and Ks, as well as a set of narrow band filters.
  • A list of available filters can be found on the ISPI filter page.


Performance

During the January 2003 engineering night weather permitted us to measure the throughput of ISPI at the telescope, see table. The ISPI throughput is very similar to the throughput of Flamingos at the Kitt Peak 4m. Exposure times and signal-to-noise ratios can be estimated using the exposure time calculator [10] on the CTIO Infrared Instrument web page [12] .

ISPI performancebased on data taken on Jan. 19, 2003
Band Background Flux
(electrons/sec)
per pixel
Integrated Flux
(electrons/sec)
for m=15 star
point source detection limit
5 sigma in
60 sec
J 310 4000 19.6
H 1560 5000 18.9
K' 2160 3300 18.3

The HAWAII-2 array in ISPI

The array, as set-up in ISPI since May 2003, has the following performance characteristics:

  • Gain: 4.25 e-/ADU
  • Well depth: 150,000 e-
  • Linearity: 1% full frame linear range is ~13,000 ADU's
  • Minimum exposure time: 3.227s

NOTE: For data taken prior to May 2003 please consult the ISPI support staff. [13]

Observing with ISPI

ISPI is run through a Graphical User Interface [9], in which all array control parameters, such as exposure time, number of co-adds, as well as filter, telescope focus and telescope off-set can be set.

The GUI can also initiate a series of observations, to conduct for example a dither sequence, a focus sequence, or a sequence of darks. These sequences are simple ASCII text files, which can be prepared in advance. Examples of such observing sequences are in the following links in H [14], J [15] and K [16].

February 13, 2007

ISPI Users Guide

    See links below, about startup, GUI, observing notes, flat fields, linearity correction, data reduction and guiding.

 

September 23, 2005

ISPI Start-Up Procedure

Normal Startup

Check all desktops on the observers console to be sure there is NOT a "TightVNC: ispi's X desktop (ctioa9:9)" already running.

  • If there IS NOT already a "TightVNC: ispi's X desktop (ctioa9:9)" running then

    • Single click on the 'ISPI VNC' button on the left of the screen. (IF button doesn't work, then in a terminal window, at the prompt type: display_ispi). This will pop up a vnc view desktop.  The desktop will have an xload window with 4 buttons. One of these is labeled 'ArcVIEW.'
    • Single click on the 'ArcVIEW' button. This will pop up the ISPI Control GUI (with the really cool cross section view of ISPI), and an 'ArcVIEW Main Application' window.  The Application window should have two round buttons, labelled 'Connected?' that appear red.
    • On the  ISPI Control GUI, there is a white arrow icon in the upper left corner of the window. Single click on the white arrow. This arrow should now turn black. The two buttons in the Application window should turn green.
    • Minimize the Application window by clicking on the underscore '_' icon on the blue bar, in the upper right corner of each window.  You are now ready to operate ISPI.
  • If there IS already a "TightVNC: ispi's X desktop (ctio:9)" running, it can be in a couple of states:

    • There is only the gray desktop with an xload window that has  4 buttons labelled ArcVIEW, IRAF, Vacuum and terminal If this is the case:

      • Single click on the 'ArcVIEW' button. This will pop up the ISPI Control GUI (with the really cool cross section view of ISPI), and an 'ArcVIEW Main Application' window.  The Application window should have two round buttons, labelled 'Connected?' that appear red.
      • On the  ISPI Control GUI, there is a white arrow icon in the upper left corner of the window. Single click on the white arrow. This arrow should now turn black. The two buttons in the Application window should turn green.
      • Minimize the Application window by clicking on the underscore '_' icon on the blue bar, in the upper right corner of each window. You are now ready to operate ISPI.
  • The ISPI Control GUI (with the really cool cross section view of ISPI) is already displayed on the desktop. In this case, Check the color of the arrow icon in the upper left corner of the ISPI Control GUI window.

    • If the arrow is white then Single click on the white arrow. This arrow should now turn black. The two buttons in the Application window should turn green.
    • If the arrow is black then you are already ready to operate ISPI, however CHECK the 'ArcVIEW Main Application' window to see that the 'Connected?' buttons are green. If you cannot see the 'ArcVIEW Main Application' window, you can find or open it by clicking on the bar at the bottom of the VNC desktop window.

July 07, 2005; Bob Blum

ISPI Observing Notes

Quick Reference [17]

Observing Sequences

  • ISPI data taking may be "scripted" and complete sequences executed through the ISPI GUI.
  • Several examples are available (H [14], J [15], K [16]). Each line of the text file corresponds to a single image written to disk. The key elements are filter changes (one filter per sequence in these examples), exposure time and coadds where exposure*coadds is the total time per image, and the resulting image is the sum of the number of coadds specified. For example, a 10 sec exposure of 3 coadds would produce one image with a total of 30 seconds. The image would be the sum of three 10 second frames. Currently Fowler samples, a type of read mode to reduce noise, are always set to 1. The offsets specify the dither pattern to use and can be either absolute or relative, the choice being made by selecting a button on the GUI itself.

    See /home2/ispi/observers/ on the ISPI computer for many more example sequences. Any of these can be loaded directly into the instrument control GUI by clicking on the load button in the Sequence Control panel of the GUI. Use the dialog to navigate to the ascii text file of choice.

  • Modify existing sequences by first making a copy. Each observer may create a sub-directory in /home2/ispi/observers/ to save their sequences. These will not be deleted.
  • Each step of the sequence needs to have the same number of fields (see the GUI) even if nothing is changed during that step (e.g. the filter). It is best to copy existing sequences and modify them as needed rather than starting from scratch.
  • Use care with the focus sequences as the f/8 focus mechanism has been known to go to the wrong position.

Focus

  • Obtaining a good focus takes practice. The main cause of focus changes is due to systematic temperature changes of the telescope truss. The change in focus due to temperature is -1000 steps/C for the f/8 focus used by ISPI (go to higher focus numbers as the temperature gets colder).
  • One can monitor the temperature by typing "temp4m" on a terminal on the observers console (not on the vnc session connected to ctioa9). Use one of the "upper truss" readouts.
  • The focus can be changed by using the TCS focus control on the ISPI System Control GUI (lower left of main GUI). Click on the button to set the value entered in the box. Do not send a "0" value to the focus. Look on the TCS status monitor (an independent CRT usually to your left). Find the current focus (something like 185,500 at 10C), and enter that value to start.
  • To find the best focus when you have no idea what it is, first use the last value from the previous night and estimate the change due to the difference in last night's temperature and the current temperature. Ask the night assisitant if you don't know these values. Next, begin a run of up to 10 images, each with a focus change of 1000 steps. The estimate for the best focus should be image five or six. The first image should be the farthest "in" and you should take successive images going "out" to higher numbers, each time entering a new focus in the System Control portion of the GUI. Analyze all the images, and set the final focus by first going "below" the target value and then moving the final small amount (say 200 steps) to the target by going out against gravity.
  • Once a rough focus is set by the above procedure, a better focus can be found by repeating the procedure with a step of 200 units and 5-6 images.
  • At the beginning of the night, the temperature can change significantly during the time it takes to set the focus. Beware.
  • Always set the desired focus by moving out against gravity to avoid backlash.
  • On the IRAF display, an out of focus image is typically elongated. For images which extend from upper left to lower right, focus should be adjusted to higher numbers. Go to lower numbers for images extended upper right to lower left. Use caution here since telescope motion or other problems could be confused with focus. At some point near focus, the image shape is not a good indicator of which way to adjust.

Detector Characteristics, Linearity, and Saturation

  • Detailed linearity behavior is shown here [18].
  • Keep ISPI counts below about 10000 ADU for data which are ~/< 1% non-linear, or 20000 ADU for data ~/< 2%.
  • Hard saturation will occur at about 40000 ADU.
  • The ISPI gain is 4.25e-/ADU.
  • The ISPI dark is ~0.1e-/sec for long exposures.
  • When changing filter or exposure time, the first one or few images will exhibit differing bias structure. This is a feature of the array. It means the first images will have poor background subtraction relative to others in the sequence. If this is important, plan to take one or more "junk" frames.

Flat Fields

  • Check here [19] for information on how to take flats.
  • The dome flat field lights are activated from control boxes in the computer room. Have the support astronomer or telops staff show you where the controls are and how to operate them. You will typically use only a single set of lamps controlled by the "bottom" box (of three).
  • Typical exposures will be 3-10 sec with the lamp control set to 3 amps.

Disk Space

  • Please write images to /home2/ispi/images/. You may make a subdirectory here of your choosing, e.g. 20040728/.
  • It may be necessary to backup some of your data and then delete it to make room for the rest of your run.
  • Your data are archived every night incase of a backup or disk failure. This archive is not user accessible. Contact telops or the support scientist if you have had a data disaster.
  • You should have the entire set of data for your run backed up when you leave the mountain. After that, it may be deleted any time.

Observing Efficiency

  • For K-band images, the background will set the maximum observing time. Depending on the ambient temperature this could range from about 15 to 30 seconds.
  • Use coadds to increase efficiency. Coadding images uses less time than writing the same number of separate images and dithering the telescope. I.e. for a given total exposure time, integrate as long as possible on each frame, take the minimum number of dither positions, and coadd as necessary.
  • Experience shows ISPI is about 60% efficient for 60 second K-band images (20 seconds by 3 coadds) and 10 dithers typically less than about 60''.

February 3, 2005

ISPI Flat Fields

The following summarizes some issues regarding flat fields taken with ISPI using the Ks (aka KA) and H (aka (HA) filters (J band results are expected to be similar to those for H).

  • These data were taken on 7 May 2004.

    • The dome flats use the lower lamp box at 3.00 amps
    • Exposure times are 3.3s, 5s, and 10s, for Ks, H, and J, respectively (15 frames lights on, 15 frames lights off)
    • Sky flats below were made from 100 data frames at 3.3 seconds (with rejection of 30 high and 20 low frames from the stack to account for stars in the field).
  • Ideally, flats should be constructed from series of "illumination on" minus "illumination off" images.
  • In all cases, flats should be made from series of frames with the exact same data taking parameters (exposure time, coadd, repeats, Fowler reads, etc).
  • Dome flats are recommended for Ks band data.
  • Sky flats made from data (i.e. observations) will probably work equally well as dome flats for J and H.
  •     Twighlight flats will probably give similar results for H and J as for dome flats. At Ks, twighlight flats with varying illumination in which the "low" frames are subtracted from the "high" frames might give similar (or even better) results as for Ks dome flats.
  •     Twighlight flats at Ks made by subtracting dark frames will probably give similar results as for the sky flats described below and these will not be as good as dome flats.

The following plots show the results of placing a star at 100 evenly spaced positions over the ISPI 2048x2048 detector under photometric conditions.

  • The individual frames were flat fielded with either a dome flat (as made in bullet 2 above) or a sky flat made by combining the data and subtracting a dark frame from the result.
  • All flats were fixed for bad pixels and normalized by the mean counts in the image.
  • Photometry is for a 6 pixel aperture with a 3 pixel sky annulus with inner radius of 6 pixels. No particular effort was made to optimize the aperture for varying image quality.
  • Normal dome flats at Ks or dome/sky flats at H (and probably J) should have large scale uncertainties of about ~/< +/- 3% over the full field of view.
  • The Ks flats might be affected by extra emission arising from the warm structure around M3 and a possible undersize of M3. The difference between "on" and "off" images appears to account for this putative extra emission. Sky flats would not normally account for such emission, and this is a possible explanation for the apparent better accuracy seen below for dome vs sky flats.
  • It is not known what produces the residual large scale flat fielding errors seen in the plots.
  • The following table contains text files of the photometry results which can be used to estimate second order corrections to H and Ks flats.
Ks phot data for dome flat [20] phot data for sky flat [21] average over x pixel, dome [22] average over x pixel, sky [23] average over y, dome [24] average over y, sky [25]
H phot data for dome flat [26] phot data for sky flat [27] average over x pixel, dome [28] average over x pixel, sky [29] average over y, dome [30] average over y, sky [31]

 


Figure 1a. Instrumental H Magnitude vs x pixel. Both sky and dome flats give similar results.


Figure 1b. Instrumental H Magnitude averaged over all y pixels for given x pixel. The linear fit suggests a full range of variation of +/- 2% due to large scale flat fielding uncertainty.


Figure 2a. Instrumental Ks Magnitude vs y pixel. The sky flat shows a much larger range of values than the dome flat. These represent systematic changes in the flat field, not variations due to photon statistics. Sky flats made by using illuminated frames and a dark frame should not be used unless the region of interest for objects and standards is restricted. The differences here are probably due to vingetting of the field at M3 combined with "off M3" background at Ks.


Figure 2b. Instrumental Ks magnitude averaged over all y pixels for a given x pixel. The fits show that the sky flat probably contains a component of emission comming in at the edges of the field around M3. The beam from M3 is "suppressed at the edges" by the flat.


Figure 3a. Same as Figure 1a, but for the y pixel vs H.


Figure 3b. Same as Figure 1b, but for y pixel vs H. The large scale variations are well fit by a 3rd order polynomial.


Figure 4a. Same as Figure 2a, but for the y pixel vs Ks.


Figure 4b. Same as Figure 2b, but for y pixel vs Ks.
 

ISPI Linearity Correction

A nonlinearity correction has been derived for ISPI data. This is an average correction based on all pixels. A series of J-band frames interleaved with minimum exposure dark ("bias") frames were used to derive the correction. The dark following each illuminated frame was subtracted from that illuminated frame.

The correction was derived by considering a linear fit to the counts vs exposure time for counts < about 7000 ADU. The residual of this linear fit, normalized by the fit value, plus one was then fit to a third order polynomial. The IRAF coefficients (see IRAF irlincor) are available through the cirred [32] data reduction package in the task "osiris.cl".

Last Update 06 July, 2005; Bob Blum

ISPI Data Reduction

Basics

ISPI data should be flat fielded and sky (and/or "dark") subtracted. These steps can be accomplished with may standard IRAF or IDL routines. One way is to use the basic reduction tool "osiris" provided in the cirred [33] set of IRAF scripts. A basic set of initial reduction steps would include:

  • Build flats from raw ISPI images, typically subtracting a combined "lights off" image from a combined "lights on" image. Use the cirred task med.cl (for "median combining").
  • Build a bad pixel mask from images with typical bad pixels (for example, use a lights on image and a lights off or dark image). Use the cirred task maskbad.cl.
  • Run the cirred meta script "osiris" to read in a series of images efficiently, divide each by a flat, and fix bad pixels (interpolates over bad pixels in x and y). For ISPI, choose flipy=yes. This puts N up and E to the left on an IRAF image display and is required to run the astrometric tools described below. If a bad pixel mask called "mask.fits" was created, and a flat called "flat.fits", then a typical osiris call might be:

    cl> osiris n1 n2 pre=r suf=n mask=mask.fits div=yes, dome=flat.fits

    n1 and n2 are command line params, the others can be set in the param file. Here pre is a "prefix" or image base name, suf will be the output image suffix and the images to be reduced have numbers n1 thru n2, inclusive. Suppose the prefix for input images is "r" for raw. If n1=1, and n2=3, then images r001.fits, r002.fits, r003.fits will be processed and images r001n.fits, r002n.fits, and r003n.fits will be output. Look over the cirred pages, the osiris param file, and the osiris.cl script itself for more information.

  • The flat fielded images should be sky subtracted. Separate dark subtraction is not needed as long as the images are not scaled during sky subtraction. If images are to be scaled significantly, and the dark and bias pattern is important compared to the sky level, then darks of the same exposure and coadds should be subtracted before sky subtraction (usually at the very beginning; this can be done in osiris.cl). Use the cirred script sky_sub. A typical sky_sub might be:

    cl> sky_sub 10 25 sky=kskyimage make+ sub+

    Where the image numbers work as in the osiris example. Prefixes and suffixes are set in the param file or on the command line. The two flags "make" and "sub" control how sky_sub behaves. Make will combine the indicated images into a sky frame, and sub will subtract that sky frame from the images. A set of independent skys could be combined with make+ ans sub- (to suppress subtraction) and applied to a separate set of object frames with make- (to suppress making of the sky frame) and sub+.

Astrometry

One of the most important tasks for ISPI data reduction is the assignment of a proper coordinate system and correction of image distortions. These tasks are conceptually straightforward, but often difficult to accomplish in practice (especially the latter). A preliminary set of tasks and scripts has been set up to allow for the determination of the world coordinate system keywords to be included in the image headers. A solution to higher order distortions is also available, which will allow for the registration and combination of dithered images.

Software.

The processing of ISPI images described here requires that one download and install several external packages.

  • WCSTOOLS [34]. These are C libraries and programs by Douglas Mink at the SAO. These routines update the headers of ISPI images with FITS keywords describing the image coordinate system, its projection on the sky, etc. Install the WCSTOOLS software, and read the information about imwcs [35] which is the program used to compute a coordinate solution and update the image headers.
  • Catalogs. WCSTOOLS needs an astrometric reference catalog. These can be existing catalogs like 2MASS or user catalogs. The best way to proceed is to download and install one or more catalogs locally. At CTIO, we have so far experimented with 2MASS, since this is very useful for near infrared observations in the Galactic plane. Other catalogs will work well for regions where the optical extinction is not high. To install the 2MASS catalog, follow the instructions given by the WCSTOOLS pages [36]. The disk space requirement is about 40 GB.
  • A program to find image centroids in the ISPI images. This can be accomplished in IRAF (for example daofind). Other well known programs are DOPHOT, DAOPHOT, and SExtractor [37]. Finding accurate centroids is necessary for certain options using WCSTOOLS.
  • A program to resample or interpolate images. We have used SWARP [38], a program by Emmanuel  Bertin.

Updating the Basic Headers for a WCS.

  • The images should be preprocessed to have a near zero rotation from some known position. For ISPI, the usual N up E left projection is easily accomplished by choosing the "flipy" flag in the osiris.cl IRAF meta script. ISPI raw images have N down and East left when displayed in IRAF. An EPOCH keyword should also be added to the header which matches the value of the ISPI EQUINOX keyword (the default is 2000.0).
  • The simplest task needed for raw ISPI images is to add appropriate world coordinate system (WCS) keywords [39]. The program imwcs should be used to do this. There are a variety of options, but a simple call would be the following:

    > imwcs -c tmc -d t150s.dao -h 200 -n 6 -vw -p 0.303 t150s.fits

    In this example, -c tmc calls for the 2MASS catalog. The location of the catalog is specified during the build process of WCSTOOLS. -h 200 calls for 200 stars to be used in the matching between the image and catalog, and -n 6 refers to the order of the fit between image and catalog positions (6 parameters in each axis). -p 0.303 gives the initial image scale in arcseconds per pixel for ISPI; this is required since this information is not in the header. Alternately, one could add SECPIX (with a value of 0.3) to the ISPI image headers. -d t150s.dao specifies that the file "t150s.dao" be used for the input x,y positions of the stars to be matched by the catalog. This file should be sorted on brightness so the brightest stars are used. The file should be one star per line, x, y, and mag. "#'s" indicate comments and are ignored in this file, as are any fields after the first three of each line. See the IMWCS Command Line Arguments [40] for more information. -vw sets the "verbose" flag along with a call (w flag) to write a new image (in this case t150sw.fits) which has the updated WCS. See the imwcs command line argument -o for more output options, including updating the input image header.

  • The result of IMWCS is an image updated (if desired) with a WCS which is a standard tangent plane projection on the sky. In the above example, a new image t150sw.fits was created from t150s.fits. The pixel data are not changed, only the header.
  • Once this step has been performed, use DS9 (SAO image) to display the updated image. The WCS coordinates should appear in the upper left of the image display. WCS specific IRAF scripts or any other program which is based on the FITS standard should also now be able to read and use the WCS information stored in the updated ISPI header.

Taking Distortion into Account

  • The above call to imwcs results in a new image header, but it is often desired to actually resample or interpolate the image to a new projection. This could just be to have the pixel to sky projection modified or to allow for correcting higher order distortions. The latter is the case for ISPI since these distortions will affect the combination of dithered images. ISPI exhibits distortion, so combining dithered images without correcting for this will result in a degraded PSF at positions away from the center of the image. The effect can be significant (of order several pixels change in FWHM.
  • The key to creating a plate solution which accounts for the ISPI distortions is to have an image or images with many stars across the field of view and a set of accurate catalog positions.
  • Once a matched catalog of pixel and sky coordinates is in hand, any program which fits the latter to the former and writes an appropriate solution can be used to further update the image header. In practice, our solution relies on the IRAF task "ccmap" because this task will do the fit and add the appropriate keywords to the header. This is done in the IRAF TNX [41] system. The advantage is that the SWARP [38] program used to resample the dithered ISPI images automatically handles the TNX system (which is not officially part of the FITS standard). TNX is a combination of the linear WCS solution plus higher order terms. The fit is a general polynomial in x and y with cross terms. Initial testing with ISPI suggests ccmap fits of 4th-6th order (in IRAF speak) with full cross terms do well.

In this image, the residuals to a linear fit are shown. The peak errors in position would be about three pixels if the distortion were not accounted for. In the next image, the residuals for a 4th order fit are shown. The RMS is less than 0.1'' which is about 1/3 of an ISPI pixel.

  • The ccmap solution is written to a database file which IRAF understands and also to the image header if the flag "update" is turned on (i.e. set to yes) and images are specified in the images parameter of the task. The image header keywords are explained in the TNX description [41].
  • The ccmap solution can easily be applied in SWARP. Refer to the SWARP pages for details on how to execute this program. A simple call to swarp would be:

    > swarp -c swarp.gc t150sw.fits t151sw.fits t152sw.fits t153sw.fits t154sw.fits

    In this example, the swarp config file is "swarp.gc" and five images are reinterpolated. Each image can have its own TNX solution, which in our example would be the combination of WCS and TNX keyword values produced by imwcs and ccmap. Unless the fit to each image is very good, it is advisable to use an average solution for the higher order terms (and rotation and scale in this case). This can easily be done in swarp by specifying an auxillary header file with the FITS keywords and values. Cutting and pasting from the above solution from the FITS header, a file of the form t15?sw.head was made for each image. The file has the following form and is identical for each image:

    WAT0_001= 'system=image'
    WAT1_001= 'wtype=tnx axtype=ra lngcor = "3. 4. 4. 1. -0.09223848232288634 0.095'
    WAT2_001= 'wtype=tnx axtype=dec latcor = "3. 4. 4. 1. -0.09223848232288634 0.09'
    CD1_1   =  -8.4729762794223E-5
    CD1_2   =  7.23448951262463E-7
    CD2_1   =  7.40672424953437E-7
    CD2_2   =  8.47678406227307E-5
    WAT1_002= '30461739012073 -0.08294874238596413 0.086260218248302 -3.81821493608'
    WAT1_003= '0989E-6 0.004201053595244364 0.004530606959289172 -0.722711400721411'
    WAT1_004= '9 -5.887103279160741E-4 0.006853642315098468 0.0766861807819218 0. -'
    WAT1_005= '0.005316483436707363 -0.4617939442868234 0. 0. 0.1156522538399243 0.'
    WAT2_002= '530461739012073 -0.08294874238596413 0.086260218248302 4.34379162931'
    WAT2_003= '0940E-6 -5.007423412969465E-4 6.003889063455503E-5 0.084985861670387'
    WAT2_004= '73 0.003355387441743075 0.003754664665905738 -0.3911461579926636 0. '
    WAT2_005= '0.001750656256554302 0.04350355226749533 0. 0. -0.5986827241716564 0'
    WAT1_006= ' 0. 0. "'
    WAT2_006= '. 0. 0. "'
    END

    These lines were copied directly from the fits header of the image used in our example to find the higher order solution. A tab follows the keyword "END."

  • Swarp will combine the dithered images using the individually resampled images (see the swarp documentation and the config file options). The resampling includes rotation, scale, distortion and the standard pixel to sky projection (in this case, the TAN or tangent plane projection).

    Bad pixels may be rejected in the combine process by using the "pixel weighting" capability in SWARP. A simple weighting which rejects bad pixels can be accomplished by using the "mask2.fits" file produced by maskbad.cl (see the cirred package). This file has a 1 for each good pixel and 0 for each bad one. The file is called in the swarp configuration file by specifying its name in the WEIGHT_IMAGE keyword and using WEIGHT_MAP as the weighting method (it appears that one must specify this file name for each input image. I.e, make a coma separated list for the value of WEIGHT_IMAGE with the same file name repeated to match the number of input images). See the swarp documentation for more information. A slightly more sophisticated weight map can be made by multiplying the normalized flatfield image by mask2. Finally, one can specify a unique weight mask for each image through the file "imagename.suffix.fits" where "suffix" is given in the swarp config file, and the default is "imagename.weight.fits."

End to End Processing of ISPI Images

The following is an example of how ISPI images may be fully processed using the tools discussed on this page (cirred, WCSTOOLS, SWARP, IRAF).

  1. Make flats: Use med.cl to combine lights on and off images into a domeon and domeff image, respectively. Subtract the off image from the on.
  2. Make bad pixel masks: Use maskbad.cl to create bad pixel masks.
  3. Fix flats: Use fixfits.f to fix the flat images. The resulting images will have bad pixels interpolated over. If the images are named "dFilterNamef.fits", then the osiris.cl script will apply the right flat according to the FITS header filter name where FilterName is th evalue of this header parameter (e.g. KA for the ISPI Ks filter).
  4. Run osiris.cl: The meta script which will preprocess images, divide by a flat, and add appropriate head key words.
  5. Sky subtract the images: Use sky_sub.cl to make and subtract sky frames.
  6. Add a WCS: Use do_wcs.cl to add a World Coordiante System to the images. (Preliminary script included in cirred package. Requires WCSTOOLS [34]. Requires a catalog implementation for the WCSTOOLS.).
  7. Compute high order distortions: Use do_ccmap.cl to compute the high order distortion coefficients and add them to the image header (Preliminary script included in cirred package. Requires WCSTOOLS. Need to incorporate TNX keywords manually as described above).
  8. Shift/Combine images: Use SWARP [38] (by Emmanuel Bertin) to reinterpolate the images and optionally shift and combine them.
  9. Make a color image: SWARP can be used to shift and orient the final J, H, and K images so they can be combined in a false color image.

Acknoledgements: We would like to acknowledge the useful comments and input of Silas Laycock regarding ISPI data reduction and astrometry. It is also a pleasure to acknowledge the very useful software and documentation of the Mink WCSTOOLS, the Bertin SWARP program, and the IRAF suite.

Please mail questions or comments to:
rblumATctio.noao.edu
nvdbliekATctio.noao.edu
bgregoryATctio.noao.edu

Guiding with ISPI

Typically, ISPI observations are short and unguided. This simplifies the execution of large mosaics. If dithering about in a restricted area for a long time, then guiding can be useful to improve or maintain image quality. The guide field is provided by the portion of the telescope field which is viewed "past" the tertiary mirror (M3) which delivers the science field to ISPI by a 45 deg reflection. M3 and the guider both live in the white rotator box fixed to the back of the primary mirror cell.

The M3 mirror as viewed from above. The mirror is actually an elliptical shape but appears round since it is mounted 45 deg to the optical axis. The guide field is to the left in this image.

This sketch shows the available guide field using the F/8 guider. The inner radius is 15.35 arcmin and the outer radius is 23.03 arcmin as measured from the center of the ISPI field. Guide stars should be chosen toward the center of the annulus if possible (i.e. between R3 and R4). This will allow for the guide stage to follow the offsets for dithers of several arcmin. If you need to offset further (e.g. to take separate sky frames), make sure to have the night assistant turn off the guiding before the offset is made. The scale is 6.58 arcsec/mm.

Last Update 21 April, 2005; Bob Blum

ISPI Troubleshooting & known problems

The most common problem when observing with ISPI, is that the ISPI GUI will loose connection with the TCS. As a result the ISPI GUI might hang. If this happens, close the GUI, stop the ArcVIEW main application and then restart the GUI by clicking on the ArcVIEW icon in the VNC window. Detailed instructions can be found under Simple Restart on the Re-starting ISPI page [42].

If simply restarting the GUI does not bring back up the connection with the TCS, ask the night assistant to re-start the TCS router. See also Telescope & TCS related problems [43]

Limit the number of processes running on ctioa9 to a minimum.
This will improve overall stability of the system.

Further information on problems one might encounter when observing with ISPI, can be found in the following pages:
    ISPI PC & ISPI GUI related problems
    Telescope & TCS related problems
    Re-starting ISPI
    Noise & noise patterns
    Condensation on entrance window

February 21, 2007

ISPI Noise & noise patterns

A simple receipe to measure the noise:

  1. Collect 10 darks at shortest exposure (3.227s)
     
  2. Subtract last 3-4 pairs from each other
    Dark10-Dark9, Dark9-Dark8, etc
     
  3. imstat with lo=-100 hi=100, to get rid of outlayers
     
  4. Divide average std number by sqrt(2) and multiply by 4.25e-/ADU (gain)

    • If number is less than 20e-rms, the noise is within specs.
    • If the number is more than 20e-rms, ask moutain support to

      • check that the cables and ground strap are well-tighed and perhaps ask them to move the cables.
      • Switching the motors for the filterwheels OFF & ON might help too.
  1. Measure the noise again, and again...

 

Site under construction

February 21,2007

ISPI Problems related to the TCS

TCS Time Outs

Occasionally, the connection between the TCS and ISPI will time out. The yellow "busy" light on the TCS section of the main GUI will light up (normal), but not turn off, indicating that the TCS command was sent but no response returned to ISPI. In this case, several things can happen. Some intervention by the user may be required.

  • The GUI will try and restart the connection automatically. This is only done if the check box "TCS auto recovery" is checked.
  • If auto recovery is off or does not respond, click on the "TCS GUI" button in the upper left area of the main GUI. This will bring up the separate TCS GUI. Note any error messages in the TCS GUI and add these to your night report. To reconnect, click on "Stop" in the upper right of the GUI. It may be necessary to click on "Stop" a second time. The GUI should flash and a white arrow appear in the upper left corner. Click on this arrow. The green TCS light should come on. The response field should indicate that the connection is back up. Click on "hide" to close the TCS GUI.
  • It may take 10 seconds or so for the "TCS auto recovery" to execute. If one does not want to wait, the TCS GUI can be restarted manually.
  • If there is ever any doubt about the status of the TCS connection, bring up the TCS GUI and click on "info." A current string of TCS info should come over the link. If the normal TCS info is reported in the GUI response window, then communication is ok.
  • As with all the popup GUI's, always open and close with the "show" and "hide" buttons. These buttons can either be clicked on the popup or the main GUI.

Other TCS Errors

Sometimes there will be a TCS error which is not a time out. That is, for whatever reason, a TCS command sent from ISPI will not be completed. In this case, ISPI will also stop, for example, during a script. The TCS error must first be solved, then the ISPI command or script restarted. Have the night assistant note any errors which may have appeared on the TCS console or on the TCS router.

TELESCOPE PROBLEMS

If you have other problems during the night, please have the night assitant call for assitance.

Detail any problems in the nightly report form [44], and ask the night assistant to fill in a trouble report (GNATS [45]) for serious problems.

It is good practice to put in the DARK filters in the filter wheel in case telops needs to turn on the lights in the dome to fix the problem. You can do this either by moving to the DARK position on the filter wheels or just clicking on the END OF NIGHT button.

February 15, 2007

Problems related to the ISPI GUI

GUI hangs because connection with TCS is lost

The most common problem when observing with ISPI, is that the ISPI GUI will loose connection with the TCS. As a result the ISPI GUI might hang. If this happens, close the GUI, stop the ArcVIEW main application and then restart the GUI. Detailed instructions can be found under Simple Restart on the Re-starting ISPI [46]page.

If simply restarting the GUI does not bring back up the connection with the TCS, ask the night assistent to re-start the TCS router. See also the pages on Telescope & TCS related problems.

Overloading ctioa9

The GUI might also hang as a result of an overload on ctioa9, the ISPI computer. Avoid opening too many terminal windows on the VNCs running on ctioa9. It is recommended to limit the number of windows to the minimum, e.g. one to run IRAF and one to edit script. Use either of these windows to check the temperature.

An apparent problem, where the GUI hangs when using too many co-adds, is most likely related to an overload on ctioa9. We have not been able to reproduce this problem. Hence, we advice to limit the number of processes run on ctioa9, rather than limiting the number of coadds.

Along the same line of thought: when transferring data to another computer, it is recommended to start the transfer at the other end, rather than running the command, e.g. scp, on ctioa9.

Other problems with ISPI GUI & ctioa9

....

NOTE: Site under construction

February 21, 2007

Re-starting ISPI

If the GUI hangs, you will need to restart ISPI, you should start from the simplest restart, building to a more full restart as neccessary.

Simple Restart:

  • Single click the red 'EXIT' button on the ISPI Control GUI (with the really cool cross section view of ISPI).
  • Quit the  'ArcVIEW Main Application' window by first clicking on the 'EXIT' button in the window and then pulling down the 'File' menu in the top, left corner of the window, and going to 'Quit'.
  • Single click on the 'ArcVIEW' button in the vnc session xload box. This will pop up the ISPI Control GUI, and an 'ArcVIEW Main Application' window.  The Application window should have two round buttons, labelled 'Connected?' that appear red.
  • On the ISPI Control GUI, there is a white arrow icon in the upper left corner of the window. Single click on the white arrow. This arrow should now turn black. The two buttons in the Application window should turn green.
  • Minimize the Application window by clicking on the underscore '_' icon on the blue bar, in the upper right corner of each window. You are now ready to operate ISPI.

Less simple:

  • Quit as above for the simple restart. Now check to see if there are any remaining labview processes. To do this, from a terminal window, connected to ctioa9 type the command "ps -ef | grep labview". Note the process# attached to the labview processes.
  • Kill all the remaining labview processes with the command "kill -9 process#" from the same terminal window.  You really only need to kill the parent process in order to kill all the subprocesses.  The parent is the process with a 1 next to the process#.
  • Verify with ps -ef | grep labview that all the labview processes are killed.
  • Now continue to re-start the GUI as above in the simple restart.

Frustration level rising, deep breath, restart

  • Do as in the less simple re-start above. Kill the "TightVNC: ispi's X desktop (ctio:9)"  by clicking on  the X in the upper, right corner of the window. THIS IS THE ONLY X IN A BOX THAT YOU SHOULD EVER CLICK ON IN  ANY PART OF THE VNC DESKTOP!
  • Click the ISPI VNC icon in the upper left of the ctioa8 desktop.
  • Single click on the 'ArcVIEW' button in the vnc session xload box. This will pop up the ISPI Control GUI, and an 'ArcVIEW Main Application' window.  The Application window should have two round buttons, labelled 'Connected?' that appear red.
  • On the  ISPI Control GUI, there is a white arrow icon in the upper left corner of the window. Single click on the white arrow. This arrow should now turn black. The two buttons in the Application window should turn green.
  • Minimize the Application window by clicking on the underscore '_' icon on the blue bar, in the upper right corner of each window. You are now ready to operate ISPI.

And finally:

  • Call the mountain electronics support to reboot ctioa9 and/or power cycle the instrument at the telescope.
  • After rebooting, follow the normal startup [47].

If in doubt, or really stuck, phone the instrument scientists:

  • Contact scientist cell phone: 09-499-6216
  • Nicole van der Bliek: home: 211869; office: 205
  • Malcom Smith: office: 217
     

February 15, 2007

ISPI contact information

For further information on ISPI, please contact

Nicole van der Bliek
Malcolm Smith

or one of the members of the ISPI project team.

February 7, 2007

ISPI Project Team

  • Project Scientist - Ron Probst
  • Project Manager - Nicole S. van der Bliek
  • Mechanical Design - Andres Montane
  • Optics - Maxim Boccas & Roberto Tighe
  • Software - Marco Bonati, Francisco Delgado & Mike Ashe
  • Electronics - Michael Warner
  • Vacuum and Cryogenics - Ramon Galvez
  • Filter Wheel Electronics - Manuel Martinez

March 23, 2003

ISPI Publications

Refereed Publications using Infrared Side Port Imager (ISPI) on the Blanco 4-m Telescope

(last updated February 2012)

  • Spectrophotometric distances to Galactic H II regions (Moises et al. 2011, MNRAS, 411, 705) [48]
  • The Ultra Cool Brown Dwarf Companion of WD 0806-661B: Age, Mass, and Formation Mechanism (Rodriguez et al. 2011, ApJ, 732, L29) [49]
  • A Wide-Field Survey of the Orion Nebula Cluster in the Near-Infrared (Robberto et al. 2010, AJ, 139, 950) [50]
  • LoCuSS: first results from strong-lensing analysis of 20 massive galaxy clusters at z = 0.2 (Richard et al. 2010, MNRAS, 404, 325) [51]
  • On a New Near-Infrared Method to Estimate the Absolute Ages of Star Clusters: NGC 3201 as a First Test Case (Bono et al. 2010, ApJ, 708, L74) [52]
  • Properties of RR Lyrae stars in the inner regions of the Large Magellanic Cloud. III. Near-infrared study (Borrisova et al. 2009, A&A, 502, 505) [53]
  • A Public, K-Selected, Optical-to-Near-Infrared Catalog of the Extended Chandra Deep Field South (ECDFS) from the Multiwavelength Survey by Yale-Chile (MUSYC) (Taylor et al. 2009, ApJS, 183, 295) [54]
  • Spitzer View of Young Massive Stars in the Large Magellanic Cloud H II Complex N 44 (Chen et al. 2009, ApJ, 695, 511) [55]
  • High- and Intermediate-Mass Young Stellar Objects in the Large Magellanic Cloud (Gruendhl & Chu, 2009, ApJS, 184, 172) [56]
  • The Multiwavelength Survey by Yale-Chile (MUSYC): Wide K-Band Imaging, Photometric Catalogs, Clustering, and Physical Properties of Galaxies at z~2 (Blanc et al. 2008, ApJ, 681, 1099 [57]
  • The First Release COSMOS Optical and Near-IR Data and Catalog (Capak et al. 2007, ApJS, 172, 99) [58]

If we have missed one of your publications which used ISPI, please contact one of the members of the ISPI instrument team [59]


Source URL (modified on 07/14/2020 - 12:01): http://www.ctio.noao.edu/noao/content/ISPI

Links
[1] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/block.png
[2] http://www.ctio.noao.edu/noao/content/ISPI-Filters
[3] http://www.astro.ufl.edu/~elston/flamingos/flamingos.html
[4] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/jspot.jpg
[5] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/hspot.jpg
[6] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/kspot.jpg
[7] http://www.janostech.com/
[8] http://www.astro.ufl.edu/~elston/flamingos/first.light.gif
[9] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/ispigui.jpg
[10] http://www.ctio.noao.edu/instruments/ir_instruments/image_cal.html
[11] http://www.ctio.noao.edu/noao/content/IR-Instruments
[12] http://www.ctio.noao.edu/noao/content/ir-instruments
[13] http://www.ctio.noao.edu/noao/content/ISPI-contact-information
[14] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/exampleH.dat
[15] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/exampleJ.dat
[16] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/exampleK.dat
[17] http://www.ctio.noao.edu/noao/content/ispi-quick-reference-guide
[18] http://www.ctio.noao.edu/noao/content/ispi-linearity-correction
[19] http://www.ctio.noao.edu/noao/content/ispi-flat-fields
[20] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/pointK.dat
[21] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/pointKsky.dat
[22] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgKx_res
[23] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgKxsky_res
[24] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgKy_res
[25] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgKysky_res
[26] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/pointH.dat
[27] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/pointHsky.dat
[28] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgHx_res
[29] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgHxsky_res
[30] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgHy_res
[31] http://www.ctio.noao.edu/noao/sites/default/files/instruments/imagers/avgHysky_res
[32] http://www.ctio.noao.edu/instruments/ir_instruments/cirred/cirred.html
[33] http://www.ctio.noao.edu/instruments/ir_instruments/datared.html
[34] http://tdc-www.harvard.edu/software/wcstools/
[35] http://tdc-www.harvard.edu/software/wcstools/imwcs/
[36] http://tdc-www.harvard.edu/software/catalogs/tmc.convert.html
[37] http://terapix.iap.fr/rubrique.php?id_rubrique=91/
[38] http://terapix.iap.fr/rubrique.php?id_rubrique=49
[39] http://tdc-www.harvard.edu/software/wcstools/wcstools.wcs.html
[40] http://tdc-www.harvard.edu/software/wcstools/imwcs/imwcs.com.html
[41] http://iraf.noao.edu/projects/ccdmosaic/tnx.html
[42] http://www.ctio.noao.edu/noao/content/re-starting-ispi
[43] http://www.ctio.noao.edu/noao/content/ISPI-Problems-related-TCS
[44] http://www.ctio.noao.edu/new/Tools/Forms/EON/Form.php?telescope=Blanco%204-m
[45] http://www3.ctio.noao.edu/cgi-bin/gnatsweb.pl
[46] http://www.ctio.noao.edu/noao/content/Re-starting-ISPI
[47] http://www.ctio.noao.edu/noao/content/ISPI-Start-Procedure
[48] http://adsabs.harvard.edu/abs/2011MNRAS.411..705M
[49] http://adsabs.harvard.edu/abs/2011ApJ...732L..29R
[50] http://adsabs.harvard.edu/abs/2010AJ....139..950R
[51] http://adsabs.harvard.edu/abs/2010MNRAS.404..325R
[52] http://adsabs.harvard.edu/abs/2010ApJ...708L..74B
[53] http://adsabs.harvard.edu/abs/2009A%26A...502..505B
[54] http://adsabs.harvard.edu/abs/2009ApJS..183..295T
[55] http://cdsads.u-strasbg.fr/abs/2009ApJ...695..511C
[56] http://adsabs.harvard.edu/abs/2009ApJS..184..172G
[57] http://adsabs.harvard.edu/abs/2008ApJ...681.1099B
[58] http://adsabs.harvard.edu/abs/2007ApJS..172...99C
[59] http://www.ctio.noao.edu/noao/content/ispi-contact-information