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
|
include <imhdr.h>
# AP_GETDATA -- Get the summed dispersion line.
# Return the IMIO pointer, pointer to image data, the aperture axis and title.
# The pointers must be freed by the calling program. Note that the value of
# line may be changed.
procedure ap_getdata (image, line, nsum, im, imdata, npts, apaxis, title)
char image[SZ_FNAME] # Image name
int line # Dispersion line to graph
int nsum # Number of dispersion lines to sum
pointer im # IMIO pointer
pointer imdata # Pointer to image data
int npts # Number of pixels
int apaxis # Aperture axis
pointer title # Title for image data
int i, j, k, l, n, dispaxis
pointer buf, medbuf
real asumr(), amedr()
pointer ap_immap(), imgs2r()
errchk ap_immap, imgs2r
begin
# Map the image
im = ap_immap (image, apaxis, dispaxis)
# Determine the dispersion and aperture axes.
if (IS_INDEFI (line))
line = IM_LEN(im, dispaxis) / 2
else
line = max (1, min (IM_LEN(im, dispaxis), line))
# Allocate memory for the image line and title.
npts = IM_LEN(im, apaxis)
call calloc (imdata, npts, TY_REAL)
call malloc (title, SZ_LINE, TY_CHAR)
# Sum the specified number of dispersion lines.
n = max (1, abs (nsum))
switch (apaxis) {
case 1:
i = max (1, line - n / 2)
j = min (IM_LEN(im, dispaxis), i + n - 1)
i = max (1, j - n + 1)
buf = imgs2r (im, 1, npts, i, j)
j = j - i + 1
if (j < 3 || nsum > 0) {
do k = 1, j
call aaddr (Memr[buf+(k-1)*npts], Memr[imdata],
Memr[imdata], npts)
call sprintf (Memc[title], SZ_LINE,
"Image=%s, Sum of lines %d-%d")
call pargstr (image)
call pargi (i)
call pargi (i+j-1)
} else {
call malloc (medbuf, j, TY_REAL)
do k = 0, npts-1 {
do l = 0, j-1
Memr[medbuf+l] = Memr[buf+l*npts+k]
Memr[imdata+k] = amedr (Memr[medbuf], j)
}
call mfree (medbuf, TY_REAL)
call sprintf (Memc[title], SZ_LINE,
"Image=%s, Median of lines %d-%d")
call pargstr (image)
call pargi (i)
call pargi (i+j-1)
}
case 2:
i = max (1, line - n / 2)
j = min (IM_LEN(im, dispaxis), i + n - 1)
i = max (1, j - n + 1)
buf = imgs2r (im, i, j, 1, npts)
j = j - i + 1
if (j < 3 || nsum > 0) {
do k = 1, npts
Memr[imdata+k-1] = asumr (Memr[buf+(k-1)*j], j)
call sprintf (Memc[title], SZ_LINE,
"Image=%s, Sum of columns %d-%d")
call pargstr (image)
call pargi (i)
call pargi (i+j-1)
} else {
do k = 1, npts
Memr[imdata+k-1] = amedr (Memr[buf+(k-1)*j], j)
call sprintf (Memc[title], SZ_LINE,
"Image=%s, Median of columns %d-%d")
call pargstr (image)
call pargi (i)
call pargi (i+j-1)
}
}
end
|
