gocr-ruby 0.0.1

data/ext/gocr/output.h

@@ -0,0 +1,37 @@
+/*
+This is a Optical-Character-Recognition program
+Copyright (C) 2000  Joerg Schulenburg
+
+This program is free software; you can redistribute it and/or
+modify it under the terms of the GNU General Public License
+as published by the Free Software Foundation; either version 2
+of the License, or (at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program; if not, write to the Free Software
+Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
+
+ see README for EMAIL-address */
+ 
+#ifndef OUTPUT_H
+#define OUTPUT_H
+
+#include <stdlib.h>
+#include <stdio.h>
+#include "pnm.h"
+#include "gocr.h"
+#include "list.h"
+
+void out_b(struct box *px, pix *b, int x0, int y0, int dx, int dy, int cs );
+void out_x(struct box *px);
+void out_x2(struct box *box1,struct box *box2);
+int output_list(job_t *job);
+int debug_img(char *fname, struct job_s *job, int opt);
+
+
+#endif

data/ext/gocr/pcx.c

@@ -0,0 +1,153 @@
+/*
+This is a Optical-Character-Recognition program
+Copyright (C) 1999  Joerg Schulenburg
+
+This program is free software; you can redistribute it and/or
+modify it under the terms of the GNU General Public License
+as published by the Free Software Foundation; either version 2
+of the License, or (at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program; if not, write to the Free Software
+Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
+
+ see README for EMAIL-address
+ */
+/* plan: use popen("ppm2pcx -packed ...","w"); for writing pcx */
+
+#include <stdio.h>
+#include <stdlib.h>
+/* #include <assert.h> */
+
+#include "pcx.h"
+
+typedef unsigned char byte;
+
+#define ERR(x) { fprintf(stderr,"ERROR "__FILE__" L%d: " x "\n",__LINE__);exit(1);}
+
+int err;
+/* --- needed for reading PCX-files */
+unsigned char read_b(FILE *f1){
+  unsigned char c=0; c=fgetc(f1); if(feof(f1) || ferror(f1))err=1; return c;
+}
+
+/* something here is wrong! */
+void readpcx(char *name,pix *p,int vvv){  /* see pcx.format.txt */
+  int page,pages,nx,ny,i,j,b,x,y,bpl,bits,pal[256][3];
+  FILE *f1;
+  unsigned char *pic,h[128],bb,b1,b2,b3;
+  err=0;
+  for(i=0;i<256;i++)for(j=0;j<3;j++)pal[i][j]=i;
+  f1=fopen(name,"rb"); if(!f1) ERR("open");
+  if(fread(h,1,128,f1)!=128)ERR("read PCX header");    /* 128 Byte lesen -> h[] */
+  if(h[0]!=10)ERR("no ZSoft sign");	/* ZSoft sign */
+  if(h[2]>  1)ERR("unknown coding");	/* run length encoding */
+  bits = h[3];		/* 1 or 8 */
+  if(bits!=1 && bits!=8)ERR("only 1 or 8 bits supported");
+  nx = h[ 9]*256+h[ 8] - h[ 5]*256-h[ 4] +1;	/* Xmax-Xmin */
+  ny = h[11]*256+h[10] - h[ 7]*256-h[ 6] +1;	/* Ymax-Ymin */
+  pages=h[65]; bpl=h[66]+256*h[67]; /* bytes per line */
+  if(vvv)
+  fprintf(stderr,"# PCX version=%d bits=%d x=%d y=%d HRes=%d VRes=%d\n"
+                 "# NPlanes=%d BytesPerLine=%d Palette=%s",
+    h[1],bits,nx,ny,h[12]+256*h[13],h[14]+256*h[15],
+    pages,bpl,((h[68]==1)?"1=color/bw":"2=gray"));
+  /* line1(NP=4): RRRRR...,GGGG....,BBBBB...,IIII...., line2: RRRR...,GGGG.... */
+  /* C4 EF = (C4&3F)*EF = EF EF EF EF */
+  fflush(stdout);
+  /* palette: for(i=0;i<16;i++) for(j=0;j<3;j++) h[16+3*i+j]  */
+  if(pages>1)for(b=0;b<16;b++) for(i=0;i<16;i++)
+             for(j=0;j< 3;j++) pal[b*16+i][j]=h[16+3*i+j]>>2;
+  if(bits>7){
+   fseek(f1,-3*256,2); if(fread(pal,3,256,f1)!=256)ERR("read palette");
+   for(i=0;i<256;i++) for(j=0;j<3;j++) pal[i][j]>>=2;
+  }
+  fseek(f1,128,0);
+  pic=(unsigned char *)malloc( nx*ny );
+  if(pic==NULL)ERR("no memory");	/* no memory */
+  x=y=0;
+  do {
+    for(page=0;page<pages;page++)    /* 192 == 0xc0 => b1=counter */
+    do {
+      b1=1; bb=read_b(f1); b2=bb; if(b1==192)fprintf(stderr,"?");
+      if((b2>=192) && (h[2]==1)){b1=b2&63;bb=read_b(f1);b2=bb;}
+      if(err){fprintf(stderr,"\nread error x=%d y=%d\n",x,y);x=nx;y=ny;break;}
+      for(b3=0;b3<b1;b3++)for(b=0;b<8;b+=bits,x++)if(x<nx){
+        bb=(b2>>(8-bits-b)) & ~((~0)<<bits);
+        if(bits==1 && bb==1) bb=240;
+        if(page==0) pic[x+nx*y] =(byte)bb;
+        else        pic[x+nx*y]|=(byte)bb<<(page*bits);
+      }
+    } while(x<(9-bits)*bpl);  x=0; y++;
+  } while(y<ny);
+  /*  */
+  fclose(f1);
+  p->p=pic;  p->x=nx;  p->y=ny; p->bpp=1;
+  if(vvv)fprintf(stderr,"\n");
+}
+
+/* -----------------------------------------------------------------------
+// write bmp 8bit palette no RLE
+//   bit 2+3 used for color coding (markers)
+//   replaced by writeppm (ppm.gz) and is obsolate now, removed later
+ */
+void writebmp(char *name,pix p,int vvv){ /* see pcx.format.txt */
+  int nx,ny,i,y,rest[4]={0,0,0,0};
+  FILE *f1;
+  /*FIXME jb static*/static unsigned char *pic, h[54+4*256];
+  long fs,fo,hs,is; /* filesize, offset, headersize, imagesize */
+
+  nx=p.x; ny=p.y; pic=p.p;
+  if (nx&3) nx+=4-(nx&3);  /* must be mod4 ? */
+  hs=40;	    /* bmi headersize fix */
+  is=nx*ny; /* imagesize */
+  fo=14+hs+4*256;
+  fs=fo+is;
+  for(i=0;i<54;i++){ h[i]=0; }
+  /* BITMAPFILEHEADER */
+  h[ 0]='B';          h[ 1]='M'; 		/* type of file BMP */
+  h[ 2]= fs     &255; h[ 3]=(fs>> 8)&255;
+  h[ 4]=(fs>>16)&255; h[ 5]=(fs>>24)&255;	/* size of file */
+  h[10]= fo     &255; h[11]=(fo>> 8)&255;
+  h[12]=(fo>>16)&255; h[13]=(fo>>24)&255;	/* offset to image data */
+  /* BITMAPINFO  (BITMAPCOREHEADER not used here) */
+  /* 14 - HEADER */
+  h[14]= hs     &255; h[15]=(hs>> 8)&255;
+  h[16]=(hs>>16)&255; h[17]=(hs>>24)&255;	/* bmi-header size */
+  h[18]= nx     &255; h[19]=(nx>> 8)&255;
+  h[20]=(0l>>16)&255; h[21]=(0l>>24)&255;	/* WIDTH/pixel */
+  h[22]= ny     &255; h[23]=(ny>> 8)&255;
+  h[24]=(0l>>16)&255; h[25]=(0l>>24)&255;	/* HIGH/pixel */
+  h[26]=1;            				/* planes */
+  h[28]=8;            				/* bits/pixel 1,4,8,24 */
+  h[30]=0;            				/* compression */
+  h[34]= is     &255; h[35]=(is>> 8)&255;
+  h[36]=(is>>16)&255; h[37]=(is>>24)&255;	/* sizeImage (can be 0 if ~RLE) */
+  h[38]=0;h[39]=1;	/* ca 100dpi, x/meter */
+  h[42]=0;h[43]=1;	/* 	      y/meter */
+  h[46]=0;h[47]=1;	/* colorused (0=maximum)  */
+  h[50]=0;h[51]=1;	/* colorimportand (0=all) */
+  /* 54 - endofheader */
+  for(i=0;i<256;i++){ 
+     h[54+4*i+0]=((~((i & 2)*64)) &                    (i & (128+64)))|63;
+     h[54+4*i+1]=((~((i & 2)*64)) &  (~((i & 4)*32)) & (i & (128+64)))|63;
+     h[54+4*i+2]=(  ((i & 2)* 8)  | ((~((i & 4)*32)) & (i & (128+64)))|63);
+  }			 /* blue-green-red */
+  f1=fopen(name,"wb"); if(!f1) fprintf(stderr," error opening file\n");
+  if(!f1)ERR("open");		/* open-error */
+  if(fwrite(h,1,54+4*256,f1)!=54+4*256)ERR("write head"); 
+  if(vvv) fprintf(stderr,"# write BMP x=%d y=%d\n",nx,ny);
+  for(y=ny-1;y>=0;y--){
+    if(((int)fwrite(pic+p.x*y,1,p.x,f1))!=p.x)ERR("write");
+    if(nx>p.x)
+    if(((int)fwrite(rest,1,nx-p.x,f1))!=nx-p.x)ERR("write");
+  }
+  fclose(f1);
+}
+
+/* ---------------------------------------------------------------------- */

data/ext/gocr/pcx.h

@@ -0,0 +1,9 @@
+
+#include "pnm.h"
+
+void readpcx(char *name,pix *p,int vvv);
+
+/* write 8bit palette no RLE, ToDo: obsolete?  */
+void writebmp(char *name,pix p,int vvv);
+
+/* ------------------------------------------------------------------------ */

data/ext/gocr/pgm2asc.c

@@ -0,0 +1,3259 @@
+/*
+This is a Optical-Character-Recognition program
+Copyright (C) 2000-2012  Joerg Schulenburg
+
+This program is free software; you can redistribute it and/or
+modify it under the terms of the GNU General Public License
+as published by the Free Software Foundation; either version 2
+of the License, or (at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program; if not, write to the Free Software
+Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
+
+ see README for EMAIL-address
+
+  sometimes I have written comments in german language, sorry for that
+
+ - look for ??? for preliminary code
+ - space: avX=22 11-13 (empirical estimated)
+          avX=16  5-7
+          avX= 7  5-6
+         
+ ToDo: - add filter (r/s mismatch) g300c1
+       - better get_line2 function (problems on high resolution)
+       - write parallelizable code!
+       - learnmode (optimize filter)
+       - use ispell for final control or if unsure
+       - better line scanning (if not even)
+       - step 5: same chars differ? => expert mode
+       - chars dx>dy and above 50% hor-crossing > 4 is char-group ?
+       - detect color of chars and background
+       - better word space calculation (look at the examples)
+          (distance: left-left, middle-middle, left-right, thickness of e *0.75)
+
+   GLOBAL DATA (mostly structures)
+   - pix   : image - one byte per pixel  bits0-2=working
+   - lines : rows of the text (points to pix)
+   - box   : list of bounding box for character 
+   - obj   : objects (lines, splines, etc. building a character)
+ */
+
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <assert.h>
+#include <string.h>
+#include <ctype.h>
+#include "config.h"
+#ifdef HAVE_WCHAR_H
+#include <wchar.h>
+#endif
+
+#include "amiga.h"
+#include "list.h"
+#include "pgm2asc.h"
+// #include "pcx.h"        /* needed for writebmp (removed later) */
+/* ocr1 is the test-engine - remember: this is development version */
+#include "ocr1.h"
+/* first engine */
+#include "ocr0.h"
+#include "otsu.h"
+#include "barcode.h"
+#include "progress.h"
+#include "unicode_defs.h" /* UNKNOWN + PICTURES + ... */
+
+#include "gocr.h"
+
+/* wew: will be exceeded by capitals at 1200dpi */
+#define MaxBox (100*200)	// largest possible letter (buffersize)
+#define MAX(a,b)			((a) >= (b) ? (a) : (b))
+
+/* if the system does not know about wchar.h, define functions here */
+#ifndef HAVE_WCHAR_H
+/* typedef unsigned wchar_t; */
+/* Find the first occurrence of WC in WCS.  */
+wchar_t *wcschr (wchar_t *wcs, wchar_t wc) {
+  int i; for(i=0;wcs[i];i++) if (wcs[i]==wc) return wcs+i; return NULL;
+}
+wchar_t *wcscpy (wchar_t *dest, const wchar_t *src) {
+  int i; for(i=0;src[i];i++) dest[i]=src[i]; dest[i]=0; return dest;
+}
+size_t wcslen (const wchar_t *s){
+  size_t i; for(i=0;s[i];i++); return i;
+}
+#endif
+#ifndef HAVE_WCSDUP
+wchar_t * wcsdup (const wchar_t *WS) {	/* its a gnu extension */
+  wchar_t *copy;
+  copy = (wchar_t *) malloc((wcslen(WS)+1)*sizeof(wchar_t));
+  if (!copy)return NULL;
+  wcscpy(copy, WS);
+  return copy;
+}   
+#endif
+
+// ------------------------ feature extraction -----------------
+// -------------------------------------------------------------
+// detect maximas in of line overlaps (return in %) and line coordinates
+// this is for future use
+#define HOR 1    // horizontal
+#define VER 2    // vertical
+#define RIS 3    // rising=steigend
+#define FAL 4    // falling=fallend
+
+/* exchange two variables */
+static void swap(int *a, int *b) { 
+  int c = *a;
+  *a = *b;
+  *b = c;
+}
+
+// calculate the overlapping of the line (0-1) with black points 
+// by recursive bisection 
+// line: y=dy/dx*x+b, implicit form: d=F(x,y)=dy*x-dx*y+b*dx=0 
+// incremental y(i+1)=m*(x(i)+1)+b, F(x+1,y+1)=f(F(x,y))
+// ret & 1 => inverse pixel!
+// d=2*F(x,y) integer numbers
+int get_line(int x0, int y0, int x1, int y1, pix *p, int cs, int ret){
+   int dx,dy,incrE,incrNE,d,x,y,r0,r1,ty,tx,
+       *px,*py,*pdx,*pdy,*ptx,*pty,*px1;
+   dx=abs(x1-x0); tx=((x1>x0)?1:-1);    // tx=x-spiegelung (new)  
+   dy=abs(y1-y0); ty=((y1>y0)?1:-1);	// ty=y-spiegelung (new)
+   // rotate coordinate system if dy>dx
+/*bbg: can be faster if instead of pointers we use the variables and swaps? */
+/*js: Do not know, I am happy that the current code is working and is small */
+   if(dx>dy){ pdx=&dx;pdy=&dy;px=&x;py=&y;ptx=&tx;pty=&ty;px1=&x1; }
+   else     { pdx=&dy;pdy=&dx;px=&y;py=&x;ptx=&ty;pty=&tx;px1=&y1; }
+   if( *ptx<0 ){ swap(&x0,&x1);swap(&y0,&y1);tx=-tx;ty=-ty; }
+   d=((*pdy)<<1)-(*pdx); incrE=(*pdy)<<1; incrNE=((*pdy)-(*pdx))<<1;  
+   x=x0; y=y0; r0=r1=0; /* dd=tolerance (store max drift) */
+   while( (*px)<=(*px1) ){ 
+     if( ((getpixel(p,x,y)<cs)?1:0)^(ret&1) ) r0++; else r1++;
+     (*px)++; if( d<=0 ){ d+=incrE; } else { d+=incrNE; (*py)+=(*pty); }
+   }
+   return (r0*(ret&~1))/(r0+r1); // ret==100 => percentage %
+}
+
+// this function should detect whether a direct connection between points
+//   exists or not, not finally implemented
+// ret & 1 => inverse pixel!
+// d=2*F(x,y) integer numbers, ideal line: ,I pixel: I@
+//   ..@  @@@  .@.  ...,@2@. +1..+3 floodfill around line ???
+//   ..@  .@@  .@.  ...,.@@@ +2..+4 <= that's not implemented yet
+//   ..@  ..@  .@.  ...,.@@@ +2..+4
+//   @.@  @..  .@.  ...,@@@. +1..+3
+//   @.@  @@.  .@.  ...I@@@.  0..+3
+//   @@@  @@@  .@.  ..@1@@..  0..+2
+//   90%   0%  100%   90%     r1-r2
+// I am not satisfied with it
+int get_line2(int x0, int y0, int x1, int y1, pix *p, int cs, int ret){
+   int dx,dy,incrE,incrNE,d,x,y,r0,r1,ty,tx,q,ddy,rx,ry,
+       *px,*py,*pdx,*pdy,*ptx,*pty,*px1;
+   dx=abs(x1-x0); tx=((x1>x0)?1:-1);    // tx=x-spiegelung (new)  
+   dy=abs(y1-y0); ty=((y1>y0)?1:-1);	// ty=y-spiegelung (new)
+   // rotate coordinate system if dy>dx
+   if(dx>dy){ pdx=&dx;pdy=&dy;px=&x;py=&y;ptx=&tx;pty=&ty;px1=&x1;rx=1;ry=0; }
+   else     { pdx=&dy;pdy=&dx;px=&y;py=&x;ptx=&ty;pty=&tx;px1=&y1;rx=0;ry=1; }
+   if( *ptx<0 ){ swap(&x0,&x1);swap(&y0,&y1);tx=-tx;ty=-ty; }
+   d=((*pdy)<<1)-(*pdx); incrE=(*pdy)<<1; incrNE=((*pdy)-(*pdx))<<1;  
+   x=x0; y=y0; r0=r1=0; ddy=3; // tolerance = bit 1 + bit 0 = left+right
+   // int t=(*pdx)/16,tl,tr;  // tolerance, left-,right delimiter
+   while( (*px)<=(*px1) ){  // not finaly implemented
+     q=((getpixel(p,x,y)<cs)?1:0)^(ret&1);
+     if ( !q ){		// tolerance one pixel perpenticular to the line
+                        // what about 2 or more pixels tolerance???
+       ddy&=(~1)|(((getpixel(p,x+ry,y+rx)<cs)?1:0)^(ret&1));
+       ddy&=(~2)|(((getpixel(p,x-ry,y-rx)<cs)?1:0)^(ret&1))*2;
+     } else ddy=3;
+     if( ddy ) r0++; else r1++;
+     (*px)++; if( d<=0 ){ d+=incrE; } else { d+=incrNE; (*py)+=(*pty); }
+   }
+   return (r0*(ret&~1))/(r0+r1); // ret==100 => percentage %
+}
+
+/* Look for dots in the rectangular region x0 <= x <= x1 and y0 <= y
+ <= y1 in pixmap p.  The two low order bits in mask indicate the color
+ of dots to look for: If mask==1 then look for black dots (where a
+ pixel value less than cs is considered black).  If mask==2 then look
+ for white dots.  If mask==3 then look for both black and white dots.
+ If the dots are found, the corresponding bits are set in the returned
+ value.  Heavily used by the engine ocr0*.cc */
+char get_bw(int x0, int x1, int y0, int y1, pix * p, int cs, int mask) {
+  char rc = 0;			// later with error < 2% (1 dot)
+  int x, y;
+
+  if (x0 < 0)        x0 = 0;
+  if (x1 >= p->x)    x1 = p->x - 1;
+  if (y0 < 0)        y0 = 0;
+  if (y1 >= p->y)    y1 = p->y - 1;
+
+  for ( y = y0; y <= y1; y++)
+    for ( x = x0; x <= x1; x++) {
+      rc |= ((getpixel(p, x, y) < cs) ? 1 : 2);	// break if rc==3
+      if ((rc & mask) == mask)
+	return mask;		// break loop
+    }
+  return (rc & mask);
+}
+
+/* more general Mar2000 (x0,x1,y0,y1 instead of x0,y0,x1,y1! (history))
+ * look for black crossings throw a line from x0,y0 to x1,y1 and count them
+ * follow line and count crossings ([white]-black-transitions)
+ *  ex: horizontal num_cross of 'm' would return 3
+ *
+ * fail for:  .a... a-to-b counts no transitions, but there is
+ *            ...#.
+ *            ..#..
+ *            .#..b
+ */
+int num_cross(int x0, int x1, int y0, int y1, pix *p, int cs) {
+  int rc = 0, col = 0, k, x, y, i, d;	// rc=crossings  col=0=white
+  int dx = x1 - x0, dy = y1 - y0;
+
+  d = MAX(abs(dx), abs(dy));
+  for (i = 0, x = x0, y = y0; i <= d; i++) {
+    if (d) {
+      x = x0 + i * dx / d;
+      y = y0 + i * dy / d;
+    }
+    k = ((getpixel(p, x, y) < cs) ? 1 : 0);	// 0=white 1=black
+    if (col == 0 && k == 1)  // found a white-black transition
+      rc++;
+    col = k;        // last color
+  }
+  return rc;
+}
+
+/* check if test matches pattern
+ *   possible pattern: "a-zA-Z0-9+--\\"  (x-y dont work for c>127)
+ *  return: 0 means dont fit, 1 means found
+ *   ToDo: wchar_t cc + matching UTF-8 pattern for nonASCII
+ */
+int my_strchr( char *pattern, wchar_t cc ) {
+  char *s1;
+  if (pattern==(char *)NULL) return 0;
+  
+  /* if (!(cc&0x80)) s1=strchr(pattern,(char)cc);  else */          
+  switch (cc) {
+    case '-':                  /* used as a special character */
+      s1=strstr(pattern,"--"); /* search string -- in pattern */
+      if (s1) return 1; break;
+    default:
+      s1=strstr(pattern,decode(cc, UTF8)); /* search string cc in pattern */
+      if (s1) return 1; /* cc simply matches */
+      /* single char not found, now check the ranges */
+      s1=pattern;
+      while (s1) {
+        s1=strchr(s1+1,'-');  /* look for next '-' */
+        if ((!s1) || (!s1[0]) || (!s1[1])) return 0; /* nothing found or end */
+        if (*(s1-1)=='-' || *(s1+1)=='-') continue; /* skip -- pattern */
+        if (*(s1-1)<=cc && *(s1+1)>=cc) return 1;  /* within range */
+      }
+  }
+  return 0;
+}
+
+/* set alternate chars and its weight, called from the engine
+    if a char is recognized to (weight) percent
+   can be used for filtering (only numbers etc)
+   often usefull if Il1 are looking very similar
+   should this function stay in box.c ???
+   weight is between 0 and 100 in percent, 100 means absolutely sure
+   - not final, not time critical (js)
+   - replace it by a string-function setaobj(*b,"string",weight)
+     and let call setac the setas function 
+ */
+
+int setas(struct box *b, char *as, int weight){
+  job_t *job=OCR_JOB;
+  int i,j;
+  if (b->num_ac > NumAlt || b->num_ac<0) {
+    fprintf(stderr,"\nDBG: There is something wrong with setas()!");
+    b->num_ac=0;
+  }
+  if (as==NULL) {
+    fprintf(stderr,"\nDBG: setas(NULL) makes no sense!"); return 0; }
+  if (as[0]==0) {
+    fprintf(stderr,"\nDBG: setas(\"\") makes no sense!"
+                   " x= %d %d", b->x0, b->y0);
+    // out_x(b);
+    return 0;
+  }
+
+  /* char filter (ex: only numbers) ToDo: cfilter as UTF-8 */
+  if (job->cfg.cfilter) {
+    /* do not accept chars which are not in the cfilter string */
+    if ( as[0]>0 && as[1]==0 )
+    if ( !my_strchr(job->cfg.cfilter,as[0]) ) return 0;
+  }
+#if 0 /* obsolete, done in setac */
+  /* not sure that this is the right place, but where else? */
+  if ( as[0]>0 && as[1]==0 )
+  if (b->modifier != SPACE  &&  b->modifier != 0) {
+    wchar_t newac;
+    newac = compose(as[0], b->modifier);
+    as = (char *)decode(newac, UTF8); /* was (const char *) */
+    if (newac == as[0]) { /* nothing composed */
+      fprintf(stderr, "\nDBG setas compose was useless %d %d",b->x0,b->y0);
+      // out_x(b);
+    }
+  }
+#endif
+  
+  /* only the first run gets the full weight */
+  weight=(100-job->tmp.n_run)*weight/100;
+  
+  /* remove same entries from table */
+  for (i=0;i<b->num_ac;i++) 
+    if (b->tas[i])
+      if (strcmp(as,b->tas[i])==0) break;
+  if (b->num_ac>0 && i<b->num_ac){
+    if (weight<=b->wac[i]) return 0; /* if found + less weight ignore it */
+    /* to insert the new weigth on the right place, we remove it first */
+    if (b->tas[i]) free(b->tas[i]); 
+    for (j=i;j<b->num_ac-1;j++){  /* shift lower entries */
+      b->tac[j]=b->tac[j+1]; /* copy the char   */
+      b->tas[j]=b->tas[j+1]; /* copy the pointer to the string */
+      b->wac[j]=b->wac[j+1]; /* copy the weight */
+    }
+    b->num_ac--; /* shrink table */
+  }
+  /* sorting and add it to the table */
+  for (i=0;i<b->num_ac;i++) if (weight>b->wac[i]) break;
+  if (b->num_ac<NumAlt-1) b->num_ac++; /* enlarge table */
+  for (j=b->num_ac-1;j>i;j--){  /* shift lower entries */
+    b->tac[j]=b->tac[j-1]; /* copy the char   */
+    b->tas[j]=b->tas[j-1]; /* copy the pointer to the string */
+    b->wac[j]=b->wac[j-1]; /* copy the weight */
+  }
+  if (i<b->num_ac) {    /* insert new entry */
+    b->tac[i]=0;        /* insert the char=0 ... */
+    b->tas[i]=(char *)malloc(strlen(as)+1);     /* ... string */
+    if (b->tas[i]) memcpy(b->tas[i],as,strlen(as)+1);
+    b->wac[i]=weight;   /* ... and its weight  */
+  }
+  if (i==0) b->c=b->tac[0];  /* char or 0 for string */
+  return 0;
+}
+
+/* ToDo: this function will be replaced by a call of setas() later */
+int setac(struct box *b, wchar_t ac, int weight){
+  int i,j;
+  job_t *job=OCR_JOB;
+  if ((!b) || b->num_ac > NumAlt || b->num_ac<0) {
+    fprintf(stderr,"\nDBG: This is a bad call to setac()!");
+    if(b && (job->cfg.verbose & 6)) out_x(b);
+    b->num_ac=0;
+  }
+  if (ac==0 || ac==UNKNOWN) {
+    fprintf(stderr,"\nDBG: setac(0) makes no sense!");
+    return 0;
+  }
+  /* char filter (ex: only numbers) ToDo: cfilter as UTF-8 */
+  if (job->cfg.cfilter) {
+    /* do not accept chars which are not in the cfilter string */
+    /* if ( ac>255 || !strchr(job->cfg.cfilter,(char)ac) ) return 0; */
+    if ( !my_strchr(job->cfg.cfilter,ac) ) return 0;
+  }
+  /* not sure that this is the right place, but where else? */
+  if (b->modifier != SPACE  &&  b->modifier != 0) {
+    wchar_t newac;
+    newac = compose(ac, b->modifier);
+    if (newac == ac) { /* nothing composed */
+      if(job->cfg.verbose & 7) 
+        fprintf(stderr, "\nDBG %s setac (%d,%d): compose was useless, wac=%d",
+                decode(ac,ASCII), b->x0, b->y0, weight);
+      /* if(job->cfg.verbose & 6) out_x(b); */
+    }
+    ac = newac;
+  }
+  
+  /* only the first run gets the full weight */
+  weight=(100-job->tmp.n_run)*weight/100;
+  
+  /* remove same entries from table */
+  for (i=0;i<b->num_ac;i++) if (ac==b->tac[i]) break;
+  if (b->num_ac>0 && i<b->num_ac){
+    if (weight<=b->wac[i]) return 0;
+    if (b->tas[i]) free(b->tas[i]); 
+    for (j=i;j<b->num_ac-1;j++){  /* shift lower entries */
+      b->tac[j]=b->tac[j+1]; /* copy the char   */
+      b->tas[j]=b->tas[j+1]; /* copy the pointer to the string */
+      b->wac[j]=b->wac[j+1]; /* copy the weight */
+    }
+    b->num_ac--; /* shrink table */
+  }
+  /* sorting it to the table */
+  for (i=0;i<b->num_ac;i++) if (weight>b->wac[i]) break;
+  if (b->num_ac<NumAlt-1) b->num_ac++; /* enlarge table */
+  for (j=b->num_ac-1;j>i;j--){  /* shift lower entries */
+    b->tac[j]=b->tac[j-1]; /* copy the char   */
+    b->tas[j]=b->tas[j-1]; /* copy the pointer to the string */
+    b->wac[j]=b->wac[j-1]; /* copy the weight */
+  }
+  if (i<b->num_ac) {    /* insert new entry */
+    b->tac[i]=ac;       /* insert the char ... */
+    b->tas[j]=NULL;     /* ... no string (?) */
+    b->wac[i]=weight;   /* ... and its weight  */
+  }
+  if (i==0) b->c=ac; /* store best result to b->c (will be obsolete) */
+
+  return 0;
+}
+
+/* test if ac in wac-table
+   usefull for contextcorrection and box-splitting
+   return 0 if not found
+   return wac if found (wac>0)
+ */
+int testac(struct box *b, wchar_t ac){
+  int i;
+  if (b->num_ac > NumAlt || b->num_ac<0) {
+    fprintf(stderr,"\n#DEBUG: There is something wrong with testac()!");
+    b->num_ac=0;
+  }
+  /* search entries in table */
+  for (i=0;i<b->num_ac;i++) if (ac==b->tac[i]) return b->wac[i];
+  return 0;
+}
+
+
+/* look for edges: follow a line from x0,y0 to x1,y1, record the
+ * location of each transition, and return their number.
+ * ex: horizontal num_cross of 'm' would return 6
+ * remark: this function is not used, obsolete? ToDo: remove?
+ */
+int follow_path(int x0, int x1, int y0, int y1, pix *p, int cs, path_t *path) {
+  int rc = 0, prev, x, y, i, d, color; // rc=crossings  col=0=white
+  int dx = x1 - x0, dy = y1 - y0;
+
+  d = MAX(abs(dx), abs(dy));
+  prev = getpixel(p, x0, y0) < cs;	// 0=white 1=black
+  path->start = prev;
+  for (i = 1, x = x0, y = y0; i <= d; i++) {
+    if (d) {
+      x = x0 + i * dx / d;
+      y = y0 + i * dy / d;
+    }
+    color = getpixel(p, x, y) < cs; // 0=white 1=black
+    if (color != prev){
+      if (rc>=path->max){
+	int n=path->max*2+10;
+	path->x = (int *) xrealloc(path->x, n*sizeof(int));
+	path->y = (int *) xrealloc(path->y, n*sizeof(int));
+	path->max = n;
+      }
+      path->x[rc]=x;
+      path->y[rc]=y;
+      rc++;
+    }      
+    prev = color;
+  }
+  path->num=rc;
+  return rc;
+}
+
+/* ToDo: only used in follow_path, which is obsolete, remove? */
+void *xrealloc(void *ptr, size_t size){
+  void *p;
+  p = realloc(ptr, size);
+  if (size>0 && (!p)){
+    fprintf(stderr, "insufficient memory");
+    exit(1);
+  }
+  return p;
+}
+
+/*
+ *  -------------------------------------------------------------
+ *  mark edge-points
+ *   - first move forward until b/w-edge
+ *   - more than 2 pixel?
+ *   - loop around
+ *     - if forward    pixel : go up, rotate right
+ *     - if forward no pixel : rotate left
+ *   - stop if found first 2 pixel in same order
+ *  go_along_the_right_wall strategy is very similar and used otherwhere
+ *  --------------------------------------------------------------
+ *  turmite game: inp: start-x,y, regel r_black=UP,r_white=RIght until border
+ * 	       out: last-position
+ * 
+ *  could be used to extract more features:
+ *   by counting stepps, dead-end streets ,xmax,ymax,ro-,ru-,lo-,lu-edges
+ * 
+ *   use this little animal to find features, I first was happy about it
+ *    but now I prefer the loop() function 
+ */
+
+void turmite(pix *p, int *x, int *y,
+	     int x0, int x1, int y0, int y1, int cs, int rw, int rb) {
+  int r;
+  if (outbounds(p, x0, y0))	// out of pixmap
+    return;
+  while (*x >= x0 && *y >= y0 && *x <= x1 && *y <= y1) {
+    r = ((getpixel(p, *x, *y) < cs) ? rb : rw);	// select rule 
+    switch (r) {
+      case UP: (*y)--; break;
+      case DO: (*y)++; break;
+      case RI: (*x)++; break;
+      case LE: (*x)--; break;
+      case ST:       break;
+      default:       assert(0);
+    }
+    if( r==ST ) break;	/* leave the while-loop */
+  }
+}
+
+/* search a way from p0 to p1 without crossing pixels of type t
+ *  only two directions, useful to test if there is a gap 's'
+ * labyrinth algorithm - do you know a faster way? */
+int joined(pix *p, int x0, int y0, int x1, int y1, int cs){
+  int t,r,x,y,dx,dy,xa,ya,xb,yb;
+  x=x0;y=y0;dx=1;dy=0;
+  if(x1>x0){xa=x0;xb=x1;} else {xb=x0;xa=x1;}
+  if(y1>y0){ya=y0;yb=y1;} else {yb=y0;ya=y1;}
+  t=((getpixel(p,x,y)<cs)?1:0);
+  for(;;){
+    if( t==((getpixel(p,x+dy,y-dx)<cs)?1:0)	// right free?
+     && x+dy>=xa && x+dy<=xb && y-dx>=ya && y-dx<=yb) // wall
+         { r=dy;dy=-dx;dx=r;x+=dx;y+=dy; } // rotate right and step forward
+    else { r=dx;dx=-dy;dy=r; } // rotate left
+    // fprintf(stderr," path xy %d-%d %d-%d %d %d  %d %d\n",xa,xb,ya,yb,x,y,dx,dy);
+    if( x==x1 && y==y1 ) return 1;
+    if( x==x0 && y==y0 && dx==1) return 0;
+  }
+  // return 0; // endless loop ?
+}
+
+/* move from x,y to direction r until pixel of color col is found
+ *   or maximum of l steps
+ * return the number of steps done */
+int loop(pix *p,int x,int y,int l,int cs,int col, DIRECTION r){ 
+  int i=0;
+  if(x>=0 && y>=0 && x<p->x && y<p->y){
+    switch (r) {
+    case UP:
+      for( ;i<l && y>=0;i++,y--)
+	if( (getpixel(p,x,y)<cs)^col )
+	  break;
+      break;
+    case DO:
+      for( ;i<l && y<p->y;i++,y++)
+	if( (getpixel(p,x,y)<cs)^col )
+	  break;
+      break;
+    case LE:
+      for( ;i<l && x>=0;i++,x--)
+	if( (getpixel(p,x,y)<cs)^col )
+	  break;
+      break;
+    case RI:
+      for( ;i<l && x<p->x;i++,x++)
+	if( (getpixel(p,x,y)<cs)^col )
+	  break;
+      break;
+    default:;
+    }
+  }
+  return i;
+}
+
+/* Given a point, frames a rectangle containing all points of the same
+ * color surrounding it, and mark these points.
+ *  ToDo:  obsolate and replaced by frame_vector
+ *
+ * looking for better algo: go horizontally and look for upper/lower non_marked_pixel/nopixel
+ * use lowest three bits for mark
+ *   - recursive version removed! AmigaOS has no Stack-OVL-Event
+ * run around the chape using laby-robot
+ * bad changes can lead to endless loop!
+ *  - this is not absolutely sure but mostly works well
+ *  diag - 0: only pi/2 direction, 1: pi/4 directions (diagonal)
+ *  mark - 3 bit marker, mark each valid pixel with it
+ */
+int frame_nn(pix *p, int  x,  int  y,
+             int *x0, int *x1, int *y0, int *y1,	// enlarge frame
+             int cs, int mark,int diag){
+#if 1 /* flood-fill to detect black objects, simple and faster? */
+  int rc = 0, dx, col, maxstack=0; static int overflow=0;
+  int bmax=1024, blen=0, *buf;  /* buffer as replacement for recursion stack */
+
+  /* check bounds */
+  if (outbounds(p, x, y))  return 0;
+  /* check if already marked (with mark since v0.4) */
+  if ((marked(p,x,y)&mark)==mark) return 0;
+
+  col = ((getpixel(p, x, y) < cs) ? 0 : 1);
+  buf=(int *)malloc(bmax*sizeof(int)*2);
+  if (!buf) { fprintf(stderr,"malloc failed (frame_nn)\n");return 0;}
+  buf[0]=x;
+  buf[1]=y;
+  blen=1;
+
+  g_debug(fprintf(stderr,"\nframe_nn   x=%4d y=%4d",x,y);)
+  for ( ; blen ; ) {
+    /* max stack depth is complexity of the object */
+    if (blen>maxstack) maxstack=blen;
+    blen--;             /* reduce the stack */
+    x=buf[blen*2+0];
+    y=buf[blen*2+1];
+    if (y < *y0) *y0 = y;
+    if (y > *y1) *y1 = y;
+    /* first go to leftmost pixel */
+    for ( ; x>0 && (col == ((getpixel(p, x-1, y) < cs) ? 0 : 1)) ; x--);
+    if ((marked(p,x,y)&mark)==mark) continue; /* already scanned */
+    for (dx=-1;dx<2;dx+=2) /* look at upper and lower line, left */
+    if (    diag && x<p->x && x-1>0 && y+dx >=0 && y+dx < p->y
+         &&  col != ((getpixel(p, x  , y+dx) < cs) ? 0 : 1)
+         &&  col == ((getpixel(p, x-1, y+dx) < cs) ? 0 : 1)
+         && !((marked(p,x-1,y+dx)&mark)==mark)
+       ) {
+       if (blen+1>=bmax) { overflow|=1; continue; }
+       buf[blen*2+0]=x-1;
+       buf[blen*2+1]=y+dx;
+       blen++;
+    }
+    if (x < *x0) *x0 = x;
+    /* second go right, mark and get new starting points */ 
+    for ( ; x<p->x && (col == ((getpixel(p, x  , y) < cs) ? 0 : 1)) ; x++) {
+      p->p[x + y * p->x] |= (mark & 7);    rc++;  /* mark pixel */
+      /* enlarge frame */
+      if (x > *x1) *x1 = x;
+      for (dx=-1;dx<2;dx+=2) /* look at upper and lower line */
+      if (     col == ((getpixel(p, x  , y+dx) < cs) ? 0 : 1)
+        && (
+               col != ((getpixel(p, x-1, y   ) < cs) ? 0 : 1)
+            || col != ((getpixel(p, x-1, y+dx) < cs) ? 0 : 1) )
+        && !((marked(p,x,y+dx)&mark)==mark) && y+dx<p->y && y+dx>=0
+         ) {
+         if (blen+1>=bmax) { overflow|=1; continue; }
+         buf[blen*2+0]=x;
+         buf[blen*2+1]=y+dx;
+         blen++;
+      }
+    }
+    for (dx=-1;dx<2;dx+=2) /* look at upper and lower line, right */
+    if (    diag  && x<p->x && x-1>0 && y+dx >=0 && y+dx < p->y
+         &&  col == ((getpixel(p, x-1, y   ) < cs) ? 0 : 1)
+         &&  col != ((getpixel(p, x  , y   ) < cs) ? 0 : 1)
+         &&  col != ((getpixel(p, x-1, y+dx) < cs) ? 0 : 1)
+         &&  col == ((getpixel(p, x  , y+dx) < cs) ? 0 : 1)
+         && !((marked(p,x,y+dx)&mark)==mark) 
+       ) {
+       if (blen+1>=bmax) { overflow|=1; continue; }
+       buf[blen*2+0]=x;
+       buf[blen*2+1]=y+dx;
+       blen++;
+    }
+  }
+  
+  /* debug, ToDo: use info maxstack and pixels for image classification */
+  g_debug(fprintf(stderr," maxstack= %4d pixels= %6d",maxstack,rc);)
+  if (overflow==1){
+    overflow|=2;
+    fprintf(stderr,"# Warning: frame_nn stack oerflow\n");
+  }
+  free(buf);
+#else /* old version, ToDo: improve it for tmp04/005*.pgm.gz */
+  int i, j, d, dx, ox, oy, od, nx, ny, rc = 0, rot = 0, x2 = x, y2 = y, ln;
+
+  static const int d0[8][2] = { { 0, -1} /* up    */, {-1, -1}, 
+				{-1,  0} /* left  */, {-1,  1}, 
+				{ 0,  1} /* down  */, { 1,  1}, 
+				{ 1,  0} /* right */, { 1, -1}};
+
+  /* check bounds */
+  if (outbounds(p, x, y))
+    return 0;
+  /* check if already marked */
+  if ((marked(p,x,y)&mark)==mark) 
+    return 0;
+
+  i = ((getpixel(p, x, y) < cs) ? 0 : 1);
+  rc = 0;
+
+  g_debug(fprintf(stderr," start frame:");)
+
+  for (ln = 0; ln < 2 && rot >= 0; ln++) {  // repeat if right-loop 
+    g_debug(fprintf(stderr," ln=%d diag=%d cs=%d x=%d y=%d - go to border\n",ln,diag,cs,x,y);)
+    
+    od=d=(8+4*ln-diag)&7; // start robot looks up, right is a wall
+    // go to right (left) border
+    if (ln==1) { 
+      x=x2;	y=y2; 
+    } 
+    /* start on leftmost position */
+    for (dx = 1 - 2*ln; x + dx < p->x && x + dx >= 0 /* bounds */ &&
+      	      	       i == ((getpixel(p, x + dx, y) < cs) ? 0 : 1) /* color */; 
+	      	       x += dx);
+
+    g_debug(fprintf(stderr," ln=%d diag=%d cs=%d x=%d y=%d\n",ln,diag,cs,x,y);)
+
+    /* robot stores start-position */
+    ox = x;	oy = y;
+    for (rot = 0; abs(rot) <= 64; ) {	/* for sure max. 8 spirals */
+      /* leftmost position */
+      if (ln == 0 && x < x2) {
+	x2 = x; 	y2 = y;
+      }	
+
+      g_debug(fprintf(stderr," x=%3d y=%3d d=%d i=%d p=%3d rc=%d\n",x,y,d,i,getpixel(p,x,y),rc);)
+
+      if ( abs(d0[d][1]) ) {	/* mark left (right) pixels */
+	for (j = 0, dx = d0[d][1]; x + j >= 0 && x + j < p->x
+	              	&& i == ((getpixel(p, x + j, y) < cs) ? 0 : 1); j += dx) {
+	  if (!((marked(p, x + j, y)&mark)==mark))
+	    rc++;
+	  p->p[x + j + y * p->x] |= (mark & 7);
+	}
+      }
+      /* look to the front of robot */
+      nx = x + d0[d][0];
+      ny = y + d0[d][1];
+      /* if right is a wall */
+      if ( outbounds(p, nx, ny) || i != ((getpixel(p,nx,ny)<cs) ? 0 : 1) ) {
+	/* rotate left */
+        d=(d+2-diag) & 7; rot-=2-diag;
+      }
+      else {	/* if no wall, go, turn back and rotate left */
+        x=nx; y=ny; d=(d+4+2-diag) & 7; rot+=2-diag+4;
+	/* enlarge frame */
+        if (x < *x0)      *x0 = x;
+	if (x > *x1)	  *x1 = x;
+	if (y < *y0)	  *y0 = y;
+	if (y > *y1)	  *y1 = y;
+      } 
+      if(x==ox && y==oy && d==od) break;	// round trip finished
+    }
+  }
+  g_debug(fprintf(stderr," rot=%d\n",rot);)
+#endif
+  return rc;
+}
+
+/*   obsolete! replaced by vectors
+ * mark neighbouring pixel of same color, return number
+ * better with neighbours of same color (more general) ???
+ * parameters: (&~7)-pixmap, start-point, critical_value, mark
+ *  recursion is removed */
+int mark_nn(pix * p, int x, int y, int cs, int r) {
+  /* out of bounds or already marked? */
+  if (outbounds(p, x, y) || (marked(p, x, y)&r)==r) 
+    return 0;
+  {
+    int x0, x1, y0, y1;
+    x0 = x1 = x;
+    y0 = y1 = y;			// not used
+    return frame_nn(p, x, y, &x0, &x1, &y0, &y1, cs, r, OCR_JOB->tmp.n_run & 1);
+    // using same scheme
+  }
+}
+
+/* ToDo: finish to replace old frame by this new one
+ *
+ *   @...........#@@@@@@@.  # = marked as already scanned black pixels
+ *   @........@@@@@@@@@@@#      only left and right border
+ *   .......#@@@@@@@@@@@@@        left side on even y
+ *   ......@@@@@@@@#.@@@@#        right side on odd y
+ *   .....#@@@@@......#@@@   no border is marked twice
+ *   ....@@@@@#......@@@#.   works also for thinn lines
+ *   ...#@@@@........#@@@. - outer loop is stored as first
+ *   ..@@@@#........@@@#.. - inner loop is stored as second
+ *   .#@@@@........#@@@@..    1st in an extra box (think on white chars)
+ *   @@@@#.......@@@@#....    2nd merge in an extra step
+ *   #@@@@@....#@@@@@.....
+ *   @@@@@@@@@@@@@@#......
+ *   .#@@@@@@@@@@@@.......
+ *
+ * run around the chape using laby-robot
+ *  - used for scanning boxes, look for horizontal b/w transitions
+ *    with unmarked black pixels and call this routine
+ *  - stop if crossing a marked box in same direction (left=up, right=down)
+ *  box  - char box, store frame_vectors and box
+ *  x,y  - starting point
+ *  mark - 3 bit marker, mark each valid pixel with it
+ *  diag - 0: only pi/2 direction, 1: pi/4 directions (diagonal)
+ *  ds   - start direction, 6=right of right border, 2=left of left border
+ *  ret  - 0=ok, -1=already marked, -2=max_num_frames_exceeded
+ *               -7=no border in direction ds
+ */
+#if 0
+#undef g_debug
+#define g_debug(x) x
+#endif
+/* grep keywords: scan_vectors frame_vector */
+int frame_vector(struct box *box1, int  x,  int  y,
+             int cs, int mark, int diag, int ds) {
+  int i1, i2, i2o,
+      new_x=1,    /* flag for storing the vector x,y */
+      steps=1,    /* steps between stored vectors, speedup for big frames */
+      d,          /* direction */ 
+      ox, oy,     /* starting point */
+      nx, ny, mx, my, /* used for simplification */
+      /* ToDo: add periphery to box (german: Umfang?) */
+      rc  = 1,    /* return code, circumference, sum vector lengths */
+      rot = 0,    /* memory for rotation, rot=8 means one full rotation */
+      vol = 0;    /* volume inside frame, negative for white inside black */
+  pix *p=box1->p;
+
+  /* translate the 8 directions to (x,y) pairs,
+   * if only four directions are used, only every 2nd vector is accessed,
+   * +1 turn left, -1 turn right
+   */
+  static const int d0[8][2] =
+    { { 0, -1}, /* up    */  {-1, -1},   /* up-le */
+      {-1,  0}, /* left  */  {-1,  1},   /* do-le */
+      { 0,  1}, /* down  */  { 1,  1},   /* do-ri */
+      { 1,  0}, /* right */  { 1, -1} }; /* up-ri */
+
+  /* check bounds */
+  if (outbounds(p, x, y))
+    return 0;
+
+  /* pixel color we are looking for, 0=black, 1=white */
+  d = ds;
+  i1 = ((getpixel(p, x,            y           ) < cs) ? 0 : 1);
+  i2 = ((getpixel(p, x + d0[d][0], y + d0[d][1]) < cs) ? 0 : 1);
+
+  g_debug(fprintf(stderr,"\nLEV2 frame_vector @ %3d %3d  d%d %2d %2d"
+    "  %d-%d pix=%3d mark=%d cs=%d",\
+    x,y,ds,d0[ds][0],d0[ds][1],i1,i2,getpixel(p,x,y),mark,cs);)
+                  
+  if (i1==i2){
+    fprintf(stderr,"ERROR frame_vector: no border\n");
+    return -7; /* no border detected */
+  }
+
+  /* initialize boxframe outside this function
+     box1->x0=box1->x1=x;
+     box1->y0=box1->y1=y;
+  */
+  
+  /* initialize boxvector outside this function
+     box1->num_frames=0
+     num_frame_vectors[0]=0 ???
+     and store start value
+   */
+  if (box1->num_frames >= MaxNumFrames) return -2;
+  /* index to next (x,y) */
+  i2o=i2=( (box1->num_frames==0)?0: 
+            box1->num_frame_vectors[ box1->num_frames ] );
+#if 0 // obsolete v0.43
+  box1->frame_vector[i2][0]=x;
+  box1->frame_vector[i2][1]=y;
+  i2++;
+  box1->num_frame_vectors[ box1->num_frames ]=i2;
+#endif
+  box1->num_frames++;
+   
+  /* robot stores start-position */
+  ox = x;  oy = y; /* look forward to white pixel */
+
+  for (;;) {	/* stop if same marked pixel touched */
+
+    g_debug(fprintf(stderr,"\nLEV3: x= %3d %3d d= %d rot= %2d  %3d",x,y,d,rot,i2);)
+    
+    /* ToDo: store max. abs(rot) ??? for better recognition */
+    if (new_x) {
+      g_debug(fprintf(stderr,"\nLEV2: markB xy= %3d %3d ", x, y);)
+      p->p[x + y * p->x] |= (mark & 7); /* mark black pixel */
+    }
+
+    /* store a new vector or enlarge the predecessor */
+    if (new_x && (rc%steps)==0) { /* dont store everything on big chars */
+      if (i2>=MaxFrameVectors) {
+        box1->num_frame_vectors[ box1->num_frames-1 ]=i2;
+        reduce_vectors(box1,1);    /* simplify loop */ 
+        i2=box1->num_frame_vectors[ box1->num_frames-1 ];
+        /* enlarge steps on big chars getting speedup */
+        steps=(box1->y1-box1->y0+box1->x1-box1->x0)/32+1;
+      }
+      /* store frame-vector */
+      if (i2<MaxFrameVectors) {
+        box1->frame_vector[i2][0]=x;
+        box1->frame_vector[i2][1]=y;
+        /* test if older vector points to the same direction */
+        if (i2>1) {
+          /* get predecessor */
+          nx=box1->frame_vector[i2-1][0]-box1->frame_vector[i2-2][0];
+          ny=box1->frame_vector[i2-1][1]-box1->frame_vector[i2-2][1];
+          mx=x                          -box1->frame_vector[i2-1][0];
+          my=y                          -box1->frame_vector[i2-1][1];
+          /* same direction? */
+          if (nx*my-ny*mx==0 && nx*mx>=0 && ny*my>=0) {
+            /* simplify by removing predecessor */
+            i2--;
+            box1->frame_vector[i2][0]=x;
+            box1->frame_vector[i2][1]=y;
+          } /* do not simplify */
+        }
+        i2++;
+        box1->num_frame_vectors[ box1->num_frames-1 ]=i2;
+      }
+      g_debug(fprintf(stderr," stored @ %3d steps= %d", i2-1, steps);)
+    } 
+    new_x=0; /* work for new pixel (x,y) done */
+
+    /* check if round trip is finished */
+    if (x==ox && y==oy && abs(rot)>=8) break;
+
+    /* look to the front of robot (turtle or ant) */
+    nx = x + d0[d][0];
+    ny = y + d0[d][1];
+    
+    /* next step, if right is a wall turn the turtle left */
+    if ( outbounds(p, nx, ny) || i1 != ((getpixel(p,nx,ny)<cs) ? 0 : 1) ) {
+      if (y==ny && nx>=0 && nx<p->x) { /* if inbound */
+        g_debug(fprintf(stderr,"\nLEV2: markW xy= %3d %3d ", nx, ny);)
+        p->p[nx + ny * p->x] |= (mark & 7); /* mark white pixel */
+      }
+      /* rotate left 90 or 45 degrees */
+      d=(d+2-diag) & 7; rot+=2-diag;
+      /* calculate volume inside frame */
+      switch (d+diag) {
+        case 2+2: vol-=x-1; break;
+        case 6+2: vol+=x;   break;
+      }
+    }
+    else { /* if no wall, go forward and turn right (90 or 45 degrees) */
+      x=nx; y=ny;
+      /* turn back and rotate left */
+      d=(d+4+2-diag) & 7; rot+=2-diag-4;
+      rc++; /* counting steps, used for speedup */
+
+      /* enlarge frame */
+      if (x < box1->x0) box1->x0 = x;
+      if (x > box1->x1) box1->x1 = x;
+      if (y < box1->y0) box1->y0 = y;
+      if (y > box1->y1) box1->y1 = y;
+      
+      new_x=1;
+    }
+  }
+
+  /* to distinguish inner and outer frames, store volume as +v or -v */
+  box1->frame_vol[ box1->num_frames-1 ] = vol;
+  box1->frame_per[ box1->num_frames-1 ] = rc-1;
+
+  /* dont count and store the first vector twice */
+  if (i2-i2o>1) { 
+    i2--; rc--; box1->num_frame_vectors[ box1->num_frames-1 ]=i2;
+  }
+  /* output break conditions */
+  g_debug(fprintf(stderr,"\nLEV2 o= %3d %3d  x= %3d %3d  r=%d v=%d",ox,oy,x,y,rot,vol);)
+  /* rc=1 for a single point, rc=2 for a two pixel sized point */
+  g_debug(fprintf(stderr," steps= %3d vectors= %3d",rc,i2);)
+  /* out_x(box1); ToDo: output only the first thousend */
+  return rc; /* return number of bordering pixels = periphery? */
+}
+
+
+
+/* clear lowest 3 (marked) bits (they are used for marking) */ 
+void clr_bits(pix * p, int x0, int x1, int y0, int y1) {
+  int x, y;
+  for ( y=y0; y <= y1; y++)
+    for ( x=x0; x <= x1; x++)
+      p->p[x+y*p->x] &= ~7;
+}
+
+/* look for white holes surrounded by black points
+ * at the moment look for white point with black in all four directions
+ *  - store position of hole in coordinates relativ to box!
+ *  ToDo: count only holes with vol>10% ???
+ * ToDo: rewrite for frame vectors (faster, no malloc)
+ *       holes are frames rotating left hand 
+ *  obsolete, do it with vectors
+ */
+int num_hole(int x0, int x1, int y0, int y1, pix * p, int cs, holes_t *holes) {
+  int num_holes = 0, x, y, hole_size;
+  pix b;			// temporary mini-page
+  int dx = x1 - x0 + 1, dy = y1 - y0 + 1;
+  unsigned char *buf;	//  2nd copy of picture, for working 
+
+  if (holes) holes->num=0;
+  if(dx<3 || dy<3) return 0;
+  b.p = buf = (unsigned char *) malloc( dx * dy );
+  if( !buf ){
+    fprintf( stderr, "\nFATAL: malloc(%d) failed, skip num_hole", dx*dy );
+    return 0;
+  }
+  if (copybox(p, x0, y0, dx, dy, &b, dx * dy))
+    { free(b.p); return -1;}
+
+  // printf(" num_hole(");
+  /* --- mark white-points connected with border */
+  for (x = 0; x < b.x; x++) {
+    if (getpixel(&b, x, 0) >= cs)
+      mark_nn(&b, x, 0, cs, AT);
+    if (getpixel(&b, x, b.y - 1) >= cs)
+      mark_nn(&b, x, b.y - 1, cs, AT);
+  }
+  for (y = 0; y < b.y; y++) {
+    if (getpixel(&b, 0, y) >= cs)
+      mark_nn(&b, 0, y, cs, AT);
+    if (getpixel(&b, b.x - 1, y) >= cs)
+      mark_nn(&b, b.x - 1, y, cs, AT);
+  }
+
+  g_debug(out_b(NULL,&b,0,0,b.x,b.y,cs);)
+  // --- look for unmarked white points => hole
+  for (x = 0; x < b.x; x++)
+    for (y = 0; y < b.y; y++)
+      if (!((marked(&b, x, y)&AT)==AT))	// unmarked
+	if (getpixel(&b, x, y) >= cs) {	// hole found
+#if 0
+	  hole_size=mark_nn(&b, x, y, cs, AT);  /* old version */
+	  if (hole_size > 1 || dx * dy <= 40)
+	    num_holes++;
+#else
+          {    /* new version, for future store of hole characteristics */
+            int x0, x1, y0, y1, i, j;
+            x0 = x1 = x;
+            y0 = y1 = y;			// not used
+            hole_size=frame_nn(&b, x, y, &x0, &x1, &y0, &y1, cs, AT, OCR_JOB->tmp.n_run & 1);
+            // store hole for future use, num is initialized with 0
+ 	    if (hole_size > 1 || dx * dy <= 40){
+ 	      num_holes++;
+              if (holes) {
+                // sort in table
+                for (i=0;i<holes->num && i<MAX_HOLES;i++)
+                  if (holes->hole[i].size < hole_size) break;
+                for (j=MAX_HOLES-2;j>=i;j--)
+                  holes->hole[j+1]=holes->hole[j];
+                if (i<MAX_HOLES) {
+                  // printf("  i=%d size=%d\n",i,hole_size);
+                  holes->hole[i].size=hole_size;
+                  holes->hole[i].x=x;
+                  holes->hole[i].y=y;
+                  holes->hole[i].x0=x0;
+                  holes->hole[i].y0=y0;
+                  holes->hole[i].x1=x1;
+                  holes->hole[i].y1=y1;
+                }
+                holes->num++;
+              }
+            }
+          }
+#endif
+	}
+  free(b.p);
+  // printf(")=%d",num_holes);
+  return num_holes;
+}
+
+/* count for black nonconnected objects --- used for i,auml,ouml,etc. */
+/* ToDo: obsolete, replaced by vectors and box.num_boxes */
+int num_obj(int x0, int x1, int y0, int y1, pix * p, int cs) {
+  int x, y, rc = 0;		// rc=num_obj
+  unsigned char *buf; // 2nd copy of picture, for working
+  pix b;
+
+  if(x1<x0 || y1<y0) return 0;
+  b.p = buf = (unsigned char *) malloc( (x1-x0+1) * (y1-y0+1) );
+  if( !buf ){
+    fprintf( stderr, "\nFATAL: malloc(%d) failed, skip num_obj",(x1-x0+1)*(y1-y0+1) );
+    return 0;
+  }
+  if (copybox(p, x0, y0, x1 - x0 + 1, y1 - y0 + 1, &b, (x1-x0+1) * (y1-y0+1)))
+    { free(b.p); return -1; }
+  // --- mark black-points connected with neighbours
+  for (x = 0; x < b.x; x++)
+    for (y = 0; y < b.y; y++)
+      if (getpixel(&b, x, y) < cs)
+	if (!((marked(&b, x, y)&AT)==AT)) {
+	  rc++;
+	  mark_nn(&b, x, y, cs, AT);
+	}
+  free(b.p);
+  return rc;
+}
+
+#if 0
+// ----------------------------------------------------------------------
+// first idea for making recognition based on probability
+//  - start with a list of all possible chars
+//  - call recognition_of_char(box *) 
+//    - remove chars from list which could clearly excluded
+//    - reduce probability of chars which have wrong features
+//  - font types list could also build
+// at the moment it is only an idea, I should put it to the todo list
+//  
+char *list="0123456789,.\0xe4\0xf6\0xfc"	// "a=228 o=246 u=252
+           "abcdefghijklmnopqrstuvwxyz"
+           "ABCDEFGHIJKLMNOPQRSTUVWXYZ";
+int  wert[100];
+int  listlen=0,numrest=0;
+// initialize a new character list (for future)
+void ini_list(){ int i;
+    for(i=0;list[i]!=0 && i<100;i++) wert[i]=0;
+    numrest=listlen=i; } 
+// exclude??? (for future) oh it was long time ago, I wrote that :/
+void exclude(char *filt){ int i,j;
+    for(j=0;filt[j]!=0 && j<100;j++)
+    for(i=0;list[i]!=0 && i<100;i++)
+    if( filt[j]==list[i] ) { if(!wert[i])numrest--; wert[i]++; } }
+// get the result after all the work (for future)
+char getresult(){ int i;
+    if( numrest==1 )
+    for(i=0;list[i]!=0 && i<100;i++) if(!wert[i]) return list[i];
+    return '_';
+ }
+#endif
+
+//  look at the environment of the pixel too (contrast etc.)
+//   detailed analysis only of diff pixels!
+//
+// 100% * "distance", 0 is ideal fit
+// = similarity of two chars for recognition of garbled (verstuemmelter) chars
+//   weight of pixels with only one same neighbour set to 0
+//   look at contours too! v0.2.4: B==H
+// changed for v0.41, Mar06
+int distance( pix *p1, struct box *box1,
+              pix *p2, struct box *box2, int cs){
+   int rc=0,x,y,v1,v2,i1,i2,rgood=0,rbad=0,x1,y1,x2,y2,dx,dy,dx1,dy1,dx2,dy2;
+   x1=box1->x0;y1=box1->y0;x2=box2->x0;y2=box2->y0;
+   dx1=box1->x1-box1->x0+1; dx2=box2->x1-box2->x0+1; dx=((dx1>dx2)?dx1:dx2);
+   dy1=box1->y1-box1->y0+1; dy2=box2->y1-box2->y0+1; dy=((dy1>dy2)?dy1:dy2);
+   if(abs(dx1-dx2)>1+dx/16 || abs(dy1-dy2)>1+dy/16) return 100;
+   // compare relations to baseline and upper line
+   if(2*box1->y1>box1->m3+box1->m4 && 2*box2->y1<box2->m3+box2->m4) rbad+=128;
+   if(2*box1->y0>box1->m1+box1->m2 && 2*box2->y0<box2->m1+box2->m2) rbad+=128;
+   // compare pixels
+   for( y=0;y<dy;y++ )
+   for( x=0;x<dx;x++ ) {	// try global shift too ???
+     v1     =((getpixel(p1,x1+x  ,y1+y  )<cs)?1:0); i1=8;	// better gray?
+     v2     =((getpixel(p2,x2+x  ,y2+y  )<cs)?1:0); i2=8;	// better gray?
+     if(v1==v2) { rgood+=8; continue; } // all things are right!
+     // what about different pixel???
+     // test overlap of 8 surounding pixels ??? bad if two nb. are bad
+     v1=-1;
+     for(i1=-1;i1<2;i1++)
+     for(i2=-1;i2<2;i2++)if(i1!=0 || i2!=0){
+       if( ((getpixel(p1,x1+x+i1*(1+dx/32),y1+y+i2*(1+dy/32))<cs)?1:0)
+         !=((getpixel(p2,x2+x+i1*(1+dx/32),y2+y+i2*(1+dy/32))<cs)?1:0) ) v1++;
+     }
+     if (v1>0) rbad+=16*v1;
+     else      rbad++;    
+   }
+   if(rgood+rbad) rc= (100*rbad+(rgood+rbad-1))/(rgood+rbad); else rc=99;
+   if(rc<10 && OCR_JOB->cfg.verbose & 7){
+     fprintf(stderr,"\n#  distance rc=%d good=%d bad=%d",rc,rgood,rbad);
+//     out_x(box1);out_x(box2);
+   }
+   return rc;
+}
+
+
+
+// ============================= call OCR engine ================== ;)
+//  nrun=0 from outside, nrun=1 from inside (allows modifications, oobsolete)
+wchar_t whatletter(struct box *box1, int cs, int nrun){
+   wchar_t bc=UNKNOWN;			// best letter
+   wchar_t um=SPACE;			// umlaut? '" => modifier
+   pix *p=box1->p;   // whole image
+   int	x,y,dots,xa,ya,x0,x1,y0,y1,dx,dy,i;
+   pix b;            // box
+   struct box bbuf=*box1;  // restore after modifikation!
+
+   if (box1->num_ac>0 && box1->wac[0]>=OCR_JOB->cfg.certainty && bc==UNKNOWN) {
+      bc=box1->tac[0];
+   }
+   // if (bc!=UNKNOWN) return bc;
+   // if whatletter() called again, only unknown chars are processed
+   // bad for splitting!
+   
+   // store box data, which can be modified for modified chars in 2nd run
+   bbuf.x0=box1->x0; bbuf.y0=box1->y0;
+   bbuf.x1=box1->x1; bbuf.y1=box1->y1;
+   
+   xa=box1->x;  ya=box1->y;
+   x0=box1->x0; y0=box1->y0;
+   x1=box1->x1; y1=box1->y1;
+   // int vol=(y1-y0+1)*(x1-x0+1);	// volume
+   // crossed l-m , divided chars
+   while( get_bw(x0,x1,y0,y0,p,cs,1)!=1 && y0+1<y1) y0++;
+   while( get_bw(x0,x1,y1,y1,p,cs,1)!=1 && y0+1<y1) y1--;
+   dx=x1-x0+1;
+   dy=y1-y0+1;	// size
+
+   // better to proof the white frame too!!! ????
+   // --- test for german umlaut and points above, not robust enough???
+   // if three chars are connected i-dots (ari) sometimes were not detected
+   //  - therefore after division a test could be useful
+   // modify y0 only in second run!?
+   // we need it here to have the right copybox
+   if (um==SPACE && dy>5 && box1->num_boxes>1)
+     testumlaut(box1,cs,2,&um); /* set box1->modifier + new y0 */
+
+   dots=box1->dots;
+   y0  =box1->y0;	// dots==2 => y0 below double dots
+   dy  =y1-y0+1;
+
+   // move upper and lower border (for divided letters)
+   while( get_bw(x0,x1,y0,y0,p,cs,1)==0  &&  y0+1<y1) y0++;
+   while( get_bw(x0,x1,y1,y1,p,cs,1)==0  &&  y0+1<y1) y1--;
+   while( get_bw(x0,x0,y0,y1,p,cs,1)==0  &&  x0+1<x1) x0++;
+   while( get_bw(x1,x1,y0,y1,p,cs,1)==0  &&  x0+1<x1) x1--;
+   dx=x1-x0+1;
+   dy=y1-y0+1;	// size
+   box1->x0=x0; box1->y0=y0;	// set reduced frame
+   box1->x1=x1; box1->y1=y1;
+
+   // set good startpoint (probably bad from division)?
+   if( xa<x0 || xa>x1 || ya<y0 || ya>y1
+     || getpixel(p,xa,ya)>=cs /* || 2*ya<y0+y1 */ || dots>0 ){
+     // subfunction? also called after division of two glued chars?
+     for(y=y1;y>=y0;y--) // low to high (not i-dot)
+     for(x=(x0+x1)/2,i=0;x>=x0 && x<=x1;i++,x+=((2*i&2)-1)*i) /* is that ok? */
+     if (getpixel(p,x,y)<cs && (getpixel(p,x+1,y)<cs
+                             || getpixel(p,x,y+1)<cs)){ xa=x;ya=y;y=-1;break; }
+     /* should box1->x,y be set? */
+   }
+
+   // ----- create char-only-box -------------------------------------
+   // ToDo: this will be obsolete if vectors are used only
+   if(dx<1 || dy<1) return bc; /* should not happen */
+   b.p = (unsigned char *) malloc( dx * dy );
+   if (!b.p) fprintf(stderr,"Warning: malloc failed L%d\n",__LINE__);
+   if( copybox(p,x0,y0,dx,dy,&b,dx*dy) ) 
+     { free(b.p); return bc; }
+   // clr_bits(&b,0,b.x-1,0,b.y-1);
+   // ------ use diagonal too (only 2nd run?) 
+   /* following code failes on ! and ?  obsolete if vectors are used
+      ToDo:
+       - mark pixels neighoured to pixels outside and remove them from &b
+         v0.40
+         will be replaced by list of edge vectors  
+       - mark accents, dots and remove them from &b
+    */
+#if 1 /* becomes obsolate by vector code */
+   if (y0>0)  // mark upper overlap
+   for ( x=x0; x<=x1; x++) {
+     if (getpixel(p,x,y0-1)<cs
+      && getpixel(p,x,y0  )<cs && (marked(&b,x-x0,0)&1)!=1)
+     mark_nn(&b,x-x0,0,cs,1);
+   }
+   if (x0>0)  // mark left overlap
+   for ( y=y0; y<=y1; y++) {
+     if (getpixel(p,x0-1,y)<cs
+      && getpixel(p,x0  ,y)<cs && (marked(&b,0,y-y0 )&1)!=1)
+     mark_nn(&b,0,y-y0,cs,1);
+   }
+   if (x1<p->x-1)  // mark right overlap
+   for ( y=y0; y<=y1; y++) {
+     if (getpixel(p,x1+1,y)<cs
+      && getpixel(p,x1  ,y)<cs && (marked(&b,x1-x0,y-y0)&1)!=1)
+     mark_nn(&b,x1-x0,y-y0,cs,1);
+   }
+   mark_nn(&b,xa-x0,ya-y0,cs,2); // not glued chars
+   for(x=0;x<b.x;x++)
+   for(y=0;y<b.y;y++){
+     if (  (marked(&b,x,y  )&3)==1 && getpixel(&b,x,y  )<cs )
+     b.p[x+y*b.x] = 255&~7;  /* reset pixel */
+   }
+#endif
+   
+   // if (bc == UNKNOWN)   // cause split to fail
+   bc=ocr0(box1,&b,cs);
+
+   /* ToDo: try to change pixels near cs?? or melt? */
+   if (box1->num_ac>0 && box1->wac[0]>=OCR_JOB->cfg.certainty && bc==UNKNOWN) {
+     bc=box1->tac[0];
+   }
+
+   if (um!=0 && um!=SPACE && bc<127) {  /* ToDo: is that obsolete now? */
+     wchar_t newbc;
+     newbc = compose(bc, um );
+     if (newbc == bc) { /* nothing composed */
+       if(OCR_JOB->cfg.verbose & 7) 
+         fprintf(stderr, "\nDBG whatletter: compose(%s) was useless (%d,%d)",
+           decode(bc,ASCII), box1->x0, box1->y0);
+       // if(OCR_JOB->cfg.verbose & 6) out_x(box1);
+     }
+     bc = newbc;
+   }
+   // restore modified boxes
+   box1->x0=bbuf.x0; box1->y0=bbuf.y0;
+   box1->x1=bbuf.x1; box1->y1=bbuf.y1;
+//   if (box1->c==UNKNOWN) out_b(box1,&b,0,0,dx,dy,cs); // test
+
+   free(b.p);
+   return bc;
+}
+
+/*
+** creates a list of boxes/frames around objects detected 
+** on the pixmap p for further work
+** returns number of boxes created.
+** - by the way: get average X, Y (avX=sumX/numC,..)
+*/
+int scan_boxes( job_t *job, pix *p ){
+  int x, y, nx, cs, rc, ds;
+  struct box *box3;
+  //  job_t *job=OCR_JOB; /* fixme */
+
+  if (job->cfg.verbose)
+    fprintf(stderr,"# scan_boxes");
+
+  cs = job->cfg.cs;
+  job->res.sumX = job->res.sumY = job->res.numC = 0;
+
+  /* clear the lowest bits of each pixel, later used as "scanned"-marker */
+  clr_bits( p, 0, p->x - 1, 0, p->y - 1);
+
+  for (y=0; y < p->y; y++)
+    for (x=0; x < p->x; x++)
+    for (ds=2; ds<7; ds+=4) { // NO - dust of size 1 is not removed !!!
+      nx=x+((ds==2)?-1:+1);
+      if (nx<0 || nx>=p->x) continue; /* out of image, ex: recframe */
+      if ( getpixel(p, x,y)>=cs || getpixel(p,nx,y)< cs)  // b/w transition?
+	continue;
+      if ((marked(p, x,y) & 1)&&(marked(p, nx, y) & 1))
+	continue;
+      /* check (and mark) only horizontal b/w transitions */
+      // --- insert new box in list
+      box3 = (struct box *)malloc_box(NULL);
+      box3->x0=box3->x1=box3->x=x;
+      box3->y0=box3->y1=box3->y=y;
+      box3->num_frames=0;
+      box3->dots=0;
+      box3->num_boxes=1;
+      box3->num_subboxes=0;
+      box3->modifier='\0';
+      box3->num=job->res.numC;
+      box3->line=0;	// not used here
+      box3->m1=0; box3->m2=0; box3->m3=0; box3->m4=0;
+      box3->p=p;
+      box3->num_ac=0;   // for future use
+
+/*  frame, vectorize and mark only odd/even horizontal b/w transitions
+ *  args: box, x,y, cs, mark, diag={0,1}, ds={2,6}
+ *  ds   - start direction, 6=right of right border, 2=left of left border
+ *  ret  - 0=ok, -1=already marked, -2=max_num_frames_exceeded
+ *               -7=no border in direction ds 
+ *  ToDo: count errors and print out for debugging
+ */
+      rc=frame_vector(box3, x, y, cs, 1, 1, ds);
+      g_debug(fprintf(stderr,"\n# ... scan xy= %3d %3d rc= %2d", x, y, rc);)
+      if (rc<0) { free_box(box3); continue; }
+      if (box3->num_frames && !box3->num_frame_vectors[0])
+        fprintf(stderr,"\nERROR scan_boxes: no vector in frame (%d,%d)",x,y);
+
+      job->res.numC++;
+      job->res.sumX += box3->x1 - box3->x0 + 1;
+      job->res.sumY += box3->y1 - box3->y0 + 1;
+
+      box3->c=(((box3->y1-box3->y0+1)
+               *(box3->x1-box3->x0+1)>=MaxBox)? PICTURE : UNKNOWN);
+      list_app(&(job->res.boxlist), box3); 	// append to list
+      // ToDo: debug
+      // if (job->cfg.verbose && box3->y0==29) out_x(box3);
+  }
+  if(job->res.numC){
+    if (job->cfg.verbose)
+      fprintf(stderr," nC= %3d avD= %2d %2d\n",job->res.numC,
+               (job->res.sumX+job->res.numC/2)/job->res.numC,
+               (job->res.sumY+job->res.numC/2)/job->res.numC);
+  }
+  return job->res.numC;
+}
+
+/* compare ints for sorting.  Return -1, 0, or 1 according to
+   whether *vr < *vs, vr == *vs, or *vr > *vs */
+int 
+intcompare (const void *vr, const void *vs)
+{
+  int *r=(int *)vr;
+  int *s=(int *)vs;
+
+  if (*r < *s) return -1;
+  if (*r > *s) return 1;
+  return 0;
+}
+
+/*
+ * measure_pitch - detect monospaced font and measure the pitch
+ * measure overall pitch for difficult lines,
+ *  after that measure pitch per line 
+ * dists arrays are limited to 1024 elements to reduce
+ *  cpu usage for qsort on images with extreme high number of objects
+ * insert space if dist>=pitch in list_insert_spaces()
+ *  ToDo: ???
+ *   - min/max distance-matrix  a-a,a-b,a-c,a-d ... etc;  td,rd > ie,el,es
+ *   - OR measuring distance as min. pixel distance instead of box distance
+ *        especially useful for italic font!
+ *   - Kerning detection? minspace<=0 ???
+ *   - iterate minMono+maxMonoWidth and count fitting and misfitting pairs
+ * Lit:
+ *  http://en.wikibooks.org/wiki/LaTeX/Formatting
+ *         #The_Space_between_Words_and_Sentences
+ *   \frenchspacing == no extra space after periods (word vs. sentences)
+ *   \sloppypar     == some spaces between words may be to large
+ *   inter word space
+ *  http://en.wikipedia.org/wiki/Space_(punctuation)
+ *   Variable-width general-purpose space == 1/5-em to 1/3-em
+ *  http://en.wikipedia.org/wiki/Em_(typography) 
+ *   em = absolute maximum high,
+ *       median cap height=0.70em,
+ *       x-height=1ex=0.45..0.48..0.5em
+ *  http://en.wikipedia.org/wiki/En_(typography) = n-width=0.5em
+ *  http://pfaedit.sourceforge.net/glossary.html#overshoot
+ *   i: left + right side bearing (character specifique, may be negative: VA)
+ *  http://en.wikipedia.org/wiki/Typeface
+ *  http://en.wikipedia.org/wiki/Letter-spacing
+ *  http://en.wikipedia.org/wiki/Tracking_(typography) # Overlap VA
+ *  http://en.wikipedia.org/wiki/Kerning  # Overlap VA AT Tx etc.
+ *    similar blank 2D-area between pairs of characters
+ *     Helvetica: ry=+30 AV=-80 units?
+ * 
+ */
+void measure_pitch( job_t *job ){     /* word spacing */
+  int numdists=0, spc=0,              /* number of stored distances */
+      pitch_p=2, pdist, pdists[1024], /* proportional distances */
+      pitch_m=10, /* monospaced em width */
+      monospaced=1, l1, char_width_min=1023, char_width_max=0,
+      mono_em_min=0,    // maximum monospace char width + 1   2010-09-25
+      mono_em_max=2047, // minimum distance left side of two chars
+      d1l, d1r; // left-left and right-right distance of 2 chars
+  int d1, d2; // temporary vars, d1l + d1r sorted
+  struct box *box2, *pre1=NULL, *pre2=NULL;
+
+  if(job->cfg.verbose){ fprintf(stderr,"# check for word pitch"); }
+  for (l1=0; l1<job->res.lines.num; l1++)
+  { /* 0 means all lines */
+    if(job->cfg.verbose){ fprintf(stderr,"\n#  line %2d\n# ...",l1); }
+    numdists = 0;  /* clear distance lists */
+    monospaced=1; mono_em_min=0;  mono_em_max=2047; // reset, 2010-09-28
+    char_width_min=1023; char_width_max=0; // reset, 2010-09-28
+    for_each_data(&(job->res.boxlist)) {
+      box2 = (struct box *)list_get_current(&(job->res.boxlist));
+      if (l1>0 && box2->line!=l1) continue; /* ignore other lines */
+      /* ignore dots and pictures (min. font is 4x6) */
+      if (box2->y1 - box2->y0 + 1 < 4 || box2->c==PICTURE) pre2=pre1=NULL;
+      if (!pre1) { pre1=box2; continue; } /* we need a predecessor */
+      /* use gap for proportional fonts */
+      pdist = box2->x0 - pre1->x1 - 1; /* do not add 1, subtract 1 ! */
+      if (pdist<0) { // new line
+        pre2=NULL; pre1=box2; continue; }
+      if ((box2->x1 - box2->x0 + 1)
+       >2*(box2->y1 - box2->y0 + 1)) { // skip long object
+        continue; }
+      if ((pre1->x1 - pre1->x0 + 1)
+       >2*(pre1->y1 - pre1->y0 + 1)) { // skip long object
+        pre1=box2; continue; }
+      // JS-2010-09 sample spaces20100910.jpg 7 chars, fix bad auto space
+      if (char_width_min > box2->x1 - box2->x0 + 1)
+          char_width_min = box2->x1 - box2->x0 + 1;
+      if (box2->x1 - box2->x0 < 4*(pre1->x1 - pre1->x0)) // ~ big lines
+      if (char_width_max < box2->x1 - box2->x0 + 1)
+          char_width_max = box2->x1 - box2->x0 + 1;
+      // may cause problems if "_" is of width em (not em-1 like mwMW etc.)
+      if (mono_em_min < char_width_max + 1)
+          mono_em_min = char_width_max + 1; // minimum monospaced width
+
+      // will fail on monospaced fonts where chars are not centered
+      if (pre1) { // 2010-09-28
+        d1l = box2->x0 - pre1->x0; // left to left distance
+        d1r = box2->x1 - pre1->x1; // right to right distance
+        if (d1l > d1r) { d1=d1r; d2=d1l; } // thinner char on the right
+        else           { d1=d1l; d2=d1r; } // thicker char on the right
+        /* d1 < 2*width && d2 < 2*width, may fail for "IIIM" d2<2*max OK */
+        if (d1>0 && d1 < 2*char_width_max && d2 < 2*mono_em_max) {
+          if (mono_em_min<d1-1) mono_em_min = d1; }
+        if (d1>0) {
+          if (mono_em_max>d2+2) mono_em_max = d2; } // not best, shifted ()
+        // 2010-10-06 examples/ocr-b add -1 +2, bad for "()"        
+#if 1
+        if ((48 & job->cfg.verbose) == 48)
+        if (monospaced && l1)  // debugging until monospaced=0
+          fprintf(stderr," L%02d DBG1 x %3d %+4d %3d %+4d  d %3d %3d"
+            "  em %2d %2d  ex %2d\n# ...",
+            l1, pre1->x0, pre1->x1-pre1->x0+1, 
+                box2->x0, box2->x1-box2->x0+1, d1, d2,
+            mono_em_min, mono_em_max, char_width_max);
+#endif
+      }
+#if 1 // needed for correct spacing of last line of tmp08/0810CSchulze_crop
+      if (pre2) {
+        d1l = box2->x0 - pre2->x0; // left to left distance
+        d1r = box2->x1 - pre2->x1; // right to right distance
+        if (d1l > d1r) { d1=d1r; d2=d1l; } // thinner char on the right 
+        else           { d1=d1l; d2=d1r; } // thicker char on the right        
+        if (d1>0 && d1 < 3*char_width_max && d2 < 3*mono_em_max) {
+          if (2*mono_em_min<d1) mono_em_min = (d1+1)/2; }
+        if (d1>0) {
+          if (2*mono_em_max>d2) mono_em_max = (d2+1)/2; }
+#if 1
+        if ((48 & job->cfg.verbose) == 48)
+        if (monospaced && l1)  // debugging until monospaced=0
+          fprintf(stderr," L%02d DBG2 x %3d %+4d %3d %+4d  d %3d %3d"
+            "  em %2d %2d  ex %2d\n# ...",
+            l1, pre2->x0, pre2->x1-pre2->x0+1,
+                box2->x0, box2->x1-box2->x0+1, d1, d2,
+            mono_em_min, mono_em_max, char_width_max);
+#endif
+      }
+#endif
+
+// the upper part does good work, we do not need this stuff ... ???
+#if 0
+      // min distance between next neighbours of pre
+      if (pre2  &&  1 < box2->x0 - pre2->x1)
+      if (mono_em_max > box2->x0 - pre2->x1)
+          mono_em_max = box2->x0 - pre2->x1;
+      // ToDo: could be a problem for " ???
+      if (pre2)
+      if (pre1->x1 - pre1->x0 >= mono_em_min) // best max mono_dx
+      if (pre1->x1 - pre1->x0 == box2->x1 - box2->x0) // best max mono_dx
+      if (mono_em_max > box2->x0 - pre1->x0)
+          mono_em_max = box2->x0 - pre1->x0;
+      /* ToDo: better take 3 instead of 2 neighbours?, smallest font 4x6 */
+      /* tmp08/gocr0801_bad5.jpg was not mono, need 2 to 3 chars */
+      /* 2010-09-27 gives precise range! 16..22 to 16..17 */
+      /* ToDo: no 2 char variant? */
+      if (pre2  &&  1 < box2->x0 - pre2->x1)
+      if (box2->x0-pre1->x1+1 < mono_em_min) // no spc between char + pre1
+      if (pre1->x0-pre2->x1+1 < mono_em_min) // no spc between pre1 + pre2
+      {
+        if (3*mono_em_min <  box2->x1 - pre2->x0)
+              mono_em_min = (box2->x1 - pre2->x0 + 2)/3;
+      }  
+#endif
+//#  tmp09/oebb_teletext_836_0001_sw.png
+//#  line 4 12 - 12  pre2 134 142           181 190
+//#                         0   8            47  56
+//#                         0    12  24  36   48                   
+      // n=2: (n-1)*min < d1 <= (n  )*max  &&  (2*n+1)*max < (2*n+2)*min
+      //      (n  )*min < d2 <= (n+1)*max  &&  (2*n+2)*max < (2*n+3)*min 
+      if (monospaced && pre1) { // check 2 chars for non mono space within
+        d1l = box2->x0 - pre1->x0; // left to left distance  (do not + 1!)
+        d1r = box2->x1 - pre1->x1; // right to right distance
+        if (d1l > d1r) { d1=d1r; d2=d1l; } // thinner char on the right
+        else           { d1=d1l; d2=d1r; } // thicker char on the right
+        if ((box2->x0 - pre1->x1 <=   mono_em_min
+          && box2->x1 - pre1->x0 >  2*mono_em_max) // crossing 1 em border?
+         || (box2->x0 - pre1->x1 >    mono_em_max 
+          && box2->x0 - pre1->x1 <= 2*mono_em_min 
+          && box2->x1 - pre1->x0 >  3*mono_em_max)) { // crossing 2 em border?
+          monospaced = 0; // can not be monospaced in that case 2010-09-25
+         if (job->cfg.verbose)
+           fprintf(stderr, " L%02d mono:=0  %d - %d  pre1 %d %d  %d %d\n# ...",
+            l1, mono_em_min, mono_em_max,
+            pre1->x0, pre1->x1, box2->x0, box2->x1);
+        }
+      }
+      // n=3: (n-1)*min < d1 <= (n  )*max  &&  (2*n+1)*max < (2*n+2)*min
+      //      (n  )*min < d2 <= (n+1)*max  &&  (2*n+2)*max < (2*n+3)*min 
+      if (monospaced && pre2 && (2*2+2)*mono_em_max < (2*2+3)*mono_em_min)
+      { // check 2 chars for non mono space within
+        d1l = box2->x0 - pre2->x0; // left to left distance
+        d1r = box2->x1 - pre2->x1; // right to right distance
+        if (d1l > d1r) { d1=d1r; d2=d1l; } // thinner char on the right
+        else           { d1=d1l; d2=d1r; } // thicker char on the right
+        if ((box2->x0 - pre2->x1 >    mono_em_max
+          && box2->x0 - pre2->x1 <= 2*mono_em_min
+          && box2->x1 - pre2->x0 >  3*mono_em_max) // crossing 2 em border?
+         || (box2->x0 - pre2->x1 >  2*mono_em_max 
+          && box2->x0 - pre2->x1 <= 3*mono_em_min  // ?????? ToDo oebb
+          && box2->x1 - pre2->x0 >  4*mono_em_max)) { // crossing 3 em border?
+          monospaced = 0; // can not be monospaced in that case 2010-09-25
+         if (job->cfg.verbose)
+           fprintf(stderr, " L%02d mono:=0  %d - %d  pre2 %d %d  %d %d\n# ...",
+            l1, mono_em_min, mono_em_max,
+            pre2->x0, pre2->x1, box2->x0, box2->x1);
+        }
+      }
+      /* fonts are expected to be 6 to 60 pixels high, which is about
+         4 to 50 pixels wide.  We allow some extra margin.
+         space > 0 2010-09-27
+         ToDo: compare left and right gap (or additional nearest 4 gaps)
+         similar to mono space detection, check min distance
+          between upper, middle and lower rightmost vector of prev char and
+          leftmost vector of right char (hight is defined by the lower char)
+          (if overlapping chars are detected! WAV,Te,...)
+      */
+      if (0 < pdist && pdist < 140) { /* better mdist < 3*Xaverage ? */
+        // ignore extra wide spaces, tmp09/gocr_screen_capture* 2010-09-28
+        if (2*pdist<5*char_width_max)
+        /* two options for overflow: 1) ignore, 2) store randomly */
+        if (numdists<1024) {   /* we do ignore here */
+          pdists[numdists] = pdist;
+          numdists++;
+        }
+      }
+      pre2 = pre1; pre1 = box2;
+    } end_for_each(&(job->res.boxlist));
+    
+    if (job->cfg.verbose)
+      fprintf(stderr, " L%02d num_gaps= %2d x_width= %2d - %2d"
+        " mono_em= %2d - %2d  mono= %d",
+        l1, numdists, char_width_min, char_width_max,
+        mono_em_min, mono_em_max, monospaced);
+    if (numdists<8) {
+      if (job->cfg.verbose && l1==0) /* only for all lines */
+        fprintf(stderr," (WARNING num_gaps<8)");
+    }
+#if 1 /* debugging */
+    if ((job->cfg.verbose&(32+16))==48) {
+      int i;
+      fprintf(stderr,"\n# ...");   
+      for (i=0;i<numdists;i++) fprintf(stderr," %2d",pdists[i]);
+      fprintf(stderr," <- pdist[%d]\n# ...",l1);
+    }
+#endif
+    if (numdists>0) {
+      int i, diff, ni_min, max, best_p, ni;
+      /* aware: takes long time for big data sets */
+      /* dilute? (german: ausduennen?) */
+      qsort (pdists, numdists, sizeof (int), intcompare);
+      /* the new method, div0? */
+      best_p=4*numdists/5;
+      /* try to find better pitch for monospaced font (ok for prop) */
+      // tolerant to 090729num*  tmp09/barcodes090916_interleaved*
+      if (mono_em_min  > mono_em_max+mono_em_min/32+1
+       || mono_em_max>=2*mono_em_min)
+        monospaced = 0;
+      else
+        pitch_m=((mono_em_max<3*mono_em_min)?
+                 (mono_em_max+3*mono_em_min)/4:mono_em_min);
+      /* try to find better pitch for proportional font */
+      // the largest diff could be the best, if diff is always 1,
+      //  take the diff with the lowest weight
+      // JS-2010-09 add numdists<8 sample spaces20100908.jpg
+      // todo: search most offen biggest gapdiff (ignore big table gaps)
+      //   mean gapdiff? gap[n-1-i]-gap[0+i] until gapdiff=0, skip table gaps
+      // 2010-09-28 check until end of table, because old bad wide gaps are 
+      //     no more added to the table
+      for (ni=ni_min=1024,max=0,i=((numdists<8)?0:numdists/2+1);
+                                i<numdists;i++) {
+        if (pdists[i]<=char_width_min/3) continue; // JS-2010-09
+        if (pdists[i]> char_width_max*2) {
+         /* set 2nd best which is numdists as default */; break; } // JS-2010-27 table gaps
+        if (numdists<16)  // single word?
+        if (pdists[i]<=char_width_max/3) continue; // JS-2010-09
+        diff=pdists[i]-pdists[i-1];
+        if (diff>max) {
+          max=diff; best_p=i-1;
+          if ((job->cfg.verbose&(32+16))==48)
+            fprintf(stderr," L%02d best_p= %3d + maxdiff=%3d\n# ...",
+              l1, pdists[best_p], max);
+          if (max>3 &&                   3*pdists[i]>=4*pdists[i-1]) { break; }
+          if (max>1 && 3*i>numdists*2 && 3*pdists[i]>=4*pdists[i-1]) { break; }
+        }
+        if (diff) {
+          if (ni<ni_min) {
+            // do not try to divide one word per line
+            ni_min=ni; if (max<=1 && numdists>16) best_p=i-1;
+            if ((job->cfg.verbose&(32+16))==48)
+              fprintf(stderr," L%02d best_p=%3d ni_min=%3d\n# ...",
+                l1, pdists[best_p], ni_min);
+          }
+          ni=1;
+        } else ni++;
+      }
+      if (numdists<16 && max<=1 && ni_min>1) best_p=numdists-1; // one word 
+#if 1 /* debugging */
+      if ((job->cfg.verbose&(32+16))==48) {
+        // fprintf(stderr,"\n# ...");   
+        for (i=0;i<numdists;i++) fprintf(stderr," %2d",pdists[i]);
+        fprintf(stderr," <- pdist[%d] sorted\n# ...",l1);
+        fprintf(stderr," L%02d maxdiff=%d min_samediffs=%d", l1, max, ni_min);
+      }
+#endif
+      /* we measure spaces in two different ways (mono, prop) */
+      /* prop: gap between boxes,   mono: distance of middle */
+      if (best_p<numdists-1) pitch_p = ((pdists[best_p]+pdists[best_p+1])/2+1);
+      else                   pitch_p = (pdists[best_p]+1  );
+      if (numdists)
+        if (   pdists[numdists-1]*2 <= pdists[0]*3
+            || pdists[numdists-1]   <= pdists[0]+3) {
+        /* line is just a single word */
+          pitch_p = pdists[numdists-1]+10;
+        }
+      if (l1>0 && job->cfg.spc==0) {
+        job->res.lines.pitch[l1]=(monospaced?pitch_m:pitch_p);
+        job->res.lines.mono[l1]=monospaced;
+      }
+      if (job->cfg.verbose) {
+        fprintf(stderr,"\n# ...");
+        fprintf(stderr," L%02d mono: num=%3d min=%3d max=%3d pitch=%3d\n# ...",
+          l1, numdists, mono_em_min,mono_em_max,pitch_m);
+        fprintf(stderr," L%02d prop: num=%3d min=%3d max=%3d pitch=%3d @ %2d%%\n# ...",
+          l1, numdists, pdists[0],pdists[numdists-1],pitch_p,best_p*100/numdists);
+        fprintf(stderr," L%02d result: mono=%d  distance >= %d considered as space\n# ...",
+          l1, monospaced, job->res.lines.pitch[l1]);
+      }
+    } /* if (not) enough spaces */
+    if (l1==0) {  /* set default spaces to each line */
+      int l2;
+      spc = job->cfg.spc;
+      if (spc==0) /* set only if not set by option */
+        spc = ((monospaced)?pitch_m:pitch_p);
+      for  (l2=0; l2<job->res.lines.num; l2++ )
+        job->res.lines.pitch[l2]=spc;
+    }
+  }  /* each line */
+  if (job->cfg.spc==0)
+    job->cfg.spc = spc;
+  if (job->cfg.verbose)
+    fprintf(stderr," overall space width is %d %s\n",
+        spc, ((monospaced)?"monospaced":"proportional"));
+
+
+}
+
+/* ---- count subboxes (white holes within black area) --------
+ *  new: count boxes lying inside another box (usually holes, ex: "aeobdg")
+ *  needed for glue_boxes, dont joining textboxes, tables and other complex
+ *    objects
+ * ToDo: count only frames of invers spin? do we need sorted list here? -> no
+ */
+int count_subboxes( pix *pp ){
+  int ii=0, num_mini=0, num_same=0, cnt=0;
+  struct box *box2,*box4;
+  job_t *job=OCR_JOB; /* fixme */
+  progress_counter_t *pc = NULL;
+  if (job->cfg.verbose) { fprintf(stderr,"# count subboxes\n# ..."); }
+  
+  pc = open_progress(job->res.boxlist.n,"count_subboxes");
+  for_each_data(&(job->res.boxlist)) {
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    box2->num_subboxes=0;
+    progress(cnt++,pc);
+    if (   (box2->x1 - box2->x0)<2
+        || (box2->y1 - box2->y0)<2) continue; /* speedup for dotted bg */
+    // holes inside box2 char, aoebdqg, 0.41
+    for_each_data(&(job->res.boxlist)) {
+      box4=(struct box *)list_get_current(&(job->res.boxlist));
+      if (box4->y0 > box2->y1) break; // faster, but boxes need to be sorted
+      // ToDo: better use binary tree (above/below x) to find near boxes?
+      if (box4==box2) continue;
+      if( box4->x0==box2->x0 && box4->x1==box2->x1
+       && box4->y0==box2->y0 && box4->y1==box2->y1)
+         num_same++; /* erroneous!? */
+      if ( box4->x0 >= box2->x0  &&  box4->x1 <= box2->x1
+        && box4->y0 >= box2->y0  &&  box4->y1 <= box2->y1
+        && box4->num_subboxes==0 ) /* box4 inside box2? */
+      {
+        box2->num_subboxes++; ii++;
+        if ((box4->x1 - box4->x0 + 1)
+           *(box4->y1 - box4->y0 + 1)<17) num_mini++;
+      }
+    } end_for_each(&(job->res.boxlist));
+#if 0
+    if (cnt < 1000 && job->cfg.verbose)
+      fprintf(stderr," %4d box %4d %4d %+3d %+3d  subboxes %4d\n# ...",
+        cnt, box2->x0, box2->y0, box2->x1-box2->x0,
+                                 box2->y1-box2->y0, box2->num_subboxes);
+#endif
+  }   end_for_each(&(job->res.boxlist));
+  close_progress(pc);
+  if (job->cfg.verbose)
+    fprintf(stderr," %3d subboxes counted (mini=%d, same=%d) nC= %d\n",
+      ii, num_mini, num_same/2 /* counted twice */, cnt);
+  return 0;
+}
+
+/* ---- join holes to chars( before step1 ) v0.42  -----------------------
+   join boxes lying inside another box (usually holes, ex: "aeobdg46890")
+   Dont add dust to a char!
+   lines are not detected yet
+*/
+int glue_holes_inside_chars( pix *pp ){
+  int ii, cs, x0, y0, x1, y1, cnt=0,
+      glued_same=0, glued_holes=0;
+  struct box *box2, *box4;
+  job_t *job=OCR_JOB; /* fixme */
+  progress_counter_t *pc = NULL;
+  cs=job->cfg.cs;
+  {
+    count_subboxes( pp ); /* move to pgm2asc() later */
+    
+    pc = open_progress(job->res.boxlist.n,"glue_holes_inside_chars");
+    if (job->cfg.verbose)
+        fprintf(stderr,"# glue holes to chars nC= %d\n# ...",job->res.numC);
+    ii=0;
+    for_each_data(&(job->res.boxlist)) {
+      // get the smaller box which may be extended by bigger boxes around it
+      box2 = (struct box *)list_get_current(&(job->res.boxlist));
+      x0 = box2->x0;  x1 = box2->x1;
+      y0 = box2->y0;  y1 = box2->y1;
+      
+      progress(cnt++,pc);
+
+      // would it better than moving vectors to build a sub-box-tree?
+
+      // do not remove chars inside pictures (car plates on photos)
+      if( box2->c == PICTURE || box2->num_subboxes > 7) continue;
+
+      // holes inside char, aoebdqg, 0.41
+      // dont merge boxes which have subboxes by itself!
+      // search boxes inside box2 
+      // if (x1-x0+1>2 || y1-y0+1>2) /* skip tiny boxes, bad for 4x6 */
+      for_each_data(&(job->res.boxlist)) {
+	box4=(struct box *)list_get_current(&(job->res.boxlist));
+        if(box4!=box2 && box4->c != PICTURE )
+	{
+	  // ToDo: dont glue, if size differs by big factors (>16?)
+	  //  box4 is of same size or smaller
+	  //if ((job->cfg.verbose & 48)==48
+	  //   && abs(box4->x0-x0)<4 && abs(box4->y0-y0)<8)
+          //  { fprintf(stderr,"\n# DBG_glue");out_x(box2);out_x(box4); }
+          if (abs(box4->frame_vol[0])
+            >=abs(box2->frame_vol[0])/512) // 2010-10 bad invalid_ogv.jpg
+          if (   (    box4->x0==x0 && box4->x1==x1
+                   && box4->y0==y0 && box4->y1==y1 ) /* do not happen !? */
+              || (    box4->x0>=x0 && box4->x1<=x1
+                   && box4->y0>=y0 && box4->y1<=y1
+                   // 2010-09 subboxes==0 to subboxes<4 for 0 with dot in it
+                   && box4->num_subboxes<2 ) )  /* no or very small subboxes? */
+          {  // fkt melt(box2,box4)
+            // same box, if very small but hollow char (4x5 o)
+            if( box4->x0==x0 && box4->x1==x1
+             && box4->y0==y0 && box4->y1==y1) glued_same++; else glued_holes++;
+            // fprintf(stderr,"\n# DEBUG merge:");
+            // out_x(box2);  // small
+            // out_x(box4);  // big
+            if ((job->cfg.verbose & 7)==7) // LEV3
+              fprintf(stderr," join hole %4d %4d %+4d %+4d %+6d"
+                                     " + %4d %4d %+4d %+4d %+6d %d\n# ...",
+                x0, y0, x1-x0+1, y1-y0+1, box2->frame_vol[0],
+                box4->x0, box4->y0,
+                box4->x1-box4->x0+1, box4->y1-box4->y0+1, 
+                box4->frame_vol[0], glued_same);
+            if ((box4->x1-box4->x0+1)< 8*(x1-x0+1)
+             || (box4->y1-box4->y0+1)<12*(y1-y0+1)) // skip dust
+            merge_boxes( box2, box4 ); // add box4 to bigger box2
+  	    //if ((job->cfg.verbose & 48)==48)
+            //  { fprintf(stderr,"\n# DBG_glue_result");out_x(box2); }
+            x0 = box2->x0; x1 = box2->x1;
+            y0 = box2->y0; y1 = box2->y1;
+            job->res.numC--;  // dont count fragments as chars
+            ii++;	// count removed
+	    list_del(&(job->res.boxlist), box4); // remove box4
+	    free_box(box4);
+	    // now search another hole inside box2
+          }
+        }
+      } end_for_each(&(job->res.boxlist));
+
+    } end_for_each(&(job->res.boxlist)); 
+
+    if (job->cfg.verbose)
+      fprintf(stderr," joined: %3d holes, %3d same, nC= %d\n",
+        glued_holes, glued_same, job->res.numC);
+    close_progress(pc);
+  }
+  return 0;
+}
+
+
+/* ---- join broken chars ( before step1 ??? )  -----------------------
+    use this carefully, do not destroy previous detection ~fi, broken K=k' g 
+    join if boxes are near or diagonally connected 
+    other strategy: mark boxes for deleting and delete in extra loop at end
+    faster: check only next two following boxes because list is sorted!
+    ToDo: store m4 of upper line to m4_of_prev_line, and check that "-points are below
+    done: join boxes lying inside another box (usually holes, ex: "aeobdg")
+    Dont add dust to a char!
+    lines should be detected already (Test it for m1-m4 unknown)
+    ToDo: divide in glue_idots, glue_thin_chars etc. and optimize it
+*/
+int glue_broken_chars( job_t *job, pix *pp ){
+  int ii, y, cs, x0, y0, x1, y1, cnt=0,
+      num_frags=0, glued_frags=0, glued_hor=0,
+      do_join=0; /* 1..n means we have a reason to join two objects to one */
+  char *(join_reason)[5]={"no","\"A\"Uij\%","!?;\%","=:;","'',,"};
+  struct box *box2, *box4;
+  // job_t *job=OCR_JOB; /* fixme */
+  progress_counter_t *pc = NULL;
+  cs=job->cfg.cs;
+  {
+    count_subboxes( pp ); /* move to pgm2asc() later */
+    
+    pc = open_progress(job->res.boxlist.n,"glue_broken_chars");
+    if (job->cfg.verbose)
+        fprintf(stderr,"# glue broken chars nC= %d avX= %d\n# ...",
+          job->res.numC, job->res.avX);
+    ii=0;
+    for_each_data(&(job->res.boxlist)) {
+      // get the box which may be extended by boxes around it
+      box2 = (struct box *)list_get_current(&(job->res.boxlist));
+      x0 = box2->x0;  x1 = box2->x1;
+      y0 = box2->y0;  y1 = box2->y1;
+      progress(cnt++,pc);
+      do_join=0;
+      // vertical broken (g965T umlauts etc.)
+      // not: f,
+      // would it better than moving vectors to build a sub-box-tree?
+      // do not remove chars inside pictures (car plates on photos)
+      if (box2->c == PICTURE || box2->num_subboxes > 7) continue;
+      /* continue loop if box is below or above line = dust */
+      if (box2->m4>0 && y0>box2->m4) continue; /* dust outside ? */
+      if (box2->m1>0 && y0<box2->m1-(box2->m3-box2->m2)) continue;
+      /* ToDo:
+       *  - check that y0 is greater as m3 of the char/line above 
+       */
+      // --- variant 1 = ij-dots umlaut-dots :;= ---
+      // check small boxes (box2) whether they belong
+      //       to near same size or bigger boxes (box4)
+      if( 2*(y1-y0) < box2->m4 - box2->m1     // care for dots etc.
+       && (   2*y1<=(box2->m3+box2->m2)       // upper fragments
+           || 2*y0>=(box2->m3+box2->m2)) ) {  // lower fragments
+        struct box *box5=NULL, *box6=NULL;    // nearest and next nearest box
+        box4=NULL;
+        num_frags++;   /* count for debugging */
+        // get the [2nd] next x-nearest box in the same line
+        for_each_data(&(job->res.boxlist)) {
+  	  box4=(struct box *)list_get_current(&(job->res.boxlist));
+          if (box4 == box2  ||  box4->c == PICTURE) continue;
+          /* 0.42 speed up for background pixel pattern, box4 to small */
+          if ( box4->x1 - box4->x0 + 1 < x1-x0+1
+            && box4->y1 - box4->y0 + 1 < y1-y0+1 ) continue;
+          // have in mind that line number may be wrong for dust 
+          if (box4->line>=0 && box2->line>=0 && box4->line==box2->line)
+          {
+             if (!box5) box5=box4;
+             if ( abs(box4->x0 + box4->x1 - 2*box2->x0)
+                 <abs(box5->x0 + box5->x1 - 2*box2->x0))
+               { box6=box5; box5=box4; }
+  	  }
+        } end_for_each(&(job->res.boxlist));
+	box4=box5; // next nearest box within the same line
+      	if (box4) {
+          // do not glue "%^" in 0811qemu2.png 2010-09-28
+      	  if (box4->x1 - box4->x0 + 1 > job->res.avX / 2
+      	   && box2->x1 - box2->x0 + 1 > job->res.avX / 2
+      	   && (  box2->x0 > box4->x1
+      	      || box4->x0 > box2->x1)) continue;
+#if 0    /* set this to 1 for debugging of melting bugs */
+          if (job->cfg.verbose & 7) {
+            fprintf(stderr,"\n# next two boxes are candidates for joining");
+            out_x(box2);
+            out_x(box4); }
+#endif
+          if ( /* umlaut "a "o "u, ij; box2 is the small dot, box4 the body */
+                       4*y1 <= 3*box2->m2 + box2->m3 // y1=box2->y1, ocr-a %
+              && 4*box4->y1 >= 3*box2->m2 + box2->m3 // dont join 2 dots
+              &&   2* y1 < box4->y1 + box4->y0  // box2 above box4
+              &&   box4->x1 + job->res.avX/2 >= x0
+              &&   box4->x0 - job->res.avX/2 <= x1
+              && (y1 < box4->y0 || x0 < box4->x1) // dont melt "d'"
+              &&   3* (      y1 - box4->y0)
+                <= 2* (box4->y1 - box4->y0)  // too far away? dust!  
+              &&   8* (      x1 -       x0 + 1)
+                >=    (box4->x1 - box4->x0 + 1)  // dot must have minimum size  
+              &&  10* (      y1 -       y0 + 1)
+                >=    (box4->y1 - box4->y0 + 1)  // dot must have minimum size  
+            ) do_join=1;
+          if ( (!do_join) /* !?; box2 is the dot, box4 the body */
+              && 2*box4->x1>=x0+x1 	/* test if box4 is around box2 */
+              && 2*box4->x0<=2*x1 /* +x0+1 Jan00 */
+              && ( x1-x0 <= box4->x1-box4->x0+2 )
+              &&   2*y0>=box2->m2+box2->m3 
+              &&   4*y1>=box2->m2+3*box2->m3 
+              &&   4*(y1-y0)<box2->m4-box2->m1
+              &&   (8*box4->y1 < box4->m2+7*box4->m3
+                   || box4->m4-box4->m1<16) /* Jan00 */
+            ) do_join=2;
+          if ( (!do_join) /* =;: box2 is the upper box, box4 the lower box */
+              &&  2*box4->x1>=x0+x1 	/* test if box4 is around box2 */
+              && 2*box4->x0<=2*x1 /* +x0+1 */
+              && ( x1-x0  <=   box4->x1-box4->x0+4 )
+              && (  4*x0  <= 3*box4->x1+box4->x0 )
+              && (( box2->m2 && box4->m2
+                &&   y1< box2->m3
+                && 2*box4->y1 >    box4->m3+box4->m2  // can be bigger than m3
+                && 4*box4->y0 >= 3*box4->m2+box4->m3 
+                && 2*box2->y0 <    box2->m3+box2->m2
+                 )
+               || ( (!box2->m2) || (!box4->m2) )
+              )
+            ) do_join=3;
+          /* '' ,, tmp08/0811qemu2 2010-10-01 */
+          if (   abs(box2->y1 - box4->y1) <= (y1-y0)/8+1 // same y1
+              && abs(box2->y0 - box4->y0) <= (y1-y0)/8+1 // same y0
+              && abs((box4->x1 - box4->x0) - (x1-x0)) <= (x1-x0)/8+1 // same dx
+              && x1-x0 <= job->res.avX/2               // small width
+              && ( abs(box4->x0 - x1 - 1) <= job->res.avX/2  // small gap
+                || abs(x0 - box4->x1 - 1) <= job->res.avX/2) // ocr-b
+              && ( 4*y1 <= 3*box2->m2 +   box2->m3     // ''
+                || 4*y0 >= 2*box2->m2 + 2*box2->m3 )   // ,,
+            ) do_join=4;
+          if (do_join>0) {  // fkt melt(box2,box4)
+            if (job->cfg.verbose & 7) // space "( " for better " x"-searching
+              fprintf(stderr," join objects  %3d %3d %+4d %+4d"
+                                         " + %3d %3d %+4d %+4d %s\n# ...",
+                x0, y0, x1-x0+1, y1-y0+1, box4->x0, box4->y0,
+                box4->x1-box4->x0+1, box4->y1-box4->y0+1,join_reason[do_join]);
+            // fprintf(stderr,"\n# DEBUG merge:");  // d=7x34 @ (109,51) ???
+            // if (job->cfg.verbose & 4) out_x(box2);
+            // if (job->cfg.verbose & 4) out_x(box4);
+            merge_boxes( box2, box4 ); // add box4 to box2
+            x0 = box2->x0; x1 = box2->x1;
+            y0 = box2->y0; y1 = box2->y1;
+#if 0
+            if (job->cfg.verbose & 7) // 
+              fprintf(stderr," join objects  %3d %3d %+4d %+4d\n# ...",
+                x0, y0, x1-x0+1, y1-y0+1);
+#endif
+            // if (job->cfg.verbose & 4) out_x(box2);
+            // 2010-09-24 hmm, correct overall hight here, later set bad???
+            // job->res.numC--;  // dont count fragments as chars
+            ii++; glued_frags++; // remove
+            // output_list(job);
+	    list_del(&(job->res.boxlist), box4); /* ret&1: error-message ??? */
+            // output_list(job);
+	    free_box(box4);
+          }
+	}
+      }
+//  continue;
+
+      // horizontally broken w' K'
+      if(     2*y1  <   (box2->m3+box2->m2) )
+      if( 2*(y1-y0) <   (box2->m3+box2->m2) )	// fragment
+      for_each_data(&(job->res.boxlist)) {
+	box4=(struct box *)list_get_current(&(job->res.boxlist));
+        if (box4!=box2 && box4->c != PICTURE)
+	{
+          if( box4->line>=0 && box4->line==box2->line
+          && box4->x1>=x0-1 && box4->x1<x0  // do not glue 6-
+          && box4->x0+3*box4->x1<4*x0)
+          if( get_bw(x0  ,x0  ,y1,y1  ,pp,cs,1) == 1)
+          if( get_bw(x0-2,x0-1,y1,y1+2,pp,cs,1) == 1)
+          {  // fkt melt(box2,box4)
+            if (job->cfg.verbose & 7)
+              fprintf(stderr," join objects  %3d %3d %+4d %+4d"
+                                         " + %3d %3d %+4d %+4d w'K'\n# ...",
+                x0, y0, x1-x0+1, y1-y0+1, box4->x0, box4->y0,
+                box4->x1-box4->x0+1, box4->y1-box4->y0+1);
+            put(pp,x0,y1+1,~(128+64),0);
+            merge_boxes( box2, box4 );
+            x0 = box2->x0; x1 = box2->x1;
+            y0 = box2->y0; y1 = box2->y1;
+            job->res.numC--; ii++;	// remove
+            glued_hor++;
+	    list_del(&(job->res.boxlist), box4);
+	    free_box(box4);
+          }
+        }
+      } end_for_each(&(job->res.boxlist));
+
+      // horizontally broken n h	(h=l_)		v0.2.5 Jun00
+      if( abs(box2->m2-y0)<=(y1-y0)/8 )
+      if( abs(box2->m3-y1)<=(y1-y0)/8 )
+      if( num_cross(x0,         x1,(y0+  y1)/2,(y0+  y1)/2,pp,cs) == 1)
+      if( num_cross(x0,         x1,(y0+3*y1)/4,(y0+3*y1)/4,pp,cs) == 1)
+      if(    get_bw((3*x0+x1)/4,(3*x0+x1)/4,(3*y0+y1)/4,y1,pp,cs,1) == 0)
+      if(    get_bw(x0,(3*x0+x1)/4,(3*y0+y1)/4,(y0+3*y1)/4,pp,cs,1) == 0)
+      if(    get_bw(x0,         x0,         y0,(3*y0+y1)/4,pp,cs,1) == 1)
+      for_each_data(&(job->res.boxlist)) {
+	box4=(struct box *)list_get_current(&(job->res.boxlist));
+      	if (box4!=box2 && box4->c != PICTURE)
+	{
+          if( box4->line>=0 && box4->line==box2->line
+          && box4->x1>x0-3 && box4->x1-2<x0
+           && abs(box4->y1-box2->m3)<2)
+      	  {  // fkt melt(box2,box4)
+            if (job->cfg.verbose & 7)
+              fprintf(stderr," join objects %3d %3d %+4d %+4d"
+                                        " + %3d %3d %+4d %+4d nh\n# ...",
+                x0, y0, x1-x0+1, y1-y0+1, box4->x0, box4->y0,
+                box4->x1-box4->x0+1, box4->y1-box4->y0+1);
+      	    y=loop(pp,x0,y0,y1-y0,cs,0,DO);if(2*y>y1-y0) continue;
+            put(pp,x0-1,y0+y  ,~(128+64),0);
+            put(pp,x0-1,y0+y+1,~(128+64),0);
+            merge_boxes( box2, box4 );  // add box4 to box2
+            x0 = box2->x0; x1 = box2->x1;
+            y0 = box2->y0; y1 = box2->y1;
+            job->res.numC--; ii++;	// remove
+            glued_hor++;
+	    list_del(&(job->res.boxlist), box4);
+	    free_box(box4);
+          }
+      	}
+      } end_for_each(&(job->res.boxlist));
+    } end_for_each(&(job->res.boxlist)); 
+    if (job->cfg.verbose)
+      fprintf(stderr," joined: %3d fragments (found %3d), %3d rest, nC= %d\n",
+        glued_frags, num_frags, glued_hor, job->res.numC);
+    close_progress(pc);
+  }
+  return 0;
+}
+
+/*
+** this is a simple way to improve results on noisy images:
+** - find similar chars (build cluster of same chars)
+** - analyze clusters (could be used for generating unknown font-base)
+** - the quality of the result depends mainly on the distance function
+*/
+  // ---- analyse boxes, compare chars, compress picture ------------
+  // ToDo: - error-correction only on large chars! 
+int find_same_chars( pix *pp){
+  int i,k,d,cs,dist,n1,dx; struct box *box2,*box3,*box4,*box5;
+  pix p=(*pp);
+  job_t *job=OCR_JOB; /* fixme */
+  cs=job->cfg.cs;
+  {
+    if(job->cfg.verbose)fprintf(stderr,"# packing");
+    i = list_total(&(job->res.boxlist));
+    for_each_data(&(job->res.boxlist)) {
+      box4 = box2 = (struct box *)list_get_current(&(job->res.boxlist));
+      dist=1000;	// 100% maximum
+      dx = box2->x1 - box2->x0 + 1;
+
+      if(job->cfg.verbose)fprintf(stderr,"\r# packing %5d",i);
+      if( dx>3 )
+      for(box3=(struct box *)list_next(&(job->res.boxlist),box2);box3;
+	  box3=(struct box *)list_next(&(job->res.boxlist),box3)) {
+        if(box2->num!=box3->num){
+          int d=distance(&p,box2,&p,box3,cs);
+          if ( d<dist ) { dist=d; box4=box3; }	// best fit
+          if ( d<5 ){   // good limit = 5% ??? 
+            i--;n1=box3->num;		// set all num==box2.num to box2.num
+	    for_each_data(&(job->res.boxlist)) {
+	      box5=(struct box *)(struct box *)list_get_current(&(job->res.boxlist));
+	      if(box5!=box2)
+              if( box5->num==n1 ) box5->num=box2->num;
+	    } end_for_each(&(job->res.boxlist));
+          // out_x2(box2,box5);
+          // fprintf(stderr," dist=%d\n",d);
+          }
+      	}
+      }
+      // nearest dist to box2 has box4
+      //    out_b2(box2,box4);
+      //    fprintf(stderr," dist=%d\n",dist); 
+    } end_for_each(&(job->res.boxlist));
+    k=0;
+    if(job->cfg.verbose)fprintf(stderr," %d different chars",i);
+    for_each_data(&(job->res.boxlist)) {
+      struct box *box3,*box4;
+      int j,dist;
+      box2=(struct box *)list_get_current(&(job->res.boxlist));
+      for(box3=(struct box *)list_get_header(&(job->res.boxlist));
+          box3!=box2 && box3!=NULL;
+	  box3=(struct box *)list_next(&(job->res.boxlist), box3))
+        if(box3->num==box2->num)break;
+      if(box3!=box2 && box3!=NULL)continue;
+      i++;
+      // count number of same chars
+      dist=0;box4=box2;
+      
+      for(box3=box2,j=0;box3;
+          box3=(struct box *)list_next(&(job->res.boxlist), box3)) {
+	if(box3->num==box2->num){
+          j++;
+          d=distance(&p,box2,&p,box3,cs);
+          if ( d>dist ) { dist=d; box4=box3; }	// worst fit
+	}
+      }
+      if(job->cfg.verbose&8){
+        out_x2(box2,box4);
+        fprintf(stderr," no %d char %4d %5d times maxdist=%d\n",i,box2->num,j,dist);
+      }
+      // calculate mean-char (error-correction)
+      // ToDo: calculate maxdist in group 
+      k+=j;
+  //    if(j>1)
+  //    out_b(box1,NULL,0,0,0,0,cs);
+      if(job->cfg.verbose&8)
+      fprintf(stderr," no %d char %4d %5d times sum=%d\n",i,box2->num,j,k);   
+    } end_for_each(&(job->res.boxlist));
+    if(job->cfg.verbose)fprintf(stderr," ok\n");
+  }
+  return 0; 
+}
+
+/*
+** call the first engine for all boxes and set box->c=result;
+**
+*/
+int char_recognition( pix *pp, int mo){
+  int i,ii,ni,cs,x0,y0,x1,y1;
+  struct box *box2;
+  progress_counter_t *pc;
+  wchar_t cc;
+  job_t *job=OCR_JOB; /* fixme */
+  cs=job->cfg.cs;
+  // ---- analyse boxes, find chars ---------------------------------
+  if (job->cfg.verbose) 
+    fprintf(stderr,"# char recognition");
+  i=ii=ni=0;
+  for_each_data(&(job->res.boxlist)) { /* count boxes */
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    /* wew: isn't this just job->res.numC? */
+    /* js: The program is very complex. I am not sure anymore
+           wether numC is the number of boxes or the number of valid
+           characters.
+           Because its not time consuming I count the boxes here. */
+    if (box2->c==UNKNOWN)  i++;
+    if (box2->c==PICTURE) ii++;
+    ni++;
+  } end_for_each(&(job->res.boxlist));
+  if(job->cfg.verbose)
+    fprintf(stderr," unknown= %d picts= %d boxes= %d\n# ",i,ii,ni);
+  if (!ni) return 0;
+  i=ii=0;
+  pc = open_progress(ni,"char_recognition");
+  for_each_data(&(job->res.boxlist)) {
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    x0=box2->x0;x1=box2->x1;
+    y0=box2->y0;y1=box2->y1;	// box
+    cc=box2->c;
+    if (cc==PICTURE) continue;
+    
+    if ((mo&256)==0) { /* this case should be default (main engine) */
+      if(cc==UNKNOWN || box2->num_ac==0 || box2->wac[0]<job->cfg.certainty)
+        cc=whatletter(box2,cs   ,0);
+    }
+
+    if(mo&2) 
+      if(cc==UNKNOWN || box2->num_ac==0 || box2->wac[0]<job->cfg.certainty)
+	cc=ocr_db(box2, job);
+
+
+    // box2->c=cc; bad idea (May03 removed)
+    // set(box2,cc,95); ToDo: is that better? 
+
+    if(cc==UNKNOWN) 
+	i++;
+    ii++;
+
+    if(job->cfg.verbose&8) { 
+      fprintf(stderr,"\n# code= %04lx %c",(long)cc,(char)((cc<255)?cc:'_')); 
+      out_b(box2,pp,x0,y0,x1-x0+1,y1-y0+1,cs);
+    }
+    progress(ii,pc); /* ii = 0..ni */
+
+  } end_for_each(&(job->res.boxlist));
+  close_progress(pc);
+  if(job->cfg.verbose)fprintf(stderr," %d of %d chars unidentified\n",i,ii);
+  return 0;
+}
+
+
+/*
+** compare unknown with known chars,
+** very similar to the find_similar_char_function but here only to
+** improve the result
+*/
+int compare_unknown_with_known_chars(pix * pp, int mo) {
+  job_t *job=OCR_JOB; /* fixme */
+  int i, cs = job->cfg.cs, dist, d, ad, wac, ni, ii;
+  struct box *box2, *box3, *box4;
+  progress_counter_t *pc=NULL;
+  wchar_t bc;
+  i = ii = 0; // ---- -------------------------------
+  if (job->cfg.verbose)
+    fprintf(stderr, "# try to compare unknown with known chars !(mode&8)");
+  if (!(mo & 8))
+  {
+    ii=ni=0;
+    for_each_data(&(job->res.boxlist)) { ni++; } end_for_each(&(job->res.boxlist));
+    pc = open_progress(ni,"compare_chars");
+    for_each_data(&(job->res.boxlist)) {
+      box2 = (struct box *)list_get_current(&(job->res.boxlist)); ii++;
+      if (box2->c == UNKNOWN || (box2->num_ac>0 && box2->wac[0]<97))
+	if (box2->y1 - box2->y0 > 4 && box2->x1 - box2->x0 > 1) { // no dots!
+	  box4 = (struct box *)list_get_header(&(job->res.boxlist));;
+	  dist = 1000;		/* 100% maximum */
+	  bc = UNKNOWN;		/* best fit char */
+	  for_each_data(&(job->res.boxlist)) {
+	    box3 = (struct box *)list_get_current(&(job->res.boxlist));
+            wac=((box3->num_ac>0)?box3->wac[0]:100);	    
+	    if (box3 == box2 || box3->c == UNKNOWN
+                             || wac<job->cfg.certainty) continue;
+	    if (box2->y1 - box2->y0 < 5 || box2->x1 - box2->x0 < 3) continue;
+	    d = distance(pp, box2, pp, box3, cs);
+	    if (d < dist) {
+		dist = d;  bc = box3->c;  box4 = box3;
+	    }
+	  } end_for_each(&(job->res.boxlist));
+	  if (dist < 10) {
+            /* sureness can be maximal of box3 */
+	    if (box4->num_ac>0) ad = box4->wac[0];
+	    else                ad = 97;
+	    ad-=dist; if(ad<1) ad=1; 
+	    /* ToDo: ad should depend on ad of bestfit */
+	    setac(box2,(wchar_t)bc,ad);
+	    i++;
+	  }			// limit as option???
+	  //  => better max distance('e','e') ???
+	  if (dist < 50 && (job->cfg.verbose & 7)) {	// only for debugging
+	    fprintf(stderr,"\n#  L%02d xy= %4d %4d best fit was %04x=%c"
+	         " dist=%3d%% i=%d", box2->line, box2->x0, box2->y0,
+	         (int)bc, (char)((bc<128)?bc:'_'), dist, i);
+	    if (box4->num_ac>0) fprintf(stderr," w= %3d%%",box4->wac[0]);
+	    if ((job->cfg.verbose & 4) && dist < 10)
+	      out_x2(box2, box4);
+	  }
+	  progress(ii,pc);
+	}
+    } end_for_each(&(job->res.boxlist));
+    close_progress(pc);
+  }
+  if (job->cfg.verbose)
+    fprintf(stderr, " - found %d (nC=%d)\n", i, ii);
+  return 0;
+}
+
+/*
+// ---- divide overlapping chars which !strchr("_,.:;",c);
+// block-splitting (two ore three glued chars)
+// division if dots>0 does not work properly! ???
+//
+// what about glued "be"?
+// what about recursive division?
+// ToDo: mark divided boxes to give the engine a chance to 
+//       handle wrong divisions
+*/
+int  try_to_divide_boxes( pix *pp, int mo){
+  struct box *box2, boxa, boxb;
+  job_t *job=OCR_JOB; /* fixme */
+  int cs=job->cfg.cs, ad=100,
+      a2[8], ar, // certainty of each part, ar = product of all certainties
+      cbest;  // best certainty, skip search of certainty<cbest-1 for speed
+  wchar_t ci[8],  // split max. 8 chars 
+          s1[]={ UNKNOWN, '_', '.', ',', '\'', '!', ';', '?', ':', '-', 
+      '=', '(', ')', '/', '\\', '\0' };	// not accepted chars, \0-terminated!
+  int x0, x1, y0, y1,
+      xi[8+1]; // cutting positions
+  int i, ii, n1, dy, dx;
+  // pix p=(*pp); // remove!
+  if (job->cfg.verbose)
+    fprintf(stderr,"# try to divide unknown chars !(mode&16)");
+  if(!(mo&16))  // put this to the caller
+  for_each_data(&(job->res.boxlist)) {
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    // don't try to split simple structures (ex: 400x30 square)
+    if ((!box2->num_frames)
+       || box2->num_frame_vectors[ box2->num_frames-1 ]<9) continue; 
+    if((box2->c==UNKNOWN || (box2->num_ac && box2->wac[0]<job->cfg.certainty))
+      && box2->x1-box2->x0>5 && box2->y1-box2->y0>4){
+      x0=box2->x0; x1=box2->x1;
+      y0=box2->y0; y1=box2->y1;
+      ad=100;
+      cbest=0;
+      
+      /* get minimum vertical lines */
+      n1 = num_cross(x0,x1,(  y1+y0)/2,(  y1+y0)/2,pp,cs);
+      ii = num_cross(x0,x1,(3*y1+y0)/4,(3*y1+y0)/4,pp,cs); if (ii<n1) n1=ii;
+      if (box2->m2 && box2->m3 > box2->m2+2)
+      for (i=box2->m2+1;i<=box2->m3-1;i++) {
+        if (loop(pp,x0+1,i,x1-x0,cs,1,RI) > (x1-x0-2)) continue; // ll
+        ii = num_cross(x0,x1,i,i,pp,cs); if (ii<n1) n1=ii;
+      } if (n1<2) continue;  // seems to make no sense to divide
+      if (n1<4) ad=99*ad/100; // not to strong because m2+m3 could be wrong
+      if (n1<3) ad=99*ad/100;
+            
+      if( 2*y1 < box2->m3+box2->m4    /* baseline char ? */
+       && num_cross(x0,x1,y1-1,y1-1,pp,cs)==1  // -1 for slopes
+       && num_cross((x0+2*x1)/3,(x0+3*x1)/4,y0,y1,pp,cs)<3  // not exclude tz
+       && num_cross((3*x0+x1)/4,(2*x0+x1)/3,y0,y1,pp,cs)<3  // not exclude zl
+       && loop(pp,x0,y1-(y1-y0)/32,x1-x0,cs,0,RI)
+         +loop(pp,x1,y1-(y1-y0)/32,x1-x0,cs,0,LE) > (x1-x0+1)/2
+        ) continue; /* do not try on bvdo"o etc. */
+        
+      // one vertical line can not be two glued chars, lc?
+      if ( num_cross(x0,x1,(y1+y0)/2,(y1+y0)/2,pp,cs)<=1 ) continue;
+      {	// doublet = 2 letters
+        // char buf[4]="\0\0\0";      // 4th byte is string end == \0
+        // buf[0]=c1;                 // c1 is wchar_t! (0xbf00 to 0) failes
+        // buf[1]=c2;
+        char buf[64]="";      // end == \0
+        if (job->cfg.verbose&2){ 
+          fprintf(stderr, "\n#\n# divide box: %4d %4d %3d %3d\n",
+                  x0, y0, x1-x0+1, y1-y0+1);
+          if (job->cfg.verbose&4) out_x(box2);
+        }
+        // it would be better if testing is only if most right and left char
+        //   is has no horizontal gap (below m2) ex: be
+        i=0; // num splittet chars
+        xi[0]=x0; xi[1]=x0+3; xi[2]=x1;
+        for ( ; ; xi[i+1]++) { // x[i] .. x[i+1], slower? but better v0.42
+          /* break if x is to near to the right border */
+          if (xi[i+1]>x1-3) { if (i==0) break; i--; xi[i+2]=x1; continue; }
+          // ToDo: skip if not a local dy-min for speedup
+        { int ymin=y1, ymax=y0, bow=0, // min max at cutting point
+              max0=y0, max1=y0, // max y on left and right side
+              min0=y1, min1=y1; // min y on left and right side
+          for (dy=0,ii=0;ii<box2->num_frame_vectors[ 0 ];ii++) {
+            int pre=ii-1, next=(ii+1)%box2->num_frame_vectors[ 0 ];
+            if (pre<0) pre=box2->num_frame_vectors[ 0 ]-1;
+            // check if vector is inside box to cut
+            if (    box2->frame_vector[ii  ][0]<=xi[i  ]) continue;
+            if (    box2->frame_vector[ii  ][0]> xi[i+2]) continue;
+            // 2nd derivation of y(x)
+            if (abs(box2->frame_vector[ii  ][0]-xi[i+1])<2) {
+              dy= 2*box2->frame_vector[ii  ][1]
+                   -box2->frame_vector[next][1]
+                   -box2->frame_vector[pre ][1];
+              dx=   box2->frame_vector[next][0]
+                   -box2->frame_vector[pre ][0];
+              // rotate 180 degree if dx<0
+              if (((dx>0)?dy:-dy)<-abs(dx)/2) { bow=1; }
+            }
+            // its not the best if we think on glued fi fo etc.
+            if ((  box2->frame_vector[pre ][0]<=xi[i+1]
+                && box2->frame_vector[next][0]>=xi[i+1])
+             || (  box2->frame_vector[pre ][0]>=xi[i+1]
+                && box2->frame_vector[next][0]<=xi[i+1])) {
+             if (    box2->frame_vector[ii  ][1]>ymax)
+               ymax= box2->frame_vector[ii  ][1];
+             if (    box2->frame_vector[ii  ][1]<ymin)
+               ymin= box2->frame_vector[ii  ][1];
+            }
+            // min and max of left and right side
+            if (    box2->frame_vector[ii  ][1]>max0
+             &&     box2->frame_vector[ii  ][0]<=xi[i+1])
+               max0=box2->frame_vector[ii  ][1];
+            if (    box2->frame_vector[ii  ][1]>max1
+             &&     box2->frame_vector[ii  ][0]> xi[i+1])
+               max1=box2->frame_vector[ii  ][1];
+            if (    box2->frame_vector[ii  ][1]<min0
+             &&     box2->frame_vector[ii  ][0]<=xi[i+1])
+               min0=box2->frame_vector[ii  ][1];
+            if (    box2->frame_vector[ii  ][1]<min1
+             &&     box2->frame_vector[ii  ][0]> xi[i+1])
+               min1=box2->frame_vector[ii  ][1];
+          }
+          if(job->cfg.verbose&2)
+            fprintf(stderr,"\n# test if to split at x%d= %2d %2d %2d"
+                    " bow,(max-min)[i,0,1] %d %3d %3d %3d"
+                    , i, xi[i]-x0, xi[i+1]-x0, xi[i+2]-x0, bow, ymax-ymin, max0-min0, max1-min1);
+          /* skip if no local minimum at xi[i+1] or if its not thin enough */ 
+          // 2010-10-11 failes for ke on tmp08/gocr0801_bad5.jpg ToDo!!!
+//          if (bow==0 || 4*(ymax-ymin)>2*(y1-y0)) continue;
+          if (bow==0) continue;
+          // cuttet parts should have about the same height (max-min)
+          // we dont want to cut an 'n' in three parts!
+          if (2*(max0-min0+1)<(y1-y0+1)) continue;  // left  height
+          if (2*(max1-min1+1)<(y1-y0+1)) continue;  // right height
+          // ToDo: thickness on xi[i+1]?
+        }
+        // try to split successive right box if left box is recognised,
+        // else shift the splitting point further to the right border 
+        // removing ->dots if dot only above one char !!! ??? not implemented
+          if(job->cfg.verbose&2)
+            fprintf(stderr,"\n# try to split, newbox[%d].x= %2d ... %2d "
+                           "dy= %d ", i, xi[i]-x0, xi[i+1]-x0, dy);
+          boxa=*box2;	// copy contents, ToDo: reset ac-list (in cut_box?)
+          boxa.x=xi[i]; boxa.y=y0;        // obsolete? mark pixel, overlap?
+          boxa.x0=xi[i];boxa.x1=xi[i+1];  // new horizontal box range
+          cut_box(&boxa); boxa.num_ac=0;  // ToDo: add box2 as src argument?
+          // out_x(&boxa);
+          // get wchar + certainty
+          ci[i]=whatletter(&boxa,cs,0); a2[i]=testac(&boxa,ci[i]);
+          if(job->cfg.verbose&2)
+            fprintf(stderr,"\n#  certainty %d  limit= %d  cbest= %d ",
+                           a2[i], job->cfg.certainty, cbest);
+	  if (a2[i]<job->cfg.certainty || a2[i]<cbest-1
+	   || wcschr(s1,ci[i]) ) { continue; }  // dont split here
+
+          for (ar=ad,ii=0;ii<=i;ii++) {
+            ar=a2[ii]*ar/100; }  // multiply all probabilities
+	  if (ar<98*job->cfg.certainty/100 || ar<cbest) {
+            continue; } // dont go deeper, no longer string
+
+	  i++; if (i==8) break; // maximum splits
+	  if (i==4) break;  // at the moment its to slow to go further 
+	  if (i+1<8) xi[i+1]=x1;  // right border of next box
+	  if (i+2<8) xi[i+2]=x1;
+
+          if(job->cfg.verbose&2)
+            fprintf(stderr,"\n try end split [%d]=%d [%d]=%d ",
+                           i, xi[i]-x0, i+1, xi[i+1]-x0);
+          boxb=*box2;  // try rest if it has to be split again
+          boxb.x=xi[i]+1; boxb.y=y0;
+          boxb.x0=xi[i]+1;boxb.x1=xi[i+1];
+          cut_box(&boxb); boxb.num_ac=0;
+          ci[i]=whatletter(&boxb,cs,0); a2[i]=testac(&boxb,ci[i]);
+	  if (a2[i]<job->cfg.certainty || a2[i]<cbest-1
+	   || wcschr(s1,ci[i]) ) { xi[i+1]=xi[i]+2; continue; } // split rest
+	  // now we have everything splittet
+
+          if(job->cfg.verbose&2) {
+            fprintf(stderr,"\n split at/to: ");
+            for (ii=0;ii<=i;ii++)
+            fprintf(stderr,"  %2d %s (%3d)", xi[ii+1]-x0,
+              decode(ci[ii],ASCII), a2[ii]);
+            fprintf(stderr,"\n");
+          }
+	  // boxa..c changed!!! dots should be modified!!!
+          // Question: cut it into boxes v0.40 or set a string v0.41?
+          // new way of building a string v0.41 (can call setas multiple)
+          // usefull if compare unknown with known strings (except barcode?)
+          // ToDo: also create alternate variants? ex: I <-> l
+          for (buf[0]=0,ar=ad,ii=0;ii<=i;ii++) {
+            ar=a2[ii]*ar/100;  // multiply all probabilities
+            if (i>0 && ci[ii]=='n' && ci[ii-1]=='r') ar--; // m == rn
+            strncat(buf,decode(ci[ii],job->cfg.out_format),20);
+          }
+
+          if (ar>cbest) cbest=ar; // best (highest) certainty found
+          // reduce, but not if we cross certainty border
+          if (99*ar/100 > job->cfg.certainty) ar=99*ar/100;
+          if (job->cfg.verbose&2)
+            fprintf(stderr,"\n split result= %s (%3d) ",buf, ar);
+          setas(box2,buf,ar); // char *, does it disturb further splitting?
+          buf[0]=0;
+          i--; xi[i+2]=x1;
+        }
+      }
+    }
+  } end_for_each(&(job->res.boxlist));
+  if (job->cfg.verbose) fprintf(stderr,", numC %d\n",job->res.numC); 
+  return 0;
+}
+
+/*
+// ---- divide vertical glued boxes (ex: g above T);
+*/
+int  divide_vert_glued_boxes( pix *pp, int mo){
+  struct box *box2,*box3,*box4;
+  job_t *job=OCR_JOB; /* fixme */
+  int y0,y1,y,dy,flag_found,dx;
+  if(job->cfg.verbose)fprintf(stderr,"# divide vertical glued boxes");
+  for_each_data(&(job->res.boxlist)) {
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    if (box2->c != UNKNOWN) continue; /* dont try on pictures */
+    y0=box2->y0; y1=box2->y1; dy=y1-y0+1;
+    dx=4*(job->res.avX+box2->x1-box2->x0+1);     // we want to be sure to look at 4ex distance 
+    if ( dy>2*job->res.avY && dy<6*job->res.avY && box2->m1
+      && y0<=box2->m2+2 && y0>=box2->m1-2
+      && y1>=box2->m4+job->res.avY-2)
+    { // test if lower end fits one of the other lines?
+      box4=box2; flag_found=0;
+      for_each_data(&(job->res.boxlist)) {
+        box4 = (struct box *)list_get_current(&(job->res.boxlist));
+        if (box4->c != UNKNOWN) continue; /* dont try on pictures */
+        if (box4->x1<box2->x0-dx || box4->x0>box2->x1+dx) continue; // ignore far boxes
+        if (box4->line==box2->line  ) flag_found|=1;    // near char on same line
+        if (box4->line==box2->line+1) flag_found|=2;    // near char on next line
+        if (flag_found==3) break;                 // we have two vertical glued chars
+      } end_for_each(&(job->res.boxlist));
+      if (flag_found!=3) continue;         // do not divide big chars or special symbols
+      y=box2->m4;  // lower end of the next line
+      if(job->cfg.verbose&2){
+        fprintf(stderr,"\n# divide box below y=%4d",y-y0);
+        if(job->cfg.verbose&6)out_x(box2);
+      }
+      // --- insert box3 before box2
+      box3= (struct box *) malloc_box(box2);
+      box3->y1=y;
+      box2->y0=y+1; box2->line++; // m1..m4 should be corrected!
+      if (box4->line == box2->line){
+        box2->m1=box4->m1;        box2->m2=box4->m2;
+        box2->m3=box4->m3;        box2->m4=box4->m4;
+      }
+      box3->num=job->res.numC;
+      if (list_ins(&(job->res.boxlist), box2, box3)) {
+          fprintf(stderr,"ERROR list_ins\n"); };
+      job->res.numC++;
+    }
+  } end_for_each(&(job->res.boxlist));
+  if(job->cfg.verbose)fprintf(stderr,", numC %d\n",job->res.numC); 
+  return 0;
+}
+
+
+/* 
+   on some systems isupper(>255) cause a segmentation fault SIGSEGV
+   therefore this function
+   ToDo: should be replaced (?) by wctype if available on every system
+ */
+int wisupper(wchar_t cc){ return ((cc<128)?isupper(cc):0); }
+int wislower(wchar_t cc){ return ((cc<128)?islower(cc):0); }
+int wisalpha(wchar_t cc){ return ((cc<128)?isalpha(cc):0); }
+int wisdigit(wchar_t cc){ return ((cc<128)?isdigit(cc):0); }
+int wisspace(wchar_t cc){ return ((cc<128)?isspace(cc):0); }
+
+/* set box2->c to cc if cc is in the ac-list of box2, return 1 on success  */
+int setc(struct box *box2, wchar_t cc){
+  int ret=0, w1, w2;
+  w1=((box2->num_ac) ? box2->wac[0] : 0);  // weight of replaced char
+  w2=testac(box2,cc);
+  if (OCR_JOB->cfg.verbose) {
+    if (box2->num_ac<2)
+      fprintf(stderr, "\n#  change [%d] %s   %3d      to %s %3d at %4d %4d",
+       box2->num_ac, decode(box2->c,ASCII), w1,
+       decode(cc,ASCII), (100+w2+1)/2, box2->x0, box2->y0);
+    else
+      fprintf(stderr, "\n#  change [%d] %s %s %3d %3d  to %s %3d at %4d %4d",
+       box2->num_ac, decode(box2->c,ASCII),
+       decode(box2->tac[1],ASCII), box2->wac[0], box2->wac[1],
+       decode(cc,ASCII), (100+w2+1)/2, box2->x0, box2->y0);
+  }
+  if (w2) { if (box2->c!=cc) { ret=1; setac(box2,cc,(100+w2+1)/2); } }
+  // if(OCR_JOB->cfg.verbose & 4) out_x(box2);
+  // ToDo: modify per setac (shift ac)
+  return ret;
+}
+
+
+/* ---- proof difficult chars Il1 by context view ----
+  context: separator, number, vowel, nonvowel, upper case ????
+  could be also used to find unknown chars if the environment (nonumbers)
+    can be found in other places!
+  ToDo:
+   - box->tac[] as set of possible chars, ac set by engine, example:
+       ac="l/" (not "Il|/\" because serifs detected and slant>0)
+       correction only to one of the ac-set (alternative chars)!
+   - should be language-settable; Unicode compatible 
+   - box2->ad and wac should be changed? (not proper yet)
+ *  ------------- */
+int context_correction( job_t *job ) {
+ // const static char
+  char *l_vowel="aeiouy";
+    // *l_Vowel="AEIOU",chars if the environment (nonumbers)
+  char *l_nonvo = "bcdfghjklmnpqrstvwxzBCDFGHJKLMNPQRSTVWXZ";
+  int  hexdigits = 0, hexdivpos = 0; // "O0lI123456789ABCDEFabcdef:"
+  struct box *box3, *box2, *prev, *next, *pre2, *pre3, *pre4;
+  //  pix *pp = &(job->src.p);
+  int nc=0, ns=0; // num corrections
+  wchar_t last_double_quotation=0; // correction of different quotations "
+  pre4=pre3=pre2=prev=next=NULL;
+
+  if (job->cfg.verbose)
+    fprintf(stderr, "# context correction Il1 0O");
+
+  for_each_data(&(job->res.boxlist)) {
+    pre4=pre3; pre3 = pre2; pre2 = prev; // 2010-10-01 tmp08/080916_JL*_150
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    prev = (struct box *)list_get_cur_prev(&(job->res.boxlist));
+    next = (struct box *)list_get_cur_next(&(job->res.boxlist));
+    // ToDo: count last_upper, lower, digits, hexdigits
+    // 2010-10-10 hex-mode tmp08/gocr0801_bad5
+    if (strchr("O0lI123456789ABCDEFabcdef",box2->c)) hexdigits++;
+    else if (strchr(": ",box2->c) && prev && prev->c!=box2->c
+         && (hexdigits-hexdivpos==2 || hexdigits-hexdivpos==4))
+           hexdivpos=hexdigits;
+      else { hexdigits=0; hexdivpos=0; }
+    if (box2->c==' ' && prev && prev->c==' ') hexdigits=0;
+    if (box2->c==':' && pre3 && pre3->c!=':') hexdigits=0; // :89:AB:CD:
+    if (strchr("O0",box2->c) && hexdigits>5)  nc+=setc(box2,(wchar_t)'0');
+    if (strchr("l1",box2->c) && hexdigits>5)  nc+=setc(box2,(wchar_t)'1');
+    // 2010-10-01 sample tmp08/0811CSchulze_crop
+    if (box2->c==DOUBLE_LOW_9_QUOTATION_MARK) {
+      last_double_quotation = box2->tac[0];
+      fprintf(stderr,"\n#  ... found DOUBLE_LOW_9_QUOTATION_MARK");
+    }
+    if (box2->c==QUOTATION_MARK // 0x22 = ""
+      && last_double_quotation == DOUBLE_LOW_9_QUOTATION_MARK) {
+      last_double_quotation = 0;
+      box2->c = box2->tac[0] = DOUBLE_HIGH_REVERSED_9_QUOTATION_MARK;
+     if (job->cfg.verbose)
+       fprintf(stderr, "\n#  change [%d] %s   %3d      to %s %3d at %3d %3d",
+       box2->num_ac, "\"", box2->wac[0],
+       decode(box2->c,ASCII), box2->wac[0], box2->x0, box2->y0);
+    }
+    
+    if (           box2->c > 0xFF ) continue; // temporary UNICODE fix 1
+    if ((prev) && (prev->c > 0xFF)) continue; // temporary UNICODE fix 2
+    if ((next) && (next->c > 0xFF)) continue; // temporary UNICODE fix 3
+    if (box2->num_ac<2) continue; // no alternatives
+    if (box2->wac[0]==100 && box2->wac[1]<100) continue;
+    if (box2->num_ac && box2->tas[0]) continue; // buggy space_remove 0.42
+
+    /* check for Il1| which are general difficult to distinguish */
+    /* bbg: not very good. Should add some tests to check if is preceded by '.',
+     spelling, etc */
+    /* ToDo: only correct if not 100% sure (wac[i]<100)
+        and new char is in wat[] */
+    if (strchr("Il1|", box2->c) && next && prev) {
+//       if( strchr(" \n",prev->c)      // SPC 
+//        && strchr(" \n",next->c) ) box2->c='I'; else // bad idea! I have ...
+      if (wisalpha(next->c) && next->c!='i' && 
+          ( prev->c == '\n' || 
+	   ( prev->c == ' ' && pre2->c == '.' ) ) ) 
+        {  nc+=setc(box2,(wchar_t)'I'); }
+      else if (
+           (   box2->c!='1'             /* lnt => Int, but 1st */
+            && strchr(l_nonvo,next->c)
+            && strchr("\" \n",prev->c))
+          ||  (prev && ((!pre2) || wisupper(pre2->c) || strchr(" \n",pre2->c))
+            && wisupper(prev->c)
+            && box2->num_frame_vectors[0]==4
+            && box2->frame_vector[0][0]==box2->x0
+            && box2->frame_vector[1][0]==box2->x0
+            && box2->frame_vector[2][0]==box2->x1
+            && box2->frame_vector[3][0]==box2->x1
+            ))  // " DI*"
+        /* do not change he'll to he'Il! */
+        { nc+=setc(box2,(wchar_t)'I'); }  // set box2->c to 'I' if 'I' is in the ac-list
+      else if (strchr(l_vowel,next->c)) /* unusual? Ii Ie Ia Iy Iu */   
+          /*  && strchr("KkBbFfgGpP",prev->c)) */ /* kle Kla Kli */
+          {  nc+=setc(box2,(wchar_t)'l'); }
+      else if (wisupper(next->c)  // ToDo: check 6 neighbours for upper+spaces
+            && !strchr("O0I123456789",next->c)
+            && !strchr("O0I123456789",prev->c)) /* avoid lO => IO (10) */
+	{  nc+=setc(box2,(wchar_t)'I'); }
+      else if (prev && wislower(prev->c))
+	{  nc+=setc(box2,(wchar_t)'l'); }
+      else if (wisdigit(prev->c)
+            || wisdigit(next->c)
+            || (next && strchr(":-",next->c) && pre2 && pre2->c==next->c
+             && prev && strchr("0123456789ABCDabcd",prev->c)) // hex 2010-10
+            || (next->c=='O' && !wisalpha(prev->c)))  /* lO => 10 */
+	{  nc+=setc(box2,(wchar_t)'1'); }
+    }
+    // JS-2010-09 (ToDo: only if I is an alternate char!?)
+    if (strchr("Il|", box2->c) && next && !prev) { // first char?
+      if (wisalpha(next->c) && next->c!='i' && !strchr(l_vowel,next->c))
+	 {  nc+=setc(box2,(wchar_t)'I'); }
+      else if (wisupper(next->c)
+            && !strchr("O0I123456789",next->c)) /* avoid lO => IO (10) */
+	{  nc+=setc(box2,(wchar_t)'I'); }
+    }
+    
+    // ToDo: count width of all "0O" to decide between wide and narrow O's
+    /* check for O0 */
+    else if (strchr("O0", box2->c)) {
+      if (((!prev) ||  wisspace(prev->c))
+          && next && ( wisalpha(next->c))
+                  && (!wisdigit(next->c)))
+	  { nc+=setc(box2,(wchar_t)'O'); } // first letter ?
+      else if (prev && (!wisdigit(prev->c)) && wisupper(prev->c)
+            && next && (!wisdigit(next->c)) && wisupper(next->c))
+	  { nc+=setc(box2,(wchar_t)'O'); } // upper case word
+      // ! "  Otto"
+      else if (((!prev && next && wisdigit(next->c))
+             || (prev && wisdigit(prev->c))
+             || (prev && pre2 && wisspace(prev->c)) // tmp10/barcode_code128_10*
+             || (prev && strchr("/:.,-+O", prev->c)) // tmp09/barcode090916*
+             || (prev && strchr("abcdABCD", prev->c)
+              && pre2 && next && pre2->c==next->c
+                      && strchr("/:- ", pre2->c)) // hex?
+             )
+            && ((next && wisdigit(next->c))
+             || (next && strchr("\n /:.,-+O",next->c))   // 2010-10-01
+             || ((next && strchr("abcdABCD", next->c)
+               && pre2 && wisdigit(pre2->c)
+               && prev && strchr("/:- ", prev->c))) // hex?
+             || !next)
+            && (!pre2
+             || (pre2 && strchr("/:.,-+O",pre2->c))   // 2010-10-01
+             || (pre2 && wisdigit(pre2->c))   // 2010-10-01
+             || (pre2 && wisspace(pre2->c)))  // 2010-10-01
+            && (  (prev && wisdigit(prev->c))
+               || (pre2 && wisdigit(pre2->c))
+               || (pre3 && wisdigit(pre3->c))
+               || (next && wisdigit(next->c))) // one near digit required
+            ) 
+	{ nc+=setc(box2,(wchar_t)'0'); }
+    }
+
+    /* check for 5S */
+    else if (strchr("5S", box2->c) && next && prev) {
+      if (wisspace(prev->c) && wisalpha(next->c)) /* initial letter */
+	{ nc+=setc(box2,(wchar_t)'S'); }
+      else if (wisalpha(prev->c) && wisalpha(next->c)
+                                 && wisupper(next->c)) /* word in upper case */
+	{ nc+=setc(box2,(wchar_t)'S'); }
+      else if (wisdigit(prev->c) || wisdigit(next->c))
+	{ nc+=setc(box2,(wchar_t)'5'); }
+    }
+
+    /* was a space not found? xXx => x Xx ??? */
+    if (wisupper(box2->c) && next && prev) {
+      if (wislower(prev->c) && wislower(next->c)
+	  && 2 * (box2->x0 - prev->x1) > 3 * (next->x0 - box2->x1)) {
+	struct box *box3 = malloc_box((struct box *) NULL);
+	box3->x0 = prev->x1 + 2;
+	box3->x1 = box2->x0 - 2;
+	box3->y0 = box2->y0;
+	box3->y1 = box2->y1;
+	box3->x = box2->x0 - 1;
+	box3->y = box2->y0;
+	box3->dots = 0;
+	box3->num_boxes = 0;
+	box3->num_subboxes = 0;
+	box3->c = ' ';
+	box3->modifier = 0;
+	setac(box3,' ',99); /* ToDo: weight depends from distance */
+	box3->num = -1;
+	box3->line = prev->line;
+	box3->m1 = box3->m2 = box3->m3 = box3->m4 = 0;
+	box3->p = &(job->src.p);
+	list_ins(&(job->res.boxlist), box2, box3);
+      }
+    }
+    
+    /* a space before punctuation? but not " ./file" */
+    if ( prev && next)
+    if (prev->c == ' ' && strchr(" \n"    , next->c)
+                       && strchr(".,;:!?)", box2->c))
+      if (prev->x1 - prev->x0 < 2 * job->res.avX) {	// carefully on tables
+	box3 = prev;
+	if ( !list_del(&(job->res.boxlist), box3) ) free_box(box3);	
+        prev = (struct box *)list_get_cur_prev(&(job->res.boxlist));
+        ns++;
+      }
+
+    /* \'\' to \" */
+    if ( prev )
+    if ( (prev->c == '`' || prev->c == '\'')
+      && (box2->c == '`' || box2->c == '\'') )
+      if (prev->x1 - box2->x0 < job->res.avX) { // carefully on tables
+        box2->c='\"';
+	box3 = prev;
+	list_del(&(job->res.boxlist), box3);
+	free_box(box3);
+      }
+  } end_for_each(&(job->res.boxlist));
+  if (job->cfg.verbose)
+    fprintf(stderr, " num_corrected= %d removed_spaces= %d\n", nc, ns);
+  return 0;
+}
+
+
+/* ---- insert spaces ----
+ *  depends strongly from the outcome of measure_pitch()
+ * ------------------------ */
+int list_insert_spaces( pix *pp, job_t *job ) { 
+  int i=0, j1, j2, i1, maxline=-1, dy=0, num_nl=0, num_spc=0, min_x0=1023;
+  char cc;
+  struct box *box2, *box3=NULL, *box4=NULL;
+
+  // measure mean line height
+  for(i1=1;i1<job->res.lines.num;i1++) {
+    dy+=job->res.lines.m4[i1]-job->res.lines.m1[i1]+1;
+    if (min_x0>job->res.lines.x0[i1])
+        min_x0=job->res.lines.x0[i1];  // 2010-09-30
+  } if (job->res.lines.num>1) dy/=(job->res.lines.num-1);
+  i=0; j2=0;
+  for(i1=1;i1<job->res.lines.num;i1++) {
+    j1=job->res.lines.m4[i1]-job->res.lines.m1[i1]+1;
+    if (j1>dy*120/100 || j1<dy*80/100) continue; // only most frequently
+    j2+=j1; i++;
+  } if (i>0 && j2/i>7) dy=j2/i;
+  if( job->cfg.verbose&1 )
+    fprintf(stderr,"# insert space between words (dy=%d) ...",dy);
+  if (!dy) dy=(job->res.avY)*110/100+1;
+
+  if (min_x0 < 4) min_x0 = 0; // tmp09/oebb_teletext* monospaced first gap
+  // ToDo: rewrite, replace cc by num_spc + num_nl
+  i=0;
+  for_each_data(&(job->res.boxlist)) {
+    int thispitch=0, thismono=0, pdist=0; // spacing paras per line
+    box2 =(struct box *)list_get_current(&(job->res.boxlist));
+    cc=0; num_nl=0; num_spc=0;
+    box3 = (struct box *)list_prev(&(job->res.boxlist), box2);
+    if (box2->line > maxline) {  // new line, lines and chars must be sorted!
+      int ydist=0, ypitch=0;
+      if (maxline>=0) {
+        // num_nl = 1; // ToDo: allow multiple newlines
+        if (box2->line>1)
+         ydist = job->res.lines.m1[ box2->line   ]
+                -job->res.lines.m1[ box2->line-1 ]; // 2010-09-26
+        ypitch = job->res.lines.m4[ box2->line ]
+                -job->res.lines.m1[ box2->line ];
+        if (ypitch>4) num_nl = ydist / (2*ypitch); // ToDo: improve it!
+        if (!num_nl) num_nl=1;
+      }
+      maxline=box2->line;
+    }
+    if (box2->line==maxline) {  // lines and chars must be sorted!
+      thispitch = job->res.lines.pitch[box2->line];
+      thismono  = job->res.lines.mono[ box2->line];
+      if (box3) pdist  = box2->x0 - box3->x1 - 1; // 2010-09-26
+      if (pdist < 0) pdist = 0; // overlap like proportional: "VA" 
+      if (num_nl || !box3)
+        pdist  = box2->x0 - min_x0; // first char of new line
+      // if (pdist >= thispitch) cc=' '; // 2010-09-24 ???
+      if (thismono) num_spc = pdist / thispitch;
+      else          num_spc = pdist*2 / (3*job->res.avX); // ToDo: use 1em!
+      if (pdist>=thispitch && !num_spc) num_spc = 1; // proportional font
+      // ToDo: multi spaces for proportional font
+    }
+
+#if 0
+    if ((job->cfg.verbose&48)==48)
+      fprintf(stderr,"\n# DBG L%02d %d mono=%d  %d pitch= %2d"
+        " pdist= %2d nl %d spc %d", maxline, box2->line, thismono,
+        job->res.lines.mono[ box2->line], thispitch, pdist, num_nl, num_spc);
+#endif
+
+    // call this multiple times
+    for (i1=0;i1<num_nl+num_spc;i1++) {
+      int mdist=0;
+      box4=(struct box *)list_prev(&(job->res.boxlist), box2);
+      if (box4) mdist  = box2->x0 - box4->x1 + 1; // 2010-09
+      else      mdist  = 0;
+      if (mdist<0) mdist=0;
+      box3=(struct box *)malloc_box(NULL);
+      box3->x0=box2->x0-2+((num_spc)?-mdist+ i1   *mdist/num_spc:0);
+      box3->x1=box2->x0-2+((num_spc)?-mdist+(i1+1)*mdist/num_spc:0);
+      box3->y0=box2->y0;
+      box3->y1=box2->y1;
+      if (i1>=num_nl &&  box4)
+        box3->x0 = box4->x1+2+((num_spc)?i1*mdist/num_spc:0);
+      if (i1< num_nl || !box4)
+        box3->x0 = job->res.lines.x0[box2->line];
+      if (i1< num_nl && box4){
+        box3->y0=box4->y1;	// better use lines.y1[box2->pre] ???
+        box3->y1=box2->y0;
+      }
+      box3->x = box3->x0; // 2010-09
+      box3->y = box2->y0;
+      box3->dots = 0;
+      box3->c = cc = ((i1<num_nl)?'\n':' ');
+      box3->num_boxes = 0;
+      box3->num_subboxes = 0;
+      box3->modifier = '\0';
+      box3->num=-1;        box3->line=box2->line;
+      box3->m1=box2->m1;   box3->m2=box2->m2;
+      box3->m3=box2->m3;   box3->m4=box2->m4;
+      box3->p=pp;
+      setac(box3,cc,100);   /* ToDo: weight depends from distance */
+      list_ins(&(job->res.boxlist),box2,box3); // insert box3 before box2
+      if( job->cfg.verbose&1 ) {
+        fprintf(stderr,"\n# insert space &%d; at %4d %4d box= %p"
+          " mono %d dx %2d pdx,mdx %2d %2d",
+          (int)box3->c, box3->x0, box3->y0, (void*)box3,
+          thismono, thispitch, pdist, mdist);
+        /* out_x(box3); */
+      }
+      i++;
+    }
+  } end_for_each(&(job->res.boxlist));
+  if( job->cfg.verbose&1 ) fprintf(stderr,"\n# ... found %d spaces\n",i);
+  return 0;
+}
+
+
+/*
+   add infos where the box is positioned to the box
+   this is useful for better recognition
+*/
+int  add_line_info( job_t *job /* , List *boxlist2 */){
+  struct tlines *lines = &job->res.lines;
+  struct box *box2;
+  int i,xx,m1,m2,m3,m4,num_line_members=0,num_rest=0;
+  if (job->cfg.verbose&1) fprintf(stderr,"# add_line_info to boxes ...");
+  for_each_data(&(job->res.boxlist)) {
+    box2 =(struct box *)list_get_current(&(job->res.boxlist));
+    for (i=1;i<job->res.lines.num;i++) /* line 0 is a place holder */
+    { // add rotated image correction dy(x)
+      if (lines->dx) xx=lines->dy*((box2->x1+box2->x0)/2)/lines->dx;
+                else xx=0;
+      m1= lines->m1[i]+xx;
+      m2= lines->m2[i]+xx;
+      m3= lines->m3[i]+xx;
+      m4= lines->m4[i]+xx;
+      // fprintf(stderr," test line %d m1=%d %d %d %d\n",i,m1,m2,m3,m4);
+      if (m4-m1==0) continue; /* no text line (line==0) */
+#if 0
+      if( box2->y1+2*job->res.avY >= m1
+       && box2->y0-2*job->res.avY <= m4 ) /* not to far away */
+#endif
+      /* give also a comma or dot behind the line a chance */
+      if ( box2->x0 >= lines->x0[i] 
+        && box2->x1 <= lines->x1[i]+job->res.avX )
+      if ( box2->y0 <= m4 + 2*job->res.avY    // 2010-10-01+09 0811qemu2
+        && box2->y1 >= m1 - job->res.avY/2  // give "a "o ... a chance
+        && box2->y1 <= m4 + 2*job->res.avY )  // 2010-10-09 ocr-b-'_'
+      if ( box2->m2==0  // already put to a line? check y-distance
+        // || abs(box2->y0 - box2->m2) > abs(box2->y0 - m2)
+        // ToDo: m3 vs m2
+        ||  (abs(box2->y1 - box2->m3) > abs(box2->y1 - m1)
+          && box2->y0 > box2->m4)
+       )
+      { /* found nearest line */
+        box2->m1= m1;
+        box2->m2= m2;
+        box2->m3= m3;
+        box2->m4= m4;
+        box2->line= i;
+      }
+    } // i=1..lines
+    if( box2->y1+2 < box2->m1
+     || box2->y0   < box2->m1 - (box2->m3-box2->m1)/2
+     || box2->y0-2 > box2->m4 + (box2->m3-box2->m2)/2 // bad m4 + ,._ ocr-b
+     || box2->y1   > box2->m3 + (box2->m3-box2->m1)
+     ) /* to far away */
+    {  /* reset */
+        box2->m1= 0;
+        box2->m2= 0;
+        box2->m3= 0;
+        box2->m4= 0;
+        box2->line= 0;
+        num_rest++;
+    } else num_line_members++;
+  } end_for_each(&(job->res.boxlist));
+  if (job->cfg.verbose&1)
+    fprintf(stderr," done, num_line_chars=%d rest=%d\n",
+            num_line_members, num_rest);
+  return 0;
+}
+
+
+/*
+ *  bring the boxes in right order
+ *  add_line_info must be executed first!
+ */
+int sort_box_func (const void *a, const void *b) {
+  struct box *boxa, *boxb;
+
+  boxa = (struct box *)a;
+  boxb = (struct box *)b;
+
+  if ( ( boxb->line < boxa->line ) ||
+       ( boxb->line == boxa->line && boxb->x0 < boxa->x0 ) )
+    return 1;
+  return -1;
+}    
+
+// -------------------------------------------------------------
+// ------             use this for entry from other programs 
+// include pnm.h pgm2asc.h 
+// -------------------------------------------------------------
+// entry point for gocr.c or if it is used as lib
+// better name is call_ocr ???
+// jb: OLD COMMENT: not removed due to set_options_* ()
+// args after pix *pp should be removed and new functions
+//   set_option_mode(int mode), set_option_spacewidth() .... etc.
+//   should be used instead, before calling pgm2asc(pix *pp)
+//   ! change if you can ! - used by X11 frontend
+int pgm2asc(job_t *job)
+{
+  pix *pp;
+  progress_counter_t *pc;
+  static int multi_image_count=0;  /* number of image within multi-image */
+  int orig_cs=0; 
+  
+  if (!multi_image_count) orig_cs = job->cfg.cs; /* save for multi-images */
+  
+  multi_image_count++;
+
+  assert(job);
+  /* FIXME jb: remove pp */
+  pp = &(job->src.p);
+
+  pc = open_progress(100,"pgm2asc_main");
+  progress(0,pc); /* start progress output 0% 0% */
+#if 0 /* dont vast memory */
+  /* FIXME jb: malloc */
+  if ( job->cfg.verbose & 32 ) { 
+    // generate 2nd imagebuffer for debugging output
+    job->tmp.ppo.p = (unsigned char *)malloc(job->src.p.y * job->src.p.x); 	
+    // buffer
+    assert(job->tmp.ppo.p);
+    copybox(&job->src.p,
+            0, 0, job->src.p.x, job->src.p.y,
+            &job->tmp.ppo,
+            job->src.p.x * job->src.p.y);
+  }
+#else
+  job->tmp.ppo=job->src.p; /* temporarely, removed later */
+#endif
+  // if (job->cfg.verbose&32) debug_img("out000.ppm",job,0);
+
+  /* ----- count colors ------ create histogram -------
+     - this should be used to create a upper and lower limit for cs
+     - cs is the optimum gray value between cs_min and cs_max
+     - also inverse scans could be detected here later */
+  if (orig_cs==0)
+    job->cfg.cs=otsu( pp->p,pp->y,pp->x,0,0,pp->x,pp->y,job->cfg.verbose & 1);
+  else  // dont set cs, output stats + do inversion if needed 2010-10-07
+    otsu( pp->p,pp->y,pp->x,0,0,pp->x,pp->y,job->cfg.verbose & 1);
+//  if (job->cfg.verbose&32) debug_img("out001.ppm",job,0);
+  /* renormalize the image and set the normalized threshold value */
+  job->cfg.cs=thresholding( pp->p,pp->y,pp->x,0,0,pp->x,pp->y, job->cfg.cs );
+  if( job->cfg.verbose ) 
+    fprintf(stderr, "# thresholding new_threshold= %d\n", job->cfg.cs);
+//  if (job->cfg.verbose&32) debug_img("out002.ppm",job,0);
+
+  progress(5,pc); /* progress is only estimated */
+
+  
+  /* this is first step for reorganize the PG
+     ---- look for letters, put rectangular frames around letters
+     letter = connected points near color F
+     should be used by dust removing (faster) and line detection!
+     ---- 0..cs = black letters, last change = Mai99 */
+  
+  progress(8,pc); /* progress is only estimated */
+
+//  if (job->cfg.verbose&32) debug_img("out008.ppm",job,8);
+  scan_boxes( job, pp );
+  if ( !job->res.numC ){ 
+    fprintf( stderr,"# no boxes found - stopped\n" );
+    if(job->cfg.verbose&32) debug_img("out01",job,8);
+    /***** should free stuff, etc) */
+    return(1);
+  }
+  // tmp10/bug100818a.pgm creates artefacts on image
+//  if (job->cfg.verbose&32) debug_img("out00",job,4+8);
+
+  progress(10,pc); /* progress is only estimated */
+  // if(job->cfg.verbose&32) debug_img("out01",job,4+8);
+  // output_list(job);  // for debugging 
+  // ToDo: matrix printer preprocessing
+
+  remove_dust( job ); /* from the &(job->res.boxlist)! */
+// if(job->cfg.verbose&32) debug_img("out02",job,4+8);
+// output_list(job);  // for debugging 
+#if 0 // ToDo 2010-10-15 destroys QR-barcodes
+  smooth_borders( job ); /* only for big chars */
+#endif
+  progress(12,pc); /* progress is only estimated */
+// if(job->cfg.verbose&32) debug_img("out03",job,4+8);
+// output_list(job);  // for debugging 
+
+  detect_barcode( job );  /* mark barcode */
+// if(job->cfg.verbose&32) debug_img("out04",job,4+8);
+// output_list(job);  // for debugging 
+
+  detect_pictures( job ); /* mark pictures */
+//  if(job->cfg.verbose&32) debug_img("out05",job,4+8);
+// output_list(job);  // for debugging 
+
+  remove_pictures( job ); /* do this as early as possible, before layout */
+//  if(job->cfg.verbose&32) debug_img("out06",job,4+8);
+// output_list(job);  // for debugging
+
+  glue_holes_inside_chars( pp ); /* including count subboxes (holes)  */
+
+  detect_rotation_angle( job );
+
+#if 1 		/* Rotate the whole picture! move boxes */
+  if( job->res.lines.dy!=0 ){  // move down lowest first, move up highest first
+    // in work! ??? (at end set dy=0) think on ppo!
+  }
+#endif
+  detect_text_lines( pp, job->cfg.mode ); /* detect and mark job->tmp.ppo */
+// if(job->cfg.verbose&32) debug_img("out07",job,4+8);
+  progress(20,pc); /* progress is only estimated */
+
+  add_line_info( job /* , &(job->res.boxlist) */);
+  if (job->cfg.verbose&32) debug_img("out10",job,4+8);
+
+  divide_vert_glued_boxes( pp, job->cfg.mode); /* after add_line_info, before list_sort! */
+//  if(job->cfg.verbose&32) debug_img("out11",job,0);
+
+  remove_melted_serifs( job, pp ); /* make some corrections on pixmap */
+  /* list_ins seems to sort in the boxes on the wrong place ??? */
+//  if(job->cfg.verbose&32) debug_img("out12",job,4+8);
+
+  glue_broken_chars( job, pp ); /* 2nd glue */
+//  if(job->cfg.verbose&32) debug_img("out14",job,4+8);
+// 2010-09-24 overall box size is correct here, but later broken
+
+  remove_rest_of_dust( job );
+//  if(job->cfg.verbose&32) debug_img("out15",job,4+8);
+
+  /* better sort after dust is removed (slow for lot of pixels) */ 
+  list_sort(&(job->res.boxlist), sort_box_func);
+
+  measure_pitch( job );
+
+  if(job->cfg.mode&64) find_same_chars( pp );
+  progress(30,pc); /* progress is only estimated */
+//  if(job->cfg.verbose&32) debug_img("out16",job,4+8);
+
+  char_recognition( pp, job->cfg.mode);
+  progress(60,pc); /* progress is only estimated */
+//  if(job->cfg.verbose&32) debug_img("out17",job,4+8);
+
+  if ( adjust_text_lines( pp, job->cfg.mode ) ) { /* correct using chars */
+    /* may be, characters/pictures have changed line number */
+    list_sort(&(job->res.boxlist), sort_box_func);
+    // 2nd recognition call if lines are adjusted
+    char_recognition( pp, job->cfg.mode);
+  }
+
+#define BlownUpDrawing 1     /* german: Explosionszeichnung, temporarly */
+#if     BlownUpDrawing == 1  /* german: Explosionszeichnung */
+{ /* just for debugging */
+  int i,ii,ni; struct box *box2;
+  i=ii=ni=0;
+  for_each_data(&(job->res.boxlist)) { /* count boxes */
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    if (box2->c==UNKNOWN)  i++;
+    if (box2->c==PICTURE) ii++;
+    ni++;
+  } end_for_each(&(job->res.boxlist)); 
+  if (job->cfg.verbose)
+    fprintf(stderr,"# debug: unknown= %d picts= %d boxes= %d\n",i,ii,ni);
+}
+#endif
+  // ----------- write out20.pgm ----------- mark lines + boxes
+  if (job->cfg.verbose&32) debug_img("out20",job,1+4+8);
+
+  compare_unknown_with_known_chars( pp, job->cfg.mode);
+  progress(70,pc); /* progress is only estimated */
+
+  try_to_divide_boxes( pp, job->cfg.mode);
+  progress(80,pc); /* progress is only estimated */
+
+  /* --- list output ---- for debugging --- */
+  if (job->cfg.verbose&6) output_list(job);
+
+  /* ---- insert spaces ---- */
+  list_insert_spaces( pp , job );
+
+  // ---- proof difficult chars Il1 by context view ----
+  if (job->cfg.verbose)
+    fprintf(stderr,"# context correction if !(mode&32)\n");
+  if (!(job->cfg.mode&32)) context_correction( job );
+  
+  store_boxtree_lines( job, job->cfg.mode );
+  progress(90,pc); /* progress is only estimated */
+
+/* 0050002.pgm.gz ca. 109 digits, only 50 recognized (only in lines?)
+ * ./gocr -v 39 -m 56 -e - -m 4 -C 0-9 -f XML tmp0406/0050002.pbm.gz
+ *  awk 'BEGIN{num=0}/1<\/box>/{num++;}END{print num}' o
+ * 15*0 24*1 18*2 19*3 15*4 6*5 6*6 6*7 4*8 8*9 sum=125digits counted boxes
+ *  9*0 19*1 14*2 15*3 11*4 6*5 5*6 6*7 4*8 8*9 sum=97digits recognized
+ * 1*1 1*7 not recognized (Oct04)
+ *  33*SPC 76*NL = 109 spaces + 36*unknown sum=241 * 16 missed
+ */
+#if     BlownUpDrawing == 1  /* german: Explosionszeichnung */
+{ /* just for debugging */
+  int i,ii,ni; struct box *box2; const char *testc="0123456789ABCDEFGHIJK";
+    i=ii=ni=0;
+  for_each_data(&(job->res.boxlist)) { /* count boxes */
+    box2 = (struct box *)list_get_current(&(job->res.boxlist));
+    if (box2->c==UNKNOWN)  i++;
+    if (box2->c==PICTURE) ii++;
+    if (box2->c>' ' && box2->c<='z') ni++;
+  } end_for_each(&(job->res.boxlist)); 
+  if(job->cfg.verbose)
+    fprintf(stderr,"# debug: (_)= %d picts= %d chars= %d",i,ii,ni);
+  for (i=0;i<20;i++) {
+    ni=0;
+    for_each_data(&(job->res.boxlist)) { /* count boxes */
+      box2 = (struct box *)list_get_current(&(job->res.boxlist));
+      if (box2->c==testc[i]) ni++;
+    } end_for_each(&(job->res.boxlist)); 
+    if(job->cfg.verbose && ni>0)
+      fprintf(stderr," (%c)=%d",testc[i],ni);
+  }
+  if(job->cfg.verbose)
+    fprintf(stderr,"\n");
+}
+#endif
+
+  // ---- frame-size-histogram
+  // ---- (my own defined) distance between letters
+  // ---- write internal picture of textsite
+  // ----------- write out30.pgm -----------
+  if( job->cfg.verbose&32 ) debug_img("out30",job,2+4);
+    
+  progress(100,pc); /* progress is only estimated */
+
+  close_progress(pc);
+  
+  return 0; 	/* what should I return? error-state? num-of-chars? */
+}