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
|
# IMALIGN - Register a set of images using the results of the IMCENTROID,
# and the IMSHIFT, and IMCOPY tasks.
procedure imalign (input, reference, coords, output)
begin
bool shifts_found, trim_found
string tmpfile, outfile, shiftfile, trimsect, tmp, junk
string l_input, l_reference, l_coords, l_output
int x1, x2, y1, y2
real xshift, yshift
struct line
# Set up some temporary files.
tmpfile = mktemp ("tmp$ia_tmp.")
outfile = mktemp ("tmp$ia_tmp.")
shiftfile = mktemp ("tmp$ia_tmp.")
# Get the required parameters.
l_input = input
l_reference = reference
l_coords = coords
l_output = output
# Write the output names to outfile.
sections (l_output, option="fullname", > outfile)
# Compute the centers and relative shifts.
imcentroid (l_input, l_reference, l_coords, shifts=shifts,
boxsize=boxsize, bigbox=bigbox, negative=negative,
background=background, lower=lower, upper=upper,
niterate=niterate, tolerance=tolerance, maxshift=maxshift,
verbose=verbose, >& tmpfile)
# Print the centering results on the screen?
if (verbose)
type (tmpfile)
# Shift the images.
if (shiftimages) {
# Read the shifts.
shifts_found = no
list = tmpfile
while (fscan (list, line) != EOF) {
tmp = substr (line, 2, 7)
if (tmp == "Shifts") {
shifts_found = yes
break
}
}
# Decode the shifts.
if (shifts_found)
while (fscan (list, junk, xshift, junk, yshift, junk) == 5)
print (xshift, " ", yshift, >> shiftfile)
else
error (1, "No shifts were calculated.")
# Shift the images.
print ("\n# Shifting images:\n")
imshift (l_input, "@"//outfile, shifts_file=shiftfile,
interp_type=interp_type, boundary_type=boundary_type,
constant=constant)
# Trim the images.
if (trimimages) {
# Check for vignetting.
trim_found = no
while (fscanf (list, "%s = [%d:%d,%d:%d]", line, x1, x2,
y1, y2) != EOF) {
tmp = substr (line, 2, 5)
if (tmp == "Vign") {
print ("Images not trimmed ! Vignetting is present.")
trim_found = no
break
} else if (tmp == "Trim") {
trim_found = yes
break
}
}
# Trim the images.
if (!trim_found) {
print ("Images not trimmed ! Trim section is undefined.")
} else {
# Correct for boundary extension "contamination".
if (interp_type == "poly3") {
x1 += 1; x2 -= 1; y1 += 1; y2 -= 1
} else if (interp_type == "poly5" ||
interp_type == "spline3") {
x1 += 2; x2 -= 2; y1 += 2; y2 -= 2
}
if (1 <= x1 && x1 <= x2 && 1 <= y1 && y1 <= y2) {
trimsect = "["//x1//":"//x2//","//y1//":"//y2//"]"
list = outfile; delete (tmpfile, ver-, >& "dev$null")
while (fscan (list, tmp) != EOF)
print (tmp//trimsect, >> tmpfile)
print ("# Trimming images: corrected section = ",
trimsect)
imcopy ("@"//tmpfile, "@"//outfile, verbose-)
} else {
print ("Images not trimmed ! No overlap region.")
}
}
}
}
list = ""
delete (tmpfile, ver-, >& "dev$null")
delete (outfile, ver-, >& "dev$null")
delete (shiftfile, ver-, >& "dev$null")
end
|
