______________________________________________________________________________ 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 ______________________________________________________________________________