This page was created by IDL
lasco_mk_html_help.pro on
Wed Aug 17 12:22:02 2005.
NAME:
CARRDATE
PURPOSE:
This function converts a carrington number and longitude into a
CDS date structure.
CATEGORY:
LASCO Synoptic
CALLING SEQUENCE:
Result = CARRDATE ( Cr, Clong)
INPUTS:
Cr: Carrington Rotation Number
Clong: Carrington Longitude
OPTIONAL INPUTS:
None
KEYWORDS:
EL: This keyword specifies that the date should be returned as an
East Limb Date
WL: This keyword specifies that the date should be returned as a
West Limb Date
CMP: This keyword specifies that the date should be returned as a
central meridian Date. This is the default.
OUTPUTS:
Result: The function result is the date as a CDS date structure
COMMON BLOCKS:
carr_long: Contains the start date of the carrington rotations
from the almanac. Generated by READ_CARR_LONG
MODIFICATION HISTORY:
Written by: RA Howard, NRL, 3/18/96
10/17/96 RAH, Added check for cday being defined
, Changed output to CDS structure
@(#)carrdate.pro 1.3 10/17/96 LASCO IDL LIBRARY
NAME:
CARRDATE
PURPOSE:
This function converts a carrington number and longitude into a
date.
CATEGORY:
LASCO Synoptic
CALLING SEQUENCE:
Result = CARRDATE ( Cr )
INPUTS:
Cr: Carrington Rotation Number
OPTIONAL INPUTS:
None
KEYWORDS:
EL: This keyword specifies that the date should be returned as an
East Limb Date
WL: This keyword specifies that the date should be returned as a
West Limb Date
CMP: This keyword specifies that the date should be returned as a
central meridian Date. This is the default.
Clong: Carrington Longitude
OUTPUTS:
Result: The function result is returned as a fractional modified
Julian date.
COMMON BLOCKS:
carr_long: Contains the start date of the carrington rotations
from the almanac. Generated by READ_CARR_LONG
MODIFICATION HISTORY:
Written by: RA Howard, NRL, 3/18/96
11/20/96 by N. RIch added read_carr_long call if cday=0; changed
cday(n-1) to cday(n) and cday(n) to cday(n+1);
made dte variable same as day; initialized clong
to 360 if not set and change clong to keyword;
changed +/- in computing el and wl day
12/30/96 by N. Rich changed output to CDS structure
@(#)carrdate.pro 1.1 10/05/96 LASCO IDL LIBRARY
NAME: CARRLONG PURPOSE: This procedure computes the carrington number and longitude given the date. CATEGORY: LASCO Synoptic CALLING SEQUENCE: CARRLONG, Date, Cr, Clong INPUTS: Date: The date in CDS date/time structure OPTIONAL INPUTS: None KEYWORD PARAMETERS: EL: Set this keyword to obtain information for East Limb date WL: Set this keyword to obtain information for West Limb date CMP: Set this keyword to obtain information for central meridian This is the default. OUTPUTS: Cr: The carrington rotation number corresponding to the input date Clong: The carrington longitude corresponding to the input date COMMON BLOCKS: carr_long: Contains the start date of the carrington rotations from the almanac. Generated by READ_CARR_LONG EXAMPLE: CARRLONG, date, cnum, clong, /EL MODIFICATION HISTORY: Written by: RA Howard, NRL, 3/18/96 10/17/95 RAH, Converted dates to CDS time structure 13/10/00 nbr, Changed "tai LT cday" to "tai LE cday" @(#)carrlong.pro 1.4 11/05/01 LASCO IDL LIBRARY
ro carrmapdisp, limb, r, RANGE=range, VERTICAL=vertical, CMAPDIR=cmapdir Purpose: Display carrington maps in a scrollable format (NOTE: previous carrmapdisp.pro renamed carrmapdisp0.pro) Optional Inputs: limb STRING 'w', 'e' or search parameter r FLOAT height of desired map KEYWORDS: RANGE = [cn0, cn1] INTARR range of carrington numbers to display /VERTICAL If set, display maps stacked, not side-by-side CMAP='/path/of/fits/files' defaults to location at NRL Written 980917 by Nathan Rich, Interferometrics/NRL 030821, nbr - add RANGE, vertical display options, rename from slide_carrmap.pro 030829, nbr - fix bug in list=findfile() @(#)carrmapdisp.pro 1.2 09/02/03 - LASCO NRL IDL Library
Procedure: CARRMAPMAKER2
PURPOSE:
Returns carrington maps at multiple radii for a passed
carrington number at radii read from user in this routine.
Program only works for maps beginning after
January,6 1996.
INPUTS:
cn INT: carrington number
cam STRING: camera, i.e. 'c2'
num_r INT: number of radii to make maps for
rad FLTARR(num_r): radii for maps
cmap FLTARR(mapsize,ht,num_r): passed empty
pol STRING: 'p','u' or ''; indicates polarizer value
limb INTARR: 0 or 1 tells which limb for carrdate procedure
level STRING: 'ql' for quick-look, 'r' for doing ratio
disp INT: 1 to display images
wlimb STRING: 'wl' or 'el'
hdr STRUCT: dummy header sent back to carrmap3 for fits header
KEYWORDS:
FULL360 If set, go all the way around, not just 180 deg.
MAP_SIZE IF set, then use it's value for horizontal size
This routine can ask the user for each radius.
OUTPUTS: rad = fltarr(num_r) : radii of maps
cmap = fltarr([360 or 720],ht,num_r): carrington maps, one per radius
ROUTINES CALLED:
qdb.pro
carrdate2.pro
resize2.pro
getc2c3norm2.pro
AUTHOR: Nathan Rich, NRL, Nov. 1996
Julia Kraemer, NRL, June 7, 1996
MODIFIED:
961129 by N. RIch each image subtends one pixel; checks for other
values of val0 if 0 when k=0
961203 by N. RIch if more than one pixel between images, use
congrid to fill space
961206 by N. Rich Create log file; rotate images from 960521;
Check median = 85 (for 1024x576 images)
961209 by N. Rich change CONGRID statement; moved carrmap_x
increment statement to after validity check
961210 by N. Rich print time0 and time1; if mapdelta GT 10, don't
do congrid
961211 by N. Rich skip bad files
961213 by N. Rich add level variable; if mapdelta GT 20, don't do
congrid; check for forward time
961217 by N. Rich add disp variable for displaying image
970102 by N. Rich make time units consistent
970103 by N. Rich rotate images for Nov. 21,22
970103 by N. Rich add white strip to top of gaps
970108 by J. Kraemer add east limb option
970123 by N. Rich changed to val0(r)=median(strip(115:135))
970213 by N. Rich changed check for Out of Order
970226 by N. Rich do normalization for whole image rather than
each strip
970311 by N. Rich make map display scale vary with camera
970402 by N. RIch add empty space to final columns if no images
970414 by N. Rich accept c1 images for Clarence
9705 by N. Rich use c1 images generally
970714 by N. Rich change criteria for data gaps
970715 by N. Rich fix final gap coverer
970722 by N. Rich check each strip against previous and reject
image if it is too different
970730 by N. Rich sort s by date_obs
970801 by N. Rich skip image if mapdelta LE 0
970807 by N. Rich scan strip ends; move imdisplay part around
970827 by N. Rich change scale; modify name of log files
970915 by N. Rich add get_im common block, compute ff_ratio; change
cx,cy for c1
971021 by N. Rich compute rconv for each image separately for c2 and c3
971103 by N. Rich prompt for mapsize
971107 by N. Rich move log print line; prompt for median skip; don't
check image median
971110 by N. Rich add skiplast variable
971114 by N. Rich use strip files if available
971216 by N. Rich make automatic (auto-enter radii, mapsize)
971217 by N. Rich if error in readfits, enter in log
980212 by N. Rich added variable boxes for normalization
980312 by N. Rich cancel box normalization; auto choose mapsize
980313 by N. Rich do not multiply image by constant (of about 0.12)
980413 by N. Rich use m = median(im(*,y21:y12))
980514 by N. Rich change upper limit of strip median, ends
980610 by N. Rich extend warpit COMMON block for fits header
980612 by N. Rich pass rad (radii) already full
980629 by N. Rich use reduce_std_size instead of resize2
981216 by N. Rich Use ht variable instead of 181 for mapheight
990114 by N. Rich Add FULL360 keyword
Feb 2000 by N. RIch Remove source from query; change GET_BKG init.;
Base level on LZ keyword, levelstr; Replace strip with
strp to avoid IDL conflict
Apr 2000, NBR, Add half a column's time to t_start; always use first image (k=0)
May 2000, NBR, Use mk_img instead of getc123; don't use avg_strip_ends
Aug 2000, NBR, Add MAP_SIZE keyword and take it out of common block; compare strip,stripends
2001.11.05, NBR - Add SCCS version tag
2002.03.14, NBR - Add skipwait, medianhist, variable sunc from mk_img, roll in mk_img
2002.09.03, nbr - Expand logging, add l1 to COMMON get_im
09/08/03 @(#)carrmapmaker2.pro 1.6 - IDL NRL LASCO Library
unction cmap2gif, ftsfile, maxd,mind, ROOT=root, CONTROL=control Converts a FITS carrington map into a byte array with scaled image and labeled axes. See front end called make_cmap_gifs.pro. INPUTS ftsfile STR The name of the carrington map FITS file. OPTIONAL INPUTS: maxd, mind FLOAT Scale min and max KEYWORDS CONTROL If set, interactively choose max/min for bytscl ROOT Returns root portion of FITS filename to named variable MODIFIED 9. 7.01, nbr - Change automax for 195A 03.05.15, nbr - Fix sf for c2/c3 05/15/03 @(#)cmap2gif.pro 1.2 : IDL LASCO NRL Library
TITLE:
COMB_FULL_EQ
PURPOSE:
This procedure combines the monthly equatorial region image
and the monthly full field image into a single monthly model.
CATEGORY:
DATA_ANAL
CALLING SEQUENCE:
COMB_FULL_EQ,Tel,Td
INPUTS:
Tel: String denoting the telescope: 'c1','c2','c3','c4'
or 'eit'
Td: String denoting the date, in one of the following
formats: YYYY/MM/DD
YYY-MM-DD
YYMMDD
KEYWORDS:
FILTER: If present this string specifies the filter to be used.
The defaults are the following:
C1 = 'FeXIV'
C2 = 'Orange'
C3 = 'Clear'
C4 = 'Clear'
POLAR: If present this string specifies the polarizer (or sector)
to be used. The defaults are the following:
C1 = 'Clear'
C2 = 'Clear'
C3 = 'Clear'
C4 = '304'
OUTPUTS:
None. The routine will store the minimum image in a file
named tmdddddd.fts, where t is the telescope number,
m is the letter 'm', and dddddd is the 6 digit date
in the format YYMMDD.
PROCEDURE:
The routine looks for file names of the form, tmxxxYYMMDDfilpol.fts
where t = '1','2','3','4' for the telescope
m = 'm' for monthly
xxx = 'equ' or 'ful' for equatorial or full images
YYMMDD = The year, month, day
fil = the filter code:
C1: 'fv', 'fx', 'na', 'ca', 'cl'
C2: 'or', 'rd', 'bl', 'ha'
C3: 'cl', 'or', 'rd', 'bl', 'dr'
C4/EIT: 'cl', 'a1', 'al4'
pol = the polarizer/sector code:
C1: 'cl', 'p0', 'pm', 'pp', 'ha'
C2: 'cl', 'p0', 'pm', 'pp'
C3: 'cl', 'p0', 'pm', 'pp', 'ha'
C4/EIT: '304', '284', '195','171'
WRITTEN:
RA Howard, NRL, 9/27/96
@(#)comb_full_eq.pro 1.1 10/21/96 LASCO IDL LIBRARY
NAME:
DIFBKGND
PURPOSE:
Difference the appropriate background image
CATEGORY:
DATA_ANAL
CALLING SEQUENCE:
Result = DIFBKGND(Img,Hdr)
INPUTS:
Img = Input Image array. It is assumed that the array is full image,
(either 1024x1024 or 512x512) so that if the original image was
a subimage, it must be put into the proper place in a full image
array.
Hdr = FITS header
KEYWORDS:
MINVAL: If set, the minimum value for the byte scaling is set to
this value. The default is the minimum value in the image.
MAXVAL: If set, the maximum value for the byte scaling is set to
this value. The default is the maximum value in the image.
LOGSCL: If set, a logarithmic scaling is applied. The default is
linear scaling.
ALTMIN: If set, the minimum value for the byte scaling when the
background image is NOT found is set to this value. The
default is the minimum value in the image.
ALTMAX: If set, the maximum value for the byte scaling when the
background image is NOT found is set to this value. The
default is the maximum value in the image.
ALTLOG: If set, a logarithmic scaling is applied when the background
image is NOT found. The default is linear scaling.
RATIO: If set the input image is divided by the model. Zeroes in
the model or input are set to 1. The default is to perform
a subtraction.
NO_SCALE:If set the output image is not byte scaled. The default is
to byte scale the image.
OUTPUTS:
The function returns a byte scaled difference image.
PROCEDURE:
The input header is examined to extract the telescope, filter and
polarizer/sector. The directory pointed to by the environment
variable $MONTHLY_IMAGES is searched for all monthly images that
exist for that combination. If only one monthly image is found it
is used. Otherwise the monthly image whose Mid-Date is closest to
the date of the input image is used.
No background subtraction is performed if the input image is a dark,
cal lamp, and continuous mode. In these cases the image is byte
scaled using histogram equalization. If the image is not 1024x1024.
or 512x512 then no background subtraction is used, but the image is
scaled according to the input keywords.
MODIFICATION HISTORY:
Written, RA Howard, NRL, 16 October 1996
24 Oct 96, RAH, Added ratio
25 Oct 96, RAH, Added save of last background in common
16 Jul 97, RAH, Modified call to OFFSET_BIAS to pass header
26 Feb 98, RAH, Modified method for determining whether a background image was found
29 Sep 99, NBR, Remove CURRENT keyword from GETBKGIMG call
@(#)difbkgnd.pro 1.6 09/29/99 LASCO IDL LIBRARY
NAME: CONGRID PURPOSE: Simulate the action of the VAX/VMS CONGRID/CONGRIDI function. Shrink or expand the size of an image. CATEGORY: Image processing. CALLING SEQUENCE: Result = CONGRID(Image, Xs, Ys [, Interp = Interp]) INPUTS: Image = 2D array to resample. Xs = desired number of columns for result. Ys = number of rows for result. KEYWORD PARAMETERS: Interp = keyword which if set causes bilinear interpolation to be used. Otherwise nearest neighbor method is used. OUTPUTS: Result = Image of same type as input, of size (Xs, Ys). COMMON BLOCKS: None. SIDE EFFECTS: None. RESTRICTIONS: Doesn't completely emulate the VAX/VMS CONGRID. The case of a rectangular grid is not implemented. This can be done using multiple calls to POLY_2D. PROCEDURE: Simple call POLY_2D with the warping coefficients. MODIFICATION HISTORY: DMS, Sept. 1988.
********************************************************************************
LASER8ew.PRO
LASER8 for printing east limb and west limb carrington maps in landscape mode
designed to be significantly more modular than its predecessors
and hopefully more comprehensible to the peruser (and myself).
Eric T Swanson
8/27/91
PROCEDURES CALLED:
makebasecarr.pro
MODIFIED:
Dec. 1996 by N. Rich make output in landscape mode
Jan. 1996 by J. Kraemer add difference and limb option
Jan. 1996 by N. Rich add east/west limb feature; do actual differencing;
display begin and end dates; display max and min (modified common
block PM1)
Feb. 14 96 by N. Rich touch up ew maps
Feb. 26 96 by N. Rich remove factor by which basecarr is multiplied
Mar. 13 97 by N. Rich add B0 angle
Mar. 14 97 by N. Rich add top sub-title; newmax=newmax-1 for row=1
Mar. 18 97 by N. Rich change B0 angle display
Mar. 28 97 by N. Rich use pb0r.pro to compute B0 angle
May 15 97 by N. Rich variable mapsize from filename
May 21 97 by N. Rich add my name to output
Nov 3 97 by N. Rich use ratio instead of difference
Apr 7 98 by N. Rich add today's date to output
Jun 26 98 by N. Rich read fits files
Mar 2000 By N. Rich Add kind to common PM1; edit PRINTMENU
Oct 2000 By N. Rich Use BYTSCL for scaling images
Nov 2001 by N. RIch Add SCCS tags
27Mar02 by N. RIch Change default Top2 title
1Apr02 by N. Rich Auto exit after save
*************************************************************************
%H% %W% - NRL LASCO IDL Library
NAME: CONGRID PURPOSE: Simulate the action of the VAX/VMS CONGRID/CONGRIDI function. Shrink or expand the size of an image. CATEGORY: Image processing. CALLING SEQUENCE: Result = CONGRID(Image, Xs, Ys [, Interp = Interp]) INPUTS: Image = 2D array to resample. Xs = desired number of columns for result. Ys = number of rows for result. KEYWORD PARAMETERS: Interp = keyword which if set causes bilinear interpolation to be used. Otherwise nearest neighbor method is used. OUTPUTS: Result = Image of same type as input, of size (Xs, Ys). COMMON BLOCKS: None. SIDE EFFECTS: None. RESTRICTIONS: Doesn't completely emulate the VAX/VMS CONGRID. The case of a rectangular grid is not implemented. This can be done using multiple calls to POLY_2D. PROCEDURE: Simple call POLY_2D with the warping coefficients. MODIFICATION HISTORY: DMS, Sept. 1988.
NAME: MAKE_ALL_DAYS PURPOSE: This procedure is an easy interface to the MK_DAILY_MED procedure CATEGORY: LASCO SYNOPTIC CALLING SEQUENCE: MAKE_ALL_DAYS,Tel,Stdte,Endte INPUTS: Tel: Telescope to be processed, (string), eg. 'C1','C2','C3' Stdte: The starting date specified as a string: YYMMDD, eg. 960601 Endte: The ending date specified as a string: YYMMDD, eg, 960630 OUTPUTS: This function generates images in $MONTHLY_IMAGES of the daily medians for the specfied telescope, and for each day specfied by the starting and ending dates. KEYWORDS: QL: If set then use quick look date, else use level 0 data FILTER: String denoting the filter position. The default depends on the telescope: C1: FeXIV C2: Orange C3: Clear C4: Clear POLAR: String denoting the polarizer position. The default depends on the telescope: C1: Clear C2: Clear C3: Clear C4: 304A WLU: String denoting the upper FP wavelength of an interval. ANYSIZE: Accept images with any dimensions GT 256 NOREBIN: Do not rebin images to 512x512 ROLLED: Save result in 'daily/rolled' directory; set equal to roll angle (CROTA) VERBOSE: Print information on rows during loop RESTRICTIONS: PROCEDURE: A loop is generated between the start and end dates, in which the routine MK_DAILY_MED is called for each day. EXAMPLE: To make all the daily median images for the month of June, 1996 for C3: MAKE_ALL_DAYS,'c3','960601','960630' MODIFICATION HISTORY: Written by: RA Howard, 20 October, 1996 9/20/99 NBR Had to add QL to mk_daily_med call 12/1999 NBR Make y2k compliant; add MK_DAILY_MED keywords 07/2000 NBR Add ROLLED keyword 5/2/2002, NBR Add keyword descriptions 07/14/03 @(#)make_all_days.pro 1.5 LASCO IDL LIBRARY
NAME: MAKE_ALL_MINDAYS PURPOSE: This procedure is an easy interface to the MK_DAILY_MIN procedure CATEGORY: LASCO SYNOPTIC CALLING SEQUENCE: MAKE_ALL_MINDAYS,Tel,Stdte,Endte INPUTS: Tel: Telescope to be processed, (string), eg. 'C1','C2','C3' Stdte: The starting date specified as a string: YYMMDD, eg. 960601 Endte: The ending date specified as a string: YYMMDD, eg, 960630 OUTPUTS: This function generates images in $MONTHLY_IMAGES of the daily minimums for the specfied telescope, and for each day specfied by the starting and ending dates. RESTRICTIONS: These routines only generate the default filter and polarizer images. PROCEDURE: A loop is generated between the start and end dates, in which the routine MK_DAILY_MIN is called for each day. EXAMPLE: To make all the daily median images for the month of June, 1996 for C3: MAKE_ALL_MINDAYS,'c3','960601','960630' MODIFICATION HISTORY: Written by: RA Howard, 20 October, 1996 9/20/99 NBR Had to add QL to mk_daily_med call 12/1999 NBR Make y2k compliant 11/05/01 @(#)make_all_mindays.pro 1.1 LASCO IDL LIBRARY
NAME: MAKE_ALL_MONEXP PURPOSE: This procedure is an easy interface to the MONITOR_EXP procedure CATEGORY: LASCO SYNOPTIC CALLING SEQUENCE: MAKE_ALL_MONEXP,Stdte,Endte INPUTS: Stdte: The starting date specified as a string: YYMMDD, eg. 960601 Endte: The ending date specified as a string: YYMMDD, eg, 960630 Lz: 0 for QL, 1 for LZ OUTPUTS: This function generates images in $MON_EXP of the statistics for each image for each day specfied by the starting and ending dates for all three LASCO telescopes. PROCEDURE: A loop is generated between the start and end dates, in which the routine MONITOR_EXP is called for each day. EXAMPLE: To compute all of the image statistics for the month of June, 1996: MAKE_ALL_MONEXP,'960601','960630' MODIFICATION HISTORY: Written by: RA Howard, 22 May 1997 @(#)make_all_monexp.pro 1.1 09/26/97 LASCO IDL LIBRARY
NAME:
MAKE_ALL_MONTHS
PURPOSE:
This procedure is an easy interface to the MK_MONTHLY_MIN procedure
CATEGORY:
LASCO SYNOPTIC
CALLING SEQUENCE:
MAKE_ALL_MONTHS,Tel,Stdte,Endte
INPUTS:
Tel: Telescope to be processed, (string), eg. 'C1','C2','C3'
Stdte: The starting date specified as a string: YYMMDD, eg. 960601
Endte: The ending date specified as a string: YYMMDD, eg, 960630
OUTPUTS:
"Monthly" minimum FITS files in $MONTHLY_IMAGES; filenames have dates
corresponding to 3m_clcl_xx'*'.fts
Logs written in $MONTHLY_IMAGES/logs
KEYWORDS:
ALL_YEARS: Use daily images from multiple years.
FWPW: String representing filter/polarizer combination; if not
set, defaults to 'clcl' for C3 and 'orcl' for C2
INTERVAL: Number of days between minimum images; default is 7
NDAYS: Number of days used to compute the min; default is 27
MANUAL: Override automatic determination of days to use and filename
TEST: Write output in $MONTHLY_IMAGES/testdir
DOREBIN: Make output half size (2x2 binned)
RESTRICTIONS:
PROCEDURE:
IF MANUAL is not set, this procedure calls mk_monthly_min.pro for dates
corresponding to the included list of dates in between the starting and
ending dates. Gaps of more that 2 days are not spanned. IF MANUAL is set,
all days in the input interval (inclusive) which exist will be used to
compute the minimum image, and it will be named according to the midpoint
of the interval of days used.
EXAMPLE:
To make all the monthly minimum images for the month of June, 1996 for
C3:
MAKE_ALL_MONTHS,'c3','960601','960630'
MODIFICATION HISTORY:
Written by: RA Howard, 20 October, 1996
Mofified:
06 Oct 1997 SEP - Added ALL_YEARS flag
1999/02/11 NBR - Added keywords FILT_POLR, DAYS_BETWN, DAYS_COMPT
1999/12 NBR - Make y2k compliant
Jan 2000 NBR - Add NOREBIN keyword
Jun 2000 NBR - Add ROLL keyword
25Feb02, NBR - Use 3m_clcl_xx'*'.fts for dates used.
10Jul02, NBR - Add TEST keyword
17Jul02, NBR - Change NOREBIN to DOREBIN
28Jan04, NBR - default nbetween = 7
%W% %H% LASCO IDL LIBRARY
NAME:
MAKE_DAILY_IMAGE
PURPOSE:
This function generates the daily image for the SOHO summary images.
CATEGORY:
LASCO_SYNOPTIC
CALLING SEQUENCE:
Result - MAKE_DAILY_IMAGE (Date, Tele)
INPUTS:
Date: Date for which the image is desired, (YYMMDD)
If undefined, then uses the current date
Tele: Telescope designator, string of 'c2','c3'
KEYWORD PARAMETERS:
TIME: If set use the image closest to the specified time (HHMM)
else uses the latest image
FITS: If set create a FITS file as output
GIF: If set create a GIF file as output
LZ: IF set, use level 0 data else use quick look data
SCREEN: If set, writes the resultant image to the current display
else writes to the Z buffer
ARCHIVE: If set, write to the SOHO Summary file area
OUTPUTS:
Result: Returns the byte-scaled image if successful, else returns -1.
RESTRICTION:
Only works for C2 and C3
PROCEDURE:
If the input parameter, date, is defined, then that date is used
otherwise the current date is used. Then reads in the img_hdr.txt
file corresponding to the desired date.
Finds the image that is closest to the input time if set, else
finds the latest image. The standard synoptic type of image is
selected (C2: orange/clear; C3: clear/clear). The image shold
be large enough and not have too many missing blocks.
Then the image is scaled by the background model.
EXAMPLE:
cimg = MAKE_DAILY_IMAGE('960810','c3')
MODIFICATION HISTORY:
Written by: S.E. Paswaters, NRL, Dec 1997
1/31/01, nbr - Change bkg model settings for C2
%H% %W% LASCO IDL LIBRARY
TITLE: MK_ALL_MIN PURPOSE: This procedure generates the minimum of the all of the monthly minimum images. CATEGORY: DATA_ANAL CALLING SEQUENCE: MK_ALL_MIN,Tel,Filpol INPUTS: Tel: String denoting the telescope: 'c1','c2','c3','c4' Filpol: String denoting the filter/polarizer configuration OUTPUTS: None. The routine will store the minimum image in a file named tm_all.fts, where t is the telescope number, m is the letter 'm', and _all is the string '_all' MODIFICATION HISTORY: WRITTEN RA Howard, NRL, 3 Oct 96 8 Nov 96 RAH, checks for bad blocks are <= 0 30 Jan 97 SEP, check filenames for YYMMDD @(#)mk_all_min.pro 1.4 11/05/01 LASCO IDL LIBRARY
NAME:
MK_DAILY_C1_MED
PURPOSE:
This procedure generates an image by taking all the files of a given
type for one day and finding the median value for each pixel.
CATEGORY:
DATA_ANAL
CALLING SEQUENCE:
MK_DAILY_C1_MED,Tel,Date
INPUTS:
Tel: String denoting the telescope, 'c1','c2','c3','c4'
Date: Gives the date to be processed in one of the following formats:
YYMMDD, 6 character string
YYYY/MM/DD, 10 character string
CDS Date Structure
Long Word of the modified julian date
OPTIONAL INPUTS:
None
KEYWORD PARAMETERS:
FILTER: String denoting the filter position. The default depends on
the telescope:
C1: FeXIV
C2: Orange
C3: Clear
C4: Clear
POLAR: String denoting the polarizer position. The default depends on
the telescope:
C1: Clear
C2: Clear
C3: Clear
C4: 304A
WLL: String denoting the lower FP wavelength of an interval.
WLU: String denoting the upper FP wavelength of an interval.
ONOFF: String denoting whether the on line or off line wavelength should be used
QL: If set uses QL data, else LZ
OUTPUTS:
None.
SIDE EFFECTS:
Writes a fits file to $MONTHLY_IMAGES/daily.
RESTRICTIONS:
Only looks for full resolution, full width, any height.
PROCEDURE:
For each image that satifies the selection conditions, (naxis1=1024),
filter and polarizer as requested), the median image is computed of
the median value of all the images for a single day after being
normalized to the median exposure time.
If the number of images is less than 7 for any pixel, that pixel is
not computed in the first pass. In the second pass, all full images
within +/- 2 days of the given day.
EXAMPLE:
To create the daily median image for 1 Sep 1996:
MK_DAILY_C1_MED,'c3','960901'
or
MK_DAILY_C1_MED,'c3','1996/09/01'
MODIFICATION HISTORY:
Written by: R.A. Howard, NRL, 9/27/96
Version 2, RAH, 10/10/96, Modified to accept all image sizes
Version 3, RAH, 10/21/96, Added filter and polarizer as inputs
Version 4, SEP, 11/12/96, Fixed problem with equatorial img boundary lines.
Added check for ADCT compressed images.
NBR, 8/ 9/01, Change output directory to $MONTHLY_IMAGES/daily
@(#)mk_daily_c1_med.pro 1.2 08/09/01 LASCO IDL LIBRARY
NAME:
MK_DAILY_MED
PURPOSE:
This procedure generates an image by taking all the files of a given
type (up to 25) for one day and finding the median value for each pixel.
CATEGORY:
DATA_ANAL
CALLING SEQUENCE:
MK_DAILY_MED,Tel,Date
INPUTS:
Tel: String denoting the telescope, 'c1','c2','c3','c4'
Date: Gives the date to be processed in one of the following formats:
YYMMDD, 6 character string
YYYY/MM/DD, 10 character string
CDS Date Structure
Long Word of the modified julian date
OPTIONAL INPUTS:
None
KEYWORD PARAMETERS:
QL: If set then use quick look date, else use level 0 data
FILTER: String denoting the filter position. The default depends on
the telescope:
C1: FeXIV
C2: Orange
C3: Clear
C4: Clear
POLAR: String denoting the polarizer position. The default depends on
the telescope:
C1: Clear
C2: Clear
C3: Clear
C4: 304A
WLL: String denoting the lower FP wavelength of an interval.
WLU: String denoting the upper FP wavelength of an interval.
ONOFF: String denoting whether the on line or off line wavelength should be used
ANYSIZE: Accept images with dimensions GT 256
NOREBIN: Do not rebin images to 512x512
SAVEDIR: Specify where to save result
ROLLED: Save result in 'rolled' directory; set equal to roll angle (CROTA)
FILES: If set to a strarr, then use files in array to compute the median image
VERBOSE: Print row information for zero rows
OUTPUTS:
None.
SIDE EFFECTS:
Writes a fits file to $MONTHLY_IMAGES/daily of the format td_fwpw_yymmdd[raaa[x]].fts
where t = telescope [1,2,3]
d = d (daily image)
fw = filter abbreviation from abbrv_filpol.pro
pw = polarizer abbreviation from abbrv_filpol.pro
yymmdd = date of the image
r = r if roll of image GT 2 degrees from north
aaa = average (dwell) roll angle during interval in degrees CCW
x = x if roll angle varies more than 5 degrees during interval
RESTRICTIONS:
PROCEDURE:
For each image that satifies the selection conditions, (default naxis1=1024,
filter and polarizer as requested), the median image is computed of
the median value of all the images for a single day after being
normalized to the median exposure time.
If the number of images is less than 7, there is a second pass.
In the second pass, images within +/- 2 days of the given day are used,
up to 15 per day.
EXAMPLE:
To create the daily median image for 1 Sep 1996:
MK_DAILY_MED,'c3','960901'
or
MK_DAILY_MED,'c3','1996/09/01'
MODIFICATION HISTORY:
Written by: R.A. Howard, NRL, 9/27/96
Version 2, RAH, 10/10/96, Modified to accept all image sizes
Version 3, RAH, 10/21/96, Added filter and polarizer as inputs
Version 4, SEP, 11/12/96, Fixed problem with equatorial img boundary lines.
Added check for ADCT compressed images.
Version 5, RAH, 6/1/98, Added keyword parameter, QL to select QL or LZ images
Force write to $MONTHLY_IMAGES
Version 6, NBR, 2/11/99, Fix 1999 bug (LE instead of LT 99)
Version 7, NBR, 8/27/99, Re-fix 1999 bug; rebin result to 512x512
Version 8, NBR, 12/1999, Update to use header structure; use READLIST
instead of HEADFITS to filter images; change
jjmax to 100; Add ANYSIZE and NOREBIN keywords; Make b INTARR
NBR, 8. 8.01 - Change output directory to $MONTHLY_IMAGES/daily
NBR, 11. 9.01 - Use IMG_HDR_TXT2STRUCT and fix exptime problem
NBR, 11.13.01 - Add ROLLED keyword
NBR, 11.16.01 - Fix findfile and ih conflict; add FILES keyword
NBR, 11.20.01 - Add CROTA to header
NBR, 12. 7.01 - Header changes
NBR, 12.18.01 - Fix bug in Pass 2
NBR, 5.10.02 - Change jjmax to 50; make date_obs of result midpoint between
date_obs of first and last images used; add N_IMAGES to header;
Use get_sc_point for CROTA; cancel Pass 2 and make n_median_min=5;
Change VERBOSE reporting; Move sdir and postd to end
NBR, 5.23.02 - add caching
NBR, 7.10.02 - Fix P[12]COL/ROW
jake 030716 fixed indentation
NBR, 030718 Automatically check ROLL, even if it changes during a day
version= '@(#)mk_daily_med.pro 1.22, 07/28/03' ; NRL LASCO IDL LIBRARY
;version = 'MK_DAILY_MED testing 2002/05/23 by NBR'
NAME:
MK_DAILY_MIN
PURPOSE:
This procedure generates an image by taking all the files of a given
type for one day and finding the minimum value for each pixel.
CATEGORY:
DATA_ANAL
CALLING SEQUENCE:
MK_DAILY_MIN,Tel,Date
INPUTS:
Tel: String denoting the telescope, 'c1','c2','c3','c4'
Date: Gives the date to be processed in one of the following formats:
YYMMDD, 6 character string
YYYY/MM/DD, 10 character string
CDS Date Structure
Long Word of the modified julian date
OPTIONAL INPUTS:
None
KEYWORD PARAMETERS:
QL: If set then use quick look date, else use level 0 data
FILTER: String denoting the filter position. The default depends on
the telescope:
C1: FeXIV
C2: Orange
C3: Clear
C4: Clear
POLAR: String denoting the polarizer position. The default depends on
the telescope:
C1: Clear
C2: Clear
C3: Clear
C4: 304A
WLL: String denoting the lower FP wavelength of an interval.
WLU: String denoting the upper FP wavelength of an interval.
ONOFF: String denoting whether the on line or off line wavelength should be used
OUTPUTS:
None.
SIDE EFFECTS:
Writes a fits file to $MONTHLY_IMAGES/daily.
RESTRICTIONS:
Only looks for full resolution, full width, any height.
PROCEDURE:
For each image that satifies the selection conditions, (naxis1=1024),
filter and polarizer as requested), the minimum image is computed of
the minimum value of all the images for a single day after being
normalized to the median exposure time.
If the number of images is less than 7 for any pixel, that pixel is
not computed in the first pass. In the second pass, all full images
within +/- 2 days of the given day.
Adapted from MK_DAILY_MED
EXAMPLE:
To create the daily minimum image for 1 Sep 1996:
MK_DAILY_MIN,'c3','960901'
or
MK_DAILY_MIN,'c3','1996/09/01'
MODIFICATION HISTORY:
Written by: R.A. Howard, NRL, 9/27/96
Version 2, RAH, 10/10/96, Modified to accept all image sizes
Version 3, RAH, 10/21/96, Added filter and polarizer as inputs
Version 4, SEP, 11/12/96, Fixed problem with equatorial img boundary lines.
Added check for ADCT compressed images.
Version 5, RAH, 6/1/98, Added keyword parameter, QL to select QL or LZ images
Force write to $MONTHLY_IMAGES
Version 6, NBR, 2/11/99, Fix 1999 bug (LE instead of LT 99)
Version 7, NBR, 8/27/99, Re-fix 1999 bug; rebin result to 512x512
NBR, 8/ 9/01, Change output directory to $MONTHLY_IMAGES/daily
@(#)mk_daily_min.pro 1.3 08/09/01 :LASCO IDL LIBRARY
TITLE:
MK_MONTHLY_MIN
PURPOSE:
This procedure generates the minimum of the daily median
images for an interval (default is 27)
CATEGORY:
DATA_ANAL
CALLING SEQUENCE:
MK_MONTHLY_MIN,Tel,Td
INPUTS:
Tel: String denoting the telescope: 'c1','c2','c3','c4'
Td: String denoting the date, in one of the following
formats: YYYY/MM/DD
YYY-MM-DD
YYMMDD
This date is the MIDPOINT of the desired interval.
OPTIONAL INPUTS:
Filpol: String denoting the filter/polarizer/sector
eg. 'clcl' for the clear/clear combination
The default is 'fvcl', 'orcl' or 'clcl' for C1, C2, C3
respectively.
KEYWORDS:
NDAYS: Number of days to use in the minimum. Must be an odd number.
The default is 27. This keyword overrides the days-available test.
OFF: The string to put after the date if C1
ALL_YEARS:Use daily images from multiple years
NOREBIN: Do not rebin result to 512x512
ROLLED: Use $MONTHLY_IMAGES/rolled directory
FILES Use files in this array to compute minimum
TEST Save in $MONTHLY_IMAGES/testdir
OUTPUTS:
None. The routine will write a FITS file in $MONTHLY_IMAGES. The
format of the filename is tm_fwpw_yymmdd.fts, where
t is the telescope number [1,2,3],
m is the letter 'm' (monthly),
fwpw is the filter/polarizer abbreviation code, and
yymmdd is the date of the MIDPOINT of the interval used to
make the file
PROCEDURE:
- Date of the file output ("DATE-OBS" and filename) is ALWAYS the midpoint of
interval used to make the file and is ALWAYS the input date
- Date of file is always the same DOY (based on 3m_clcl_xx*.fts)
- Searches daily median images in $MONTHLY_IMAGES/daily for Td +/- 13 days.
- Minimum interval is Td +/- 7 days (15 days); if interval is less, no file is
created, unless NDAYS is set.
- Does not span gaps more than 2 days unless NDAYS keyword is set.
MODIFICATION HISTORY:
Written RA Howard, NRL, 9/27/96
15 October 1996 RAH, mods for full images (eliminated regions)
22 Sep 1997 RAH, corrected definition of MID_DATE and MID_TIME
to be mid date and time of the entire interval
06 Oct 1997 SEP, added ALL_YEARS option.
10 Dec 1998 NBR, for non-clear filters, save result as 512x512 image;
divide by 4 for summed images
25 Mar 1999 NBR, insert number of images used in the header
20 Sep 1999 NBR, Save all images as 512x512
Jan 2000 NBR, Add NOREBIN keyword
Jun 2000 NBR, Add ROLL keyword for using *r.fts daily median files
27 Mar 2001 NBR, change mjd0 computation
9 Aug 2001 NBR, Find daily images in $MONTHLY_IMAGES/daily
20 Nov 2001 NBR, Add CROTA keyword to header
25 Feb 2002 NBR, Modify to make files with equal number of days before
and after date of file name; write log file
10 Jul 2002 NBR, Fix PxCOL/ROW, RxCOL/ROW, add keyword TEST
18 Jul 2003 NBR, Fix endutc bug for no days found; add CROTA1/2
10 Oct 2003 NBR, Make ALL_YEARS for 2000+ only
er= '@(#)mk_monthly_min.pro 1.14 04/12/05' ; LASCO IDL LIBRARY
NAME:
READ_CARR_LONG
PURPOSE:
This procedure reads in the carrington number and date information file.
It is called from other routines as needed, and thus doesn't need to be
called by the user explicitly.
CATEGORY:
LASCO Synoptic
CALLING SEQUENCE:
READ_CARR_LONG, Date, Cr, Clong
INPUTS:
None
OPTIONAL INPUTS:
None
OUTPUTS:
The only output is the storing of the carrington information into a common block.
COMMON BLOCKS:
carr_long: Contains the start date of the carrington rotations from
the almanac. Generated by this routine.
EXAMPLE:
MODIFICATION HISTORY:
Written by: RA Howard, NRL, 3/18/96
@(#)read_carr_long.pro 1.4 11/07/96 LASCO IDL LIBRARY
unction rt_carrmapmaker, yymmdd, num_r, rad, limb, wlimb, hdr, disp, saveset, nextrot
PURPOSE:
Restores partially finished Carrington Maps for a given limb/camera
and appends the passed day's worth of images.
INPUTS:
yymmdd STRING: date
num_r INT: number of radii to make maps for-passed empty
rad FLTARR(num_r): radii for maps-passed empty
cmap FLTARR(mapsize,181,num_r): passed undefined
limb INTARR: 0 or 1 tells which limb for carrdate procedure
disp INT: 1 to display images
wlimb STRING: 'wl' or 'el'
hdr STRUCT: dummy header sent back to carrmap3 for fits header
saveset STRING: filename of saved set of carrmapmaker
nextrot INT 1 if going on to the next rotation
OUTPUTS: rad = fltarr(num_r) : radii of maps
num_r : number of radii
RESULT = fltarr(mapsize,181,num_r): carrington maps, one per radius
ROUTINES CALLED:
carrdate2.pro
getc2c3norm2.pro
savestrips.pro
AUTHOR: Nathan Rich, NRL, Nov. 1996
Julia Kraemer, NRL, June 7, 1996
MODIFIED:
11/12/98 NBR modified from carrmapmaker2.pro
12/09/98 NBR do cn=cn+1 in beginning
04/16/99 NBR make into a function
06/08/99 NBR Put check for m LE 0 (median)
11/23/99 NBR Change strip end skip criteria
02/2000 NBR GET_BKG=2 for any_year
12/17/01 NBR Use GET_BKG=1 for current year
12/17/01, @(#)rt_carrmapmaker.pro 1.2
NAME:
SYNSCAN
PURPOSE:
This routine scans an image file and creates a
constant radius scan for synoptic maps.
CATEGORY:
REDUCTION
CALLING SEQUENCE:
SYNSCAN,Filename
INPUTS:
Filename: Ascii string of the name of the
FITS file
OPTIONAL INPUTS:
None
KEYWORD PARAMETERS:
None
OUTPUTS:
None
OPTIONAL OUTPUTS:
None
COMMON BLOCKS:
None
SIDE EFFECTS:
Creates a FITS file of the constant radius
scans in the current directory
RESTRICTIONS:
The directory containing the FITS file should
be the local directory or !imgdir should be set
to point to the directory
PROCEDURE:
EXAMPLE:
MODIFICATION HISTORY:
Written by: RA Howard, NRL, 26 Nov 1995.
Version 1 26 Nov 95 Initial Release
Version 2 13 Apr 96 Modified for keywords
@(#)synscan.pro 1.2 01/19/00 LASCO IDL LIBRARY
NAME: WRITE_SUMMARIES
PURPOSE: Create FITS and GIF image of the summary image
CATEGORY: REDUCTION
CALLING SEQUENCE: WRITE_SUMMARIES,Img,Hdr
INPUTS: Img = Input Image array
Hdr = FITS header
OPTIONAL INPUTS:
Img2 = Input image array for C1 ase image
Hdr2 = FITS header for Img2 array
OUTPUTS: FITS and GIF files are written
RESTRICTIONS: This must be done with IDL having an X window
display.
PROCEDURE: Scales the image to not larger than 512 x 512,
and creates a FITS image with LASCO logo in
image Derived from WRITE_LAST_IMG
EXAMPLE:
MODIFICATION HISTORY: Written, RA Howard, NRL
VERSION 1 rah 12 Apr 1996
VERSION 2 rah 15 Apr 1996
VERSION 3 rah 29 Jun 1996
VERSION 4 rah 29 Oct 1996, Img2 array added to subtract C1 base
VERSION 5 rah 21 Jul 1997, changed scaling on C2/C3 to ratio model
@(#)write_summaries.pro 1.7 11/20/97 :NRL Solar Physics
NAME: WRITE_SUMMARIES
PURPOSE: Create FITS and GIF image of the summary image
CATEGORY: REDUCTION
CALLING SEQUENCE: WRITE_SUMMARIES,Img,Hdr
INPUTS: Img = Input Image array
Hdr = FITS header
OPTIONAL INPUTS: None
OUTPUTS: FITS and GIF files are written
RESTRICTIONS: This must be done with IDL having an X window
display.
PROCEDURE: Scales the image to not larger than 512 x 512,
and creates a FITS image with LASCO logo in
image Derived from WRITE_LAST_IMG
EXAMPLE:
MODIFICATION HISTORY: Written, RA Howard, NRL
VERSION 1 rah 12 Apr 1996
VERSION 2 rah 15 Apr 1996
VERSION 2.1 SPP 28 May 1996 (Replaced DATE-OBS and TIME-OBS
keywords with DATE_OBS and TIME_OBS).
@(#)write_summaries2.pro 1.1 10/05/96 LASCO IDL LIBRARY
lasco_mk_html_help.pro on
Wed Aug 17 12:22:02 2005.