; TVMAXFIT
;---------------------------------------------------------------
;! displays fit pixel positions and intensity at maxima on TV
;# Procedure TV Analysis
;-----------------------------------------------------------------------
;;  Copyright (C) 1995, 2003, 2009
;;  Associated Universities, Inc. Washington DC, USA.
;;
;;  This program is free software; you can redistribute it and/or
;;  modify it under the terms of the GNU General Public License as
;;  published by the Free Software Foundation; either version 2 of
;;  the License, or (at your option) any later version.
;;
;;  This program is distributed in the hope that it will be useful,
;;  but WITHOUT ANY WARRANTY; without even the implied warranty of
;;  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
;;  GNU General Public License for more details.
;;
;;  You should have received a copy of the GNU General Public
;;  License along with this program; if not, write to the Free
;;  Software Foundation, Inc., 675 Massachusetts Ave, Cambridge,
;;  MA 02139, USA.
;;
;;  Correspondence concerning AIPS should be addressed as follows:
;;         Internet email: aipsmail@nrao.edu.
;;         Postal address: AIPS Project Office
;;                         National Radio Astronomy Observatory
;;                         520 Edgemont Road
;;                         Charlottesville, VA 22903-2475 USA
;-----------------------------------------------------------------------
TVMAXFIT  LLLLLLLLLLLLUUUUUUUUUUUU CCCCCCCCCCCCCCCCCCCCCCCCCCCCC
TVMAXFIT: Procedure to find coordinates and values on the TV
IMSIZE          3.0        16.0    Search area for peak value.
PRTLEV         -2.0        10.0    <= -1 => no display results
         @         Output adverbs
INNAME   @                         Image name (name)
INCLASS  @                         Image name (class)
INSEQ    @                         Image name (seq. #)
INTYPE   @                         Image physical type ('MA')
INDISK   @                         Disk number
PIXXY    @                         Pixel location (X,Y,...) in
         @                         image - last chosen
PIXVAL   @                         Intensity at PIXXY.
COORDINA @                         Coordinates at pixel
ERROR    @                         > 0 => fit failed, else okay
TVBUT    @                         Value of last button = 8.

          Instructions will also appear on terminal.
----------------------------------------------------------------
TVMAXFIT
Type: procedure
Use:  TVMAXFIT is a simple procedure that loops reading TV
      cursor positions with IMXY and returning the image value
      and coordinates at an maximum near that location to the
      terminal and message file with MAXFIT.  Hit button A, B,
      or C to select a location.  Hit button D to exit.

      MAXFIT fits a quadratic function to a 3x3 map array to
      determine the position and strength of an extremum.  In
      this procedure, the cursor position when the button is
      hit is used as an estimate of the center of the 3x3
      array.  The fit will search over an area of IMSIZE pixels
      for the actual peak value before fitting.
Adverbs: (inputs)
  IMSIZE......(X,Y) search area in pixels.  Verb uses 11 if
              the values are outside the range 3 to 16.
  PRTLEV......Print level: if <= -1, then do not put result info
              in message file, otherwise do the coordinate and
              flux displays.
Adverbs: (outputs)
  INNAME......Image name(name).
  INCLASS.....Image name(class).
  INSEQ.......Image name(seq. #).
  INDISK......Disk drive # of image.
  PIXXY.......Pixel location (X,Y,...) in the image.  The last
              one selected is returned.
  PIXVAL......Intensity at this pixel.
  COORDINA....The X and Y coordinates are found as:
              Xpos = abs(CO(1)) + abs(CO(2))/60 + abs(CO(3))/3600
                 if any of CO(1), CO(2), CO(3) < 0: Xpos = -Xpos
              Ypos = abs(CO(4)) + abs(CO(5))/60 + abs(CO(6))/3600
                 if any of CO(4), CO(5), CO(6) < 0: Ypos = -Ypos
              Note that, although these are most suited to RA, Dec in
              sexagesimal notation, they can be used for any type of
              coordinate.  The units are standard FITS units (e.g
              degrees, m/sec, Hz, sec, etc.) except that right
              ascensions are in hours of time.
  ERROR.......If the verb is not aborted for error, but the fit
              fails, ERROR will be TRUE (1), else it will be
              FALSE (-1).  Procedure exits if ERROR is TRUE.
----------------------------------------------------------------
See EXPLAIN MAXFIT.
----------------------------------------------------------------
