summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat (limited to 'leptonica/prog/binmorph4_reg.c')
-rw-r--r--leptonica/prog/binmorph4_reg.c544
1 files changed, 544 insertions, 0 deletions
diff --git a/leptonica/prog/binmorph4_reg.c b/leptonica/prog/binmorph4_reg.c
new file mode 100644
index 00000000..a9023c56
--- /dev/null
+++ b/leptonica/prog/binmorph4_reg.c
@@ -0,0 +1,544 @@
+/*====================================================================*
+ - Copyright (C) 2001 Leptonica. All rights reserved.
+ -
+ - Redistribution and use in source and binary forms, with or without
+ - modification, are permitted provided that the following conditions
+ - are met:
+ - 1. Redistributions of source code must retain the above copyright
+ - notice, this list of conditions and the following disclaimer.
+ - 2. Redistributions in binary form must reproduce the above
+ - copyright notice, this list of conditions and the following
+ - disclaimer in the documentation and/or other materials
+ - provided with the distribution.
+ -
+ - THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ - ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ - LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ - A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL ANY
+ - CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+ - EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+ - PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
+ - PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
+ - OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
+ - NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+ - SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+ *====================================================================*/
+
+/*
+ * binmorph4_reg.c
+ *
+ * Regression test for dwa brick morph operations
+ * We compare:
+ * (1) morph composite vs. morph non-composite
+ * (2) dwa non-composite vs. morph composite
+ * (3) dwa composite vs. dwa non-composite
+ * (4) dwa composite vs. morph composite
+ * (5) dwa composite vs. morph non-composite
+ * The brick functions all have a pre-allocated pix as the dest.
+ */
+
+#ifdef HAVE_CONFIG_H
+#include <config_auto.h>
+#endif /* HAVE_CONFIG_H */
+
+#include "allheaders.h"
+
+void TestAll(L_REGPARAMS *rp, PIX *pixs, l_int32 symmetric);
+
+l_int32 DoComparisonDwa1(L_REGPARAMS *rp,
+ PIX *pixs, PIX *pix1, PIX *pix2, PIX *pix3,
+ PIX *pix4, PIX *pix5, PIX *pix6, l_int32 isize);
+l_int32 DoComparisonDwa2(L_REGPARAMS *rp,
+ PIX *pixs, PIX *pix1, PIX *pix2, PIX *pix3,
+ PIX *pix4, PIX *pix5, PIX *pix6, l_int32 isize);
+l_int32 DoComparisonDwa3(L_REGPARAMS *rp,
+ PIX *pixs, PIX *pix1, PIX *pix2, PIX *pix3,
+ PIX *pix4, PIX *pix5, PIX *pix6, l_int32 isize);
+l_int32 DoComparisonDwa4(L_REGPARAMS *rp,
+ PIX *pixs, PIX *pix1, PIX *pix2, PIX *pix3,
+ PIX *pix4, PIX *pix5, PIX *pix6, l_int32 isize);
+l_int32 DoComparisonDwa5(L_REGPARAMS *rp,
+ PIX *pixs, PIX *pix1, PIX *pix2, PIX *pix3,
+ PIX *pix4, PIX *pix5, PIX *pix6, l_int32 isize);
+void PixCompareDwa(L_REGPARAMS *rp,
+ l_int32 size, const char *type, PIX *pix1, PIX *pix2,
+ PIX *pix3, PIX *pix4, PIX *pix5, PIX *pix6);
+
+#define TIMING 0
+
+ /* Note: the symmetric case requires an extra border of size
+ * approximately 40 to succeed for all SE up to size 64. With
+ * a smaller border the differences are small, and most of the
+ * problem seems to be in the non-dwa code, because we are doing
+ * sequential erosions without an extra border, and things aren't
+ * being properly initialized. To avoid these errors, add the border
+ * in advance for symmetric b.c.
+ * Note that asymmetric b.c. are recommended for document image
+ * operations, and this test passes for asymmetric b.c. without
+ * any added border. */
+
+int main(int argc,
+ char **argv)
+{
+PIX *pixs;
+L_REGPARAMS *rp;
+
+ if (regTestSetup(argc, argv, &rp))
+ return 1;
+
+ pixs = pixRead("feyn-fract.tif");
+ TestAll(rp, pixs, FALSE);
+ TestAll(rp, pixs, TRUE);
+ pixDestroy(&pixs);
+ return regTestCleanup(rp);
+}
+
+void
+TestAll(L_REGPARAMS *rp,
+ PIX *pixs,
+ l_int32 symmetric)
+{
+l_int32 i;
+PIX *pix1, *pix2, *pix3, *pix4, *pix5, *pix6;
+
+ if (symmetric) {
+ /* This works properly with an added border of 40 */
+ resetMorphBoundaryCondition(SYMMETRIC_MORPH_BC);
+ pix1 = pixAddBorder(pixs, 40, 0);
+ pixTransferAllData(pixs, &pix1, 0, 0);
+ lept_stderr("Testing with symmetric boundary conditions\n");
+ } else {
+ resetMorphBoundaryCondition(ASYMMETRIC_MORPH_BC);
+ lept_stderr("Testing with asymmetric boundary conditions\n");
+ }
+
+ pix1 = pixCreateTemplateNoInit(pixs);
+ pix2 = pixCreateTemplateNoInit(pixs);
+ pix3 = pixCreateTemplateNoInit(pixs);
+ pix4 = pixCreateTemplateNoInit(pixs);
+ pix5 = pixCreateTemplateNoInit(pixs);
+ pix6 = pixCreateTemplateNoInit(pixs);
+
+ for (i = 2; i < 64; i++) {
+
+#if 1
+ /* Compare morph composite with morph non-composite */
+ DoComparisonDwa1(rp, pixs, pix1, pix2, pix3, pix4, pix5, pix6, i);
+#endif
+
+#if 1
+ /* Compare DWA non-composite with morph composite */
+ if (i < 16)
+ DoComparisonDwa2(rp, pixs, pix1, pix2, pix3, pix4, pix5, pix6, i);
+ /* Compare DWA composite with DWA non-composite */
+ if (i < 16)
+ DoComparisonDwa3(rp, pixs, pix1, pix2, pix3, pix4, pix5, pix6, i);
+ /* Compare DWA composite with morph composite */
+ DoComparisonDwa4(rp, pixs, pix1, pix2, pix3, pix4, pix5, pix6, i);
+ /* Compare DWA composite with morph non-composite */
+ DoComparisonDwa5(rp, pixs, pix1, pix2, pix3, pix4, pix5, pix6, i);
+#endif
+ }
+ lept_stderr("\n");
+ pixDestroy(&pix1);
+ pixDestroy(&pix2);
+ pixDestroy(&pix3);
+ pixDestroy(&pix4);
+ pixDestroy(&pix5);
+ pixDestroy(&pix6);
+}
+
+
+ /* Morph composite with morph non-composite */
+l_int32
+DoComparisonDwa1(L_REGPARAMS *rp,
+ PIX *pixs,
+ PIX *pix1,
+ PIX *pix2,
+ PIX *pix3,
+ PIX *pix4,
+ PIX *pix5,
+ PIX *pix6,
+ l_int32 isize)
+{
+l_int32 fact1, fact2, size;
+
+ selectComposableSizes(isize, &fact1, &fact2);
+ size = fact1 * fact2;
+
+ lept_stderr("..%d..", size);
+
+ if (TIMING) startTimer();
+ pixDilateCompBrick(pix1, pixs, size, 1);
+ pixDilateCompBrick(pix3, pixs, 1, size);
+ pixDilateCompBrick(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixDilateBrick(pix2, pixs, size, 1);
+ pixDilateBrick(pix4, pixs, 1, size);
+ pixDilateBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "dilate", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixErodeCompBrick(pix1, pixs, size, 1);
+ pixErodeCompBrick(pix3, pixs, 1, size);
+ pixErodeCompBrick(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixErodeBrick(pix2, pixs, size, 1);
+ pixErodeBrick(pix4, pixs, 1, size);
+ pixErodeBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "erode", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixOpenCompBrick(pix1, pixs, size, 1);
+ pixOpenCompBrick(pix3, pixs, 1, size);
+ pixOpenCompBrick(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixOpenBrick(pix2, pixs, size, 1);
+ pixOpenBrick(pix4, pixs, 1, size);
+ pixOpenBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "open", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixCloseSafeCompBrick(pix1, pixs, size, 1);
+ pixCloseSafeCompBrick(pix3, pixs, 1, size);
+ pixCloseSafeCompBrick(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixCloseSafeBrick(pix2, pixs, size, 1);
+ pixCloseSafeBrick(pix4, pixs, 1, size);
+ pixCloseSafeBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "close", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ return 0;
+}
+
+
+ /* Dwa non-composite with morph composite */
+l_int32
+DoComparisonDwa2(L_REGPARAMS *rp,
+ PIX *pixs,
+ PIX *pix1,
+ PIX *pix2,
+ PIX *pix3,
+ PIX *pix4,
+ PIX *pix5,
+ PIX *pix6,
+ l_int32 isize)
+{
+l_int32 fact1, fact2, size;
+
+ selectComposableSizes(isize, &fact1, &fact2);
+ size = fact1 * fact2;
+
+ lept_stderr("..%d..", size);
+
+ if (TIMING) startTimer();
+ pixDilateBrickDwa(pix1, pixs, size, 1);
+ pixDilateBrickDwa(pix3, pixs, 1, size);
+ pixDilateBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixDilateCompBrick(pix2, pixs, size, 1);
+ pixDilateCompBrick(pix4, pixs, 1, size);
+ pixDilateCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "dilate", pix1, pix2, pix3, pix4, pix5, pix6);
+
+/* pixDisplay(pix1, 100, 100); */
+/* pixDisplay(pix2, 800, 100); */
+
+ if (TIMING) startTimer();
+ pixErodeBrickDwa(pix1, pixs, size, 1);
+ pixErodeBrickDwa(pix3, pixs, 1, size);
+ pixErodeBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixErodeCompBrick(pix2, pixs, size, 1);
+ pixErodeCompBrick(pix4, pixs, 1, size);
+ pixErodeCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "erode", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixOpenBrickDwa(pix1, pixs, size, 1);
+ pixOpenBrickDwa(pix3, pixs, 1, size);
+ pixOpenBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixOpenCompBrick(pix2, pixs, size, 1);
+ pixOpenCompBrick(pix4, pixs, 1, size);
+ pixOpenCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "open", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixCloseBrickDwa(pix1, pixs, size, 1);
+ pixCloseBrickDwa(pix3, pixs, 1, size);
+ pixCloseBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixCloseSafeCompBrick(pix2, pixs, size, 1);
+ pixCloseSafeCompBrick(pix4, pixs, 1, size);
+ pixCloseSafeCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "close", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ return 0;
+}
+
+
+ /* Dwa composite with dwa non-composite */
+l_int32
+DoComparisonDwa3(L_REGPARAMS *rp,
+ PIX *pixs,
+ PIX *pix1,
+ PIX *pix2,
+ PIX *pix3,
+ PIX *pix4,
+ PIX *pix5,
+ PIX *pix6,
+ l_int32 isize)
+{
+l_int32 fact1, fact2, size;
+
+ selectComposableSizes(isize, &fact1, &fact2);
+ size = fact1 * fact2;
+
+ lept_stderr("..%d..", size);
+
+ if (TIMING) startTimer();
+ pixDilateCompBrickDwa(pix1, pixs, size, 1);
+ pixDilateCompBrickDwa(pix3, pixs, 1, size);
+ pixDilateCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixDilateBrickDwa(pix2, pixs, size, 1);
+ pixDilateBrickDwa(pix4, pixs, 1, size);
+ pixDilateBrickDwa(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "dilate", pix1, pix2, pix3, pix4, pix5, pix6);
+
+/* pixDisplay(pix1, 100, 100); */
+/* pixDisplay(pix2, 800, 100); */
+
+ if (TIMING) startTimer();
+ pixErodeCompBrickDwa(pix1, pixs, size, 1);
+ pixErodeCompBrickDwa(pix3, pixs, 1, size);
+ pixErodeCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixErodeBrickDwa(pix2, pixs, size, 1);
+ pixErodeBrickDwa(pix4, pixs, 1, size);
+ pixErodeBrickDwa(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "erode", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixOpenCompBrickDwa(pix1, pixs, size, 1);
+ pixOpenCompBrickDwa(pix3, pixs, 1, size);
+ pixOpenCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixOpenBrickDwa(pix2, pixs, size, 1);
+ pixOpenBrickDwa(pix4, pixs, 1, size);
+ pixOpenBrickDwa(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "open", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixCloseCompBrickDwa(pix1, pixs, size, 1);
+ pixCloseCompBrickDwa(pix3, pixs, 1, size);
+ pixCloseCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixCloseBrickDwa(pix2, pixs, size, 1);
+ pixCloseBrickDwa(pix4, pixs, 1, size);
+ pixCloseBrickDwa(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "close", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ return 0;
+}
+
+
+ /* Dwa composite with morph composite */
+l_int32
+DoComparisonDwa4(L_REGPARAMS *rp,
+ PIX *pixs,
+ PIX *pix1,
+ PIX *pix2,
+ PIX *pix3,
+ PIX *pix4,
+ PIX *pix5,
+ PIX *pix6,
+ l_int32 isize)
+{
+l_int32 fact1, fact2, size;
+
+ selectComposableSizes(isize, &fact1, &fact2);
+ size = fact1 * fact2;
+
+ lept_stderr("..%d..", size);
+
+ if (TIMING) startTimer();
+ pixDilateCompBrickDwa(pix1, pixs, size, 1);
+ pixDilateCompBrickDwa(pix3, pixs, 1, size);
+ pixDilateCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixDilateCompBrick(pix2, pixs, size, 1);
+ pixDilateCompBrick(pix4, pixs, 1, size);
+ pixDilateCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "dilate", pix1, pix2, pix3, pix4, pix5, pix6);
+
+/* pixDisplay(pix1, 100, 100); */
+/* pixDisplay(pix2, 800, 100); */
+
+ if (TIMING) startTimer();
+ pixErodeCompBrickDwa(pix1, pixs, size, 1);
+ pixErodeCompBrickDwa(pix3, pixs, 1, size);
+ pixErodeCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixErodeCompBrick(pix2, pixs, size, 1);
+ pixErodeCompBrick(pix4, pixs, 1, size);
+ pixErodeCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "erode", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixOpenCompBrickDwa(pix1, pixs, size, 1);
+ pixOpenCompBrickDwa(pix3, pixs, 1, size);
+ pixOpenCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixOpenCompBrick(pix2, pixs, size, 1);
+ pixOpenCompBrick(pix4, pixs, 1, size);
+ pixOpenCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "open", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixCloseCompBrickDwa(pix1, pixs, size, 1);
+ pixCloseCompBrickDwa(pix3, pixs, 1, size);
+ pixCloseCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixCloseSafeCompBrick(pix2, pixs, size, 1);
+ pixCloseSafeCompBrick(pix4, pixs, 1, size);
+ pixCloseSafeCompBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "close", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ return 0;
+}
+
+ /* Dwa composite with morph non-composite */
+l_int32
+DoComparisonDwa5(L_REGPARAMS *rp,
+ PIX *pixs,
+ PIX *pix1,
+ PIX *pix2,
+ PIX *pix3,
+ PIX *pix4,
+ PIX *pix5,
+ PIX *pix6,
+ l_int32 isize)
+{
+l_int32 fact1, fact2, size;
+
+ selectComposableSizes(isize, &fact1, &fact2);
+ size = fact1 * fact2;
+
+ lept_stderr("..%d..", size);
+
+ if (TIMING) startTimer();
+ pixDilateCompBrickDwa(pix1, pixs, size, 1);
+ pixDilateCompBrickDwa(pix3, pixs, 1, size);
+ pixDilateCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixDilateBrick(pix2, pixs, size, 1);
+ pixDilateBrick(pix4, pixs, 1, size);
+ pixDilateBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "dilate", pix1, pix2, pix3, pix4, pix5, pix6);
+
+/* pixDisplay(pix1, 100, 100); */
+/* pixDisplay(pix2, 800, 100); */
+
+ if (TIMING) startTimer();
+ pixErodeCompBrickDwa(pix1, pixs, size, 1);
+ pixErodeCompBrickDwa(pix3, pixs, 1, size);
+ pixErodeCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixErodeBrick(pix2, pixs, size, 1);
+ pixErodeBrick(pix4, pixs, 1, size);
+ pixErodeBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "erode", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixOpenCompBrickDwa(pix1, pixs, size, 1);
+ pixOpenCompBrickDwa(pix3, pixs, 1, size);
+ pixOpenCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixOpenBrick(pix2, pixs, size, 1);
+ pixOpenBrick(pix4, pixs, 1, size);
+ pixOpenBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "open", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ if (TIMING) startTimer();
+ pixCloseCompBrickDwa(pix1, pixs, size, 1);
+ pixCloseCompBrickDwa(pix3, pixs, 1, size);
+ pixCloseCompBrickDwa(pix5, pixs, size, size);
+ if (TIMING) lept_stderr("Time Dwa: %7.3f sec\n", stopTimer());
+ if (TIMING) startTimer();
+ pixCloseSafeBrick(pix2, pixs, size, 1);
+ pixCloseSafeBrick(pix4, pixs, 1, size);
+ pixCloseSafeBrick(pix6, pixs, size, size);
+ if (TIMING) lept_stderr("Time Rop: %7.3f sec\n", stopTimer());
+ PixCompareDwa(rp, size, "close", pix1, pix2, pix3, pix4, pix5, pix6);
+
+ return 0;
+}
+
+
+void
+PixCompareDwa(L_REGPARAMS *rp,
+ l_int32 size,
+ const char *type,
+ PIX *pix1,
+ PIX *pix2,
+ PIX *pix3,
+ PIX *pix4,
+ PIX *pix5,
+ PIX *pix6)
+{
+l_int32 same;
+
+ pixEqual(pix1, pix2, &same);
+ regTestCompareValues(rp, TRUE, same, 0);
+ if (!same)
+ lept_stderr("%s (%d, 1) not same\n", type, size);
+ pixEqual(pix3, pix4, &same);
+ regTestCompareValues(rp, TRUE, same, 0);
+ if (!same)
+ lept_stderr("%s (1, %d) not same\n", type, size);
+ pixEqual(pix5, pix6, &same);
+ regTestCompareValues(rp, TRUE, same, 0);
+ if (!same)
+ lept_stderr("%s (%d, %d) not same\n", type, size, size);
+}
+