1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
|
include <imhdr.h>
include "ms.h"
# T_MSEXTRACT -- General MULTISPEC extraction task.
#
# The general task parameters are obtained and the desired extraction
# procedure is called. The input database and image are accessed and
# the output image is created.
procedure t_msextract ()
# User parameters:
char image[SZ_FNAME] # Image
char output[SZ_FNAME] # Output image file
real lower # Lower limit of strip
real upper # Upper limit of strip
int spectra[3, MAX_RANGES] # Spectra to be extracted
int lines[3, MAX_RANGES] # Lines to be extracted
bool ex_model # Extract model or data
bool integrated # Extract integrated spectra?
bool unblend # Correct for spectra blending
bool clean # Correct for bad pixels
int nreplace # Maximum number pixels replaced
real sigma_cut # Threshold for replacing bad pixels
int model # Model type: gauss5, profile
bool ex_spectra
int nlines
int nspectra
pointer ms, im_in, im_out
int clgeti(), ms_model_id(), clgranges()
bool clgetb()
real clgetr()
pointer msmap(), immap()
begin
# Access input and output files.
call clgstr ("image", image, SZ_FNAME)
ms = msmap (image, READ_ONLY, 0)
im_in = immap (image, READ_ONLY, 0)
call clgstr ("output", output, SZ_FNAME)
im_out = immap (output, NEW_IMAGE, 0)
# Determine extraction limits.
lower = clgetr ("lower")
upper = clgetr ("upper")
nlines = clgranges ("lines", 1, IM_LEN(im_in, 2), lines, MAX_RANGES)
nspectra = clgranges ("spectra", 1, MS_NSPECTRA(ms), spectra,
MAX_RANGES)
# Determine type of extraction.
ex_spectra = TRUE
ex_model = clgetb ("ex_model")
integrated = clgetb ("integrated")
# Determine whether to clean data spectra and the cleaning parameters.
clean = clgetb ("clean")
if (clean) {
nreplace = clgeti ("nreplace")
sigma_cut = clgetr ("sigma_cut")
} else
nreplace = 0
# Determine whether to apply blending correction.
if (!ex_model)
unblend = clgetb ("unblend")
# Set type of model to be used. If a blending correction is desired
# the model must GAUSS5 otherwise the user selects the model.
model = NONE
if (unblend)
model = GAUSS5
else if (ex_model || clean)
model = ms_model_id ("model")
# Set verbose output.
call ex_set_verbose (clgetb ("verbose"))
call ex_prnt1 (MS_IMAGE(ms), output)
# Set image header for output extraction image file.
IM_NDIM(im_out) = 3
if (integrated)
IM_LEN(im_out, 1) = 1
else
IM_LEN(im_out, 1) = nint (upper - lower + 1)
IM_LEN(im_out, 2) = nlines
IM_LEN(im_out, 3) = nspectra
IM_PIXTYPE(im_out) = TY_REAL
call strcpy (IM_TITLE(im_in), IM_TITLE(im_out), SZ_IMTITLE)
# Select extraction procedure based on model.
switch (model) {
case GAUSS5:
call set_fit_and_clean (clgeti ("niterate"), nreplace, sigma_cut,
clgeti ("fit_type"), ex_model)
call ex_gauss5 (ms, im_in, im_out, spectra, lines, lower, upper,
ex_spectra, ex_model, integrated)
case SMOOTH:
call set_fit_smooth (nreplace, sigma_cut)
call ex_smooth (ms, im_in, im_out, spectra, lines, lower, upper,
ex_spectra, ex_model, integrated)
default:
call ex_strip (ms, im_in, im_out, spectra, lines, lower, upper,
ex_spectra, ex_model, integrated)
}
# Close files.
call imunmap (im_in)
call imunmap (im_out)
call msunmap (ms)
end
|
