; IRING ;--------------------------------------------------------------- ;! integrates intensity / flux in rings / ellipses ;# TASK ANALYSIS SPECTRAL PLOT ;----------------------------------------------------------------------- ;; Copyright (C) 1995-1996, 1999-2001, 2008-2009, 2013, 2015 ;; 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 ;----------------------------------------------------------------------- IRING LLLLLLLLLLLLUUUUUUUUUUUU CCCCCCCCCCCCCCCCCCCCCCCCCCCCC IRING: Task to integrate in confocal ellipses. USERID -32000.0 32000.0 User ID - ignored INNAME Image name (name) INCLASS Image name (class) INSEQ 0.0 9999.0 Image name (seq. #) INDISK 0.0 9.0 Disk drive # BLC 0.0 4096.0 Bottom left corner of image 0=>entire image TRC 0.0 4096.0 Top right corner of image 0=>entire image APARM Orientation parameters: 1: x - coordinate center wrt reference point in pixels 2: y - coordinate center wrt reference point in pixels 3: position angle major axis, degrees 4: inclination of rings to plane of sky, deg. 5: full range of azimuthal angles to consider, deg. 0 => 360 degrees 6: central azimuthal angle of the range specified in APARM(5), deg. CPARM Specification of rings: 1 : radius inner ring (") 2 : radius outer ring (") 3 : width of each ring (") 4 : 1=>treat blanks as zero 5 : 1=>plot radial intensity 2=>plot integrated radial intensity, from inner ring outwards. 6 : 1 => CPARM(1-3) in pixels 7 : scale X axis by CPARM(7) 8 > 0 => y limits in 9, 10 KEYSTRNG X axis type if CPARM(7) > 0 SYMBOL -24.0 24.0 Plot symbol number: 1 - 24 FACTOR -20.0 20.0 Scale plot symbols by FACTOR <0 => connect points too DOTV -1.0 1.0 > 0 Do plot on the TV, else make a plot file GRCHAN 0.0 8.0 Graphics channel 0 => 1. OUTTEXT Not blank => write list x,y plot points to text file = 'NONE' -> no text file and no message file of values ---------------------------------------------------------------- IRING Type: Task Use: Integrates a map in concentric annulli projected onto the sky according to the orientation parameters in APARM. As well the mean value per ring is given, as the total cumulative flux in the right units. The latter is only possible if information about the (Gaussian) beam is present in the header. The RA-dec coordinates should be the first two axes. Error bars are given using the rms of the pixels in the annulus divided by the square root of the number of pixels counted in the annulus scaled by an estimate of the correlation of those pixels. If the annulus is wider than the Clean beam minimum then that estimate is the beam area in pixels. If the annulus is smaller, then that estimate is the beam area divided by the ratio of the Clean beam minimum to the width of the annulus. The beam area and the correlation correction in pixels are given. If there is no Clean beam in the header, the errors in the mean values are computed assuming the noise in each pixel is independent. Adverbs: USERID......User ID of owner of image. Ignored INNAME......Image name(name). blank=>any INCLASS.....Image name(class). blank=>any INSEQ.......Image name(seq. #). 0=>any INDISK......Disk drive # of image. 0=>any BLC.........The Bottom Left-hand pixel of the subarray of the image to be analysed. The value (0,0) means (1,1). TRC.........The Top Right-hand pixel of the subarray of the image to be analysed. The value (0,0) means the top right hand corner of the entire image. APARM.......Orientation parameters : (1) : x pixel central position relative to the image reference point in pixels. (2) : y pixel central position relative to the image reference point in pixels. APARM(1) AND (2) are added to the pixel position of the reference point (use IMHEADER to see those pixel coordinates) to obtain the pixel coordinate of the center of the ellipses. (3) : position angle (degrees) receding major axis; East from North. (4) : inclination (degrees) of orbital plane to plane of sky. A face on galaxy is thus at inclination 0 degrees. (5) : full-width of range of azimuthal angles (degrees) over which distributions should be computed 0 => use all azimuthal angles (6) : central azimuthal angle (degrees) for the range specified in APARM(5); East from North. CPARM.......(1) : Radius inner circle of first annulus ("), default is 0 (2) : Radius outer circle of last annulus ("), defaults to min(mapsize x,mapsize y) (3) : Width of each annulus, defaults to 10 annuli (4) : 0 => disregard blanks, 1 => treat blanks as zero (5) : 1 => plot radial intensity distribution. : 2 => plot summed intensity from first annulus outwards. (6) : 1 => CPARM(1-3) in pixels. Use this for non-astronomical images although it could be used also for those images that have no Clean beam. The pixel values are assumed independent when computing the errors in the mean. (7) : >0 => scale the X axis by CPARM(7) and use KEYSTRNG as the label (8) : >0 => use CPARM(9) and CPARM(10) as the plot value range on the Y axis. (9) : Minimum value in Y plotted (CPARM(8) > 0) (10) : Maximum value in Y plotted (CPARM(8) > 0) KEYSTRNG...Label to use on X axis when CPARM(7) > 0 SYMBOL.....1: Plus sign (default) 12: Five pointed star 2: Cross (X) 13: Star of David 3: Circle 14: Seven-pointed star 4: Box 15: Eight-pointed star 5: Triangle 16: Nine-pointed star 6: Diamond 17: Ten-pointed star 7: Pentagon 18: 11-pointed star 8: Hexagon 19: 12-pointed star 9: Septagon 20: 13-pointed star 10: Octagon 21: 14-pointed star 11: Nine-gon 22: Plus with gap 23: vertical line 24: Cross with gap SYMBOL <= 0 => do not plot error bars. FACTOR.....Scale plot symbols by FACTOR 0 -> 1. If FACTOR < 0, draw lines between the points and use abs(FACTOR) to scale the points plotted. DOTV........> 0 => plot directly on the TV device, otherwise make a plot file for later display on one or more devices (including the TV if desired). GRCHAN......Graphics channel (1 - 8) to use for line drawing. 0 => 1 OUTTEXT.....Not blank => write a list of the X (radius) and Y (mean or sum) plot points and the uncertainties in the mean or sum to a text file specified, in the usual manner, by OUTTEXT. If OUTTEXT='NONE', no text file is written and the values found are not written to the message file. Any other value of OUTTEXT will have the ring radii, number of non-zero pixels, mean, sigma, and cumulative values displayed to the message file for every anulus. ----------------------------------------------------------------