______________________________________________________________________________
CO
COMPARSLINES computes the positions of a given set of spectral
comparison lines for a given set of rows, by nonlinear least squares or
by a centroid calculation (faster but less accurate). The calculation
can be done automatically, or in an interactive mode where the user is
asked to give, for each comparison line, the left and right boundaries
for the fit. This can be done via cursor from a plot device or a
pseudocolor display (preferred), or manually. The user is then asked to
give the set of input rows to use, and the number of rows to bin together
for each solved position. An output SAD file is created which has a
number of columns equal to the number of comparison lines selected and
a number of rows is equal to that of the original map. The value at each
location is the X position of the comparison line. Rows not selected for
calculation have zeroes written to them. The results can be inspected
by plotting columns in option PL.
______________________________________________________________________________
Enter an option: 7-MAY-1995 21:54:33.09
OPTION=.CO
COMPARSLINES, Version 1
Open the file to have the position of comparison lines measured
FILE=CT002028
There is one map in the file
Reading from the file
400 rows of 576 pixels in R4 format
Choose a method for selecting the comparison lines:
A = Automatic
D = Cursor with pseudocolor display (default)
G = Cursor with plot on graphics device
M = Manual input of position
SELECTION_METHOD=A
Choose a mode for calculating the line positions:
C = Centroid
N = Nonlinear least squares (default)
MODE=C
Enter the start row, end row, and number of rows in each solution group:
ROW_PARAMETERS=1 400 25
LINE_WIDTH= 3
FACTOR= 3.00000
19 lines were chosen to be fit between the limits of
132 and 558.
TRY_AGAIN=N
Open the output file to store the line positions
NMAPS= 1 NPLNS= 1 NROWS= 400 NELS= 19
FILE_POSITIONS=TEMP
Writing to an existing file
There is one map in the file
400 rows of 19 pixels in R4 format
MAPS_OUT=1
Please wait for the solution
For the first group of rows:
X = 137.7604
X = 165.5240
X = 305.6377
X = 315.0695
X = 328.2583
X = 342.5175
X = 355.5219
X = 365.0461
X = 383.2852
X = 396.6978
X = 420.5465
X = 439.8323
X = 449.3805
X = 469.6818
X = 478.5776
X = 501.7251
X = 522.8957
X = 539.3678
X = 552.0583
Map number 1 row group 2 is done.
.
.
.
Map number 1 row group 16 is done.
Interpolating and writing to output file; this may take some time
CALIBRATE=T
Open the file of comparsion line wavelengths
FILE_WAVELENGTHS=NE-HG
An initial fit of a small number of lines in the central row will
establish which wavelengths to use.
DEGREE= 3
Enter a sufficient number of column wavelength pairs for the initial
fit: (enter 0 0 to continue)
COLUMN_WAVELENGTH=165.500 5460.700
COLUMN_WAVELENGTH=305.600 5769.600
COLUMN_WAVELENGTH=501.700 6217.300
COLUMN_WAVELENGTH=552.000 6334.400
COLUMN_WAVELENGTH=0 0
Column Wavelength Wavelength(calc.) Difference
165.500 5460.700 5460.700 0.000
305.600 5769.600 5769.599 0.000
501.700 6217.300 6217.299 0.000
552.000 6334.400 6334.400 0.000
The mean error = 0.0000 Angstroms.
The reciprocal dispersion = 2.0583 Angstroms per pixel.
A reference value of the mean error equal to 0.2 Angstroms will be used below.
We will now automatically include as many lines as possible in a second fit.
Enter the degree of the polynomial fit and the multiple of sigma, beyond
which wavelengths are to be rejected: (defaults = 3 2; maximum degree = 7)
DEGREE_MULTIPLE=
3.00000
2.00000
Column Wavelength Wavelength(calc.) Difference
328.423 5820.200 5820.200 0.000
449.736 6096.200 6096.199 0.001
470.063 6143.100 6143.116 -0.016
478.904 6163.600 6163.571 0.029
502.075 6217.300 6217.323 -0.023
523.206 6266.500 6266.502 -0.002
539.608 6304.800 6304.775 0.025
552.288 6334.400 6334.414 -0.014
The mean error = 0.0197 Angstroms.
The reciprocal dispersion = 1.9732 Angstroms per pixel.
Enter a choice for what to do next:
A = Try another fit to all the lines
S = Try another fit to a small number of lines
P = Proceed to fit all the rows with the above degree and
multiple of sigma (default)
LINE_FIT_CHOICE=A
Enter the degree of the polynomial fit and the multiple of sigma, beyond
which wavelengths are to be rejected: (defaults = 3 2; maximum degree = 7)
DEGREE_MULTIPLE=4 2
Column Wavelength Wavelength(calc.) Difference
328.423 5820.200 5820.200 0.000
449.736 6096.200 6096.203 -0.003
470.063 6143.100 6143.113 -0.013
478.904 6163.600 6163.567 0.033
502.075 6217.300 6217.322 -0.022
523.206 6266.500 6266.506 -0.006
539.608 6304.800 6304.778 0.022
552.288 6334.400 6334.410 -0.010
The mean error = 0.0201 Angstroms.
The reciprocal dispersion = 2.2458 Angstroms per pixel.
Enter a choice for what to do next:
A = Try another fit to all the lines
S = Try another fit to a small number of lines
P = Proceed to fit all the rows with the above degree and
multiple of sigma (default)
LINE_FIT_CHOICE=A
Enter the degree of the polynomial fit and the multiple of sigma, beyond
which wavelengths are to be rejected: (defaults = 3 2; maximum degree = 7)
DEGREE_MULTIPLE=5 2
Column Wavelength Wavelength(calc.) Difference
328.423 5820.200 5820.200 0.000
449.736 6096.200 6096.196 0.004
470.063 6143.100 6143.126 -0.026
478.904 6163.600 6163.572 0.028
502.075 6217.300 6217.306 -0.006
523.206 6266.500 6266.505 -0.005
539.608 6304.800 6304.793 0.007
552.288 6334.400 6334.402 -0.002
The mean error = 0.0170 Angstroms.
The reciprocal dispersion = -18.3326 Angstroms per pixel.
Enter a choice for what to do next:
A = Try another fit to all the lines
S = Try another fit to a small number of lines
P = Proceed to fit all the rows with the above degree and
multiple of sigma (default)
LINE_FIT_CHOICE=A
Enter the degree of the polynomial fit and the multiple of sigma, beyond
which wavelengths are to be rejected: (defaults = 3 2; maximum degree = 7)
DEGREE_MULTIPLE=5 3
Column Wavelength Wavelength(calc.) Difference
306.215 5769.600 5769.627 -0.027
315.602 5790.700 5790.841 -0.141
328.423 5820.200 5819.759 0.441
343.039 5852.500 5852.714 -0.214
356.027 5881.900 5882.037 -0.137
383.717 5944.800 5944.813 -0.013
397.160 5975.500 5975.444 0.056
420.948 6030.000 6029.906 0.094
440.256 6074.300 6074.319 -0.019
449.736 6096.200 6096.179 0.021
470.063 6143.100 6143.155 -0.055
478.904 6163.600 6163.618 -0.018
502.075 6217.300 6217.341 -0.041
523.206 6266.500 6266.464 0.036
539.608 6304.800 6304.733 0.067
552.288 6334.400 6334.450 -0.050
The mean error = 0.1671 Angstroms.
The reciprocal dispersion = 9.3850 Angstroms per pixel.
Enter a choice for what to do next:
A = Try another fit to all the lines
S = Try another fit to a small number of lines
P = Proceed to fit all the rows with the above degree and
multiple of sigma (default)
LINE_FIT_CHOICE=
P
A fit will now be made to each row of the file of line positions.
Open a file to store the polynomial parameters for each row
NMAPS= 1 NPLNS= 1 NROWS= 400 NELS= 6
Enter the filename: (<CR> for default)
FILE_POLYNOMIAL=TEMPP
Creating new file
______________________________________________________________________________