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
|
include <imhdr.h>
define SZ_REGION 4 # Size of region structure
define C1 Memi[$1] # Minimum column
define C2 Memi[$1+1] # Maximum column
define L1 Memi[$1+2] # Minimum line
define L2 Memi[$1+3] # Maximum line
# T_TEXT2MASK -- Create a text file description (FIXPIX) from a mask.
procedure t_mask2text ()
pointer mask # Pixel mask
pointer text # Text file
int i, fd, nc, nl, c1, c2, l, nalloc, nregions
pointer sp, regions, p, pmatch, im, bp
pointer immap(), imgl2s()
int open()
errchk immap, open
begin
call smark (sp)
call salloc (text, SZ_FNAME, TY_CHAR)
call salloc (mask, SZ_FNAME, TY_CHAR)
# Get task parameters.
call clgstr ("mask", Memc[mask], SZ_FNAME)
call clgstr ("text", Memc[text], SZ_FNAME)
# Open the files and abort on an error.
im = immap (Memc[mask], READ_ONLY, 0)
fd = open (Memc[text], NEW_FILE, TEXT_FILE)
nc = IM_LEN(im,1)
nl = IM_LEN(im,2)
nalloc = 0
nregions = 0
do l = 1, nl {
bp = imgl2s (im, l) - 1
for (c1=1; c1<=nc && Mems[bp+c1]==0; c1=c1+1)
;
while (c1 <= nc) {
for (c2=c1+1; c2<=nc && Mems[bp+c2]!=0; c2=c2+1)
;
c2 = c2 - 1
pmatch = NULL
for (i=0; i<nregions; i=i+1) {
p = Memi[regions+i]
if (c1 <= C2(p) && c2 >= C1(p)) {
if (pmatch == NULL) {
L2(p) = l
C1(p) = min (C1(p), c1)
C2(p) = max (C2(p), c2)
pmatch = p
} else {
L1(pmatch) = min (L1(pmatch), L1(p))
C1(pmatch) = min (C1(pmatch), C1(p))
C2(pmatch) = max (C2(pmatch), C2(p))
Memi[regions+i] = Memi[regions+nregions-1]
Memi[regions+nregions-1] = p
nregions = nregions - 1
i = i - 1
}
}
}
if (pmatch == NULL) {
if (nregions == nalloc) {
nalloc = nregions + 1
if (nalloc == 1)
call malloc (regions, nalloc, TY_STRUCT)
else
call realloc (regions, nalloc, TY_STRUCT)
call salloc (Memi[regions+nregions], SZ_REGION,
TY_STRUCT)
}
p = Memi[regions+nregions]
L1(p) = l
L2(p) = l
C1(p) = c1
C2(p) = c2
nregions = nregions + 1
}
for (c1=c2+1; c1<=nc && Mems[bp+c1]==0; c1=c1+1)
;
}
for (i=0; i<nregions; i=i+1) {
p = Memi[regions+i]
if (L2(p) != l) {
call fprintf (fd, "%4d %4d %4d %4d\n")
call pargi (C1(p))
call pargi (C2(p))
call pargi (L1(p))
call pargi (L2(p))
Memi[regions+i] = Memi[regions+nregions-1]
Memi[regions+nregions-1] = p
nregions = nregions - 1
i = i - 1
}
}
}
for (i=0; i<nregions; i=i+1) {
p = Memi[regions+i]
call fprintf (fd, "%4d %4d %4d %4d\n")
call pargi (C1(p))
call pargi (C2(p))
call pargi (L1(p))
call pargi (L2(p))
}
call close (fd)
call imunmap (im)
call mfree (regions, TY_POINTER)
call sfree (sp)
end
|
