Leptonica 1.85.0
Image processing and image analysis suite
Loading...
Searching...
No Matches
dwacomb.2.c
1/*====================================================================*
2 - Copyright (C) 2001 Leptonica. All rights reserved.
3 -
4 - Redistribution and use in source and binary forms, with or without
5 - modification, are permitted provided that the following conditions
6 - are met:
7 - 1. Redistributions of source code must retain the above copyright
8 - notice, this list of conditions and the following disclaimer.
9 - 2. Redistributions in binary form must reproduce the above
10 - copyright notice, this list of conditions and the following
11 - disclaimer in the documentation and/or other materials
12 - provided with the distribution.
13 -
14 - THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
15 - ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
16 - LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
17 - A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL ANY
18 - CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
19 - EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
20 - PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
21 - PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
22 - OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
23 - NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
24 - SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 *====================================================================*/
26
34#ifdef HAVE_CONFIG_H
35#include <config_auto.h>
36#endif /* HAVE_CONFIG_H */
37
38#include <string.h>
39#include "allheaders.h"
40
41PIX *pixMorphDwa_2(PIX *pixd, PIX *pixs, l_int32 operation, char *selname);
42PIX *pixFMorphopGen_2(PIX *pixd, PIX *pixs, l_int32 operation, char *selname);
43l_int32 fmorphopgen_low_2(l_uint32 *datad, l_int32 w,
44 l_int32 h, l_int32 wpld,
45 l_uint32 *datas, l_int32 wpls,
46 l_int32 index);
47
48static l_int32 NUM_SELS_GENERATED = 76;
49static char SEL_NAMES[][80] = {
50 "sel_comb_4h",
51 "sel_comb_4v",
52 "sel_comb_5h",
53 "sel_comb_5v",
54 "sel_comb_6h",
55 "sel_comb_6v",
56 "sel_comb_7h",
57 "sel_comb_7v",
58 "sel_comb_8h",
59 "sel_comb_8v",
60 "sel_comb_9h",
61 "sel_comb_9v",
62 "sel_comb_10h",
63 "sel_comb_10v",
64 "sel_comb_12h",
65 "sel_comb_12v",
66 "sel_comb_14h",
67 "sel_comb_14v",
68 "sel_comb_15h",
69 "sel_comb_15v",
70 "sel_comb_16h",
71 "sel_comb_16v",
72 "sel_comb_18h",
73 "sel_comb_18v",
74 "sel_comb_20h",
75 "sel_comb_20v",
76 "sel_comb_21h",
77 "sel_comb_21v",
78 "sel_comb_22h",
79 "sel_comb_22v",
80 "sel_comb_24h",
81 "sel_comb_24v",
82 "sel_comb_25h",
83 "sel_comb_25v",
84 "sel_comb_27h",
85 "sel_comb_27v",
86 "sel_comb_28h",
87 "sel_comb_28v",
88 "sel_comb_30h",
89 "sel_comb_30v",
90 "sel_comb_32h",
91 "sel_comb_32v",
92 "sel_comb_33h",
93 "sel_comb_33v",
94 "sel_comb_35h",
95 "sel_comb_35v",
96 "sel_comb_36h",
97 "sel_comb_36v",
98 "sel_comb_39h",
99 "sel_comb_39v",
100 "sel_comb_40h",
101 "sel_comb_40v",
102 "sel_comb_42h",
103 "sel_comb_42v",
104 "sel_comb_44h",
105 "sel_comb_44v",
106 "sel_comb_45h",
107 "sel_comb_45v",
108 "sel_comb_48h",
109 "sel_comb_48v",
110 "sel_comb_49h",
111 "sel_comb_49v",
112 "sel_comb_50h",
113 "sel_comb_50v",
114 "sel_comb_52h",
115 "sel_comb_52v",
116 "sel_comb_54h",
117 "sel_comb_54v",
118 "sel_comb_55h",
119 "sel_comb_55v",
120 "sel_comb_56h",
121 "sel_comb_56v",
122 "sel_comb_60h",
123 "sel_comb_60v",
124 "sel_comb_63h",
125 "sel_comb_63v"};
126
146PIX *
147pixMorphDwa_2(PIX *pixd,
148 PIX *pixs,
149 l_int32 operation,
150 char *selname)
151{
152l_int32 bordercolor, bordersize;
153PIX *pixt1, *pixt2, *pixt3;
154
155 if (!pixs)
156 return (PIX *)ERROR_PTR("pixs not defined", __func__, pixd);
157 if (pixGetDepth(pixs) != 1)
158 return (PIX *)ERROR_PTR("pixs must be 1 bpp", __func__, pixd);
159
160 /* Set the border size */
161 bordercolor = getMorphBorderPixelColor(L_MORPH_ERODE, 1);
162 bordersize = 32;
163 if (bordercolor == 0 && operation == L_MORPH_CLOSE)
164 bordersize += 32;
165
166 pixt1 = pixAddBorder(pixs, bordersize, 0);
167 pixt2 = pixFMorphopGen_2(NULL, pixt1, operation, selname);
168 pixt3 = pixRemoveBorder(pixt2, bordersize);
169 pixDestroy(&pixt1);
170 pixDestroy(&pixt2);
171
172 if (!pixd)
173 return pixt3;
174
175 pixCopy(pixd, pixt3);
176 pixDestroy(&pixt3);
177 return pixd;
178}
179
180
204PIX *
205pixFMorphopGen_2(PIX *pixd,
206 PIX *pixs,
207 l_int32 operation,
208 char *selname)
209{
210l_int32 i, index, found, w, h, wpls, wpld, bordercolor, erodeop, borderop;
211l_uint32 *datad, *datas, *datat;
212PIX *pixt;
213
214 if (!pixs)
215 return (PIX *)ERROR_PTR("pixs not defined", __func__, pixd);
216 if (pixGetDepth(pixs) != 1)
217 return (PIX *)ERROR_PTR("pixs must be 1 bpp", __func__, pixd);
218
219 /* Get boundary colors to use */
220 bordercolor = getMorphBorderPixelColor(L_MORPH_ERODE, 1);
221 if (bordercolor == 1)
222 erodeop = PIX_SET;
223 else
224 erodeop = PIX_CLR;
225
226 found = FALSE;
227 for (i = 0; i < NUM_SELS_GENERATED; i++) {
228 if (strcmp(selname, SEL_NAMES[i]) == 0) {
229 found = TRUE;
230 index = 2 * i;
231 break;
232 }
233 }
234 if (found == FALSE)
235 return (PIX *)ERROR_PTR("sel index not found", __func__, pixd);
236
237 if (!pixd) {
238 if ((pixd = pixCreateTemplate(pixs)) == NULL)
239 return (PIX *)ERROR_PTR("pixd not made", __func__, NULL);
240 }
241 else /* for in-place or pre-allocated */
242 pixResizeImageData(pixd, pixs);
243 wpls = pixGetWpl(pixs);
244 wpld = pixGetWpl(pixd);
245
246 /* The images must be surrounded, in advance, with a border of
247 * size 32 pixels (or 64, for closing), that we'll read from.
248 * Fabricate a "proper" image as the subimage within the 32
249 * pixel border, having the following parameters: */
250 w = pixGetWidth(pixs) - 64;
251 h = pixGetHeight(pixs) - 64;
252 datas = pixGetData(pixs) + 32 * wpls + 1;
253 datad = pixGetData(pixd) + 32 * wpld + 1;
254
255 if (operation == L_MORPH_DILATE || operation == L_MORPH_ERODE) {
256 borderop = PIX_CLR;
257 if (operation == L_MORPH_ERODE) {
258 borderop = erodeop;
259 index++;
260 }
261 if (pixd == pixs) { /* in-place; generate a temp image */
262 if ((pixt = pixCopy(NULL, pixs)) == NULL)
263 return (PIX *)ERROR_PTR("pixt not made", __func__, pixd);
264 datat = pixGetData(pixt) + 32 * wpls + 1;
265 pixSetOrClearBorder(pixt, 32, 32, 32, 32, borderop);
266 fmorphopgen_low_2(datad, w, h, wpld, datat, wpls, index);
267 pixDestroy(&pixt);
268 }
269 else { /* not in-place */
270 pixSetOrClearBorder(pixs, 32, 32, 32, 32, borderop);
271 fmorphopgen_low_2(datad, w, h, wpld, datas, wpls, index);
272 }
273 }
274 else { /* opening or closing; generate a temp image */
275 if ((pixt = pixCreateTemplate(pixs)) == NULL)
276 return (PIX *)ERROR_PTR("pixt not made", __func__, pixd);
277 datat = pixGetData(pixt) + 32 * wpls + 1;
278 if (operation == L_MORPH_OPEN) {
279 pixSetOrClearBorder(pixs, 32, 32, 32, 32, erodeop);
280 fmorphopgen_low_2(datat, w, h, wpls, datas, wpls, index+1);
281 pixSetOrClearBorder(pixt, 32, 32, 32, 32, PIX_CLR);
282 fmorphopgen_low_2(datad, w, h, wpld, datat, wpls, index);
283 }
284 else { /* closing */
285 pixSetOrClearBorder(pixs, 32, 32, 32, 32, PIX_CLR);
286 fmorphopgen_low_2(datat, w, h, wpls, datas, wpls, index);
287 pixSetOrClearBorder(pixt, 32, 32, 32, 32, erodeop);
288 fmorphopgen_low_2(datad, w, h, wpld, datat, wpls, index+1);
289 }
290 pixDestroy(&pixt);
291 }
292
293 return pixd;
294}
295
PIX_CLR
#define PIX_CLR
Definition pix.h:447
PIX_SET
#define PIX_SET
Definition pix.h:448