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
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
|
include "apertures.h"
define NRANGES 50
# AP_RECENTER -- Recenter apertures.
procedure ap_recenter (image, line, nsum, aps, naps, apedit)
char image[SZ_FNAME] # Image name
int line # Image dispersion line
int nsum # Number of dispersion lines to sum
pointer aps[ARB] # Aperture pointers
int naps # Number of apertures
int apedit # Called by apedit?
pointer ranges # Apertures to select
int npeaks # Number of bright peaks to select
bool shift # Shift instead of center?
real center, delta
int i, j, k, na, npts, apaxis
pointer sp, str, im, imdata, title, index, peaks, deltas
int decode_ranges()
real apgetr(), ap_center(), ap_cveval(), asokr()
bool clgetb(), ap_answer(), apgetb(), is_in_range()
errchk ap_getdata
begin
# Check if apertures are defined.
na = 0
do i = 1, naps
if (AP_SELECT(aps[i]) == YES)
na = na + 1
if (na < 1)
return
# Query user.
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
if (apedit == NO) {
call sprintf (Memc[str], SZ_LINE, "Recenter apertures for %s?")
call pargstr (image)
if (!ap_answer ("ansrecenter", Memc[str])) {
call sfree (sp)
return
}
if (clgetb ("verbose"))
call printf ("Recentering apertures ...\n")
}
# Get parameters
delta = apgetr ("npeaks")
shift = apgetb ("shift")
if (IS_INDEFR (delta))
npeaks = naps
else if (delta < 1.)
npeaks = max (1., delta * naps)
else
npeaks = delta
# Map the image and get the image data.
call ap_getdata (image, line, nsum, im, imdata, npts, apaxis, title)
if (npeaks == naps && !shift) {
na = 0
do i = 1, naps {
if (AP_SELECT(aps[i]) == NO)
next
center = AP_CEN(aps[i], apaxis) +
ap_cveval (AP_CV(aps[i]), real (line))
center = ap_center (center, Memr[imdata], npts)
if (!IS_INDEF(center)) {
AP_CEN(aps[i], apaxis) = center -
ap_cveval (AP_CV(aps[i]), real (line))
na = na + 1
}
}
} else {
call salloc (ranges, 3*NRANGES, TY_INT)
call salloc (index, naps, TY_REAL)
call salloc (peaks, naps, TY_REAL)
call salloc (deltas, naps, TY_REAL)
call apgstr ("aprecenter", Memc[str], SZ_LINE)
if (decode_ranges (Memc[str], Memi[ranges], NRANGES, i) == ERR)
call error (0, "Bad aperture list")
j = 0
do i = 1, naps {
if (!is_in_range (Memi[ranges], AP_ID(aps[i])))
next
center = AP_CEN(aps[i], apaxis) +
ap_cveval (AP_CV(aps[i]), real (line))
delta = ap_center (center, Memr[imdata], npts)
if (!IS_INDEF(delta)) {
k = max (1, min (npts, int (delta+0.5)))
Memr[index+j] = i
Memr[peaks+j] = -Memr[imdata+k-1]
Memr[deltas+j] = delta - center
j = j + 1
}
}
if (j > 0 && npeaks > 0) {
if (npeaks < j) {
call xt_sort3 (Memr[peaks], Memr[deltas], Memr[index], j)
j = npeaks
}
if (shift) {
if (mod (j, 2) == 0)
delta = (asokr (Memr[deltas], j, j/2) +
asokr (Memr[deltas], j, 1+j/2)) / 2
else
delta = asokr (Memr[deltas], j, 1+j/2)
na = 0
do i = 1, naps {
if (AP_SELECT(aps[i]) == NO)
next
center = AP_CEN(aps[i], apaxis) + delta
AP_CEN(aps[i], apaxis) = center
na = na + 1
}
} else {
na = 0
do k = 1, j {
delta = Memr[deltas+k-1]
i = Memr[index+k-1]
if (AP_SELECT(aps[i]) == NO)
next
center = AP_CEN(aps[i], apaxis) + delta
AP_CEN(aps[i], apaxis) = center
na = na + 1
}
}
}
}
# Log the operation, write the apertures to the database,
# unmap the image and free memory.
if (shift) {
call sprintf (Memc[str], SZ_LINE,
"RECENTER - %d apertures shifted by %.2f for %s.")
call pargi (na)
call pargr (delta)
call pargstr (image)
} else {
call sprintf (Memc[str], SZ_LINE,
"RECENTER - %d apertures recentered for %s")
call pargi (na)
call pargstr (image)
}
if (apedit == NO)
call ap_log (Memc[str], YES, YES, NO)
else
call ap_log (Memc[str], YES, NO, NO)
call appstr ("ansdbwrite1", "yes")
call mfree (imdata, TY_REAL)
call mfree (title, TY_CHAR)
call imunmap (im)
call sfree (sp)
end
|
