*** io ***
cvtByteOrder - equiv to swap_endian
OPENR_ - openr, /SWAP_IF_LITTLE_ENDIAN
OPENU_ - openu, /SWAP_IF_LITTLE_ENDIAN
OPENW_ - openw, /SWAP_IF_LITTLE_ENDIAN
SGK_PICKFILE - return a picked file
*** color ***
NOBPTT - do NOT bypass translation table
PC8 - Set device, pseudo_color = 8
PICKCOLOR - help pick a color, using HSV model (best w/ 24 bit disp)
PLCT - Load a `plotting' color map (16 predefined colors)
SGKLCT - Load a color table, using SGK's format
SGKLCT - Load a color table, using SGK's format
TC24 - Set device, true_color = 24
UDC - use decomposed color model for 24 bit dispays
XSGKLCT - GUI interface to SGKLCT
XSGKLCT - GUI to SGKLCT
*** math ***
CMP_AVG - computes average of argument
CMP_CC - compute cross-correlation
CMP_COVAR - compute covariance from Hessian matrix (use num. der.)
CMP_COVAR_A - compute covariance from Hessian matrix (analytical)
CMP_RMS - computes root-mean-squares of argument
CMP_WAVG - computes weighted average of arguments
DLSQ_FIT - perform non linear least-square fit
FCC - compute cross-correlation using FFT
LZ - Lorentzian((x-x0)/w)
FCN - c(0)*lz(x, c(1), c(2)) + c(3)
DLZ - 1st derivative wrt kth param of lz()
DFCN - 1st derivative wrt kth param of fcn()
D2FCN - 2nd derivative wrt kth & kkth param of fcn()
FRSTDER - computes the 1st derivative
FSMOOTH - frequency smoothing, over `n' pixels
FSMOOTH2D - 2d frequency smoothing, over `n' pixels
GEN_POLY - Generalized Polynomial Computation
getFWHM - computes FWHM of y(x)
GETMAX - compute the maximum, using parabolic fit
INTERP_1D - 1D linear interpolation
INTERP_2D - 2D bi-cubic interpolation
LGFIT - Fit a Legendre polynomial expansion
LGPOLY - evaluate Legendre polynomial expansion
LINFIT - perform a linear fit
LSQ - perform a Least-Squares fit
POLY2D - estimate a 2D polynomial
POLY_SOLVE - solve for x the polynomial expansion Sum_i Pi(x) = y
POLYFIT - perform a polynomial fit, include x-sigma rejection
POLYFIT2D - perform a 2D polynomial fit, include x-sigma rejection
POLYX - evaluates extended polynomial
QEQROOTS - solve quadratic equation a x^2 + b x + c
SANMBA - simulated annealing amoeba minimization
SCNDER - computes second derivative
SVDINVERT - invert matrix using single value decomposition
TZ_INTEGRATE - trapeze rule integration
TZ_INTEGRATE2D - 2D trapeze rule integration
XBIN - bin a data set xybin = xbin(x,y,n)
XSMOOTH - smoothing up to the edges
ZBIN2D - bin a data set zbin = zbin2d(x,y,z,nx,ny)
*** misc ***
@ADD-ALL-ASTRON-LIBS - add ALL pro/* for astron lib
addal - add ASTRON lib to the path
ADDSLASH - add a slash (/) to end of a string, unless ...
CTRL - return control char
cvt_d2HMS - return 'HH:MM:SS.SS' string of a value
cvt_HMS2D - return double value of a 'HH:MM:SS' string
cvt_jd2t - return 'YYYY.MM.DD_HH:MM:SS_UT' string for JD
cvt_s2hms - return 'HH:MM:SS.SSS' or '+D:HH:MM' string for s
cvt_s2t - return 'YYYY.MM.DD_HH:MM:SS' string for seconds
cvt_T2JD - convert 'YYYY.MM.DD_HH:MM:SS' string to JD
DEG2RAD - convert degree to radians
EPLAPSED - measured elapsed time since last call to ELAPSED
FILEEXIST - check if a file exists
FLIP - flip an array
FTOA - return string(f), using fixed formating rule
GETVAL - extract the value of a string in the form 'key: value'
IMAg - return imaginary part of value
ISAPEAK - return list of peaks
ISDOUBLE - return 1 if argument is of type DOUBLE
ISLONG - return 1 if argument is of type LONG
ISSTRING - return 1 if argument is of type STRING
ISSTRUC - return 1 if argument is of type STRUCTURE
ITOA - return strcat(string(round(i), format='(i0)'), sep)
NEEDSWAP - check OS byte order
NLAST - n_elements()-1
NLINES - return the no of lines in a file
@NO-WARNINGS - turn off all warnings
OLH - wrapper to widget_olh
PAUSE - pause until a CR is pressed
PCURSOR - print the current value of the cursor
PLOT_CLOCK - plot the time as a clock
PTELL - print a number to tell progress of a computation
RAD2DEG - convert radians to degrees
RAMP - return a ramp of N values for a given range
RCALL - return & close,/all
READ_RDB - read a RDB table, ie tab separated
READ_STRUCT - read a struct array from a flat ASCII table
READ_TABLE - read a flat ASCII table
REAL - return real part of value
REWIND, lu - point_lun, lu, 0
RHIST - [cumul] histogram w/ range
STRCAP - capitalize arg
STRCAT - catenate string array in one string w/ given separator
STRSEP - calls strsplit(string, separator, /extract)
WRITE_MGIFS - hacked version of write_gif, w/ LOOP keyword
WRITE_RDB - write a RDB table, ie tab separated
WRITE_STRUCT - write a struct array as a flat ASCII table
*** plot ***
AUTO_RANGE - return a padded plotting range
bMenu - button selection menu
BWIN - wrapper around window to open a big window
COLBAR - Add a color bar to a plot
DEMO_MODE - Returns true if IDL is in Demo Mode.
EB_PLOT - overplot error-bars
@EPS-PRINT - print the postscript plot
@EPS-PRINT-C - print the postscript plot
@EPS-PRINT-C - print the postscript plot
EPSCLOSE - call ps_close for EPS files (noprint, noid)
EPSOPEN - call ps_open for EPS files
FSDSP - Scaled display of a 2D array in a FITS file
GET_CPOS - get the cursor postion
GV - run cmd=ghostview on plot file
HHASH - plot hashes in the given yRange
hist_oplot - oplot histogram using polyfill
HLINE - overplot an horizontal line
IDENTIFY - identify points on a plot
LABEL - LaTeX alike string parser for labeling
MAKE_LEGEND - put a legend on a plot
MAKE_USERSYM - create the user symbol of given name
PLOT3D - 3d point plot
plot_vs_date - wrapper to plot, t, y, to labe x-axis as time/date spec
@PS-PRINT-C - print the postscript plot
@PS-PRINT-CP - print the postscript plot: color & portrait
@PS-PRINT-GV - print the postscript plot
@PS-PRINT-P - print the postscript plot, portrait mode
@PS-PRINT-Q - queue the postscript plot to be printed
@PS-PRINT - print the postscript plot
PS_CLOSE - Close the PostScript plotting device opened with PS_Open
PS_OPEN - Set the plotting device to be PostScript
PS_OPTS - Set options for plotting device if set to be PostScript
PUT_ID - Put an ID (time stamp, user, & hostname) on a plot
RELSCALE - Relative scaling, convert relative unit to data unit
RXYOUTS - "Relative" version of xyouts
SDSP - Scaled display of a 2D array (ie image)
SDSPOS - position equiv to sdsp, when using !p.multi
setRelPos - set relpos parameters for getrelpos()
SHOW3D - displays volumetric data using cuts
SHOW_FONT - Show all the characters in a given font
SPOSITION - Return dev. pos. for proper scaling of plots on images
SQUARE_PLOT - Restrict the plotting area to a relative square position
STV - Performs a "scaled" TV, see SWINDOW
STVSCL - Performs a "scaled" TVSCL, see SWINDOW
SWINDOW - Open a "scaled" window, w/ proper scaling whether X or PS
SXYOUTS - "scaled" version of XYOUTS, see SWINDOW
VEC2D - Plot 2D vector
VEC3D - Make a 3d vector plot
VECOL2D - Plot 2D vector, using color coding
VHASH - plot hashes in the given xRange
VLINE - overplot a vertical line
XBOPLOT - over-plot binned values
XFONTSEL - Select an X font
XGetScreenSize - return the screen(s) size (using xdpyinfo)
XPLOT_IO - Log plot x vs y, works when y is not always positive
XYLABEL - wrapper to xyouts, x, y, label(string), but add /rel
ZOOM - Zoom an array by factor(s) n1, [n2]
ZOOMTO - Zoom an array to be N1 x N2
*** widgets ***
getIDsFromCheckList - get widget IDs from checkList (see topMBW)
PlayAnyMovie - Play movie from list of any images.
PlayPPMMovie - Play movie from list of PPM images.
PlayPPMMovie - Play movie from list of PPM images.
setIDsFromCheckList - (re)set check marks from checkList (see topMBW)
SHOWPERCENT - Creates a widget to indicate fraction completed
TOPMBW - Creates a top level base widget w/ menu bar
XBLINK - blink two images
XWDEL - wdel for xwin
XWIN - widget wrapper around window to open a big window
XWSET - wset for xwin
; covar = cmp_covar_a(x, y, c, nfree)
;
; compute the covariance using the Hessian matrix, and analytical expressions
; of the derivatives, for the DLSQ problem y = fcn(x, c)
;
; keywords:
; weigths: weigths associated to (x,y) == 1/sigma
; (as for obsolete cmp_covar_aw)
; var: return variance
; hmat: return the Hessian matrix
; no2d: don't include 2nd deriv terms
; debug: if set, will print stuff, and stop inside routine
;
; calls dfcn(x, c, k) and d2fcn(x, c, k, kk), for 1st and 2nd derivatives
;
; if NFREE == 0 indicates maximum likelyhood, use
; NFREE = n_elements(c) otherwise
;
; <- Last updated: Fri Dec 4 16:51:48 1998 -> SGK
;
; NAME:
; CWX_ANIMATE
;
; PURPOSE:
; This widget displays an animated sequence of images using
; X-windows Pixmaps. This is a compound widget, based on the
; XINTERANIMATE procedure, with the following advantages:
; - It can be included in other applications.
; - Multiple copies can be run simultaneously.
;
; The speed and direction of the display can be adjusted using
; the widget interface.
;
; CATEGORY:
; Image display, compound widgets.
;
; CALLING SEQUENCE:
; To initially create:
; widget = CWX_ANIMATE(PARENT, SIZEX, SIZEY, NFRAMES)
;
; To reinitialize when another animation is loaded:
; CWX_ANIMATE_INIT, ANIMATEBASE, SIZEX, SIZEY, NFRAMES
;
; To load a single image:
; CWX_ANIMATE_LOAD, WIDGET, IMAGE = IMAGE, FRAME = FRAME_INDEX
;
; To load a single image that is already displayed in an existing window:
;
; CWX_ANIMATE_LOAD, WIDGET, FRAME = FRAME_INDEX, $
; WINDOW = [WINDOW_NUMBER [, X0, Y0, SX, SY]]
;
; (This technique is much faster than reading back from the window.)
;
; To display the animation after all the images have been loaded:
;
; CWX_ANIMATE, WIDGET [, RATE]
;
; To get a copy of the vector of Pixmaps being used by the widget.
; If this routine is called, the widget does not destroy the pixmaps
; when it is destroyed. The user can then provide them to a later
; call to CWX_ANIMATE to re-use them while skipping the Pixmap creation
; and rendering step:
;
; CWX_ANIMATE_GETP, widget, PIXMAPS
;
; INPUTS:
; CWX_ANIMATE:
; PARENT: The ID of the parent widget.
; SIZEX: The width of the displayed image.
; SIZEY: The height of the displayed image.
; NFRAMES: The number of frames in the animation sequence.
;
; CWX_ANIMATE_INIT:
; ANIMATEBASE: The ID of the base animation widget.
; SIZEX: The width of the displayed image.
; SIZEY: The height of the displayed image.
; NFRAMES: The number of frames in the animation sequence.
;
; CWX_ANIMATE_LOAD:
; WIDGET: The ID of the widget (as created with CWX_ANIMATE)
; into which the image should be loaded.
;
; CWX_ANIMATE_RUN:
; WIDGET: The ID of the widget (as created with CWX_ANIMATE)
; into which the image should be loaded.
; RATE: A value between 0 and 100 that represents the
; speed of the animation as a percentage of the
; maximum display rate. The fastest animation has
; a value of 100 and the slowest has a value of 0.
; The default animation rate is 100.
; STOP: If this keyword is set, the animation is stopped.
; NFRAMES: Specifies the number of frames to animate, must
; <= the number specified in CWX_ANIMATE().
;
; KEYWORD PARAMETERS:
; CWX_ANIMATE:
; PIXMAPS: This keyword provides the new widget with a vector
; of pre-existing pixmap (off screen window) IDs.
; This vector is usually obtained from a call to
; CWX_ANIMATE_GETP applied to a previous animation
; widget.
; UVALUE: A user supplied value to be stored in the widget's
; user value field.
; UNAME: A user supplied string name to be stored in the
; widget's user name field.
; NO_KILL: If NOT set, an "End Animation" button is added to the
; animation base. If set the button is not added.
; OPEN_FUNC: A user supplied string that specifies a callback
; function name. When a value is specified for this
; keyword, an "Open..." pushbutton is added to the
; window. When the "Open..." pushbutton is clicked
; the OPEN_FUNC function is called to load new
; animation data.
; INFO_FILE: A filename containing text to be displayed by
; XDISPLAYFILE when user selects the help button.
;
; CWX_ANIMATE_INIT:
; PIXMAPS: This keyword provides the new widget with a vector
; of pre-existing pixmap (off screen window) IDs.
; This vector is usually obtained from a call to
; CWX_ANIMATE_GETP applied to a previous animation
; widget.
;
; CWX_ANIMATE_LOAD:
; CYCLE: If set, cycle. Normally, frames are displayed
; going either forward or backwards. If CYCLE is
; set, reverse direction after the last frame in
; either direction is displayed.
; FRAME: The frame number to be loaded. This is a value
; between 0 and NFRAMES. If not supplied, frame 0
; is loaded.
; IMAGE: The image to be loaded.
; ORDER: Set this keyword to display images from the top
; down instead of the default bottom up. This keyword
; is only used when loading images with the IMAGE
; keyword.
; TRACK: If set, the frame slider tracks the current frame.
; Default is not to track.
; WINDOW: When this keyword is specified, an image is copied
; from an existing window to the animation pixmap.
; When using X windows, this technique is much faster
; than reading from the display and then loading with
; the IMAGE keyword.
;
; The value of this parameter is either an IDL window
; number (in which case the entire window is copied),
; or a vector containing the window index and the
; rectangular bounds of the area to be copied. For
; example:
; WINDOW = [Window_Number, X0, Y0, Sx, Sy]
;
; XOFFSET: The horizontal offset, in pixels from the left of
; the frame, of the image in the destination window.
;
; YOFFSET: The vertical offset, in pixels from the bottom of
; the frame, of the image in the destination window.
;
; OUTPUTS:
; No explicit outputs.
;
; SIDE EFFECTS:
; If the widget is realized before calls to CWX_ANIMATE_LOAD, the frames
; are displayed as they are loaded. This provides the user with an
; indication of how things are progressing.
;
; When the widget is destroyed, it destroys the pixmaps used in the
; animation, unless they were previously obtained via CWX_ANIMATE_GETP
; and the KILL_ANYWAY keyword was not set.
;
; The only event returned by this widget indicates that the user
; has pressed the DONE button. The parent application should use
; this as a signal to kill the animation widget via WIDGET_CONTROL.
;
; RESTRICTIONS:
; If more than one animation widget is running at a time, they
; will fight for resources and run slower.
;
; PROCEDURE:
; When initialized, this procedure creates pixmaps containing the
; frames of the animation sequence. Once the images are loaded,
; they are displayed by copying the images from the pixmap or buffer
; to the visible draw widget.
;
; EXAMPLE:
; Assume the following event handler procedure exists:
; PRO EHANDLER, EV
; WIDGET_CONTROL, /DESTROY, EV.TOP
; end
;
; Enter the following commands to open the file ABNORM.DAT (a series
; of images of a human heart) and load the images it contains into
; an array H:
;
; OPENR, 1, FILEPATH('abnorm.dat', SUBDIR = 'images')
; H = BYTARR(64, 64, 16)
; READU, 1, H
; CLOSE, 1
; H = REBIN(H, 128, 128, 16)
;
; Create an instance of the animation widget at load the frames:
;
; base = widget_base()
; animate = CWX_ANIMATE(base, 128, 128, 16)
; WIDGET_CONTROL, /REALIZE, base
; for i=0,15 do CWX_ANIMATE_LOAD, animate, FRAME=i, IMAGE=H[*,*,I]
;
; Start the animation:
;
; CWX_ANIMATE_RUN, animate
; XMANAGER, "CWX_ANIMATE Demo", base, EVENT_HANDLER = "EHANDLER"
;
; Pressing the DONE button kills the application.
;
; MODIFICATION HISTORY:
; AB, June 1992 Heavily based on the XINTERANIMATE procedure.
; SR, September 1992 Fixed a problem when a paused animation's
; frame selection was moved and the resulting
; frame change ended up in another animation.
; SR, November 1992 Fixed a problem when a single paused animation
; would fail when the frame selection slider
; event tried to set do a bad drawing window.
; DMS/AB, March, 1993 Got rid of state caching. Got rid of
; XMANAGER background tasks in favor of new
; "WIDGET_CONTROL,timer=" feature.
; ACY, October 1993 Set RETAIN=2 for draw widget to prevent
; clipping by an overlapping window when
; loading frames.
; DMS, Dec, 1993 Added STOP and NFRAMES keywords to CWX_ANIMATE_RUN.
; Added KILL_ANYWAY keyword to CWX_ANIMATE_GETP.
; WSO, Jan, 1995 Added OPEN_FUNC keyword and updated UI.
; ACY, Jan, 1997 Added INFO_FILE keyword to allow user-supplied
; files for help text
; JLP, Jan, 2000 Allow TrueColor images as input to CWX_ANIMATE_LOAD.
; Name:
; dlsq_fit
; Purpose:
; perform non linear least-square fit to
; Y = FCN(X, C)
;
; Calling Sequence:
; niter = dlsq_fit(x, y, c, eps, niterx)
; Keywords:
; silent: if set, perform fit silently
; debug: if set, print the 'A' and 'B' matrices
; Inputs:
; c(): first guess of fitting coefs
; eps(): relative accuracy required for each coef
; niterx: max no of iterations
; weights: weights associated to (x,y) == 1/sigma
; useNumDer: flag, if set, uses numerical derivative
; (instead of DFCN)
; delta(): array of relative changes to compute numerical
; derivatives (default: eps/2.)
;
; Output:
; niter: no if iteration required to achieve required accuracy
; (negative if error
; or did not converge in niterx iterations)
; c: fitting coefs
;
; dc: keyword w/ values of delta-c at last iteration
; hasConverged: keyword w/ index of coefs that have converged
;
; Notes:
; uses the function FCN(x,c) to fit,
; and DFCN(x,c,k) as derivative of FCN vs c(k)
;
; adjust only c( where eps gt 0 ) until abs(d-c/c) < eps,
; and # iter <= niterx
; Name:
;
; EB_PLOT
;
; Purpose:
;
; overplot error-bars (in y-direction by default)
;
; Calling sequence:
;
; EB_PLOT, X, Y, S
;
; Inputs:
;
; X: vector of absices
; Y: vector of ordinates
; S: error-bar size (y-s to y+s, or x-s to x+s if XAXIS is set)
;
; Input Keyword:
;
; XAXIS: if set, specify abcisses style error bars (x-axis)
; SERIF: add serifs
; LENGTH: relative length of the serifs, default : 0.005
; POSITIVE: draw only postive (y+s, or x+s) error bars
; NEGATIVE: draw only negative (y-s, or x-s) error bars
;
; Support all oplot keywords, as well
;
; Example:
;
; plot, x, y, psym=4
; eb_plot, x, y, s, color = 2
;
; <- Last updated: Tue Mar 6 19:49:36 2018 -> SGK
;
; Name:
; Gen_Poly
; Purpose:
; Generalized Polynomial Computation
; y(*) = Sum c(i).p(*, i)
; i = 0,#(c)-1
;
; Calling Sequence:
; y = Gen_Poly(c, p)
;
; Example:
; Computes Legendre polynomials to 5th order, for x = -1.,(0.1),1., then
; value for y = Sum_i c_i P_i(x)
;
; x = findgen(21)/10.-1.
; p =fltarr(21, 6)
; p(*,0) = 1.0
; p(*,1) = x
; for i = 2, 5 do $ ; Legendre polynomials recursion relation
; p(*,i) = ((2*i-1)*x*p(*, i-1) - (i-1)*p(*, i-2))/i
;
; c = [0.1, 420., 0.1, 20.0, 0.1, -5.0]
; y = gen_poly(c, p)
; plot, x, y
;
; See Also: LGPOLY LSQ
; index = identify(x, y, [/markit], [psym=], /loop)
;
; identify points on a plot, ie the closest value in x[], y[] selecte with the
; cursor; return the index/indices of the identified point(s)
;
; input keywords:
;
; markit = flag, of set will mark the identified point
; loop = flag, if set loop until click outside plot boundary
; psym = value of psym for marking identified points
;
; support all oplot keywords, like color, etc...
;
; example:
;
; plot, x, y, /psym
; idx = identify(x, y, /loop, /mark, color=2, symsize=2)
;
; lsq, y, sy, x, m, a, sa, r, rr, s, ss, $
; nonerror = nonerror, scaleSA = scaleSA
;
; least square solution to:
;
; y(k) +/- sy(k) = Sum a(i).x(k,i), for k=1,n
; i=1,m
;
; m< don't use that point
; r is the average residual, rr its stdev,
; s is the stdev of the prob, ss=std(s).
;
; keyword:
;
; nonerror: if set disable the 'on_error, 2'
; scaleSA: how to svcale SA
; if <0 sa = sqrt(cinv[i,i])
; if =0 sa = sqrt(cinv[i,i]) * sqrt(s)
; if >0 sa = sqrt(cinv[i,i]) * (sqrt(s) > scaleSA)
;
; Ref : Linnik, Y. V., Method of Least Squares and Principle of the Theory
; of Observation, Pergamon Press, Mc Millan Co., 1961, ch. 6, pp 135-187.
;
; <- Last updated: Wed Mar 7 13:20:57 2018 -> SGK
;
; Name:
;
; PLCT
; Purpose:
;
; Load a `plotting' color map:
; red, green, blue, yellow, magenta, indigo, etc...
; up to 14 colors.
;
; Calling Sequence:
;
; PLCT
;
; Keywords:
; n = - no of colors [3-14],
; default 14, min 3, max 14.
; offset = - offset at which to load the colmap,
; default 1, min 1.
; invertBW = flag, of set invert black and white levels
;
; Side Effects:
; Modify the color map, using TVLCT.
;
; Example:
;
; IDL> plct
; IDL> plot, findgen(3), thick=4, yr=[0, 6]
; IDL> for i = 1, 14 do oplot,findgen(3) + i/4., color = i, thick = 4
;
; Modification History:
; Jul 7 1995 Written by SGK, CfA.
; Jun 29 1998 modified to create a full color map, and deal w/ b/w
; levels, SGK
; May 13 2000 no action if TEK window system
; May 26 2002 call UDC if !D.N_COLORS NE !D.TABLE_SIZE
; PLOT3D, Z, X, Y
;
; makes a 3d point plot using a dummy
; SURFACE, /SAVE
; followed by a
; PLOTS, /T3D
;
; if [xyz]ranges are not set, uses SGK's AUTO_RANGE()
;
; The following keywords are supported:
; AX, AZ, XSTYLE, YSTYLE, ZSTYLE, XRANGE, YRANGE, ZRANGE,
; XTITLE, YTITLE, ZTITLE, TITLE, SUBTITLE,
; PSYM, SYMSIZE, CHARSIZE
; Name:
; c = polyfit(x, y, m)
; Purpose:
; perform a polynomial fit to
;
; y = Sum a P (x ) for i=0,m and j=0,N-1 w/ n >> m
; j i i i j
;
; by default P_i(x) = x^ii where ii = i-xZero
;
; Input Keywords:
;
; sy specify uncertainties for y,
; if sy(i) <= 0.0 this point is NOT used in the fit
; Legendre use legendre polynomials
; xZero offset to exponent def. 0 (none)
; noDC flag, if set does not include dc term (fit must go thru [0,0])
; usePX flag, if set uses the values passed in px= for P_i(x_j)
; nRej no of rejection passes, def. 0 (none)
; xRej x-sigma rejection factor, def. 3.0
; print flag, if set print the fitting coefs
; plot flag, if set plot
; pResid flag, if set plot plot residuals (sets plot)
; psym psym for oplot def. 1
; sSA see LSQ scaleSA (0, <0 or >0)
; chiFit rms of reduced residuals (used by scaleSA)
; noLabel flag, if set don't label values of coefs on plot
; cz character size for labels, def: 1.
; xtitle xtitle for the plot
; ytitle ytitle for the plot
; title subtitle for the plot
; noerase /noerase set to plot
; nonerror flag, if set disable the 'on_error, 2'
;
; Output:
; polynonial coefficients
;
; Output Keywords:
;
; rms RMS of plotResiduals to the fit
; cnt no of points included in the fit
; idi list of included indices
; idx list of excluded indices
; px values of P_i(x_j) (for gen_poly)
; sc uncertainties on coefs
; extra [r, rr, s, ss] as defined by LSQ, ie:
; r is the average residual, rr its stdev,
; s is the stdev of the prob, ss = std(s).
;
; <- Last updated: Thu Mar 8 08:17:03 2018 -> SGK
;
; Name:
;
; PS_Close
;
; Purpose:
;
; Closes the PostScript plotting device (opened with PS_Open),
; adds an optional time-stamp in the lower left corner,
; sends the PostScript file to the default printer
; and restore the plotting device saved by PS_Open.
;
; the default b/w printer is $PRINTER, or `ssplw',
; the default color printer is $COLORPRINTER, or `psc1',
;
; the default b/w print command is $LPCMD, or `lp',
; the default color print command is $COLORLPCMD, or `lp',
;
; Calling sequence:
;
; PS_Close
;
; Keywords:
;
; NOID - If set, the time-stamp ID is not added to the plot.
;
; NOPRINT - If set, the PostScript plotting device is closed, and
; the plotting device saved by PS_Open is restore,
; but the PostScript file is NOT sent to the printer.
;
; GV - set NOPRINT and run GV (ghostview previewer routine)
;
; SAVEAS = . Save the PostScript file by copying it to the
; specified file. The plot is still printed, unless NOPRINT is set.
;
; PRINTER = . Send the plot the specified printer.
; By default the enviroment variable PRINTER (or COLORPRINTER) is used to
; determine the printer to use. If it is not set the defaul values are
; ssplw (or psc1).
;
; TRANSPARENCY - If set, requests a transparency.
; Not all printers support this option (-SInputSlot=Transparency).
; (Only the color printers do)
;
; OPTIONS = . Optional extra options for the pslpr command.
; Example: OPTIONS = '-SDuplex' to use duplex printing.
;
; LPR = . Specify the command to use to print the file; by defaut
; pslpr is used. Note that if this keyword is specified, the keywords
; PRINTER, TRANSPARENCY and OPTIONS are ignored.
; Example: LPR = '/usr/bin/lp'
;
; SILENT - If set, disable warning messages
;
; QUEUE - If set, save the plot to print them all later
;
; NOEXEC - flag, if set won execute the spawn, cmd (for debug/testing mode)
;
; See also: PS_Open
; Name:
;
; PS_Open
;
; Purpose:
;
; Set the plotting device to be PostScript, and save the current plotting
; device (to restore it when the plot is done).
;
; Calling sequence:
;
; PS_Open
;
; Keywords:
;
; PORTRAIT - Normally, PS_Open defines the plotting area as landscape. If
; this keyword is set, the plotting area will be portrait.
;
; FILENAME = - Normally, PS_Open uses '/tmp/idl-$USER.ps' for the
; PostScript file, the keyword FILENAME can be used to override that value.
; ($USER is retrieved with a getenv('USER')).
;
; SIZES = <2-elements array> - Specify the plot size.
; The default plotting area is 10" x 7.5".
;
; OFFSETS = <2-elements array> - Specify the plot margins.
; The default plotting area has 0.5", 0.5" margins.
;
; SQUARE - if set, sets the plotting area to be square.
; Since the plotting area includes room for labels to axis, this option
; alone does not gives you a square plot (cf: !P.REGION).
;
; FLUSH - flush (ie ignore) all previously queued plots (see PS_Close)
;
; FONT - if numeric, set !P.FONT to value of FONT and turn on IsoLatin1
; - if string, set !P.FONT=0 and use the given PostScript font name as
; default (!3), checks the font name against a list, and if doesn't
; match use 'Helevetica-Narrow-Bold' as font
; (ie font=' ' equivalent to font='Helevetica-Narrow-Bold')
;
; PCZ - set !P.charsize to value PCZ (otherwise set it to zero)
; saves old value and ps_close restores it
;
; ENCAPSULATED - produces an encapsulated PostScript plot.
;
; BPP = - Set the number of bits per pixel; default value is 8,
; valid values are 8, 4, OR 2. A higher no produces a better
; resolution but also a bigger PostScript file.
;
; COLOR - If set, produces a color PostScript plot.
;
; TABLOID - If set, selects tabloid/ledger format (17" x 11" paper).
; The default plotting area is 16" x 10", with 0.5" margins.
; Not all printers can print tabloid format.
;
; LEDGER = TABLOID, /landscape
;
; PAGE_LENGTH = - Specify the actual page length.
; The default lenght is 17" for tabloid format, 11" otherwise.
; PAGE_WDITH = Specify the actual page width.
; The default lenght is 11" for tabloid format, 8.5" otherwise.
; This is needed because IDL's implementation of the X/YOFFSET keywords to
; DEVICE in landscape mode. The offset must be specified
; relative to the page width or length. This way PS_Open() offsets are more
; intuitive. (cf: IDL ref. guide, section: PostScript Positioning).
;
; NOFLIP - do not flip landscape mode (leave it as seascape)
;
; OKFontList - string array of valid font names to check font=
; against
;
; SILENT - If set, disable warning messages
;
; See also: PS_Close
; data = read_table(fn)
;
; read an ASCII file organized with a constant no. of col. per row, and return
; a float array. Uses the first line to deduce the number of columns (based
; on the column separator), unless nCols, format or struct are set.
;
; Input Keywords:
;
; nSkip: no of (header) lines to skip, def.: 0
; nCols: no of columns in the file, (don't try to guess using first line)
; colSep: column separator, def.: ' '
; format: format specification, def.: none
; [must supply nCols when using format, unless struct is given]
; [if format is supplied, colSep is ignored]
; autoSkip: flag, if set skip as many line as needed,
; skip lines that start w/ the comDel or are empty
; comDel: the comments delimitor def.: #
; double: flag, if set return double precision array
; verbose: flag, if set print out some information about the file
; struct: structure to use to read the data
; noOnerror: flag, if set won't do `on_error, 2'
;
; Output Keywords:
;
; header: header, ie the lines skipped
;
; <- Last updated: Mon Dec 4 11:12:24 2017 -> SGK
;
; sdsp, img, wno = wno, xrange = xr, yrange = yr, $
; xstyle = xstyle, ystyle = ystyle, zrange = zrange, $
; xsize = xsize, ysize = ysize, noScale = noScale, swFill = swFill, $
; xOffset = xOff, yOffset = yOff, noPmulti = noPmulti, $
; nlevels = nlevels, levels = levels, c_lines = c_lines, $
; zoom = zoomIn, sample = sample, $
; wtitle = wtitle, xtitle = xtitle, ytitle = ytitle, title = title, $
; margin1 = w1, margin2 = w2, tickf = tickf, $
; useSaveImgSz = useSaveImgSz, useIndivImgSz = useIndivImgSz,$
; absOffset = absOffset, noXC = noXC, _extra = _extra
;
; scaled display of a 2D array (ie image)
;
; zoomIn can be either a scalar or a 2-elements array
; if nlevels > 0, coutours are superposed
; if wno < 0, use current window
;
; uses