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
|
include <imhdr.h>
include <fset.h>
include "iralign.h"
# T_IRMATCHD2 -- Align the individual subraster elements in the input image.
# In order to run this program the user should have created the output image
# and the database file with the IRMOSAIC task. In addition the user should
# supply a coordinate list consisting of pairs of coordinates of identical
# objects or features in two adjacent subrasters.
procedure t_irmatchd2 ()
int cl, interp, align, verbose, nimages, nmatch
pointer sp, inimage, outimage, database, coords, matchlist, trimlimits, ranges
pointer str, ir, im, outim, dt
bool clgetb()
int open(), clgwrd(), btoi(), ir_links(), ir_clinks(), ir_flinks()
int decode_ranges()
pointer immap(), dtmap()
real clgetr()
begin
# Allocate space for the ir strucuture.
call ir_init (ir)
# Set the standard output to flush on a new line.
call fseti (STDOUT, F_FLUSHNL, YES)
# Allocate temporary working space.
call smark (sp)
call salloc (inimage, SZ_FNAME, TY_CHAR)
call salloc (outimage, SZ_FNAME, TY_CHAR)
call salloc (coords, SZ_FNAME, TY_CHAR)
call salloc (database, SZ_FNAME, TY_CHAR)
call salloc (matchlist, SZ_LINE, TY_CHAR)
call salloc (trimlimits, SZ_FNAME, TY_CHAR)
call salloc (ranges, 3 * MAX_NRANGES + 1, TY_INT)
call salloc (str, SZ_LINE, TY_CHAR)
# Get the input and output images and the coordinate list.
call clgstr ("input", Memc[inimage], SZ_FNAME)
call clgstr ("output", Memc[outimage], SZ_FNAME)
call clgstr ("database", Memc[database], SZ_FNAME)
call clgstr ("coords", Memc[coords], SZ_FNAME)
align = clgwrd ("alignment", Memc[str], SZ_LINE, ",coords,shifts,file,")
call clgstr ("match", Memc[matchlist], SZ_LINE)
call clgstr ("trimlimits", Memc[trimlimits], SZ_FNAME)
# Open the images and files.
im = immap (Memc[inimage], READ_ONLY, 0)
outim = immap (Memc[outimage], NEW_COPY, im)
dt = dtmap (Memc[database], READ_ONLY)
# Get the data base parameters.
call ir_dtrparams (dt, Memc[inimage], ir)
call ir_params (ir, im, outim)
interp = clgwrd ("interpolant", Memc[str], SZ_LINE,
",nearest,linear,poly3,poly5,spline3,")
verbose = btoi (clgetb ("verbose"))
# Decode the list of input images to be intensity matched.
nimages = IR_NXSUB(ir) * IR_NYSUB(ir)
if (Memc[matchlist] == EOS) {
Memi[ranges] = NULL
} else if (Memc[matchlist] == '*') {
Memi[ranges] = 1
Memi[ranges+1] = nimages
Memi[ranges+2] = 1
Memi[ranges+3] = NULL
} else if (decode_ranges (Memc[matchlist], Memi[ranges], MAX_NRANGES,
nmatch) == ERR) {
call error (0,
"Cannot decode list of rasters to be intensity matched.")
}
# Allocate working space.
call ir_arrays (ir, nimages)
# Compute the shifts for each subraster.
switch (align) {
case IR_COORDS:
cl = open (Memc[coords], READ_ONLY, TEXT_FILE)
if (ir_links (cl, Memr[IR_XRSHIFTS(ir)], Memr[IR_YRSHIFTS(ir)],
Memr[IR_XCSHIFTS(ir)], Memr[IR_YCSHIFTS(ir)],
Memi[IR_NRSHIFTS(ir)], Memi[IR_NCSHIFTS(ir)],
IR_NCOLS(ir), IR_NROWS(ir), IR_NXRSUB(ir), IR_NYRSUB(ir),
IR_NXSUB(ir), IR_NYSUB(ir), IR_NXOVERLAP(ir), IR_NYOVERLAP(ir),
IR_ORDER(ir)) > 0) {
call ir_shifts (ir, im, outim, Memr[IR_XRSHIFTS(ir)],
Memr[IR_YRSHIFTS(ir)], Memr[IR_XCSHIFTS(ir)],
Memr[IR_YCSHIFTS(ir)], Memi[IR_IC1(ir)], Memi[IR_IC2(ir)],
Memi[IR_IL1(ir)], Memi[IR_IL2(ir)], Memi[IR_OC1(ir)],
Memi[IR_OC2(ir)], Memi[IR_OL1(ir)], Memi[IR_OL2(ir)],
Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)])
call ir_m2match (ir, im, Memi[ranges], Memi[IR_IC1(ir)],
Memi[IR_IC2(ir)], Memi[IR_IL1(ir)], Memi[IR_IL2(ir)],
Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
Memr[IR_DELTAI(ir)])
} else
call error (0, "There are no legal shifts in the coords file.")
call close (cl)
case IR_SHIFTS:
if (ir_clinks (Memr[IR_XRSHIFTS(ir)], Memr[IR_YRSHIFTS(ir)],
Memr[IR_XCSHIFTS(ir)], Memr[IR_YCSHIFTS(ir)], IR_NXRSUB(ir),
IR_NYRSUB(ir), IR_NXSUB(ir), IR_NYSUB(ir), clgetr ("xshift"),
clgetr ("yshift")) > 0) {
call ir_shifts (ir, im, outim, Memr[IR_XRSHIFTS(ir)],
Memr[IR_YRSHIFTS(ir)], Memr[IR_XCSHIFTS(ir)],
Memr[IR_YCSHIFTS(ir)], Memi[IR_IC1(ir)], Memi[IR_IC2(ir)],
Memi[IR_IL1(ir)], Memi[IR_IL2(ir)], Memi[IR_OC1(ir)],
Memi[IR_OC2(ir)], Memi[IR_OL1(ir)], Memi[IR_OL2(ir)],
Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)])
call ir_m2match (ir, im, Memi[ranges], Memi[IR_IC1(ir)],
Memi[IR_IC2(ir)], Memi[IR_IL1(ir)], Memi[IR_IL2(ir)],
Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
Memr[IR_DELTAI(ir)])
} else
call error (0, "There are no legal shifts in the coords file.")
case IR_FILE:
cl = open (Memc[coords], READ_ONLY, TEXT_FILE)
if (ir_flinks (cl, Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
Memr[IR_DELTAI(ir)], nimages) >= nimages) {
call ir_fshifts (ir, im, outim, Memr[IR_DELTAX(ir)],
Memr[IR_DELTAY(ir)], Memi[IR_IC1(ir)], Memi[IR_IC2(ir)],
Memi[IR_IL1(ir)], Memi[IR_IL2(ir)], Memi[IR_OC1(ir)],
Memi[IR_OC2(ir)], Memi[IR_OL1(ir)], Memi[IR_OL2(ir)])
} else
call error (0, "There are fewer shifts than input subrasters.")
call close (cl)
default:
call error (0, "T_IRALIGN: Undefined alignment algorithm")
}
# Fill the output image with the unknown value.
call ir_imzero (outim, int (IM_LEN(outim,1)), int (IM_LEN(outim, 2)),
IR_OVAL(ir))
# Shift and match all the subrasters.
call ir_subalign (ir, im, outim, Memc[trimlimits], Memi[IR_IC1(ir)],
Memi[IR_IC2(ir)], Memi[IR_IL1(ir)], Memi[IR_IL2(ir)],
Memi[IR_OC1(ir)], Memi[IR_OC2(ir)], Memi[IR_OL1(ir)],
Memi[IR_OL2(ir)], Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
Memr[IR_DELTAI(ir)], YES, interp, verbose)
# Close up files.
call imunmap (im)
call imunmap (outim)
call dtunmap (dt)
call sfree (sp)
call ir_free (ir)
end
|
