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>
include <imset.h>
include "ccdred.h"
# SET_TRIM -- Set the trim parameters.
#
# 1. Return immediately if the trim correction is not requested or
# if the image has been previously corrected.
# 2. Determine the trim section. This may be specifed directly or
# indirectly through the image header or symbol table.
# 3. Parse the trim section and apply it to the output image.
# 4. If the image is trimmed then log the operation and reset the output
# image size.
procedure set_trim (ccd)
pointer ccd # CCD structure
int xt1, xt2, yt1, yt2
int nc, nl, c1, c2, l1, l2
pointer sp, str, image
bool clgetb(), ccdflag()
begin
# Check if the user wants this operation or it has been done.
if (!clgetb ("trim") || ccdflag (IN_IM(ccd), "trim"))
return
# Check trim section.
nc = IM_LEN(IN_IM(ccd),1)
nl = IM_LEN(IN_IM(ccd),2)
c1 = TRIM_C1(ccd)
c2 = TRIM_C2(ccd)
l1 = TRIM_L1(ccd)
l2 = TRIM_L2(ccd)
if ((c1 < 1) || (c2 > nc) || (l1 < 1) || (l2 > nl)) {
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
call salloc (image, SZ_LINE, TY_CHAR)
call imstats (IN_IM(ccd), IM_IMAGENAME, Memc[image], SZ_FNAME)
call sprintf (Memc[str], SZ_LINE,
"Error in trim section: image=%s[%d,%d], trimsec=[%d:%d,%d:%d]")
call pargstr (Memc[image])
call pargi (nc)
call pargi (nl)
call pargi (c1)
call pargi (c2)
call pargi (l1)
call pargi (l2)
call error (0, Memc[str])
}
# If no processing is desired print trim section and return.
if (clgetb ("noproc")) {
call eprintf (" [TO BE DONE] Trim section is [%d:%d,%d:%d].\n")
call pargi (c1)
call pargi (c2)
call pargi (l1)
call pargi (l2)
return
}
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
xt1 = max (0, c1 - IN_C1(ccd))
xt2 = min (0, c2 - IN_C2(ccd))
yt1 = max (0, l1 - IN_L1(ccd))
yt2 = min (0, l2 - IN_L2(ccd))
CCD_C1(ccd) = CCD_C1(ccd) + xt1
CCD_C2(ccd) = CCD_C2(ccd) + xt2
CCD_L1(ccd) = CCD_L1(ccd) + yt1
CCD_L2(ccd) = CCD_L2(ccd) + yt2
IN_C1(ccd) = IN_C1(ccd) + xt1
IN_C2(ccd) = IN_C2(ccd) + xt2
IN_L1(ccd) = IN_L1(ccd) + yt1
IN_L2(ccd) = IN_L2(ccd) + yt2
OUT_C1(ccd) = IN_C1(ccd) - c1 + 1
OUT_C2(ccd) = IN_C2(ccd) - c1 + 1
OUT_L1(ccd) = IN_L1(ccd) - l1 + 1
OUT_L2(ccd) = IN_L2(ccd) - l1 + 1
IM_LEN(OUT_IM(ccd),1) = c2 - c1 + 1
IM_LEN(OUT_IM(ccd),2) = l2 - l1 + 1
CORS(ccd, TRIM) = YES
COR(ccd) = YES
call sprintf (Memc[str], SZ_LINE, "Trim data section is [%d:%d,%d:%d]")
call pargi (c1)
call pargi (c2)
call pargi (l1)
call pargi (l2)
call timelog (Memc[str], SZ_LINE)
call ccdlog (IN_IM(ccd), Memc[str])
call hdmpstr (OUT_IM(ccd), "trim", Memc[str])
call sfree (sp)
end
|
