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
|
include <mach.h>
include "epix.h"
# EP_SEARCH -- Search input data for maximum or minimum pixel in search radius.
# Return the new aperture positions. The magnitude of the search radius
# defines the range to be searched (bounded by the raster dimension) and
# the sign of the radius determines whether a minimum or maximum is sought.
procedure ep_search (ep, data, nx, ny, ap, xa, ya, xb, yb)
pointer ep # EPIX pointer
real data[nx,ny] # Subraster
int nx, ny # Subraster size
int ap # Aperture type
int xa, ya, xb, yb # Aperture (initial and final)
real xc, yc, search2, dj2, r2, dmax
int i, j, i1, i2, j1, j2, imax, jmax
begin
if (EP_SEARCH(ep) == 0.)
return
search2 = abs (EP_SEARCH(ep))
xa = xa - EP_X1(ep) + 1
xb = xb - EP_X1(ep) + 1
xc = (xa + xb) / 2.
i1 = max (1., xc - search2)
i2 = min (real(nx), xc + search2)
imax = nint (xc)
ya = ya - EP_Y1(ep) + 1
yb = yb - EP_Y1(ep) + 1
yc = (ya + yb) / 2.
j1 = max (1., yc - search2)
j2 = min (real(ny), yc + search2)
jmax = nint (yc)
dmax = data[imax,jmax]
switch (ap) {
case 1:
search2 = EP_SEARCH(ep) ** 2
do j = j1, j2 {
dj2 = (j - yc) ** 2
do i = i1, i2 {
r2 = dj2 + (i - xc) ** 2
if (r2 > search2)
next
if (EP_SEARCH(ep) > 0.) {
if (data[i,j] > dmax) {
dmax = data[i,j]
imax = i
jmax = j
}
} else {
if (data[i,j] < dmax) {
dmax = data[i,j]
imax = i
jmax = j
}
}
}
}
default:
do j = j1, j2 {
do i = i1, i2 {
if (EP_SEARCH(ep) > 0.) {
if (data[i,j] > dmax) {
dmax = data[i,j]
imax = i
jmax = j
}
} else {
if (data[i,j] < dmax) {
dmax = data[i,j]
imax = i
jmax = j
}
}
}
}
}
xa = xa + (imax - xc) + EP_X1(ep) - 1
xb = xb + (imax - xc) + EP_X1(ep) - 1
ya = ya + (jmax - yc) + EP_Y1(ep) - 1
yb = yb + (jmax - yc) + EP_Y1(ep) - 1
end
|
