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
|
include <imio.h>
include <error.h>
include "hdiff.h"
# HDIFF -- Show the differences between two header files
#
# B.Simon 23-Apr-90 Original
procedure hdiff ()
#--
pointer input1 # First image template
pointer input2 # Second image template
pointer keywords # List of keywords to compare
bool nofile # Ignore error if files do not exist
bool negate, bjunk
int nlist, fd, ic, nc, ilist, junk
int ngroup1, ngroup2, igroup, len1, len2, ndx1, ndx2
pointer sp, cluster1, cluster2, image1, image2, keynames
pointer im1, im2, tp1, tp2, buffer1, buffer2, index1, index2
string badimcount "The number of images in both lists are not the same"
string badgcount "The number of groups in both iamges are not the same"
bool clgetb(), tp_fetch()
int imtlen(), open(), getline(), imtgetim(), word_count()
pointer imtopenp(), tp_open(), immap()
begin
# Allocate dynamic memory for strings
call smark (sp)
call salloc (cluster1, SZ_FNAME, TY_CHAR)
call salloc (cluster2, SZ_FNAME, TY_CHAR)
call salloc (image1, SZ_FNAME, TY_CHAR)
call salloc (image2, SZ_FNAME, TY_CHAR)
call salloc (keywords, SZ_FNAME, TY_CHAR)
call salloc (keynames, SZ_COMMAND, TY_CHAR)
# Read the task parameters
input1 = imtopenp ("input1")
input2 = imtopenp ("input2")
call clgstr ("keywords", Memc[keywords], SZ_FNAME)
nofile = clgetb ("nofile")
# Make sure that the number of images is the same in both lists
nlist = imtlen (input1)
if (imtlen (input2) != nlist)
call error (1, badimcount)
# Create a list of keyword names
if (Memc[keywords] != '@') {
call strcpy (Memc[keywords], Memc[keynames], SZ_COMMAND)
} else {
fd = open (Memc[keywords+1], READ_ONLY, TEXT_FILE)
for (ic = 0; ic < SZ_COMMAND; ic = ic + nc) {
nc = getline (fd, Memc[keynames+ic])
if (nc <= 0)
break
}
Memc[keynames+ic] = EOS
call close (fd)
}
call strupr (Memc[keynames])
if (Memc[keynames] == '~') {
negate = true
keynames = keynames + 1
} else if (word_count (Memc[keynames]) == 0) {
negate = true
Memc[keynames] = EOS
} else {
negate = false
}
# Compare each image
do ilist = 1, nlist {
junk = imtgetim (input1, Memc[cluster1], SZ_FNAME)
junk = imtgetim (input2, Memc[cluster2], SZ_FNAME)
# Compare each group in each image
iferr {
tp1 = tp_open (Memc[cluster1], 0, ngroup1)
tp2 = tp_open (Memc[cluster2], 0, ngroup2)
} then if (nofile) {
next
} else {
call erract (EA_ERROR)
}
if (ngroup1 != ngroup2)
call error (1, badgcount)
do igroup = 1, ngroup1 {
bjunk = tp_fetch (tp1, Memc[image1])
bjunk = tp_fetch (tp2, Memc[image2])
im1 = immap (Memc[image1], READ_ONLY, 0)
im2 = immap (Memc[image2], READ_ONLY, 0)
len1 = IM_LENHDRMEM(im1) * SZ_STRUCT
len2 = IM_LENHDRMEM(im2) * SZ_STRUCT
ndx1 = len1 / SZ_RECORD + 1
ndx2 = len2 / SZ_RECORD + 1
call malloc (index1, ndx1, TY_INT)
call malloc (index2, ndx2, TY_INT)
call malloc (buffer1, len1, TY_CHAR)
call malloc (buffer2, len2, TY_CHAR)
call getkeybuf (im1, igroup, negate, Memc[keynames],
Memc[buffer1], Memi[index1], ndx1)
call getkeybuf (im2, igroup, negate, Memc[keynames],
Memc[buffer2], Memi[index2], ndx2)
call pr_hdiff (im1, im2, Memc[buffer1], Memi[index1], ndx1,
Memc[buffer2], Memi[index2], ndx2)
call mfree (index1, TY_INT)
call mfree (index2, TY_INT)
call mfree (buffer1, TY_CHAR)
call mfree (buffer2, TY_CHAR)
call imunmap (im1)
call imunmap (im2)
}
call tp_close (tp1)
call tp_close (tp2)
}
call sfree (sp)
end
|
