;+ ; froi.pro IDL procedure ; ; PURPOSE ; Read FITS images and define a region of interest. ; ; SYNTAX ; froi, 'fits_file' ; ; PROCEDURES USED ; readfits ; writefits ; fxaddpar ; ; HISTORY ; Author: Andrew L. Stanger 2 June 1998 ;- pro froi, fits_file loadct, 0 lct,'/home/stanger/color/bwy.lut' red = indgen (256) green = indgen (256) blue = indgen (256) xdim = 0 ydim = 0 xcen = 0 ycen = 0 roll = 0.0 mag = 1.0 ;--- Read FITS image. img = readfits (fits_file, hdu) xdim = fxpar (hdu, "NAXIS1") ydim = fxpar (hdu, "NAXIS2") nbypix = fxpar (hdu, "BITPIX") * 8 xcen = fxpar (hdu, "CRPIX1") ycen = fxpar (hdu, "CRPIX2") pixrs = fxpar (hdu, "PIXRS") ;--- Display image. tvscl, img ;--- Define Region-of-interest (ROI). polypix = defroi (xdim, ydim, xv, yv) size_roi = size (polypix) nloc = size_roi (1) xloc = intarr (nloc) yloc = intarr (nloc) FOR i = 0, nloc-1 DO $ BEGIN yloc (i) = polypix (i) / xdim xloc (i) = polypix (i) - yloc (i) * xdim END size_xv = size (xv) nv = size_xv (1) print, 'size (xv): ', size (xv) print, 'size (yv): ', size (yv) FOR i = 0, nv - 1 DO $ print, 'i, xv, yv: ', i, xv (i), yv (i) print, 'size (polypix): ', size (polypix) FOR i = 0, nloc - 1 DO $ print, 'polypix, xloc, yloc: ', polypix (i), xloc (i), yloc (i) xmass, imgref, imgcme, xdim, ydim, nbypix, $ xcen, ycen, pixrs, xloc, yloc, nloc END