Generate a local FITS image (named "rad.fits") of elliptical radii computed by projecting a (tilted) circle.

%  ellradius.sh 1992 1992 123.0 300.0 20.0 0.75 25.0   

 arg1 - X dimension of output image
 arg2 - Y dimension of output image
 arg3 - X center in pixels 
 arg4 - Y center in pixels 
 arg5 - circle radius (asemi) in pixels) 
 arg6 - axis ratio (b/a) 
 arg7 - pa in degrees (GALPHOT convention)

FIGURE 1: Radius maps from ellradius. In every case above an ellipse with b/a=0.35 and a semi-major axis length of asemi=150 pixels was generated with different pa (Position Angle) values. The value of each pixel in the output radius map (rad.fits) is the radius in the deprojected (original) circle. The color scale at the bottom of the figure shows that the maximum radius values run from 0 to about 150 pixels. The pa values are measured from the POSITIVE-Y axis in a CLOCKWISE direction. For comparison we also show the corresponding ds9 Angle value used to plot a vector representing the semi-major axis of the projected ellipse.

Use the gal2dsubs_mark2.f (SCOLIB) routine called pa_2_angleds9.f to convert the pa (in degrees) to ds9 angle value:

 
 
      call pa_2_angleds9(pa,angleds9)