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
|
include <imhdr.h>
include <imio.h>
include <mwset.h>
define SYSTEMS "|equispec|multispec|physical|image|world|linear|"
# SMW_OPENIM -- Open the spectral MWCS for various input formats.
pointer procedure smw_openim (im)
pointer im #I Image pointer
pointer mw #O MWCS pointer
pointer sp, system, mw_openim()
bool streq()
int i, wcsdim, sys, strdic(), mw_stati()
errchk mw_openim, smw_oldms, smw_linear
begin
call smark (sp)
call salloc (system, SZ_FNAME, TY_CHAR)
# Workaround for truncation of header during image header copy.
IM_HDRLEN(im) = IM_LENHDRMEM(im)
# Force higher dimensions to length 1.
do i = IM_NDIM(im) + 1, 3
IM_LEN(im,i) = 1
mw = mw_openim (im)
call mw_seti (mw, MW_USEAXMAP, NO)
wcsdim = mw_stati (mw, MW_NDIM)
call mw_gwattrs (mw, 0, "system", Memc[system], SZ_FNAME)
sys = strdic (Memc[system], Memc[system], SZ_FNAME, SYSTEMS)
# Set various input systems.
switch (sys) {
case 1:
call smw_equispec (im, mw)
case 2:
call smw_multispec (im, mw)
default:
if (sys == 0) {
call eprintf (
"WARNING: Unknown coordinate system `%s' - assuming `linear'.\n")
call pargstr (Memc[system])
} else if (sys == 3)
call mw_newsystem (mw, "image", wcsdim)
# Old "multispec" format.
ifnoerr (call imgstr (im, "APFORMAT", Memc[system], SZ_FNAME)) {
if (streq (Memc[system], "onedspec"))
call smw_onedspec (im, mw)
else
call smw_oldms (im, mw)
# Old "onedspec" format or other 1D image.
} else if (wcsdim == 1) {
call smw_onedspec (im, mw)
# N-dimensional image.
} else
call smw_nd (im, mw)
}
call sfree (sp)
return (mw)
end
|
