sp - Make a scatter plot comparing two images RWW & MWP
EXAMPLES:
sp 1,2,1 Make a scatter plot of image 2 vs image 1 within
area 1. The results will be in the file "scatter"
in your home directory.
sp an:3 fn:"scatter3" yfact:1E-4 minx:.5
Make a scatter plot with file named scatter3 in the
current directory using area 3 multiplying the
y values by 3 and suppressing data with x < .5
(unscaled).
sp logx: logy: miny:0 minx:0
Make a log-log plot of the images suppressing negative
values.
sp ps:2 fill: ovl: Plot with filled triangles instead of dots and overlay
on the current scatter plot. See 'gm ?!' for a list of
plot symbols and their figure numbers. See also note on
plot symbols below.
DESCRIPTION:
Sp will produce a plot comparing two images using comb graphics.
The default is a linear-linear plot with dots as symbols. Log-linear and
log-log plots are also possible (this does not affect the printed output
described below.) Other symbols as described in "gm ?!" are also available.
Sp produces a four column file suitable for plotting with mongo or
other package. The four columns are x and y offsets from the image
center, the x image value, and the y image value. The points will all
be from within the specified area. The beginning of such a file
follows:
-11 -21 0.3192 2.6e-07
-10 -21 0.3849 2.571e-07
-9 -21 0.3784 2.483e-07
-8 -21 0.4007 2.424e-07
-7 -21 0.4447 2.306e-07
-6 -21 0.3818 2.189e-07
After making the graph, the integral of both images over the area is
shown on the standard output. The value shown is the sum of the points
times the area (in steradians) represented by each point in the map, e.g.:
Integral in x image is 0.000106336 [K-KM/S]*[SR] over 3.59618e-05 SR
Integral in y image is 7.77671e-12 [W/M2/SR]*[SR] over 3.59618e-05 SR
The integral is multiplied by the scale factor for that image.
Note for an integral with units [JY/BEAM]*[SR], you must divide
by the beam area in SR to get total JY. (Hint: use that number for
the scale factor.)
RELATED MACROS:
bmsr in .GMACROS will print out the beam area in steradians given
a FWHM in arcsec.
A NOTE ON PLOT SYMBOLS:
A list of figures and corresponding numbers is listed at the end
of "gm ?!". If the given figure width and height (fw:, fh:) are
zero (which is the initial value the first time you use "sp"),
then comb will scale them to the plot. Otherwise, the previous
value is used so that overlay plots (ovl:) come out looking nice.
So if you want comb to re-compute the plot symbol sizes, just
give "fh:0" and "fw:0".
This document was last updated on
Tue Oct 17 16:24:05 EDT 2000
by Marc W. Pound.