adjust-map
adjust-map
adjust-map uses a comparison arc image to make small corrections to the mapping created by map-spectra.
| USAGE | adjust-map -m mapfile -f framename[s] [-d interval]
| ||||||||||||||||||||||||||||||
| INPUT |
| ||||||||||||||||||||||||||||||
| OUTPUT |
| ||||||||||||||||||||||||||||||
| PARAMETERS |
|
Details:
Using an existing spectral map, adjust-map searches in the region of each spectral line to calculate an offset. All offsets for one spectrum are fit to a polynomial which is used to correct the polynomial solutions in the spectral map. If the parameter iterations is greater than 1, the order of the fits is reduced by 1 on the first iteration, and points which deviate from the fit by more than siglimit times the standard deviation are excluded from the next iteration.
The order of the fits must be no greater than the order of the original map:
which is 7 for ord_disp and 5 for ord_sag. Take care if using high order fits:
unless you have very many comparison lines you may make things worse rather
than better. Also, the parameters search_height and search_width should be set
with care. If set too large, the line search may find the incorrect line-
either a neighboring comparison line in the same spectrum, or a line in an
adjacent spectrum. If the initial map errors are comparable to the spacing
between comparison lines and/or spectra, begin with a value of the search
parameters which is sufficiently small so that no misidentifications will
occur, and iterate the solution several times until it converges. Considerable
care should also be devoted to selected a good set of comparison lines. Lines
near to other, brighter lines are particularly to be avoided. If, as
occasionally happens, the ends of the slit are badly cut, edge may be set to a
non-zero value to exclude the slit ends in the adjustments in the wavelength
direction.
The quality of the fits may be tested by setting the diagnostic parameter to a
value other than zero, in which case each nth fit is plotted. Green points are
those which were used on the final iteration; red points were excluded.
Inspection of these plots will help the user set the optimal order for the
fits, and also look for any troublesome comparison lines. In diagnostic mode no
output files are written. If not run in diagnostic mode, plots of the fits to
all spectra are output into a file named mapfile.ps
If the order of the adjustment to either fit is set to a negative number, no correction of that fit is made.
It is often useful to interpolate between two comparison arc exposures which, for example, bracket an object spectrum. If two framenames are provided, the map adjustment is the mean of that of the two frames. The output mapfile name is a composite of the two framenames; e.g. from ccd0029 and ccd0030, a mapfile ccd0029-30.map is produced.
Suggested fit order is 3 to 4 for ord_disp, 1 to 3 for ord_sag, but the order
for ord_disp may need to be lowered if there are insufficient comparison lines.
Two options are available to improve the quality of the fits:
- In the line list file, if a second number follows the wavelength, it is used as a weight for the line in the fitting process. Lines without a weight are assigned weights of unity.
- If
outlier-wtis set to y, a nearest neighbor algorithm is used to identify outlying points, which might be due to misidentified lines (usually caused by interference with a neighboring spectrum). Points which are far from other points are weighted inversely by the square of their distance to neighbors.
| ADVANCED PARAMETERS |
|
Important note: It is strongly recommended that bias-subtracted frames be used with adjust-map. Slit positions of weak lines in frames without bias subtraction may be less accurate.