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
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
|
/* libpbm1.c - pbm utility library part 1
**
** Copyright (C) 1988 by Jef Poskanzer.
**
** Permission to use, copy, modify, and distribute this software and its
** documentation for any purpose and without fee is hereby granted, provided
** that the above copyright notice appear in all copies and that both that
** copyright notice and this permission notice appear in supporting
** documentation. This software is provided "as is" without express or
** implied warranty.
*/
/* See pmfileio.c for the complicated explanation of this 32/64 bit file
offset stuff.
*/
#define _FILE_OFFSET_BITS 64
#define _LARGE_FILES
#include <stdio.h>
#include "pm_c_util.h"
#include "mallocvar.h"
#include "shhopt.h"
#include "pbm.h"
bit *
pbm_allocrow(unsigned int const cols) {
bit * bitrow;
MALLOCARRAY(bitrow, cols);
if (bitrow == NULL)
pm_error("Unable to allocate space for a %u-column bit row", cols);
return bitrow;
}
void
pbm_init(int * const argcP,
char ** const argv) {
pm_proginit(argcP, (const char **)argv);
}
void
pbm_nextimage(FILE *file, int * const eofP) {
pm_nextimage(file, eofP);
}
void
pbm_check(FILE * const fileP,
enum pm_check_type const checkType,
int const format,
int const cols,
int const rows,
enum pm_check_code * const retvalP) {
if (rows < 0)
pm_error("Invalid number of rows passed to pbm_check(): %d", rows);
if (cols < 0)
pm_error("Invalid number of columns passed to pbm_check(): %d", cols);
if (checkType != PM_CHECK_BASIC) {
if (retvalP)
*retvalP = PM_CHECK_UNKNOWN_TYPE;
} else if (format != RPBM_FORMAT) {
if (retvalP)
*retvalP = PM_CHECK_UNCHECKABLE;
} else {
pm_filepos const bytesPerRow = (cols+7)/8;
pm_filepos const needRasterSize = rows * bytesPerRow;
pm_check(fileP, checkType, needRasterSize, retvalP);
}
}
static unsigned int
bitpop8(unsigned char const x) {
/*----------------------------------------------------------------------------
Return the number of 1 bits in 'x'
-----------------------------------------------------------------------------*/
static unsigned int const p[256] = {
0, 1, 1, 2, 1, 2, 2, 3, 1, 2, 2, 3, 2, 3, 3, 4,
1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,
1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,
2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,
1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,
2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,
2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,
3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,
1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5,
2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,
2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,
3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,
2, 3, 3, 4, 3, 4, 4, 5, 3, 4, 4, 5, 4, 5, 5, 6,
3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,
3, 4, 4, 5, 4, 5, 5, 6, 4, 5, 5, 6, 5, 6, 6, 7,
4, 5, 5, 6, 5, 6, 6, 7, 5, 6, 6, 7, 6, 7, 7, 8 };
return p[x];
}
static int
bitpop(const unsigned char * const packedRow,
unsigned int const cols,
unsigned int const offset) {
/*----------------------------------------------------------------------------
Return the number of 1 bits in 'packedRow', ignoring 0 to 7 bits
at the row start (= on the left edge), indicated by offset.
-----------------------------------------------------------------------------*/
unsigned int const fullLength = cols + offset;
unsigned int sum;
if (fullLength <= 8) {
/* All bits are in a single byte */
sum = bitpop8((packedRow[0] << offset ) & (0xff << (8 - cols)));
} else {
unsigned int const colByteCnt = pbm_packed_bytes(fullLength);
unsigned int const fullByteCnt = fullLength/8;
unsigned int i;
/* First byte, whether it is full or not */
sum = bitpop8(packedRow[0] << offset );
/* Second byte to last full byte */
for (i = 1; i < fullByteCnt; ++i)
sum += bitpop8(packedRow[i]);
/* Partial byte at the right end */
if (colByteCnt > fullByteCnt)
sum += bitpop8(packedRow[i] >> (8 - fullLength%8));
}
return sum;
}
bit
pbm_backgroundbitrow(unsigned const char * const packedBits,
unsigned int const cols,
unsigned int const offset) {
/*----------------------------------------------------------------------------
PBM version of pnm_backgroundxelrow() with additional offset parameter.
When offset == 0, produces the same return value as does
pnm_backgroundxelrow(promoted_bitrow, cols, ...)
-----------------------------------------------------------------------------*/
const unsigned char * const row = &packedBits[offset/8];
unsigned int const rs = offset % 8;
unsigned int const last = pbm_packed_bytes(cols + rs) - 1;
unsigned int retval;
bool const firstbit = (row[0] >> (7-rs)) & 0x01;
bool const lastbit = (row[last] >> (7- (cols+rs-1)%8)) & 0x01;
if (firstbit == lastbit)
retval = firstbit;
else {
if (bitpop(row, cols, rs) >= cols/2)
retval = PBM_BLACK;
else
retval = PBM_WHITE;
}
return retval;
}
|
