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
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
|
include <error.h>
include <imhdr.h>
include <smw.h>
# SFLIP -- Flip data and/or coordinate system in spectra.
procedure t_sflip ()
pointer inlist # Input list
pointer outlist # Output list
bool coord_flip # Flip coordinates?
bool data_flip # Flip data?
bool in_place
int i, j, k, n, axis
pointer sp, input, output, temp, a, b, c
pointer in, out, smw, mw, tmp, inbuf, outbuf
bool clgetb(), streq()
int imtgetim(), imtlen()
pointer imtopenp(), immap(), smw_openim(), imgl3r(), impl3r()
errchk immap, smw_openim
begin
call smark (sp)
call salloc (input, SZ_FNAME, TY_CHAR)
call salloc (output, SZ_FNAME, TY_CHAR)
call salloc (temp, SZ_FNAME, TY_CHAR)
call salloc (a, 3*3, TY_DOUBLE)
call salloc (b, 3, TY_DOUBLE)
call salloc (c, 3, TY_DOUBLE)
# Get task parameters.
inlist = imtopenp ("input")
outlist = imtopenp ("output")
coord_flip = clgetb ("coord_flip")
data_flip = clgetb ("data_flip")
# Loop over all input images.
in = NULL
out = NULL
smw = NULL
while (imtgetim (inlist, Memc[input], SZ_FNAME) != EOF) {
if (imtlen (outlist) > 0) {
if (imtgetim (outlist, Memc[output], SZ_FNAME) == EOF)
break
} else
call strcpy (Memc[input], Memc[output], SZ_FNAME)
if (streq (Memc[input], Memc[output])) {
if (data_flip) {
in_place = false
call mktemp ("temp", Memc[temp], SZ_FNAME)
} else
in_place = true
} else {
in_place = false
call strcpy (Memc[output], Memc[temp], SZ_FNAME)
}
iferr {
# Map the images and WCS.
if (in_place) {
tmp = immap (Memc[input], READ_WRITE, 0); in = tmp
out = in
} else {
tmp = immap (Memc[input], READ_ONLY, 0); in = tmp
tmp = immap (Memc[temp], NEW_COPY, in); out = tmp
}
tmp = smw_openim (in); smw = tmp
# Flip coordinates.
if (coord_flip) {
mw = SMW_MW(smw,0)
n = SMW_PDIM(smw)
axis = SMW_PAXIS(smw,1) - 1
call mw_gltermd (mw, Memd[a], Memd[b], n)
Memd[a+axis*(n+1)] = -Memd[a+axis*(n+1)]
Memd[b+axis] = SMW_LLEN(smw,1) - Memd[b+axis] + 1
call mw_sltermd (mw, Memd[a], Memd[b], n)
call smw_saveim (smw, out)
}
# Flip data.
if (data_flip) {
n = IM_LEN(in,1)
do j = 1, IM_LEN(in,3) {
do i = 1, IM_LEN(in,2) {
inbuf = imgl3r (in, i, j)
switch (SMW_FORMAT(smw)) {
case SMW_ND:
switch (SMW_LAXIS(smw,1)) {
case 1:
outbuf = impl3r (out, i, j) + n - 1
do k = 0, n-1
Memr[outbuf-k] = Memr[inbuf+k]
case 2:
outbuf = impl3r (out, IM_LEN(in,2)-i+1, j)
call amovr (Memr[inbuf], Memr[outbuf], n)
case 3:
outbuf = impl3r (out, i, IM_LEN(in,3)-j+1)
call amovr (Memr[inbuf], Memr[outbuf], n)
}
case SMW_ES, SMW_MS:
outbuf = impl3r (out, i, j) + n - 1
do k = 0, n-1
Memr[outbuf-k] = Memr[inbuf+k]
}
}
}
} else if (!in_place) {
n = IM_LEN(in,1)
do j = 1, IM_LEN(in,3) {
do i = 1, IM_LEN(in,2) {
inbuf = imgl3r (in, i, j)
outbuf = impl3r (out, i, j)
call amovr (Memr[inbuf], Memr[outbuf], n)
}
}
}
} then {
if (!in_place && out != NULL) {
call imunmap (out)
call imdelete (Memc[temp])
}
call erract (EA_WARN)
}
if (smw != NULL)
call smw_close (smw)
if (!in_place && out != NULL) {
call imunmap (out)
call imunmap (in)
if (streq (Memc[input], Memc[output])) {
call imdelete (Memc[input])
call imrename (Memc[temp], Memc[output])
}
} else if (in != NULL)
call imunmap (in)
}
call imtclose (inlist)
call imtclose (outlist)
call sfree (sp)
end
|
