Category

Development Platform
Home > SourceCode/Document > Graph program > View Code
HOUGH.C
Package: Algorithms for Image Processing and Computer Visio　[view]
Upload User: gzsenex
Upload Date: 2013-02-24
Package Size: 25984k
Code Size: 3k
Category:

Graph program

Development Platform:

C/C++

HOUGH.C：Code Content
  1. /* Estimate the skew angle: Hough transform of low resolution smoothed image */
  3. #define MAX
  4. #include "lib.h"
  5. #define OBJECT 1
  6. #define BACK 0
  8. struct lrec {
  9. int rs, re, cs, ce;
  10. } LINE[400];
  11. int Nlines = -1;
  12. int Space = 1;
  13. int Meanwidth = 0;
  15. typedef unsigned char ** GIMAGE;
  17. struct grec { /* Information concerning a glyph */
  18. char value;
  19. short int nr, nc;
  20. GIMAGE ptr;
  21. struct grec * next;
  22. };
  23. typedef struct grec * GPTR;
  24. GPTR database[256];
  26. void FDisplay (float ** x, int nr, int nc);
  27. void hough (IMAGE x, float *theta);
  29. void main (int argc, char *argv[])
  30. {
  31. IMAGE im;
  32. char text[128];
  33. int i, j;
  34. FILE *f;
  35. float theta;
  37. im = newimage (64, 64);
  38. for (i=0; i<64; i++)
  39.   for (j=0; j<64; j++)
  40. im->data[i][j] = BACK;
  41. im->data[2][2] = OBJECT;
  42. im->data[22][22] = OBJECT;
  43. im->data[55][55] = OBJECT;
  44. hough (im, &theta);
  45. printf ("Approximate skew angle is %10.5fn", theta-90);
  46. }
  48. void hough (IMAGE x, float *theta)
  49. {
  50. float **z;
  51. int center_x, center_y, r, omega, i, j, rmax, tmax;
  52. double conv;
  53. double sin(), cos(), sqrt();
  54. float tmval;
  56. conv = 3.1415926535/180.0;
  57. center_x = x->info->nc/2; center_y = x->info->nr/2;
  58. rmax = 
  59.  (int)(sqrt((double)(x->info->nc*x->info->nc+x->info->nr*x->info->nr))/2.0);
  61. /* Create an image for the Hough space - choose your own sampling */
  62. z = f2d (180, 2*rmax+1);
  64. for (r = 0; r < 2 * rmax+1; r++)
  65.    for (omega = 0; omega < 180; omega++)
  66.     z[omega][r] = 0;
  68. tmax = 0; tmval = 0;
  69. for (i = 0; i < x->info->nr; i++)
  70.   for (j = 0; j < x->info->nc; j++)
  71. if (x->data[i][j])
  72.    for (omega = 0; omega < 180; ++omega) 
  73.    {
  74.     r = (i - center_y) * sin((double)(omega*conv)) 
  75.    + (j - center_x) * cos((double)(omega*conv));
  76. /* if (r == 0) continue;  */
  77.     z[omega][rmax+r] += 1;
  78.    }
  80. for (i=0; i<180; i++)
  81.   for (j=0; j<2*rmax+1; j++)
  82.     if (z[i][j] > tmval)
  83.     {
  84. tmval = z[i][j];
  85. tmax = i;
  86.     }
  87. *theta = tmax;
  88. FDisplay (z, 180, 2*rmax+1);
  89. free (z[0]); free (z);
  90. }
  92. void FDisplay (float ** x, int nr, int nc)
  93. {
  94. float xmax, xmin, z, rng;
  95. int count=0, k=0;
  96. int i,j;
  97. FILE *f;
  99. xmax = xmin = x[0][0];
  100. for (i=0; i<nr; i++)
  101.   for (j=0; j<nc; j++)
  102.   {
  103.     if (xmax < x[i][j]) xmax = x[i][j];
  104.     if (xmin > x[i][j]) xmin = x[i][j];
  105.   }
  106. fprintf (stderr, "DISPLAY: xmax = %f xmin = %fn", xmax, xmin);
  107. rng = xmax - xmin;
  109. f = fopen ("/tmp/z", "w");
  110. fprintf (f, "P2n%d %dn255n", nc, nr);
  111. for (i=0; i<nr; i++)
  112.   for (j=0; j<nc; j++)
  113.   {
  114.     z = x[i][j];
  115.     k = (unsigned char) ((z-xmin)/rng * 255);
  116.     fprintf (f, "%3d ", k);
  117.     count++;
  118.     if (count > 30)
  119.     {
  120. count = 0;
  121. fprintf (f, "n");
  122.     }
  123.   }
  124. fclose (f);
  125. system ("xv /tmp/z");
  126. }